亚星游戏官网-www.yaxin868.com

专利名称：一种卫星光学遥感相机内方元素在轨检校方法
技术领域：
本发明属于遥感卫星在轨几何检校领域，特别是涉及一种卫星光学遥感相机内方 元素在轨检校方法。
背景技术：
随着国产遥感卫星空间分辨率的提高，遥感影像的内部几何精度成为一项重要的 考察指标。线阵推扫式光学遥感相机是一种主要的光学遥感卫星载荷，其内方位元素检校 能够提取光学畸变模型参数，减小畸变，提高影像内部几何精度。目前国内卫星遥感领域，在光学遥感相机的内方元素几何检校方面，除了在三线 阵立体成像相机的在轨动态检校方面有一定的研究外，其他还有对于面阵相机地面检校的 研究，其中三线阵相机动态检校的重点是研究三线阵相机的主点和焦距及三个相机之间的 几何关系；面阵相机地面检校则要借助地面测试设备如转台和光源发生器，并校正的是面 阵模型，但是对于单线阵推扫式相机在轨检校方面则缺少相应的研究。

发明内容
本发明的目的在于克服现有技术的上述不足，提供一种卫星光学遥感相机内方元 素在轨检校方法，该方法通过在轨检校光学遥感相机的几何参数，提高了遥感卫星图像的 内部几何精度，减小了光学畸变误差。本发明的上述目的是通过如下技术方案予以实现的一种卫星光学遥感相机内方元素在轨检校方法，包括如下步骤建立0XYZ三维坐标系，其中相机在坐标系中的位置为Oi，地球中心在坐标系中的 位置为02，坐标系中X向为线阵推扫方向，Y向为垂轨方向，Z向为遥感相机视轴，设地面某 控制点P，向量0P与地球相交于P点，0P与焦平面相交于P2点，由于光学部件设计与加工 缺陷的存在，视向量Of在焦平面的实际成像位置为Pi点，设P2点的Y坐标为y2，Pi点的Y 坐标为，则根据径向畸变模型公式得出y2与yi的关系式y2 = c0+c1y1+c2y12+c3y13计算模型参数C(l、Cl、c2和c3，将计算得到的模型参数C(l、Cl、c2和c3代入上式，得 到己点修正后的Y坐标y2'，进一步得到拟合残差Ay2 = y' 2_y2，调整相机焦距f，使得 拟合残差Ay2最�。�完成在轨检校，其中计算模型参数C(l、Cl、c2和c3的步骤如下(1)选择几何检校区，卫星在几何检校区成像后得到1级图像，同时选择与1级图 像具有相同分辨率或更高分辨率的正射影像和DEM高程图作为参考图，在1级图像上选取 k个控制点，同时在所述参考图上选取k个同名点，记录1级图像上k个控制点的坐标(m， n)，m表示列号，n表示行号，记录参考图上k个同名点的经纬度坐标和高度值(Lon，Lat, h)，Lon表示经度，Lat表示纬度，h表示高度，其中k为正整数；(2)根据列号m，计算0XYZ三维坐标系下控制点的y坐标，y = c(m-M/2)，其中M表示一行的像元个数，c表示像元尺寸，
得到控制点线阵向量《(|(0，7，小进一步归一化得到控制点线阵向量& ；(3)根据行号n，计算控制点所在行的成像时刻秒计数t，t = n d_tQ，其中d为积分时间，t0为1级图像第一行秒计数，根据秒计数t，利用卫星下传的GPS数据，计算成像时刻卫星在地固系的位置坐标 向量Pt，并假设卫星和相机的坐标重合，根据参考图上同名点的经纬度坐标和高度值(Lon， Lat, h)，计算得到地固系下控制点的坐标向量Pj，控制点坐标向量Pj减去卫星位置坐标向 量Pt得到控制点视向量％，进一步归一化得到控制点视向量Vi ；(4)在k个控制点的线阵向量Wi中选择最接近相机视轴的向量《。.，计算向量Wij 与视轴的夹角a = arctan(^)，Yo表示离相机视轴最近的控制点的y坐标；分别以最接近相机视轴的向量wu、ViJ为基准，计算其它k-1个向量与所述向量 wu、ViJ的夹角，得到k-1个控制点线阵向量间夹角AWi和k-1个控制点视向量间夹角AVi， 每个夹角都加上一个夹角a，得到Aw’，和Av’ x ；本步骤中以最接近相机视轴的向量Wij、ViJ为基准，并假设该向量的畸变可以忽 略，得到其它向量与基准向量wu、 之间的夹角后，两组夹角数据的差异反映了光学几何 畸变的特征，这是本方法中能够解耦卫星外方位元素影响的关键。(5)计算线阵坐标BWi和BVi，计算公式为Bwj = f*tan(Aw'BVi = f*tan (Av' ，f 为相机焦距；线阵坐标BWi和BVi分别对应公式y2 = c0+c1y1+c2y12+c3y13中的和y2，将BWi和 BVi代入公式y2 = (^(^…—(^/—(仏3，用最小二乘法拟合得到模型参数cQ、Cl、C2和c3。步骤(5)中的公式是将两组夹角数据转换为线阵Y坐标方向的数值，使得到的畸 变模型系数单位和Y坐标单位一致，方便应用该模型；对于直接利用角度进行几何校正的 程序，则可以直接进行角度畸变计算C(l、Cl、c2和c3，而模型中的yi和y2也对应于角度值。在上述卫星光学遥感相机内方元素在轨检校方法中，步骤(1)中几何检校区选择 人工地物目标丰富，道路交通发达的地区。在上述卫星光学遥感相机内方元素在轨检校方法中，步骤(1)中在1级图像上选 取沿Y方向均勻分布，沿X方向较窄区域中的k个控制点，并且20 < k < 40。本发明与现有技术相比具有如下优点(1)本发明在轨检校方法根据通常的光学镜头畸变模型建立了 3阶内方元素模 型，并基于相对角度误差受外方元素误差影响小的原理，从视向量之间的相对角度关系中 提取了内方元素的模型参数，实现了内外方元素解耦；(2)本发明在轨检校方法通过对光学高阶畸变建模，并通过最小二乘法拟合得到 模型参数，从而得到修正后的线阵坐标，并进一步得到拟合残差，通过调整焦距，使得拟合 残差最�。�完成几何检校，本发明方法能够达到较高的检校精度；(3)本发明在轨检校方法采用3阶模型，具有相当的柔韧性，对于不同的焦距误差 均能予以吸收，并通过模型参数平衡能够达到总体模型的高精度；(4)本发明在选择控制点时，选取沿像元排列方向均勻分布的控制点，并且控制点的选择区域为沿X推扫方向的较窄区域，并且越窄约好，目的是使得时间引起的姿态误差 达到最小。


图1为本发明单线阵推扫相机在轨成像模型；图2为本发明内方元素几何检校流程图；图3为本发明内方元素几何检校过程中控制点选取示意图；图4为本发明实施例中相机畸变拟合前后散点图；图5为本发明实施例中相机畸变拟合后残差图。具体实时方式下面结合附图通过具体实施例对本发明进行进一步详细的描述如图1所示为本发明单线阵推扫相机在轨成像模型，建立0XYZ三维坐标系，其中 航天相机在坐标系中的位置为，地球中心在坐标系中的位置为02，图1中左图表示线阵推 扫的情况，右图表示相机视向量和地球相交的情况。相机焦距为f，x为线阵推扫方向，Y为 垂轨方向，Z为相机视轴。设地面某控制点P，向量Of与地球相交于P点，与焦平面XY相 交于P2点，由于光学部件设计和加工的缺陷，光学相机总是存在一定的几何畸变，使得视向 量Of在焦平面的实际成像位置为Pi点。通常光学畸变为径向畸变，提出建立一个3阶多项式畸变模型来表达该径向畸 变。r' = c0+c1r+c2r2+c3r3(1)其中，r为畸变后实际成像距离即OPpr'为共线径向距离0P2。地球中心为02，卫星和地球的距离为OiOy向量01己在(》￥2相机坐标系下构成视向量《(0，72，0，己点的￥坐标为72。向量 在0XYZ相机坐标系下构成视向量v(0，yi，f)，Pi点的Y坐标为yi。根据公式(1)的径向畸变模型，和y2的关系为y2 = c0+c1y1+c2y12+c3y13(2)计算模型参数C(I、C1、C2和c3，将计算得到的模型参数C(I、C1、C2和c3代入公式⑵， 得到P2点修正后的Y坐标y2'，进一步得到拟合残差Ay2 = y' 2_y2，调整相机焦距f，使 得拟合残差A ^最�。�完成在轨检校。其中计算模型参数C(l、Cl、c2和c3的步骤如下步骤一控制点数据采集。选择一景人工地物目标丰富，道路交通发达的地区作为几何检校区，卫星在该区 域成像后得到1级图像，同时选择与1级图像同分辨率或更高分辨率的正射影像和DEM高 程图作为参考图。在经过了辐射校正的1级图像上选取控制点，同时在参考图上选取同名 点。记录1级图像上控制点坐标(m，n)，m表示列号，n表示行号，记录参考图上的经纬度坐 标和高度值(Lon, Lat, h)，Lon表示经度，Lat表示纬度，h表示高度。如图3所示为本发明内方元素几何检校过程中控制点选取示意图，每一个“ + ”表 示一个控制点，选择清晰交叉点，个数k个，20 < k < 40。控制点沿像元排列方向选择，均 勻分布；控制点的选择区域为沿X推扫方向的较窄区域，并且越窄约好，使得时间引起的姿
6态误差达到最小。步骤二 数据处理首先，根据列号m，计算相机坐标系下控制点的y坐标y = c (m-M/2)(3)这里，M表示一行的像元个数，c表示像元尺寸。得到控制点线阵向量％(0，y，f)， 进一步归一化得到控制点线阵向量^。如图2所示为本发明内方元素几何检校流程图。然后，根据行号n，可以计算控制点所在行的成像时刻秒计数t = n d-t0(4)其中d为积分时间，t0为该景第一行秒计数。根据卫星秒计数可以推算出成像时刻的卫星位置在地固系下的向量Pt，并假设卫 星质心和相机位置Oi重合。根据地面控制点的经纬度和高度值(Lon，Lat, h)，计算得到地固系下控制点坐标
向量Pj。控制点坐标向量？」减去卫星位置向量Pt得到控制点视向量V(l，并对该向量归一化 得到控制点视向量Vi。步骤三角度计算上述步骤二分别得到了归一化线阵向量&和归一化视向量Vi，对于k个控制点， 有线阵向量Wil k和视向量vn k。首先，在一系列线阵向量wn k中选择最接近相机视轴的向量Wij，计算其与视轴的 夹角a，该夹角计算公式为 = arctan(^)
J(5)其中y(1表示离相机视轴最近的控制点的y坐标，分别以最接近相机视轴的的向量wu、vu为基准，计算其他k_l个向量与所述向量 wu、ViJ的夹角，得到k-1个控制点线阵向量间夹角AWi和k-1个控制点视向量间夹角AVi， 每个夹角都加上一个常值角度a，得到新的Aw’ i*Av’ 10根据三角关系计算在0XYZ相机坐标系下的各个y坐标，即线阵坐标，得到BWi和 BVi，计算公式为Bwj = f*tan(Aw' ^Bvj = f*tan(Av'(6)线阵坐标BWi和BVi分别对应公式y2 = c0+c1y1+c2y12+c3y13中的Yl和y2，将BWi和 BVi代入y2 = Ca+CA+Cj^+CA3，用最小二乘法拟合得到模型参数cQ、Cl、c2和c3。基于模型参数C(l、Cl、c2和c3，可以得到修正后的P2点的Y坐标y2'，进一步得到 拟合残差Ay2 = y' 2-y2。调整相机焦距f，使得拟合残差A y2最�。�完成在轨检校。下面列举一个具体的实施例以我国环境减灾卫星1B的(XD2数据为参考，抽取2009年8月29日位于安徽和江 苏北部地区的一景遥感影像，进行试验计算。相机(XD像元尺寸0. 065mm,对应地面30米。 影像参考图选择landsat的ETM正射影像，高程数据为STRM90。控制点共有21个，平面精 度约为50米，高程精度约为5米。
如图4所示为本发明实施例相机畸变拟合前后散点图，由图4可知拟合前畸变误 差非常明显，其中‘☆’表示控制点畸变误差坐标，‘ + ’表示控制点误差拟合修正坐标，本实 施例分别选择了 5个焦距值，得到5组结果，见表1为畸变拟合结果图，图4的相机畸变拟 合前后散点图为选取表1中序号1数据得到的结果图。表1 从结果可以看出，3阶畸变模型得到的残差大约在0. 015mm,等于2. 3个像元，残差 分布见图5，说明残差已经没有明显规律。考虑参考点精度也大约在2个像元左右，因此说 明本发明在轨检校方法达到了预期精度，同时，本发明方法采用的3阶模型，具有相当的柔 韧性，对于不同的焦距差异，均能予以吸收，图5相机畸变拟合后各个控制点的残差图。以上所述，仅为本发明最佳的具体实施方式
，但本发明的保护范围并不局限于此， 任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换， 都应涵盖在本发明的保护范围之内。本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
权利要求
一种卫星光学遥感相机内方元素在轨检校方法，其特征在于包括如下步骤建立OXYZ三维坐标系，其中相机在坐标系中的位置为O1，地球中心在坐标系中的位置为O2，坐标系中X向为线阵推扫方向，Y向为垂轨方向，Z向为遥感相机视轴，设地面某控制点P，向量O1P与地球相交于P点，O1P与焦平面相交于P2点，由于光学部件设计与加工缺陷的存在，视向量O1P在焦平面的实际成像位置为P1点，设P2点的Y坐标为y2，P1点的Y坐标为y1，则根据径向畸变模型公式得出y2与y1的关系式y2＝c0+c1y1+c2y12+c3y13计算模型参数c0、c1、c2和c3，将计算得到的模型参数c0、c1、c2和c3代入上式，得到P2点修正后的Y坐标y2′，进一步得到拟合残差Δy2＝y′2-y2，调整相机焦距f，使得拟合残差Δy2最�。�完成在轨检校，其中计算模型参数c0、c1、c2和c3的步骤如下(1)选择几何检校区，卫星在几何检校区成像后得到1级图像，同时选择与1级图像具有相同分辨率或更高分辨率的正射影像和DEM高程图作为参考图，在1级图像上选取k个控制点，同时在所述参考图上选取k个同名点，记录1级图像上k个控制点的坐标(m，n)，m表示列号，n表示行号，记录参考图上k个同名点的经纬度坐标和高度值(Lon，Lat，h)，Lon表示经度，Lat表示纬度，h表示高度，其中k为正整数；(2)根据列号m，计算OXYZ三维坐标系下控制点的y坐标，y＝c(m-M/2)，其中M表示一行的像元个数，c表示像元尺寸，得到控制点线阵向量w0(0，y，f)，进一步归一化得到控制点线阵向量wi；(3)根据行号n，计算控制点所在行的成像时刻秒计数t，t＝n·d-t0，其中d为积分时间，t0为1级图像第一行秒计数，根据秒计数t，利用卫星下传的GPS数据，计算成像时刻卫星在地固系的位置坐标向量Pt，并假设卫星和相机的坐标重合，根据参考图上同名点的经纬度坐标和高度值(Lon，Lat，h)，计算得到地固系下控制点的坐标向量Pj，控制点坐标向量Pj减去卫星位置坐标向量Pt得到控制点视向量v0，进一步归一化得到控制点视向量vi；(4)在k个控制点的线阵向量wi中选择最接近相机视轴的向量wij，计算向量wij与视轴的夹角αy0表示离相机视轴最近的控制点的y坐标；分别以最接近相机视轴的向量wij、vij为基准，计算其它k-1个向量与所述向量wij、vij的夹角，得到k-1个控制点线阵向量间夹角Awi和k-1个控制点视向量间夹角Avi，每个夹角都加上一个夹角α，得到Aw’i和Av’l；(5)计算线阵坐标Bwi和Bvi，计算公式为Bwi＝f*tan(Aw′i)Bvi＝f*tan(Av′i)，f为相机焦距；线阵坐标Bwi和Bvi分别对应公式y2＝c0+c1y1+c2y12+c3y13中的y1和y2，将Bwi和Bvi代入公式y2＝c0+c1y1+c2y12+c3y13，用最小二乘法拟合得到模型参数c0、c1、c2和c3。FSA00000137999500021.tif
2.根据权利要求1所述的一种卫星光学遥感相机内方元素在轨检校方法，其特征在 于所述步骤(1)中几何检校区选择人工地物目标丰富，道路交通发达的地区。
3.根据权利要求1所述的一种卫星光学遥感相机内方元素在轨检校方法，其特征在于所述步骤(1)中在1级图像上选取沿Y方向均勻分布，沿X方向较窄区域中的k个控制 点，并且20 < k< 40。
全文摘要
本发明涉及一种卫星光学遥感相机内方元素在轨检校方法，该方法根据通常的光学镜头畸变模型建立了3阶内方元素模型，并基于相对角度误差受外方元素误差影响小的原理，从视向量之间的相对角度关系中提取了内方元素的模型参数，并得到拟合残差，通过调整相机焦距f，使得拟合残差最�。�完成在轨检校，本发明方法实现了内外方元素解耦，并通过在轨检校光学遥感相机的几何参数，提高了遥感卫星图像的内部几何精度，减小了光学畸变误差。
文档编号G01C11/02GK101852623SQ20101019488
公开日2010年10月6日 申请日期2010年6月2日 优先权日2010年6月2日
发明者徐建艳, 郝雪涛 申请人:中国资源卫星应用中心