解決python gdal投影座標系轉換的問題
阿新 • • 發佈:2020-01-18
要將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",0.0174532925199433]],PROJECTION["Gauss_Kruger"],PARAMETER["False_Easting",37500000.0],PARAMETER["False_Northing",PARAMETER["Central_Meridian",111.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",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) **想想之前也沒發現過這個問題呀,最後發現之前投影座標系:** 'PROJCS["Xian_1980_3_Degree_GK_Zone_37",PROJECTION["Transverse_Mercator"],1.0]]' **不是高斯-克呂格投影而是橫軸墨卡託,我想原因可能是gdal庫裡只能支援高斯-克呂格投影吧。**
以上這篇解決python gdal投影座標系轉換的問題就是小編分享給大家的全部內容了,希望能給大家一個參考,也希望大家多多支援我們。