97和98遥感图像的几何处理.docVIP

§9.7 线阵列传感器影像的共线方程纠正 如前所述，与缝隙摄影机一样，线阵列传感器采用推扫式成像，在每一成像时刻是中心投影，得到一条“线中心投影”影像，一幅图像有若干条线影像构成，形成一个多中心投影方式。其外方位元素随时间变化，不同线影像外方位元素之间是相关的，不可能把每条线影像的6个外方位元素都作为独立参数来求解。通常采用“中心投影加改正”的方法，即首先将一幅影像当作框幅式中心投影影像来看待，这时它只有一组外方位元素，实际上这组外方位元素相应于一个起算时刻（通常选择像幅的中心）的线影像；然后，考虑到传感器通常是装在卫星上的，具有很平稳的轨道和姿态变化率很小的特点，因此在某一范围内，可以认为外方位元素是近似随时间线性变化的。而且，当选取运行方向为影像上轴时，用像点在整幅图像上的坐标值表示该像点的成像时刻，这样，外方位元素随时间变化的影响就可以用下面的关系式来描述： （9.39） 式中： 、、、、、为像幅中心对应的线影像的外方位元素，下标表示起算时刻。 、、、、、为某一线影像的外方位元素，见图：线阵列图像的处理方法。例如为图中时刻的外方位元素。 、、、、、为外方位元素随时间线性变化率。 由以上分析可知，对于阵列传感器影像，可采用与框幅式中心投影影像相类似的方法，利用线中心投影的构像关系，并且考虑到外方位元素的线性变化特性（9.39）式，即可实现几何处理：包括（1）三维数字微分纠正，制作正射影像图。（2）双像解析摄影测量，立体定位获取地面三维坐标。 现写出线中心投影的构像关系： 若是旁向倾斜一个角摄影获得立体影像（公式推导详见§9.4的有关章节），在此给出其共线方程的全部内容 （9.25） 若是航向倾斜一个角摄影获得立体影像，（公式推导详见§9.4的有关章节），其共线方程的具体表达式为 （9.29） 一.几何纠正 1.影像坐标系的选择 对于由若干条线阵列影像构成的像幅，影像坐标系的原点设在像幅的中心，该点所在的扫描线方向为轴，垂直方向（飞行方向）为轴，构成右手直角坐标系，现以法国在1986年发射的卫星为例，每幅影像由6000???扫描线组成，每条线上又有6000个像元。这样，坐标原点就在第3000条线的第3000个像元上，第3000条线就是影像的轴，而各线上第3000个像元的连线就是轴。如图：卫星阵列影像坐标系定义。 这样，在由影像坐标与扫描坐标之间的转换（即内定向）时只有坐标原点的平移，而无旋转。 2.像点坐标的计算 由于线中心投影的特殊性，任意地面点所对应的外方位元素是未知的，它需由所对应的像点坐标按（9.39）式来计算，而像点坐标又是待定的。因此，不可能如同框幅式中心投影那样，直接由共线方程解算出像点的影像坐标，而应是一个迭代解算的过程。 迭代解算的基本思想是：先给定一近似的值（一般为0），代入（9.39）式计算出该点的外方位元素，然后将外方位元素代入（9.25）或（9.29）式的第一式，如果此时的是正确的，则（9.25）或（9.29）式的第一式，其等式应成立；否则利用等式右边算出的修正原近似值，即 （9.40） 式中，为迭代次数，重复此过程，直到不大于给定的限差，迭代计算结束。根据迭代计算出的以及由此计算的外方位元素、、、、、，由（9.25）或（9.29）式的第二式的右边计算出相应的改化坐标，即 由（9.25）式的第二式的右边计算： 可写出： (9.41) 或 由（9.29）式的第二式的右边计算： 可写出： （9.42） 由（9.41）或（9.42）式可分别解得，因此得到像点的真正影像坐标： （9.43） 或 （9.44） 3.确定纠正影像的灰度 上述计算得到的像点坐标除以像元尺寸（如卫星影像每个像元的大小为），并作坐标原点的平移后，作为双线性内插进行灰度重采样时的依据。从而得到应用间接法进行三维数字微分纠正的正射影像。 综上所述，当“影像坐标系的选择”和“像点坐标的计算”这两项工作完成之后，就可以对线阵列传感器影像进行几何纠正。通常采用以共线方程作为几何纠正的变换函数，在此就是（9.25）和（9.29）式，即为间接法数字微分纠正的变换关系式。共线方程（9.25）和（9.29）式右边给出了地面坐标，是三维几何纠正。此时需结合数字高程模型、影像的外方位元素、、、、、（包括起算时刻的6个元素、、、、、和6个元素的变化率、、、、、）、等效主距和倾角，由提供的每个格网中心点的地面坐标，按共线方程计算其相应于原始影像上的像点坐标，最后由像点坐标按双线性内插方法进行灰度重采样，求得该点的灰度值，并赋给纠正后的正射影像的相应像元。 二.CCD线阵列立体影像的立体定位 这里讨论