CN106403902B - 一种星地协同的光学卫星在轨实时几何定位方法及系统 - Google Patents
一种星地协同的光学卫星在轨实时几何定位方法及系统 Download PDFInfo
- Publication number
- CN106403902B CN106403902B CN201610786005.0A CN201610786005A CN106403902B CN 106403902 B CN106403902 B CN 106403902B CN 201610786005 A CN201610786005 A CN 201610786005A CN 106403902 B CN106403902 B CN 106403902B
- Authority
- CN
- China
- Prior art keywords
- calibration
- satellite
- model
- positioning
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000003287 optical effect Effects 0.000 title claims abstract description 44
- 238000012545 processing Methods 0.000 claims abstract description 41
- 238000003384 imaging method Methods 0.000 claims abstract description 37
- 238000013461 design Methods 0.000 claims abstract description 11
- 239000000523 sample Substances 0.000 claims description 29
- 239000011159 matrix material Substances 0.000 claims description 24
- 238000009434 installation Methods 0.000 claims description 18
- 238000012937 correction Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000002347 injection Methods 0.000 claims description 8
- 239000007924 injection Substances 0.000 claims description 8
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 6
- 239000000243 solution Substances 0.000 claims description 6
- 238000012795 verification Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000007711 solidification Methods 0.000 claims description 3
- 230000008023 solidification Effects 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 abstract 1
- 230000005540 biological transmission Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000003672 processing method Methods 0.000 description 4
- 238000001444 catalytic combustion detection Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000013144 data compression Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012634 optical imaging Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种星地协同的光学卫星在轨实时几何定位方法及系统，包括构建适用于星上实时处理单元的光学卫星成像定位模型，由地面实验室对相机内部及平台安装关系的检校参数或设计参数获取星上定位模型参数的初始值，卫星在轨后获取地面定标场影像数据并下传，利用下传的定标场影像在地面处理系统中完成定标匹配及定标参数解算；定标结果验证与模型参数上注更新。本发明针对星上处理环境的局限性，结合地面处理系统，形成星地协同的处理模式，实现了光学遥感卫星在轨高精度实时几何定位，提高了对地观测系统的应用效能和时效性，为实现智能化、高效的遥感数据星上实时处理技术提供必要的基础。
Description
技术领域
本发明的技术方案属于卫星遥感数据处理技术领域，特别涉及一种星地协同的光学卫星星上几何定位处理的方法及系统。
背景技术
随着影像分辨率的提高，光学遥感卫星获取的数据量呈几何级数增长，远远超出了数据压缩传输能力的发展，导致星上实时获取的数据无法及时下传。以高分二号遥感卫星为例，该卫星原始获取数据率达到7Gb/s，若采用常规压缩传输方法，通过2*450Mbps数传链路，也无法完成全部获取数据的实时下传，严重制约了用户获取遥感数据和信息的时效性。因此，面向高分辨率光学成像所获取的实时海量遥感数据，现有的星地数据处理模式、数据压缩方法、数据处理方法无法满足各类用户及时准确获取信息的要求。从而迫切需要研究新的数据处理模式、自动化和高时效性的数据处理方法，提升海量遥感影像的实时数据处理能力和信息提取水平，充分发挥对地观测系统的应用效能。
光学遥感卫星星上在轨实时处理打破传统“星上成像-影像下传-地面处理”的数据处理模式，可在星上从海量遥感影像数据中，针对兴趣任务目标进行提取，实时对目标影像块进行数据处理，进而转化为有效信息快速分发至地面用户，以任务驱动角度出发，大大提高了遥感数据应用的时效性及自动化、智能化程度。其中，具有高精度、实时性的星上几何定位技术是实现以任务驱动的光学遥感卫星星上在轨实时处理的必要环节，精确、可靠的地理位置信息是卫星获取的空间数据得以实时提取、转化为有效信息的基础。
光学卫星高精度几何定位依赖于精确的成像几何模型参数，受卫星发射时应力释放、卫星运行时空间热环境及力学环境等因素的影响，地面实验室对于光学相机内部及平台安装参数的检校无法满足卫星定位的精度需求。目前，地面系统处理一般采用基于地面定标场的几何定标方法，以卫星在轨运行时获取的几何定标场影像数据进行模型参数的精化解算，而星上几何定位技术受制于存储环境及处理环境的约束，无法在星上进行几何定标。
发明内容
本发明针对光学卫星成像高精度在轨几何定位问题，提出了一种星地协同的光学卫星在轨实时几何定位方法及系统。
本发明提供的技术方案为一种星地协同的光学卫星在轨实时几何定位方法，包括以下步骤：
步骤1，定位模型构建及算法固化，构建适用于星上实时处理单元的光学卫星成像定位模型，并将相应定位求解方法固化于星上硬件环境，保留模型参数更新上注接口；
所述光学卫星成像定位模型，是采用线阵CCD探元指向角的内定向模型，建立基于严密共线方程模型与地球椭球面模型相交的在轨定位模型；
步骤2，初值确定，由地面实验室对相机内部及平台安装关系的检校参数或设计参数获取星上定位模型参数的初始值；
步骤3，定标数据获取，卫星在轨运行后对地面定标场进行成像，获取适宜几何定标的影像数据并下传至地面系统；
步骤4，地面系统几何定标，包括在地面处理系统中完成定标控制点的密集匹配、定标参数解算；
步骤5，定标结果验证与模型参数上注更新，对定标精度评价，确定定标结果的正确性后，更新星上定位的相应参数。
而且，步骤1中，所述线阵CCD探元指向角的内定向模型如下，
其中，(Vimage)cam为像元在像空间坐标系下的指向矢量，x、y分别为像元在像平面坐标系下垂轨和沿轨方向的坐标，ψx(s)、ψy(s)为探元s指向矢量在沿轨和垂轨方向的角度分量，f为相机主距；
对于具有多片CCD拼接成像的相机，设有m片CCD，则对于m片CCD分别采用多组三次多项式(ψxj(s)ψyj(s))进行描述，
其中，s为探元号，ψxj(s)、ψyj(s)为各片上探元指向角在沿轨和垂轨方向的角度分量，j表示CCD的标号，(ax0j,ax1j,ax2j,ax3j,ay0j,ay1j,ay2j,ay3j)为内定标参数。
而且，步骤1中，所述用基于严密共线方程模型与地球椭球面模型相交的在轨定位模型如下，
其中，为卫星本体坐标系到相机坐标系的旋转矩阵，为J2000坐标系到本体坐标系的旋转矩阵，为WGS84坐标系到J2000坐标系的旋转矩阵，(Xs,Ys,Zs)为成像时刻的卫星在WGS84坐标系下的位置，(X,Y,Z)为目标点在WGS84坐标系下的坐标，λ为比例因子，aWGS84和bWGS84分别为WGS84椭球的长半轴和短半轴，h为目标定位点在物方的高程。
而且，步骤1中，所述定位求解方法采用基于DEM数据的高程迭代定位求解方法，实现方式如下，
以目标点高程初值h0＝0，在DEM数据的支持下，进行高程迭代求解，包括执行以下步骤，
a.令i＝1，目标点高程h＝h0，即令h＝0，代入椭球面模型；
b.共线方程与椭球面方程联立，求得目标点物方坐标，获取光线与椭球高h处的交点Mi；
c.若i＞1，则判断此次获得的交点Mi与上次计算坐标Mi-1的修正量d(Mi-1,Mi)是否小于阈值d；
d.若修正量小于阈值，则输出定位结果，若i＝1或修正量大于阈值，则由目标点物方坐标Mi在DEM上内插更新高程值h＝h(Mi)，令i＝i+1，返回步骤b，重复b、c、d步骤直至收敛。
而且，步骤2中，确定星上模型参数的初值，实现方式如下，
对于相机在平台上安装关系的实验室检校，获取三安装角参数，与构成的三个旋转角一致，直接作为初值；
对于相机内部参数的实验室检校，按照严格物理模型，测量相机主距f、各片CCD首像元在相机坐标系下的坐标(x0j,y0j)，设定内定向参数的初始值，如下式，
其中，pixelsize为CCD像元大小设计值。
而且，步骤4中构建结合内部探元指向角模型与外部安装矩阵补偿的在轨几何定标模型，实现方法如下：
作为外定标参数，用于恢复相机坐标系在空间中的位置和姿态；(ax0j,ax1j,ax2j,ax3j,ay0j,ay1j,ay2j,ay3j)(j＝1,2,...,m)作为内定标参数，用于确定相机内部CCD各探元在相机坐标系下的坐标。
而且，步骤4中，基于步骤1所得在轨定位模型，执行以下步骤，
a.设在待定标影像上量测了K个高精度地面控制点作为定向点，控制点的WGS84地心直角坐标为(Xi Yi Zi)，像点坐标为(si li)，i＝1,2,3...k；
b.令：
其中，为像点光线在本体坐标系下的矢量，(pitch,roll,yaw)为相机在本体坐标系上的三个安装偏置角，为本体坐标系到相机坐标系的旋转矩阵
设外定标参数XE、内定标参数XI为自变量，令F()、G()分别为像空间坐标系下像点沿轨及垂轨方向的矢量残差函数，则有：
c.对外定标参数XE、内定标参数XI赋初值
d.将当前的内定标参数XI视为“真值”，将外定标参数XE视为待求的未知参数，将相应的当前值代入矢量残差函数，对每个定向点，进行线性化处理，建立误差方程式，利用最小二乘平差计算更新外定标参数XE的当前值，
e.重复步骤d，迭代计算直至外定标参数改正数均小于预设阈值时停止，进入步骤f；
f.将当前的外定标参数XE视为“真值”，将内定标参数XI视为待求的未知参数，将相应的当前值代入矢量残差函数，对每个定向点，进行线性化处理，建立误差方程式，利用最小二乘平差计算更新内定标参数XI的当前值，
g.重复步骤f，迭代计算直至内定标参数改正数均小于预设阈值时停止，完成几何内外定标参数求解。
本发明相应提供一种星地协同的光学卫星在轨实时几何定位系统，包括以下模块：
定位模型构建及算法固化模块，用于构建适用于星上实时处理单元的光学卫星成像定位模型，并将相应定位求解方法固化于星上硬件环境，保留模型参数更新上注接口；
所述光学卫星成像定位模型，是采用线阵CCD探元指向角的内定向模型，建立基于严密共线方程模型与地球椭球面模型相交的在轨定位模型；
初值确定模块，用于由地面实验室对相机内部及平台安装关系的检校参数或设计参数获取星上定位模型参数的初始值；
定标数据获取模块，用于卫星在轨运行后对地面定标场进行成像，获取适宜几何定标的影像数据并下传至地面系统；
地面系统几何定标模块，用于在地面处理系统中完成定标控制点的密集匹配、定标参数解算；
定标结果验证与模型参数上注更新模块，用于对定标精度评价，确定定标结果的正确性后，更新星上定位的相应参数。
本发明提供一种以地面几何定标和星上实时几何定位协同处理的星上高精度实时几何定位技术方案，满足了星上高精度实时几何定位的需求，解决了光学卫星在轨实时处理的一个关键技术问题。该技术方案针对星上处理环境的局限性，结合地面处理系统，形成星地协同的处理模式，实现了光学遥感卫星在轨高精度实时几何定位，提高了对地观测系统的应用效能和时效性，为实现智能化、高效的遥感数据星上实时处理技术提供必要的基础，具有重要的市场价值。
发明附图
图1为本发明的光学卫星在轨几何定位方法流程图。
图2为本发明的探元指向角模型示意图。
图3为本发明的基于DEM的单点定位迭代计算流程图。
图4为本发明的基于定标场影像的地面系统在轨几何定标流程图。
具体实施方式
以下结合附图和实施例详细说明本发明技术方案。
本发明技术方案可采用计算机软件方式支持自动运行流程。本发明实施例星地协同的光学卫星在轨实时几何定位方法如图1所示，分以下步骤进行详细描述。
(1)定位模型构建及算法固化：构建适用于星上实时处理单元的光学卫星成像定位模型，并将算法固化于星上硬件环境，保留模型参数更新上注接口。
进一步地，步骤(1)中考虑成像过程中各种误差的几何特性、统计特性及变形规律，建立了优化的光学卫星星上在轨成像定位模型。其中采用探元指向角模型代替相机的严格内检校模型，避免模型的过度参数化，并可在同等精度水平下消除相机镜头畸变、CCD畸变及内方位元素标定误差；相机外部误差通过由相机与平台间安装角度构成的旋转矩阵进行补偿。将算法固化于星上硬件时，保留探元指向角模型参数及安装角参数的更新接口。
对于光学卫星相机，可以将影响几何定位的因素分为两类，一类是内部误差，包括镜头畸变、线阵CCD变形因素导致的相机内部光轴指向的改变，另一类是外部误差，包括相机的安装误差、热变形导致的相机安装相对关系的改变和外方位元素的观测误差。其中镜头畸变、线阵CCD的变形、相机的安装误差和热变形属于静态误差，具有很强的系统性，可以通过在轨几何标定的方式进行标校和补偿。
相机严格内检校模型中由于包含众多的物理畸变参数，某些参数之间具有强相关性，因此存在过度参数化的问题，难以分别精确标校，并不适合作为相机的内定向模型。因此设计CCD探元指向角模型，如图2，X1、Y1、Z1为像空间坐标系的三轴，O1为投影中心，Vimage为像元在像空间坐标系下的指向矢量，ψx、ψy为Vimage分别在沿轨和垂轨方向的角度分量。对各探元的指向角进行描述，其本质上就是确定各探元的像方矢量在单位主距下相机焦平面上的投影平面坐标，即对相机主距进行了归一化处理，采用线阵CCD探元指向角的内定向模型，如下式：
其中，(Vimage)cam为像元在像空间坐标系下的指向矢量，x、y分别为像元在像平面坐标系下垂轨和沿轨方向的坐标，ψx(s)、ψy(s)为探元s指向矢量在沿轨和垂轨方向的角度分量，f为相机主距。
由于相机严格物理模型本质上就是一个三次多项式模型，为了确定各探元的指向角，采用一个三次多项式对相机CCD上各探元在相机坐标系下的指向角进行拟合，作为相机的内定向模型，对于具有多片CCD拼接成像的相机，设有m片CCD，则对于m片CCD分别采用多组三次多项式(ψxj(s)ψyj(s))进行描述：
其中，s为探元号，ψxj(s)、ψyj(s)为各片上探元指向角在沿轨和垂轨方向的角度分量，j表示CCD的标号，(ax0j,ax1j,ax2j,ax3j,ay0j,ay1j,ay2j,ay3j)为内定标参数。
则代入光学卫星成像严密共线方程模型，与地球椭球面模型联立，利用像点光束与地面物方高程面相交的几何关系，建立的光学成像在轨几何定位模型如下：
其中，为卫星本体坐标系到相机坐标系的旋转矩阵，由相机在平台上的安装角求得；为J2000坐标系到本体坐标系的旋转矩阵，由成像时刻的卫星平台三姿态角求得；为WGS84坐标系到J2000坐标系的旋转矩阵，由成像时刻及国际地球自转服务(IERS)发布的地球自转参数求得；(Xs,Ys,Zs)为成像时刻的卫星在WGS84坐标系下的位置；(X,Y,Z)为目标点在WGS84坐标系下的坐标；λ为比例因子；aWGS84和bWGS84分别为WGS84椭球的长半轴和短半轴；h为目标定位点在物方的高程。
基于构建的定位模型，在星上定位解算过程中，首先组成共线方程：输入目标点的像方坐标，则可获取其在相机坐标系下的矢量；根据像点的成像时刻，可由GPS轨道测量星历及星敏姿态测量星历内插出对应的卫星位置(WGS84坐标系)及姿态(J2000坐标系)；根据像点的成像时刻及地球自转参数可求得WGS84坐标系与J2000坐标系的转换关系；由相机安装角计算本体坐标系到相机坐标系的旋转矩阵。然后，将共线方程与地球椭球面方程联立，以目标点高程初值h0＝0，在DEM数据的支持下，进行高程迭代求解，如图3所示，具体步骤为：
a.令i＝1，目标点高程h＝h0，即令h＝0，代入椭球面模型；
b.共线方程与椭球面方程联立，求得目标点物方坐标，即光线与椭球高h处的交点Mi；
c.若i＞1，则判断此次获得的交点Mi与上次计算坐标Mi-1的修正量d(Mi-1,Mi)是否小于阈值d；
d.若修正量小于阈值，则输出定位结果，若i＝1或修正量大于阈值，则由目标点物方坐标Mi在DEM上内插更新高程值h＝h(Mi)，令i＝i+1，返回步骤b，重复b、c、d步骤直至收敛。
在卫星地面建设阶段，完成该定位模型及算法在星上处理模块(如DSP)上的固化以及DEM数据在星上存储模块的注入，并将模型中存在系统误差、需在轨几何标定的参数(ax0j,ax1j,ax2j,ax3j,ay0j,ay1j,ay2j,ay3j)(j＝1,2,...,m)保留上注更新的接口。
(2)初值确定：由地面实验室对相机内部及平台安装关系的检校参数或设计参数获取星上定位模型参数初始值。
卫星发射前，采用实验室检校的手段可获取相机内部参数及安装矩阵的地面检校值，作为在轨精确标定前的模型初值。
对于相机在平台上安装关系相应卫星本体坐标系到相机坐标系的旋转矩阵的实验室检校，一般获取三安装角参数，与构成的三个旋转角一致，直接作为初值；
对于相机内部参数的实验室检校，一般按照严格物理模型，测量其相机主距f、各片CCD首像元在相机坐标系下的坐标(x0j,y0j)(j＝1,2,...,m)，此时可忽略相机高阶畸变，结合CCD像元大小设计值pixelsize，设定内定向参数的初始值，如下式。
(3)定标数据获取：卫星在轨运行后对地面定标场进行成像，获取适宜几何定标的影像数据并下传至地面系统。
进一步地，步骤(3)中定标影像数据应该选择地物清晰无云、成像角度较小、控制点在影像内均匀分布的影像。
卫星发射后在轨运行过程中，规划成像任务时将定标场目标列为优先，对定标场进行较小成像角度的拍摄(星下点成像为最优选方案)，并将成像数据下传至地面处理系统，选取成像天气晴朗无云、地面控制数据覆盖均匀的影像数据作为待定标影像数据。
(4)地面系统几何定标：包括在地面处理系统中完成定标控制点的密集匹配、定标参数解算。
基于定标场影像的在轨几何定标，处理流程如图4所示，基于实验室标定参数、轨道和姿态数据(由GPS和星敏观测获得)和待定标影像数据进行。
在轨几何定标需要密集控制点匹配：获取待定标影像数据后，需一定数量、在影像上均匀分布的控制点作为控制信息。目前国内卫星地面几何定标场均提供了高精度数字正射影像(DOM)和数字高程模型(DEM)等参考数据，利用影像高精度匹配技术，将待定标影像直接和定标场的DOM和DEM参考数据进行影像匹配，从而实现控制点的自动量测，获取大量同名像点，为后续的平差解算提供必要可靠的控制信息。
同步骤(1)，将采用探元指向角模型的严密共线方程作为在轨几何定标模型，作为外定标参数，用于精确恢复相机坐标系在本体坐标系下的安装偏置关系；(ax0j,ax1j,ax2j,ax3j,ay0j,ay1j,ay2j,ay3j)(j＝1,2,...,m)作为内定标参数，用于确定相机内部CCD各探元在相机坐标系下的坐标。
由控制点匹配结果，基于建立的内外几何定标模型，进行在轨几何内外定标参数求解，具体解算公式及流程如下：
a.假设在待定标影像上量测了K个高精度地面控制点作为定向点，控制点的WGS84地心直角坐标为(Xi Yi Zi)，像点坐标为(si li)，i＝1,2,3...k；
b.令：
其中，为像点光线在本体坐标系下的矢量，(pitch,roll,yaw)为相机在本体坐标系上的三个安装偏置角，为本体坐标系到相机坐标系的旋转矩阵。
设外定标参数XE、内定标参数XI为自变量，令F()、G()分别为像空间坐标系下像点沿轨及垂轨方向的矢量残差函数，则有：
c.对外定标参数XE、内定标参数XI赋初值这里初值为实验室检校值或初始设计值。
d.将当前内定标参数XI视为“真值”，将外定标参数XE视为待求的未知参数。将它们的当前值代入上式，对每个定向点，对其进行线性化处理，建立误差方程式：
Vi＝AiX-Li Pi
其中
式中，Li是利用内外定标参数当前值代入公式计算得到的常数向量；Ai是误差方程式的系数矩阵；X代表外定标参数改正数dXE；Pi是观测值的权；Vi是像点残差向量；(dpitch,droll,dyaw)是相机平台安装角改正数；Fi、Gi与步骤b中公式一致，为各像点在像空间坐标系下的矢量残差函数。
计算法方程系数矩阵，
其中，L是常数向量，A是系数矩阵，P为权值向量。
利用最小二乘平差计算X，如下，
X＝(ATPA)-1(ATPL)
更新外定标参数XE的当前值：
e.重复步骤d，迭代计算直至外定标参数改正数均小于阈值(本领域技术人员可自行预设，优选地取10-12)时停止，进入步骤f。
f.同样，将当前的外定标参数XE视为“真值”，将内定标参数XI视为待求的未知参数，将相应的当前值代入矢量残差函数，对每个定向点，进行线性化处理，建立误差方程式，利用最小二乘平差计算更新内定标参数XI的当前值，具体实现如下，
将外定标参数的当前值视为“真值”，内定标参数视为待求的未知参数，代入公式对每个定向点构建误差方程式：
Vi＝BiY-Li Pi
其中，
式中，Li是利用内外定标参数当前值计算得到的常数向量；Bi是误差方程式的系数矩阵；Y代表内定标参数改正数dXI；Pi是观测值的权；Vi是像点残差向量；(dax0,dax1,dax2,dax3,day0,day1,day2,day3)是相机内定标参数改正数；Fi、Gi与步骤b中公式一致，为各像点在像空间坐标系下的矢量残差函数。
计算法方程系数矩阵；
其中，L是常数向量，B是系数矩阵，P为权值向量。
利用最小二乘平差计算Y，如下式；
Y＝(BTPB)-1(BTPL)
更新内定标参数XI的当前值。
g.重复步骤f，迭代计算直至内定标参数改正数均小于阈值(本领域技术人员可自行预设，优选地取10-12)时停止，完成几何内外定标参数求解。
(5)定标结果验证与模型参数上注更新：对定标精度评价，确定定标结果的正确性后，更新星上定位算法的相应参数。
定标完成后，需要对定标效果进行评价验证，一般通过对定标后生产的产品数据进行产品几何精度测试以验证其正确性。
进一步地，步骤(5)中定标结果的正确性评价是指对定标后生产的产品数据进行产品几何精度的评价，利用地面检查点评价影像的绝对定位精度，对比定标前后的定位精度指标，评估参数精化效果。
针对星上单点定位，验证定标后产品的绝对几何精度即可，通过利用地面检查点参考信息，将检查点的地面坐标采用几何检校后的模型反算得到相应的图像坐标，求出沿轨道/垂直轨道方向真实的图像坐标与计算的图像坐标之间的差值，并统计其数学期望，由影像的几何分辨率可换算为物方精度。对比定标前后产品的绝对几何精度提升水平，评价定标结果参数的正确性和可用性。
最后，将地面处理系统在轨定标所获取的精确参数，通过星上保留的参数上注接口，提供给固化在星上硬件环境中的实时几何定位模型，更新不精确的实验室检校参数，从而实现星上高精度实时定位。
具体实施时，本发明所提供方法可基于软件技术实现自动运行流程，也可采用模块化方式实现相应系统。
本发明实施例相应提供一种星地协同的光学卫星在轨实时几何定位系统，包括以下模块：
定位模型构建及算法固化模块，用于构建适用于星上实时处理单元的光学卫星成像定位模型，并将相应定位求解方法固化于星上硬件环境，保留模型参数更新上注接口；
所述光学卫星成像定位模型，是采用线阵CCD探元指向角的内定向模型，建立基于严密共线方程模型与地球椭球面模型相交的在轨定位模型；
初值确定模块，用于由地面实验室对相机内部及平台安装关系的检校参数或设计参数获取星上定位模型参数的初始值；
定标数据获取模块，用于卫星在轨运行后对地面定标场进行成像，获取适宜几何定标的影像数据并下传至地面系统；
地面系统几何定标模块，用于在地面处理系统中完成定标控制点的密集匹配、定标参数解算；
定标结果验证与模型参数上注更新模块，用于对定标精度评价，确定定标结果的正确性后，更新星上定位的相应参数。
各模块具体实现可参见相应步骤，本发明不予赘述。
本文中所描述的具体实例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (8)
1.一种星地协同的光学卫星在轨实时几何定位方法，包括以下步骤：
步骤1，定位模型构建及算法固化，构建适用于星上实时处理单元的光学卫星成像定位模型，并将相应定位求解方法固化于星上硬件环境，保留模型参数更新上注接口；所述光学卫星成像定位模型，是采用线阵CCD探元指向角的内定向模型，建立基于严密共线方程模型与地球椭球面模型相交的在轨定位模型；
步骤2，初值确定，由地面实验室对相机内部及平台安装关系的检校参数或设计参数获取星上定位模型参数的初始值；
步骤3，定标数据获取，卫星在轨运行后对地面定标场进行成像，获取适宜几何定标的影像数据并下传至地面系统；
步骤4，地面系统几何定标，包括在地面处理系统中完成定标控制点的密集匹配、定标参数解算；
步骤5，定标结果验证与模型参数上注更新，对定标精度评价，确定定标结果的正确性后，更新星上定位的相应参数。
2.根据权利要求1所述星地协同的光学卫星在轨实时几何定位方法，其特征在于：步骤1中，所述线阵CCD探元指向角的内定向模型如下，
其中，(Vimage)cam为像元在像空间坐标系下的指向矢量，x、y分别为像元在像平面坐标系下垂轨和沿轨方向的坐标，ψx(s)、ψy(s)为探元s指向矢量在沿轨和垂轨方向的角度分量，f为相机主距；
对于具有多片CCD拼接成像的相机，设有m片CCD，则对于m片CCD分别采用多组三次多项式(ψxj(s) ψyj(s))进行描述，
其中，s为探元号，ψxj(s)、ψyj(s)为各片上探元指向角在沿轨和垂轨方向的角度分量，j表示CCD的标号，(ax0j,ax1j,ax2j,ax3j,ay0j,ay1j,ay2j,ay3j)为内定标参数。
3.根据权利要求2所述星地协同的光学卫星在轨实时几何定位方法，其特征在于：步骤1中，所述基于严密共线方程模型与地球椭球面模型相交的在轨定位模型如下，
其中，为卫星本体坐标系到相机坐标系的旋转矩阵，为J2000坐标系到本体坐标系的旋转矩阵，为WGS84坐标系到J2000坐标系的旋转矩阵，(Xs,Ys,Zs)为成像时刻的卫星在WGS84坐标系下的位置，(X,Y,Z)为目标点在WGS84坐标系下的坐标，λ为比例因子，aWGS84和bWGS84分别为WGS84椭球的长半轴和短半轴，h为目标定位点在物方的高程。
4.根据权利要求3所述星地协同的光学卫星在轨实时几何定位方法，其特征在于：步骤1中，所述定位求解方法采用基于DEM数据的高程迭代定位求解方法，实现方式如下，
以目标点高程初值h0＝0，在DEM数据的支持下，进行高程迭代求解，包括执行以下步骤，
a.令i＝1，目标点高程h＝h0，即令h＝0，代入椭球面模型；
b.共线方程与椭球面方程联立，求得目标点物方坐标，获取光线与椭球高h处的交点Mi；
c.若i＞1，则判断此次获得的交点Mi与上次计算坐标Mi-1的修正量d(Mi-1,Mi)是否小于阈值d；
d.若修正量小于阈值，则输出定位结果，若i＝1或修正量大于阈值，则由目标点物方坐标Mi在DEM上内插更新高程值h＝h(Mi)，令i＝i+1，返回步骤b，重复b、c、d步骤直至收敛。
5.根据权利要求1或2或3或4所述星地协同的光学卫星在轨实时几何定位方法，其特征在于：步骤2中，确定星上模型参数的初值，实现方式如下，
对于相机在平台上安装关系的实验室检校，获取三安装角参数，与构成的三个旋转角一致，直接作为初值；
对于相机内部参数的实验室检校，按照严格物理模型，测量相机主距f、各片CCD首像元在相机坐标系下的坐标(x0j,y0j)，设定内定向参数的初始值，如下式，
其中，pixelsize为CCD像元大小设计值。
6.根据权利要求1或2或3或4所述星地协同的光学卫星在轨实时几何定位方法，其特征在于：步骤4中构建结合内部探元指向角模型与外部安装矩阵补偿的在轨几何定标模型，实现方法如下：
作为外定标参数，用于恢复相机坐标系在空间中的位置和姿态；(ax0j,ax1j,ax2j,ax3j,ay0j,ay1j,ay2j,ay3j)(j＝1,2,...,m)作为内定标参数，用于确定相机内部CCD各探元在相机坐标系下的坐标。
7.根据权利要求1或2或3或4所述星地协同的光学卫星在轨实时几何定位方法，其特征在于：步骤4中，基于步骤1所得在轨定位模型，执行以下步骤，
a.设在待定标影像上量测了K个高精度地面控制点作为定向点，控制点的WGS84地心直角坐标为(Xi Yi Zi)，像点坐标为(si li)，i＝1,2,3...k；
b.令：
其中，为像点光线在本体坐标系下的矢量，(pitch,roll,yaw)为相机在本体坐标系上的三个安装偏置角，为本体坐标系到相机坐标系的旋转矩阵
设外定标参数XE、内定标参数XI为自变量，令F()、G()分别为像空间坐标系下像点沿轨及垂轨方向的矢量残差函数，则有：
c.对外定标参数XE、内定标参数XI赋初值
d.将当前的内定标参数XI视为“真值”，将外定标参数XE视为待求的未知参数，将相应的当前值代入矢量残差函数，对每个定向点，进行线性化处理，建立误差方程式，利用最小二乘平差计算更新外定标参数XE的当前值，
e.重复步骤d，迭代计算直至外定标参数改正数均小于预设阈值时停止，进入步骤f；
f.将当前的外定标参数XE视为“真值”，将内定标参数XI视为待求的未知参数，将相应的当前值代入矢量残差函数，对每个定向点，进行线性化处理，建立误差方程式，利用最小二乘平差计算更新内定标参数XI的当前值，
g.重复步骤f，迭代计算直至内定标参数改正数均小于预设阈值时停止，完成几何内外定标参数求解。
8.一种星地协同的光学卫星在轨实时几何定位系统，包括以下模块：
定位模型构建及算法固化模块，用于构建适用于星上实时处理单元的光学卫星成像定位模型，并将相应定位求解方法固化于星上硬件环境，保留模型参数更新上注接口；所述光学卫星成像定位模型，是采用线阵CCD探元指向角的内定向模型，建立基于严密共线方程模型与地球椭球面模型相交的在轨定位模型；
初值确定模块，用于由地面实验室对相机内部及平台安装关系的检校参数或设计参数获取星上定位模型参数的初始值；
定标数据获取模块，用于卫星在轨运行后对地面定标场进行成像，获取适宜几何定标的影像数据并下传至地面系统；
地面系统几何定标模块，用于在地面处理系统中完成定标控制点的密集匹配、定标参数解算；
定标结果验证与模型参数上注更新模块，用于对定标精度评价，确定定标结果的正确性后，更新星上定位的相应参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610786005.0A CN106403902B (zh) | 2016-08-31 | 2016-08-31 | 一种星地协同的光学卫星在轨实时几何定位方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610786005.0A CN106403902B (zh) | 2016-08-31 | 2016-08-31 | 一种星地协同的光学卫星在轨实时几何定位方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106403902A CN106403902A (zh) | 2017-02-15 |
| CN106403902B true CN106403902B (zh) | 2019-01-04 |
Family
ID=58000841
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610786005.0A Active CN106403902B (zh) | 2016-08-31 | 2016-08-31 | 一种星地协同的光学卫星在轨实时几何定位方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106403902B (zh) |
Families Citing this family (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107085856B (zh) * | 2017-04-06 | 2020-01-31 | 上海航天测控通信研究所 | 一种基于光学影像的在轨高精度实时定位方法 |
| CN107644435B (zh) * | 2017-09-19 | 2020-09-22 | 武汉大学 | 顾及姿态校正的敏捷光学卫星无场地几何标定方法及系统 |
| CN107745830B (zh) * | 2017-09-29 | 2020-08-18 | 浙江大学 | 一种基于微小卫星的高分遥感定标系统及定标方法 |
| CN107679748B (zh) * | 2017-09-30 | 2020-10-30 | 哈尔滨工业大学 | 面向星群观测任务自主规划的星地联合运行方法 |
| CN107945235B (zh) * | 2017-10-17 | 2022-02-01 | 许昌学院 | 一种高轨大面阵静止卫星影像几何定位仿真方法 |
| CN109039419A (zh) * | 2018-06-15 | 2018-12-18 | 上海卫星工程研究所 | 天基全球实时连续监视网络系统 |
| CN109143287A (zh) * | 2018-08-15 | 2019-01-04 | 上海航天电子通讯设备研究所 | 定位解算处理系统 |
| CN109612447B (zh) * | 2018-12-29 | 2020-12-15 | 湖南璇玑信息科技有限公司 | 遥感影像地图数据的增强定位变换模型的构建方法、增强定位方法及增强定位服务器 |
| CN110046430B (zh) * | 2019-04-17 | 2022-04-15 | 武汉大学 | 基于椭球分层大气折射的光学卫星影像精确对地定位方法 |
| CN110211054B (zh) * | 2019-04-28 | 2021-01-15 | 张过 | 一种星载推扫式光学传感器无畸变影像制作方法 |
| CN110503692A (zh) * | 2019-07-31 | 2019-11-26 | 武汉大学 | 基于小范围定标场的单线阵光学卫星几何定标方法及系统 |
| CN110986998B (zh) * | 2019-10-28 | 2021-09-14 | 武汉大学 | 一种基于有理函数模型的卫星视频相机在轨几何定标方法 |
| CN111508028A (zh) * | 2020-04-09 | 2020-08-07 | 武汉大学 | 光学立体测绘卫星相机的自主在轨几何定标方法及系统 |
| CN111524196B (zh) * | 2020-04-30 | 2023-04-07 | 中国科学院微小卫星创新研究院 | 一种摆扫大幅宽光学卫星在轨几何标定方法 |
| CN111561935B (zh) * | 2020-05-19 | 2022-03-15 | 中国科学院微小卫星创新研究院 | 旋转大幅宽光学卫星在轨几何标定方法及系统 |
| CN111798523B (zh) * | 2020-06-19 | 2022-09-02 | 武汉大学 | 星相机在轨定标定姿及遥感影像几何定位方法、系统 |
| CN111912430B (zh) * | 2020-08-13 | 2022-07-01 | 中国人民解放军海军航空大学 | 高轨光学卫星的在轨几何定标方法、装置、设备及介质 |
| CN112212833B (zh) * | 2020-08-28 | 2021-07-09 | 中国人民解放军战略支援部队信息工程大学 | 机械拼接型tdi ccd推扫相机整体几何平差方法 |
| CN112435301A (zh) * | 2020-11-18 | 2021-03-02 | 中国科学院上海技术物理研究所 | 一种基于恒星轨迹的遥感相机在轨几何定标方法 |
| CN113311421A (zh) * | 2021-05-24 | 2021-08-27 | 北京市遥感信息研究所 | 一种目标高精度星上实时定位解算系统 |
| CN113393571B (zh) * | 2021-06-15 | 2023-06-27 | 北京市遥感信息研究所 | 一种无云卫星图像生成方法及装置 |
| CN113536485B (zh) * | 2021-07-20 | 2022-12-06 | 中国科学院西安光学精密机械研究所 | 电离层成像探测器图像地理坐标解算方法 |
| CN113900122B (zh) * | 2021-09-27 | 2024-05-14 | 武汉大学 | 星地联合的面阵成像遥感卫星全自主几何定标方法及系统 |
| CN113900123B (zh) * | 2021-09-27 | 2024-07-26 | 武汉大学 | 一种基于rpc的面阵成像遥感卫星在轨几何定标方法及系统 |
| CN114383632B (zh) * | 2021-12-23 | 2023-09-29 | 北京市遥感信息研究所 | 一种基于均方根误差的光学星上处理目标定位精度评估方法 |
| CN114241064B (zh) * | 2022-02-24 | 2022-05-17 | 中国科学院空天信息创新研究院 | 一种遥感卫星内外方位元素实时几何定标方法 |
| CN118129798B (zh) * | 2024-04-30 | 2024-07-30 | 中国科学院空天信息创新研究院 | 天回归轨道遥感卫星的几何定标方法、装置、设备及介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4652884A (en) * | 1984-07-20 | 1987-03-24 | Deutsche Forschungs-Und Versuchsanstalt Fur Luft-Und Raumfahrt E.V. | Satellite navigational system and method |
| US7613565B2 (en) * | 2005-01-07 | 2009-11-03 | Mitac International Corp. | Voice navigation device and voice navigation method |
| CN101833090A (zh) * | 2010-03-12 | 2010-09-15 | 中国科学院遥感应用研究所 | 利用全球卫星定位系统信号源的机载海洋微波遥感系统 |
| CN103968838A (zh) * | 2014-05-09 | 2014-08-06 | 哈尔滨工程大学 | 一种基于极坐标系的auv曲线运动状态下的协同定位方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7917155B2 (en) * | 2006-10-27 | 2011-03-29 | Roundtrip Llc | Location of cooperative tags with personal electronic device |
-
2016
- 2016-08-31 CN CN201610786005.0A patent/CN106403902B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4652884A (en) * | 1984-07-20 | 1987-03-24 | Deutsche Forschungs-Und Versuchsanstalt Fur Luft-Und Raumfahrt E.V. | Satellite navigational system and method |
| US7613565B2 (en) * | 2005-01-07 | 2009-11-03 | Mitac International Corp. | Voice navigation device and voice navigation method |
| CN101833090A (zh) * | 2010-03-12 | 2010-09-15 | 中国科学院遥感应用研究所 | 利用全球卫星定位系统信号源的机载海洋微波遥感系统 |
| CN103968838A (zh) * | 2014-05-09 | 2014-08-06 | 哈尔滨工程大学 | 一种基于极坐标系的auv曲线运动状态下的协同定位方法 |
Non-Patent Citations (1)
| Title |
|---|
| C-GPS 协同式卫星定位技术研究;敖富江，龚国辉，李思昆;《测控技术》;20041231;第23卷(第4期);61-63页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106403902A (zh) | 2017-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106403902B (zh) | 一种星地协同的光学卫星在轨实时几何定位方法及系统 | |
| CN103674063B (zh) | 一种光学遥感相机在轨几何定标方法 | |
| CN104897175B (zh) | 多相机光学推扫卫星在轨几何定标方法及系统 | |
| CN110500995B (zh) | 利用rpc参数建立高分辨率卫星影像等效几何成像模型的方法 | |
| CN105091906B (zh) | 高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 | |
| CN103557841B (zh) | 一种提高多相机合成影像摄影测量精度的方法 | |
| CN106780321B (zh) | 一种cbers-02卫星hr传感器影像整体严密定向与纠正拼接方法 | |
| CN105698764B (zh) | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 | |
| CN107644435B (zh) | 顾及姿态校正的敏捷光学卫星无场地几何标定方法及系统 | |
| CN110006452B (zh) | 高分六号宽视场相机相对几何定标方法及系统 | |
| CN103129752B (zh) | 一种基于地面导航的光学遥感卫星姿态角误差动态补偿方法 | |
| CN103697864B (zh) | 一种基于大虚拟相机的窄视场双相机影像拼接方法 | |
| CN102901519B (zh) | 一种基于探元指向角光学推扫卫星在轨分步几何定标方法 | |
| CN107564057B (zh) | 顾及大气折光校正的高轨面阵光学卫星在轨几何标定方法 | |
| CN106885585B (zh) | 一种基于光束法平差的星载摄影测量系统一体化检校方法 | |
| CN104807477B (zh) | 一种基于靶标控制点的卫星ccd阵列影像几何检校方法 | |
| CN109696182A (zh) | 一种星载推扫式光学传感器内方位元素定标方法 | |
| CN102519433B (zh) | 一种利用rpc反演星载线阵传感器几何定标参数方法 | |
| CN106873004B (zh) | 基于太阳高度角自适应的高轨面阵相机在轨几何检校方法 | |
| CN111473802A (zh) | 一种基于线阵推扫的光学传感器内方位元素定标方法 | |
| CN104820984A (zh) | 一种卫星遥感立体影像处理系统及方法 | |
| CN115326025B (zh) | 一种用于海浪的双目影像测量与预测方法 | |
| CN114838740A (zh) | 一种考虑不同经纬度区域的卫星图像几何定标方法 | |
| Pi et al. | On-orbit geometric calibration using a cross-image pair for the linear sensor aboard the agile optical satellite | |
| CN109029379B (zh) | 一种高精度小基高比立体测绘方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |