python 如何在地图上绘制矩形框

python 如何在地图上绘制矩形框

通过官网的一些示例，对于处理海洋上的数据时，在地图上绘制矩形框，只需要增加一下投影，然后根据经纬度设置你的范围即可。

matplotlib.patches.Rectangle官网说明

matplotlib.patches.Rectangle(xy, width, height, angle=0.0, **kwargs)[source]

从函数可以看出来，主要就是设置（起始点，宽，高），所谓的起始点就是从左下角开始绘图的点，投影到地图上就是我们的（经度第一个值，纬度第一个值），宽度和高度就简单了，选取的经度范围和纬度范围罢了。下面这个图可以很形象的表示出来：

:                +------------------+
:                |                  |
:               lat                 |
:                |                  |
:     (lon[0]，lat[0])---- lon -----+

首先导入所需的包：

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

然后，随便导入你要进行绘制的数据经纬度，我这里使用的是SST（海表温度）:

path='D:\\sst.nc'
ds=xr.open_dataset(path).sortby("lat", ascending= True)
sst=ds.sel(lat=slice(-20,20),time=slice('2000','2010'))
# 选择矩形框的区域，矩形框的经纬度范围
region=sst.sel(lat=slice(-15,15),lon=slice(180,180+90)).sst
lon=region.lon.data
lat=region.lat.data
#定义一个绘图函数，避免以后绘制多张图时，重复输入代码
def make_map(ax, title):
    # set_extent  set crs
    ax.set_extent(box, crs=ccrs.PlateCarree())
    land = cfeature.NaturalEarthFeature('physical',
                                        'land',
                                        scale,
                                        edgecolor='grey', 
                                        facecolor='grey',
                                        zorder=2)
    ax.add_feature(land)  # set land color
    ax.coastlines(scale)  # set coastline resolution
    # set coordinate axis
    ax.set_xticks(np.arange(0, 360+45, 45),crs=ccrs.PlateCarree())
    ax.set_yticks(np.arange(-20, 30, 10),crs=ccrs.PlateCarree())
    ax.xaxis.set_major_formatter(LongitudeFormatter(zero_direction_label =False))#经度0不加标识
    ax.yaxis.set_major_formatter(LatitudeFormatter())
    # plt.tick_params(labelsize=25)
    lon_formatter = LongitudeFormatter(zero_direction_label=False)  
    # zero_direction_label set 0 without E or W
    lat_formatter = LatitudeFormatter()
    ax.xaxis.set_major_formatter(lon_formatter)
    ax.yaxis.set_major_formatter(lat_formatter)
    ax.set_title(title, fontsize=25, loc='center')
    ax.tick_params(which='major', 
                    direction='out', 
                    length=8,
                    width=0.99, 
                    pad=0.2, 
                    labelsize=20,
                    bottom=True, left=True, right=False, top=False)
    return ax
# 设置绘图函数所需的参数
box = [0, 361, -20, 20]  
scale = '50m'            
xstep, ystep = 45, 10   
cmap=plt.get_cmap('jet',9)#'RdYlBu_r'  'viridis'
proj=ccrs.PlateCarree(central_longitude=180)
box_top = 45
#设置矩形框的起始点
start_point = (lon[0], lat[0])
# 进行绘图
# # #===============================开始绘图==============================
fig=plt.figure(figsize=(20,12))
ax= fig.add_axes([0.1, 0.2, 0.5, 0.5],projection = proj)
make_map(ax,' draw rectangle')
#  #  # select 矩形区域
RE=Rectangle(start_point,lon[-1]-lon[0],lat[-1]-lat[0],linewidth=2,linestyle='-' ,zorder=2,\
              edgecolor='red',facecolor='none', transform=ccrs.PlateCarree()) #25表示纬度间隔25
ax.add_patch(RE)  #添加矩形框，如果这一句求掉了就没有框显示了

运行一下代码可以发现绘制成功啦~~非常简单

#当然，你可以进行下面的操作，使框内增加颜色
re=Rectangle(start_point,lon[-1]-lon[0],lat[-1]-lat[0],linewidth=2,linestyle='-' ,zorder=2,\
              edgecolor='red',facecolor='blue',alpha=0.5, transform=ccrs.PlateCarree()) 
ax.add_patch(re)  

ok，that’s all。
感兴趣的小伙伴们可以尝试。

										    		一个努力学习python的海洋菜鸡
										                 水平有限，欢迎指正！！！
										              欢迎关注、点赞、评论、收藏。