文章/答案/技术大牛

发布

如何为添加Cartopy shapely功能的GeoAxes添加图例？

为添加Cartopy shapely功能的GeoAxes添加图例，可以按照以下步骤进行操作：

首先，确保已经安装了Cartopy和shapely库，并导入所需的模块：

import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
from shapely.geometry import Point

创建一个GeoAxes对象，并设置投影方式：

fig = plt.figure()
ax = fig.add_subplot(1, 1, 1, projection=ccrs.PlateCarree())

添加地图背景和要素：

ax.add_feature(cfeature.LAND)
ax.add_feature(cfeature.OCEAN)
ax.add_feature(cfeature.COASTLINE)
ax.add_feature(cfeature.BORDERS)

创建一个shapely的Point对象，并设置其位置和属性：

point = Point(lon, lat)
attributes = {'name': 'City', 'population': 1000000}

将Point对象转换为GeoAxes坐标系下的点，并绘制在地图上：

ax.plot(point.x, point.y, 'ro', transform=ccrs.PlateCarree())

创建图例，并设置图例的标签和样式：

legend_labels = ['City']
legend_styles = ['ro']
ax.legend(legend_styles, legend_labels)

显示地图和图例：

plt.show()

这样，就可以为添加了Cartopy shapely功能的GeoAxes添加图例了。

关于腾讯云相关产品和产品介绍链接地址，由于要求不能提及具体品牌商，建议在腾讯云官方网站上查找与地理信息处理、地图可视化相关的云服务产品和解决方案。

相关·内容

Cartopy 系列教程：安装并绘制地图

Cartopy 也是一个 python 地图绘制包，同样能完成很多 Basemap 能实现的功能，而且旨在使数据分析及可视化尽可能简单。...其利用了强大的 PROJ.4，numpy，shapely库，提供了简单直观的绘图接口，可以创建满足出版质量的地图。...添加数据到地图一旦你创建好了地图，添加数据的方式就和常规的 matplotlib 绘图方法相同。默认情况下，添加到 GeoAxes 的任何数据的坐标系统和 GeoAxes 本身的坐标系统是相同的。...为了控制给定数据的坐标系，你可以添加 transform 关键词，并传递一个合适的 cartopy.crs.CRS 实例。...因为 cartopy 使用的 GeoAxes 类，这相当于询问地图的坐标轴范围。有时候会得到预期效果，但有时候并不会。

10.7K1 0

Python气象绘图教程（十三）—Cartopy_4

本节提要：关于子图的一些问题、使用path添加示意框线、Cartopy台风实例本土化一、关于子图的一些问题在某些时候，我们需要展示某个地区在整个地图中的位置，常规的方法是绘制两幅地图，比如一张为全国地图...唯有一点希望读者注意，在此时ax1与ax2已经不是一类子图了，因为ax2在使用了projection命令之后，已经转变为cartopy中的GeoAxes。...但是当我们在GeoAxes中使用这个命令时，他是无效的。（然而并不会出现报错的情况），这是因为在GeoAxes中，设置了新的命令来管理框线——ax.outline_patch。...三、Cartopy台风实例本土化在官网上，有一个台风行径影响图，我在学习时将其中国化并实用化，添加了一些中文注释，并增加了一个本土化例子，希望能帮助到学习的同学： import numpy as np...+10, 10)) gl.ylocator = mticker.FixedLocator(np.arange(extent[2], extent[3]+10, 10)) #######添加图例

9.7K8 2

【ProPlot库】ProPlot3兰伯特投影-可添加刻度(三)

提供的 cartopy.mpl.geoaxes.GeoAxes.gridlines 方法添加 grid： import proplot as plot proj = plot.Proj('lcc',...添加自定义刻度虽然 cartopy 提供了 gridlines 的方法，但经度的标签出现在了 y 轴上（80°E），并且和matplotlib 的 matplotlib.axes.Axes.set_xticks...网上的解决方案并不多，自己尝试了很多办法。下面提供一种可以给 cartopy 兰伯特投影添加 major、minor ticks 的方法(略微有些复杂)。...在proplot中,目前我只找到了用ax.tick_params添加geoaxes 的tickmarker 的例子，在format 里给地图设置刻度标签的方法可能还未实现。...，为geoAxes 添加一个笛卡尔坐标系的轴，再使用 lambert_xticks 函数添加标签。

2K2 1

Python气象绘图教程（十五）—Cartopy_5

最后，通过前面几节提到的添加图例的方法，完善图形。...最后添加图例时候，我用的是逐个添加，实际上可以用循环的方式添加。...其中，mpatches.Rectangle这个语句就是添加图例上的色条，rectangle就是矩形的意思，延伸的，你也可以使用其他形状添加图例（包括自定义形状）。...首先说crs，这个是GeoAxes的基础，只有在projection设置了投影之后，才能添加地图。...在cartopy=0.17中，不能使用draw_labels=True来为除PlateCarree、Mercator之外的投影添加经纬标签，不过据说在0.18版本中已经优化，读者可以试试。

11.5K5 2

这份Cartopy绘图入门指南，请查收

，pyecharts在地学可视化中功能实在过于简陋;geopandas是基于pandas的，一般用于商业数据分析，所以我更推荐大家学习Cartopy，虽然Cartopy气象专业用得很多哈哈哈 Cartopy...GeoAxes 类为特定于绘制地图的轴添加了额外的功能。创建一个 GeoAxes 对象的办法是，在创建 axes（或 subplot）时，通过参数 projection 指定一个 ccrs 中的投影。...因此用 Cartopy 画地图的基本流程并不复杂：创建画布。通过指定 projection 参数，创建 GeoAxes 对象。调用 GeoAxes 的方法画图。...，仅有海岸线，尝试一下添加更多地理信息添加预定义要素首先需要导入一个cartopy.feature 常量，为了简化一些非常常见的情况，如大陆海洋国界等cartopy都已经进行了预定义，但需要注意的是直接导入的中国国界线并不是标准的.../generated/cartopy.mpl.geoaxes.GeoAxes.html#cartopy.mpl.geoaxes.GeoAxes.add_feature cartopy更多地图 https

5.3K2 1

想知道世界大坝分布情况吗？这篇可视化作品告诉你

图2 「象形柱状图部分」原作品中下方部分的象形柱状图也是非常的形象生动切合主题，通过观察可以发现每个大坝logo代表数值200，而末尾不足200的部分就按照余数/200的「透明度」进行渲染，配合右下角的图例...2.2 半球部分的制作「正射投影」部分，我选择使用cartopy.ccrs内置的Orthographic()，通过传入中央经纬度，即可得到期望的半球面： ?...图3 我们利用mpl_toolkits.axes_grid1.inset_locator中的inset_axes()将两个半球各自对应的axe对象插入到主体axe中，再利用cartopy的add_geometries...0.58, 0.26, 0.26), bbox_transform=ax.transAxes, axes_class=cartopy.mpl.geoaxes.GeoAxes...图6 而右下角图例中的第一个logo上下渐变的效果其实是配合numpy数组，从上往下线性地降低rgba第四个通道的透明度值得到了，非常的容易~ ?

6694 0

（在模仿中精进数据可视化07）星球研究所大坝分布可视化

图2 象形柱状图部分　　原作品中下方部分的象形柱状图也是非常的形象生动切合主题，通过观察可以发现每个大坝logo代表数值200，而末尾不足200的部分就按照余数/200的透明度进行渲染，配合右下角的图例...摸清楚要做的内容之后，下面让我们开始吧！ 2.2 半球部分的制作正射投影部分，我选择使用cartopy.ccrs内置的Orthographic()，通过传入中央经纬度，即可得到期望的半球面： ?...图3 　　我们利用mpl_toolkits.axes_grid1.inset_locator中的inset_axes()将两个半球各自对应的axe对象插入到主体axe中，再利用cartopy的add_geometries...0.58, 0.26, 0.26), bbox_transform=ax.transAxes, axes_class=cartopy.mpl.geoaxes.GeoAxes...图6 　　而右下角图例中的第一个logo上下渐变的效果其实是配合numpy数组，从上往下线性地降低rgba第四个通道的透明度值得到了，非常的容易~ ?

6212 0

在模仿中精进数据可视化07：星球研究所大坝分布可视化

图2 「象形柱状图部分」原作品中下方部分的象形柱状图也是非常的形象生动切合主题，通过观察可以发现每个大坝logo代表数值200，而末尾不足200的部分就按照余数/200的「透明度」进行渲染，配合右下角的图例...2.2 半球部分的制作「正射投影」部分，我选择使用cartopy.ccrs内置的Orthographic()，通过传入中央经纬度，即可得到期望的半球面：图3 我们利用mpl_toolkits.axes_grid1....inset_locator中的inset_axes()将两个半球各自对应的axe对象插入到主体axe中，再利用cartopy的add_geometries进行矢量元素的叠加和色彩映射即可，我已经在数据中算好了归一化数值方便色彩映射...0.58, 0.26, 0.26), bbox_transform=ax.transAxes, axes_class=cartopy.mpl.geoaxes.GeoAxes...，即可在前面的基础上，得到下面的图像：图6 而右下角图例中的第一个logo上下渐变的效果其实是配合numpy数组，从上往下线性地降低rgba第四个通道的透明度值得到了，非常的容易~ 图7 而图中其他的小元素譬如各种文字就不赘述了

6931 0

Python可视化 | 三维地图可视化实例

首先必须明确一点，matplotlib的axes3D这个投影中，是不能用add_geometry这个功能来直接将读取到的shp文件添加上去的。...add_geometry这个功能是cartopy下的geoaxes才能使用，同理add_feature也不能再3d图中使用。但是这个功能确实又是比较常用，而且酷炫的。...所以我在s站上查到了cartopy库包的开发人员直接给出的回答。我们不直接开讲怎么绘制，而先回忆在不久之前的推文中，我们使用过的一个功能。...而是使用了shapely的一个功能，构造了一个geometry import shapely.geometry as sgeom x2=[105,115,115,105,105] y2=[25,25,35,35,25...恰巧，matplotlib的axes3D投影中，允许我们使用polygon功能。那么我们是不是可以将shp文件中的geometry读取出来，转变成polygon，然后添加到三维图中呢。

4.3K2 1

GPM卫星 IMERG 降水产品的简单可视化

xarray as xr #可视化 import matplotlib.pyplot as plt import cartopy.crs as ccrs import cartopy.feature..., LATITUDE_FORMATTER import cmaps from shapely import geometry as sgeom # 辅助模块 from glob import glob...1.2 IMERG降水可视化加入修正cartopy无法正常显示经纬度的函数 ''' 参考Andrew Dawson 提供的解决方法: https://gist.github.com/ajdawson/...def find_side(ls, side): """ Given a shapely LineString which is assumed to be rectangular, return...Use GeoAxes.spines['geo'] or the default Axes properties instead. 1.3 GMI(hdf5) # 读取GPM GMI降水产品 import

1602 0

使用Matplotlib & Cartopy绘制我国台风路径图

主要第三方库用到的主要工具包有pandas、numpy、matplotlib、cartopy、shapely，前三个库大家可能都熟悉，下面介绍下后两个库的使用场景。...❝cartopy：基于matplotlib的python地理数据处理和可视化库，本文会用来展示地图 shapely: 是一个对几何对象进行操作和分析的Python库，本文用来处理点线数据 ❞ cartopy...文档：https://scitools.org.uk/cartopy/docs/latest/ shapely文档：https://shapely.readthedocs.io/en/stable/...方式有很多种，既可以用离线的GeoJson数据，也可以用JPG图片，或者第三方库提供的地图。我这里用的是cartopy内置的地图数据，可以很方便的修改配置属性。...首先导入本次会用到的所有库： # cartopy：用来获取地图 import cartopy.crs as ccrs import cartopy.feature as cfeature # matplotlib

3.1K2 0

地图1-图例式地图

图例式地图图例式地图往往是论文中的第一张图片，用来展示研究区域的某种特征。一般有讲述研究地区的地貌图、展示数据站点分布的站点图、展示某些特定图注信息的地图，或者这些地图的叠加形式。...图1即是典型的图例式地图，主要介绍研究区域的情况，没有涉及到太多具体的研究数据和结论。图2则是一张综合式的图例式地图，主图为地貌图，附图是介绍地质结构和海拔的图像。...引入这种理论的原因是matplotlib自带的legend命令绘制图例的时候不是非常方便，那么我们需要使用其他的技法来实现图例的绘制。...as ccrs import cartopy.mpl.ticker as cmt import cartopy.mpl.gridliner as cmg import cartopy.mpl.patch...在添加图例的过程中，我们并没有直接使用legend命令，而是新建了一个小ax作为载具来承载图例。ax的绘图方法是比legend丰富的。

680 0

基于geopandas的精美地图绘制：geoplot

温馨提示由于可视化代码过长隐藏，可点击运行Fork查看若没有成功加载可视化图，点击运行可以查看 ps:隐藏代码在【代码已被隐藏】所在行，点击所在行，可以看到该行的最右角，会出现个三角形，点击查看即可...前言项目目标本项目旨在介绍geoplot的地图绘图示例。...lib/python3.9/site-packages (from contextily>=1.0.0->geoplot) (1.3.6) Requirement already satisfied: shapely...(2022.1) Requirement already satisfied: pyshp>=2.3 in /opt/conda/lib/python3.9/site-packages (from cartopy...，当然还有更多功能有待发掘

1141 0

python计算与绘制WRF降水量

as ccrs import cartopy.feature as cfeature import cartopy.io.shapereader as shpreader from cartopy.mpl.ticker...import shapely.geometry as sgeom import cmaps 简单计算与绘图(累计值) from netCDF4 import Dataset from wrf import...cmap=cmaps.radar, transform=ccrs.PlateCarree()) # 为降雨量添加...shrink=1) cbar.set_label('Rainfall (mm)', fontdict={'size':20}) cbar.ax.tick_params(labelsize=20) # 添加经纬度网格线...ax.gridlines(color='black', linestyle='dotted') # 设置标签大小 plt.tick_params(labelsize=15) # 添加标题 plt.title

6351 1

气象绘图——白化杂谈

库包下的clip白化； ④公众号DataCharm推送的Fiona、shapely库包赋nan值白化； ⑤气象家园Masterpiece提出的shp转path判别法。...在上一期提到的三维地图的添加时我们就提到了各种绘图集合图形之间是可以相互转化的，我们还是得使用那个geos_to_path的命令，将shp里的集合形状geometry转化成path，然后用set_clip_path...不过大佬为了丰富白化功能，增添了一个region参数，使得我们能够凭关键词筛选要白化的省份，并且不使用cartopy中的geos_to_path命令，而是使用的shapefile库包来读取边界。...如果看过气象家园的文章，可以知道作者提供了一个名为maskout.py的Python程序。但是不知道怎么引入。这里我提供了一个懒办法——添加指定的读取路径。...fiona和shapely的方法可以原汁原味的返回原来网格。

1.3K3 2

GPM 卫星三级产品的简单可视化

as ccrs import cartopy.feature as cfeat import cartopy.io.shapereader as shpreader from cartopy.mpl.gridliner...import LONGITUDE_FORMATTER, LATITUDE_FORMATTER import cmaps from shapely import geometry as sgeom #...1.2 IMERG降水可视化加入修正cartopy无法正常显示经纬度的函数 ''' 参考Andrew Dawson 提供的解决方法: https://gist.github.com/ajdawson/...dd536f786741e987ae4e ''' # 给出一个假定为矩形的LineString，返回对应于矩形给定边的直线。...def find_side(ls, side): """ Given a shapely LineString which is assumed to be rectangular, return

1531 0

库|分享一些python库

GeoPandas 这不算是一个单纯的地理数据分析库。这个库类似pandas，但是为了能够更容易的进行地理数据处理而设计。GeoPandas扩展了pandas的数据类型，从而允许对几何类型的空间操作。...GeoPandas利用了shapely库来执行几何操作。说到地理数据操作，也应该提及xarray库，尤其是对于地球科学中的高维数据表示来说，这个库是非常适宜了。...xarray也对pandas有很好的支持。之前也专门写过xarray库的介绍。...和cartopy。...基于HoloViews构建而成，而且基于Cartopy库GeoViews添加了一部分地理图形，绘图部分实用matplotlib或Bokeh完成。能够很好的支持geopandas数据结构。

1K1 0

一刻也没有为圣诞的结束而悲伤，下一刻赶来的是ERA5数据计算700hPa水汽通量散度

绘制水汽通量散度分布图将计算结果可视化，展示水汽通量散度的空间分布特征。逐步分析水汽通量散度的单位对每一个计算步骤的物理单位进行剖析，确保物理意义清晰且计算结果准确。项目功能 1....绘制水汽通量散度图将计算结果可视化，使用 cartopy 库绘制水汽通量散度的空间分布图。...代码示例如下： import matplotlib.pyplot as plt import cartopy.crs as ccrs import cartopy.feature as cfeature...ax.pcolormesh(div_q.longitude, div_q.latitude, div_q * 1e6, cmap='RdBu', vmin=-1, vmax=1) # 添加色标.../mpl/geoaxes.py:1797: MatplotlibDeprecationWarning: shading='flat' when X and Y have the same dimensions

1851 0

WRF domain 绘制改进

as mtransforms import cartopy.crs as ccrs from cartopy.io.shapereader import BasicReader import cmaps...from matplotlib.path import Path from cartopy.mpl.patch import geos_to_path from mpl_toolkits.axes_grid1....inset_locator import TransformedBbox, BboxPatch, BboxConnector import shapely.geometry as sgeom from.../ajdawson/dd536f786741e987ae4e ''' # 给出一个假定为矩形的LineString，返回对应于矩形给定边的直线。...def find_side(ls, side): """ Given a shapely LineString which is assumed to be rectangular, return

1591 0

（数据科学学习手札83）基于geopandas的空间数据分析——geoplot篇(下)

2 geoplot进阶　　上一篇文章中的pointplot()、polyplot以及webmap()帮助我们解决了在绘制散点、基础面以及添加在线地图底图的问题，为了制作出信息量更丰富的可视化作品，我们需要更强的操纵矢量数据与映射值的能力...：list型，用于自定义图例显示的各个具体数值对应的文字标签，与legend_values搭配使用 legend_kwargs：字典，在legend参数设置为True时来传入更多微调图例属性的参数...legend_labels：list型，用于自定义图例显示的各个具体数值对应的文字标签，与legend_values搭配使用 legend_kwargs：字典，在legend参数设置为True时来传入更多微调图例属性的参数...子模块来源于cartopy，这一点我跟geoplot的主要开发者聊过，他表示geoplot暂时不支持geopandas中那样自定义任意投影或使用EPSG投影，而是内置了一系列常用的投影，譬如我们上文中绘制美国区域时频繁使用到的...AlbersEqualArea()即之前我们在geopandas中通过proj4自定义的阿尔伯斯等面积投影，其他常见投影譬如Web Mercator、Robinson，或者直接绘制球体地图，如本文开头的图

1.8K3 0

点击加载更多

扫码

添加站长进交流群

领取专属 10元无门槛券

手把手带您无忧上云