1、摄影测量及遥感实习报告HEBEI UNITED UNIVERSITY摄影测量与遥感实习报告学生：学 号：专业班级：学 院：轻工学院指导教师：谷海红2011年12月30日1实习目的摄影测量与遥感实习是摄影测量学和遥感技术及应用两门课程的综合实习课。本课程的任务是使学生通过编写程序实现单影像空间前方交会计算，掌握非线性方程线性化的过程、相应数据读入与存储的方法以及迭代计算的特点。掌握遥感影像预处理、信息获取、实地调绘及制图的方法等作业程序。使学生进一步稳固遥感上机课的知识。从而更系统地掌握摄影测量与遥感技术。通过实习使学生熟练地掌握摄影测量及遥感的原理，信息获取的途径，数字处理系统和应用处理方法。
2、进一步稳固和深化理论知识，理论与实践相结合。2实习容1)控制点的采集2)遥感影像的几何校正。3)遥感影像的野外调绘、属性调查。4)遥感影像的室清绘；5)遥感影像图像矢量化；6)属性信息的录入。7)摄影测量单像空间前方交会程序窗体设计8)编写代码9)运行程序，得出结果3实习设备与资料硬件设备：计算机的根本配置CPU型号、存大小，硬盘空间软件：ENVI或ArcGIS的软件名称及版本空间前方交会原理：根据影像覆盖围一定数量的分布合理的地面控制点其像点和地面点的坐标，利用共线条件方程求解像片外方位元素。计算过程：获取数据m, x0 , y0 , f , Xtp, Ytp, Ztp量测控制点像点坐标x,
3、y确定未知数初值Xs0, Ys0, Zs0, 0, 0, 0组成误差方程式并法化解求外方位元素改正数检查迭代是否收敛实习数据：遥感实习数据：市高分辨率遥感影像摄影测量实习数据：航摄仪的方位元素：f=153.24mm，x0=y0=0.0mm，摄影比例尺为1:50000； 4个地面控制点的地面坐标及其对应像点的像片坐标如下表1 左像片点号像片坐标(mm) 地面点坐标(m) x y X Y Z 1 -86.15 -68.99 36589.41 25273.32 2195.17 2 -53.40 82.21 37631.08 31324.51 728.69 3 -14.78 -76.63 39100.
4、97 24934.98 2386.50 4 10.46 64.43 40426.54 30319.81 757.31 表2 右像片点号像片坐标(mm) 地面点坐标(m) x y X Y Z 1 57.67 35.54 36589.41 25273.32 2195.17 2 67.33 -23.12 37631.08 31324.51 728.69 3 32.21 72.21 39100.97 24934.98 2386.50 4 -44.12-66.21 40426.54 30319.81 757.31 且dX dY dZ 0.005限差d d dOpen Image File，如图1所示，弹
5、出翻开图像文件窗口。图1 翻开图像文件建立两个显示窗口，分别翻开拼接好的市地形图与遥感影像图。选择校正控制点。按如下步骤操作：在主菜单下，Map Registration Select GCPs:Image to Image Regi 在弹出对话框中选择拼接后的市地形图为Base Image, 遥感影像图为Warp Image OK 弹出对话框。然后在两幅图上找到一样的控制点，点击Add Piont添加更多的控制点至少3个。将控制点数据保存。同样的步骤，继续校正至少十个控制点数据。控制点找完保存后进展校正。在Ground Control Points Selecti.对话框中，Options
6、Warp File 在“Input Warp Image对话框中选择遥感影像图 OK Choose，选择保存路径，然后保存。建立连接查看校正后的影像图和市图间的连接，查看校正情况。显示窗口下，Tools Link Link Displays OK。5.3 遥感图像的野外调绘、属性调查在进展野外调绘之前，将调绘图平放在大图板上，然后再将比调绘图稍大一些的透明纸盖于调绘图上，用胶带粘好，连同调绘图用夹子固定于大图板。图像固定好后，首先在室直接判读地物、地貌，对于清晰的可以先用铅笔与尺子按照国家标准地形图图式将其边界勾勒在透明纸上。对于在室无法直接判读的，均应带上遥感调绘图到实地对各种地物进展比照。
7、完成准备工作后，那么开场进展野外调绘。首先，选定调绘路线。一般在实际调绘工作中，均应按行调绘或者按列调绘，即地物地貌的调绘要连续进展，防止调绘不连贯。然后，依照选好的路线，按行或按列进展调绘。调绘时，调绘点选择在地势较高的地方，主要调查属性，比方建筑物的建筑材料、层数、名称，道路的名称，各种绿地的分类等。对于已经在室绘好的地物地貌，那么在实地检查是否有改变，有改变的按照地形图图式将它们修改绘制在对应的透明纸上。对于原先不确定的地物地貌，按照实地情况在透明纸上按图式勾勒出它们的边界并注明属性。直至将整图所包括的容全部调绘完毕。调绘的同时，要将调绘的容特别是与实地有所出入的地物地貌及其属性填写在调
8、绘手簿中。5.4遥感影像的室清绘像片调绘可采用先室判读清绘，后野外检核和调查，再室修改和补充清绘的方法。野外调绘完成后，当天进展分色清绘。根据要求，用红、蓝、绿、黑四种颜色进展清绘。其中应注意：当地理名称注记过密时，可适当取舍。调绘工作应按照国家标准的地形图图式进展，说明性质的注记应采用“简注表，不得任意命名。实际调绘工作中，均应按行调绘或者按列调绘，不得随意调绘。调绘要按照实地情况严格进展，不得伪造，篡改。在调绘好的透明纸上，图名注于调绘片正上方，像片号注于调绘片右上角，调绘者及调绘日期等信息在调绘航片的右下角。在清绘时，要遵守“片片清、天天清、幅幅清。5.5 遥感影像图像矢量化5.5.1
9、遥感影像翻开翻开ENVI4.6，点击File-Open Image File，选择影像如图2，图2 选择影像在Available Bands List窗口中点击Load RGB按钮，在#1窗口中加载影像，翻开后，显示遥感影像如图3，图3 显示遥感影像5.5.2 图层的建立在图3显示的窗口上，选择Overlay vectors，弹出对话框#1 Vector Parameters：，如图4图4 矢量参数设置在此窗体上选择File Create New Layer.在弹出的对话框中填写新建图层的名字与保存位置，点击确定。再次弹出如图4窗口，且新建的图层名字以显示在文本框中。按照这种方法我新建了原有建
10、筑、新建地物、公路、居民区及办公区道路、草坪及绿化带、变更建筑、辅助7个图层。5.2.3 矢量化按照清绘完成的纸图，选择相应图层的颜色与线条等。在#1 Vector Parameters：中，选中所要设置的图层并右击，在弹出的菜单中选择Edit Layer Properties,弹出对话框如图5。图5 设置图层样式在我组图层设置时选择原有建筑：黑色实线；新建地物：红色点线，3Thick；公路：黑色实线，2Thick；居民区及办公区道路：黑色实线，1Thick；草坪及绿化带：边框为绿色实线，用绿色点线填充；变更建筑：红色实线；辅助：蓝色实线。设置完成后，按照着色后的纸图在如图4显示的窗口中，选择
11、Mode菜单下的Add New Vectors ,并选择相应的点、线、多边形、圆、方，在Zoom窗口中以遥感影像图为底图添加矢量元素。5.6 属性信息的录入当完成一个矢量元素后，在#1 Vector Parameters：窗体中选择Mode菜单下的Edit Existing Vectors,然后选择Edit菜单下的Add New Attributes,显示如图6图6添加属性信息选项按照图6显示选择，设置完成后，点击确定。5.6.1 属性字段的添加完成上步后，弹出窗体 如图7图7 属性字段设置设置完成后，点击确定，弹出窗口，如图8图8图层属性信息在图7显示的窗口中选择Options菜单下的Add
12、 Record Column，即可再次显示图7窗口，然后就可以按照需求添加相应的字段了。5.6.2 属性信息的录入在图8显示的窗口中双击想要添加的字段，添加相应的属性信息，如图9图9 属性信息的录入5.6.3 属性信息的修改、删除想要修改属性信息时，可以在#1 Vector Parameters：中选中要修改的属性信息所在的图层并选择Mode菜单下的Edit Existing Vectors,然后选择Edit菜单下的View/Edit/Query Attrbutes,在Zoom窗口单击选中要修改的矢量元素，此时在属性窗口中该元素的属性行也被高亮显示出来，然后就可以对选中的元素进展节点的移动、添
13、加、删除，属性信息的修改，同时也可以再选中元素后在右击显示的菜单中选择删除此元素，在相应的属性信息列表中的属性信息也会被删除。要删除某一图层的属性字段时，可在图7显示的窗口中选中该列，然后选择Options菜单下的Delete Record Column,那么该列属性字段就会被删除了。6 摄影测量实习过程6.1 窗体设计6.1.1计算过程分析根据实习要求及试验数据实现左右两影像空间前方交会的计算，并将最后计算得到的旋转矩阵及外方位元素输出，得到如图10的程序流程图。图10 空间前方交会程序流程图6.1.2设计窗体根据设计要求将已有的控制点像点坐标及物方坐标可通过输入界面进展输入，已可定义成文本
14、文件由程序读取文本文件；定义变量和数组来存储数据；最终计算结果的输出可输出为文本文件。设计出主要计算窗体如下：图11 数据输入界面图12 结果显示界面6.2 编写代码private void button1_Click(object sender, EventArgs e) try for (int j = 0; j 5; j+) for (int i = 0; i 3; i+) Xi = Convert.ToDouble(dataGridView10, i.Value); Yi = Convert.ToDouble(dataGridView11, i.Value); Zi = Convert
15、.ToDouble(dataGridView12, i.Value); xi = (Convert.ToDouble(dataGridView13, i.Value); yi = (Convert.ToDouble(dataGridView14, i.Value); #region/ 初始化 for (i = 0; i 3; i+) S1 += Xi; S2 += Yi; S3 += Zi; /确定外方位元素的初始值 = 0; = 0; = 0; if (textBox1 .Text = | textBox2 .Text = ) MessageBox.Show(请输入比例尺和焦距数据); re
16、turn; m = Convert.ToDouble(textBox1 .Text ); / m = 50000; f = Convert.ToDouble(textBox2 .Text ); / f = 0.15324; Xs0 = S1 / 3.0; Ys0 = S2 / 3.0; Zs0 = m * f + S3 / 3.0; /像片坐标的单位换算 for (i = 0; i 3; i+) xi /= 1000.0; yi /= 1000.0; catch (Exception ex) MessageBox.Show(this, Error01 + ex.Message, 初始化错误);
17、 #endregion / try do K+; cal_R(); cal_xy(); cal_A(); WindowsApplication1.Matrix.Transpose(A, AT); WindowsApplication1.Matrix.Matrix6_6(AT, A, AT_A); E = WindowsApplication1.Matrix.ReverseMatrix(AT_A, 6);/将A的转置与A的乘积存放在临时矩阵E中 WindowsApplication1.Matrix.Matrix6_6(E, AT, ATA_ni_AT); WindowsApplication1.
18、Matrix.Matrix6_1(ATA_ni_AT, L, det); /计算外方位元素 Xs0 += det0, 0; Ys0 += det1, 0; Zs0 += det2, 0; += det3, 0; += det4, 0; += det5, 0; listBox1.Items.Add(String.Format(第0次迭代, K); listBox1.Items.Add(.); listBox1.Items.Add(String.Format( Xs0 = 0, Ys0=1, Zs0=2, Xs0, Ys0, Zs0); listBox1.Items.Add(String.Form
19、at( = 0, =1, =2, , , ); listBox1.Items.Add(.); if (Math.Abs(det0, 0) = 0.0005 & Math.Abs(det1, 0) = 0.0005 & Math.Abs(det2, 0) = 0.0005 & Math.Abs(det3, 0) = 0.00005 & Math.Abs(det4, 0) = 0.00005 & Math.Abs(det5, 0) 500 & Iscontiue = true) if (MessageBox.Show(this, 已经超过500次n是否继续, 提示, MessageBoxButto
20、ns.YesNo, MessageBoxIcon.Question, MessageBoxDefaultButton.Button2) = DialogResult.No) Result1(); break; else if (K 1000) MessageBox.Show(迭代超过1000次); Result1(); break; else Iscontiue = false; while (true); catch (Exception ex) MessageBox.Show(this, Error02 + ex.Message, 迭代错误); /左片计算6.3 运行程序并得出结果输入数据
21、，点击计算，此时可查看迭代过程与计算结果，如图13和图14。图13迭代过程图14 计算结果7实习心得通过本次实习，我不仅加深了对课本上抽象知识的理解，通过实际操作明白了摄影测量和遥感应用的原理、方法、原那么要求与实际操作等，而且更深刻地明白了摄影测量和遥感对于现代土地资源调查和环境规划有着重要的意义。 同时在这次实习中我获得了很多启示，并且我相信这次实习将会成为我记忆深处的一段美好回忆。首先，小组成员的分工合作，默契配合是高效完成任务的重要途径。在前几次实习的过程中，我们从陌生走向熟悉，从彼此的拘谨走向大家相互调侃的融洽关系，这为我们高效完成这项任务提供了很好的条件。在实习的前期，我们认真地分
22、配各自任务，规划好每天的工作安排，保证每个人都能够充分地利用时间完成任务。其次，在遥感实习过程中，不可防止的遇到了一些困难，通过教师的提点、理论知识与课上讲解的结合分析以及我们的认真地对待最终还是很好的克制了，使我深刻地体会了将学以致用的重要性。同时充分开动脑筋，利用所学知识解决外业工作中的问题，对于有效提高数据的准确性起了很大的作用。在摄影测量实习中要根据选择的编程语言C#结合计算公式与其他实例，设计出符合要求的程序，让我们对课本上的繁琐的计算过程有了更深的理解，并加深了印象。最后，端正的态度对于实习结果是十分重要的。我们组员的相互配合，以一个严谨的态度完成每一项工作，这使得我们得到一个客观的有经济价值的结果起到了很大的作用。时光如白驹过隙，我会十分珍惜我们每一次实习的时机，我十分感教师的殷勤教导。
