Agisoft Metashape 照片高程改正 附python源码

Agisoft Metashape 照片高程改正 附python源码

---

前言

Agisoft Metashape可以直接读取照片中的经纬度和高程信息，无人机获取的照片中保存的都是WGS84经纬度坐标和大地高，而我们实际测量中用到的为国家2000高斯投影坐标和85高程，可以通过我往期博文（Metashape 坐标系选择 坐标转换）实现WGS84经纬度坐标转投影坐标，本文讲解通过python脚本，通过高程改正数实现大地高转85高程（小面积区域），大面积可以通过七参数的方式实现大地高转85高程，从而实现将照片的WGS84经纬度坐标和大地高转到地面测量坐标系的投影坐标和85高程与控制点坐标系和高程一致，方便后期的刺点，同时这对于提高成果精度很有必要。

一、计算高程改正数

可以通过测量获取的像控点的原始大地高和转换后的85高程求解对应的高程改正数，在面积不大（小于10平方千米）的平原地区，一般情况下高程差值（高程异常）不大，航测地形图项目可以通过加减高程改正数得的方式实现大地高转85高程。

测区有控制点的情况下，可以使用求七参数的方法，实现大地高转85高，没有控制点则需要去当地的省/市CORS中心进行转化，RTK测得的高程为大地高程，1985国家高程基准属正常高，CORS不能直接测85高程，需要转换。

二、python源码

将下面代码复制、粘贴到文本文件中，将后缀名改成.py,即可在A