利用已有像控点的.doc
《利用已有像控点的.doc》由会员分享,可在线阅读,更多相关《利用已有像控点的.doc(9页珍藏版)》请在三一文库上搜索。
1、精品论文利用已有像控点的多期航空影像光束法区域网联合平差1王建梅 1,朱紫阳 21 同济大学测量与国土信息工程系,上海(200092)2 广东省国土资源厅测绘院,广东广州(510500)E-mail:摘要:为了充分利用已有航空影像与控制成果,节约作业成本,缩短作业周期,本文提出 了一种利用已有像控点联合多期航空影像进行整体平差加密的数学模型,并采用我国南方某 地区的数据资料进行了试验证明。结果表明,联合平差能将已有像控点的控制作用通过多期 影像间的同名点传递到最新航空影像上,联合平差加密点的平面精度与常规光束法区域网平 差加密精度基本一致,而高程精度则取决于多期影像间同名点的观测值个数,同名点
2、观测值 越多,加密点高程精度越高,因此只要多期影像间存在足够多的同名点观测值,则联合平差 的加密精度完全能够满足摄影测量地形图测绘的生产需要。关键词:像控点加密精度空中三角测量联合光束法平差区域网平差 中图分类号:1引言航空摄影测量是利用二维对地观测影像提取三维地表空间信息的重要技术手段,其关键 是快速而准确地恢复影像获取时的空间方位。长期以来,这一目标是借助大量合理分布的地 面控制点通过空中三角测量间接实现的。随着空间定位技术、传感器技术、计算机技术等的 飞速发展,航空摄影测量几何定位方法正朝着无需地面控制点的方向迈进1。20世纪70年代, 随着GPS的出现,人们开始尝试采用载波相位差分GP
3、S 动态定位技术来确定摄影瞬间摄站 的空间位置(即像片的3个外方位线元素)2,3,4,5;进入20世纪90年代,人们又开始研究 GPS/IMU(简称POS系统)来获取像片摄影时的空间方位(即利用GPS确定摄站的空间位置, 利用IMU惯性测量装置获取影像的姿态角),以直接用于航测内业的像片定向6,7,8,这些新 技术的引入确实使摄影测量作业大量减少了地面控制点,降低了生产成本,缩短了航测成图 周期。但大量研究和生产实践同时也表明,不论采用GPS还是POS系统,都不能完全满足大 比例尺地形图测绘的精度要求,仍需要一定数量的地面控制点改正系统误差的影响。近年来随着我国经济建设的高速发展,地表变化日新
4、月异,快速更新基础地理信息是时 代所需,航空摄影测量因其技术成熟、速度快、精度高依旧是基础地理数据更新的主要方式。 鉴于我国大部分地区已经完成了基本比例尺的基础测绘工作,拥有了大量的像控点与航空摄 影成果,因此考虑利用已有像控点完成对同一地区新获影像的空间定位便成为一种自然而然 的想法。但在实际工作中受多种因素的影响,同一地区不同时相的单张航空影像覆盖的地面 范围通常并不一致,致使前期布设的像控点对于后期的航空摄影测量加密不能满足作业规范 要求,更为常见的是由于实地发生各种变化,某些前期布设的像控点已经不复存在。这就提 出了一个问题,能否在不增加野外控制测量的情况下,合理利用已有的像控点和航片
5、资料进 行同一地区新获影像的定位,通过提高作业技术含量来降低成本,提高效率。针对此类问题1本课题得到国家高技术研究发展计划(项目编号:2007AA12Z178)的资助。-9-已有一些研究成果,代表性的如袁修孝教授等研究了基于已知定向参数影像的光束法区域网平差进行后期影像定位9,10,但该方法就实质而言属于二次加密,某种程度上降低了成果的 精度。为了克服此方法的不足,本论文提出了一种利用已有像控点,联合多期影像资料整体 平差解算后期影像外方位元素及其加密点坐标的方法,该方法能更深层次挖掘现有成果资 料,使新获航片达到其原始设计级别的成图精度。由于多期影像存在地区重叠,对于没有发 生变化的区域其同
6、名像点坐标观测值数量增加很多,而对于发生变化的区域,其像点坐标观 测值数也与单期平差一致,因此联合多期影像资料整体平差观测值数量将大大增加,不但增 强了整个区域网的可靠性,而且克服了控制点的漏区,因此具有很好的实用价值。下面首先 介绍该方法的理论基础和数学模型,然后以我国南部某地区为例,验证该方法的有效性与可 行性。2联合平差的数学模型由文献11,12可知,单期自检校光束法区域网平差的误差方程式可写为V = AX + BY + Ct LP(1)式中,V 为像点坐标观测值的改正数向量,X = X YZ T 为加密点坐标未知数改正数向量,Y = XY ZT 为像片外方位元素未知数改正数向量, t
7、为自检校参数向量, L = x ( x)S S Sy ( y)T 为像点坐标观测值的残差向量,其中( x , y )为像点坐标观测值,( x ),( y )为按共线条件方程计算的像点坐标值, A、 B 、 C 为对应的系数矩阵,P 为像点坐标观测值的权矩阵,其中 ( P )ii = 1 2 , 为像点坐标观测值的中误差。0 0由公式(1)我们不难推导出利用已有像控点,联合前(I)后(II)两期影像像点坐标 观测值的自检校光束法区域网平差的误差方程式为V1 = A11 X1=V2+ A13 X 3 + B1Y1A22 X 2 + A23 X 3+ B2Y2+C1t1 L1 P1+C2 t2 L2
8、 P2(2)这里,V1 、V2 分别为期与期的像点坐标观测值改正数向量, X1 、 X 2 分别为期与期独立加密点坐标未知数向量, X 3 为两期共有加密点坐标未知数向量,Y1 、Y2 分别为期与 期的外方位元素未知数改正数向量,t1 、t2 分别为期与期的自检校参数向量,L1 、L2 分别为期与期的像点坐标观测值的残差向量, A11 、 A13 、 A22 、 A23 、B1 、B2 、C1 、C2分别为对应的系数矩阵, P1 、 P2 分别为期与期的像点坐标观测值的权矩阵,其中( P )= 1 2 、 ( P )= 1 2 , 与 分别为期与期像点坐标观测值的中误差。1 i , i 12
9、i , i 2 1 2T对于每一像点坐标观测值根据式(2)建立误差方程式,组成法方程,得到 AT P AOAT P AAT P BOAT P CO X AT P L11 1 11 11 1 13 11 1 1 11 1 11 11 1 1OAT P AAT P AOAT P BOAT P C X AT P L22 2 22 22 2 23 22 2 2 22 2 2 2 22 2 2 AT P AAT P AAT P A+ AT P AAT P BAT P BAT P CAT P C X AT P L + AT P L 13 1 11 23 2 22 13 1 13 23 2 2313 1 1
10、23 2 2 13 1 1 23 2 2 3 13 1 123 2 2 (3)1 1 11 1 1 13 1 1 1 1 1 1B P L BT P AOBT P ABT P BOBT P CO Y1 = 1 1 1OBT P ABT P AOBT P BOBT P C Y BT P L2 2 22 2 2 23 2 22 2 2 2 2 2 2 21 1 11 1 1 13 1 1 1 1 1 1 C T P AOC T P AC T P BOC T P CO t1 TC P L1 1 1OC T P AC T P AOC T P BOC T P C t C T P L2 2 22 2 2
11、23 2 2 2 2 2 2 2 2 2 2如果考虑到外业控制点坐标误差的影响,将控制点坐标值作为观测值看待,则可得到修正后 的两期自检校光束法区域网平差的误差方程式为V1 = A11 X1V2=+ A13 X 3 + B1Y1A22 X 2 + A23 X 3+ B2Y2+C1 t1+ D1 s L1+C2 t2 + D2 s L2P1P2 (4)Vs =+ Es Ls Ps这里V 为控制点虚拟观测值的坐标改正数向量, s = XY Z 为控制点坐标改正数,s c c c2Ls 为控制点坐标虚拟观测值的残差向量(取控制点坐标为初值时, Ls = 0 ), D1 、 D2 为对 A P AOA
12、 P AA P BOA P COA P D应的系数阵, E 为单位权矩阵, Ps 为控制点坐标虚拟观测值的权矩阵,其中 ( Ps )i , i = 1 s , s 为控制点坐标中误差,其它字符含义同式(2)。由式(4)可得到对应的法方程为T T T T T11 1 11 11 1 13 11 1 1 11 1 1 11 1 1 X1 OAT P AAT P AOAT P BOAT P CAT P D X B P AOB P AB P B22 2 22 22 2 23 22 2 2 22 2 2 22 2 2 2 AT P AAT P AAT P A+ AT P AAT P BAT P BAT
13、P CAT P CAT P D + AT P D X 1 1 1 1 1 1 13 1 11 23 2 22 13 1 13 23 2 23 13 1 123 2 2 13 1 1 23 2 2 13 1 1 23 2 2 3 T T T1 1 11 1 1 13 1 1 1OBT P COBT P D Y1 OBT P ABT P AOBT P BOBT P CBT P D Y 2 2 22 2 2 23 2 2 2 2 2 2 2 2 2 2 C T P AOC T P AC T P BOC T P COC T P D t 1 1 11 1 1 13 1 1 1 1 1 1 1 1 1OC
14、 T P AC T P AOC T P BOC T P CC T P D1 t (5)A P L2 2 22 2 2 23 2 2 2 2 2 2 2 2 2 2 DT P ADT P ADT P A+ DT P ADT P BDT P BDT P CDT P CDT PD + DT P D + P s 1 1 11 2 2 22 1 1 13 2 2 23 1 1 1 2 2 2 1 1 1 2 2 2 11 1 2 2 2s T11 1 1TA P LA P L + A P L22 2 2 T T B P LT13 1 123 2 2 =B P L1 1 1 T C P L2 2 2 T1
15、 1 1C P LT 2 2 2 DT P L + DT P L + P L 1 1 12 2 2s s 当两期影像上的同名像点足够多时,利用式(5)就可以求解期影像的外方位元素及 其加密点的坐标。3联合平差实例为了验证所建数学模型的正确性以及上述方法的可行性,本文选择了我国南方某地区的 两期航空摄影成果进行试验,具体摄影技术参数见表 1 所示。试验区位于城市与山区之间, 中间有一条地区主干河道穿越,地形为丘陵地。摄影时间间隔四到五年,在此期间由于经济建设和产业结构的调整,实地地物变化较快, 为期生产布设的很大一部分控制点在期影像上已不复存在。另外由于两期影像摄影比例 尺相差较大,即使保存完好
16、的控制点,从重叠度、像片位置条件等方面考虑,也不能满足 期影像空中三角测量加密的规范要求,为此在期的航空摄影生产作业中根据项目的实际需 要,重新布设了一次控制点。其中期影像航向每隔 3 条基线布设一个平高点、旁向每隔一 条航线布设一个平高点,共 37 个控制点;期影像航向每隔 6 条基线布设一个平高点、旁 向每隔一条航线布设一个平高点,共 52 个控制点。具体像主点和控制点分布见图 1 所示。试验采用 VirtuoZo AAT/PATB 自动空中三角测量软件进行了同名像点坐标量测与粗差探 测,共有像点坐标观测值 13579 个,其中期 3937 个,平均每片 57.9 个,期 9642 个,
17、平均每片 70.9 个,同名像点观测值 6701 个,期独立观测值 1329 个,占 33.8%,期独立观测值 5549 个,占 57.6%。两期像点观测值具体分布见图 2 所示。表 1 两期试验航空摄影影像技术参数Tab1. Technical Data of two sets Aerial Imageries项目影像技术参数期影像(黑白)期影像(真彩色)摄影时间2001 年 10 月2005 年 12 月-2006 年 12 月 航摄仪型号 WILD RC-10WILD RC-30 摄影仪焦距(mm)89.000153.500像幅(cm)23232323 摄影比例尺1:250001:150
18、00-1:20000扫描分辨率() 2521航向重叠度60%60%旁向重叠度30%30% 航线数4 6总像片数68136期像主点 期像主点 期控制点 期像主点图 1 试验区两期影像与控制点分布图Fig1. Distribution Figure of Two Sets Aerial Imageries and Ground Control Points观测值数 (个)600550500450400350后期影像 重叠度300109250 8200 76150 541003250 11 2 3 4 5 6 7 8 9 10前期影像 重叠度图 2 两期加密点观测值分布Fig2. Distribut
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 利用 已有像控点
链接地址:https://www.31doc.com/p-3623996.html