modis数据预处理

1、MODIS数据预处理1 .波段设置表IMODIS波段分布特征基本用途波段序号波段宽度nrft光谱灵徼度VV.;ni2信曝出陆地与云的界限i020-670_2LB_】2界同上2R-1I-S7624.7201陆地与云的性膜3459-47V35.3243同上4545-565290琢同上51230-12503.474同上6162S-165275邱同上72105-21551J0110海洋颜色.水体立层性成生物化学R405-2044.9期。同上g438-44841.9R3Sr同上HI483-49332.1X0?同上I1526536之7754同上12546-556ZLO750同上13662-6729591。

2、同上14573-6X38,7I0N7同上1574375310.25脑同上16S62-877516大气水分17K90-920I0.G167同上IS931-0413657同上1991596515.0250衣2MQD临波段分布特征（社，某本用途波段序号渡段宽度pm光漕昆通度V7m2-pm-sr等效噪声温度基地表/云温度2(13661)-3X4()0.45(300K)preprocessing-data_specificutilities-viewHDFdatasetattributes可以打开数据集里每个要素的属性表，在里面选中需要的HDF文件中的数据集，就会打开其属性表，波段内容如下Attribu

3、te3-5:band_names8,9,10,11,12,13lo,13hi,14lo,14hi,15,16,17,18,19,26”对应打开的HDF文件里1KM辐亮度文件的波段数，一共16个波段.其中13/14波段比较特殊，都有hi和lo两组数据，它们是传感器高敏感度和低敏感度两种状态下获取到的DN值,分别对应于较暗地物和较亮地物，使用哪个文件根据需要而定.但是在太湖湖区,13/14波段大部分区域效果都不太好.值会很大，出现溢出.可能是由于太湖的高浑浊度.2 .几何校正几何校正有三种方法：1）用envi自带模块进行几何校正，通过菜单栏中的Map-GeoreferencesMODIS选中env

4、i中已经打开的需要校正的数据集，输入研究区的地理位置，如下图左，投影用UTM,基准面用WGS-84,区域根据经纬度确定。输入完成，envi会自动校正,GLT2）用GLT即是查找表法对图像进行几何校正Map-GeoreferencefrominputGeometry-buidGLT并执行去蝴蝶结效应算法，有点是能对我们需要的那些波段进行校正。缺点也很明显。如下图右，校正结束的图像会失去原始图像四个角的信息，这样就无法和校正的图像很好的匹配起来，不利于一些后续的处理。用来建立查找表。在弹出的对话框中选择查找表的XY信息，其中X对应图像经度信息，Y对应纬度信息。然后只需要规定投影、基准面和区位信息，

5、就可以生成一个查找表文件。这个查找表文件的实质也是两幅图像，分别在每个像元上保存着经纬度值，并且像元位置是拉伸到我们规定的输出投影上面去了,那么剩下的矫正工作就只是把想要矫正的信息和查找表一一匹配起来，因此速度也很快。Map-GeoreferencefrominputGeometry-GeoreferencefromGLT就是上面所说的，把查找表运用到需要矫正的每一个像元上。这里只需要选择想要矫正的影响就好了。有一点需要注意：这里查找表的分辨率比辐亮度信息要低，因此如果直接用这个查找表对原始的辐亮度图像进行矫正，只能得到原始图像一小部分的信息。因为矫正时是按照像元行列号一一对应查询的。解决方案

6、有两个：可以把原始图像分辨率降低，重采样成和查找表一样大小的图像，再用查找表进行几何校正，这样运行速度快但是损失原始图像信息；也可以将查找表重采样成和原始的辐亮度图像一样大小的图像，再用这个新的查找表对图像进行矫正。这样就不会损失原始图像信息，但是计算速度会大大降低。视应用选择。3 )IDL批处理;forward_functionenvi_proj_createPROModis_gef_batchenvi,/restore_base_save_files恢复ENVIsav文件envi_batch_init,log_file=zbatch.tXt模式;inpath=DIALOG_PICKFILE

7、(/DIRECTORYTITLE=selectMODISfilespath)CD,inpathfilename=FILE_SEARCH(*.HDF);print,resultn=N_ELEMENTS(filename)outpath=DIALOG_PICKFILE(/DIRECTORYTITLE=selectMODISoutpath)FORi=0,n-1DOBEGINin_name=inpath+filenameiout_name=ReGeo+filenamei;设置校正方法;0=RadianceEmissivity,1=ReflectanceEmissivitycalib_method=1;

8、设置输出方法;0=Standard,1=Projected,2=StandardandProjectedout_method=1;设置输出投影output_projection=envi_proj_create(/geographic);在输出时设置去除蝴蝶效应convert_modis_data,in_file=in_name,out_path=outpath,$out_root=out_name,/l1b,out_method=out_method,$out_proj=output_projection,calib_method=calib_method,/bowtie,$sd_pos=1

9、,3,/no_msg,background=0.0而且是逐像元的拉伸。ENDFORenvi_batch_exitEND只能生成envi标准格式打开的一些文件，无法对具体的所有辐亮度信息进行矫正。3 .辐射定标辐射定标信息也通过打开HDF数据集属性信息查看basic_tools-preprocessing-data_specificutilities-viewHDFdatasetattributesAttribute3-7:radiance_offsets316.97219849,316.97219849,316.97219849,316.97219849,316.97219849,316.97

10、219849,316.97219849,316.97219849,316.97219849,316.97219849,316.97219849,316.97219849,316.97219849,316.97219849,316.97219849Attribute3-8:radiance_unitsWatts/mA2/micrometer/steradianAttribute3-9:reflectance_scales0.00002352,0.00001261,0.00000720,0.00000538,0.00000462,0.00000227,0.00000168,0.00000255,0

11、.00000141,0.00000232,0.00000292,0.00002401,0.00003249,0.00002493,0.00002598Offsets和scales规定了辐射定标的系数。Modis辐射定标公式是R=scales*（DN-offsets）其中R代表辐射定标之后的辐亮度,scales和offsets分别代表列表中的定标系数,DN代表图像原始数码值。4 .大气校正6s大气校正，最关键的就是获取卫星、太阳几何信息。Envi操作稍显繁琐，大致分为以下步骤：1）用上面的批处理程序对modis进行几何校正（modis处理工具包也可以），关键是获取到没有边角信息损失的几何校正影像

12、。因为后面要和几何信息进行叠加；2）用GLT方法对太阳几何、卫星几何信息进行波段合成（要用另存为envi格式影像。因为没有进行过几何校正，图像不包含经纬度信息，layerstacking时会报错，只能saveas成一幅没有经纬度信息的多波段图像，这一步是为了减少几何校正步骤）、几何校正；3）对几何校正之后的几何信息（四个波段的图像）进行重采样，得到和1）矫正得到的modis校准影像等行等列的图像；4）对3）、4）得到的结果进行波段叠加，就可以在湖心区，用查看z_profile的方式，查看该点的太阳几何、卫星几何信息。为了简化该步骤，可以用IDL辅助查找这些信息。大致思路和envi操作相似，不过

13、在查找湖心像元几何信息的部分做了一些改进，直接用经纬度进行匹配，可以减少一些步骤。（envi也可以这么操作）functionGLT_regeo,filename;,pos1,pos2,pos3,pos4;获取经纬度信息和需要进行几何校正的波段；输出文件6个波段信息;0:longtitude;1:latitude;2:sensor_zenith;3:sensor_azimuth;4:solar_zenith;5:solarazimuthoutname=strsplit(filename,.,/extract)ENVI_OPEN_DATA_FILE,filename,r_fid=fid_x,/hd

14、f_sd,$hdfsd_dataset=1,hdfsd_interleave=0ENVI_OPEN_DATA_FILE,filename,r_fid=fid_y,/hdf_sd,$hdfsd_dataset=0,hdfsd_interleave=0;ENVI_OPEN_DATA_FILE,filename,r_fid=fid_sensor_zenith,/hdf_sd,$hdfsd_dataset=13,hdfsd_interleave=0ENVI_OPEN_DATA_FILE,filename,r_fid=fid_sensor_azimuth,/hdf_sd,$hdfsd_dataset=1

15、4,hdfsd_interleave=0ENVI_OPEN_DATA_FILE,filename,r_fid=fid_solar_zenith,/hdf_sd,$hdfsd_dataset=16,hdfsd_interleave=0ENVI_OPEN_DATA_FILE,filename,r_fid=fid_solar_azimuth,/hdf_sd,$hdfsd_dataset=17,hdfsd_interleave=0;查询图像信息envi_file_query,fid_sensor_zenith,dims=dims;将四个有用信息的图像合成一个图像，便于后续的几何校正工作envi_doi

16、t,cf_doit,$fid=fid_sensor_zenith,fid_sensor_azimuth,fid_solar_zenith,fid_solar_azimuth,$pos=0,0,0,0,dims=dims,$remove=0,$;out_name=out_name,$r_fid=info_fid,/in_memory;将四个单独的几何信息文件释放envi_file_mng,id=fid_sensor_zenith,/removeenvi_file_mng,id=fid_sensor_azimuth,/removeenvi_file_mng,id=fid_solar_zenith,

17、/removeenvi_file_mng,id=fid_solar_azimuth,/remove;建立查找表i_proj=envi_proj_create(/geographic,datum=WGS-84)o_proj=envi_proj_create(/utm,datum=WGS-84,zone=50)envi_glt_doit,i_proj=i_proj,$o_proj=o_proj,$rotation=0,$;保证图像北面朝上r_fid=glt_fid,$x_fid=fid_x,y_fid=fid_y,$out_name=outname0+_GLT,$x_pos=0,y_pos=0;,

18、/in_memory;根据查找表对数据进行几何校正envi_doit,envi_georef_from_glt_doit,fid=info_fid,$glt_fid=glt_fid,out_name=outname0+_regeo,pos=0,1,2,3,$subset=dims,r_fid=regeo_fidenvi_file_query,regeo_fid,dims=dimsenvifilemng,id=gltfid,/removeenvi_file_mng,id=info_fid,/removeenvi_file_mng,id=fid_x,/removeenvi_file_mng,id=

19、fid_y,/remove;将最后生成数据的id号，和最后图像的维度信息，存成一个数组返回print,regeo_fidreturn,regeo_fid,dimsend;proget_modis_infofile_name=DIALOG_PICKFILE(TITLE=selectaHDFfile)ifstrlen(file_name)eq0thenreturninfo=GLT_regeo(file_name);info0:函数输出图像的id值;info1:5:函数输出图像的纬度信息;用记录下的太湖图像中心部分的经纬度，转化成投影坐标，再转化成图像坐标,对返回的图像做行列号匹配，即是最佳位置。;地理坐标转化为投影坐标center_lat=31.22723889center_lon=120.2057722datum=WGS-84iproj=ENV