Python完美白化

kkkk9090

求助大神，bou2_4p这个shp你有没有，能否共享一下，谢谢，我只有全国的底图才用得起，就是想把四川以外的都白化了

Mr.Gao

平流层的萝卜 发表于 2019-4-25 04:12
背景的颜色应该是可以设置的

请问是在maskout里面吗

hailonghelen

神贴啊，我要好好学习

jinhe

希望能学习一下插值边界处理能够下载一下文件

平流层的萝卜

Mr.Gao 发表于 2019-4-25 14:48
请问是在maskout里面吗

不再maskout里，应该是pyplot，或者basemap

dongfangxing

学习和互鉴提高，白化是画图基本功，{:eb502:}{:eb502:}

sam_doggy

楼主，以后可能专用cartopy，最近从basemap改学，然后这个clip文件，如果等经纬不用动，如果是其他投影，原来传递basemap的参数的，现在传递什么啊

妮妮

用了楼主的maskout,但是一直 出现这样的问题，不知道怎么回事，可以麻烦楼主帮忙看一下吗

平流层的萝卜

妮妮 发表于 2019-6-3 10:51
用了楼主的maskout,但是一直 出现这样的问题，不知道怎么回事，可以麻烦楼主帮忙看一下吗

你好。
提示错误是index超出了shapfile的范围。如果你是用其他shp文件的话，需要重新确认一下，想要mask地区的特殊标识符是什么，换句话说，就是用那个属性等于哪个值来表示mask地区。

xc123

请问为什么我出现如下错误呢？
local variable 'clip1' referenced before assignment

代码如下：
fig = plt.figure(figsize=[10, 8])
    # Set projection and plot the main figure
    ax = plt.axes(projection=ccrs.LambertConformal(central_latitude=30.5,
                                                   central_longitude=118.5))
    border = r"d:\forcast\shp\宣城边界.shp"
    shape_feature = cfeature.ShapelyFeature(Reader(border).geometries(),ccrs.PlateCarree(),edgecolor='k',facecolor='none')
    ax.add_feature(shape_feature)
    ax.gridlines(linestyle='--')
    # Set figure extent
    ax.set_extent([117.8, 119.8, 29.8, 31.5])
   
   
    #==============================================================================
    #set colormap
    levels=[0,0.1,10,25,50,100,250,500]#雨量间隔
    colors=['#FFFFFF','#A5F38D','#3DB93F','#63B8F9','#0000FE','#F305EE','#810040']#雨量标准色标

    cf=plt.contourf(lons,lats,f,15, levels=levels, colors=colors,transform=ccrs.PlateCarree())
    plt.colorbar()
    shijian=str(d1.year)+'年'+str(d1.month)+'月'+str(d1.day)+'日'+str(d1.hour)+'时'+'--'+str(d2.year)+'年'+str(d2.month)+'月'+str(d2.day)+'日'+str(d2.hour)+'时'+'雨量预报'
    plt.title(shijian,fontsize=20,color='blue')
    d1=d2
    clip2=maskout.shp2clip(cf,ax,'d:\\forcast\\shp\\xc.shp',[341881]) #白化
    plt.savefig(outfile)
    plt.show()