投影坐标与地理坐标的相互转化（Python）

        昨天遇到了个小问题，就是对方发来几个投影坐标，坐标系是大地2000，给了X和Y值，然后尝试用Arcgis做，但参照了几个教程之后，使用Arcgis Pro无论怎么转化得到的都是空值，鉴于之前的作业做过坐标转化，就截取一段代码，正好可以解决这个问题。这样不用安装Arcgis这类专业工具也可以完成相应任务。

1.1 数据格式

        首先将得到的坐标输入到表格中，生成csv文件，文件格式如图。第一列为X，第二列为Y。注意不要把X，Y顺序弄错。

1.2 运行代码

简单介绍下导入的两个库，pandas和pyproj。

第一个是读取数据需要的库，第二个则是用于坐标系换算的库，下面是官网地址：

Pandas官网：

pandas - Python Data Analysis Library (pydata.org)

pyproj官网：

pyproj 3.6.0 documentation (pyproj4.github.io)

在运行之前，需要进行环境配置，即安装这两个库。没安装的话可以pip install一下。

得到表格之后，我们在python中输入下列代码：

import pandas as pd
from pyproj import Transformer

df1 = pd.read_csv(r'E:\Desktop\postion2000.csv')
transformer = Transformer.from_crs(4550, 4490)
df1['lat'], df1['lon'] = transformer.transform(df1['x'],df1['y'])
df1.to_csv
print(df1)

在运行之前，需要确保库的安装成功和文件地址正确。此外，转换坐标系的WKID也应填写正确。

这里有个小插曲，对方只给了大地2000并没有给具体的投影坐标系，所以只能在中央经线最接近现实地理位置的坐标系中选择，大家根据自己的具体位置来查询应采用什么坐标系。

WKID即4550和4490，即坐标系的通用编号。查询方式有很多种，可以直接搜索也可以通过软件查询，笔者参考了Arcgis中的信息，具体查询方式为：

打开Arcgis，在工具箱中搜索坐标表转点，点击地球图标，就可以选择想要的坐标系进行查询：

点击详细信息，就可以查到对应的WKID。

可以看到，投影坐标系CGCS2000_3_Degree_GK_CM_123E的WKID为4550。同样的你也可以查询到你想要的坐标系信息。下面回到代码对应的位置改成你所需要的坐标系就可以直接运行。结果如下：