GDAL命令-坐标快速转换
cs2cs 是PROJ的命令行工具，如果需要解决很多坐标转换的内容，建议使用这个方式进行快速处理。
这个命令使用有一个前提：配置好PROJ4的配置文件，通常情况下载Python中安装好proj文件后就可以进行使用。找到proj对应的nad文件，执行如下命令，将nad文件放到计算机环境中，这个文件主要是坐标系文件存放位置。
set PROJ_LIB=C:\PROJ\NAD
假设您在名为Coordinates.txt的文件中的EPSG：4326（WGS84）中拥有这些疯狂的准确超坐标（例如，从Excel中复制列）。
108.8196257 36.8196257
108.4590673 36.45906729
108.6955076 36.69550756
108.2114871 36.21148709
108.3992153 36.39921532
108.78723 36.78722996
108.1194116 36.1194116
108.1644091 36.16440907
108.2315863 36.23158627
108.0148562 36.01485621
并且您希望它们位于非常优越的公制EPSG：4545中。您只需使用+init=sourceCRS +to +init=targetCRS。如果您不知道坐标所在的CRS，请在http://projfinder.com/查看。
cs2cs + init = epsg：4326 + to + init = epsg：4545 test.txt它会打印
573128.93 4076807.25 0.00
541150.11 4036580.24 0.00
562154.72 4062945.42 0.00
519017.43 4009030.22 0.00
535812.52 4029914.72 0.00
570268.05 4073187.77 0.00
510750.34 3998799.20 0.00
514792.91 4003798.11 0.00
520819.47 4011264.62 0.00
501339.24 3987191.20 0.00
不过，您可能想要更多的小数位，为此，您可以使用-f指定输出格式。
cs2cs + init = epsg：4326 + to + init = epsg：4545 -f “％.17f” test.txt
573128.93235548702068627 4076807.24640083173289895 0.00000000000000000
541150.11231174739077687 4036580.23916309839114547 0.00000000000000000
562154.71878067729994655 4062945.42327193915843964 0.00000000000000000
519017.42718655971111730 4009030.21620059525594115 0.00000000000000000
535812.52077009994536638 4029914.71937049925327301 0.00000000000000000
570268.05317815206944942 4073187.77400729805231094 0.00000000000000000
510750.33795329823624343 3998799.19788714684545994 0.00000000000000000
514792.90707002923591062 4003798.11148226726800203 0.00000000000000000
520819.47146004205569625 4011264.61779326666146517 0.00000000000000000
501339.23850690264953300 3987191.20464842254295945 0.00000000000000000
同时可以通过添加输出重定向来将这些转换后的坐标定向到一个名为MyCoolCoordinates.txt的新文件中：
cs2cs + init = epsg：4326 + to + init = epsg：4545 Coordinates.txt > CoolCoordinates.txt
