arcpy投影(一)

2023-12-29 14:20| 来源: 网络整理| 查看: 265

arcpy投影这一个专题从文件位置、文件含义、空间参照获取、转换关系查询、投影定义、自定义转换关系、投影变换这几个角度上系统的进行了介绍，整理出了：

arcpy投影(一)——prj、gtf文件定义、路径及解析（arcmap,arcpro）

arcpy投影(二)——基准面变换概念及参数、空间参考对象获取、变换关系获取方法梳理与解析（Spatial Reference、ListTransformations）

arcpy投影(三)——定义投影、地理变换关系自定义和投影变换Project_managemen（含基准面/椭球体转换参数使用方法，arcpro/arcmap）三篇文章，本文为第一篇：

一、投影坐标prj文件，转换gtf文件

ArcGIS的prj文件是一个专门用来描述坐标系投影的文件，里面定义了坐标系的名称、长短半轴、基准面、投影参数等一些其他重要的信息。 gtf文件是专门用来描述坐标转换关系的文件，在你自定义转换关系（Create Custom Geographic Tramsformation）后均会在默认路径下生成，该文件中包含了转换关系名称、原始坐标系名称及相应参数、目标坐标系名称及相应参数，以及两个坐标系统转换的方法与参数。

二、文件路径

无论是ArcMap还是Arcpro，这两种文件都是放在C盘用户文件夹在的AppData文件夹中，要注意的是AppData文件夹为隐藏文件夹，如果找不到的话可以勾选文件夹“查看”中的“隐藏的项目”选项。

（一）ArcMAP路径

ArcMap中自定义坐标系和投影关系分别放在Coordinate System和CustomTransformations文件夹中

(1)自定义坐标系

若Coordinate System文件夹不存在，则在第一行路径下进行创建，若存在prj文件直接放入到文件夹中

C:\User\Administrator\AppData\Roaming\Esri\Desktop10.6\ArcMap\ C:\User\Administrator\AppData\Roaming\Esri\Desktop10.6\ArcMap\Coordinate System （2）转换关系

同上

C:\User\Administrator\AppData\Roaming\Esri\Desktop10.6\ArcToolbox\ C:\User\Administrator\AppData\Roaming\Esri\Desktop10.6\ArcToolbox\CustomTransformations （二）ArcPro路径 (1)自定义坐标系

在pro中进行了测试，新建Coordinate System文件夹等方式都不可以识别到，通过收藏的方式可以进行识别，自己弄的投影可以放到这个路径下，但是经测试目前pro在这个路径下只能识别出来全英文名的prj文件，如果使用Arcpy进行处理时，可以直接通过路径进行指认。

C:\User\Administrator\AppData\Local\ESRI\ArcGISPro\Favorites （2）转换关系

同ArcMap,注意如果要使用Arcpro的Arcpy进行投影关系的变换，一定要检查该路径下是否已经存在对应的转换关系，没有的话要进行生成或从他处拷贝，否则函数会报参数错误，这部分后面的文章中会提到。

C:\User\Administrator\AppData\Roaming\Esri\ArcGISPro\ArcToolbox \ C:\User\Administrator\AppData\Roaming\Esri\ArcGISPro\ArcToolbox\CustomTransformations 三、文件内容解析（一）投影坐标系prj文件 PROJCS[ "Ain_el_Abd_UTM_Zone_38N" ,GEOGCS[ "GCS_Ain_el_Abd_1970" ,DATUM["D_Ain_el_Abd_1970" ,SPHEROID["International_1924",6378388.0,297.0]] ,PRIMEM["Greenwich",0.0] ,UNIT["Degree",0.0174532925199433] ] ,PROJECTION["Transverse_Mercator"] ,PARAMETER["False_Easting",500000.0] ,PARAMETER["False_Northing",0.0] ,PARAMETER["Central_Meridian",45.0] ,PARAMETER["Scale_Factor",0.9996] ,PARAMETER["Latitude_Of_Origin",0.0] ,UNIT["Meter",1.0] ,AUTHORITY["EPSG",20438] ]

这是一个投影坐标系。 PROJCS：代表的是其投影坐标系。投影名称是：Ain_el_Abd_UTM_Zone_38N 基准面信息：D_Ain_el_Abd_1970 椭球体：SPHEROID[“International_1924”,6378388.0,297.0]-长半轴和反扁率单位：UNIT[“Degree”,0.0174532925199433] PARAMETER：代表投影重要的参数信息

（二）地理坐标prj文件 GEOGCS[ "GCS_Beijing_1954" ,DATUM[ "D_Beijing_1954" ,SPHEROID["Krasovsky_1940",6378245.0,298.3] ] ,PRIMEM["Greenwich",0.0] ,UNIT["Degree",0.0174532925199433] ,AUTHORITY["EPSG",4214] ]

这是一个经纬度坐标系。 GEOGCS：代表是其地理坐标系基准面：D_Beijing_1954 椭球体：Krasovsky_1940 长半轴：6378245.0 反扁率：298.3 中央经线：PRIMEM[“Greenwich”,0.0] - 格林威治单位：UNIT[“Degree”,0.0174532925199433]

（三）变换关系gtf文件 GEOGTRAN["CGS2000_TO_BJ1954", #转换关系名称 GEOGCS["GCS_China_Geodetic_Coordinate_System_2000", #数据原始坐标系名称 DATUM["D_China_2000", #原始坐标系椭球体名称及参数 SPHEROID["CGCS2000",6378137.0,298.257222101]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433] ], GEOGCS["GCS_Beijing_1954", #数据目标坐标系名称 DATUM["D_Beijing_1954", #目标坐标系椭球体名称及参数 SPHEROID["Krasovsky_1940",6378245.0,298.3]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433] ], METHOD["Geocentric_Translation"], #转换方法及参数 PARAMETER["X_Axis_Translation",0.0], PARAMETER["Y_Axis_Translation",0.0], PARAMETER["Z_Axis_Translation",0.0] ] 四、其他说明

1.以上是ArcGIS中投影文件和转换关系文件的含义、所在路径、内容结构和含义解释。 2.由于ArcMap所带的Arcpy为2.7版本，存在编码问题，同时很多三方库无法兼容，使用起来不是很友好，而Arcpro所带的Arcpy为python3.6版本，编码问题较少，同时能够更好的使用第三方库文件，Esri将大部分函数在老版本中的参数设定及命名直接搬到了新版本上，给用户带来很多便利，但仍存在部分细节性问题，目前网上资料较少，博主结合目前在做的工作，对投影转换部分进行梳理，欢迎大家一同交流进步。 3.文件参数解析部分参考了文章ArcGIS的.prj文件生成proj4格式的字符串，对博主的整理表示感谢。

【本文地址】

公司简介

联系我们