成都网站建设设计

将想法与焦点和您一起共享

如何解决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投影坐标系转换的问题”这篇文章的所有内容,感谢各位的阅读!希望分享的内容对大家有帮助,更多相关知识,欢迎关注创新互联行业资讯频道!


新闻标题:如何解决pythongdal投影坐标系转换的问题
URL网址:http://chengdu.cdxwcx.cn/article/pejgoc.html