||
好久没更新了,这次准备公开共享自己博士期间独立研发的全球叶绿素产品生产代码。在去年写博士论文间隙,挤出时间帮导师独立研发了全球植物叶片叶绿素产品的生产代码。所用数据源如下:
数据源:
(1) MERIS 7天地表反射率数据,全年共52期数据
文件名:ESACCI-LC-L3-SR-MERIS-300m-P7D-h**v**-yyyymmdd-v2.0.nc
分辨率:(0.002778°*0.002778°)
投影:经纬度投影
(2) MCD12Q1 地表覆盖产品,一年一期
文件名:MCD12Q1.A2010001.h**v**.006.2018054195151.hdf
分辨率:500m*500m
投影:正弦投影
目标是需要同时用到MERIS和MODIS数据来进行反演。但考虑到MERIS数据与MODIS数据的投影和格网分幅管理方式各不相同,如何有效率的将其进行统一,并将其全部变为MODIS的格网管理方式尤为困难,这也是全球叶绿素产品反演生产的第一步也是最重要一步。经过几夜查询与思考,最终决定采用IDL平台来进行研发,并且巧妙利用其中的ENVISpatialGridRaster函数来统一两类数据的投影和格网管理格式。
由于当时只有4,5天时间去研发,其他时间需要完成博士论文,所以代码稍显凌乱,也没做注释。其中包含有相当多遥感数据应用方面的细节,例如MODIS的.hdf类型的数据读取,Meris的.nc类型数据的读取,如何给此类型数据加入投影,以及如何调用ENVISpatialGridRaster函数来统一投影和格网分幅管理等操作,可能会对一些初接触该类数据的人有较大帮助。现将代码全部公开,供各位需要的同行使用。
PS:此代码需要一个辅助的MODIS全球数据读取其投影信息(即代码部分的ref_file变量),此辅助数据可以是MODIS任何产品数据,但需为全球尺度。截图是其中部分代码内容,全部代码见附件。代码分为两部分:meris_retile.pro为MERIS和MODIS数据分幅管理格式的统一,globechlorophyll_retrieval.pro为叶绿素反演,已统一放入下面压缩包中。
chlorophyll_retrieval_v1.0.rar
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-1 16:21
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社