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


Posted in Python onJanuary 17, 2020

要将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库里只能支持高斯-克吕格投影吧。**

以上这篇解决python gdal投影坐标系转换的问题就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
python字符串连接的N种方式总结
Sep 17 Python
分享给Python新手们的几道简单练习题
Sep 21 Python
python书籍信息爬虫实例
Mar 19 Python
Python 读取指定文件夹下的所有图像方法
Apr 27 Python
浅谈python 读excel数值为浮点型的问题
Dec 25 Python
Python图像处理之gif动态图的解析与合成操作详解
Dec 30 Python
使用python实现抓取腾讯视频所有电影的爬虫
Apr 15 Python
关于python中plt.hist参数的使用详解
Nov 28 Python
Python模块的定义,模块的导入,__name__用法实例分析
Jan 07 Python
python数据预处理方式 :数据降维
Feb 24 Python
python中sys模块是做什么用的
Aug 16 Python
Django项目配置Memcached和Redis, 缓存选择哪个更有优势
Apr 06 Python
python3 Scrapy爬虫框架ip代理配置的方法
Jan 17 #Python
Python计算机视觉里的IOU计算实例
Jan 17 #Python
python实现百度OCR图片识别过程解析
Jan 17 #Python
Python实现图像的垂直投影示例
Jan 17 #Python
基于python实现语音录入识别代码实例
Jan 17 #Python
Python操作MySQL数据库实例详解【安装、连接、增删改查等】
Jan 17 #Python
python实现高斯投影正反算方式
Jan 17 #Python
You might like
我的论坛源代码(八)
2006/10/09 PHP
深入理解PHP中的empty和isset函数
2016/05/26 PHP
php利用imagemagick实现复古老照片效果实例
2017/02/16 PHP
Prototype1.5 rc2版指南最后一篇之Position
2007/01/10 Javascript
jQuery对JSON数据进行排序输出的方法
2015/06/24 Javascript
基于javascript实现根据身份证号码识别性别和年龄
2016/01/22 Javascript
JavaScript中的时间处理小结
2016/02/24 Javascript
jQuery+CSS实现一个侧滑导航菜单代码
2016/05/09 Javascript
jQuery1.9+中删除了live以后的替代方法
2016/06/17 Javascript
使用JavaScript判断手机浏览器是横屏还是竖屏问题
2016/08/02 Javascript
js基础之DOM中元素对象的属性方法详解
2016/10/28 Javascript
ES6概念 Symbol toString()方法
2016/12/25 Javascript
jQuery中Datatables增加跳转到指定页功能
2017/02/08 Javascript
使用mint-ui实现省市区三级联动效果的示例代码
2018/02/09 Javascript
微信小程序实现banner图轮播效果
2020/06/28 Javascript
python单链表实现代码实例
2013/11/21 Python
python实现矩阵乘法的方法
2015/06/28 Python
Python中对象迭代与反迭代的技巧总结
2016/09/17 Python
解决python matplotlib imshow无法显示的问题
2018/05/24 Python
Django多数据库配置及逆向生成model教程
2020/03/28 Python
对Keras中predict()方法和predict_classes()方法的区别说明
2020/06/09 Python
基于PyTorch的permute和reshape/view的区别介绍
2020/06/18 Python
Python 中如何使用 virtualenv 管理虚拟环境
2021/01/21 Python
利用html5 file api读取本地文件示例(如图片、PDF等)
2018/03/07 HTML / CSS
HTML5 Convas APIs方法详解
2015/04/24 HTML / CSS
世界排名第一的运动鞋市场:Flight Club
2020/01/03 全球购物
最新会计专业求职信范文
2014/01/28 职场文书
聚美优品励志广告词
2014/03/14 职场文书
市场拓展计划书
2014/05/03 职场文书
优秀管理者事迹材料
2014/05/22 职场文书
小学三好学生事迹材料
2014/08/15 职场文书
心术观后感
2015/06/11 职场文书
心得体会该怎么写呢?
2019/06/27 职场文书
python库sklearn常用操作
2021/08/23 Python
详解JSON.parse和JSON.stringify用法
2022/02/18 Javascript
Python调用腾讯API实现人脸身份证比对功能
2022/04/04 Python