航空摄影测量中的立体模型重建外文翻译.doc

航空摄影测量中的立体模型重建摘要 本文描述的是现代航空摄影测量的操作问题和基本的技术需要。当立体模型重建时，利用航空摄影测量中的外方位元素决定摄影测量点的精度和在对应的模型点中的y-视差分析。真正的航空摄影，在图像的比例，由1:2 500至1:6 0000，与dgps/imu的数据来源于各种地形，在中国由我们的pos-支持的大型区域网平差计划wucaps处理。实证结果证实来源于大型区域网平差的外方位元素的精度符合地形勘测规范的要求。然而，通过pos确定的外方位元素的精度不能满足地形勘测规范的要求。关键词 空中三角测量（at）；gps（全球定位系统）；pos（定位和定向系统）； 立体模型重建；地面控制点（gcps）；精度导言航空摄影测量是从空中影像获得关于地球表面的三维空间信息的科学和技术。摄影点的决定，其中通过使用图像找出地面对象，是依据识别物体的遥感。并且问题的关键是迅速和准确地确定图像的位置和行为上的即时影像。通过基于分布式地面控制点的空中三角测量满足这一目标。随着空间定位技术的发展，遥感技术和计算机科学，以及空中三角测量的演变和发展走向没有地面控制点的数字化勘测。早在1950年，摄影科学家就开始研究如何利用各种辅助数据，以减少地面控制点的需要。然而，由于技术的局限性，方法没有变成现实。直到20世纪70年代，出现了美国的全球定位系统（gps），在航空摄影过程中人们仅得到通过载波相位差分全球定位系统（dgps）技术来确定曝光驻地的位置(即航摄照片的三个线性元素)，用于执行空中三角测量（简称gps-支持 at）可以减少摄影对地面控制点的依赖，缩短测绘周期；并降低生产成本，在摄影测量的领域触发革命。然而，gps-支持 at在空中摄影测量的操作是有利的，主要是在浩大和困难的区域，在中小型的比例尺，而不是带状区域和城市大比例尺测图。在20世纪90年代，人们开始探讨采用gps/lns集成系统(也称pos)获取照片的位置和姿态(即利用gps获得曝光驻地的位置，由imu获得图像姿态元素)，目的是照片的定向，最终目标是取代区域空中。三角测量程序。现代数字摄影测量学在4d产品(dem，dom，dlg，drg)的自动化的生产和空间数据库的更新中将扮演一个重要角色。本文将介绍航空摄影测量学和相关的技术需要在当前的操作应用，特别是，摄影信息链的几何定位精度可从照片方向到立体模型重建获得，旨在探讨4d产品生产的可实行性。人民希望这项研究的研究结果可以在国土普查，地图测绘和基础地理信息采集方面为航空摄影测量的操作提供指导。1 现代航空摄影测量的模式现今，航空摄影测量主要有三种模式，即标准的航空摄影测量、gps-支持航空摄影测量和pos-支持航空摄影测量。它们主要程序如图1所示。从图1，我们可以了解到，区别这三种模式的方法主要是如何获取的航摄照片以及照片的方向。对于标准空中三角测量，它是通过区域空中三角测量的大量地面控制点获得模型定向点的坐标来完成图像的定向。对于gps-支持 at，在航空照片获得的过程中，动态gps定位是用来代替地面控制点以确定曝光中心的位置和获得该模型定向点的坐标，然后用于纠正图像的方向。对于pos-支持 at，图像和它们对应的方位元素(图像的六个外方位元素)都是已获取的，用于了解在曝光时刻几何反演摄影存储的空间位置和姿态。2 相关技术要求2.1 空中摄影在现代航空摄影，为了提高获得的图像质量，除了新增飞行控制系统到空中摄影机以外(例如ascot，ccns4，空中跟踪系统)，当采取gps空中摄影和在照相机上安置pos系统进行dgps/imu空中摄影时，方法还包括牢固黏附与照相机的一台gps接收器。根据空中摄影的不同的模式，我们可以制定一个负责计划如图2所示。2.2 地面控制计划在数字摄影测量工作站，空中三角测量进行了理论上的最严格的大型区域网平差，但为了获得照片最佳的传输点的坐标和方向的外方位元素，地面控制计划的设计应如图3所示，即不同模式的空中摄影。2.3 数字映射从理论上说，在得到准确的内外方位元素的图像之后，可衡量的立体模型可利用模型重建恢复，其中我们可以做地形的测绘以及物体的自动运行。然而，目前的四维产品的生产工艺是：单张照片的内定向立体像对的相对定向单一模型的绝对定向立体模型的测绘。该方法的模型只有通过pos-支持的航空摄影测量直接地理参考恢复。3 实验和分析航摄定位有两种方法。其中之一被称作区域空中三角测量，关于图像点的坐标，地面控制点的坐标和（或）图像的外方位元素加权观测值，并结合大型区域网平差来解决图像定向参数和目标点的空间坐标，来作为方向控制点的立体模型绘图和做高度精确的几何定位的应用。为不同尺度和地形类型的航空摄影测量， 航摄照片办公室操作的地形图规格定义了各自空中三角测量方法，地面控制计划，以及传输点精度的具体标准。这种方法已被建立并得到了广泛的应用。另一种是所谓的直接地理参考，假定高精确的图像外方位元素是可以得到的，在立体像对中通过使用图像坐标系统的同名像点的坐标，利用空间交会计算出对应的目标点物体的空间坐标。这种方法直接地确定对象的位置，因此4d产品可以被生产。然后本文将主要讨论当利用各种方式获得图像的外方位元素时，如何定位精度可以完成立体模型的y视差。3.1 数据如表1所示，来自不同领域的4个组的实际图像的实验被实施。所有的底片被扫描，其分辨率为21米，为了得到连接点，pos-支持的大型区域网平差软件wucaps被用于图像的测试1，测试2和测试4，自制的jx- 4数字摄影测量工作站用于图像的测试3。根据相对定向的结果与由wucaps去除严重错误的功能，地面控制点全部在立体镜方式下手动地被测量，并且所有图像点的统计精度(rms)均优于6.0米。在那以后，我们使用applanix pos/av系统的后处理软件pospac做测试场校准和dgps和imu数据的整合，然后通过应用坐标系变换和系统误差改正每个图像的六个外方位元素，这是由pos系统规定的，可以获得。表1 实验计划中的图像参数测试1测试2测试3测试4时间2004.112005.12005.92005.10飞机yun-12yun-12yun-8aviation-照相机leica rc-30leica rc-30leica rc-30leica rc-30班次控制系统航空轨道航空轨道ccns 4ccns 4pos系统pos av 510pos av 510pos av 510pos av 510gps接收机astechtrimble 5 700trimble 5 700trimble 5 700底片kodak 2 444kodak 2 044kodak 2 402kodak 2 402距离原则/mm153.84303.64154.06153.53框架/cmcm2323232323232323相片比例1:2 5001:3 0001:32 0001:60 000前部重叠度/%61636464边缘重叠度/%32333330航带数量91094控制带2220相片25537724448gcps731603429密集点3 6315 4422 951712区域/kmkm455847524057最大地带波动/m38.6（平坦地区）181.6（山地）729.3（高山地区）109.3（低地）gps更新比率/s20.511gps时间设定/min101055静态gps/min5555人工gps/m0.303，0.110，-2.0290.303，-0.110，-2.002-2.015，-0.030，3.1022.034，-0.520，1.320人工imu/m0.000，0.200，-0.5590.000，0.200，-0.7100.000，-0.201，0.427-0.006，-0.202，0.4303.2 外方位元素的性能为了得到通过不同的方法分析外方位元素的性能，标准at与边缘密集地面控制点和gps-支持 at与四组完全地面控制点在基础领域被首先执行，可以获得每幅图像的六个外方位元素，并且它们的理论精度可以估计。然后我们假设标准的结果作为“真理”并估计通过pos提供的外方位元素的性能。结果如表2所示。从表1我们看到的图像的测试1，测试2可用于生产规模为1:5001:2 000的四维产品，图像的测试4可用于生产规模为1:50001:10 000的四维产品。在空中摄影测量办公运行的地形图规格原则中，测试1，测试2，测试3和测试4分别属于平坦的土地，山地，高山区的土地和低地。如表2的结果，一些结论可归纳如下。表2 不同方法的外方位元素精度的确定图像方法/gcps控制 pts控制pts残差的rmseo的理论精度水平垂直水平垂直xyxyz测试1std.5.7233949330.090.060.1040.0790.0280.0300.019gps7.04467670.100.090.1370.1050.0300.0340.02pos0.1230.1120.09362.840.732.4测试2std.4.93969116860.060.060.0870.1280.0970.1040.03922.624.44.1gps6.74415115430.1530.0800.1380.06821.325.812.1pos0.2240.2940.16553.243.945.1测试3gps7.64430300.740.761.0610.5030.2030.2400.2329.710.99.9pos1.0641.4141.78135.931.831.9测试4std.7.6151910141.301.291.8301.4540.9200.9460.67518.818.56.9gps7.04425251.552.332.7981.2750.8780.9880.65817.919.66.5pos1.3242.8492.81761.857.167.3注：1）std.，gps和pos取于通过边缘密集gcps标准的方法获得的外方位元素，各自的gps-支持与四种正式的基础gcps和pos系统（与下表相同）。 2）因为gcps的分配不满足标准的要求，我们不能计算出边缘密集gcps大型区域网平差（与下表相同）。 3）控制点残差的rms计算来源于坐标平差和n个控制点的地面测量坐标之间的误差，即；。 4）eo的理论精度从根据误差的增值规律计算的未知量的分散矩阵获得，。1)对于测试1，密集点的水平精度和垂直精度优于0.15米，完全满足精度要求：平面精度0.25米和垂直精度0.30米的1:500地形图规格的空中摄影测量办公操作平坦土地。2) 对于测试2，密集点的水平精度优于0.15米，垂直精度优于0.20米，完全满足精度要求：平面精度0.35米和垂直精度0.40米1:500地形图规格的空中摄影测量办公操作的山地。3) 对于测试3，密集点的水平精度优于1.1米，垂直精度优于1.0米，完全满足精度要求:平面精度2.5米和垂直精度2.0米的1:5000地形图规格的空中摄影测量办公操作的山地。4) 对于测试4，密集点的水平精度优于3.0米，垂直精度优于1.5米,完全满足精度要求：平面精度17.5米和垂直精度3.0米的1:50000地形图规格的空中摄影测量办公操作的高山区的土地。从表2可以看出，在图像的不同土地类型不同规模下，从标准 at和gps-支持at获得的致密点均满意的四维产品的需求，而且获得外方位元素的精度的这两种方法通常是相似的。相片的比例越大，我们获得的相片的线性元素的精度就越高；但是相片的角元素的精度与相片比例无关，而与相机的焦距有关，因此焦距越短，相片的角元素的精度越高。此外，通过比较表2的数据，可以发现，通过pos系统提供的外方位元素的性能不如解析空中三角测量，也明显低于pos的理论精度 ，。3.3 直接地理参考的精度目前，四维产品普遍使用通过空中三角测量获得的致密点作为模型定向点取代直接从图像的外方位元素重建的立体模型。因此，外方位元素的精度要求没有定义在当前的规范。一般来说，只要我们可以获得足够的致密点满足每个模型的误差极限，通过执行绝对取向，可以重建一个可测量的模型。然后我们可以获取令人满意空间信息。通过使用不同的外方位元素，我们分析了直接地理参考的精度，首先使用前方交会的计算方法，通过使用在上述部分得到的不同精度的六个外方位元素，计算地面坐标，并且通过比较地面坐标与多数地面控制点的地面坐标，计算目标点地面坐标的精度(详细参见表3)。表3 直接地理参考的精度图像获得eq方法控制点数量最大差值/m最小差值/mrms/mxyxy高程xyxy高程xyxy高程测试1std.1880.18-0.250.295-0.2760.000.000.0090.0000.080.070.1070.081gps196-0.28-0.280.3340.2540.000.000.014-0.0010.130.080.1560.100pos196-0.310.390.3900.3100.000.010.0550.000100.127测试2std.412-0.22-0.140.222-0.3170.000.000.007-0.0010.060.050.0830.121gps419-0.290.210.322-0.4280.000.000.0060.0000.110.060.1310.166pos4190.26-0.490.4970.4870.000.000.0030.0010.090.240.2570.182测试3gps682.341.922.540-2.6600.020.020.146-0.0301.000.831.2991.325pos64-2.95-2.433.3992.330-0.030.040.2980.0821.200.781.4351.051测试4std.46-3.633.013.8844.374-0.030.050.1740.0281.401.331.9301.880gps464.434.754.8324.3590.000.020.1080.0331.451.502.2151.912pos465.294.475.4475.947-0.14-0.020.584-0.0582.551.793.1153.711从表3，这样的结论可归纳如下：1) 对于测试1，直接地理参考的水平精度优于0.25米，垂直精度优于0.15米，这也符合要求指明的1:500的平坦的土地的地形测绘标准：平面0.3米，海拔0.2米。2) 对于测试2，直接地理参考的水平精度优于0.3米，垂直精度优于0.2米，这也符合要求指明的1:500山区地形测绘标准：平面0.4米，海拔0.5米。3) 对于测试3，直接地理参考的水平精度和垂直精度均优于1.5米，这也符合要求指明的1:5000地形测绘标准：平面3.75米，海拔2.5米。4) 对于测试4，直接地理参考的水平精度优于3.15米，垂直精度优于3.75米，这也符合要求指明的1:50000的丘陵地形测绘标准：平面25.0米，海拔4.0米。从表3可以了解到，当直接地理参考应用到图像的不同土地类型和规模时，空中三角测量的方法优于利用pos获得的外方位元素。这个结果表明将通过空中三角测量的方法取得的外方位元素运用到直接地理参考中可以得到满意的地形测量的精度标准。因此可以推断出，只要空中三角测量符合精度标准，来源于此的外方位元素的内容绝对可以用于生产四维产品。3.4 重建立体模型中的y-视差在摄影测量工作中，我们应该注意的另一问题是利用外方位元素直接重建地形测绘的立体模型的可行性，即模型点的y-视差不超过20米。因此描述不同地形类型的三个立体像对利用得到的外方位元素从四个重建立体模型测试中分别获得，表4是y-视差在立体模型中的对应点。从表4可以看出，当使用空中三角测量的方法得到的外方位元素重建立体模型时，无论何种类型的地形和应用多大的比例，该模型点的最大y-视差位于一个像素，精度的最大y-视差低于半个像素，每个模型的能满足精度要求的地形测量垂直视差的有效值不应该超过20米。但是，当我们使用pos提供的外方位元素时，每个模型点的垂直视差都偏大，较小比例的相片，较大的垂直视差，完全不符合要求。表4 利用外方位元素的立体模型重构的y-视差图像模型pts最大高程不同 /my-视差/最大值最小值rms规定gpspos规定gpspos规定gpspos测试1282/281230.578.08.812.9274/2733217.8711.214.512.9343/3442736.3813.115.718.0测试214/132626.4115.011.718.129/284768.2710.917.819.60.07.210.3275/27434105.7414.310.012.8测试31017/10183126.0018.744.91.04.810.022.91013/10142291.6017.57.917.8234/23354192.9018.08.212.9238/23723375. 7617.08.117.7测试41076/1075327.4818.71.80.612.312.921.31115/11142817.4514.817.239.06.97.322.01105/110428109.5511.612.91.06.56.314.04 结论从实验可以看出，如果从空中三角测量获得的外方位元素符合精度标准，它们可以直接用于图像的定向和立体模型重建。然而，由于存在pos外方位元素的系统误差，当前很难符合摄影测量的标准，特别是当提取三维空间信息时。发现，在数字摄影测量时期，工作可以由计算机自动地完成，对地面控制点的信任逐渐减少，从而简化了摄影测量的操作。从总体上看，标准，是获取图像定向参数最确定的和广泛使用的方法，仍然是摄影测量的主体；gps-支持at是易于操作和低成本的方法，相应的标准已经被起草：pos直接地理参考在摄影测量中是一个重要的尖端技术。基本空间信息采集应利用这一点，设计良好的计划获得最大的经济利益。我们建议，对于交通状况良好的平坦地区的大比例测图，应该主要采用标准空中三角测量；对于困难区域，非指明的地区或领域无法访问，没有地面控制点的gps-支持at可以为国民生产的基础地图获取基本的空间信息；pos摄影测量可用于在小区域生产正射影像和更新四维产品。然而，pos在大比例的城市测图，lidar，数字航空摄影测量的领域中是有前景的。我们应该通过开展大规模的实验推动pos系统的整合技术及其传感器，从而为经济和迅速聚集的地理空间信息提供技术支持。鸣谢数据收集的实验由遥感应用中国科学院，中非通用航空有限公司，辽宁经纬测绘科技有限公司，大连市测绘设计研究院，四维航空遥感有限公司，西安全国测绘航空遥感有限公司等支持。对这些支持表示感谢。作者希望对参与部分实验的付建宏，谢筹，季顺平和杨明表现出他热诚的谢意。作者要感谢杨明和张井雄教授为他们润饰英语。参考文献1 李徳仁，单杰（1989）导航数据大型区域网平差的质量分析j。工程摄影测量与遥感，55（12）：1743-17462 ackermann f（1994) practical experience with gps-supported aerial triangulation jphotogrammetric record, 16(84):861-8743 袁修孝（2001）gps-支持空中三角测量的原理及应用m。北京：测绘出版社（中文）4 cannon m e, sun h (1996) experimental assessment of a non-dedicate gps receiver system for airborne attitude determination jisprs journal of photogrammetry & remote sensing, 51(2): 99-1085 cramer m, stallmann d, haala n (2000) direct georeferencing using gps/inertial exterior orientations for photogrammetric applications c. international archives of photogrammetry and remote sensing, 33(b3):198-2056 gruen a, baer s (2001）aerial mobile mapping-georeferencing without gps/ins cthe 3rd intemational symposium on mobile mapping technology, cairo7 jacobsen k (2002) calibration aspects in direct georeferencing of frame imagery cpecora 15/land satellite information, isprs commission i/iv8 mostafa m r (2002) camera/imu boresight calibration:new advances and performance analysis c.the asprs annual meeting, washington d.c.9 袁修孝，谢筹，王书亘(2004)通过与固定图像组合的大型区域网平差进行图像定向j。摄影测量与遥感的国际文件，35(b2)：34-3810 applanix产品概述ol。/products/ pospac_airborne_index.php11 国家统计局的标准（1998）1:500，1:1000，1:2000空中摄影测量办公操作的地形图规格s。北京：中国标准报(中文)12 国家统计局技术监督（1993）1:5000，1:10000航摄相片办公操作的地形图规格s。北京：中国标准报(中文)13 国家统计局技术监督（1991）1：25000，1:50000，1:100 000航摄相片办公操作的地形图规格s。北京：中国标准报(中文)14 王之卓（1990）摄影测量原理（遥感）m。武汉：武汉测绘学院学报；北京：测绘出版社15 袁修孝，付建宏，阎佩利等(2007)通过使用从定位和定向系统获得的外方位元素分析立体模型重建的垂直视差j。武汉大学测绘和信息科学，32（6）：489-493 on stereo model reconstitution in aerial photogrammetryyuan xiuxiaoabstract this paper describes the operational issues and basic technical requirements of modern aerial photogrammetry.the accuracy of photogrammetric point determination and the y-parallax at corresponding model points is analyzed when stereo models are reconstituted by using the exterior orientation elements of aerial images. real aerial photographs, at image scales from 1:2 500 to 1:6 0000，with dgps/imu data taken from various topographies in china were processed by our pos-supported bundle block adjustment program wucaps. the empirical results verified that the accuracy of the exterior orientation elements from bundle block adjustment meets the requirements of the specifications of topographic mapping.however, the accuracy of the exterior orientation elements determined by pos fails to meet the requirements of the specifications of topographic mapping.keywords aerial triangulation (at); gps (global positioning system）；pos（position and orientation system）；stereo model reconstitution；ground control points (gcps)；accuracyintroduction aerial photogrammetry is the science and technology for obtaining 3-dimensional spatial information about the earths surface from aerial imagesphotogrammetric point determination, which locates ground objects by using images, is the basis for object recognition in remote sensing. and the key point of this issue is the rapid and accurate determination of an images position and behavior at the instant of imaging. this goal was met by aerial triangulation based on well distributed gcps. with the development of spatial positioning technology, remote sensing technology, and computer science, aerial triangulation evolved and progressed towards digital mapping without gcpsin the early 1950s, photogrammetric scientists began studying how to utilize various auxiliary data to reduce the number of gcps required. however, the methods havent become practical due to technological limitationsuntil 1970s, with the emergence of american global positioning system (gps), people got to adopt carrier phase differential gps（dgps）technology to determine an exposure stations positions(that is three linear elements of aerial photos) during aerial photographic process, which was used to perform aerial triangulation (called cjps-supported at for short) that can decrease photogrammetric reliance on gcps, shorten the mapping cycle; and reduce production costs, triggering the revolution in the field of photogrammetry.nevertheless, gps-supported at is advantageous for aerial photogrammetric operation primarily over vast and difficult areas, at small and medium mapping scales, not for strip-like zone and urban large-scale mappingin the 1990s,people started to investigate employing gps/lns integrated system (also called pos) to acquire a photos position and attitude (i.e.，to obtain exposure stations position by gps，and imagesattitude elements by imu), for the purpose of photo orientation, and the final goal is to replace block aero triangulation procedure modern digital photogrammetry will play an important role in automated productions of 4d products(dem, dom, dlq drg) and updating of spatial databases. this paper will introduce current operational applications of aerial photogrammetry and related technical requirements, in particular, geometric positioning accuracy obtainable in the photogrammetric information chain from photo orientation to stereo-model reconstitution, aiming to investigate their practicability for 4d products production. it is hoped that findings from this study will provide guidance for operational aerial photogrammetry in the context of national, land surveying, mapping, and fundamental geographic information acquisition.1 current patterns of the modern aerial photogrammetry nowadays, there are primarily three patterns for aerial photogrammetry, namely, standard aerial photogrammetry gps-supported aerial photogrammetry and pos-supported aerial photogrammetry. their main procedures are shown as fig. l. from fig.l，we can learn that the main difference between these three patterns lies in the ways of aerial photo acquisition and photo orientation. for standard at, it is through block aerotriangulation with a large number of gcps to get a model orientation points coordinates to complete image orientation. for gps-supported at, in aerial photo acquisition process,dynamic gps positioning is used instead of gcps to determine the positions of exposure center and meanwhile obtain the models orientation pointscoordinates, which are then used to rectify the images orientation. for pos-supported at, images and their corresponding orientation elements (six exterior orientation elements of images) are both acquired, in order to realize geometric inversion of photography by storing their spatial positions and attitude at the moment of exposure2 related technologic requirements2 .1 aerial photography in modern aerial photography, in order to improve the quality of obtained images, besides adding flight control systems to aerial camera (such as ascot,ccns4, aerial tracker system), the methods include sticking a gps receiver with the camera firmly when adopting gps aerial photography and mounting pos system on the camera in dgps/imu aerial photography. according to the different patterns of aerial photography, we can formulate an answerable plan as shown in fig.2.2.2 ground control plan in digital photogrammetry workstations, aerotriangulation is carried out by the most theoretically rigorous procedure of bundle block adjustment, but for the sake of obtaining the best pass points coordinates and the exterior orientation elements of photos,ground control plan should be designed, as shown in fig.3，for different patterns of aerial photogrammetry.2.3 digital mapping theoretically, after getting the accurate exterior orientation elements of images, measurable stereo models can be reconstructed using model restoration,by which we can do surveying and mapping of terrain and objects automatically. however, the current process of producing 4d product is：single photo interior orientationrelative orientation of stereo pairsingle model absolute orientationsurveying and mapping on stereo modelsthe method of model restoration is only adopted in the direct georeferencing of pos-supported aerial photogrammetry.3 experiments and analysis there are two ways of aerial photogrammetric positioning. one is called block aerotriangulation, regarding image pointscoordinates, gcpscoordinates and/or the exterior orientation elements of images as weighted observed values, and combined bundle block adjustment is performed to solve the images orientation parameters and target pointsspatial coordinates, so as to supply orientation control points for stereo model mapping and