MODIS图像下载及处理.docxVIP

[转]MODIS下载与处理 ? 转自丫丫博客 从USGS网站  HYPERLINK / 网站可以浏览并下载需要的MODIS数据，不过另一些提供下载的网站还有：  HYPERLINK /wist-bin/api/ims.cgi/u228616 /wist-bin/api/ims.cgi/u228616 (WIST 搜索和下载MODIS产品的一个网站） （FTP下载方式的MODIS产品网站） 其中MODIS产品的介绍可以从网站  HYPERLINK /lpdaac/products/modis_products_table /lpdaac/products/modis_products_table 了解。 无论何种方式下载，都需要等待很长时间，与网速有很大关系，因为需要一景一景的下载，文件大小根据所选数据产品的类型有关。 一、MODIS产品数据的打开 数据下载后，在ENVI、ERDAS IMAGINE等遥感图像处理软件中都可以打开。如果下载的是HDF格式，在ENVI中直接使用 [File] ---- [Open Image File]命令打开，选择HDF格式就行。打开后，产品自动加载入“Available Band List”中，自动识别地理参考信息。在ERDAS IMAGINE中可以使用[Import]工具选择相应类型打开（这里不作介绍）。 二、MODIS产品数据的感兴趣提取 在ENVI中打开后，需要对影像进行感兴趣区域的提取。这里只介绍一种方法，即使用已有的矢量多边形数据提取影像中的感兴趣区域。一般这种情况是根据某国家（省、市、地区等）行政边界提取某地区的影像数据。首先在选中的产品数据的打开窗口（Image Display）中添加矢量边界层（[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文件。利用生成的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，不符合我们正常使用数据的查看方式，因此需要进行重新投影。一般我选择UTM，WGS84形式的。以河南省为例，我通过上一步的感兴趣提取操作后，对感兴趣区域进行重新投影。在ENVI中的[Map]菜单下选择[Convert Map Projection]命令，进入对话框，点击[Chang Proj...]按钮，选择UTM，[Datum]为WGS84，[Zone]为50。[OK]后，再[OK]，完成重新投影。如果打开的影像发现是全白色的，只有一些区域是黑色的，那么你将红色框调整到某个位置，进行图像增强就可以了，如线性拉伸、直方图均衡等。 如果输出的重投影影像看上去还是感兴趣小，图像大的时候，你可以再针对这个数据做一次MASK操作。记得[Spatial Subset]一定要定义好，和上同。 今天就说到这里了，附上河南省的MODIS产品图（EVI的）。 四、MRT批处理MODIS产品数据 Modis Reprojection Tool 可以批处理Modis产品，如mosaic，resample，reprojection。在windows平台下可以在Dos命令框中运行自己创建的批处理*.bat文件实现批处理。值得注意的是：MRT工具读取文件的路径不可以包含中文字符，因此，请将文件保存在英文字母的路径下。这里只以Global 0.05Deg MOD11C3产品为例，提取2001~2007年北京地区的图像。使用的命令是resample。 Resample -p *.prm -I *.hdf -o *.