解决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进阶教程之文本文件的读取和写入
Aug 29 Python
Python中用pycurl监控http响应时间脚本分享
Feb 02 Python
详解Python Socket网络编程
Jan 05 Python
使用python编写简单的小程序编译成exe跑在win10上
Jan 15 Python
python画出三角形外接圆和内切圆的方法
Jan 25 Python
对python抓取需要登录网站数据的方法详解
May 21 Python
Python高斯消除矩阵
Jan 02 Python
python在回调函数中获取返回值的方法
Feb 22 Python
Python hashlib加密模块常用方法解析
Dec 18 Python
Python3爬虫带上cookie的实例代码
Jul 28 Python
Python发送邮件实现基础解析
Aug 14 Python
详解Python+Selenium+ChromeDriver的配置和问题解决
Jan 19 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
音乐朗读剧《MARS RED》2021年TV动画化决定!
2020/03/06 日漫
php 服务器调试 Zend Debugger 的安装教程
2009/09/25 PHP
PHP设计模式 注册表模式
2012/02/05 PHP
Thinkphp5框架实现图片、音频和视频文件的上传功能详解
2019/08/27 PHP
PHPstorm启用自动换行的方法详解(IDE)
2020/09/17 PHP
ExtJs Excel导出并下载IIS服务器端遇到的问题
2011/09/16 Javascript
转义字符(\)对JavaScript中JSON.parse的影响概述
2013/07/17 Javascript
火狐下input焦点无法重复获取问题的解决方法
2014/06/16 Javascript
javascript记录文本框内文字个数检测文字个数变化
2014/10/14 Javascript
JavaScript中的this关键字使用方法总结
2015/03/13 Javascript
js中substr,substring,indexOf,lastIndexOf,split,replace的用法详解
2015/11/09 Javascript
详谈for循环里面的break和continue语句
2017/07/20 Javascript
vue如何实现动态加载脚本
2020/02/05 Javascript
Vue +WebSocket + WaveSurferJS 实现H5聊天对话交互的实例
2020/11/18 Vue.js
python中的reduce内建函数使用方法指南
2014/08/31 Python
深入理解python中的闭包和装饰器
2016/06/12 Python
Python 提取dict转换为xml/json/table并输出的实现代码
2016/08/28 Python
TensorFlow实现Softmax回归模型
2018/03/09 Python
用python统计代码行的示例(包括空行和注释)
2018/07/24 Python
对python3新增的byte类型详解
2018/12/04 Python
利用python求积分的实例
2019/07/03 Python
200行python代码实现贪吃蛇游戏
2020/04/24 Python
python实现数据结构中双向循环链表操作的示例
2020/10/09 Python
python实现扫雷游戏的示例
2020/10/20 Python
美国知名日用品连锁超市:Dollar General(多来店)
2017/01/14 全球购物
现金出纳岗位职责
2014/03/15 职场文书
体育口号大全
2014/06/18 职场文书
党员干部群众路线个人整改措施
2014/09/18 职场文书
地道战观后感400字
2015/06/04 职场文书
费城故事观后感
2015/06/10 职场文书
关于感恩的歌曲整理(8首)
2019/08/14 职场文书
redis通过6379端口无法连接服务器(redis-server.exe闪退)
2021/05/08 Redis
mysql5.7使用binlog 恢复数据的方法
2021/06/03 MySQL
Win11 S Mode版本泄露 正式上线后叫做Windows 11 SE
2021/11/21 数码科技
如何利用Python实现n*n螺旋矩阵
2022/01/18 Python
vue修饰符.capture和.self的区别
2022/04/22 Vue.js