SAR geotiff 图像中的 Landmasking

Question

我正在尝试在卫星 (SAR) 灰度 geotiff 图像中掩盖土地。该功能在 rsgislib 中可用，但它适用于 Linux，我正在使用 conda python 3.5 (Windows) 并且无法找到可能的出路。 请指导如何在图像中掩盖土地。

Answer 1

我找到了出路：

首先，我们必须下载我们想要遮罩的区域的适当形状文件， 然后在 gdal 中有一个漂亮的功能，称为gdalwarp。我们只需要打开 anaconda prompt 并从那里输入：`

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

现在，带有土地边界的图像将保存在文件output_filename.tif
这是包含陆地部分和海洋的文件。

然后，通过从原始图像中减去output_filename.tif 图像，该过程变得相当简单。

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