数字产品生产实训报告.doc

数字产品生产实训报告

数字产品生产实训

实 习 报 告

学 院 测绘科学与技术学院
专业班级 测绘工程1203
姓 名 长孙建坤
学 号
指导教师 邱春霞 张静 娄宁

日 期 01月15日

目录
1. 实习目的 1
2. 实习概况 1
3. 1:10000 DRG制作 2
3.1. 基础数据 2
3.2. 应用软件 2
3.3. 制作步骤 2
3.3.1. 图像预处理 2
3.3.2. 数据格式转换 3
3.3.3. 图像校正 4
3.3.4. 图像裁剪 8
3.3.5. 图像拼接 8
3.3.6. 质量控制 9
4. 1:10000 DLG制作 10
4.1. 基础数据 10
4.2. 应用软件 10
4.2.1. MAPGIS6.7输入编辑模块 10
4.2.2. MAPGIS6.7中的快捷键 10
4.2.3. MAPGIS6.7中几个重要概念 11
4.3. 制作步骤 11
4.3.1. 读图分析 11
4.3.2. 新建工程 11
4.3.3. 设置工作目录 12
4.3.4. 创立工程文件 12
4.3.5. 创立图例板 13
4.3.6. 地物信息采集 14
4.4. 质量控制 19
5. 1:10000 DEM制作 20
5.1. 基础数据 20
5.2. 应用软件 20
5.3. 制作步骤 20
6. 1:10000 DOM制作 22
6.1. 基础数据 22
6.2. 应用软件 22
6.3. 制作步骤 22
7. 基于VirtuoZo的4D产品生产 26
7.1. 操作步骤 26
7.1.1. 新建测区设置测区参数文件 26
7.1.2. 模型生成 28
7.1.3. 模型定向与生成核线影像 29
7.2. 单模型DEM的生成 33
7.3. DOM数据制作 34
7.4. DLG数据制作 34
8. 实习心得 36

1. 实习目的
随着科学技术的不断更新,计算机技术和3S技术也得到了飞速的发展。数字城市、数字地球等一系列空间描述的概念在不断的提出。显然,传统的测绘已经不能满足社会发展的需求,必须在其基础上进行革新。在新科学、新理念的支撑下,测绘与地理信息系统两门学科的整合将成为一种趋势,单纯的测绘人或者地信人也将不复存在。21世纪,全球地理信息产业将会得到飞速的发展,地球村的伟大目标也将会变成现实。
数字产品生产实习既是地信专业的一门专业实习,也是测绘工程专业的一次大型综合类实习,它结合了地理信息系统原理与应用、摄影测量学、数字摄影测量学,以及地图学、遥感原理等课程。
本次实习涉及的内容多、范围广、综合性强。按照生产实践过程要求,使学生在理论学习的基础上,真正掌握基于GIS软件、遥感图像处理软件、数字摄影测量系统等的1:1万4D(DRG、DLG、DEM、DOM)产品的主要制作方法和技术流程,以及4D产品之间的关系、4D产品制作的质量控制措施。
2. 实习概况
本次实习以西安市周边地区1:1万扫描地形图为数据源如下表1-1所示,以及西安市周边地区资源三号全色卫星影像,分辨率为2.1米。借助地理信息软件MapGIS/ArcGIS、遥感图像处理软件ENVI/ERDAS、数字摄影测量系统Virtuo-Zo以及Photoshop等辅助软件完成4D产品生产的方案设计、制作方法选定和质量控制等环节。

表格 2-1 西安市周边地区地形图
I49G04
陕棉十厂
I49G04
西安西站
I49G04
西安市
I49G043014
鱼化寨
I49G043015
电子科技大学
I49G043016
大雁塔
I49G044014
傅村
I49G044015
杜　城
I49G044016
曲江池

3. 1:10000 DRG制作
3.1. 基础数据
本次1:10000DRG制作的基础数据采用的是1998年出版的西安市周边的九幅1:10000扫描地形图,数据的格式为标准TIFF文件。数据采用新版地图标准分幅,扫描分辨率为300DPI。平面坐标系统采用西安80坐标系,高程系统采用1985国家高程系统。
3.2. 应用软件
当前,市场上商用的4D产品生产软件种类繁多,功能也比较丰富,但4D产品生产的流程及原理都基本相同。本环节所能涉及到的软件有MAPGIS、PHOTOSHOP、ARCGIS。
3.3. 制作步骤
3.3.1. 图像预处理
由于扫描地形图在扫描过程中存在灰度失真的情况,加之,地形图年代较久,因此扫描后的地形图色彩显示比较模糊。这在图像校正的过程中对于控制点点位中心的判断造成影响,从而增大校正误差。为了改进图像的视觉效果,提高图像的清晰程度,将图像转换为一种更适合于人机交互处理和分析的形式。
本次采用PHOTOSHOP图像处理软件对西安市周边的九幅扫描地形图进行灰度变换,以增强图像的显示效果。经过图像处理后,我们发现图像的清晰程度得到了有效的改进,如图3-1所示。
图表 3-1 图像增强对比

图像未处理前
图像处理后

3.3.2. 数据格式转换
本次DRG的制作采用国产软件MAPGIS。MAPGIS是由中地数码科技有限公司研发的具有独立自主知识产权的大型地理信息系统平台,主要功能包括数字制图、数据库管理、空间分析等。该软件在国内外取得较高的知名度。
由于本次采用的基础数据是西安市周边的九幅扫描TIFF文件,该文件格式MAPGIS不能进行直接编辑,因此需要将其转换为MAPGIS的标准格式*.MSI格式。具体操作流程如下:
运行软件狗,打开软件,设置工作环境目录,在图形处理模块中选择图像分析,点击打开该模块后显示界面如图3-1。
图 3-1

在文件菜单下点击输入数据,弹出如下图3-2所示的对话框,在数据格式中选择TIFF文件格式,然后点击添加目录,在对应的数据路径下找到目标文件,然后点击打开,对应的西安市九幅TIFF文件就被加载上来。
选中所有文件,然后点击转换按钮,所有的文件将被转换成MAPGIS的标准文件格式,同时输出到原来的文件路径下。就此完成图像格式的转换工作。
图 3-2

3.3.3. 图像校正
DRG的制作实质是将纸质地形图转换为电子的栅格地形图,可是由于纸质地形图在扫描的过程中不可避免的存在几何畸变。因此,需要利用图框点理论坐标以及其对应的图面坐标,经过内插控制点,形成密集的图面控制网,然后经过平差的方法求出扫描图像的最或然坐标。具体的操作流程如下:
打开图像处理模块,点击文件,打开影像,选择第一幅影像i49g04 打开,如图3-3所示。
图 3-3

在镶嵌融合菜单选择DRG生产,在三级菜单下选择图幅生成控制点,根据该幅影像所在标准分幅的图幅号可计算出图像四个脚点的理论坐标值,弹出窗口如图3-4所示。
图 3-4

根据扫描地形图的实际情况设置图框信息时,应将坐标系统选则为西安80坐标系,网格间距设为1公里,图框类型选则四点图框。大地坐标系一栏如果选则,则生成的坐标系统以米为单位,否则,以毫米为单位。在这里我们选则以毫米为单位。
根据图幅号生成的四个脚点坐标分别在图上选则对应的实际图面坐标位置,然后一次添加。再添加坐标时,先将光标移到图框坐标上点击鼠标左键激活该点。然后在图像上对应位置寻找内图廓点坐标。如图3-5所示。为了保证控制点具有较高的精度,以便为后续工作提供可靠的控制基础,在寻找点位过程中应当尽量将图像放大到像元。同时应根据像元的灰度变化规律,准确的判断出点位。
图 3-5

根据以上方法添加完左上角控制点后,采用同样的方法添加完剩余的三个控制点。然后点击生成GCP,软件会根据四个脚点的控制点坐标在图像上以1公里为间隔生成一系列的公里网格交点和公里网格与内图框交点。在镶嵌融合菜单下点击控制点浏览和控制点信息能够浏览控制点的图面位置和大致坐标信息。如图3-6所示。
图 3-6

由于生成的控制点的校正坐标并没有精确地与实际图面位置相对应,因此需要将内插的所有控制点逐点进行精确定位。如图3-7所示。完成所有控制点的修改后,点击保存所有控制点信息。在镶嵌融合菜单下设置校正参数。由于图像经过校正后原来的坐标位置及灰度值都将发生改变,因此需要设计合适的采样和量化方法,以期达到较高的校正精度和效果。本次实习数据量大小适中,因此采样方法选择三次多项式,量化方法采用双线性内插法。
图 3-7

在二级菜单DRG制作下点击逐网格校正,然后软件会根据生成的控制点逐网格对图像进行校正。此时会生成一个新的校正后文件Ji49g04 .msi。
完成校正后,软件会生成一个校正后的质量评估报告。报告中显示了校正后图像点位中误差、图廓边长及对角线误差信息。如图3-8所示。本次实习规定图像校正后点位中误差不得大于0.8米。横向、纵向误差不得大于2米,对角线误差不得大于3米。由报告可知本次校正后的图像点位中误差都基本达到厘米级精度。
图 3-8

图像i49g04 DRG制作至此则已经完成。此时图像已经大幅度的消除了几何畸变。采用相同的方法分别对剩余的8幅图像进行处理,便可得到相应的DRG产品。
3.3.4. 图像裁剪
图像完成单幅图像的DRG制作后并不能满足对整个西安市及其周边4D产品的制作要求,因此需要对9幅图像进行裁剪,去除图框,然后拼接在一起形成一副西安市及其周边地区完整的DRG产品。本次裁剪采用MAPGIS的标准图框裁剪功能。由于9幅扫描图像是按国家1:10000标准分幅编制的,因此利用标准图框裁剪功能对校正后的DRG进行裁剪可达到零误差的精度。对图像i49g04 的裁剪结果如图3-9所示。
图 3-9

3.3.5. 图像拼接
MAPGIS中采用影像镶嵌的功能可实现两幅图像之间的高精度拼接。可是由于盗版软件的原因,在拼接的过程中只能完成三幅影像之间的拼接,形成一条带状DRG,当两条带状DRG进行拼接时,软件则会报错。针对该问题,经过网上查询,得出结论是由于盗版软件的原因。
因此,在拼接过程中采用PHOTOSHOP进行拼接。先在MAPGIS中将三张单独的DRG拼接成一条带状DRG,最后采用PHOTOSHOP完成三条带状DRG的拼接任务,最终得到结果图如图3-10所示。
图 3-10

3.3.6. 质量控制
DRG制作的质量控制主要从三个方面进行控制:点位定位精度、采样量化方法选择、校正限差控制。
点位定位精度:控制点的定位精度直接影响到图像校正的精度,以及后续4D产品生产的精度。因此必须进行精确地确定点位位置。在点位定位时应尽量放大图像,使得图像显示到像素级别。同时应结合周围像元的灰度值进行综合判断来确定控制点的图面坐标位置。
采样量化方法选择:采样量化选择不同的方法将达到不同的校正效果和精度。可是方法的复杂程度与图像的数据量及电脑处理器的配置高低有着密切的关系,因此应该综合图像的数据量选择合适的采样和量化方法。本次实习采用的数据基本在300M-400M之间,数据量算适中,计算机处理器为酷睿i3。因此,综合考虑,采样方法选择三次多项式,量化方法采用双线性内插法。
校正限差控制:校正限差是DRG生产质量的最终保证条件。满足精度要求则为合格,否则,应进行重新制作。本次实习规定图像校正后点位中误差不得大于0.8米。横向、纵向误差不得大于2米,对角线误差不得大于3米。由报告可知本次校正后的图像点位中误差都基本达到厘米级精度。

4. 1:10000 DLG制作
4.1. 基础数据
DLG制作是利用已经完成的DRG产品进行生产的。由于DLG生产过程中是纯人工进行的,因此工作量比较大,因此我小组将测区进行分块分派给每个组员进行处理。最后组员之间进行拼接、组与组之间进行拼接。我负责的是数字栅格图像i49g043016的右上角区域。
4.2. 应用软件
4.2.1. MAPGIS6.7输入编辑模块
本次DLG制作采用MAPGIS6.7软件。在这个阶段主要利用的是MAPGIS6.7中的图形处理模块下的输入编辑模块。这部分几乎是整个矢量化的主体部分,几乎90％的工作在这个部分完成。在该模块下能够实现多种类型数据的输入、输出、编辑、统计、分析和管理等功能。同时,也能创立和编辑几何属性和拓扑关系。
4.2.2. MAPGIS6.7中的快捷键
Shift＋鼠标左键:捕捉线上点或者线的端点;
Ctrl＋鼠标右键:闭合线;
F5:图面放大;
F6:图面移动;
F7:图面缩小;
F8: 加点(如果你选择交互式矢量化,必须用这个键,不能用鼠标左键加点);
F9:撤销上一段所画的线;
F11: 画线方向反向;
F12:选择造线捕捉方式(如下图,用的比较多的是”靠近线(母线加点)”);
Ctrl:按下CTRL的同时选择图元,能够多个选择;
Ctrl＋C:拷贝;
Ctrl＋X:剪切;
Ctrl＋V:粘贴。
4.2.3. MAPGIS6.7中几个重要概念
参数:主要是指其外形,属性则是其本质。在MAPGIS6.7中,参数就是线的类型。
属性:就是所代表的地物。
地图参数:是指对整个矢量化界面的参数设置。
属性结构:是指其结构主要包括多少项的一个分类。
4.3. 制作步骤
4.3.1. 读图分析
测区位于西安市碑林区与新城区交界地带,面积约9平方公里。测区包括城区和郊区两部分,以城区为主。城区内居民地密集,地物种类繁杂。在测区的南边为郊区地带,主要分布着一些村庄和果园,地物类型较为简单。
经过大致了解,测区内地物类型能够分:测量控制点、居民地及设施、公共服务及设施、名胜古迹、宗教设施、道路及附属物、管线设施、城市绿地、等高线九大类。
4.3.2. 新建工程
在MAPGIS6.7图像编辑模块下选择输入编辑子模块,单击左键打开。就会弹出如下图4-1所示界面。选择新建工程,然后确定。弹出设置工程地图参数界面,在该界面设置坐标系统为西安80坐标系,设置高程系统为1985国家高程系,设置地图单位为毫米,投影类型为高斯-克吕格投影。然后点击确定,完成新建工程。
图 4-1

4.3.3. 设置工作目录
在设置菜单下选择修改目录环境,此时选择
工作目录路径为:D:\PROGRAM FILES\MAPGIS67\SLIB5000
系统字库目录为:D:\PROGRA~1\MAPGIS67\CLIB
系统库目录为:D:\PROGRAM FILES\MAPGIS67\SLIB5000
系统临时目录为:D:\PROGRA~1\MAPGIS67\TEMP
4.3.4. 创立工程文件
在显示工作台单击鼠标右键,选择添加文件,在工作目录下将已经做好的DRG i49g043016.msi加载上来。
在编辑窗口点击右键,选择复位窗口,便可看到加载的图像已经被显示出来。如图4-2所示。
根据读图结果,对图中的地理要素进行分类,然后查阅<1:10000DLG地形要素分类代码表>,将所有的地图要素分类成点、线和面三种形式。然后,根据每个要素的几何属性特征建立相应的点、线或面文件。
图 4-2

4.3.5. 创立图例板
在编辑窗口空白处点击鼠标右键,选择新建工程图例,随即弹出工程图例编辑窗口。如图4-3所示。
图 4-3

根据要素几何特征设置图例类型为点类型图例、线类型图例或类型图例;
在图例信息栏中填写图例名称,图例编码,图例描述以及图例分类码信息。图例编码是指图例在<1:10000DLG地形要素分类代码表>中所对应的编号;
点击图例参数,寻找对应的地理要素符号,并根据地形图上要素的大小调节符号的X,Y系数,直至符合要求为止;
点击添加按钮,将设置好的图例添加到图例列表中,至此图例板创立完成。如图4-4所示
图 4-4

4.3.6. 地物信息采集
4.3.6.1. 关联工程图例
将光标放在左窗口空白的地方,按右键,在弹出的菜单中,选择打开图例板的功能,图例板关联路径必须是”C:\USERS\ADMINISTRATOR\DESKTOP\DLG制作\工程图例.CLN”,点击确定以后,再在左窗口空白的地方按右键,选择”打开图例文件”弹出图例板。
如果要对工程文件中的图例参数进行修改,可利用图例板进行统改。在任意一个工程文件上右击,选择”修改地图参数”,”选取全部文件”,”进行设置”,统一修改比例参数。
4.3.6.2. 矢量化前参数设置
第一,参数设置:选择”设置”菜单下面的”参数设置”,修改参数如下图:勾选”坐标点可见”和”显示线首尾点”以及”拓扑重建时搜子区。
第二,置系统参数:选择”设置”菜单下面的”置系统参数”, 设置结点/裁剪搜索半径为0.0001,坐标点间最小距离为0.0001。
4.3.6.3. 点状地物矢量化
在点编辑菜单下点击输入点按钮,弹出输入点图形对话框,在框内设置点参数。以输入工矿点为例,输入类型选择子图,输入方式选择光标定角参数缺省。
先在工程文件窗口新建点文件,命名为工矿点。然后,在工程文件窗口打开图例板,在图例板中添加工矿点图例。完成后保存设置并将图例板与文件相关联。然后点击确定,在图例板中选中工矿点,右键编辑图例参数,设置图例高度为50,宽度为40,旋转角度为0,图层为工业设施。然后单击确定,在地图上找到对应的图标位置,用鼠标制定其位置。如图4-5所示。其它点状地物矢量化方法类似。
图 4-5

输入点状地物后,如果点状地物要素的符号不符合要求能够利用点参数修改命令对符号的参数进行修改,直至满足要求即可。多个要素的属性修改能够利用快捷键Ctrl+左键,依次选中需要编辑的对象,然后进行参数修改。也能够利用图例板进行统一修改,方法在上文中已经介绍过。
4.3.6.4. 线状地物矢量化
在线编辑窗口下点击输入线按钮,弹出输入线图形对话框,在框内设置线参数。以输入陡坎为例,介绍一下线状地物要素的矢量化方法。
先在工程文件窗口新建线文件,命名为陡坎。然后,在工程文件窗口打开图例板,在图例板中添加陡坎图例。完成后保存设置并将图例板与文件相关联。在线编辑菜单下或者直接在工具条中点击输入线按钮,软件会弹出如图4-6的窗口。在窗口中根据地物要素的几何特征设置相关参数。
绘制陡坎时根据地形图上的陡坎特征,设置线型为37号图形,线形为折线,线宽为0.05mm,X,Y系数为10,颜色为蓝色,图层为坡坎层。
图 4-6

设置完相关参数后便可进行陡坎矢量化工作。如图4-7所示。其它线状地形要素的录入方法同陡坎的录入方法相同。
线状地形要素的属性修改能够在线编辑菜单下找到,也能够在工具条中直接选择。属性的修改能够经过图例板进行统一修改或者利用快捷键进行多项修改,或者进行单一修改。方法同点状地形要素的修改方法相同。
图 4-7

4.3.6.5. 面状地物矢量化
面状地物主要是指以区域的形式分布的一些地形要素。本次测区中主要有居民区、植被、公共服务及设施等。在MAPGIS中面的构成是由弧段构成的,而线和弧段在MAPGIS中是两个不同的概念。但线和弧段之间能够相互转换。因此,在MAPGIS中构面的方法有两种。一种是利用弧段直接构面,另一种,是先利用线构成面,再利用线转弧段的方法进行构面。可是两种方法都要求形成闭合的图形结构,否则,无法构成面。此时,经常要用到Ctrl+右键快捷键形成封闭图形结构。
下面以构成房子为例,介绍面状地形要素的矢量化方法。
先在工程文件窗口新建线文件,命名为房子。在该窗口打开图例板,编辑图例板,新建图例线型,设置相关参数。在这里设置线形为1号线形,线型为折线,颜色为5号色,宽度为0.05mm,图层为普通房屋层等。设置完成后保存设置,并调用图例板,将图例板与文件关联。
在工具条中点击输入线按钮,设置线型为折线。在图例板中打开设置好的线,然后根据房屋的形状进行矢量化得到房屋的轮廓如图4-8所示。
使房屋所在线文件处于编辑状态,在其它菜单下点击线转弧段命令,然后会弹出另存为一个窗口,在该窗口设置转换后的弧段文件名及路径。然后单击确定,保存文件。
在工程文件窗口添加刚才保存的文件,完成相关的拓扑错误检查。然后在其它菜单下点击拓扑重建,此时系统会给该区域赋上不同的颜色。如图4-9所示。至此,面状地形要素的矢量化工作已经完成。
图 4-8

另一种方法,是利用弧段直接构面。同样以构建房子为例,介绍用弧段构成面状要素的方法。首先,在工程文件窗口新建线文件,命名为房子。在该窗口打开图例板,编辑图例板,新建图例线型,设置相关参数,参数值设置同上。
图 4-9

在区编辑菜单下点击输入弧段,设置弧段的线型为折线,在图中找到要矢量化的房子,然后利用弧段依次描绘出房子的轮廓。其形状同上方法形同,如图4-8所示。
然后在区编辑中点击输入区,然后用鼠标在已经绘好的房子轮廓内点击,此时该区域会被激活并弹出区参数修改窗口。然后,根据需要,设置相关参数,如图4-10所示。

图 4-10

区参数的修改方法同点、线参数修改方法相同。能够统改也能够进行逐个修改。按照本次实习的要求,对所有的房屋进行统一参数修改,设置颜色为3号色,填充图案为1号图案,X,Y的系数为40,图层为3号普通房屋层。修改后的图案如图4-11所示。
图 4-11

4.4. 质量控制
为了保证DLG制作的精度及质量,矢量化过程中必须遵循一下几个重要原则:
第一、所有图上的地物都必须矢量化,不能遗漏,特殊地物要做特殊处理;
第二、光栅图面上的线条只能矢量化一次;
第三、在矢量化的同时必须打开图例板,而且要保证只能有一个图层处于可编辑状态;
第四、矢量化中以现有扫描图为标准,如果图面有错误或者表示不清楚的时候,参考其它资料;
第五、矢量化时,必须保证图面光栅文件的放大倍数为10倍以上。

5. 1:10000 DEM制作
5.1. 基础数据
DEM制作的基础数据是根据上一个步骤中生成的DLG中的等高线来制作的。确保每条等高线已经赋值。
5.2. 应用软件
DEM制作采用的软件为MAPGIS6.7。
5.3. 制作步骤
在 MAPGIS6.7 DTM 分析中,支持 Grd 和 TIN 两种类型的高程模型,根据已有等高线创立两类高程数据使用 DTM 分析 -> 处理点线菜单 -> 高程点/线三角化 或 高程点/线栅格化,高程点/线栅格化步骤如下: 图 5-1
启动 DTM 分析模块,经过文件菜单 -> 打开数据文件... -> 线数据文件,打开已有的 等高线数据,如图5-1所示。
选择处理点线菜单 -> 高程点/线栅格化,在打开的如下对话框中分别设置如下参数,
设置完毕确定退出,系统根据设置创立 Grd 数据,处理完毕提示结束。
图 5-2
网格 DX/DY:要生成的 Grd 网格数据的横向/纵向分辨率,为单个网格对应的坐标范围。与 X/Y 轴点数对应,数据范围不变的情况下,网格 DX/DY 的值越小,则 X/Y 轴点数越大, 生成 Grd 数据越大。如图5-2所示。
等高线高程属性项:选择高程字段。
栅格化保存:设置生成 Grd 文件的保存路径及名称。
查看生成的高程数据:如上图所示,经过工具栏的”打开高程数据文件”打开上一步生成的 Grd 文件,如图5-3所示。
图 5-3

6. 1:10000 DOM制作
6.1. 基础数据
DOM制作的基础数据是根据上一个步骤中生成的DEM和扫描影像制作的。
6.2. 应用软件
DEM制作采用的软件为ERDAS8.4。
6.3. 制作步骤
1、加载航空影像文件。
2、启动几何校正模块
(1)在Viewer窗口菜单条选择Raster|Geometric Correction命令,打开Set Geometric Model对话框。
(2)选择Camera选项,单击OK按钮,关闭Set Geometric Model对话框,同时启动Geometric Correction Tools对话框和航摄模式属性(Camera Model Properties)对话框。
3、输入航摄模式参数
(1)在Camera Model Properties窗口中,输入高程模型文件(Elevation Units)为(Meters);
(2)输入像主点(Principal Point)坐标为X值为-0.004,Y值为0.000;
(3)输入像机焦距(Focal Length)为152.8204;
(4)像机焦距单位(Units)为毫米(Millimeters);
(5)不选择地球曲率(Account Earth’s Curvature);
(6)空间后方交会计算迭代次数(Number of Iterations)为5;
(7)设置完以上参数后,点击Apply按钮。

4、确定内定向参数
在Camera Model Properties对话框中,选择Fiducials选项卡,进入航空摄影内定向参数的设定窗口,定义Fiducials选项卡中的类型和框标点的位置。
(1)框标类型(Fiducials Type)选择第一种,即四个角点。
(2)定义框标位置(Viewer Fiducials Locator)。单击Toggle Image Fiducials Input按钮开,,打开Viewer Selection Instruction对话框。
(3)单击ps_napp.ing窗口,图像窗口中出现一个关联框,同时打开局部放大窗口Viewer #2。
(4)在ps_napp.ing窗口左下角寻找框标点,在航空影像内定向窗口中单击Place Image Fiducial按钮,进入框标定位状态。
(5)在ps_napp.ing窗口或Viewer #2窗口中的框标点位置单击,输入第一个框标点位置,该点的图像坐标(Image X,Image Y)显示在框标数据表中。
(6)按顺时针方向重复步骤(4)和(5),直到其它三个框标点都输入到数据列表中为止。
(7)在框标数据列表中输入四个框标点的实际坐标值:
Point # Film X Film Y
1 -106.000 106.000
2 105.999 105.994
3 105.998 -105.999
4 -106.008 -105.999
当所有框标点的坐标都输入以后,Status变为Solved,同时自动计算误差(误差不得大于1)
(8)再次单击Toggle Image Fiducials Input按钮,关闭局部放大窗口。

5、设置投影参数
(1)在Camera Model Properties窗口选择Projection选项卡,进入投影设置窗口。
(2)选择Add/Change Projection,打开Projection Chooser对话框,单击自定义(Custom)选项卡。
(3)投影类型(Projection Type)为UTM。
(4)参考椭球体(Spheroid Name)为Clarke 1866。
(5)基准面名称(Datum Name)为NAD27。
(6)UTM投影分带(UTM Zone)为11。
(7)南北半球(North to South)选择North。
(8)单击OK,关闭Projection Chooser对话框,返回Projection标签窗口。
(9)选择地图坐标单位(Map Units)为Meters,同时激活Apply按钮,单击保存设置。
(10)单击Save As按钮,打开Geometric Model Name对话框,确定文件名为geomodel.gms,单击OK,关闭Geometric Model Name对话框。

6、读取地面控制点
读取地面控制点的目的是用于空间后方交会求解航片的外方位元素。由于有六个未知数,因此至少需要知道三个已知的地面控制点对,为了能够平差并有较好的精度,一般在像片四周选取6个或更多的地面控制点。
(1)在Geo Correction Tools对话框中,单击Star GCP Editor按钮,打开GCP Tools Reference Setup对话框。
(2)在Collect Reference Points Form选项中选择GCP File(.gcc),单击OK,打开Reference GCP File对话框。
(3)选择ps_camera.gcc,单击OK按钮并关闭Reference GCP File对话框。
(4)在GCP Tool对话框中,单击Solve Geometric Model图标,系统自动计算求解模型,计算中误差(RMS)、残差(Residuals)及控制点X、Y坐标值误差。
(5)在相机模型属性(Camera Model Properties)窗口中,单击Save按钮保存。
7、重采样
(1)在Geo Correction Tools对话框中,单击Image Resample按钮,打开重采样窗口。
(2)输出文件(Output Files)为geomodel.img。
(3)采样方法(Resample Method)选择三次卷积插值法(Cubic Convolution)。
(4)设置输出格网大小(Output Cell Sizes)为X值10,Y值为10。
(5)选中忽略零值(Ignore Zero in Status)复选框。
(6)单击OK,结束航空影像的正射校正。
8、数据格式转换
DOM成果的格式一般为*.tif格式,因此要将*.img格式转换为*.tif格式,并生成*.tfw文件,或者将*.img格式转换为Geotif格式。
利用ERDAS的”Import”模块进行数据格式转换。
ERDAS8.4软件中(灰度图像):

7. 基于VirtuoZo的4D产品生产
7.1. 操作步骤
本次实习以影像像对164-165、165-166为原始材料,经过全数字摄影测量系统VirtuoZo实现测区DEM、DOM、DLG的生产。现就各个产品的生产过程作详细描述。
7.1.1. 新建测区设置测区参数文件
在VirtuoZo系统中,在用户输入影像之前需要用户自定义一个文件夹,用于新建测区和存放该测区所有参数文件及中间结果、成果等。
运行VirtuoZo系统,进入到系统系统主界面,如图7-1所示。
图 7-1
在主菜单中选择”设置”、”测区参数项”,屏幕显示”打开或创立一个测区”文件对话框,在对应位置输入测区名称后进入到参数设置界面,如图7-2所示。
图 7-2
在参数设置界面设置好控制点、相机参数、加密点文件的路径。点击保存返回到系统主界面。
在设置中点击相机参数,屏幕弹出相机检校参数界面,由已知相机数据,在输入处点击鼠标左键,在VirtuoZo安装目录下找到对应的相机校检参数文件。点击确定按钮,将参数存盘。测区文件参数设置如图7-3所示。
图 7-3
在文件中点击引入,输入影像,屏幕弹出输入影像界面,在界面上点击增加,弹出添加界面,在VirtuoZo安装目录下找到对应的影像文件164和165。点击确定按钮,将影像加载上来。
在输入影像界面将像素大小参数设为:0.045mm,影像类型设置为:测量。然后点击处理按钮,系统将会对像对进行相关处理。最后点击退出,返回到系统主界面。如图7-4所示。
图 7-4

在设置中点击地面控制点,屏幕弹出地面控制点文件选择界面,由已知控制点文件,在输入处点击鼠标左键,在VirtuoZo安装目录下找到对应的控制点参数文件。点击确定按钮,将参数存盘。测区文件参数设置如图7-5所示。
图 7-5

7.1.2. 模型生成
在文件中单击打开或创立模型,新建一个模型,命名为164-165。如图7-6所示。
图 7-6

模型创立完成后,会弹出模型参数设置界面,在该界面我们需要设置左右影像。我们选择将164影像作为模型的右影像,将165影像作为模型的左影像,将核线参数设置为水平方向。如图7-7所示。点击保存,返回到系统主界面。
图 7-7

7.1.3. 模型定向与生成核线影像
7.1.3.1. 模型内定向
(1)启动内定向模块。
(2)像片内定向
内定向中,左右像片要分别进行内定向。下图中,左边窗口的中心是按钮面板,每个方块按钮对应于一个框标。单击其中一个按钮,则右边微调窗口中将放大显示其对应的框标影像。左边窗口的四周是框标影像窗口,每个小窗口显示一个框标。如图7-8所示。
图 7-8
右边窗口的上边是IO参数显示/修改窗口,可在此微调框标坐标。上半部的参数显示窗口,用来显示各框标的像片坐标、残差、内定向变换矩阵和中误差。下半部显示当前框标的放大影像。
为使内定向的精度满足作业要求,应尽量使白色的十字丝对准框标的中心。这时,要使用到框标的放大影像。具体操作是经过方块按钮选择第一个框标,然后利用右边窗口中的按钮,进行微调,直到框标放大影像中的白色十字丝对准像机的框标中心。对其它框标采用同样的方法进行调整。并保存、退出,完成左右片的内定向。
7.1.3.2. 模型相对定向
像片的内定向完成后,就要进行相对定向。
(1)启动相对定向模块。
(2)自动相对定向
相对定向时,采用自动相对定向。程序自动进行相对定向,寻找左右影像上合适的同名点,自动相对定向完成后,所有找到的同名点均以红色的”+”分别显示在左右影像上。
(3)检查和调整同名点
在界面的定向结果窗中,显示相对定向的中误差等。拉动定向结果窗的滚动条可看到所有相对定向点的上下视差。若某点误差过大,可进行调整,或将该点删除或微调。
选中某对同名点时,下面的点位放大窗口会显示出该同名点在左右影像上的图像。如果点位没有对准相同的地面位置,能够经过点位放大窗口下的按钮进行调整。选中要微调的点(将光标置于定向结果窗口该点的误差行再击鼠标左键)后,分别选择界面右下方的左影像或右影像按钮,分别点击向上、向下、向左、向右按钮,使左、右影像的十字丝中心位于同一影像点上。如图7-9所示。
图 7-9
自动相对定向完成后,在定向结果窗口检查同名点的上下视差,若上下视差较大,就把该点删除。具体操作是用鼠标左键将该点选中,然后再用鼠标左键单击窗口中的删除点按钮,即可删除该点。
编辑完成后,在影像显示窗口中的任意位置单击鼠标右键,选择弹出的右键菜单中的保存,再选择退出。
至此,模型相对定向完成。
7.1.3.3. 模型绝对定向
绝对定向前,以手工的方式在当前立体模型的左右影像上准确地定位一定数量的控制点。打开hammerIndex.html文件,查看测区的控制点分布情况,在影像上点击控制点,会出现放大图像,是量测控制点的依