Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.------------------------------------------author------------------------------------------dateMODIS数据处理11. 数据下载MODIS 1B数据下载 L1B数据下载地址：http://ladsweb.nascom.nasa.gov/data/search.html 打开网页，选中需要的数据源类型图 1数据源类型选择注意：其中，MOD03数据是用于对1KM,QKM,HKM数据进行几何纠正所用，如果采用Modis Swath Tool等工具进行几何纠正，需下载MOD03数据接下来选择需下载的数据的时间图 2时间选择其中，日期类型为：月/日/年 时:分:秒需要注意的是，网页中显示的时间为UTC时间在‘spatial selection’选项中选择“latitude/longtitude”，按经纬度形式选择影像范围。
图 3范围选择点击’search’查到需要的数据然后点击’View RGB’预览所选的数据 图 4预览数据勾选所需要的数据，点击‘order files now’输入接收信息的邮箱，点’order’开始订购该数据如果要搜索多天数据，可以选‘add files to shopping cart’继续搜索其他日期的数据所订购数据的存放位置信息： 点击‘Data->Track Orders ’可以查看所有已订购的数据的状态如果’state’显示‘avalable’即可开始下载使用FTP下载软件下载如CUTEFTP，首先新建站点：ladsweb.nascom.nasa.govusername: anonymouspassword: seline808@ （自己申请数据的邮箱）连接成功后，开始下载数据申请的数据存放在orders>’订单号’下每个订单中的数据只会在FTP上存放5天，需要及时下载2. 几何校正2.1 Modis Swath Tool安装Modis Swath Tool是NASA网站提供的对HDF格式的1B数据进行几何精校正的工具，该软件使用MOD03数据对影像进行纠正，处理速度快且使用简单方便。
其安装步骤如下：(一) 获取MRTSwath软件包通过USGS EROS MODIS Reprojection Tool 网站（http://lpdaac.usgs.gov/landdaac/tools/mrtswath/index.asp）获取对应操作平台的Modis Swath Tool软件包二) 安装Java至少为java2运行环境1.3/1.4以上版本或者java2 SDK1.3以上版本三) 开始正式安装打开命令行窗口，进入安装文件所在目录四) 安装目录下运行install脚本文件，进入安装命令界面输入y，即同意安装五) 开始进行安装，按照所提示问题进行输入1）输入安装目标位置，如果想安装在当前目录下，直接回车即可，否则需要输入绝对路径（如：f:/MRTSwath)2）如果当前已经存在所输入目录则会提示是否替换，否则提示是否创建根据实际情况确定是否替换或创建3）解压缩完成后会提示：unzip executed successfully.并要求选择操作系统，选择对应数字4）输入java的bin文件夹所在路径(如：f:/java/bin)然后回车5）按回车键完成安装6）重启计算机六) 创建桌面快捷方式。
找到ModisSwathTool.bat发送到桌面快捷方式更改图标样式，右键图标＞属性＞更改图标＞浏览到ModisSwathTool图标需要注意的是：无论是java bin的路径还是安装目标路径中都不要出现空格如(c:/program files)，否则会造成GUI不能正常使用另外，在路径中以“/”代替“\”所有路径中最好都为英文，否则会出现不可预期的错误2.2 应用Modis Swath Tool进行校正（1） 打开ModisSwathTool，进入GUI界面图 5ModisSwathTool GUI界面（2） 点击“Open Input File…”，出现文件选择对话框，打开待校正的文件图 6打开待校正文件（3） 将Selected Bands中不需要的波段移至Available Bands中，保留需要的波段这里保留“EV_500_RefSB_b0”，“EV_500_RefSB_b1”，“EV_250_Aggr500_RefSB_b0”这三个波段（143真彩色合成）图 7选择波段（4） 点击“Open Geolocation File…”，打开相应的地理定位文件，即MOD03文件图 8加载地理定位文件（5） 在destination模块中单击“Specify Output File…”按钮，指定输出目录和文件名，这里选择原文件，输出文件名以原文件名为基准。
图 9指定输出目录及文件名（6） 分别在Output File Type和Resampling Type的下拉菜单中设置输出的文件格式和重采样的方法选择GEOTIFF格式和Nearest Neighbor（最近邻采样法）图 10设置输出文件类型及重采样方法（7） 分别在Output Projection Type和Output Data Type中设置输出的投影类型和数据类型选择Polar Stereographic投影方式图 11设置投影类型点击Edit Projection Parameters，弹出Projection Parameters对话框，输入投影参数这里所采用的投影是Antarctic Polar Stereographic，TrueScale参数应设为-71（标准纬线在南纬71度），SMajor = 6378137.0, Sminor = 6356752.314245179, LongPol = 0.0, FE = 0.0, FN = 0.0最后，在Output Data Type的下拉菜单中选择“Same As Input Data Types”，设置输出的数据类型与输入的数据类型相同。
图 12设置投影参数（8） 设置输出的像元大小这里根据输入影像的分辨率设置为500米9） 设置完毕后，单击“Run”开始校正，此时会弹出Status文本框，显示运行的状态图 13运行状态（10） 合并各个波段的校正结果影像校正完毕后，可以在前面给定的输出目录下看到各个选定的波段的校正结果影像在ENVI中打开这些文件，应用Basic Tools > Layer Stacking选项，按照“EV_250_Aggr500_RefSB_b0”，“EV_500_RefSB_b1”，“EV_500_RefSB_b0”的顺序将三个波段叠合成一个影像文件（ENVI格式）图 14校正结果图 15波段叠合3. 相对辐射校正待研究4. 影像拼接在ENVI中打开经过几何校正和合并的结果影像，打开Basic Tools > Mosaicking > Georeferenced模块单击菜单栏中的Import > Import Files and Edit Properties，选择需要加载的文件之后，在弹出的Entry对话框中设置Data Value to Ignore为65535（ModisSwathTool校正结果的默认背景值为65535），其它值均为默认。
载入后可以在对话框下部的列表中调整各个文件的叠合顺序调整完毕后单击File > Apply进行镶嵌，弹出Mosaic Parameters对话框，在此可以设置像元大小、重采样方法、输出目录及文件名以及背景值（设为0）单击OK运行，影像拼接过程完成最后，选择ENVI主菜单下File > Save file as > Erdas Image，将拼接结果影像存为Erdas Image格式以便在ArcMap中被加载图 16打开Mosaick模块图 17设置输入参数图 18加载影像文件后的对话框图 19 设置镶嵌参数图 20 执行镶嵌5. 影像裁剪5.1 应用ArcGIS ExtractByMask命令进行影像裁剪（1） 打开ArcMap，选择Windows > Command Line，显示Command Line窗口2） 加载前面影像拼接得到的结果影像（ERDAS IMAGE格式）以及所需用到的网格文件（Shapefile格式，共112个）3） 在Excel中编辑一系列所需的命令如下所示：ExtractByMask_sa 2008001_all.img 001 D:\data\Default\001\20080101.tifExtractByMask_sa 2008001_all.img 002 D:\data\Default\002\20080101.tifExtractByMask_sa 2008001_all.img 003 D:\data\Default\003\20080101.tifExtractByMask_sa 2008001_all.img 004 D:\data\Default\004\20080101.tifExtractByMask_sa 2008001_all.img 005 D:\data\Default\005\20080101.tifExtractByMask_sa 2008001_all.img 006 D:\data\Default\006\20080101.tif……ExtractByMask_sa 2008001_all.img 110 D:\data\Default\110\20080101.tifExtractByMask_sa 2008001_all.img 111 D:\data\Default\111\20080101.tifExtractByMask_sa 2008001_all.img 112 D:\data\Default\112\20080101.tif该命令的句法为：ExtractByMask_sa 。
在这里，in_raster为上面得到的拼接结果影像，in_mask_data为各个网格对应的shapefile（（2）中已加载）4） 将上述一系列命令复制到Command Line Window中，回车开始执行5.2 格式转换（TIFF to BMP）格式转换在IDL中实现具体代码如下所示：pro tiff_to_bmp;;convert tiff files to bmp files;envi,/restore_base_save_filesenvi_batch_init, log_file='batch.txt'for i=1,112 do begin ;判断并设置三种情况下的路径 if (i le 9) then begin path = 'D:\data\Default\。