Python cartopy 绘图 显著性打点

import xarray as xr
import cartopy.crs as ccrs
import matplotlib.pyplot as plt
import numpy as np
tif1 = r'.tif'
data = xr.open_dataarray(tif1)
fig = plt.figure(figsize=(5,5))
ax11 = fig.add_axes([0.05,0.05,0.4,0.42], projection =ccrs.PlateCarree())
trend = data[0].plot.imshow(ax=ax11,cmap='RdYlGn',add_colorbar=False,
                                  transform=ccrs.PlateCarree(),
                                  vmax=1,vmin=-1,
                                  zorder=1)
#0.26为一个阈值,根据自己的需求设置即可
dot_area = np.asarray(np.abs(data.values)>0.26).reshape(279,618)   #279,618 矩阵行列号
mesh_lon, mesh_lat = np.meshgrid(data['x'], data['y'])
ax11.scatter(mesh_lon[dot_area],mesh_lat[dot_area],transform=ccrs.PlateCarree(), color='k',marker='x',s=0.05,linewidth=0.1,zorder=2)
ax11.set_title('')
plt.show()

喜欢的点个赞哦!关注一下

  • 2
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值