Python 3.x SAR-geotiff图像中的地形遮挡

我试图在卫星（SAR）灰度geotiff图像中遮罩陆地。该功能在rsgislib中可用，但它可以在Linux上运行，我正在使用conda python 3.5（Windows），无法找到可能的解决方法。 请指导如何在图像中掩盖土地

我找到了出路：

首先，我们必须下载一个合适的形状文件的区域，我们希望遮罩， 然后，gdal中有一个漂亮的功能，称为

gdalwarp

。我们只需打开anaconda提示符，然后从那里输入：`

gdalwarp -cutline shapefile_name.shp original_image.tif output_filename.tif

现在，带有土地边界的图像将保存在文件

output\u filename.tif

这是一个包含陆地部分和海洋部分的文件

然后，通过从原始图像中减去

output_filename.tif

图像，该过程变得相当简单

我们将得到海洋部分的图像，陆地部分为黑色，然后我们可以将陆地部分制作为

NaN