python用cartopy包画地图_python – 使用Cartopy在地图上显示图像时的投影问题
我有一些我想用Cartopy显示的卫星图像数据.我已成功按照图像示例详细介绍了here.导致此代码:
import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
fig = plt.figure(figsize=(12,12))
img_extent = (-77,-59,9,26)
ax = plt.axes(projection=ccrs.PlateCarree())
# image data coming from server,code not shown
ax.imshow(img,origin='upper',extent=img_extent)
ax.set_xmargin(0.05)
ax.set_ymargin(0.10)
# mark a known place to help us geo-locate ourselves
ax.plot(-117.1625,32.715,'bo',markersize=7)
ax.text(-117,33,'San Diego')
ax.coastlines()
ax.gridlines()
plt.show()
此代码生成以下图像
我的问题是卫星图像数据不在PlateCarree投影中,而是墨卡托投影.
但是当我得到轴对象时
ax = plt.axes(projection=ccrs.Mercator())
我失去了海岸线.
我看到了here报告的问题.但是
ax.set_global()
得到这张图片的结果:
数据不存在,圣地亚哥位置错误.纬度/经度范围也发生了变化.我究竟做错了什么?
发布讨论更新
主要问题是我没有使用transform_points方法在目标投影中正确指定图像范围.我还必须具体说明菲尔建议的imshow方法中的坐标参考系统.这是正确的代码:
import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
proj = ccrs.Mercator()
fig = plt.figure(figsize=(12,12))
extents = proj.transform_points(ccrs.Geodetic(),np.array([-77,-59]),np.array([9,26]))
img_extents = (extents[0][0],extents[1][0],extents[0][6],extents[1][7] )
ax = plt.axes(projection=proj)
# image data coming from server,extent=img_extents,transform=proj)
ax.set_xmargin(0.05)
ax.set_ymargin(0.10)
# mark a known place to help us geo-locate ourselves
ax.plot(-117.1625,markersize=7,transform=ccrs.Geodetic())
ax.text(-117,'San Diego',transform=ccrs.Geodetic())
ax.coastlines()
ax.gridlines()
plt.show()
导致这个正确地理投影的卫星图像:
python用cartopy包画地图_python – 使用Cartopy在地图上显示图像时的投影问题相关推荐
- python用cartopy包画地图_python绘制地图的利器Cartopy使用说明
python绘制地图一般使用Basemap绘图包,但该包配置相对较繁琐,自定义性不强,这里介绍一个绘制地图的利器Cartopy,个人认为该工具方便.快捷,附上一些自己写的程序. 准备工作,工欲善其事, ...
- python的turtle怎么画曲线_python怎么画曲线图,
如何用python turtle 画n阶希尔伯特曲线 工大少年你好 Python如何画函数的曲线 输下导入我们用到数库. >>> import numpy as np >> ...
- python 引入同一路径的类_python小课堂15 - 史上最详细的包和模块import讲解篇
python小课堂15 - 史上最详细的包和模块import讲解篇 前言 在大量的代码设计中,我们不可能将所有代码都写在一个.py文件,所以有了包.模块,而为了代码可以重复利用(复用性),就有了类.函 ...
- python函数实例与类实例_python – 在类和实例上表现不同的函数
好问题!您使用描述符可以轻松完成您所寻求的内容. Descriptors是实现描述符协议的Python对象,通常以__get __()开头. 它们主要存在于不同类中的类属性.在访问它们时,将调用它们的 ...
- python散点地图_Python交散点地图
原标题:Python交散点地图 leftlet给R语言提供了很好用的交互式动态地图接口,其在Python中得API接口包名为folium(不知道包作者为何这样起名字,和leaflet已经扯不上关系了) ...
- python circle函数如何画圆_Python练习实例56 | 画图,学用circle画圆形
作者还记得在抖音上,有不少的Python学习爱好者还能够用Python里面自带的工具,来画各种各样的图形. 一句话,秀的简直不行不行的. 小黄人.jpg 那么今天,我们也来开始用Python来画出第一 ...
- python用turtle库画正方形_Python也可以这么“文雅”,用turtle库画一幅高傲的梅花...
文末有惊喜!!! 前言: "梅花"我们最早认识它是在小学的课本诗词中,在诗人心中是那么的文雅.坚韧,梅花是花中四君子之一,傲而不俗.好梅花就说到这里,下面讲怎样用Python画一幅 ...
- python恶搞表情包下载大全_Python:恶搞,将你朋友照片做成熊猫人表情包
本文的文字及图片来源于网络,仅供学习.交流使用,不具有任何商业用途,版权归原作者所有,如有问题请及时联系我们以作处理 以下文章来源于腾讯云 作者:龙哥 ( 想要学习Python?Python学习交流群 ...
- python读取csv数据画直方图_Python数据可视化(Pygal、svg、csv、Bar、Line)
一.pygal(图表类型Bar) 将使用Python可视化包Pygal来生成可缩放的矢量图形文件 pygal官方文档:[www.pygal.org/en/stable/](http://www.pyg ...
最新文章
- 扩增子统计绘图6韦恩图:比较组间共有和特有OTU或分类单元
- 二三代基因组混合组装流程的搭建与香菇基因组精细图谱的获得
- Teaching Machines to Converse--阅读笔记1 摘要及介绍
- java 首字母小写_java实现将字符串中首字母转换成大写,其它全部转换成小写的方法示例...
- gan简介_GAN简介
- BZOJ 2440 完全平方数(莫比乌斯-容斥原理)
- 企业管理软件随想-透视-包容,无形思想-有形方便
- php-fpm进程利用CPU不均问题的优化过程
- ES6 数组高频使用方法
- python语句示例_Python学习笔记之if语句的使用示例
- 疑难杂症篇(一)--安装Visio与已安装的office冲突的解决方案
- 浅谈逆向——从案例谈OD的使用(OD的使用2)
- catia圆管焊接焊接_CATIA焊接设计实例教程
- 手把手教你制作炫酷的PCB板3D效果图
- 以计算机为主题的英语作文范例,英语作文范例,一篇万能英语作文范文。
- OpenStack虚拟云桌面在携程呼叫中心的应用
- jarvis OJ web babyphp
- java端微信公众号发送模板信息乱码解决
- 君子喻于义,小人喻于利
- 软件系统六大设计原则
热门文章
- java语言中如何表示素数,使用Java语言求素数的几个方法
- HyperLynx(十六)PCI-E的设计与仿真
- 自动化测试框架[Cypress命令行执行测试详解]
- 安兔兔苹果html5排行榜,iPhone8Plus最强?9月安兔兔手机性能排行榜出炉
- icloud连接服务器时出现问题_iphone8 icloud连接服务器时出现问题
- gRPC大数据量消息传递方法
- 「链节点活动年度总结」2019年区块链行业会议回顾
- Mac OS X中virtualenv里python shell无法使用光标键问题的解决
- C语言求1到100的和
- win10修改ntp服务器地址,修改win10 ntp服务器地址