ICESat激光高程点辅助的天绘一号卫星影像立体区域网平差

无地面控制点（简称无控）区域网平差是实现卫星影像无控测图的一项重要技术,对于境外和外业测控困难区域的测图具有重要意义。然而,无控区域网平差的定位精度一般难以满足对应比例尺测图规范要求。利用公开、可稳定获取的公众地理信息数据辅助区域网平差,是提高卫星影像无控定位精度的有效途径,其中ICESat激光高程点便是一种良好的高程控制数据。为了提高天绘一号卫星影像无控定位精度,本文提出ICESat激光高程点辅助的卫星影像模型法立体区域网平差方法。首先,以30 m分辨率SRTM估算的地形坡度作为限制条件,结合激光高程点自身质量评价信息,自动提取高质量ICESat激光高程点;其次,利用自动匹配的连接点进行模型法自由网平差,实现卫星影像几何定位精度的相对一致性（内部一致性）;最后,将激光高程点自动量测至卫星影像作为控制点,其平面坐标根据自由网平差结果前方交会计算而得,高程坐标取自激光点高程,再次进行区域网平差精化定向参数,提高卫星影像的绝对高程精度。最后本文利用山东全省的天绘一号卫星影像进行试验,验证了本文方法的有效性和可行性。     

基于已知定向参数影像的光束法区域网平差

从自检校光束法区域网平差模型出发,建立了利用同一地区已知定向参数的影像解求新获取影像外方位元素并进行目标定位的数学模型。以此为基础,对来自某地区相隔两年的三期不同摄影比例尺航空影像两两组合地进行联合光束法区域网平差。结果表明,当两期影像比例尺相近时,所解求的新影像外方位元素和加密点坐标精度基本达到常规光束法区域网平差的精度;当两期影像比例尺不同时,利用已知定向参数的大比例尺影像解求出的小比例尺影像的外方位元素和加密点坐标可满足现行规范精度要求,利用已知定向参数的小比例尺影像去解求大比例尺影像的外方位元素和加密点坐标精度明显下降且不可用。本文的研究证实,就同一地区的多时相航摄影像而言,利用已知定向参数的影像进行新影像的定向方法是可行的。     

附加视线向量修正的卫星影像区域网平差

针对传统有理函数模型(RFM)区域网平差方法局限于姿态和轨道测量误差小、相机视场角小及影像交会角良好的情况,提出了附加视线向量修正的卫星影像区域网平差方法。首先利用影像附带的有理多项式系数(RPC)计算出像元视线向量,其次根据该视线向量恢复成像时刻虚拟位置和姿态信息,然后对恢复的虚拟位置和姿态构建误差补偿模型,最后通过最小二乘方法整体解算模型参数和连接点物方坐标。该方法从系统误差产生的原因构建补偿模型,可以规避传统区域网平差方法的近似假设和条件限制。通过对模拟数据以及多套测绘卫星和非测绘卫星数据进行试验的结果表明,该方法处理大姿态角误差、大视场角以及弱交会角等各种严苛条件下的卫星影像能达到比传统方法更好的效果。     

激光测高数据辅助卫星成像几何模型精化处理

资源三号02星搭载了中国首个对地观测试验性激光测高载荷。借鉴目前较成熟的卫星影像区域网平差理论的基础上,结合近年来激光测高数据精度的大幅提升以及资源三号02星激光测高数据的特点,首次提出了激光测高数据辅助卫星立体影像进行成像几何模型精化处理的通用理论。首先,利用传统的区域网平差算法对所处理影像进行高精度连接点匹配处理,并对其进行无约束的自由网平差处理,获得高精度相对精度及不亚于原始成像几何模型的绝对精度;其次,根据激光测高数据3维坐标和精化后参考影像成像几何模型获取激光数据参考影像坐标;而后将参考影像坐标通过几何模型映射获取目标影像上待匹配影像坐标,通过连接点匹配算法,对待匹配目标影像坐标进行精化获取高精度像方同名点;最后,以同名点作为高程控制进行区域网平差计算,对影像成像几何模型进一步处理,获取高精度补偿参数。通过湖北、青海两测区的试验,以激光测高数据辅助卫星影像几何模型精化精度可分别达到1.97 m、3.23 m,结果表明本文提出的方法可有效提高卫星立体数据测图精度。     

Complete Automation of Digital Aerial Triangulation

Automatic digital aerial triangulation is reviewed. It depends essentially on the automatic selection and transfer of tie points by digital image processing. The MATCH-AT system is briefly sketched. Its basic design parameters are feature based matching, use of tie point clusters and integrated block adjustment. The system runs almost completely in batch mode, with as little interactive interference as possible. Some test results are discussed which confirm the general expectation that automatic digital aerial triangulation is faster, more precise and less costly than conventional aerial triangulation.     

Solution of the weighted symmetric similarity transformations based on quaternions

A new method through Gauss–Helmert model of adjustment is presented for the solution of the similarity transformations, either 3D or 2D, in the frame of errors-in-variables (EIV) model. EIV model assumes that all the variables in the mathematical model are contaminated by random errors. Total least squares estimation technique may be used to solve the EIV model. Accounting for the heteroscedastic uncertainty both in the target and the source coordinates, that is the more common and general case in practice, leads to a more realistic estimation of the transformation parameters. The presented algorithm can handle the heteroscedastic transformation problems, i.e., positions of the both target and the source points may have full covariance matrices. Therefore, there is no limitation such as the isotropic or the homogenous accuracy for the reference point coordinates. The developed algorithm takes the advantage of the quaternion definition which uniquely represents a 3D rotation matrix. The transformation parameters: scale, translations, and the quaternion (so that the rotation matrix) along with their covariances, are iteratively estimated with rapid convergence. Moreover, prior least squares (LS) estimation of the unknown transformation parameters is not required to start the iterations. We also show that the developed method can also be used to estimate the 2D similarity transformation parameters by simply treating the problem as a 3D transformation problem with zero (0) values assigned for the z-components of both target and source points. The efficiency of the new algorithm is presented with the numerical examples and comparisons with the results of the previous studies which use the same data set. Simulation experiments for the evaluation and comparison of the proposed and the conventional weighted LS (WLS) method is also presented.     

From point-based to feature-based aerial triangulation

This paper is concerned with using linear features in aerial triangulation. Without loss of generality, the focus is on straight lines with the attempt to treat tie lines in the same fashion as tie points. The parameters of tie lines appear in the block adjustment like the tie points do. This requires a unique representation of lines in object space. We propose a four-parameter representation that also offers a meaningful stochastic interpretation of the line parameters. The proposed line representation lends itself to a parameterized form, allowing use of the collinearity model for expressing orientation and tie line parameters as a function of points measured on image lines. The paper describes in detail the derivation of the extended collinearity model and discusses the advantages of this new approach compared to the standard coplanarity model that is used in line photogrammetry. The intention of the paper is to make a contribution to feature-based aerial triangulation on the algorithmic level.     

影像连接点均衡化高精度自动提取

针对空中三角测量影像连接点提取中存在的误匹配、点位多而分布不均和点位定位精度低等问题,提出了一种影像连接点均衡化高精度自动提取方法。首先采用分块SIFT技术进行特征的提取与匹配,并利用并查集数据结构进行特征点的高效多视追踪;然后采用提出的物方分块点位筛选算法对点位进行了相对均衡化的择优挑选;最后采用最小二乘匹配技术对得到的SIFT连接点坐标位置进行精化。选取中国嵩山遥感定标场的有人机影像和沙漠地区无人机影像作为试验数据,通过目视检查、像方反投影误差和检查点精度等3个指标进行了分析,结果表明本文方法有效克服了弱纹理和重复纹理导致的连接点提取与匹配困难,并改善了连接点分布的均匀性和提高了连接点的定位精度。     

基于附加参数模型的自检校光束法区域网平差算法运行结果对比

像点偏移是影响光束法区域网平差结果精度的主要因素之一。克服平差过程中像点偏移的影响,以提高区域网平差精度,一直受到业界高度重视。当前,克服像点偏移的常规做法是借助特定的参数模型对其进行描述,构建基于附加参数模型的自检校光束法区域网平差模型,在区域网平差的同时求解引入的参数,进而实现对系统误差的补偿。Bauer模型、Ebner模型、Brown模型、光学畸变模型是现有算法中描述像点偏移的常用模型。为了对上述4种模型的系统误差进行验证,本文首先在经典的光束法区域网平差模型基础上分别设计了与4种模型相对应的自检校光束法区域网平差模型,并予以编程实现;然后通过实例对上述4种模型的算法运行结果进行了验证、比较及分析。试验结果表明,Ebner模型和Brown模型的引入可有效提高光束法区域网平差的精度;Bauer模型的引入对小型区域网平差精度的提升更加有效;光学畸变模型的引入对区域网平差结果的精度提升作用相对较小。     

THE USE OF AIRBORNE AUXILIARY DATA IN THE RIGOROUS LEAST SQUARES ADJUSTMENT OF A BLOCK OF AERIAL TRIANGULATION

The problem of incorporating airborne auxiliary data into an aerial triangulation block adjustment is considered. The conclusion is reached that this problem can be most readily solved by the use of a fully rigorous adjustment method. The first order adjustment equations for all types of airborne auxiliary control are developed and the method of incorporating these equations into the block solution is demonstrated. The general computing problems of the aerial triangulation block adjustment are briefly discussed and a method outlined for the incorporation of tie strips, at any angle across the block, into the adjustment without unduly increasing the storage requirements for the solution. The problem of correctly weighting the different types of available data is discussed. The paper concludes with a demonstration of the improved accuracy possible if precise aerial triangulation data are combined with apparently crude absolute airborne data.     

Nonlinear analysis of the three-dimensional datum transformation [conformal group ℂ7(3)]

 The weighted Procrustes algorithm is presented as a very effective tool for solving the three-dimensional datum transformation problem. In particular, the weighted Procrustes algorithm does not require any initial datum parameters for linearization or any iteration procedure. As a closed-form algorithm it only requires the values of Cartesian coordinates in both systems of reference. Where there is some prior information about the variance–covariance matrix of the two sets of Cartesian coordinates, also called pseudo-observations, the weighted Procrustes algorithm is able to incorporate such a quality property of the input data by means of a proper choice of weight matrix. Such a choice is based on a properly designed criterion matrix which is discussed in detail. Thanks to the weighted Procrustes algorithm, the problem of incorporating the stochasticity measures of both systems of coordinates involved in the seven parameter datum transformation problem [conformal group ℂ₇(3)] which is free of linearization and any iterative procedure can be considered to be solved. Illustrative examples are given. Received: 7 January 2002 / Accepted: 9 September 2002 Correspondence to: E. W. Grafarend     

GIS结点捕捉的广义算法及误差传播模型

根据最小二乘原理，本文提出了结点捕捉的一种广义算法，并建立了伴随的误差传播模型，针对位于模糊公差范围内待捕捉点组中各点坐标误差统计我的各种可能的特殊情况，进一步导出了相应的简化算法以及其误差传播模型，并从理论和数值模拟两方面系统地分析了结点坐标的相准确性一对捕捉结果的影响，理论推导表明，现有的GIS结点捕捉算法属广义算法的一种特殊情况，最后通过算例说明了广义算法和误差传播模型的实际应用前景。     

On symmetrical three-dimensional datum conversion

A 3-D similarity transformation is frequently used to convert GPS-WGS84-based coordinates to those in a local datum using a set of control points with coordinate values in both systems. In this application, the Gauss-Markov (GM) model is often employed to represent the problem, and a least-squares approach is used to compute the parameters within the mathematical model. However, the Gauss–Markov model considers the source coordinates in the data matrix (A) as fixed or error-free; this is an imprecise assumption since these coordinates are also measured quantities and include random errors. The errors-in-variables (EIV) model assumes that all the variables in the mathematical model are contaminated by random errors. This model may be solved using the relatively new total least-squares (TLS) estimation technique, introduced in 1980 by Golub and Van Loan. In this paper, the similarity transformation problem is analyzed with respect to the EIV model, and a novel algorithm is described to obtain the transformation parameters. It is proved that even with the EIV model, a closed form Procrustes approach can be employed to obtain the rotation matrix and translation parameters. The transformation scale may be calculated by solving the proper quadratic equation. A numerical example and a practical case study are presented to test this new algorithm and compare the EIV and the GM models.     

一种结合已有DOM和DEM的高精度自动刺像控点方法

针对基础地理数据更新需求，为了免去外业布设控制点，减少内业刺点工作量，提高控制点精度、密度，提出了一种基于已有DOM和DEM数据，利用IPS软件和辅助软件，实现高精度自动刺像控点的方法。将已有DOM根据具体匹配需求裁切，裁切后的DOM和像片同时加载到IPS中进行同名像点匹配，根据DOM的地理参考、匹配点像平面坐标、TFW头文件及DEM共同解算出匹配点的物方坐标。通过平差探测和目视检查剔除错误点，改变匹配点属性为控制点，即完成了控制点刺入工作。之后可在IPS中继续完成空三处理、地理信息产品制作等，也可导出、转换为其他软件的文件格式，供后续处理使用。     

Principle, software and experiment of GPS-supported aerotriangulation

1 IntroductionAs is now well known, the high accurate point de-termination with airborne remOe sensing data hasalways ben one of the most fundaxnental prObletns..in aerial photOgrammtry. According to the princi-ple of the geOmtry reversal in photOgramrntry,the interior and exterior orientation elements ofaerial phOtOgraphs must first be known in order toreconstnJct the measuring stereo geometric medels.For the past 60 years, however, the interior orienta-tion parameters of carnera were main…     

基于史赖伯规则的机载InSAR区域网平差

采用区域网平差的方法进行合成孔径雷达干涉测量（interferometric synthetic aperture radar,InSAR）干涉参数定标,是减小大区域、多个干涉像对重叠处反演高程差异的有效方法。为了提高InSAR区域网平差的答解效能,设计一种利用史赖伯规则的机载InSAR区域网平差干涉参数定标方法。该方法利用史赖伯规则,由连接点误差方程式构建等效误差方程式,可有效减小法方程系数矩阵大小,降低对计算机配置的要求,提高平差的答解效能。采用我国机载双天线InSAR数据进行干涉参数定标试验,改化前后两种方案,干涉参数定标结果没有明显变化,改化后的平差耗时均值减小了,验证了利用史赖伯规则的机载InSAR区域网平差干涉参数定标的正确性和有效性。     

数字航空摄影三维重建理论与技术发展综述

航空摄影测量作为摄影测量学最重要的分支之一,近年来得到了长足的发展。倾斜航空摄影和无人机摄影测量等多种新作业模式的出现,给传统航空摄影测量带来新的挑战的同时也催生出了诸多新的解决方案。此外,人工智能领域计算机视觉技术和深度学习技术中的新理论、新方法不断融入航空摄影测量中,推动航空摄影测量向智能化、自动化方向发展。当代航空摄影测量学已经是多种传感器融合、多种数据采集方式结合、传统摄影测量和人工智能技术交叉的产物。三维重建是航空摄影测量的核心问题之一。本文阐述了当代航空摄影三维重建技术的发展趋势和存在的问题,着重从航空影像的同名连接点自动提取与匹配、区域网平差、密集匹配和单体化建模4个方面对当前的研究现状进行了总结讨论,给出了当前国内外主流的航空影像摄影测量处理框架。     

控制网自动平差软件的研制

阐述控制网自动平差软件设计的基本思想,研究了点号的自动最优编排、近似座标的自动推算、变带宽压缩存储等关键算法,据此在微机上设计了一个自动平差软件。该软件能直接处理电子手薄、磁卡和人工手薄记录的数据,用户不需绘制计算略图、编号编码、数据摘录等,处理过程完全自动进行,使平差过程变得非常简单。软件可用于平差导线网约800个点,水准网约2000个点。     

Orientation of MOMS-02/D2 and MOMS-2P/PRIRODA imagery

This paper deals with the orientation of three-line imagery which has been taken during the MOMS-02/D2 experiment in spring 1993, and during the MOMS-2P/PRIRODA mission since April 1996. The reconstruction of the image orientation is based on a combined adjustment of the complete image, ground control, orbit and attitude information. The combined adjustment makes use of the orientation point approach or the orbital constraints approach. In the second case, the bundle adjustment algorithm is supplemented by a rigorous dynamical modeling of the spacecraft motion. The results of the combined adjustment using MOMS-02/D2 imagery and control information of orbit #75b are presented. Based on the orientation point approach an empirical height accuracy of up to 4 m (0.3 pixel) is obtained. In planimetry the empirical accuracy is limited to about 10 m (0.7 pixel), since the ground control points (GCP) and check points could not be identified in the imagery with the required accuracy. Computer simulations on MOMS-2P/PRIRODA image orientation based on realistic input information have shown that good accuracies of the estimated exterior orientation parameters and object point coordinates can be obtained either with a single strip and a few precise GCP or even without ground control information, if a block of several overlapping and crossing strips with high geometric strength (q≈60%) is adjusted.