python使用外部.shp进行水陆掩膜

本人纯属初学者，如有表达不当请在评论区纠正我

一、掩摸之前需要知道：读tif、读shp、掩摸、掩摸取反

1、读tif、shp

python中gdal与rasterio库均可以读取tif，但是rasterio库中的操作更简单一些。

下面链接是rasterio包的简单学习python 栅格处理利器之Rasterio - 知乎 (zhihu.com)

同时附上rasterio学习文档栅格：访问地理空间栅格数据 — rasterio 文档 (osgeo.cn)

读取tif：src_tif是原始影像，dst_tif是最后掩摸之后的影像

src_img = r'C:\Users\Dell\Desktop\COMS_GOCI_L1B_GA_20200429001641.he5\res1_tif.tif'#原始tif图
dst_img = r'F:\data\goci\goci-2016\201606\06-result-mask\L2_res.tif'#最终掩摸之后的文件

读取shp及其操作：shapes是关于shp文件要素的集合，此处注意一定要将shp和tif投影参数设置一致。

mask = r'F:\常用数据\2019全国行政区划\2019全国行政区划\市.shp'
maskdata = GeoDataFrame.from_file(mask)
with fiona.open(mask, 'r') as shapefile:
    shapes = [feature["geometry