PDF《使用ERDAS软件处理ADS80数据》

使用ERDAS软件处理ADS80数据 生成DEM和DOM 北京天图科技有限公司 I 目录

一、 前期准备.

1

二、 提取 DEM.2

三、 正射纠正.8

四、 裁切、降位.9

五、 影像镶嵌.11

六、 遇到的问题和建议.16

1

一、 前期准备

1、 确保 sup 所对应路径等正确, 使用批处理创建影像金字塔(3*3)及进行统计运算(统计选 项中 Skip Factor X、Skip Factor Y 务必设为 1),在viewer 中打开浏览确保金字塔无误;

2、 将计算好的七参数和2000 椭球信息,加进…\Intergraph\ERDAS IMAGINE 2014\etc\spheroid.tab 文件;

七参数为加入高程值计算所得,因此无需定义水准精化 bin 文件.

2

二、 提取 DEM

1、 点击 ERDAS IMAGINE 主界面 Toolbox 下IMAGINE Photogrammetry, 创建工程, 选择 Mixed Sensor 模型;

设置水平和垂直投影坐标信息;

2、 导入全色前后视 sup 文件,保存工程;

3、 启动 eATE 模块, 点击 设置 Overlap 标签下 Maximum Rays 为2, 不匹配航带间重叠 区;

3

4、 点击 设置 General 标签下 Stop at Pyramid Level 为1,Point Sampling Density 为4, Pixel Block Size 为1000,Threads 为8;

5、 点击 设置 General 标签下工程路径, Output Files 标签下 Output Points in Raster Format

4 路径(栅格格式 DEM,可选,后续可用 TPT 手动生成)和Nongridded Format 路径(点云 格式 DEM,可选,供编辑)

6、 点击 生成处理处理配置文件

5

7、 点击 进入批处理界面,点击 Submit,根据计算机 CPU 线程数设置并行数;

8、 匹配过程进度可通过 eATE 主界面和 ERDAS IMAGINE 的Process List 看到

6

9、 DEM 提取完成后,如果合并生成栅格文件失败,使用 ERDAS IMAGINE 主界面 Terrain 标 签下的 TPT(Terrain Prep Tool)进行合并,点击 TPT 中 图标加载匹配好的 las 文件 点击 下进入 Surface DTM 界面,设置 Rasterization 选项下栅格 DEM 的文件名、路径和分 辨率,点击 OK 执行.

7 8

三、 正射纠正

1、 因为本次仅需要 ADS 的RGB 正射成果,因此使用 IMAGINE Photogrammetry,创建一个 新的工程,按照全色工程类似设置将 ADS 的RGB 数据加入,点击 进入重采样界面,选择DTM Source 为DEM,设定为上面提取的 DEM,设定分辨率,点击 Add Multiple 设定路 径,将所有条带加入,点击 Batch 执行正射重采样并行批处理.

9

四、 裁切、降位

1、 根据作业区裁切正射好的条带数据,使用 Raster 选项卡下 subset 工具,可根据坐标范 围、AOI 等方式进行裁切.

2、 点击 Raster 选项卡下 Radiometric 中的 Rescale 进行降位,设置输入输出数据,算 法为 Minimum-Maximum,输出统计忽略 0,点击 OK 即可,如果需要批处理,点击 Batch.

10 11

五、 影像镶嵌

1、 点击 ERDAS IMAGINE 主界面 Toolbox 下MosaicPro,选择降好位的数据,点击 Use Entire Image 按钮(如果影像背景区域较大可选择 Compute Active Area) ,点击 OK 关闭 Add Image 对话框,同时弹出 MosaicPro 视窗,数据被加载进去.

2、 点击 ,打开 Seamline Generation Options 对话框.选择一种拼接线生成方法,如Weighted Seamline,点击 OK 生成拼接线.

3、 点击 按钮进行拼接线编辑.

12

4、 点击 ,可以选择直方图匹配算法,设置提前调好颜色的模板数据进行匹配,达到匀 色目的.

5、 点击 ,设置五到十个象元平滑、羽化,使接边处过渡自然.

13

6、 点击 按钮,可以设置输出影像范围,数据类型,分辨率,波段信息. 如果需要进行标准分幅,可以选择 Method 下拉条中的 ,选 择标准分幅坐标文件,进行分幅.坐标文件格式为文本,每条记录格式如下: 图幅号 X1 Y1 X2 Y2(中间以空格分隔,坐标值分别代表左下角、右上角坐标) 如:K45E001055