微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

使用shapefile剪辑或遮罩多波段图像

如何解决使用shapefile剪辑或遮罩多波段图像

下面的代码允许用户使用简单的多边形shapefile剪切图像(1个带)。但是,即使我尝试了几种变体,也无法将其应用于由多个波段组成的图像。我想做的是读取具有多个条带的Geotiff图像并将其全部裁剪,而无需将它们分成不同的图像。

from Rasterio.mask import mask
import geopandas as gpd

img = Rasterio.open("..\Image.tiff")
aoi = gpd.read_file("..\polyGON.shp")
geometry = aoi["geometry"]
clipImg,clipImg_Transform = mask(img,geometry,crop=True,invert=False)
clipImg_Meta = img.Meta.copy()
clipImg_Meta.update({"driver": "GTiff","height": clipImg.shape[1],"width": clipImg.shape[2],"transform": clipImg_Transform})

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。