《MODIS数据处理.doc》由会员分享，可在线阅读，更多相关《MODIS数据处理.doc（39页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、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等工具进行几何纠正，需下载
2、MOD03数据。接下来选择需下载的数据的时间。图 2时间选择其中，日期类型为：月/日/年 时:分:秒。需要注意的是，网页中显示的时间为UTC时间。在spatial selection选项中选择“latitude/longtitude”，按经纬度形式选择影像范围。图 3范围选择点击search查到需要的数据。然后点击View RGB预览所选的数据。 图 4预览数据勾选所需要的数据，点击order files now。输入接收信息的邮箱，点order开始订购该数据。如果要搜索多天数据，可以选add files to shopping cart继续搜索其他日期的数据。所订购数据的存放位置信息： 点击
3、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数据进行几何精校正的工具，该软件使用M
4、OD03数据对影像进行纠正，处理速度快且使用简单方便。其安装步骤如下：(一) 获取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，即同意安
5、装。(五) 开始进行安装，按照所提示问题进行输入。（1）输入安装目标位置，如果想安装在当前目录下，直接回车即可，否则需要输入绝对路径（如：f:/MRTSwath)。（2）如果当前已经存在所输入目录则会提示是否替换，否则提示是否创建。根据实际情况确定是否替换或创建。（3）解压缩完成后会提示：unzip executed successfully.并要求选择操作系统，选择对应数字。（4）输入java的bin文件夹所在路径(如：f:/java/bin)然后回车。（5）按回车键完成安装。（6）重启计算机。(六) 创建桌面快捷方式。找到ModisSwathTool.bat发送到桌面快捷方式。更改图标样式
6、，右键图标属性更改图标浏览到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中不需要的波
7、段移至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和Resampl
8、ing Type的下拉菜单中设置输出的文件格式和重采样的方法。选择GEOTIFF格式和Nearest Neighbor（最近邻采样法）。图 10设置输出文件类型及重采样方法（7） 分别在Output Projection Type和Output Data Type中设置输出的投影类型和数据类型。选择Polar Stereographic投影方式。图 11设置投影类型点击Edit Projection Parameters，弹出Projection Parameters对话框，输入投影参数。这里所采用的投影是Antarctic Polar Stereographic，TrueScale参数应设为
9、-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） 合并各个波段的校正结果影像。校正完毕后，可以在前面
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 Propert
11、ies，选择需要加载的文件之后，在弹出的Entry对话框中设置Data Value to Ignore为65535（ModisSwathTool校正结果的默认背景值为65535），其它值均为默认。载入后可以在对话框下部的列表中调整各个文件的叠合顺序。调整完毕后单击File Apply进行镶嵌，弹出Mosaic Parameters对话框，在此可以设置像元大小、重采样方法、输出目录及文件名以及背景值（设为0）。单击OK运行，影像拼接过程完成。最后，选择ENVI主菜单下File Save file as Erdas Image，将拼接结果影像存为Erdas Image格式以便在ArcMap中被加载
12、。图 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:dataDefault0012
13、0080101.tifExtractByMask_sa 2008001_all.img 002 D:dataDefault00220080101.tifExtractByMask_sa 2008001_all.img 003 D:dataDefault00320080101.tifExtractByMask_sa 2008001_all.img 004 D:dataDefault00420080101.tifExtractByMask_sa 2008001_all.img 005 D:dataDefault00520080101.tifExtractByMask_sa 2008001_all.
14、img 006 D:dataDefault00620080101.tifExtractByMask_sa 2008001_all.img 110 D:dataDefault11020080101.tifExtractByMask_sa 2008001_all.img 111 D:dataDefault11120080101.tifExtractByMask_sa 2008001_all.img 112 D:dataDefault11220080101.tif该命令的句法为：ExtractByMask_sa 。在这里，in_raster为上面得到的拼接结果影像，in_mask_data为各个网格
15、对应的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.txtfor i=1,112 do begin ;判断并设置三种情况下的路径 if (i le 9) then begin path = D:dataDefault00+
16、strtrim(string(i),1) endif else if (i le 99) then begin path = D:dataDefault0+strtrim(string(i),2) endif else begin path = D:dataDefault+strtrim(string(i),3) endelse ;打开第i个分块的TIFF格式影像 fname = path +20080101.tif envi_open_data_file, fname, /tiff, r_fid=fid if (fid eq -1) then begin envi_batch_exit re
17、turn endif ;获取文件信息 envi_file_query, fid, dims=dims, nb=nb pos = lindgen(nb) ;设置输出的临时文件名 out_name = path +tempimg ;对影像进行2%线性拉伸 envi_doit, stretch_doit, $ fid=fid, pos=pos, dims=dims, $ method=1, out_name=out_name, $ i_min=2.0, i_max=98.0, range_by=0, $ out_min=0, out_max=255, out_dt=12, $ r_fid=r_fid
18、 ;打开经过拉伸的影像 envi_open_file, out_name, r_fid=fid envi_file_query, fid, dims=dims, nb=nb, ns=ns, nl=nl ;获取各个波段的数据 dataR = ENVI_GET_DATA(fid=fid, dims=dims, pos=0) dataG = ENVI_GET_DATA(fid=fid, dims=dims, pos=1) dataB = ENVI_GET_DATA(fid=fid, dims=dims, pos=2) ;将三个波段的数据分别存入待写出的数组中 image = uintarr(3,nl,ns) image0,*,* = rotate(dataR,7) image1,*,* = rotate(dataG,7) image2,*,* = rotate(dataB,7) ;将image数组写出为bmp文件 write_image, path +20080101.bmp, bmp, image,/append ;删除临时文件 file_delete, path +tempimg, path +tempimg.hdrendforenvi_batch_exitend-