关于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()

edwardli · 发表于 2021-8-20 16:29:00

把Cartopy降级到0.17.0即可。

BrandNewMe · 发表于 2021-8-12 16:01:35

本帖最后由 BrandNewMe 于 2021-8-12 16:04 编辑

数据做了经纬度截取后成了这样，图片不是最大，只是框外的填图没有了。

BrandNewMe · 发表于 2021-8-23 15:14:57

edwardli 发表于 2021-8-20 16:29
把Cartopy降级到0.17.0即可。

谢谢！我试一试

mikado · 发表于 2021-9-15 16:23:23

本帖最后由 mikado 于 2021-9-16 09:49 编辑

Cartopy降级到0.17.0能画了。。。

BrandNewMe · 发表于 2021-9-16 11:03:12

mikado 发表于 2021-9-15 16:23
Cartopy降级到0.17.0能画了。。。

谢谢，确实能用了~

		自动登录	找回密码
密码			立即注册

[求助] 关于maskout白化之后放大区域的问题

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