解决python gdal投影坐标系转换的问题

解决python gdal投影坐标系转换的问题,第1张

概述解决python gdal投影坐标系转换的问题 要将xian80地理坐标系转换成投影坐标系: xian1980 = """ GEOGCS["GCS_Xian_1980", DATUM["Xian_1980", SPHEROID["Xian_1980",6378140.0,298.257]], PRIMEM["Greenwich",0.0], UNIT["Degree",0.0174532925199433]]&qu ...

要将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投影坐标系转换的问题就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持我们。

总结

以上是内存溢出为你收集整理的解决python gdal投影坐标系转换的问题全部内容,希望文章能够帮你解决解决python gdal投影坐标系转换的问题所遇到的程序开发问题。

如果觉得内存溢出网站内容还不错,欢迎将内存溢出网站推荐给程序员好友。

欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/langs/1199391.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2022-06-04
下一篇 2022-06-04

发表评论

登录后才能评论

评论列表(0条)

保存