一种基于密集控制点的高分辨率卫星姿态颤振探测方法技术

本发明专利技术涉及一种基于密集控制点的高分辨率卫星姿态颤振探测方法，包括以下步骤：1)参考影像预处理；2)预处理后的参考影像和实时获取的卫星影像的密集匹配；3)卫星姿态颤振的解算。与现有技术相比，本发明专利技术利用外部参考数据(数字正射影像DOM和数字高程模型DEM)的支持，能够准确地探测出卫星平台的姿态颤振现象，从而还原卫星在轨运动的真实状态，具有准确性高、探测精度高等优点。

【技术实现步骤摘要】
一种基于密集控制点的高分辨率卫星姿态颤振探测方法
本专利技术涉及卫星姿态颤振探测方法，尤其是涉及一种基于密集控制点的高分辨率卫星姿态颤振探测方法。
技术介绍
从90年代开始，对地观测卫星对精度和稳定性的要求越来越高，国内外开始重视和研究由各种原因造成的卫星平台的波动。与其他成像方式的传感器相比，线阵CCD传感器在推扫成像的过程中更容易受到卫星运动和旋转的影响。卫星本身的机械运动如太阳帆板、动量轮和天线等，星体的受热扭曲以及姿态控制等操作，都会引起姿态的变化。当卫星的测姿设备的测量频率低于姿态变化的频率时，无法从姿态观测数据中反演其高频的变化，导致卫星影像的几何处理精度降低。因此，在遥感影像的几何处理之前，必须先去除姿态波动造成的影响。
技术实现思路
本专利技术的目的就是为了克服上述现有技术存在的缺陷而提供一种能够准确地探测出卫星平台的姿态颤振现象的基于密集控制点的高分辨率卫星姿态颤振探测方法。本专利技术的目的可以通过以下技术方案来实现：一种基于密集控制点的高分辨率卫星姿态颤振探测方法，包括以下步骤：1)参考影像预处理；2)预处理后的参考影像和实时获取的卫星影像的密集匹配；3)卫星姿态颤振的解算。步骤1)中，所述预处理包括高斯滤波和降采样处理。步骤2)中，所述密集匹配具体为：201)从预处理后的参考影像中提取出待匹配特征点，并根据相应的DEM获得待匹配特征点的地面坐标，所述待匹配特征点的采点间隔小于5*5像素区域；202)将所述待匹配特征点通过所述卫星影像的严格成像模型向后投影至卫星影像上，获得匹配点的初始值；203)根据匹配点初始值使用最小二乘匹配的方法，得到卫星影像和预处理之后的参考影像之间的密集控制点。所述密集匹配还包括：采用相关系数阈值法和基于严格成像模型和仿射变换模型的几何约束条件的剔除策略剔除误匹配点。所述严格成像模型为：其中，(X,Y,Z)为地面点在物方空间坐标系中的坐标，(XS,YS,ZS)为卫星在在物方空间坐标系中的位置，λ为比例因子，(x,y)为像方空间坐标，f为相机焦距，是由地球惯性系到WGS84坐标系的旋转正交矩阵，是由卫星轨道坐标系到地球惯性系的旋转正交矩阵，是由卫星本体坐标到卫星轨道坐标系的旋转正交矩阵，是由相机坐标系到卫星本体坐标系的旋转正交矩阵。所述步骤201)中，采用特征提取算子从参考影像中提取待匹配特征点。所述步骤203)具体为：设参考影像上的点g1(x1,y1)和卫星影像上的点g2(x2,y2)之间存在线性的辐射畸变和仿射几何畸变，则有g1(x1,y1)+n1＝h0+h1g2(x2,y2)+n2＝h0+h1g2(a0+a1x1+a2y1,b0+b1x1+b2y1)+n2其中，(h0,h1)为辐射畸变参数，(n1,n2)为影像所包含的随机噪声，(a0,a1,a2,b0,b1,b2)为几何畸变参数，将上式线性化，并根据像点匹配初值，迭代计算几何形变和辐射畸变参数，得到最佳匹配点位。所述步骤3)具体为：301)通过严格成像模型，将密集控制点后向投影至影像像方，迭代计算其像方像素坐标；302)提取像方残差中周期性的残差分量(Δx′,Δy′)，并利用正弦函数拟合所述残差分量：式中，分别表示拟合得到的沿轨方向的周期性像方残差的幅值、频率和相位，分别表示拟合得到的垂轨方向的周期性像方残差的幅值、频率和相位；303)通过以下公式获得姿态颤振变化幅值：式中，Aroll为翻滚角姿态颤振变化幅值，Apitch为俯仰角姿态颤振变化幅值，n为线阵CCD探元的数量，persize为线阵CCD探元的大小，f为相机的焦距。与现有技术相比，本专利技术具有以下优点：(1)本专利技术利用外部参考数据(数字正射影像DOM和数字高程模型DEM)的支持，能够准确地探测出卫星平台的姿态颤振现象，从而还原卫星在轨运动的真实状态。(2)本专利技术在对参考影像进行特征提取前对其进行预处理，减少影像中的高频成分，使得与卫星影像的特性较为接近，进一步提高了卫星平台姿态颤振探测的准确性。(3)本专利技术在获取密集控制点时，采用相关系数阈值法和基于严格成像模型和仿射变换模型的几何约束条件的剔除策略进行误匹配点剔除，提高了方法的精度和卫星平台姿态颤振探测的准确性。附图说明图1为本专利技术的流程示意图；图2为国产测绘卫星立体影像；其中，(2a)为前视影像，(2b)为后视影像，(2c)为下视影像；图3为图2中影像的密集控制点分布示意图；其中，(3a)为前视影像的密集控制点分布，(3b)为后视影像的密集控制点分布，(3c)为下视影像的密集控制点分布；图4为下视影像的像方残差示意图；其中，(4a)为下视第2片CCD像方残差，(4b)为下视第3片CCD像方残差；图5为前视影像的像方残差示意图；其中，(5a)为前视第3片CCD像方残差，(5b)为前视第4片CCD像方残差；图6为后视影像的像方残差示意图；其中，(6a)为后视第3片CCD像方残差，(6b)为后视第4片CCD像方残差。具体实施方式下面结合附图和具体实施例对本专利技术进行详细说明。本实施例以本专利技术技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本专利技术的保护范围不限于下述的实施例。如图1所示，本实施例提供一种基于密集控制点的高分辨率卫星姿态颤振探测方法，其特征在于，包括以下步骤：步骤S1：参考影像预处理。本实施例使用的参考影像为DOM影像，由于DOM参考影像分辨率高、图像细节信息丰富，在提取特征点之前，首先要对DOM影像进行预处理即低通滤波，尽量降低影像的细节层次，使得与卫星影像的特性较为接近。对滤波后的图像进行降采样，使其和卫星影像的空间分辨率保持一致。高斯滤波用正态分布计算图像中每个像素的变换。在二维空间定义：其中，σ是正态分布的标准偏差。由这个公式生成的曲面的等高线是一组从中心开始呈正态分布的同心圆。因此，中心即原始像素的值具有最大的权重，相邻像素距离中心像素越远权重就越小。这样得到的图像更高地保留了边缘效果。步骤S2：预处理后的参考影像和实时获取的卫星影像(三线阵影像)的密集匹配，具体为：201)采用特征提取算子从预处理后的参考影像中提取出大量密集的待匹配特征点，并根据相应的DEM获得待匹配特征点的地面坐标。202)将所述待匹配特征点通过所述卫星影像的严格成像模型向后投影至卫星影像上，获得匹配点的初始值。严格成像模型为：式中，(X,Y,Z)为地面点在物方空间坐标系中的坐标，(XS,YS,ZS)为卫星在在物方空间坐标系中的位置，λ为比例因子，(x,y)为像方空间坐标，f为相机焦距，是由地球惯性系到WGS84坐标系的旋转正交矩阵，是由卫星轨道坐标系到地球惯性系的旋转正交矩阵，是由卫星本体坐标到卫星轨道坐标系的旋转正交矩阵，是由相机坐标系到卫星本体坐标系的旋转正交矩阵。203)根据匹配点初始值使用最小二乘匹配的方法，得到卫星影像和预处理之后的参考影像之间的密集控制点，具体为：设参考影像上的点g1(x1,y1)和卫星影像上的点g2(x2,y2)之间存在线性的辐射畸变和仿射几何畸变，则有g1(x1,y1)+n1＝h0+h1g2(x2,y2)+n2(3)＝h0+h1g2(a0+a1x1+a2y1,b0+b1x1+b2y1)+n2其中，(h0,h1)为辐射畸变参数，(n1,n2)为影像所包含的随机本文档来自技高网...

【技术保护点】
一种基于密集控制点的高分辨率卫星姿态颤振探测方法，其特征在于，包括以下步骤：1)参考影像预处理；2)预处理后的参考影像和实时获取的卫星影像的密集匹配；3)卫星姿态颤振的解算。

【技术特征摘要】
1.一种基于密集控制点的高分辨率卫星姿态颤振探测方法，其特征在于，包括以下步骤：1)参考影像预处理；2)预处理后的参考影像和实时获取的卫星影像的密集匹配，具体为：201)从预处理后的参考影像中提取出待匹配特征点，并根据相应的数字高程模型获得待匹配特征点的地面坐标，所述待匹配特征点的采点间隔小于5*5像素区域；202)将所述待匹配特征点通过所述卫星影像的严格成像模型向后投影至卫星影像上，获得匹配点的初始值；203)根据匹配点初始值使用最小二乘匹配的方法，得到卫星影像和预处理之后的参考影像之间的密集控制点；3)卫星姿态颤振的解算，具体为：301)通过严格成像模型，将密集控制点后向投影至影像像方，迭代计算其像方像素坐标；302)提取像方残差中周期性的残差分量(Δx′,Δy′)，并利用正弦函数拟合所述残差分量：式中，Axjitter、ωx、分别表示拟合得到的沿轨方向的周期性像方残差的幅值、频率和相位，Ayjitter、ωy、分别表示拟合得到的垂轨方向的周期性像方残差的幅值、频率和相位；303)通过以下公式获得姿态颤振变化幅值：式中，Aroll为翻滚角姿态颤振变化幅值，Apitch为俯仰角姿态颤振变化幅值，n为线阵CCD探元的数量，persize为线阵CCD探元的大小，f为相机的焦距。2.根据权利要求1所述的基于密集控制点的高分辨率卫星姿态颤振探测方法，其特征在于，步骤1)中，所述预处理包括高斯滤波和降采样处理。3.根据权利要求1所述的基于密集控制点的高分辨...

【专利技术属性】
技术研发人员：童小华，王凤香，刘世杰，李凌云，叶真，谢欢，陈鹏，张松林，
申请(专利权)人：同济大学，
类型：发明
国别省市：上海;31