SAR地理图像中的陆地掩蔽

Question

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

Answer 1

我找到了出路：

首先，我们必须下载我们想要屏蔽的区域的相应shapefile， 然后在gdal中有一个漂亮的功能，称为gdalwarp。我们需要打开 anaconda提示，然后输入：`

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

现在，带有边界线的图像将保存在文件output_filename.tif中 这是包含陆地部分的文件，海洋被掩盖了。

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

我们将获得海洋部分的图像，其中陆地部分为黑色，之后我们可以将陆地部分设为NaN。