专利名称：一种gps/ins数据直接定向精度评估方法
技术领域：
本发明涉及摄影测量技术领域，尤其是涉及一种GPS/INS数据直接定 向精度评估方法。
背景技术：
获取影像的外方位元素是摄影测量的关键问题之一，很多摄影测量工 作者一直为解决如何快速高效地获取影像的外方位元素的问题而努力着。 长期以来，人们主要用空中三角测量技术，通过大量的外业控制点来间接 解求外方位元素，上述方法既浪费经费，又延长了成图周期，严重影响国 民经济建设所急需的资料。20世纪90年代以来，国际摄影测量界开始研 究釆用GPS/INS组合导航系统(即GPS/INS系统GPS系统和INS惯性导 航系统的组合系统)来获取像片摄影时的空间方位，以完全取代摄影测量 加密，从而实现无地面控制的空中三角测量航空遥感已成为国内航空摄影 测量的一种趋势。
而如何评估GPS/INS系统所得到的外方位元素精度成为目前迫切需要 解决的问题，常用的评估方法有两种 一种是将GPS/INS系统所得的外方 位元素数据即GPS/INS数据导入立体测图仪建立立体模型，人工立体观测 获取地面检查点的物方坐标，再与实测检查点坐标对比来评估；另一种是 利用高精度量测的检查点的像点坐标，经点投影系数法两片前交计算得到 检查点的物方坐标，再与实测检查点坐标对比来评估。
第一种方法从常规立体测图入手，也是生产中常容易想到的评估外方 位元素的方法。立体测图仪是成熟的用于数字摄影测量的系统，将外方位 元素数据导入立体测图仪，即可在两张像片组成的立体模型中，人工观测 地面控制点而获取地面控制点物方坐标。但是这种方法存在两个不足1、由于GPS/INS系统获取的外方位元素存在误差，在立体测图仪中存在很大
视差，很难进行立体测图，在视差的影响下，就不能正确量测到所需的地 面控制点，导致像点坐标的量测存在较大误差，这样量测的地面控制点坐 标的误差中有很大一部分来自像点坐标误差的影响，不能严密地评估外方
位元素的数据质量。2、地面一个检查点会出现在多个立体像对中，需要
人工逐个模型地观测，再将多个模型观测得到的坐标值平差得到最后的坐 标值，同实测点坐标对比，才能比较客观地评估外方位元素质量，这是一 个非常费工费时的工作，不仅浪费经费，而且也增加了质量评估的时间。 正是因为上述问题的存在，导致本方法的弱点不可避免，很难很好地应用
到实际的GIS/INS数据评估工作中。从实际用立体测图仪人工量测的
GPS/INS数据直接定向的精度分析表可看出，由于人工成分太多，在存在
视差的条件下，量测的像点坐标误差太多，无法正确评估直接定向精度，
分析所得数据无法用于评估工作。
另 一种方法是从摄影测量理论入手，进行摄影测量立体像对空间前方
交会中的点投影系数法空间前方交会。用空三加密所得到的高精度的量测
像点坐标，联合GPS/INS系统所得的外方位元素数据，代入点投影系数的
空间前方交会公式，求出地面控制点坐标。空间前方交会公式如下 义=IS1 +碼=X'引+ & + iVX2
r = & + m; = ^ + & + tvt2 =4+| (M +群2 + A)
Z = ZS1 =ZS1 +5Z +iV'Z2 式(1 )
详细内容可参考摄影测量学教材中立体像对空间前方交会一节。这种 方法有效避免了量测像点坐标的误差，理论上应该符合数学模型。但存在
以下两种缺点1、经过实验证明这种方法不能符合实际摄影测量生产 情况，只是数学理想化的双线前交，实际加密生产中的光東法平差是实现 多个同名光東最佳地交会到一起，并不能真是交会到一起，只有由同一个 地面控制点和左右片的外方位元素通过共线方程计算出来的理论的像点 坐标数据，再和左右片的外方位元素数据计算前方交会才符合该方法的数
7学模型。即使使用高精度的空三加密成果数据，通过该方法前方交会计算， 得到的地面控制点物方坐标的误差也很大，尤其是高程误差大，基本不符 合实际情况，无法反映外方位元素数据的质量，所以不能正确评估GPS/INS
的数据质量，算法不够严密。2、点投影系数法只限于两片前方交会，不
能实现多光线的交会，受影像数约東。正是由于该方法的以上不足，所以
只能作为理论的前方交会模型方法，不能真正用于实际的GPS/INS数据直 接定向精度评估工作。据实际分析得出的常规空中三角测量所得的像片外 方位元素定向的精度统计结果，其空三精度足以满足丘陵地区的规范要 求；但用常规空中三角测量所得的外方位元素数据和像点坐标数据，通过 点投影系数法两片前方交会所得的精度统计结果，其误差大大超过了规范 要求，甚至不合常理，可见，用点投影系数法两片前方交会的方法不能正 确评估外方位元素数据质量，不能用于GPS/INS所得外方位元素直接定向 的精度评估工作。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一 种GPS/INS数据直接定向精度评估方法，其评估精度高、量测误差小，评 估结果能直接应用于生产实践，能有效避免人工立体量测中增加的像点坐 标量测误差，并且其算法不受像片数约東。
为解决上述技术问题，本发明釆用的技术方案是 一种GPS/INS数据直 接定向精度评估方法，其特征在于该方法包括以下步骤
步骤一、通过航空摄影方法按照航空摄影规范摄取待测区域的多张摄
影像片并通过像片量测方法获取各张像片中地面控制点的像点坐标观测 数据，同时通过GPS/INS系统获取各张像片摄影时的外方位元素，所述外 方位元素包括外方位角元素 p 、 《和以及摄影中心点坐标(Xs ， & ， zs ); 并通过常规摄影测量方法获取所述各张像片的内方位元素;c。、 y。和/;所
述多张摄影像片的数量为地面每个控制点在多张像片上成像的像片数，设
8为N张；
釆用常规实地测量方法，测量出待测区域内所有控制点的实测坐标数
据且将所有实测坐标数据传送至中央处理器；
步骤二、将所获取的N张摄影像片对应的像点坐标观测数据、内外方 位元素、每个地面控制点对应的摄影像片数N均传送至中央处理器进行分 析处理求取待测区域内所有地面点的最佳物方坐标及其坐标理论误差，对 于待测区域内一地面控制点Ml即物点而言，所述中央处理器的分析处理 过程如下
第一步、获取地面控制点M1所对应的摄影像片数N;
第二步、利用N张摄影像片对应的外方位元素^、^和K空间方位元素，
计算各张像片的旋转矩阵，第i张像片的旋转矩阵《=
^ a2 a3
4 & ^ A c2 c3
,其中，
i=l、 2、 3…N;
第三步、求出N张像片中N个同名点对应的像点坐标计算值和误差方 程式的常数项
首先，利用共线方程公式列出N张像片中N个同名点对应的N组共线 方程和N组误差方程式(因此，对于任何地面控制点而言，对应N(N22) 个同名像点可列出2N ( 2"^4 )个共线方程和误差方程，但未知数只有三 个，所以有2N-3个多余观测)，对于第i张像片而言，其同名点对应的
共线方程和误差方程式分别为
a)
"3(x-^)+63(y-}g+C3(z-zs)
"2(us)+Z)2(r-ys)+c2(z-zs)
和<
v =(小x + 一AX + — + —AZ \ , 汉 W
3 (2)，式(1)和(2)中(X， Y, Z)
、" az w 况
为地面控制点Ml对应的未知坐标数据，(x， y)在第i张像片中为地面控制点Ml对应像点的像点坐标观测数据；为像点坐标数据改正
数；
假定(X, Y， Z)的初值为(0， 0, 0),将(0, 0, 0)作为地面控 制点Ml坐标的近似值代入式(l)，求出第i张像片中同名点对应的像点坐
"K)化(r-&)+c2(z-Zs)，并相应计算出
标计算值(x)和(y)
第i张像片中同名点对应的误差方程式的常数项/.和"
|X = X - CO
第四步、逐点列出N张像片中N个同名点对应的N个误差方程式 对于第i张像片而言，其误差方程式V产AiQ-Li,其中R二
V V
/ / L y」 a12"13
—"21a22a23 —
,Q = [AX， AY， AZ];
Ai中的各偏导数为
a,，
a12 =
a13 =
a21 =
22
3x —
i_ —
5_ ，= 双
"l/ + "3(X —X。)
a3(X-
V + 63(H0)
-A)+,rs)+C3(z
q/ + c3(x_x0)
-zs)"3(r-}g+C3(z-
a2/+"30;->;o)
"3(x-xs)+63(y-^)+C3(z-zs)
"3(Us)+,l^)+C3(Z-Zs) c2/ + c3(>;->;0)
总误差方程式 V= V2 ... VN
况a3(ns)+63(n;)+c3(z-zj，
第五步、将逐点列出的N个误差方程式合成为 V=AQ-L(3),其中Q:[AX， AY， AZ], A= [A! A2 ... AN: L= [Li L2 LN]T;
第六步、根据最小二乘间接平差原理，列出式(3)的法方程式 ATPAQ=ATPL(4)，式(4)中，P为单位矩阵；求出式(4)解的表达式为
10而可求出地面控制点Ml坐标近似值的改正数 AX，AY!和AZ"
第七步、判断改正数AX,, AY,和AZ,是否小于限差，若小于限差，计算 此时地面控制点M1改正后的坐标(X!， Y1; ZJ ，其中X产0+A&， Y产0+AY, 和Z严0+AZ"坐标(Xi, Y" ZJ即为地面控制点Ml的最佳物方坐标，地 面控制点M1的分析处理过程完毕，转入第九步；
若改正数AX,， A乂和AZ,不小于限差，则转入第八步；
第八步、以第七步改正后的坐标作为地面控制点Ml坐标的近似值代 入式(l)，重复第三步至第七步，经多次迭代运算后，直至第六步中所计 算得出的改正数AXk， AYk和AZk小于限差，此时地面控制点Ml改正后的坐 标为(Xk， Yk， Zk)，其中Xk-0+AXJ…+ AXk， Y^O+AY, +…+AYk, Z^O+AZ, +… + AZk, k为迭代的次数；坐标(Xk， Yk， Zk)即为地面控制点M1的最佳物 方坐标，地面控制点M1的分析处理过程完毕，转入第九步；
第九步、通过第六步中所述法方程式中未知数矩阵Q的系数矩阵的逆 矩阵(^4'来求解地面控制点Ml坐标的理论精度即理论误差，此时视地面 控制点M1的像点坐标观测数据为等精度不相关观测值，由于矩阵t^力-'中
第f个主对角线上的元素2#就是法方程式中第f个未知数的权倒数，则 第f个未知数的中误差 二^/^m。，其中f-l、 2和3, m。为单位权中误差
且附。"J^P^，其中—]-iX+《，则地面控制点Ml坐标的理论误差
附迴=附i,附少組=附2和附固=附3;
第十步、将计算得出的地面控制点mi的最佳物方坐标与地面控制点
Ml的实测坐标数据进行作差比较，得出相应差值AX腿，AY^和AZm,;
步骤三、重复步骤二，计算下一地面控制点M2的最佳物方坐标及其 坐标理论误差，直至求出待测区域内所有地面控制点的最佳物方坐标及其
坐标理论误差附一,附貞和m一；并且将待测区域内所有地面控制点的最佳物
方坐标与各自的实测坐标数据进行作差比较，并得出相应多组差值
11AXMj， AYmj和AZmj,其中j = l、 2…t， t为待测区域内所有地面控制点的数
步骤四、计算直接定向精度结果，即所有地面控制点坐标的实际中误 差mx、 my 、 ms 、 mh、实际最大误差AXmax， AYm 、 AS,和AH歸和理论中误 差mlx、 mly、 mls、 mlh:
首先，计算多组差值AXMj， AYMj和AZMj的坐标值中误差 jAXM12+... + AXMt2=mx、 |AYM12—. + AYMt2 =my和卜Zj +... + AZMt2 =mz=mh同
Vt t ~ t
时，相应计算得出平面中误差ms = VTOc2+，2 ;
从AXM」'AYMj中计算所有地面控制点的平面误差值asmj = V《+A《， 从AXm」'aymj, ASMj和AZM」筛选出最大的坐标值差值即坐标值最大误差
△xmax， AYmax， ASmax^ AZmiK=AHmax ;
从所有地面控制点坐标的理论误差，计算得出总的理论中误差
mlx= mxM1 +'" + mxMt , mly= myM1 +.■. +,mlh= mhM1 +". + mhMt 、 V t ， "V t ， "\( t 、
ml S = iJot/x2 +m/_y2 。
步骤一中所述的常规实地测量方法为常规野外控制点测量方法。 本发明与现有技术相比具有以下优点l、设计新颖，本发明以摄影 测量的重点方程共线方程为基础，以光束法多片前方交会为理论依据，结 合最小二乘原理平差计算得到地面控制点物方坐标，再经过与地面控制点 的实测坐标比较，统计得到坐标的误差值来评估GPS/INS系统所得外方位 元素直接定向的精度。2、自动化程度高，不需要人工参与量测；而传统 立体测图仪人工量测方法自动化程度低，人工参与成分大，手工工作量大， 且手工工作难度较大，效率较低。3、本发明釆用光東法原理进行计算处 理，算法严密，能正确反映外方位元素质量；传统点投影系数法两片前交 方法算法不严密，不能正确反映外方位元素质量。4、本发明釆用多片最 小二乘计算，即采用多像空间前方交会方法计算得出地面控制点的物方坐 标，因而不受像片数约東；传统立体测图仪人工量测方法为两片立体模型量测，因而受像片数约東；传统点投影系数法两片前交方法为两片前方交 会，因而也受像片数约束。5、本发明用加密的像点坐标，能有效避免像
点坐标误差问题；传统立体测图仪人工量测方法存在视差，不能有效避免
像点坐标误差。综上所述，本发明从摄影测量的核心方程共线方程入手， 结合最小二乘原理，利用一点多度重叠的多余观测条件，实现地面控制点
最佳值的计算，完成GPS/INS数据直接定向精度的评估，其不仅能够避免
人工立体量测中增加的像点坐标量测误差，而且能充分利用一个地面控制
点的多度重叠的条件和数据，增加平差条件，其算法不受像片数约束；同 时，其也能结合实际生产现实，实现最小二乘光東法平差，算法理论严密、 计算出地面控制点坐标的最佳值，能真正实现GPS/INS数据直接定向精度 的评估。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
图1为本发明的流程框图。
具体实施例方式
如图1所示，本发明所述的GPS/INS数据直接定向精度评估方法，包 括以下步骤
步骤一、通过航空摄影方法按照航空摄影规范摄取待测区域的多张摄 影像片并通过像片量测方法获取各张像片中地面控制点的像点坐标观测 数据，同时通过GPS/INS系统获取各张像片摄影时的外方位元素，所述外 方位元素包括外方位角元素P、 w和/c以及摄影中心点坐标(&， Zs); 并通过常规摄影测量方法获取所述各张像片的内方位元素x。、 y。和/;所
述多张摄影像片的数量为地面每个控制点在多张像片上成像的像片数，设 为N张。
釆用常规实地测量方法，测量出待测区域内所有控制点的实测坐标数
13据且将所有实测坐标数据传送至中央处理器。
本实施例中，所述常规实地测量方法为常规野外控制点测量方法。求 取内方位元素x。、 y。和/的方法为相机检校和内定向等方法。
步骤二、将所获取的N张摄影像片对应的像点坐标观测数据、内外方
位元素、每个地面控制点对应的摄影像片数N均传送至中央处理器进行分
析处理求取待测区域内所有地面点的最佳物方坐标及其坐标理论误差，对
于待测区域内一地面控制点Ml即物点而言，所述中央处理器的分析处理 过程如下
第一步、获取地面控制点M1所对应的摄影像片数N;
第二步、利用N张摄影像片对应的外方位元素^、w和K空间方位元素，
计算各张像片的旋转矩阵，第i张像片的旋转矩阵^:
(32 (73 & ~ ^
Cj c
,其中，
i-l、 2、 3…N;
第三步、求出N张像片中N个同名点对应的像点坐标计算值和误差方 程式的常数项
首先，利用共线方程公式列出N张像片中N个同名点对应的N组共线 方程和N组误差方程式，对于第i张像片而言，其同名点对应的共线方程
和误差方程式分别为
少— V _
和<
一;c十一AZ + —AF + —AZ
汉
^3 (2),式(1)和(2)中(X, Y， Z)
L y 、,/ , 汉 ay 况
为地面控制点Ml对应的未知坐标数据，(x, y)在第i张像片中为地面 控制点Ml对应像点的像点坐标观测数据；为像点坐标数据改正
数；
假定(X， Y， Z)的初值为(0， 0, 0),将(0, 0, 0)作为地面控
14制点Ml坐标的近似值代入式(1)，求出第i张像片中同名点对应的像点坐
(T)_ "(Us) + W"s) + q(Z-Zs)
"3^_A)+63(,_ys)+e3(，_zs)并相应计算出
标计算值(x)和(y):
"2(X-Xs)+62(W)+c2(Z-Zs)
"3(x-xs)《-rs)+C3(z-zs)
第i张像片中同名点对应的误差方程式的常数项/.和乙
|7义=X - (X)
lW-(力'
第四步、逐点列出N张像片中N个同名点对应的N个误差方程式 对于第i张像片而言，其误差方程式V严AiQ-Li,其中K"、 、f,<formula>formula see original document page 15</formula>
第五步、将逐点列出的N个误差方程式合成为一总误差方程式 V=AQ-L(3),其中Q-[AX， AY， AZ]， A= [A! A2 ... AN]T, V= A= ^ V2 VN] T， L= [Li L2 … LN] t;
第六步、根据最小二乘间接平差原理，列出式(3)的法方程式 ATPAQ=ATPL(4)，式(4)中，P为单位矩阵；求出式(4)解的表达式为 Q = (#^—^7 (4),从而可求出地面控制点Ml坐标近似值的改正数 AXP AY和AZ"
第七步、判断改正数AX,， AY,和AZ,是否小于限差，若小于限差，计算此时地面控制点M1改正后的坐标(Xi, Y丄,ZJ ，其中X产0+AX" Y产0+AY' 和Z产0+AZ。坐标(Xn Y" ZJ即为地面控制点Ml的最佳物方坐标，地 面控制点M1的分析处理过程完毕，转入第九步；
若改正数AX,， AY,和AZ,不小于限差，则转入第八步；
第八步、以第七步改正后的坐标作为地面控制点Ml坐标的近似值代
入式(l),重复第三步至第七步，经多次迭代运算后，直至第六步中所计 算得出的改正数AXk， AYk和AZk小于限差，此时地面控制点Ml改正后的坐 标为(Xk， Yk, Zk)，其中X^O+AX, +…+ AXk， Y^O+A^ +…+AYk, Z^O+AZJ… + AZk， k为迭代的次数；坐标(Xk, Yk， Zk)即为地面控制点M1的最佳物 方坐标，地面控制点M1的分析处理过程完毕，转入第九步；
第九步、通过第六步中所述法方程式中未知数矩阵Q的系数矩阵的逆 矩阵(Z々'来求解地面控制点Ml坐标的理论精度即理论误差，此时视地面 控制点Ml的像点坐标观测数据为等精度不相关观测值，由于矩阵(^^广中 第f个主对角线上的元素2#就是法方程式中第f个未知数的权倒数，则 第f个未知数的中误差柳,二^m。，其中f=l、 2和3， m。为单位权中误差
且附±,_H_，其中["=|^+/;，则地面控制点Ml坐标的理论误差
附祐==附2和附顧=附3 ;
同时，需注意的是此处的/ 和~为第i张像片经k次迭代后的误差
方程式的常数项，具体是第三步中计算得出的第i张像片中同名点对应的 误差方程式的常数项t和/y。
第十步、将计算得出的地面控制点Ml的最佳物方坐标与地面控制点 Ml的实测坐标数据进行作差比较，得出相应差值AX組，AY腿和AZm"
步骤三、重复步骤二，计算下一地面控制点M2的最佳物方坐标及其
坐标理论误差，直至求出待测区域内所有地面控制点的最佳物方坐标及其 坐标理论误差附谢,m洲.和附鄉.；并且将待测区域内所有地面控制点的最佳物
方坐标与各自的实测坐标数据进行作差比较，并得出相应多组差值
16AXM」，AYmj和AZmj，其中j = l、 2…t, t为待测区域内所有地面控制点的数
步骤四、计算直接定向精度结果，即所有地面控制点坐标的实际中误 差mx、 my、 ms、 mh、实际最大误差AX,， AY腿、AS,和AH隱和理论中误 差mlx、 mly、 mls、 mlh:
首先，计算多组差值AXMj， AYMj和AZMj的坐标值中误差 jAXM12+... + AXMt2=mx、 |AYM12—. + AYMt2 =my和+AZMt2 =mz=mh同
Vt Vt V t
时，相应计算得出平面中误差ms = >w2+mj;2 ;
从AXM AYMj中计算所有地面控制点的平面误差值ASm」=V《+《， 从AXm' aym」，ASMj和AZMj筛选出最大的坐标值差值即坐标值最大误差
AXmax， AYmax ， ASmaxAZ咖x = AH咖x ;
从所有地面控制点坐标的理论误差，计算得出总的理论中误差
mlx= m顺2+…+ 2 , mly= myM1 + " +，mlh=卜m卩+…+ mhM,2 、 V t ， t ， "\( t 、
综上，要进行GPS/INS系统的数据直接定向的精度评估工作，就要将 GPS/INS数据直接定向的结果和常规空中三角测量的定向结果进行比较，所 以先进行常规的空中三角测量工作并且达到高精度要求，除了得到加密精度 成果之外，从常规空中三角测量中，可以得到高精度的地面控制点的量测像 点坐标，首先保证了前方交会中像点坐标的量测精度。
因而当将本发明用于GPS/INS应用的实验研究中时，先由GPS/INS组合 导航系统来获取像片摄影时的空间方位，即像片的外方位元素^、 &、 Zs、 ^、"、、因精度对比所需，在试验区进行了常规摄影测量即常规空中三角 测量得到了 一套满足精度要求的外方位元素数据和控制点像点坐标数据。从 外方位元素数据和像点坐标数据这些已知数据出发，经上述步骤即步骤一到 步骤四(在Microsoft Visual C++ 6. 0平台上用0++语言开发具体实用程序)， 可得到地面控制点的物方坐标值。
17用本发明统计的常规空中三角测量所得外方位元素定向精度，如表1所
示
表1 本发明统计的常规空中三角测量成果直接定向精度
统计中误差(cm)最大误差(cm)所用mxmymsmhM歸max"^"""1 max(ax
点数
4908.54468.706412. 166. 522019.948521,298029. 734720. 6727
利用常规空中三角测量系统统计的常规空中三角测量所得的像片外方位
元素定向的精度统计结果如表2所示 表2
-f平 高定句点残差(cm)多余控制点较差(cm)中误差最大误差中误差最大误差定检mymsmh"隨厶h眼mxmymsmh△ x隱A h隨
向査
点，点
544808. 27.711. 25. 519. 820. 621. 518. 67.78. 811. 77. 516. 419. 922. 1IS. 7
将表1和表2对比分析可知，表l与表2的统计结果基本一致，可见本 发明能够正确客观地评估外方位元素的定向精度。
选用一个测区的574个控制点的地面控制点坐标，与实测地面控制点数 据对比，即用本发明统计GPS/INS系统所得外方位元素直接定向精度结果如 表3所示
表3 本发明统计的GPS/INS数据直接定向精度
统计 所用实际中误差(cm)实际最大误差(cm)理论中误差(cm)点数raxmymsmhmax《W隨mixmlymismlh
49019. 417. 325. 9935.858. 573. 874.481.95. 56. 08. 212. 7
从表3可以看出，本发明能够完成GPS/INS数据的直接定向精度评估工
18釆用现有的立体测图仪人工量测的GPS/INS直接定向的精度分析结果， 和通过点投影系数法两片前方交会统计的常规空中三角测量所得的像片外方
位元素定向的精度结果分别如表4和表5所示
表4:立体测图仪统计GPS/INS数据直接定向精度结果
统计所 用点数中误差(m)最大误差(m)mxmymsmhAX """^ i maxmax認腿
5100. 2540. 2070. 3280. 1470. 7390. 6690. 8200, 617
表5:点投影系数法统计的常规空中三角测量成果直接定向精度结果
统计中误差(m)最大误差(m)所用 点数mxmymsmhAX max"丄max認隨
50477. 4 297152.877 8584.0448740.2591877. 968483.2799918. 87 27838. 2643
将表4和表5分别与表2对比分析可知，现有的立体测图仪人工量测和 通过点投影系数法两片前方交会方法均不能正确客观地评估外方位元素的定 向精度。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是 根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构 变化，均仍属于本发明技术方案的保护范围内。
19
权利要求
1. 一种GPS/INS数据直接定向精度评估方法，其特征在于该方法包括以下步骤步骤一、通过航空摄影方法按照航空摄影规范摄取待测区域的多张摄影像片并通过像片量测方法获取各张像片中地面控制点的像点坐标观测数据，同时通过GPS/INS系统获取各张像片摄影时的外方位元素，所述外方位元素包括外方位角元素ω和κ以及摄影中心点坐标(XS，YS，ZS)；并通过常规摄影测量方法获取所述各张像片的内方位元素x0、y0和f；所述多张摄影像片的数量为地面每个控制点在多张像片上成像的像片数，设为N张；采用常规实地测量方法，测量出待测区域内所有控制点的实测坐标数据且将所有实测坐标数据传送至中央处理器；步骤二、将所获取的N张摄影像片对应的像点坐标观测数据、内外方位元素、每个地面控制点对应的摄影像片数N均传送至中央处理器进行分析处理求取待测区域内所有地面点的最佳物方坐标及其坐标理论误差，对于待测区域内一地面控制点M1即物点而言，所述中央处理器的分析处理过程如下第一步、获取地面控制点M1所对应的摄影像片数N；第二步、利用N张摄影像片对应的外方位元素ω和κ空间方位元素，计算各张像片的旋转矩阵，第i张像片的旋转矩阵其中，i＝1、2、3...N；第三步、求出N张像片中N个同名点对应的像点坐标计算值和误差方程式的常数项首先，利用共线方程公式列出N张像片中N个同名点对应的N组共线方程和N组误差方程式，对于第i张像片而言，其同名点对应的共线方程和误差方程式分别为和式(1)和(2)中(X，Y，Z)为地面控制点M1对应的未知坐标数据，(x，y)在第i张像片中为地面控制点M1对应像点的像点坐标观测数据；(vx，vy)为像点坐标数据改正数；假定(X，Y，Z)的初值为(0，0，0)，将(0，0，0)作为地面控制点M1坐标的近似值代入式(1)，求出第i张像片中同名点对应的像点坐标计算值(x)和(y)并相应计算出第i张像片中同名点对应的误差方程式的常数项lx和ly第四步、逐点列出N张像片中N个同名点对应的N个误差方程式对于第i张像片而言，其误差方程式Vi＝AiQ-Li，其中Vi＝[vx vy]T，Li＝[lx ly]T，Q＝[ΔX，ΔY，ΔZ]；Ai中的各偏导数为第五步、将逐点列出的N个误差方程式合成为一总误差方程式V＝AQ-L(3)，其中Q＝[ΔX，ΔY，ΔZ]，A＝[A1A2...AN]T，V＝A＝[V1V2...VN]T，L＝[L1L2...LN]T；第六步、根据最小二乘间接平差原理，列出式(3)的法方程式ATPAQ＝ATPL(4)，式(4)中，P为单位矩阵；求出式(4)解的表达式为Q＝(ATA)-1ATL(4)，从而可求出地面控制点M1坐标近似值的改正数ΔX1，ΔY1和ΔZ1；第七步、判断改正数ΔX1，ΔY1和ΔZ1是否小于限差，若小于限差，计算此时地面控制点M1改正后的坐标(X1，Y1，Z1)，其中X1＝0+ΔX1，Y1＝0+ΔY1和Z1＝0+ΔZ1，坐标(X1，Y1，Z1)即为地面控制点M1的最佳物方坐标，地面控制点M1的分析处理过程完毕，转入第九步；若改正数ΔX1，ΔY1和ΔZ1不小于限差，则转入第八步；第八步、以第七步改正后的坐标作为地面控制点M1坐标的近似值代入式(1)，重复第三步至第七步，经多次迭代运算后，直至第六步中所计算得出的改正数ΔXk，ΔYk和ΔZk小于限差，此时地面控制点M1改正后的坐标为(Xk，Yk，Zk)，其中Xk＝0+ΔX1+…+ΔXk，Yk＝0+ΔY1+…+ΔYk，Zk＝0+ΔZ1+…+ΔZk，k为迭代的次数；坐标(Xk，Yk，Zk)即为地面控制点M1的最佳物方坐标，地面控制点M1的分析处理过程完毕，转入第九步；第九步、通过第六步中所述法方程式中未知数矩阵Q的系数矩阵的逆矩阵(ATA)-1来求解地面控制点M1坐标的理论精度即理论误差，此时视地面控制点M1的像点坐标观测数据为等精度不相关观测值，由于矩阵(ATA)-1中第f个主对角线上的元素Qff就是法方程式中第f个未知数的权倒数，则第f个未知数的中误差其中f＝1、2和3，m0为单位权中误差且其中则地面控制点M1坐标的理论误差mxM1＝m1，myM1＝m2和mhM1＝m3；第十步、将计算得出的地面控制点M1的最佳物方坐标与地面控制点M1的实测坐标数据进行作差比较，得出相应差值ΔXM1，ΔYM1和ΔZM1；步骤三、重复步骤二，计算下一地面控制点M2的最佳物方坐标及其坐标理论误差，直至求出待测区域内所有地面控制点的最佳物方坐标及其坐标理论误差mxMj，myMj和mhMj；并且将待测区域内所有地面控制点的最佳物方坐标与各自的实测坐标数据进行作差比较，并得出相应多组差值ΔXMj，ΔYMj和ΔZMj，其中j＝1、2...t，t为待测区域内所有地面控制点的数量；步骤四、计算直接定向精度结果，即所有地面控制点坐标的实际中误差mx、my、ms、mh、实际最大误差ΔXmax，ΔYmax、ΔSmax和ΔHmax和理论中误差mlx、mly、mls、mlh首先，计算多组差值ΔXMj，ΔYMj和ΔZMj的坐标值中误差和同时，相应计算得出平面中误差从ΔXMj，ΔYMj中计算所有地面控制点的平面误差值从ΔXMj，ΔYMj，ΔSMj和ΔZMj筛选出最大的坐标值差值即坐标值最大误差ΔXmax，ΔYmax，ΔSmax和ΔZmax＝ΔHmax；从所有地面控制点坐标的理论误差，计算得出总的理论中误差
2.按照权利要求1所述的一种GPS/INS数据直接定向精度评估方法， 其特征在于步骤一中所述的常规实地测量方法为常规野外控制点测量方 法。
全文摘要
本发明公开了一种GPS/INS数据直接定向精度评估方法，其以摄影测量的重点方程共线方程为基础，以光束法多片前方交会为理论依据，结合最小二乘原理平差计算得到地面控制点物方坐标，再经过与地面控制点的实测坐标比较，统计得到坐标的误差值来评估GPS/INS系统所得外方位元素直接定向的精度。本发明评估精度高、量测误差小，评估结果能直接应用于生产实践，能有效避免人工立体量测中增加的像点坐标量测误差，并且其算法不受像片数约束。
文档编号G01C21/28GK101509784SQ20091002162
公开日2009年8月19日 申请日期2009年3月20日 优先权日2009年3月20日
发明者任为民, 吴晓燕, 丽 汤, 苗小利 申请人:西安煤航信息产业有限公司