脚本之家 服务器常用软件
快捷导航
您的位置：首页 → 脚本专栏 → python → python Cartopy

python Cartopy的基础使用详解

 更新时间：2020年11月01日 10:02:20   作者：小龙人_tao  
这篇文章主要介绍了python Cartopy的基础使用详解，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧

前言

常用地图底图的绘制一般由Basemap或者cartopy模块完成，由于Basemap库是基于python2开发的一个模块，目前已经不开发维护。故简单介绍cartopy模块的一些基础操作。 一、基础介绍

首先导入相关模块。

import numpy as np
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature
from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter

首先介绍参数projection，该命令可以配合ccrs设置投影类型，此处以方形投影命令为示例。其中central_longitude参数为投影中心位置。其中心设置与Basemap设置规则一样，详情可以看上一篇文章。

ax=plt.axes(projection=ccrs.PlateCarree(central_longitude=0))

在设置好绘制类型后，绘制地图各特征量。其代码如下:

#ax.add_feature(cfeature.LAKES.with_scale(scale))
ax.add_feature(cfeature.OCEAN.with_scale(scale))
#ax.add_feature(cfeature.RIVERS.with_scale(scale))
#ax.add_feature(cfeature.LAND.with_scale(scale),lw=0.5)
ax.add_feature(cfeature.COASTLINE.with_scale(scale),lw=2)

参数scale为地图分辨率，目前支持10m,50m,110m,参数lw为线条粗细。此处绘制海岸线和海洋，效果图如下:

在这里插入图片描述

在绘制结束后，作为地图。经纬度自然是必不可少的，在该模块中，引进同时设置坐标轴标签改变该标签刻度的表示，具体形式如下:

ax.set_xticks(np.arange(0,361,40), crs=ccrs.PlateCarree())
ax.set_yticks(np.arange(-90,90+30,30), crs=ccrs.PlateCarree())
#zero_direction_label用来设置经度的0度加不加E和W
lon_formatter = LongitudeFormatter(zero_direction_label=False)
lat_formatter = LatitudeFormatter()
ax.xaxis.set_major_formatter(lon_formatter)
ax.yaxis.set_major_formatter(lat_formatter)

可以看到效果图如下:

在这里插入图片描述

当然如果想对坐标轴粗细变化可以引入一下命令。

ax.outline_patch.set_visible(False)
ax.spines['bottom'].set_visible(True)
ax.spines['left'].set_visible(True)
ax.spines['right'].set_visible(True)
ax.spines['top'].set_visible(True)
ax.spines['bottom'].set_linewidth(2.5);###设置底部坐标轴的粗细
ax.spines['left'].set_linewidth(2.5);####设置左边坐标轴的粗细
ax.spines['right'].set_linewidth(2.5);###设置右边坐标轴的粗细
ax.spines['top'].set_linewidth(2.5);####设置上部坐标轴的粗细

应该在该模块下，控制坐标轴的命令已经和常规不一样。因此先关闭该控制，然后开启常规坐标轴设置。

二、区域地图的绘制

当我们在某一小块区域研究时，需要绘制区域地图。此时我们可以引入命令：

ax.set_extent(box,crs=ccrs.PlateCarree())

其中box为绘制区域，crs为投影类型。其他命令基本不变。设置box为[40,180,0,90],可得到效果图如下:

在这里插入图片描述

总结

为方便各位读者，我书写了绘制地图的函数，大家在使用时可直接调用。此处示例为方形投影，若希望绘制其他投影。只需要修改函数部分参数即可。代码如下:

def map_make(scale,box,xstep,ystep):
  ax=plt.axes(projection=ccrs.PlateCarree(central_longitude=180))
  a = (box[1]-box[0])//xstep
  x_start = box[1] - a*xstep
  a = (box[3]-box[2])//ystep
  y_start = box[3] - a*ystep
  ax.set_extent(box,crs=ccrs.PlateCarree())
  #ax.add_feature(cfeature.LAKES.with_scale(scale))
  #ax.add_feature(cfeature.OCEAN.with_scale(scale))
  #ax.add_feature(cfeature.RIVERS.with_scale(scale))
  #ax.add_feature(cfeature.LAND.with_scale(scale),lw=0.5)
  ax.add_feature(cfeature.COASTLINE.with_scale(scale),lw=2)
  
  ax.set_xticks(np.arange(x_start,box[1]+xstep,xstep), crs=ccrs.PlateCarree())
  ax.set_yticks(np.arange(y_start,box[3]+ystep,ystep), crs=ccrs.PlateCarree())
  #zero_direction_label用来设置经度的0度加不加E和W
  lon_formatter = LongitudeFormatter(zero_direction_label=False)
  lat_formatter = LatitudeFormatter()
  ax.xaxis.set_major_formatter(lon_formatter)
  ax.yaxis.set_major_formatter(lat_formatter)
  #添加网格线
  ax.grid()
  
  ax.outline_patch.set_visible(False)
  ax.spines['bottom'].set_visible(True)
  ax.spines['left'].set_visible(True)
  ax.spines['right'].set_visible(True)
  ax.spines['top'].set_visible(True)
  ax.spines['bottom'].set_linewidth(2.5);###设置底部坐标轴的粗细
  ax.spines['left'].set_linewidth(2.5);####设置左边坐标轴的粗细
  ax.spines['right'].set_linewidth(2.5);###设置右边坐标轴的粗细
  ax.spines['top'].set_linewidth(2.5);####设置上部坐标轴的粗细
  
  return ax

到此这篇关于python Cartopy的基础使用详解的文章就介绍到这了,更多相关python Cartopy内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！

您可能感兴趣的文章:

相关文章

  • 一文教会你用Python实现pdf转word

    一文教会你用Python实现pdf转word

    python实现pdf转word,支持中英文转换,转换精度高,可以达到使用效果,下面这篇文章主要给大家介绍了关于用Python实现pdf转word的相关资料,需要的朋友可以参考下
    2023-01-01
  • pytorch加载自己的数据集源码分享

    pytorch加载自己的数据集源码分享

    这篇文章主要介绍了pytorch加载自己的数据集源码分享,标准的数据集流程梳理分为数据准备以及加载数据库–>数据加载器的调用或者设计–>批量调用进行训练或者其他作用，需要的朋友可以参考下
    2022-08-08
  • Django 权限认证(根据不同的用户，设置不同的显示和访问权限)

    Django 权限认证(根据不同的用户，设置不同的显示和访问权限)

    这篇文章主要介绍了Django 权限认证(根据不同的用户，设置不同的显示和访问权限)，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
    2019-07-07
  • Python使用扩展库pywin32实现批量文档打印实例

    Python使用扩展库pywin32实现批量文档打印实例

    这篇文章主要介绍了Python使用扩展库pywin32实现批量文档打印实例，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
    2020-04-04
  • Python封装解构以及丢弃变量

    Python封装解构以及丢弃变量

    这篇文章主要介绍了Python封装解构以及丢弃变量，文章围绕主题展开详细的内容介绍，具有一定的参考价值，需要的朋友可以参考一下
    2022-09-09
  • 浅析python连接数据库的重要事项

    浅析python连接数据库的重要事项

    这篇文章主要介绍了python连接数据库的重要事项,本文给大家介绍的非常详细，对大家的学习或工作具有一定的参考借鉴价值，需要的朋友可以参考下
    2021-02-02
  • opencv实现图像几何变换

    opencv实现图像几何变换

    这篇文章主要为大家详细介绍了opencv实现图像几何变换，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
    2021-03-03
  • Python顺序结构语句详解

    Python顺序结构语句详解

    这篇文章主要为大家介绍了Python顺序结构语句，具有一定的参考价值，感兴趣的小伙伴们可以参考一下，希望能够给你带来帮助
    2022-01-01
  • python的常见矩阵运算(小结)

    python的常见矩阵运算(小结)

    这篇文章主要介绍了python的常见矩阵运算(小结)，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
    2019-08-08
  • Python tkinter padx参数详解

    Python tkinter padx参数详解

    这篇文章主要介绍了tkinter padx参数，帮助大家更好的理解和学习,感兴趣的朋友可以了解下,希望能够给你带来帮助
    2021-10-10

最新评论