import regionmask shp = gp.read_file("shp/Uzb.shp") mask = regionmask.mask_geopandas(shp, lons, lats) def makemask(data): return np.ma.masked_array(data, mask=mask) fig2 = plt.figure() proj = ccrs.PlateCarree(central_longitude=lons[int((len(lons)+1)/2)]) #设置一个圆柱投影坐标...
mask = regionmask.defined_regions.natural_earth_v5_0_0.land_50.mask(precip) p_land = precip.where((mask == 0)) 对于海洋进行掩膜 maskocean = regionmask.defined_regions.natural_earth_v5_0_0.ocean_basins_50.mask(precip) p_ocean = precip.where((maskocean >= 1)) 绘图展示掩膜结果 fig =...
import geopandas as gp import regionmask shp = gp.read_file("shp/Uzb.shp") mask = regionmask.mask_geopandas(shp, lons, lats) def makemask(data): return np.ma.masked_array(data, mask=mask) 1. 2. 3. 4. 5. 6. 7. 8. fig2 = plt.figure() proj = ccrs.PlateCarree(central_longit...
方法1：global_land_mask 该方法使用了global_land_mask库。通过安装并导入此库，定义了一个函数mask_land，将输入数据集进行陆地掩膜。该函数通过globe对象的is_ocean函数判断经纬度网格点是否在海洋上。适用于Netcdf格式的文件，其他格式可能需调整。方法2：regionmask regionmask库提供了海洋区域轮廓，从Na...
savefig("Himalayas_mask.png", dpi=500, bbox_inches='tight') plt.show() 出图: 4.中国区域白化(包含南海子图) 代码语言:javascript 复制 #https://regionmask.readthedocs.io/en/stable/defined_landmask.html #coding=utf-8 from mpl_toolkits.basemap import Basemap from matplotlib.path import Path ...
Natural Earth掩膜海陆 利用regionmask预设的natural erath陆地信息提取相关值,达到掩膜海洋的作用,但这里的分辨率只有1:110,比较粗糙。 利用regionmask预设的natural erath海盆信息提取相关值,达到掩膜陆地的作用,默认分辨率为1:50。 自定义shp进行掩膜
# 需要导入模块: import geopandas [as 别名]# 或者: from geopandas importGeoSeries[as 别名]deftest_network_false_nodes(self):test_file_path2 = mm.datasets.get_path("tests") self.false_network = gpd.read_file(test_file_path2, layer="network") ...
- global-land-mask==1.0.0 - munch==4.0.0 - pillow==9.1.1 - pooch==1.8.1 - pyogrio==0.7.2 - pyproj==3.3.1 - rasterio==1.2.10 - regionmask==0.12.1 - salem==0.3.5 - shapely==1.8.2 prefix: D:\anaconda\ANACONDA\envs\py311...