遥感图像预处理汇总.docx

遥感图像预处理汇总 1.自定义坐标系 1.1.坐标系概括 常用坐标系分为地理坐标系和投影坐标系 1.1.1.地理坐标系 地理坐标系中含椭球体(spheroid)和大地基准面(datum)两个重要部分 1)椭球体 地球并不是一个规则球体，因此我国常用以下椭球体模拟地球，见表1 2)大地基准面 大地基准面指参考椭球与WGS84参考椭球的相对位置关系（3个平移，3个旋转，1个缩放），可用3个、4个或者7个参数进行描述 每个椭球体都对应一个或多个大地基准面 1.1. 2.投影坐标系 投影坐标系是利用一定数学法则将地球表面上的经纬线网表示到平面上，属于平面坐标系一个投影坐标系包含椭球体、投影方法和大地基准面 我国普遍采用的投影方法是高斯-克吕格投影，在欧美国家称为横轴莫卡托投影(Transverse Mercator)为保证地图精度，采用分带投影法，有3°分带和6°分带两种对于每分带的投影，需指定中央经线以及坐标轴原点X、Y坐标我国在北半球，为避免Y坐标出现负值，规定统一将各带Y值加500km。

为避免各带坐标重复，可在不同带间Y值冠以带号，即通用坐标 商业软件坐标系分为标准坐标系和自定义坐标系，我国实际应用中往往需要自定义坐标系，以北京54坐标系和西安80坐标系为例讲解如何在Envi中自定义坐标系 1.2.Envi自定义坐标系 1.2.1.定义椭球体 使用记事本打开HOME\ITT\IDL80\products\envi48\map_proj\ellipse.txt，将“Krasovsky, 6378245.0, 6356863.3”(Beijing-54坐标系使用)和 “IAG-75, 6378140.0, 6356755.3”(Xian-80坐标系使用)加入文本末端 1.2.2.定义基准面 使用记事本打开HOME\ITT\IDL80\products\envi48\map_proj\datum.txt，将“D_BEIJING_1954, Krasovsky, -12, -113, -41”(Beijing-54坐标系使用)和“D_Xian_1980, IAG-75, 0, 0, 0”(Xian-80坐标系使用)加入文本末端。

注：此处的D_BEIJING_1954和D_Xian_1980在Envi 4.7版本以上不可更改，否则会转换失败 1.2.3.定义投影系 在Envi任何用到投影坐标的功能模块中，都可以新建坐标系（在任何地图投影选择对话框中，点击New按钮），也可通过主菜单→Map→Customize Map Projection，弹出Customized Map Projection Definition对话框，如图1 Beijing 1954 6 Degree 111E 图 1 定义投影系 依序定义如下参数： 1)Projection Name：投影坐标系名称； 2)Projection Type：选择投影类型，如Transverse Mercator； 3)Projection Datum：选择基准面； 4)False easting：东偏距离500000，若需X轴含带号，添加相应带号于500000 前； 5)False northing：北偏距离，填写0； 6)Latitude：中央纬线，填写0； 7)Longitude：中央经线，填写某带的中央经线，如117； 8)Scale factor：中央经线长度比，填写0.9996，精度与高斯-克吕格等同。

选择Projection→Add New Projection，将投影添加到Envi使用的列表中 选择File→Save Projection，存储新添或更改过的投影信息 打开HOME\ITT\IDL71\products\envi48\map_proj\map_proj.txt，即可看到新建的坐标信息，如图2至此，投影系已添加完毕，可在Envi中使用 图 2 添加完成的投影坐标系 1.3.图像投影转换 描述一个栅格文件地理位置信息由两部分组成：坐标信息（map）和投影信息（projection）坐标信息由起始点像素坐标以及对应的地理（投影）坐标和像素大小组成；投影信息就是坐标系信息一般来说，如果坐标信息丢失，这个文件将会失去坐标；投影信息是用来描述坐标信息，如果投影信息丢失，可以重新设定 下面演示如何将北京54坐标系的栅格图像转化为西安80坐标系 1.3.1.定义初始投影系 1)打开.\data\ 0-坐标定义\ f49e011021.img文件； 2)在Available Bands List中，Map Info节点右击，单击Edit Map Information， 弹出对话框； 3)单击Change Projection按钮，在弹出对话框中选择目标投影系，如图3； 图 3 设定影像投影系 4)点击OK，Available Bands List中的Map Info下拉节点坐标信息发生了相 应改变，表明投影信息已成功赋予图像，如图4。

图 4 设定影像投影系前后对比 1.3. 2.投影转换 选择主菜单→Map→Convert Map Projection，在弹出的对话框中设置如下参数，如图5： 1)Change Projection：选择投影目标系名称； 2)Conversion Parameters：Method选用Polynomial，Polynomial Degree为2； 3)Output Result to：选择Memory 图 5 投影转换参数设置单击OK后，执行投影转换，转换结果如图6 图 6 投影转换后坐标信息 2.图像几何校正 几何校正是利用地面控制点和几何校正数学模型来矫正非系统因素产生的误差，同时也是将图像投影到平面上，使其符合地图投影系统的过程；由于校正过程中会将坐标系统赋予图像数据，所以此过程包括了地理编码(geo-coding) 2.1.基于自带定位信息的几何校正 对于重返周期短、空间分辨率较低的卫星数据，如A VHRR、MODIS、SeaWiFS 等，地面控制点的选择有相当难度。

因此，可以用卫星传感器自带的地理定位文件进行几何校正 在Envi中，基本方法为主菜单→Map→Georeference传感器名称本部分以MODIS Level 1B级数据介绍基于自带定位信息的几何校正方法 2.1.1.打开数据文件 打开.\data\1-Modis\ MOD02HKM.A2022248.0345.005.2022348121959.hdf文件 2.1.2.选择校正模型 选择主菜单→Map→Georeference MODIS，弹出Input MODIS File对话框（图7），单击文件名，选择校正的文件单击OK，进入下一步 图7 选择要处理的MODIS波段 2.1. 3.设置输出参数 a)在Georeference MODIS Parameters对话框中（图4），设置输出坐标系； b)在Number Warp Points中，键入X、Y方向校正点的数量在X方向的 校正点数量应该小于等于51，在Y方向的校正点数应该小于等于行数； c)Enter Output GCP Filename：设置控制点输出文件； d)Perform Bow Tie Correction：用来消除MODIS的“蝴蝶效应”，默认Yes。

如图8，点击OK，进入下一步 图8 选择投影坐标系 2.1.4.设置校正参数 a)系统自动计算起始点坐标值、像元大小、图像行列数据，也可根据要求 更改； b)Background设置为0； c)设置输出影像路径 图9 校正参数及输出路径 如图9，点击OK，开始执行校正，经过图10两个步骤后，最终校正结果 如图11所示 图10 影像校正 图11 校正前(左)后(右)影像 2.2.基于GLT的FY3数据校正 GLT几何校正法利用输入的集合文件生成一个GLT(geographic lookup table)地理位置查找表文件它是一个二维图像文件，包含地理校正图像的行列两个波段，文件对应的灰度值表示原始图像每个像素对应的地理位置坐标信息，用符号整型存储，符号为正时，说明使用了真实的像元位置值，符号为负时，说明使用了邻近像元的位置值，值为0说明周围7个像元内没有邻近像元位置值 GLT文件包含初始图像每个像元的地理定位信息，它的校正精度是很高的， 避免了通过地面控制点利用二次多项式几何校正法对低分辨率图像数据的处理。

下面以我国(FY3)气象卫星可见光红外扫描辐射计(Visible and InfraRed Radiometer, VIRR)数据为例，进行GLT几何校正 2.2.1.安装补丁 安装HDF5补丁，将.\data\2-FY3中的open_hdf5_event.sav文件拷贝到HOME\ITT\IDL80\products\envi48\save_add文件夹中，重启ENVI 2.2.2.打开文件 a)主菜单→File→Open External File→Generic Formats→HDF5，选择文件 Z_SATE_C_BAWX_20220104070730_P_FY3A_VIRRX_GBAL_L1_2022 0104_0510_1000M_MS.HDF（文件类型选择*.*） b)如图12，选择EV_RefSB图像数据，点击Import to ENVI 图12 HDF5文件查看 c)同上步，将定位经纬度文件打开(Latitude和Longitude) 2.2. 3.生成GLT文件 a)选择主菜单→Map→Georeference from Imput Geometry→Build GLT。

b)在Input X Geometry Band对话框中，选择精度Longitude作为X波段， 如图13 图13 选择X波段 c)由于X波段左边边缘为0值，因此有必要对边缘进行掩膜处理，点击 Spatial Subset按钮，在Select Spatial Subset对话框中，设置Samples为3 （去掉开始3。