《mrt处理modis数据资料》由会员分享，可在线阅读，更多相关《mrt处理modis数据资料（7页珍藏版）》请在金锄头文库上搜索。
1、转MODIS 产品下载与数据处理小结 2011-03-15 15:15原文链接：http:/ MODIS 产品下载与数据处理方面的一些内容，供想要入门的朋友参考。从 USGS 网站 http:/glovis.usgs.gov 网站可以浏览并下载需要的 MODIS 数据，不过另一些提供下载的网站还有：http:/e4eil01.cr.usgs.gov:22000/WebAccess/drill?attrib=home&next=group （最佳下载渠道）http:/nsidc.org/data/modis/order_data.html （拥有几乎全部方式的下载链接）https:/wist.e
2、cho.nasa.gov/api/ (WIST 搜索和下载 MODIS 产品的一个网站）ftp:/e4ftl01u.ecs.nasa.gov （FTP 下载方式的 MODIS 产品网站）其中 MODIS 产品的介绍可以从网站 https:/lpdaac.usgs.gov/lpdaac/products/modis_products_table 了解。无论何种方式下载，都需要等待很长时间，与网速有很大关系，因为需要一景一景的下载，文件大小根据所选数据产品的类型有关。一、MODIS 产品数据的打开数据下载后，在 ENVI、ERDAS IMAGINE 等遥感图像处理软件中都可以打开。如果下载的是HD
3、F 格式，在 ENVI 中直接使用 File - Open Image File命令打开，选择 HDF 格式就行。打开后，产品自动加载入“Available Band List”中，自动识别地理参考信息。在 ERDAS IMAGINE 中可以使用Import工具选择相应类型打开（这里不作介绍）。二、MODIS 产品数据的感兴趣提取在 ENVI 中打开后，需要对影像进行感兴趣区域的提取。这里只介绍一种方法，即使用已有的矢量多边形数据提取影像中的感兴趣区域。一般这种情况是根据某国家（省、市、地区等）行政边界提取某地区的影像数据。首先在选中的产品数据的打开窗口（Image Display）中添加矢量
4、边界层（Overlay-Vectors），然后选择Basic Tools中的Masking菜单下的Build Mask工具，打开后，首先选择你打开影像数据的窗口（Select Input Display），如果你影像在 Display #1 中打开，就选择选项中的 Display #1 即可。点击OK后跳出的对话框中的Options菜单下选择Import EVFs命令（当然，这个根据你定义感兴趣区域的数据选择，如果是 ROI，则选择Import ROIs），选中矢量边界文件（.evf 格式），再确定输出结果的方式File或Memory。一般我选后者。点击Apply后生成 MASK 文件。利用生
5、成的 MASK 文件使用Apply Mask工具进行掩膜操作，即只提取感兴趣区域。在Apply Mask Input File窗口中，选择原始影像文件，即打开的影像。选中文件后，进入Spatial Subset中，点击ROI/EVF按钮，选择感兴趣矢量文件，和刚才定义 MASK 所用的矢量 evf 文件相同。OK后推出，再点击Select Mask Band按钮，选择刚刚生成的MASK 文件，完成后，OK即可完成感兴趣影像的提取。忘记提到的是，MASK 掉的部分默认值为 0.0。三、MODIS 产品数据的重新投影MODIS 产品数据的投影一般是 Sinusoidal，不符合我们正常使用数据的查
6、看方式，因此需要进行重新投影。一般我选择 UTM，WGS84 形式的。以河南省为例，我通过上一步的感兴趣提取操作后，对感兴趣区域进行重新投影。在 ENVI 中的Map菜单下选择Convert Map Projection命令，进入对话框，点击Chang Proj.按钮，选择 UTM，Datum为WGS84，Zone为 50。OK后，再OK，完成重新投影。如果打开的影像发现是全白色的，只有一些区域是黑色的，那么你将红色框调整到某个位置，进行图像增强就可以了，如线性拉伸、直方图均衡等。如果输出的重投影影像看上去还是感兴趣小，图像大的时候，你可以再针对这个数据做一次 MASK 操作。记得Spatia
7、l Subset一定要定义好，和上同。今天就说到这里了，附上河南省的 MODIS 产品图（EVI 的）。四、MRT 批处理 MODIS 产品数据Modis Reprojection Tool 可以批处理 Modis 产品，如mosaic，resample，reprojection。在 windows 平台下可以在 Dos 命令框中运行自己创建的批处理*.bat 文件实现批处理。值得注意的是：MRT 工具读取文件的路径不可以包含中文字符，因此，请将文件保存在英文字母的路径下。这里只以 Global 0.05Deg MOD11C3 产品为例，提取 20012007 年北京地区的图像。使用的命令是
8、resample。Resample -p *.prm -I *.hdf -o *.tif*.prm 文件如下： INPUT_FILENAME = C:TDDOWNLOADMOD11C3MOD11C3.A2001001.005.2006360120227.hdfSPECTRAL_SUBSET = ( 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 )SPATIAL_SUBSET_TYPE = INPUT_LAT_LONGSPATIAL_SUBSET_UL_CORNER = ( 41.06 115.41 )SPATIAL_SUBSET_LR_CORNER = ( 39.42 1
9、17.5 ) OUTPUT_FILENAME = C:TDDOWNLOADMOD11C3outputMOD11C3.A2001001.005.tifRESAMPLING_TYPE = BILINEAROUTPUT_PROJECTION_TYPE = GEOOUTPUT_PROJECTION_PARAMETERS = ( 0.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.0 )DATUM = WGS84OUTPUT_PIXEL_SIZE = 0.005批处理命令可在 EXCEL 中制作完成，如图 五、使用 MCTK 进行 MODIS 数据
10、的批处理由于 MODIS 数据量比较大，有时候批处理使用第四部分的方法还是比较麻烦。可以使用MCTK 工具（MODIS CONVERSION TOOLKIT)在 IDL 编程环境下进行批处理。网站上有 MCTK 工具的下载链接。下面提供本人最近处理一批 MOD13Q1 数据的批处理代码（仅供参考）：PRO BATCH_MOD13Q1_PROCESSCOMPILE_OPT IDL2ENVI, /RESTORE_BASE_SAVE_FILESENVI_BATCH_INIT, LOG_FILE=BATCH.LOGPRINT, START : ,SYSTIME()ROOT_DIR = G:MODISH
11、ENAN_RICE_SEASONMOD13Q1FNS = FILE_SEARCH(ROOT_DIR,*.HDF,COUNT = COUNT)PRINT, There ara totally, COUNT, images.OUTPUT_LOCATION = G:MODISHENAN_RICE_SEASONMOD13Q1OUTPUT_LOCATIONGRID_NAME = MODIS_Grid_16DAY_250m_500m_VISD_NAME_EVI = 250m 16 days EVISD_NAME_QA = 250m 16 days VI Quality SD_NAME_NIR = 250m
12、 16 days NIR reflectanceSD_NAME_MIR = 250m 16 days MIR reflectance;OUTPUT_METHOD = 1 ;REPROJECTEDOUTPUT_METHOD = 0 ;Standard;投影转换设定;UNITS = ENVI_TRANSLATE_PROJECTION_UNITS(Meters);OUTPUT_PROJECTION = ENVI_PROJ_CREATE(/UTM,ZONE=49,UNITS=UNITS);OUTPUT_PS_X = 250;OUTPUT_PS_Y = 250;INTERPOLATION_METHOD
13、= 8 ;TRIANGULATION WITH NEAREST NEIGHBOR;导出 EVI, MIR, NIR, QAFOR i = 0, COUNT-1 DO BEGINFILENAME = FNSiA = STRPOS(FILENAME,.)OUTPUT_ROOT_NAME_NIR = NIR_+ STRMID(FILENAME,A+1,8)CONVERT_MODIS_DATA, IN_FILE = FILENAME, $OUT_PATH = OUTPUT_LOCATION, OUT_ROOT=OUTPUT_ROOT_NAME_NIR, $/HIGHER_PRODUCT, /GRID,
14、 GD_NAME=GRID_NAME,SD_NAME = SD_NAME_NIR, $OUT_METHOD = OUTPUT_METHOD, $;OUT_PROJ = OUTPUT_PROJECTION, $OUT_PS_X = OUTPUT_PS_X, OUT_PS_Y = OUTPUT_PS_Y, $NUM_X_PTS = 50, $NUM_Y_PTS=50, $;INTERP_METHOD = INTERPOLATION_METHOD, $BACKGROUND=0, FILL_REPLACE_VALUE=0,$R_FID_ARRAY=R_FID_ARRAY, R_FNAME_ARRAY=
15、R_FNAME_ARRAY, /NO_MSGOUTPUT_ROOT_NAME_MIR = MIR_ + STRMID(FILENAME,A+1,8)CONVERT_MODIS_DATA, IN_FILE = FILENAME, $OUT_PATH = OUTPUT_LOCATION, OUT_ROOT=OUTPUT_ROOT_NAME_MIR, $/HIGHER_PRODUCT, /GRID, GD_NAME=GRID_NAME,SD_NAME = SD_NAME_MIR, $OUT_METHOD = OUTPUT_METHOD, $;OUT_PROJ = OUTPUT_PROJECTION, $OUT_PS_X = OUTPUT_PS_X, OUT_PS_Y = OUTPUT_PS_Y, $NUM_X_PTS = 50, $NUM_Y_PTS=50, $;INTERP_METHOD = INTERPOLATION_METHOD, $BACKGROUND=0, FILL_REPLACE_VALUE=0,$R_FID_ARRAY=R_FID_ARRAY, R_FNAME_ARRAY=R_FNAME_ARRAY, /NO_MSGOUTPUT_ROOT_NAME_EVI = EVI_ + STRMID(FIL
