Python遥感图像处理应用篇(十九)：GDAL +numpy批量对遥感图像外围背景值进行处理

Posted 2022-06-29 空中旋转篮球

1.问题描述

最近下载了一些遥感影像数据，这些数据都包含大量的外围背景数据，如下图所示：

外围背景值都为0值。

本文描述的是采用python批量处理外围背景，只保留最小外围背景区域。

如下图：

2.实现思路

基本思路如下：

首先获取遥感影像有数据区域四个角的坐标，这里的坐标我们用行列号表示，即左上角（x1,y1),右上角（x2,y2），左下角（x3,y3），右下角（x4,y4）。我们需要通过遍历遥感数据的方法获取到四个角点的行列号。

然后，根据四角坐标计算四至范围，即left，top，right，bottom值。

方法一：之后，分别获取左上角坐标和右下角坐标值:(left,top),(right,bottom)。最后读取(left,top),(right,bottom)范围内的数据，重新写入并保存到新的影像中。这样计算出来的结果有个问题，影像的坐标范围会发生变化。因此不可取，应该计算为经纬度坐标范围，然后根据经纬度范围裁剪。

方法2：根据以上四至点坐标值计算经纬度为单位的边界范围