实习三地形图的配准ppt课件

1、遥感图像的几何精纠正遥感图像的几何精纠正一、实验目的及内容一、实验目的及内容1.学会利用学会利用ERDAS软件对地形图进行配准。软件对地形图进行配准。2.学会利用学会利用ERDAS软件对遥感图像进行几何软件对遥感图像进行几何精纠正。精纠正。几何精纠正几何精纠正为什么要进行几何精纠正：遥感器本身引起的畸变 外部因素引起的畸变 处理过程中引起的畸变 二、地形图配准原理二、地形图配准原理 数字图象几何纠正：通过计算机对离散结构的数字图象几何纠正：通过计算机对离散结构的数字图像中的每一个像元逐个进行纠正处理的方法。数字图像中的每一个像元逐个进行纠正处理的方法。 这种方法能够精确地改正动态扫描图像的误差

2、。这种方法能够精确地改正动态扫描图像的误差。 基本原理：利用图像坐标和地面坐标另一图像基本原理：利用图像坐标和地面坐标另一图像坐标、地图坐标等之间的数学关系，即输入图像坐标、地图坐标等之间的数学关系，即输入图像和输出图像间的坐标转换关系实现。和输出图像间的坐标转换关系实现。 直接纠正方法:从原始图像，依次对每个像元根据变换函数 F（），求得它在新图像中的位置。并将灰度值付给新图像的对应位置上。间接纠正法:从新图像中依次每个像元，根据变换函数 f () 找到它在原始图像中的位置，并将图像的灰度值赋予新图像的像元。灰度重采样纠正后的新图像的每一个像元，根据变换函数，可以得到它在原始图像上的位置。如

3、果求得的位置为整数，则该位置处的像元灰度就是新图像的灰度值。如果位置不为整数，则有几种方法：1) 最近邻法2) 双线性内插法3三次卷积法 灰度重采样1) 最近邻法：距离实际位置最近的像元的灰度值作为输出图像像元的灰度值；灰度重采样2双线性法：以实际位置临近的4个像元值，确定输出像元的灰度值。公式为： 4141432144332211,iiiiinmpgpppppgpgpgpgpg）（ G2 G3 G4 G1 灰度重采样3三次卷积法以实际位置临近的16个像元值，确定输出像元的灰度值。公式为： 161161,iiiiinmpgpg）（ PSF 三次样条函数 sinc函数灰度重采样几种采样方法的优缺

4、点：1最近邻法：计算简单，一般情况下可接受，灰度级细微变化，保留原始数据，这种方法最大可产生半个像元的位置偏移，可能造成输出图像中某些地物不连续。2) 双线性插值：计算量中等，效果较好，细节丧失，具有平均化的滤波效果，边缘受到平滑作用，原来的像元值等到破坏，在以后的波谱识别分类分析中，会引起一些问题；3)三次卷积法:计算量较大，但它可既避免了采用最近邻近点法所带来地物不连续的现象，提供了比双线性内插法更加清晰的图像，具有均衡化和清晰化的效果。总上所述： 纠正的函数可有多种选择：多项式方法、共线方程纠正的函数可有多种选择：多项式方法、共线方程方法、随机场内插方法等等。其中多项式方法的应用最方法、

5、随机场内插方法等等。其中多项式方法的应用最为普遍。为普遍。 准准备备 工工作作 输输入入原原始始图图象象 建建立立纠纠正正函函数数 确确定定输输出出图图象象的的范范围围 逐逐个个像像元元进进行行几几何何变变化化 灰灰度度的的重重采采样样 输输出出纠纠正正后后的的图图象象 效效果果 评评价价 多项式纠正的基本思想：图像的变性规律可以看作是平移、缩放、旋转、仿射、偏扭、弯曲等形变的合成。一般的公式为： .54321022),(iiiiiiiixiYcXcYXcYcXccYXfx.54321022),(iiiiiiiiyiYdXdYXdYdXddYXfy.54321022),(iiiiiiiixiy

6、cxcyxcycxccyxFX.54321022),(iiiiiiiiyiydxdyxdydxddyxFY2) 直接法 1) 间接法 利用有限的控制点的已知坐标，解求多项式的系数，确定变换函数。然后讲各个像元带入多项式进行计算，得到纠正后的坐标。1) 1) 表征空间位置的可靠性，道路交叉点，标志物，水域表征空间位置的可靠性，道路交叉点，标志物，水域的边界，山顶，小岛中心，机场等。的边界，山顶，小岛中心，机场等。同名控制点要在图像上均匀分布；同名控制点要在图像上均匀分布；清楚辨认；清楚辨认；数量应当超过多项式系数的个数（数量应当超过多项式系数的个数（(n+1)(n+1)* *(n+2)/2(n+

7、2)/2）。当）。当控制点的个数超过多项式的系数个数时，采用最小控制点的个数超过多项式的系数个数时，采用最小2 2乘乘法进行系数的确定，使得到的系数最佳。法进行系数的确定，使得到的系数最佳。 控制点的选择原则控制点的选择原则:三、地形图配准操作步骤三、地形图配准操作步骤（一单幅图的配准步骤 经扫描后的地形图，并没有投影信息和坐标系统，在ERDAS软件下，利用公里格网或四个角点的经纬度座标给地形图赋上投影和坐标信息，建议利用公里格网点，因为一幅地形图的公里格网点较多，可以很好的纠正由于纸张变形，扫描等引起的误差。1、将扫描后的地形图转换为image格式、启动配准程序：Maindata preap

8、rationimage geometric correction3、选择要配准的地形图从下窗口中选择）4、选择几何校正模型 一般选择多项式变换模型polynomial）；单幅图多项式次方数可选1，两幅图多项式的次方数可数3或者3以上5、选择投影类型、填写投影参数在上图中选择add/change projection在上图中注意几点：1projection type:选择transverse mercator2)Spheroid name和datum:都选择:krasovsky3)Scale factor at central meridian:选择1.0004)Longitude of cen

9、tral meridian:中央子午线的经度需要计算5)False easting:北半球填写为50,000m(500公里)其他不变、选择几何校正采点模式 1单击set projection GCP tool2)单幅图的配准选择通过键盘方式输入控制点：即keyboard only 两幅图的配准选择exsiting viewer,即用一幅有坐标的图给另一幅没有坐标的图赋予坐标单位选择米、在图中选点，在此X residual 和Y residual输入坐标值，注意输入的单位为要换算成米，如图上的坐标为(19757,3571)，其中19757中19是带号，所以输入值为757000,3571应输入3571000单位为米）8、再采集其他控制点 最少控制点的个数可按下列公式计算：(t+1)*(t+2)/2，当次方数为时最少需要个控制点，次方需要个控制点，次方需要个控制点。 当输入3个控制点后本次实习），每采集一个地面控制点，系统自动产生一个ref.gcp9、其他知识：x residual y和 residual分别表示单个GCP控制点的X、Y残差， RMS error 和contribution是单个GCP的