如何解決pythongdal投影坐標系轉換的問題

這篇文章主要介紹如何解決pythongdal投影坐標系轉換的問題，文中介紹的非常詳細，具有一定的參考價值，感興趣的小伙伴們一定要看完！

要將xian80地理坐標系轉換成投影坐標系：

xian1980 = """GEOGCS["GCS_Xian_1980",  DATUM["Xian_1980",    SPHEROID["Xian_1980",6378140.0,298.257]],  PRIMEM["Greenwich",0.0],  UNIT["Degree",0.0174532925199433]]"""prjreference = 'PROJCS["Xian_1980_3_Degree_GK_Zone_37",GEOGCS["GCS_Xian_1980",DATUM["D_Xian_1980",SPHEROID["Xian_1980",6378140.0,298.257]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Gauss_Kruger"],PARAMETER["False_Easting",37500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",111.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]'

轉換后會報錯：

latlong = transform.TransformPoint(x,y) File "/usr/lib/python3/dist-packages/osgeo/osr.py", line 711, in TransformPoint  return _osr.CoordinateTransformation_TransformPoint(self, *args)NotImplementedError: Wrong number or type of arguments for overloaded function 'CoordinateTransformation_TransformPoint'. Possible C/C++ prototypes are:  OSRCoordinateTransformationShadow::TransformPoint(double [3])  OSRCoordinateTransformationShadow::TransformPoint(double [3],double,double,double)**想想之前也沒發現過這個問題呀，最后發現之前投影坐標系：**'PROJCS["Xian_1980_3_Degree_GK_Zone_37",GEOGCS["GCS_Xian_1980",DATUM["D_Xian_1980",SPHEROID["Xian_1980",6378140.0,298.257]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",37500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",111.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]'**不是高斯-克呂格投影而是橫軸墨卡托，我想原因可能是gdal庫里只能支持高斯-克呂格投影吧。**

以上是“如何解決pythongdal投影坐標系轉換的問題”這篇文章的所有內容，感謝各位的閱讀！希望分享的內容對大家有幫助，更多相關知識，歡迎關注億速云行業資訊頻道！