Python白化截取研究区域

zzw9701

仿照Python完美白化-编程作图-气象家园_气象人自己的家园 (06climate.com)的方法对地图进行白化，只想对中国的部分区域进行白化，但是最后好像把整个wrf文件包含的区域给画出来了，请问怎样才能只画矩形里的图呢？

#绘图
ax1 = fig.add_axes([0.1, 0.3, 0.4, 0.4],projection = proj) #左 底 宽 高
m = Basemap(projection='cyl', llcrnrlat=lowerlat, urcrnrlat=upperlat, llcrnrlon=rightlon, urcrnrlon=leftlon, resolution='l')
m.readshapefile('D:\\Compressed\\country\\country1', 'country1,shp', color='black',linewidth=0.3)
m.drawcoastlines(linewidth=0.3, color='black')
rh=ax1.contourf(lon,lat,tot_rain,transform=proj,levels=bounds,norm=norm,extend='both',cmap=plt.cm.GnBu)
cb=fig.colorbar(rh,pad=0.07,shrink=0.7,aspect=25)#方向
#地图白化
clip = maskout.shp2clip(rh,ax1,m,'D:\\Compressed\\country\\country1',["China"])
xstep, ystep = 2, 2
ax1.set_extent(img_extent,crs=proj)
ax1.set_xticks(np.arange(leftlon, rightlon+xstep,xstep), crs=proj)
ax1.set_yticks(np.arange(lowerlat, upperlat+ystep,ystep), crs=proj)
lon_formatter = LongitudeFormatter(zero_direction_label=False)
lat_formatter = LatitudeFormatter()
ax1.xaxis.set_major_formatter(lon_formatter)
ax1.yaxis.set_major_formatter(lat_formatter)
ax1.set_title('(a) Total Rain',loc='left',fontsize =15,fontname ='Times New Roman')

werewolf

数据做个经纬度范围选取或者画图范围改一下

likailing2008

同楼上方法，或者你可以做一个经纬度范围的shp