介绍一种绘图白化的方法

po_po1

Byby 发表于 2018-11-25 17:23
楼主楼主 我最近有个脚本也想用白化。我的数据是中国及周边近海的，我想把陆地的数据掩盖掉，只显示海洋 ...

或者：

1. import numpy as np
   
2. import matplotlib.pyplot as plt
   
3. import cartopy as cart
   
4. from mpl_toolkits.basemap import Basemap
   
5. from netCDF4 import Dataset
   
6. 
   
7. f = Dataset('sst.mnmean.nc')
   
8. sst = f.variables['sst'][0,:,:]
   
9. lats = f.variables['lat'][:]
   
10. lons = f.variables['lon'][:]
    
11. 
    
12. ax = plt.axes(projection=cart.crs.PlateCarree())
    
13. ax.coastlines()
    
14. ax.add_feature(cart.feature.LAND, zorder=100, edgecolor='k')
    
15. ax.set_global()
    
16. plot = ax.contourf(lons, lats, sst, 60, transform=cart.crs.PlateCarree())
    
17. cb = plt.colorbar(plot)
    
18. plt.show()
    

复制代码

改变cart.feature.LAND 为cart.feature.OCEAN，mask 海洋

Byby

po_po1 发表于 2018-11-26 14:28
或者：

excellent 可行  通过画图顺序进行掩盖 多谢楼主

sres

这一步的时候显示错误，这是为啥呢
main_us_geom = sorted(us_multipoly.geoms, key=lambda geom: geom.area)[-1]

AttributeError: 'Polygon' object has no attribute 'geoms'

chongzika

po_po1 发表于 2014-6-5 11:18
附带楼主上次提到的那个图和程序：

下面是程序：

最近又重新试了下这个mask方法，但是发现，为什么海南岛和台湾的mask后没有值？是shp文件的问题么？

po_po1

chongzika 发表于 2019-5-18 23:18
最近又重新试了下这个mask方法，但是发现，为什么海南岛和台湾的mask后没有值？是shp文件的问题么？

不是ｓｈｐ文件有问题，只是这个方法仅限于简单区域

zhongyon

楼主，有个问题请教下您。我裁剪之后，想把其他区域设置成背景透明，只显示我需要的那一部分怎么弄呢。按照您的步骤，我裁剪之后其他区域是白色，怎么设置成透明呢？

汐年

您好，我想请问下，我想对印度洋地区进行白化，但是在用您的代码时，这一句总是报错country.geometry for country in countries  if country.attributes['NAME'] == 'China']，就是其中的NAME，我想问下这个是什么原因呢？

tanyueye

sres 发表于 2019-1-9 11:20
这一步的时候显示错误，这是为啥呢
main_us_geom = sorted(us_multipoly.geoms, key=lambda geom: geom.ar ...

同问，是不是地图数据一定要Micaps的呢？
还有个问题是如果shp里面只有一个区域，us_multipoly, = [province.geometry for province in provinces  if province.attributes['NAME']==‘China’] 这里提示报错（如果是多个区域就可以）。

prn

必须称赞，感谢分享。

novelbean

如果 country.attributes['NAME'] == 'China'类似的属性没有怎么办，只有Id，并且Id的值全是0，那则么整？