2023摄影测量实习报告.docx

2023摄影测量实习报告_摄影测量实习报告 这篇摄影测量实习报告总结：测量实习报告是由小编整理供应的，请大家参考！ 摄影测量实习报告总结：测量实习报告 一、 实习任务 利用自己所熟识的一种编程语言，实现单像空间前方交会，解求此张像片的6个外方位元素 ， ， ， ， ，范文之实习报告:摄影测量实习报告。 二、 实习目的 1、 深刻理解单张像片空间前方交会的原理与意义；2、 在存在多余观测值时，利用最小二乘平差方法，经过迭代，求的外方位元素的值；3、 熟识VC编程方法，利用编程实现计算。 三、 实习原理 以单幅影像为根底，从该影像所掩盖地面范围内若干掌握点的已知地面坐标和相应点的像坐标量测值动身，依据共线条件方程，求解该影象在航空摄影时刻的像片外方位元素 ， ， ， ，共线条件方程如下：x-x0=-f*a1(X-Xs)+b1(Y-Ys)+c1(Z-Zs)/a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)y-y0=-f*a2(X-Xs)+b2(Y-Ys)+c2(Z-Zs)/a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)其中：x,y为像点的像平面坐标； x0,y0,f为影像的外方位元素；， ，为摄站点的物方空间坐标；X,Y,Z为物方点的物方空间坐标；旋转矩阵R为 ;由于此共线条件方程是非线性方程，先对其进展线性化，利用泰勒绽开得：=(x)-x+=(y)-y+像点观测值一般视为等权，即p=I；矩阵形式：V=AX-L，p=I；通过间接平差，为提高精度，增加多余观测方程，依据最小二乘平差原理，可计算出外方位元素的改正数。经过迭代计算，每次迭代用未知数的近似值与上次迭代计算的改正数之和作为新的近似值，重复计算，求出新的改正数，这样反复趋近，直到改正数小于某个限值为止。 四、 程序框图 输入原始数据 归算像点坐标x，y 计算并确定初值 ， ， ， ， 组成旋转矩阵R 计算（x）（y）和逐点组成误差方程式并法化全部点完否？解法方程，求未知数改正数计算改正后的外方位元素未知数改正数Create pass point，如图VirtuoZoNT 3.5.0软件试验步骤：(一)数据预备：1.启动 软件2.翻开测区3.翻开模型4.设置模型参数：(二)定向操作：1.内定向：2.自动相对定向：3.一般方式的肯定定向：(1)半自动量测：依次量测3个点，然后点击“猜测掌握点”。(2)肯定定向计算添加各掌握点，并调准各掌握点，使其误差小于0.03。4.定义作业区此处定义的作业区应大于自动定义的作业区5.自动生成核线影像：自动生成核线影像，单击鼠标右键弹出菜单，选择“生成核线影像”“非水平核线”，程序依次对左、右影像进展核线重采样，生成模型的核线影像。单击鼠标右键弹出菜单，选择“保存”，然后再弹出菜单，选择“退出”，然后答复界面上的提示，程序退出相对定向的界面，回到系统主界面。(三)、同名核线影像的采集与匹配1.影像匹配在VirtuoZo NT主菜单中，选择菜单“处理”“影像处理”，消失影像匹配计算的进程显示窗口，自动进展影像匹配。2.匹配结果编辑对选中区域编辑运算：(1)平滑算法：选择编辑区域后，选择平滑档次(轻、中、重);再单击“平滑算法”按钮，即对当前编辑区域进展平滑运算。(2)拟合算法：选择编辑区域后，选择外表类型(曲面、平面);再单击“拟合算法”按钮，即对当前编辑区域进展拟合运算。(四)生成DEM、等高线、正射影像及等高线叠合正射影像的操作：1.生成数字高程模型DEM在系统主菜单中，选择“产品”“生成DEM”“生成DEM(M)”项，屏幕显示计算提示界面，计算完毕后，即建立了当前模型的DEM。2.显示DEM，观看DEM是否与实际地形相符在系统主菜单中，选择“显示”“立体显示” “透视显示”项，进入显示界面，屏幕显示当前模型的数字地面模型。3.生成数字正射影像在系统主菜单中，选择“产品”“生成正射影像”项，自动制作当前模型的正射影像，屏幕显示计算提示界面，计算完毕后，自动生成当前模型的正射影像。4. 显示正射影像，观看正射影像是否有变形正射影像生成后，在系统主菜单中，选择“显示”“正射影像”项，屏幕显示当前模型的正射影像。将光标移至影像中，按鼠标右键弹出菜单，供选择不同的比例，可对影像进展缩放。5. 质量报告内定向信息：(H:GIS06hammer-165_50mic.vz_02-166_50mic.vz)-左原始影像 ( H:GIS06hammerimages-165_50mic.vz )：RMS:Mx = 0.006 My = 0.555残差: 点号 dx dy1 -0.016 0.0042 0.002 0.0043 -0.007 0.0074 0.012 0.0065 0.003 -0.0096 0.009 0.0017 0.005 -0.0108 -0.008 -0.004残差: 点号 dx dy1 -0.016 0.0042 0.002 0.0043 -0.007 0.0074 0.012 0.0065 0.003 -0.0096 0.009 0.0017 0.005 -0.0108 -0.008 -0.004右原始影像 ( H:GIS06hammerimages-166_50mic.vz )：RMS: Mx = 0.005 My = 0.555残差: 点号 dx dy1 -0.001 -0.0012 -0.001 -0.0013 0.002 -0.0074 0.004 0.0075 0.001 -0.0046 0.000 -0.0017 0.002 0.0108 -0.008 -0.002-相对定向信息：(H:GIS06hammer-165_50mic.vz_02-166_50mic.vz)-相对定向信息：左旋转矩阵：0.99995100 0.00873200 0.00467500-0.00873200 0.99996197 0.00000000-0.00467500 -0.00004100 0.99998897 右旋转矩阵：0.99996698 -0.00751000 0.003105000.00743000 0.99965900 0.02504200-0.00329200 -0.02501900 0.99968201右片旋转角 (rad)：phi = -0.00467500Omiga = 0.00000000Kappa = -0.00873200左片旋转角 (rad)：phi = -0.00310600Omiga = -0.02504500Kappa = 0.00743200残差： 点号 dq0002 0.011000166011 -0.004000166042 -0.005000166044 -0.003000166035 0.001000166037 0.0000001660310 0.004000166022 0.016000166023 -0.007000166025 0.007000165131 -0.003000165035 -0.003000165038 -0.006000165048 -0.003000165021 0.004000165049 -0.0020231650410 -0.004000165022 -0.0060006156 0.0020236155 -0.0010002265 0.0010002266 -0.0010006265 0.0000006266 0.0000001 -0.0020232 0.0050003 0.0140004 0.0020235 -0.0030006 0.0040007 0.0070008 -0.0040009 0.00300010 -0.01202311 -0.00202312 0.00100013 0.00000014 -0.00202315 -0.00800016 0.00400017 0.00100018 0.00600019 0.00202320 0.01000021 0.00202322 0.00000023 -0.00500024 0.00000025 -0.00100026 0.00800027 -0.00100028 -0.01400029 0.00000030 -0.00202331 -0.00900032 0.00400033 0.00800034 -0.00100035 -0.00100036 -0.00400037 -0.00300038 0.00500039 0.00600040 0.00000041 0.00900042 -0.00100043 0.00202344 -0.00100045 -0.00400046 -0.00500047 -0.00202348 0.00000049 0.01600050 0.00400051 -0.00100052 0.00500053 0.00202354 0.00202355 0.00500056 -0.00500057 -0.00600058 0.00600059 0.00500060 0.00000061 0.00000062 -0.00900063 0.00900064 0.00500065 0.00202366 0.00202367 -0.00600068 -0.00202369 -0.00300070 -0.00300071 -0.00202372 -0.00400073 -0.00700074 0.00000075 0.00500076 -0.00800077 0.00700078 -0.00202379 -0.00900080 0.00100081 0.00400082 -0.00300083 0.00100084 -0.00600085 -0.00500086 0.01300087 -0.00600088 0.00202389 0.00300090 0.00202391 -0.01100092 -0.00202393 0.00202394 -0.00100095 -0.01300096 -0.00202397 -0.00100098 -0.00500099 -0.003000100 0.000000101 -0.008000102 0.001000103 -0.001000104 0.004000105 0.002023106 0.000000107 -0.005000108 0.005000109 -0.009000110 -0.006000111 -0.004000112 -0.006000113 0.002023114 -0.012023115 0.014000116 -0.004000117 0.005000118 -0.001000119 -0.001000120 0.000000121 -0.001000122 0.004000123 -0.006000124 0.008000125 -0.004000126 0.008000127 -0.002023128 -0.005000129 0.006000130 0.005000131 -0.007000132 0.000000133 0.003000134 -0.002023135 -0.005000136 -0.004000137 -0.001000138 0.001000139 0.000000140 -0.005000141 0.003000142 -0.002023143 0.010000144 0.011000145 -0.001000146 -0.008000147 -0.004000148 -0.004000149 -0.009000150 -0.001000151 -0.009000152 0.002023153 0.007000154 -0.008000155 0.004000156 0.009000157 0.003000158 0.000000159 0.007000160 -0.006000161 0.001000162 0.011000163 0.003000164 -0.001000165 0.008000166 0.002023167 -0.005000168 0.013000169 -0.004000170 0.000000171 0.011000172 0.011000173 0.005000174 0.005000175 0.001000176 0.001000177 0.010000178 -0.003000179 -0.001000180 0.000000181 -0.001000182 0.000000183 0.002023184 0.013000185 0.001000186 0.009000187 0.003000188 -0.008000189 0.003000190 0.001000191 -0.006000192 0.012023193 -0.005000194 -0.009000195 0.000000196 -0.003000197 -0.003000198 -0.006000199 0.008000200 -0.015000201 -0.001000202 0.001000203 -0.001000204 0.012023205 0.009000206 0.008000207 -0.004000208 0.001000209 0.001000210 -0.011000211 -0.016000212 -0.005000213 0.003000214 -0.006000215 0.012023216 0.000000217 0.007000218 -0.007000219 -0.005000220 -0.009000221 -0.001000222 -0.008000223 -0.009000224 0.015000225 -0.004000226 -0.008000227 -0.009000228 0.016000229 -0.001000230 0.007000231 -0.008000232 -0.005000233 -0.016000234 0.001000235 0.015000236 0.006000237 -0.008000RMS: Mq = 0.006000-肯定定向信息： (H:GIS06hammer-165_50mic.vz_02-166_50mic.vz)-肯定定向信息：左旋转矩阵：0.99972469 0.01942904 0.01315499-0.01934092 0.99978989 -0.00679305-0.01328421 0.00653675 0.99989039 右旋转矩阵：0.99992472 0.00299169 0.01189982-0.00320877 0.99982804 0.01826570-0.01184313 -0.01830251 0.99976236左片摄站坐标：Xs=14867.338 ,Ys=9093.649, Zs=3275.462右片摄站坐标：Xs=16247.774 ,Ys=9078.950 , Zs=3263.634残差：6156 0.001087 -0.003654 0.0039126155 0.004336 0.014669 -0.005152265 -0.004715 -0.012243 -0.0078972266 0.000646 -0.010288 0.0097966265 0.003212 0.010743 0.0050896266 -0.004567 0.000773 -0.005747RMS: mx=0.00350 my= 0.010002mxy=0.010599 mz= 0.006571-影像匹配信息： (H:GIS06hammer-165_50mic.vz_02-166_50mic.vz)_Initia parameters_left image: rows =4320 columns =2580right image: rows =4320 columns =2580Match window width = 21Match window length = 21Searching range = 5Match grid X_interval = 21Match grid Y_interval = 21_MATCH_BLOCK = 1MATCH_LEVEL = 3MATCH_AREAS = 31 X 19589 : 0 0.0 %589 : 453 76.9 %589 : 520 88.3 %589 : 544 92.4 %589 : 547 92.9 %589 : 548 93.0 %589 : 548 93.0 %_MATCH_BLOCK = 1MATCH_LEVEL = 2MATCH_AREAS = 95 X 57_5415 : 0 0.0 %5415 : 4543 83.9 %5415 : 4980 92.0 %5415 : 5090 94.0 %5415 : 5115 94.5 %5415 : 5130 94.7 %5415 : 5143 95.0 %5415 : 5152 95.1 %5415 : 5159 95.3 %5415 : 5163 95.3 %_MATCH_BLOCK = 1MATCH_LEVEL = 1MATCH_AREAS = 205 X 122_25010 : 0 0.0 %25010 : 21733 86.9 %25010 : 23037 92.1 %25010 : 23293 93.1 %25010 : 23442 93.7 %25010 : 23537 94.1 %25010 : 23595 94.3 %25010 : 23652 94.6 %25010 : 23739 94.9 %25010 : 23787 95.1 %25010 : 23812 95.2 %25010 : 23832 95.3 %25010 : 23855 95.4 %25010 : 23873 95.5 %_MATCH_BLOCK = 1MATCH_LEVEL = 0MATCH_AREAS = 205 X 122_25010 : 0 0.0 %25010 : 20350 81.4 %25010 : 21419 85.6 %25010 : 21842 87.3 %25010 : 22083 88.3 %25010 : 22244 88.9 %25010 : 22372 89.5 %25010 : 22498 90.0 %25010 : 22595 90.3 %25010 : 22663 90.6 %25010 : 22722 90.9 %25010 : 22757 91.0 %25010 : 22787 91.1 %25010 : 22819 91.2 %25010 : 22853 91.4 %25010 : 22867 91.4 %25010 : 22888 91.5 %25010 : 22900 91.6 %-DEM 检查点中误差:-DEM 文件： H:GIS06hammer-165_50mic.vz_02-166_50mic.vzproduct-165_50mic.vz_02-166_50mic.vz.dem检查点(掌握点)文件：H:GIS06hammerhammer.pas-点号 X Y Z dZ2265 14787.371 9101.982 786.751 -3.4382266 16327.646 9002.483 748.470 -0.4116155 16340.235 10314.228 751.178 -7.0376156 14947.986 10435.860 765.182 0.0026265 14888.312 7769.835 707.615 -4.7496266 16232.309 7741.696 703.121 -0.614点数 = 6均值 = -2.7肯定均值 = 2.7均方根 = 3.8点号 百分比dZ <= 1.0 : 3 50.01.0<= 2.0 : 0 0.02.0<= 3.0 : 0 0.03.0 < dZ < 4.0 : 1 16.74.05.06.010.020.0-VirtuoZo 影像文件信息 (正射影像):-影像文件名： H:GIS06hammer-165_50mic.vz_02-166_50mic.vzproduct-165_50mic.vz_02-166_50mic.vz.orl行列数 行数X列数： 3191 X 1791模型颜色： 24位彩色影像X-方向像素大小： 0.100000 毫米地理信息：-影像比例尺： 1 : 10000旋转角： 0.00000 度X-方向地面分解率： 1.000000Y-方向地面分解率： 1.000000左下角坐标 x,y：14630.000 7450.000右下角坐标 x,y：16420.000 7450.000左上角坐标 x,y：14630.000 10640.000右上角坐标 x,y：16420.000 10640.000-VirtuoZo 影像文件信息 (等高线影像):-影像文件名： H:GIS06hammer-165_50mic.vz_02-166_50mic.vzproduct-165_50mic.vz_02-166_t行列数 行数X列数： 3191 X 1791模型颜色： 8位单色影像X-方向像素大小： 0.100000 毫米地理信息：-影像比例尺： 1 : 10000旋转角： 0.00000 度X-方向地面分解率： 1.000000Y-方向地面分解率： 1.000000左下角坐标 x,y：14630.000 7450.000右下角坐标 x,y：16420.000 7450.000左上角坐标 x,y：14630.000 10640.000右上角坐标 x,y：16420.000 10640.000-VirtuoZo影像文件信息(等高线叠加正射影像):-影像文件名： H:GIS06hammer-165_50mic.vz_02-166_50mic.vzproduct-165_50mic.vz_02-166_50mic.vz.orc行列数 行数X列数： 3191 X 1791模型颜色： 24位彩色影像X-方向像素大小： 0.100000 毫米地理信息：-影像比例尺： 1 : 10000旋转角： 0.00000 度X-方向地面分解率： 1.000000Y-方向地面分解率： 1.000000左下角坐标 x,y：14630.000 7450.000右下角坐标 x,y：16420.000 7450.000左上角坐标 x,y：14630.000 10640.000右上角坐标 x,y：16420.000 10640.000(五)生成数字影像图1.进入测图界面在VirtuoZo NT系统主菜单中，选择“测图” “IGS数字测图”项，调用测图模块，屏幕弹出测图界面。2.新建或翻开测图文件新建一个测图文件：选择FileNew Xyz File项，屏幕弹出文件查找对话框，输入一个新的xyz文件名，弹出测图参数对话框。在对话框中输入各项测图参数：成图比例尺(分母);高程注记的小数位数;流数据压缩容限(单位：毫米);图廓坐标：Xtl、Ytl(左上角)、Xtr、Ytr(右上角)、Xbl、Ybl(左下角)、Xbr、Ybr(右下角)。选择“保存”按钮后，将创立一个新的测图文件。此时屏幕弹出矢量图形窗并显示其测图的图廓范围。翻开一个测图文件：选择“文件”“翻开”项，此时弹出文件查找对话框，选择一个已有“*.xyz文件”，翻开后，屏幕显示当前的矢量图形文件。3.装入立体模型当翻开测图文件后，方可翻开立体模型。在菜单栏中选择“文件”“翻开”项，在文件查找对话框中，选择一个模型*.mod(或*.set)文件，翻开后，屏幕弹出影像窗显示立体影像。4.影像贴图与矢量图形的层掌握(1)矢量贴图：按下 图标，可将测量的结果(矢量图形)显示在立体影像上，便于检查遗漏和所测地物的精度。(2)层掌握：在数字化测图中，同一种地物为一层，每一层都有一个属性码(或层号)。所测的地物都被分层治理，层掌握就是对地物分层治理的工具。5.影像显示(1)左右影像分屏显示，由立体反光镜观测立体;(2)立体显示双影像：通过硬件的支持，左右影像交替显示，戴上相应的立体眼镜，可以进展立体观测。(六)实习中遇到的问题1、新建数据文件夹时应留意哪些问题?答：给文件夹命名时应当留意不要含有中文和空格，用字母为宜，以免程序出错。2、使用测区参数界面下的重置模型参数功能应当留意哪些问题?答：首先测区参数界面的参数项中不能有空白项;在填入参数时掌握点文件、加密点文件、相机参数文件的名字可任意命名，但是要切记不能使得这三者的名字重名，否则，可能会导致文件的冲突，影响到内定向、肯定定向的成果，甚至无法正常的采集核线影像。3、有多个相机的测区如何处理?答：分别在测区名目下建立多个相机参数文件(留意影像参数要与之对应)，分别作内定向即可。4、定义核线影像范围应留意那些事项?答：首先：定义核线范围以将掌握点划在范围内为宜，但不能超控过多;其次：应结合实际地形状况，如高山地或大比例尺城区，由于左右片视差较大，就应适当将核线范围划大些。(七)实习心得通过这次试验，我系统学习了VirtuoZo的使用操作流程，了解了VirtuoZo的根本功能，一般作业过程及主要产品的制作过程。其中主要有把握创立/翻开测区及测区参数文件的设置。把握参数文件的数据录入。通过对模型定向的作业，了解数字影像立体模型的建立方法及全过程，并能较娴熟地应用定向模块进展作业，满意定向的根本精度要求。把握核线影像重采样，生成核线影像对。把握匹配窗口及间隔的设置，运用匹配模块，完成影像匹配。把握匹配后的根本编辑，能依据等视差曲线(立体观看)发觉粗差，并对不行靠区域进展编辑，到达最根本的精度要求。把握DEM格网间隔的正确设置，生成单模型的DEM。把握正射影像辨别率的正确设置，制作单模型的数字正射影像。通过DEM及正射影像的显示，检查是否有粗差。把握拼接区域的选定及确定拼接产品的路径。把握DEM拼接及自动正射影像镶嵌。分析拼接精度。理解数据格式输出的意义。了解VirtuoZo NT系统的数据格式输出的详细操作。虽然这次做的成果并不是很完善，与试验指导中的还存在肯定的差距。但通过对实习成果的分析，了解数字产品的根本质量要求，对其进展肯定的改善，使其精度有肯定的提高。而且通过这次实习，我对数字摄影测量数据猎取有了更深刻的了解，同时对摄影测量课程有了更深更详细的体会。1.正射影像2.等高线影像3.等高线叠合正射影像4.镶嵌 5.DEM拼接数字摄影测量实习报告模板 小编为大家整理的数字摄影测量实习报告模板，供大家参考。更多阅读请查看本站实习报告网频道。 实 习 报 告 实习名称： 班 姓 级: 名: 数字摄影测量实习 09 测绘一班 刘胜 试验室 X5504 实习地点： 实习指导教师: 龚涛 实习时间： 20xx.920xx.10 西南交大地学学院 数字摄影测量实习报告： 一 LpS 简述 LpS 工程治理器是一个基于 WINDOWS 的综合数字摄影测量软件包，可 以对来 自不同类型的航空遥感相机及卫星传感器的图像进展快速和准确地 三角测量和 正射校正，与传统的三角测量和正射校正相比，可以极大的削减费用 和时间可 以处理各种各样的图像数据，诸如来自不同的摄影相机、不同的卫星 传感器、 不同的航空 GpS 数据等，处理过程涉及许多不同类型的几何模型。 二、数字摄影测量处理过程 1 创立 LpS 工程文件 2 向 LpS 工程加载图像 3 刺点 4 自动量测图像同名点 5 执行航空三角测量 DTM+等高线 DTM+等高线 6 图像正摄校正处理 处理前 处理后 掌握点坐标 三、LpS 数字摄影测量系统的应用 Leica photogrammetry SuiteLpS 是徕卡公司最新推出的数字摄影测量及遥 感处理软件系列。LpS 为影像处理及摄影测量供应了高精度及高效能的生产工 具、它可以处理各种航天(最常用的包括卫星影像 QuickBird、IKONOS、SpOT5 及 LANDSAT 等等)及航空(扫描航片、ADS40 数字影像)的各类传感器影像定向 及空三加密，处理各种数字影像格式，黑/白、彩色、多光谱及高光谱等各类数 字影像。LpS 的应用还包括矢量数据采集、数字地模生成、正射影像镶嵌及遥感 处理，它是第一套集遥感与摄影测量在单一工作平台的软件系列。 LpS 制作 DOM 的全过程如下： LpS 数字摄影测量系统制作 DOM 详细制作过程如下： 首先创立工程文件，选择相机类型，设置投影参数，输入相片参数，创立相 机参数，导入外方为元素;其次数据处理，内定向，人工选择一个点后，自动完 成内定向。建立金字塔影像，加载掌握点文件，并在图上刺出相应的点!一般说 来，选择 6 个匀称分布的点作为掌握点，其他的设为检查点。同名点自动匹配， 三角测量，直接进展空三解算，再接着生成 T