李国春

GF7 DEM/DSM制作让我纠结了好一阵子了,总希望通过不进行人工干预达到一些网站展示的处理结果。后来我发现原来人家那是手工编辑后的结果。因为这种技术路线总体上就有它的局限性,想不进行人工编辑直接拿出来就有完美的结果不太现实。大家发表出来的结果大都是手工编辑后的,也难怪,人们总是会倾向于展示自己认为满意的结果。:)

RSD制作DEM/DSM数据的方法在《高分7(GF7)卫星数据制作DEM》已经介绍过,这里再介绍一下近期针对平原区域求精完善的过程。还有就是将一些主要计算的中间结果和处理编辑过程都交给用户自行处理,用户可以根据本地地形地物特点编辑处理以得到更为精确的本地DEM产品。

RSD 3.3.1针对平原地区地貌,重新设计了DEM制作过程。其实只是增加了一些参数给用户处理更大的灵活性,并没有太大的本质上的改进。山区DEM制作容后续分析评价。

说明:本文介绍的处理结果均无人工编辑。本人时间比较紧,一个人的团队又没有其它成员协助,没时间顾及这些重复劳动的工作。以前有介绍基于对象的高程编辑,又有了一点改进,后续再介绍如何进行手工交互编辑。

一、原理

根据GF7前后视探测成像原理,在核线上检测像点位移得到前后视视差,构造出视差图,再使用视差图计算得到逐像点的高程。像点位移是发生在沿核线方向的,理论上在与核线垂直的方向上前后视没有视差。但是观察GF7数据,前后视卫星方位角并不在180°直线上,再加上二者卫星探测天顶角不同,需要纠正成直线核线形成核线图,然后进行核线匹配得到视差图。

这样需要对前后视图像先刺点得到一系列的同名点,依据这些同名点建立一个核线方程,生成核线图像。在核线图像上就基本能够满足像点位移只发生在沿核线方向而垂直核线的方向上基本没有视差的条件。考虑到实际存在探测角度的不一致和核线图的误差,在与核线垂直方向给与一定的容差,而不是使用单一的一条核线数据,以提高计算结果的稳定性。

RSD计算输出的高程使用基准高程进行校准。实际上计算连接点时,所有连接点用一个二维平面进行了拟合,所有视差都是在这个平面上的位移,根据视差所算得的高程也都是相对这个平面的相对高程,有正有负。这个平面暂且称之为相对高程平面。转换成真实高程时,需要根据该范围真实高程在该范围拟合另一个平面,称为绝对高程平面。以两个平面均值点为基准,计算两个平面逐像点的高程差。这个方法在平原地区由于两个平面基本是平行的,这个高程差可能只用一个常数就可以。在山区由于TIE POINT 采集分布原因会造成相对高程平面倾斜,与绝对高程平面不平行发生相交,需要逐点计算高程差,然后逐点修正高程。

实际上高程在正射过程使用了基准高程进行了正射校正,这样参考基准高程平面换成了参考高程的高程曲面。高程计算结果换成了在该曲面上的校正结果。由于是平原地区,可以只使用比较权威的通用高程(例如SRTM90等)做基准,然后用实测高程与计算结果进行线性拟合修正。用户如果有自己的高程数据当然也可以用来做基准高程修正计算结果。

二、DSM制作在RSD平台的交互操作

打开RSD 3.3.1,找到图1的菜单

图1 启动GF7 DEM制作菜单

菜单中 DEM生成-核线约束NCC匹配DEM生成-核线约束Census匹配 仍然是早期介绍的菜单选项,这里推荐使用 DEM生成-NCC匹配平原模式。虽然这里强调是平原模式,但除了高程处理等细节外山区使用也是一样的,具体我后面再分析评价。

点击DEM生成-NCC匹配平原模式后出现图2。点击输入编辑框后面的小按钮,依次输入GF7的前视和后视文件。注意,后视数据使用单通道PAN数据,不要输入MUX。

图2 输入前视和后视文件的对话框

点击确定后,出现图3的对话框。这个对话框比前文《高分7(GF7)卫星数据制作DEM》有了一些变化。

图3 平原低目标高程参数设置对话框

图3中,左上和左中两个框内分别为GF7后视和前视基本信息,这些信息多是从卫星数据集的属性数据中自动获取的,或者使用的是已知数据,通常不需要用户更改。

左下框中为输出选项,其中高程是我们需要的最终结果。高程分第一次匹配高程和第二次匹配高程,这是为了提高配准精度设计的。如果不需要二次匹配,可以在右侧框中第二搜索长度、第二匹配长度、第二容差或第二步长中任何一个输入0,就禁止了二次匹配。

右侧框的前4个参数是位置信息,是本次处理的数据在后视数据中的位置。例如,本例中后视数据有35818行×40001列数据,本次处理从第3500行、4000列开始,处理9000行×9000列的数据。可以调整处理数据的起点和范围来处理数据集中任意位置的数据。注意:一些未注册用户不能任意扩大处理范围。实际上RSD推荐分块处理,逐渐更改处理起点和范围遍历整个数据集处理效率更高。

右侧框的其它值是控制处理过程和精度的参数,后面要对这些参数分别进行详细讨论。

在此图3对话框中点击确定,出现图4的分裂窗。

图4 选择 TIE POINT 点的对话框

必须在后视和前视图像上创建一些连接点(TIE POINT),通过两张图上足够数量的同名连接点来建立核线方程和创建核线图像。

图中分裂窗左侧是后视图像,右侧是前视图像。一对连接点需要在这两张图象上获取相同位置的点。方法是先在左侧图像选取位置,光标放在该位置右击鼠标,出现图4中间的弹出菜单。第一行菜单命令出现光标的坐标位置,点击该菜单项,窗口中出现一个小十字,表明该点被选中。在右窗找到相同位置,同样右击鼠标,也点击菜单,窗口中也出现一个小十字。此时在图4最左侧的列表中,出现了一对坐标数据。表明一对连接点选择成功。

重复上述过程,选择足够的连接点。使用连接点建立的核线方程是一个二元一次的线性式,不需要很多的连接点。通常在图像四角和中间散布选择几个连接点就够用了,但是为了稳定起见,多选择十几二十几个连接点可能更好。

图像比较大时,在两张图找到同名点有时需要滚动图像,不容易找到。一个小窍门是在弹出菜单选择 Full Image Display,这时显示整个图像的全貌,再右击需要显示在窗口的点,选择 1:1 Real Pix Display,这时你感兴趣的位置显示在窗口。再用同样的方法在右图显示同一个感兴趣的位置。这样两个窗口都显示了同一个位置,便于选择连接点。

连接点准备完毕以后,可以按 保存数据 按钮将数据保存为一个文本文件,下次可以打开这个文件继续使用。

开始 按钮,软件开始自动建立核线方程、创建核线图像和视差图等工作。

更改图3中的位置信息,遍历整个图像,得到整个图像的分块处理结果。

三、输出的结果/中间结果

一个块的高程计算一共10层输出,这里禁止了两层像点位移的中间结果,仅显示8层。

第1层,第一次配准高程。绝对高程。

图5 第一匹配高程

第2层,第一次配准可信度。这是与第1层高程数据对应的可信度,高可信度达0.95以上,低值仅0.2左右。在图3设置可信度参数时已经将可信度低于0.4的对应高程删除掉了(用0.0001填充)。用户也可以根据需要自行设置可信度门限,来决定保留或者删除不可靠的高程数据。注意:高可信度也仍然可能存在像点错误匹配,所以,这个可信度也有可能不可信,仅供参考。见图6。

图6 第一匹配可信度

第3层是第二次配准结果的高程,图7。

图7 第二次配准高程

第4层是第二次配准结果的可信度,图8。

图8 第二次配准可信度

第5层后视正射图像,图9。

图9 后视正射图像

第6层核线图正射图像,图10。

图10 核线图正射图像

第7层,合成可信度,图11。

图11 合成可信度

第8层,合成高程,图12。

图12合成高程结果

加企鹅758461012,原来的满了。

高分7(GF7)卫星数据制作平原地区DSM/DEM(一)相关推荐

  1. 高分7(GF7)卫星数据制作平原地区DSM/DEM(三)——可信度分割

    李国春 (接上) 五.高程可信度分割 很多人纠结高程数据的空间分辨率大小.高空间分辨率固然是一个重要指标,但是可靠数据所占比例也同样重要.RSD可以根据不同的可信度标准对高程数据进行分割. 首先勾选高 ...

  2. 高分7(GF7)卫星数据制作DEM

    RSD平台新增利用GF7卫星前视PAN数据(FWDPAN)和后视PAN数据(BWDPAN)制作数字高程模型/数字表面模型(DEN/DSM).实际上多光谱后视数据(BWDMUX)的参与能进一步提高前后视 ...

  3. 高分1(GF1)、高分2(GF2)卫星数据大气校正

    KEYWORDS:GF1,GF2,RSD,大气校正,遥感软件 0. RSD大气校正 RSD是李国春教授团队开发的一款遥感数处理软件.其大气校正模块是参照USGS LaSRC大气校正流程,使用VC++重 ...

  4. 使用RSD对高分1号卫星数据进行批量大气校正

    高分数据处理可以有下述5中选项 MSS数据的批量正射校正 MSS数据批量正射校正和大气校正 PAN 数据的批量正射校正 MSS与PAN批量正射校正和融合 MSS与PAN批量正射校正.MSS的大气校正和 ...

  5. 使用RSD批量正射校正和批量融合高分1号卫星数据

    这里提供菜命令的方式启动批量处理.菜单批量处理进行以下几种处理工作. MSS数据的批量正射校正 MSS数据批量正射校正和大气校正 PAN 数据的批量正射校正 MSS与PAN批量正射校正和融合 MSS与 ...

  6. 资源3(ZY3)号卫星三线阵数据制作DEM

    资源3号卫星有正视和前.后视3线阵数据,应用这些数据可以制作高程数据DEM/DSM.上一篇<资源3(ZY3)号卫星NAD-FWD数据制作DEM>介绍了在RSD使用正视(NAD)和前视(FW ...

  7. 资源3(ZY3)号卫星NAD-FWD数据制作DEM

    资源3号卫星有正视和前.后视3线阵数据,应用这些数据可以制作高程数据DEM/DSM.使用同轨3个线阵数据制作DEM请参阅<资源3(ZY3)号卫星3线阵数据制作DEM>.这里仅介绍在RSD平 ...

  8. 高分二号(GF-2)号卫星数据的查询下载地址和方法

    高分二号卫星是我国自主研制的首颗空间分辨优于1米的民用光学遥感卫星可在平台中查询到,搭载有两台高分辨率1米全色.4米多光谱相机,具有亚米级空间分辨率.高定位精度和快速姿态机动能力等特点,有效地提升了卫 ...

  9. 遥感影像云检测-现有光学遥感卫星数据参数合集

    现有光学遥感卫星 1.MODIS卫星介绍 1.1.MODIS传感器简介 1.2.MODIS传感器基本参数 1.3.MODIS 数据特点及波段分布特征 1.4.MODIS 数据产品产品级别及类型 1.4 ...

最新文章

  1. Oracle如何代码编辑,配置UltraEdit为Oracle PL/SQL代码编辑器
  2. IBatis.net动态SQL语句
  3. MediaCodec在Android视频硬解码组件的应用
  4. PHPCMS V9.3.2用户注册模板中的一个低级Bug
  5. 使用idea编写消费者,接收生产者的持续日志输出【小案例】(三)
  6. [导入]相片: 23853741.jpg
  7. Java项目开发全流程实践
  8. 深入ASP.NET 2.0的提供者模型
  9. Template Method (模板方法模式)
  10. 计算机网络 — 网络层
  11. 遗传算法原理案例及MATLAB代码
  12. node中封装MongoDB
  13. 有效 TCP RST
  14. git更换主分支master
  15. 京东架构师讲解京东是如何优化列表页的
  16. 神经网络和深度神经网络,深度神经网络类型包括
  17. [iOS开发]Status Bar Style
  18. 怎样把其他格式的音频文件转换成MP3格式的
  19. 使用cephadm部署单节点ceph集群,后期可扩容(基于官方文档,靠谱,读起来舒服)
  20. 利用存储过程分页显示查询数据

热门文章

  1. Linux正则表达式使用方法详解
  2. 疫情大考:钉钉等移动办公们的短板与变局
  3. linux mysql开远程访问,Linux mysql开启远程访问
  4. macmini和imac哪个好?
  5. Android L Camera2 API sample ver1 - startPreview
  6. c语言左移右移运算符详解
  7. 捷豹改激光大灯 捷豹XF氙气灯升级LED灯作业
  8. 软件测试人员的职业规划 西安尚学堂
  9. 数据可视化之Echarts
  10. flash在线剪切缩放批量修改程序