基于和软件进行ads80数据生产

1、基于Xpro和ERDAS软件进行ADS80数据生产流程北京天图科技目录第 1 章处理流程1第 2 章数据准备与处理2数据路径修改22.1辅助数据整理22.2第 3 章数据预处理4第 4 章新建工程7第 5 章空三11自动生产同名点115.1量测点125.2执行空三135.3第 6 章输出 L1 级数据18第 7 章生成DEM22自动匹配点云格式高程数据227.1将点云转化为栅格 DEM247.2第 8 章输出 L2 级数据26第 9 章影像降位、镶嵌与标准分幅29影像降位299.1影像镶嵌与标准分幅299.2I第1章 处理流程整个处理按如程进行：ADS 80 RAW数据数据预处理自动点匹配与点

2、量测工程创建生产L1级数据与测图软件对接进行数字测图点云数据拼接转换匹配点云格式DSM生成DEM降位、调色、镶嵌、分幅生产标准分幅DOM生产L2级数据其中点云数据拼接转换、拼接生成 DEM，L2 级数据降位、调色、镶嵌、分幅使用 ERDASIMAGINE，其余部分使用Xpro。L1 级成果大部分国内及国外测图软件可以直接对接使用。1空三计算第2章 数据准备与处理2.1 数据路径修改ADS80 数据包含很多内容，处理前应把数据整理好，放到指定的目录。一般情况下由于拷贝移动等，使得 ADS80 数据处理时和数据时的文件路径不同，需要修改 sup，info 中的路径。用文本编辑器打开所有相关文件，批

3、量修改 IMAGE_ FILE NAME 的存放路径，使其指向正确的数据路径。2.2 辅助数据整理1. 相机文件： 将相机文件放到一个固定路径后， 需在C:Program FilesLeicaGeosystemsXProXpro.ini 中将该路径定义好，本例将相机文件夹放在 C:ADScalibrations目录下，则对 Xpro.ini 的修改如下图所示。22.资料：点坐标文件，点之记等，点坐标文件以“IDXYZ”形式的文本文件即可。3.航线信息（非必须）：航线轨迹图，飞行方向及架次等。3第3章 数据预处理1、点击开始菜单Leica GeosystemsXProLeica XPro，启动

4、XPro 程序打开 XPro 软件2、点击数据预处理图标，打开数据预处理界面3、点击 Browse图标，选择保存 ADS80 数据目录（包含 gps-imu、raw 的那一级目录），点击 Open 按钮开始验证文件，注意查看提示，如果出现 error 要根据报错信息检查数据准备阶段是否有问题。44、检查完毕后，在 Data Preparation框上勾选中 Georeferencing with 选项框，并勾选上 Pyramids，Statistics框5、在 Data Preparation框上，在Take下面勾选中需要进行数据预处理的那些ADS80 条带，在 All Bands 右侧的图标

5、点击，选择要处理的波段。5如果需要设定footprint 的路径，可以在shapefile 处选择。6、设置完成后，在Data PreparationYes框上，点击 Process 按钮，在弹出的框点击在弹出框选择Run Locally，点击 Run 执行，对ADS80 数据进行预处理。6第4章 新建工程1)在 XPro 界面上，点击新建工程图标，打开新建工程框2)点击File/ New Block，打开 New Block M框3)在 New Block M入工程名称框上点击 Browse 按钮选择新建工程的路径，在 Name 后面输4)点击Next 按钮，设置需要新建工程的投影信息75)

6、在 Coordinate System 后面选择 Projection System6)在 Choose Projection System框上选择与所做的 ADS80 数据对应的投影，点击OK，关闭 Choose Projection System框7)在 Block Preparation框上点击Block/ Add Take 按钮，加载 stp 文件88)9)选择需要处理的.stp 文件，单击 open 加载选中的数据如果有点，点击 Block > Edit Ground Control.点击左上角 Import Control Points 导入点。根据点实际投影坐标设定相关信息

7、9导入之后，可根据每个点情况修改 Type 列，修改完毕点击Save and Close10) 点击File/Save Block 按钮，保存新建的工程10第5章 空三5.1自动生产同名点1)，打开 Triangulation点击空三图标框。2)在 Triangulation框上点击File/Open Block 按钮，打开 Block Selection框3)在 Block Selection框上点击 Browse选择需要加载的工程文件114) 在 Triangulation框上点击Run APM 图标，打开 Run APM框5) 在 Run APM框上点击 OK 按钮，关闭Run APM框

8、，自动生成同名点。5.2量测点1) 在空三界面中可以看到点以三角形显示， 在需要量测的点上点击右键选择Measurement，会弹出与该点相关的所有影像122)在 Measurement 窗口处根据需要使用相关工具调节缩放和颜色，使用刺点图标添加控制点3)点击主界面的Window > Point List 可以打开Point List for GCP，查看点列表，可以看到哪些点已经刺过，也可以通过这个列表选择某个点，点击下图的图标直接跳转。5.3执行空三1) 在 Triangulation框上，点击 Create scenario 按钮，打开 New Scenario Name 对13话

9、框2) 在 New Scenario Name框上输入新的名字，点击 OK 按钮，关闭 New ScenarioName框。3) 空三参数设定，一般调节基本参数就能完成计算-基本参数: sigma0 可使用默认 3.5m；如果有点参与运算，GPS 参数可设定为Datum，若是多架次则设为 Datum per Session，在极少数情况下如果有明显的系统误差，可尝试设为 Drift； IMU，若是多架次则设为 Session，在极少数情况下如果有明显的系统误差，设为 Take，当有至少两个从相反方向飞行的重叠条带时，可设为 Sensor。4)以上设置完毕后，在 Triangulation框上，

10、点击Run Adjustment 按钮，运行空三5)空三完毕可以看到如下图的整体误差信息，一般红域为需要进行修正的区域，可以通过调整连接点或是匹配点的方式进行修改。146)在 Triangulation框上，选中粗差较大的点，单击右键选择 Measurement，在该界面的下方会显示该点所在影像的多个视窗，调整位置偏差较大的影像的点，也可以选中粗差较大的点，单击右键选择Remove Point 按钮，删除粗差较大的点。7)也可通过点列表浏览，查看红色的点，进行重新量测或调整。调节完可以再次执行空三，查看误差。158)若点调节完后，误差仍难以满足要求，可以尝试使用 Triangulation-A

11、nalysis-APM loop功能，在基本参数设定界面，勾选 Densify Points Threshold，设定循环次数，并根据需要设定高级参数中的 Densify Points Threshold 阈值。执行空三即可。9)若空三误差已经满足要求，点击 Write to Block 按钮，将空三后的结果写入工程中。10) 若需要导出点，点击Scenario 下的 Export Ground Points16两个文件会被生成：adjustedGroundControlPoints.txt、adjustedGroundTiePoints.txt.17第6章 输出 L1 级数据1)点击Prod

12、uct Generation 图标，打开 Product Generation框，点击File/OpenBlock 按钮，选择空三后的工程文件2)在 Product Generation框上点击 Generate 图标，打开 Generate Product框3)在 Generate Product框上选中 L1_RGB，点击 Next 按钮184)在 Generate Product框上选中Atmospheric + BRDF 选项，点击Next 按钮5)6)在 Generate Product在 Generate Product框上选择 Standard 选项，点击Next 按钮框上点击O

13、utput Root 后面的按钮，设置输出文件的路径，在 Perfix 后面设置输出影像的前缀。197)在 Generate Product框上，点击 Finish 按钮，关闭 Generate Product框8)在框上点击 Execute jobs 图标，打开 Process Locally or Remotely框9)在 Process Locally or Remotely框上点击Run Locally 按钮，运行程序，此成果输出后即可得到ADS80 的 1 级RGB 数据。10) 如果需要 1 级 PAN 数据，则需在Edit 点击 Products2011) 在 Product T

14、ype Editor 中新建 L1_PAN，参考 L1_RGB 设定相关参数，保存。在前面数据生产中选择新建的 L1_PAN 数据类型即可获得 1 级 PAN 数据。21第7章 生成 DEM7.1自动匹配点云格式高程数据1)点击Product Generation 图标，打开 Product Generation框，点击File/OpenBlock 按钮，选择空三后的工程文件2)在 Product Generation框上点击 Generate 图标，打开 Generate Product框223)在 Generate Product框上选中 DSM_RGB，点击 Next 按钮4)在 Gen

15、erate Product框上选中Atmospheric + BRDF 选项，点击Next 按钮5)6)在 Generate Product在 Generate Product框上选择Standard 选项，点击Next 按钮框上点击Output Root 后面的按钮，设置输出文件的路径，在 Perfix 后面设置输出影像的前缀。237)在 Generate Product框上，点击 Finish 按钮，关闭 Generate Product框8)在框上点击 Execute jobs 图标，打开 Process Locally or Remotely框9)在 Process Locally o

16、r Remotely出后即可得到点云格式的高程数据。框上点击Run Locally 按钮，运行程序，此成果输7.2将点云转化为栅格 DEM1) DEM 提取完成后，需要使用 ERDAS IMAGINE 将点云转换为栅格格式 DEM，才能继续进行。242)点击 ERDAS IMAGINE 主界面 Terrain下的 TPT(Terrain Prep Tool)进行合并，点击 TPT 中图标加载匹配好的 las 文件3)点击下进入 Surface DTM 界面，设置 Rasterization 选项下栅格 DEM 的文件名、路径和分辨率，点击 OK 执行。4)输出的结果就是栅格格式的 DEM。25

17、第8章 输出 L2 级数据1)点击Product Generation 图标，打开 Product Generation框，点击File/OpenBlock 按钮，选择空三后的工程文件2)在 Product Generation框上点击 Generate 图标，打开 Generate Product框3)在 Generate Product框上选中 L2_RGB，点击 Next 按钮264)在 Generate Product框上选中Atmospheric + BRDF 选项，点击Next 按钮5)6)在 Generate Product在 Generate Product框上选择Standa

18、rd 选项，点击Next 按钮框上设定 GSD 地面分辨率；点击 Elevation Source 中 Selecta DEM，选择上一步转换好的栅格 DEM；点击 Output Root 后面的按钮，设置输出文件的路径，在 Perfix 后面设置输出影像的前缀。277)在 Generate Product框上，点击 Finish 按钮，关闭 Generate Product框8)在框上点击 Execute jobs 图标，打开 Process Locally or Remotely框9)在 Process Locally or Remotely框上点击Run Locally 按钮，运行程序，

19、此成果输出后即可得到ADS80 的 2 级正射的 RGB 数据。28第9章 影像降位、镶嵌与标准分幅9.1影像降位1)由于默认输出的为 16 位的数据，为了生产需要将数据进行降位处理2)点击 ERDAS IMAGINE 主界面 Raster 选项卡下 Radiometric中的 Rescale 进行降位，设置输入输出数据，算法为 Minimum-要批处理，点击Batch 进行。um，输出统计忽略 0，点击 OK 即可，如果需9.2影像镶嵌与标准分幅1)点击 ERDAS IMAGINE 主界面 Toolbox 下 MosaicPro，选择降好位的数据，点击 Use Entire Image 按钮（如果影像背景区域较大可选择 Compute Act