求助各位大神，只想画出相对湿度大于等于70%的区域，代码部分该怎么写啊

Tobias

求助各位大神，只想画出相对湿度大于等于70%的区域，代码部分该怎么写啊。
尝试了levels=arange(70,101,10),rh=contourf(rh,levels),但这样似乎没用，有没有大神给指导一下啊

MeteoInfo

contourf 函数中要加上设置的 levels 参数

Tobias

MeteoInfo 发表于 2021-10-26 11:12
contourf 函数中要加上设置的 levels 参数

老师，我刚把levels加进去了，画出来的图70%以下都是蓝色，有没有可能向micaps一样，让70%以下的区域不填色，显示白色

气象小白01

把70%的mask掉

气象小白01

或者手动把70%以下的值设置成NaN

Tobias

气象小白01 发表于 2021-10-26 13:20
或者手动把70%以下的值设置成NaN

好的，我试试

MeteoInfo

Tobias 发表于 2021-10-26 13:28
好的，我试试

可以将第一个颜色设为白色 [255,255,255]， 或者颜色第四个参数alpha设为 0 (全透明)。参考下面的脚本：

1. fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')
   
2. f = addfile(fn)
   
3. ps = f['PS'][0,'10:60:8','60:140:10']
   
4. 
   
5. #Plot
   
6. geoshow('country', edgecolor=(0,0,255))
   
7. levs = arange(600, 1051, 30)
   
8. cols = makecolors(len(levs)+1)
   
9. cols[0] = [255,255,255,0]
   
10. layer = contourf(ps, levs, colors=cols)
    
11. title('Pressure')
    
12. colorbar(layer, orientation='horizontal', extendrect=False, shrink=0.8,
    
13.     aspect=20, label='hPa', labelloc='right')

fn = os.path.join(migl.get_sample_folder(), 'GrADS', 'model.ctl')<br /> f = addfile(fn)<br /> ps = f['PS'][0,'10:60:8','60:140:10']<br /> <br /> #Plot<br /> geoshow('country', edgecolor=(0,0,255))<br /> levs = arange(600, 1051, 30)<br /> cols = makecolors(len(levs)+1)<br /> cols[0] = [255,255,255,0]<br /> layer = contourf(ps, levs, colors=cols)<br /> title('Pressure')<br /> colorbar(layer, orientation='horizontal', extendrect=False, shrink=0.8, <br />     aspect=20, label='hPa', labelloc='right')

Tobias

MeteoInfo 发表于 2021-10-26 16:51
可以将第一个颜色设为白色 [255,255,255]， 或者颜色第四个参数alpha设为 0 (全透明)。参考下面的脚本：
...

收到，谢谢老师