标签：2021.5 21 Pix4d EPSG AUTHORITY filename DSM WGS
航摄像片生成DSM（以Pix4d软件为例）2021.5.21
1、常用无人机数据处理软件
1.1、ContextCapture
1.2、Pix4d
2、Pix4d安装步骤
2.1 Pix4d安装包下载
到百度网盘下载链接中下载Pix4d安装包文件，提取码：mw66
下载到本地文件夹后如下图所示：
2.2 Pix4d安装方法
首先双击Pix4Dmapper.msi安装包文件，弹出安装界面，点击Next
选择安装路径为 D:\Pix4Dmapper\ 后点击Next
之后点击Install开始安装，安装结束后取消勾选 Launch Pix4Dmapper后点击Finish即可。
下面利用Pix4d文件夹中的破解工具进行破解，双击硬件获取.exe后弹出窗口坐下角输入123后点击右侧的Generate HWID后在上方的文本框中生成了硬件ID码，复制下来。
之后打开注册机文件夹下的注册.exe，将HWID码粘贴到左侧文本框中后在右侧Owner中输入123，Days输入9999后点击Decode后会在右侧生成log日志，最终点击最右侧的Save REG两次即可在注册机文件夹生成123_pix4d.reg注册表文件。
现在双击123_pix4d.reg注册表文件，看到123_pix4d.reg中包含的项和值已成功添加到注册表中后点击确定即可，至此Pix4d破解完成。
最后双击Pix4d文件夹下的联网工具.exe后在右下角弹出Pix4d Tool v4.0，点击AutoStart即可，启动桌面或开始菜单的Pix4Dmapper.exe便可打开Pix4D。
进入上述Pix4d界面表示破解成功，可以新建项目了。（注：如果没正常启动该界面，需要多次尝试注册机破解后并注意打开联网公工具中的AutoStart后打开Pix4Dmapper.exe）
3、准备无人机倾斜摄影测量数据（航摄像片）
外业无人机倾斜摄影测量数据如下图所示。
4、利用Pix4d生成DOM和DSM
4.1 新建项目开始处理
打开Pix4Dmapper后，点击新项目，项目名称输入TestUAVGerateDOMandDSM后点击Next
在弹出的新项目窗口中点击添加图像，在图像浏览窗口选择无人机倾斜摄影测量的数据（54张航摄像片）后点击Next
在图片属性窗口会显示图像的地理位置、基准和相机型号等信息，直接点击Next
在选择输出坐标系窗口可以看到自动检测到WGS 84/ UTM zone 44N，而非纯粹的地理坐标系WGS 84，处理的结果如果要用WGS 84地理坐标系的话需要进行坐标转换。
选择处理选项模板为 3DMaps，右侧可以看到为应用实例生成DSM和正射影像图，选中开始处理后点击Finish。
下面进入数据处理阶段，只需等待
4.2 Pix4d处理航摄像片生成的DSM和DOM
最后项目处理结果文件夹如下图所示，数字表面模型DSM和数字正射影像DOM结果都被保存在3_dsm_ortho文件夹下，其中DSM被保存在1_dsm文件夹，DOM被保存在2_mosaic文件夹。
DSM正射影像对应的空间参考投影文件TestUAVGerateDOMandDSM_dsm.prj内容为：
PROJCS["WGS 84 / UTM zone 44N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",81],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32644"]]
DOM正射影像对应的空间参考投影文件TestUAVGerateDOMandDSM_transparent_mosaic_group1.prj内容为：
PROJCS["WGS 84 / UTM zone 44N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",81],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","32644"]]
4.3 ArcGIS打开DSM和DOM结果查看坐标系
打开ArcMap后，点击添加数据，选择DSM和DOM后确定，显示效果如下图所示，同时在图层上右键->属性->源可查看对应的空间参考（坐标系）为WGS 84 / UTM zone 44N。
5、DSM结果tif投影坐标转换为地理坐标WGS84
在ArcGIS中导出地理坐标系WGS84对应的空间参考文件WGS 1984.prj
WGS 1984.prj内容如下：
GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433],AUTHORITY["EPSG",4326]]
5.1 使用ArcGIS中的ToolBox（ArcMap手动操作）
在ArcGIS的工具箱ToolBox中找到数据管理工具下的投影与变换，进入栅格中的投影栅格，打开投影栅格窗口，选择要重投影的栅格文件和输出坐标系WGS84地理坐标系后确定即可。
在ArcMap中查看重投影后DSM84.tif和DOM84.tif的坐标系
DSM84影像中通过识别工具查看影像像素值发现黑色的边框部分对应的RGB的值均为0，可在图层上右键选择属性，在符号系统下的RGB合成中选中显示背景值左侧的方框，并选择背景颜色为白色则可去除黑色边框。
5.2 后台调用ArcGIS中的Arcpy实现栅格重投影并处理NoData（python脚本）
5.2.1 栅格重投影坐标系arcpy脚本所用工具介绍
这里主要利用ArcGISToolBox中的投影栅格和复制栅格两个工具，通过查看工具帮助里的arcpy调用API来实现python脚本对DSM和DOM的重投影，最终编写了一个ConvertTifToWGS84.py脚本文件，放到ArcGIS对应的python编译器路径（D:\ArcGIS\Python27\ArcGIS10.4）下打开cmd窗口输入python.exe ConvertTifToWGS84.py即可运行。
5.2.2 ConvertTifToWGS84.py调用示例
如上图所示，要转换的栅格影像Tif所在的文件夹为C:\Users\Administrator\Documents\pix4d\TestUAVGerateDOMandDSM\3_dsm_ortho\DSMandDOMToWGS84，需要手动在此文件夹下新建一个result文件夹，这个result文件夹用来保存转换为WGS 84坐标系后的Tif文件，同时在D:\ArcGIS\Python27\ArcGIS10.4文件夹下修改ConvertTifToWGS84.py代码的第三行为rootdir = "C:/Users/Administrator/Documents/pix4d/TestUAVGerateDOMandDSM/3_dsm_ortho/DSMandDOMToWGS84"后内容如下所示。
ConvertTifToWGS84.py（利用ArcGIS 的python2.7 arcpy）
import os
import arcpy
rootdir = "C:/Users/Administrator/Documents/pix4d/TestUAVGerateDOMandDSM/3_dsm_ortho/DSMandDOMToWGS84"
flag = 0
for (dirpath, dirnames, filenames) in os.walk(rootdir):
for filename in filenames:
suffex = filename[len(filename) - 4] + filename[len(filename) - 3] + filename[len(filename) - 2] + filename[len(filename) - 1]
if (suffex == '.tif'):
inputTif = rootdir + '/' + filename
middleTif = rootdir + '/result/' + '' + 'wgs1984' + filename
outputTif = rootdir + '/result/' + '' + 'Noda_wgs1984' + filename
print(inputTif)
filename_prj = inputTif.replace('.tif','.prj')
#sourceSR = arcpy.SpatialReference(filename_prj)
print(filename_prj)
targetSR = arcpy.SpatialReference("C:/Users/Administrator/Desktop/WGS 1984.prj")
arcpy.ProjectRaster_management(inputTif, middleTif, targetSR, "BILINEAR")
# arcpy.ProjectRaster_management(inputTif, middleTif, targetSR, "BILINEAR", "#", "#", "#", sourceSR)
arcpy.CopyRaster_management(middleTif,outputTif,"#", "#", 0);
print('Result:'+outputTif)
print('successed!')
# python.exe ConvertTifToWGS84.py
GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433],AUTHORITY["EPSG",4326]]
下面需要将WGS 1984.prj复制到桌面，在D:\ArcGIS\Python27\ArcGIS10.4文件夹下打开cmd命令行窗口，并输入命令python.exe ConvertTifTiWGS874.py回车开始执行坐标转换，执行完成后会在result文件夹得到转换后的DSM和DOM。
result文件夹下转换为WGS 84坐标系后的DSM和DOM结果如下，可以在ArcGIS中打开查看坐标系。
5.2.3 成功叠加WGS84坐标系的DSM和DOM到三维地球OSGEarth显示
标签：2021.5,21,Pix4d,EPSG,AUTHORITY,filename,DSM,WGS 来源： https://blog.csdn.net/jing_zhong/article/details/116789576
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。