一种无地面控制点的遥感影像定位方法技术

本发明专利技术公开了一种无地面控制点的遥感影像定位方法，其整体步骤为：采用外定标获取定轨定姿数据，即用于描述摄影中心的空间坐标值和像片面的空间姿态；构建严密几何模型以及二次多项式外方位元素模型；前方交会获取连接点地面坐标值；后方交会结合自适应权值结算外方位元素；拟合RPC参数；精化外方位和RPC参数。本发明专利技术在没有地面控制点条件下，仅利用卫星下传数据即可满足目标定位精度要求的目的，具有操作更加简便、工作量明显降低的优点。

A Remote Sensing Image Location Method without Ground Control Points

The invention discloses a remote sensing image positioning method without ground control points, whose overall steps are: obtaining orbit and attitude data by external calibration, i.e., describing the space coordinate value of the photographic center and the space attitude of the image plane; constructing a rigorous geometric model and a quadratic polynomial exterior orientation element model; obtaining the ground coordinate value of the connecting point by forward intersection; and combining the rear intersection with self-adaptation. External orientation elements should be settled by weight; RPC parameters should be fitted; exterior orientation and RPC parameters should be refined. In the absence of ground control points, the method can satisfy the target positioning accuracy requirement only by downloading data from satellite, and has the advantages of simpler operation and lower workload.

【技术实现步骤摘要】
一种无地面控制点的遥感影像定位方法
本专利技术涉及一种定位方法，尤其涉及一种无地面控制点的遥感影像定位方法。
技术介绍
为了利用遥感方法识别地面目标，并确定其地理空间坐标，需要将获取到的遥感影像进行正射纠正。通常，遥感影像正射纠正需要依赖已知大地坐标的地面控制点(GroudControlPoints，简称GCP)。地面控制点是目前对航空像片和卫星遥感影像进行各种几何纠正和地理定位的重要数据源，其数量、质量和分布等指标对影像纠正的精确性和可靠性都具有直接影响。然而，地面控制点的选取会占用遥感影像处理的大量工作量，而且获取足够数量的地面控制点通常是十分困难的，能够替代野外测量的大比例尺地形图也很难得到，因此，研究无地面控制点条件下的卫星遥感影像的高精度几何处理就成为一项紧迫的任务。
技术实现思路
为了解决上述技术所存在的不足之处，本专利技术提供了一种无地面控制点的遥感影像定位方法。为了解决以上技术问题，本专利技术采用的技术方案是：一种无地面控制点的遥感影像定位方法，其整体步骤为：S1、采用外定标获取定轨定姿数据，即用于描述摄影中心的空间坐标值和像片面的空间姿态(XSi,YSi,ZSi,ωi,κi)，其中i表示第i扫描行；通过像点坐标量测获取线阵影像连接点坐标值(x,y)；S2、模型构建Ⅰ、构建严密几何模型通过推扫式成像的线阵CCD传感器获取连续的影像条带；每一扫描行影像与被摄物体之间具有严格的中心投影关系，并且都有各自的外方位元素；以扫描行方向为x、飞行方向为y建立瞬时像平面坐标系，设第i扫描行的外方位元素为XSi,YSi,ZSi,ωi,κi，则瞬时构像方程式为：其中，a1,a2,a3,b1,b2,b3,c1,c2,c3为像点坐标对应的姿态角元素计算的旋转矩阵元素；f为摄影中心到像片的垂距；x0,y0为像主点；b1＝cosωisinκib2＝cosωicosκib3＝-sinωiⅡ、构建二次多项式外方位元素模型线阵CCD影像在同一摄影时刻具有相同的外方位元素，由于星载CCD传感器所在的航天器受外界阻力小，飞行轨道平稳，姿态变化小，因此在某一范围内近似认为外方位元素是时间的低阶函数；采用如下数学模型对摄影时刻t的角元素加以表示：其中，φ0,φ1,φ2,ω0,ω1,ω2,k0,k1,k2为二次多项式系数；利用行时信息以及外定标获取的定姿数据，根据公式2最小二乘法拟合二次多项式系数，构建出二次多项式模型，为后续提供初始姿态角元素；S3、前方交会根据构建的模型，利用连接点坐标值(x,y)，读取对应的行时信息，结合已经获取的二次多项式外方位元素模型，求解连接点对应的姿态角元素ωi,κi，联合定轨参数XSi,YSi,ZSi，组成六个外方位元素；将瞬时构像方程式进行整理，得：式中：l1＝fa1+(x-x0)a3,l2＝fb1+(x-x0)b3,l3＝fc1+(x-x0)c3lx＝fa1XSi+fb1YSi+fc1ZSi+(x-x0)a3XSi+(x-x0)b3YSi+(x-x0)c3ZSil4＝fa2+(y-y0)a3,l5＝fb2+(y-y0)b3,l6＝fc2+(y-y0)c3ly＝fa2XSi+fb2YSi+fc2ZSi+(y-y0)a3XSi+(y-y0)b3YSi+(y-y0)c3ZSi其中，a1,a2,a3,b1,b2,b3,c1,c2,c3为根据公式2结合连接点坐标对应的行时信息求解初始姿态角元素ωi,κi构建的旋转矩阵元素；若n幅影像含有同一空间点，则总共有n个如公式3的方程组求解地面坐标值(X,Y,Z)；根据最小二乘原理，求解连接点地面坐标值；S4、后方交会将前方交会获取的所有连接点地面坐标值(X,Y,Z)视为真值，保持不变，而把相应的像点坐标值(x,y)视为观测值，加入相应的改正数vx,vy，代入共线方程式公式1，并结合二次多项式外方位元素模型公式2，泰勒展开取一次项，使之线性化得到误差方程式的一般形式：用矩阵形式表示为：V＝AX-l公式5式中：v＝[vxvy]Tl＝[lxly]TA矩阵为公式1对二次多项式外方位元素模型公式2系数求偏导；其中，每一对连接点的像方坐标值均可列出如公式4的方程组；根据最小二乘间接平差原理，可构成法方程式，为：ATPAX＝ATPL公式6式中，P为观测值的权矩阵；x方向权值通过连接点像方x方向的残差值vx与所有残差vx的中误差值σ0确定，原理为：残差大权值小，残差小权值大；同理，y方向权值采用残差值vy和同上述原理确定；式中k0＝σ0；k1＝3σ0，设初始权矩阵P＝E，E为单位矩阵，v为各像点的残差，由后方交会解算得到并作为下一次求权的依据，以达到剔除粗差的目的；根据公式6式，求得：X＝(ATPA)-1ATPL公式8根据公式8，利用逐渐迭代趋近的方法，解算改正数X；结合公式2中的模型，获取修正后的外方位姿态角元素。本专利技术在没有地面控制点条件下，仅利用卫星下传数据即可满足目标定位精度要求的目的，具有操作更加简便、工作量明显降低的优点。附图说明图1为本专利技术的整体流程示意图。具体实施方式下面结合附图和具体实施方式对本专利技术作进一步详细的说明。实施例一、图1所示的一种无地面控制点的遥感影像定位方法，其整体流程为：S1、采用外定标(即热红外波段定标法)获取定轨定姿数据，即获取用于描述摄影中心的空间坐标值和像片面的空间姿态(XSi,YSi,ZSi,ωi,κi)，其中i表示第i扫描行；通过像点坐标量测获取线阵影像连接点坐标值(x,y)。本实施例以WGS-84坐标系为坐标基准。S2、模型构建Ⅰ、构建严密几何模型线阵CCD传感器采用推扫式成像，获取连续的影像条带。每一扫描行影像与被摄物体之间具有严格的中心投影关系，并且都有各自的外方位元素。以扫描行方向为x方向，飞行方向为y方向，建立瞬时像平面坐标系，设第i扫描行的外方位元素为XSi,YSi,ZSi,ωi,κi，则瞬时构像方程式为：其中，a1,a2,a3,b1,b2,b3,c1,c2,c3为像点坐标对应的姿态角元素计算的旋转矩阵元素；f为摄影中心到像片的垂距；x0,y0为像主点。b1＝cosωisinκib2＝cosωicosκib3＝-sinωiⅡ、构建二次多项式外方位元素模型线阵CCD影像在同一摄影时刻具有相同的外方位元素，由于星载CCD传感器所在的航天器受外界阻力小，飞行轨道平稳，姿态变化小，因此在某一范围内，可以近似认为外方位元素是时间的低阶函数。可采用如下数学模型对摄影时刻t的角元素加以表示。其中，φ0,φ1,φ2,ω0,ω1,ω2,k0,k1,k2为二次多项式系数。本专利技术利用行时信息以及外定标获取的定姿数据，根据公式2最小二乘法拟合二次多项式系数，构建出二次多项式模型，为后续提供初始姿态角元素，无需对外定标获取的定轨数据进行修正。S3、前方交会根据构建的模型，利用连接点坐标值(x,y)，读取对应的行时信息，结合已经获取的二次多项式外方位元素模型(公式2)，求解连接点对应的姿态角元素ωi,κi，联合定轨参数XSi,YSi,ZSi，组成六个外方位元素。将瞬时构像方程式(公式1)进行整理，得：式中：l1＝fa1+(x-x0)a3,l2＝fb1+(x-x0)b3,l3＝fc1+(x-x0)c3lx＝fa1XSi+fb1YSi+fc1ZSi+(本文档来自技高网...

【技术保护点】
1.一种无地面控制点的遥感影像定位方法，其特征在于：所述方法的整体步骤为：S1、采用外定标获取定轨定姿数据，即用于描述摄影中心的空间坐标值和像片面的空间姿态

【技术特征摘要】
1.一种无地面控制点的遥感影像定位方法，其特征在于：所述方法的整体步骤为：S1、采用外定标获取定轨定姿数据，即用于描述摄影中心的空间坐标值和像片面的空间姿态其中i表示第i扫描行；通过像点坐标量测获取线阵影像连接点坐标值(x,y)；S2、模型构建Ⅰ、构建严密几何模型通过推扫式成像的线阵CCD传感器获取连续的影像条带；每一扫描行影像与被摄物体之间具有严格的中心投影关系，并且都有各自的外方位元素；以扫描行方向为x、飞行方向为y建立瞬时像平面坐标系，设第i扫描行的外方位元素为则瞬时构像方程式为：其中，a1,a2,a3,b1,b2,b3,c1,c2,c3为像点坐标对应的姿态角元素计算的旋转矩阵元素；f为摄影中心到像片的垂距；x0,y0为像主点；b1＝cosωisinκib2＝cosωicosκib3＝-sinωiⅡ、构建二次多项式外方位元素模型线阵CCD影像在同一摄影时刻具有相同的外方位元素，由于星载CCD传感器所在的航天器受外界阻力小，飞行轨道平稳，姿态变化小，因此在某一范围内近似认为外方位元素是时间的低阶函数；采用如下数学模型对摄影时刻t的角元素加以表示：其中，φ0,φ1,φ2,ω0,ω1,ω2,k0,k1,k2为二次多项式系数；利用行时信息以及外定标获取的定姿数据，根据公式2最小二乘法拟合二次多项式系数，构建出二次多项式模型，为后续提供初始姿态角元素；S3、前方交会根据构建的模型，利用连接点坐标值(x,y)，读取对应的行时信息，结合已经获取的二次多项式外方位元素模型，求解连接点对应的姿态角元素联合定轨参数XSi,YSi,ZSi，组成六个外方位元素；将瞬时构像方程式进行整理，得：式中：l1＝fa1+(x-x0)a3,l2＝fb1+(x-x0)b3,l3＝fc1+(x-x0)c3lx＝fa1XSi+fb1YSi+fc1ZSi+(x-x...

【专利技术属性】
技术研发人员：吴方才，靳笑琳，张丽，牛东，
申请(专利权)人：航天星图科技北京有限公司，
类型：发明
国别省市：北京,11