关于maskout白化之后放大区域的问题

BrandNewMe · 发表于 2021-8-12 15:52:13

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 BrandNewMe 于 2021-8-12 15:54 编辑

我用中国区域做了maskout，但是想关注环渤海，因此将图片范围聚焦于环渤海。没有maskout的时候如下图是正常的，mask之后就成上图这样了，请问有没有解决的方法？

画图部分脚本如下：

fig = plt.figure(figsize=(15, 12))
ax = plt.axes(projection=ccrs.PlateCarree())
ax.set_extent([100, 135, 30, 45])
contours = ax.contourf(lon_grid, lat_grid, first_guess, levels=range(0,30001,1000),transform=ccrs.PlateCarree(),cmap=cmaps.MPL_terrain_r)
clip=maskout.shp2clip(contours,ax,'country1.shp','China')

province = shpreader.Reader(r'bou2_4l.shp')
ax.add_geometries(province.geometries(), crs=ccrs.PlateCarree(), linewidths=0.5,edgecolor='k',facecolor='none')
gl = ax.gridlines(draw_labels=True, linewidth=1, color='k', alpha=0.5, linestyle='--')
gl.top_labels = False  #关闭顶端标签
gl.right_labels = False  #关闭右侧标签
gl.xformatter = LONGITUDE_FORMATTER  #x轴设为经度格式
gl.yformatter = LATITUDE_FORMATTER  #y轴设为纬度格式
cbar = fig.colorbar(contours,shrink=1)
cbar.set_ticks(np.arange(0,30001,1000))
ax.set_title('OBS_Visability  valid at 2019-01-12T00:00:00(UTC)',fontsize=15)
plt.show()