摄影测量后方交会、近景摄影测量DLT直接线性变换(代码下载见介绍)

近景摄影测量实习报告班级：07021班学号：200732590062姓名：方毅日期：2010年12月1日第一部分摄影、像点量测以及DLT和单像后方空间交会解算实习的目的和意义1.1数码相机摄影：为后续的摄影测量解析处理提供质量合格的数字影像，了解所用数码相机的特点及使用，学习获取数字影像的方法。1.2像点量测：量测所拍摄的高精度室内三维控制场中控制点的像点坐标，为后续摄影测量解析处理准备计算数据。1.3直接线性变换（DLT）与单像空间后方交会解算：加深理解近景摄影测量直接线性变换与单像空间后方交会的理论，学习准备数据和调试程序的方法。2.实习原理2.1DLT直接线性变换直接线性变换解法是建立像点坐标仪坐标和相应物点物方空间坐标直接的线性关系的算法。它的基本关系式如下： （2.1）展开可得到以为未知数的方程： （2.2）当有个控制点时，即可列出个方程式，写为矩阵的形式如下： （2.3） 即。其中阵为由像点坐标和物方坐标构成的大小的系数矩阵，阵为由构成的大小为的未知数矩阵。为由像点坐标构成的大小为的观测值矩阵。根据最小二乘法，可求得未知数矩阵为。由外方位元素的线元素满足 （2.4）根据公式（2.4）解方程组可得外方位元素的线元素。根据可求得相片的旋转矩阵，继而根据公式（2.5）可求得外方位元素的角元素 （2.5）2.2单片空间后方交会单片空间后方交会是基于共线条件方程的： （2.6）将（2.6）式看成外方位元素的函数时，方程（2.6）是非线性的。将像点坐标视为观测值，将（2.6）式线性化并展开可得误差方程的矩阵形式为（2.7）当内方位元素已知时，可视。则（2.7）式子可变形为 （2.8）若已知个控制点，可列出个方程式，写成总误差方程为：其中： （2.9）要解6个未知数，至少需要3个控制点。由于后方交会求解是非线性的，故求解过程需要迭代。3.实习过程和步骤3.1数码相机摄影与像点量测使用柯达NIKOND2Xs型单反数码相机，配Nikon28mm定焦镜头，CCD芯片大小为23.584mm*15.664mm，影像分辨率设定为4288*2848，像元大小为对遥感信息工程学院近景摄影实验室内的高精度三维控制场进行摄影，建议使用28mm焦距，以便使视场较大，拍摄范围广。选择合适的摄站位置，按交向摄影方式或正直摄影方式拍摄室内高精度三维控制场，另外还要注意不要使吊杆互相遮挡。使用Photoshop软件。采用数码相机所拍摄的控制场影像，使用Photoshop的坐标量测功能，人工逐点量测30个同名控制点的影像坐标，并记录。步骤如下：（1）运行Photoshop软件，打开要量测的影像；（2）点击菜单项Window，选中Info（信息）选项；右键点击显示坐标的“+”符号，选中Pixel（像素）选项。此时，量测单位已选定为像素。（3）将鼠标左键移动到要量测的点上去，在Info（信息）选项中读出该点的像点坐标，并记录下点的编号和坐标。（4）依上步骤在左右片上量测30个同名点的坐标，按照相应格式保存为文档文件。3.2编写算法程序 按照算法原理，利用VC6.0平台编写解算程序。程序可以按照文件格式导入数据。编写程序界面和功能如图所示：图3.1 左图为主界面图、右图为导入数据后界面图图3.2 DLT解算结果图图3.3 后方交会结果图图3.4 前方交会结果与误差图图3.5 批量选取组合点进行处理运算界面图4.实验结果与分析评价在左右片量测40个同名像点的像素坐标后，运行软件选择“数据准备”→“导入数据”，依次导入数据文件。选择“数据解算”，依次进行“DLT解算”、“后方交会”、“前方交会”。运算结果如下：DLT直接解算结果：x0y0fxXsYsZsphiomegakapa

0.0281-0.056128.44121314.6669-115.6885483.35170.1055-0.0849-0.0866

0.0047-0.036128.43533782.6437-114.5162682.3636-0.1684-0.0832-0.0872 后方交会结果：*****后方交会外方位元素及单位权中误差************************************************左片**Xs=1314.637200±0.517828mm**Ys=-114.766504±0.572956mm**Zs=481.954069±0.315933mm**Φ=0.104550±0.000087rad**ω=0.086731±0.000096rad**κ=0.005592±0.000049rad**σ=1.536274e-003****右片**Xs=3782.915294±0.745427mm**Ys=-114.039851±0.836396mm**Zs=681.804932±0.385380mm**Φ=-0.168606±0.000116rad**ω=0.084374±0.000129rad**κ=0.003071±0.000053rad**σ=1.843557e-003*************************************************************************************方毅200732590062前方交会结果：点名X'Y'Z'△X△Y△Z4203603.52621385.2196-4933.7788-2.7968-5.20112.69654204603.2884936.1517-4933.4445-2.6344-4.85522.96114205602.6082479.8230-4932.3016-2.7157-4.54892.49454206602.7378196.4443-4931.2174-2.6540-4.17282160319.1461-4932.9895-2.6028-4.51913.75084208602.0186-184.6048-4931.1170-3.2669-4.35062.27564209602.4185-387.1557-4931.9077-3.0084-4.29422.99574210602.7104-625.5547-4931.3964-3.1973-4.43572.35522203602.7262927.7681-5934.7715-3.6001-5.06890.54752204602.8841629.4455-5934.9909-3.5174-5.09421.11802205601.9474289.8247-5935.9648-3.5225-4.88481.80662206602.6785-294.3308-5935.8088-3.5074-4.90062.031643031618.17971320.2521-4917.8507-3.5393-5.18631.126943041619.15841005.8658-4917.6622-3.3748-4.50810.971543051620.6429694.2870-4916.7398-3.4457-4.2469-0.329943061619.1293391.7027-4917.4459-3.1962-4.21070.211543071621.9768149.8209-4918.1213-3.0310-4.32990.314328022155.49931667.3320-5924.3092-5.0200-6.25522.708728032153.97631460.9532-5920.4991-4.8706-5.89572.743328042157.43651210.0222-5921.7221-4.9185-5.44270.404728052154.3354963.5929-5922.0160-4.9021-5.06900.673928062158.2612663.8233-5921.4079-4.8503-4.7973-0.059228072155.1431441.4144-5921.1157-4.5074-4.93520.184838042142.69601135.6523-5424.1222-3.8762-5.0298-0.728238052144.9852844.7984-5423.0634-3.8217-4.3598-1.287438062139.6365528.2201-5424.3787-4.0491-4.2408-1.614638072142.2686210.8295-5423.0929-4.1699-4.3258-0.799238082141.151468.0598-5424.0179-3.7598-4.47671.418534042668.6369917.4369-5416.0460-4.0835-4.5503-1.885634052668.3196475.6684-5418.5402-4.1770-4.6853-0.453334062667.2598174.2287-5418.1364-4.2007-4.3354-0.305034072666.2516-125.7956-5416.8454-4.0372-4.39601.150634082663.5026-425.9057-5418.2525-3.8805-4.35372.347725053697.8780513.8848-5896.6424-4.2134-4.6980-0.845025073696.0782-109.4291-5897.9393-4.1111-5.09000.060845043690.39691044.9818-4901.0362-3.6872-4.76460.193845063691.3125382.4256-4900.3078-3.7024-4.29860.212116024746.64721006.0282-6993.2420-3.3256-5.4945-1.848915013719.97971876.1954-7005.5435-5.3158-5.8617-1.812715033718.291222.4742-7007.5606-4.7525-5.5608-0.4808根据运算结果可见，后方交会中左右片外方位元素的线元素中误差大小均小于1mm，角元素中误差小于0.0002rad。基本满足误差要求。通过前方交会反算对比原控制点物方坐标，可得、不超过6mm，不超过3mm。程序中还有一项批处理计算数据的功能，其原理是在已知量测的个控制点中，任选个控制点进行DLT——后交运算。即算出个控制点任选个点的全部组合方式，将进行次运算。通过观察批处理运算的结果可以发现，左右片外方位元素计算结果与点位分布以及刺点精度相关。当某一点为粗差点时，后交所得的线元素误差可能大于1mm。而当选点点位分布不具有层次性（仅在一个或几个深度平面中选点），后交所得精度将受到严重影像。总得来说，要提高精度，选点点位分布应该①控制点在空间分布均匀，具有一定深度层次性；②在像片上的构象范围大。控制点不能不设在空间任意同一平面内，引用为控制点布设在空间的起伏不大，则解算结果也不稳定。同时对控制点像素坐标量测时，应该尽可能地瞄准控制点刺点中心。若某一控制点像素信息较为模糊，应舍去该点重新选点。5.实习心得与体会在本次实习的过程中，我了解到了近景摄影测量与传统摄影测量的异同。它们的联系是联系：①基本原理相同，即我们可以用相同的原理实现算法设计；②模拟处理方法、解析处理方法、数字影像处理方法相同，即量测图片像点坐标，完成像点坐标到像素坐标的转换等处理方法是相同的；他们的区别：①测量的目的不同；②被测量目标物不同；③目标物纵深尺寸与摄影距离比不同，一个明显的区别是近景摄影测量的相片倾角比较大；④摄影方式不同；⑤影像获取设备不同；⑥控制方式不同；⑦近景摄影测量适合动态目标；由于传统航空摄影大致符合竖直摄影，平均航高和摄影比例尺都知道，摄影相片倾角较小，所以其外方位元素的初始值可以比较近似的确定。而近景摄影测量由于其摄影方式决定了摄影相片倾角较大，摄影距离较小，因此观测点坐标应该包含畸变差改正，同时其外方位元素无法用任意初值确定。必须要用近似值，否则后方交会解算时无法迭代。直接线性变换能够在不需要外方位元素的初始值的情况下求出像点的仪器坐标与物方坐标的对应关系也即外方位元素的值，因此DLT直接解算可以提供后方交会的初始值。在具体编写程序代码过程中，前后花了4天的时间，共完成了4个版本的程序。前2个版本的程序仅能实现实习的要求，在具体操作上、人机交互上略显复杂。因此后续2个版本程序重新设计了程序框架与算法格式，将程序的数据与具体算法完全分开。通过定义的接口调用函数实现相应功能。重新设计的程序框架极大地提高了函数的重用率，并且引入线程机制，在大量重复运算过程中，程序不会过多征用CPU资源而导致假死。程序的稳定性也大大提高。同时，线程函数的引入可以控制运算过程，可以随时中断或暂停运算。4天的编程实习，收获颇多。本次实习进一步理解了摄影测量的基本原理，同时强化了程序设计的方法与流程，希望以后再接再厉！第二部分“南海观音”模型摄影测量内业处理。实习的目的和意义1）学习使用VirtuoZo软件处理近景摄影测量立体像对；2）学习摄影测量处理时P31相机影像内定向的方法；3）学习相对定向过程中人工加点、人工剔除点的方法；实习原理建立模型—设置相机参数—半自动内定向—人工选点相对定向—半自动绝对定向—生成非水平核线影像—影像匹配—匹配点的编辑—生成DEM—生成正射影像实习流程打开VirtuoZo软件建立测区文件、控制点文件、加密点文件、相机文件设置—相机参数像主点x0:0.004mm像主点y0:-0.011mm焦距：100.64mm+4+0+4+0+31++21）主距（mm）：99.582）主点（mm）：0.011 0.0043）主距改正值（mm）：1.064）框标理论坐标（mm）： 0 -0.002 0.0001 57.497 0.0002 -0.001 -27.5013 -57.496 0.0004 0.000 57.49857.498+10++10++42++3+3+0+24++1图3-1影像格式转换文件-引入-影像文件界面中像素大小为：0.02mm注意：1）左影像转换用左相机文件，右影像转换用右相机文件；2）输出路径名改在自己的工作目录下的images文件夹中。建立模型文件-打开模型内定向1