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：根据以上四至点坐标值计算经纬度为单位的边界范围

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

Python遥感图像处理应用篇(十二)：Arcpy批量计算比值植被指数RVI
1.比值植被指数RVI 公式:RVI=TM4/TM3(Landst5-7) RVI=TM5/TM4(Landst8-9) RVI=NIR/R NIR是近红外波段,R是红光波段的反射率.对于浓密植物 ...
Python遥感图像处理应用篇(十六)：GDAL 将归一化处理csv数据转化为遥感影像
1.使用数据将上一篇文章中得到的计算结果作为转换数据.链接如下: Python遥感图像处理应用篇(十五):GDAL 读取多光谱数据做归一化处理_空中旋转篮球的博客-CSDN博客 2.实现代码基本思 ...
Python遥感图像处理应用篇(十八)：GDAL +numpy对遥感图像归一化处理
1.使用数据所用输入数据为遥感影像数据,可以是单波段也可以是多波段数据. 运行环境:windows10 pycharm python3.7.7 GDAL-3.2.3-cp37-cp37m-wi ...
Python遥感图像处理应用篇(十四)：GDAL 读取多光谱数据为二维数组并存入csv文件
1.使用数据之前写过一遍文章是采用Arcgis读取多光谱数据为csv文件的,相对来说有一些麻烦,本文有个更简单一点的方法,使用python直接读取. 本文使用Landsat05数据,采用GEE下载的 ...
Python遥感图像处理应用篇(十)(续):使用EVI指数批量计算叶面积指数LAI
1.叶面积指数概念叶面积指数(leaf area index, LAI)不同资料显示定义有所不同,有的说是指单位地面上的绿叶面积,是植被冠层最显著的特征之一,具体是指一定土地面积上植物叶面面积总和与 ...
Python遥感图像处理应用篇(十)：使用NDVI指数批量计算叶面积指数LAI
1.叶面积指数概念叶面积指数(leaf area index, LAI)不同资料显示定义有所不同,有的说是指单位地面上的绿叶面积,是植被冠层最显著的特征之一,具体是指一定土地面积上植物叶面面 ...
Python遥感图像处理应用篇(二十四)：Python绘制遥感图像各波段热力图（相关系数矩阵）
给多光谱遥感图像各个波段绘制热力图,首先需要计算波段之间的相关系数矩阵,而计算遥感图像波段相关系数矩阵有不同的方法,常用的我们可以采用遥感图像处理软件计算,比如ENVI软件就可以计算相关系数矩阵,使用 ...
Python遥感图像处理应用篇(二十二)：Python+GDAL 批量等距离裁剪影像-续
之前写过一篇按照指定行列号数量来进行影像等距离裁剪的博客,链接如下: Python遥感图像处理应用篇(二十二):Python+GDAL 批量等距离裁剪影像_空中旋转篮球的博客-CSDN博客_pytho ...
Python遥感图像处理应用篇(二十五)：Python+GDAL 波段组合
1.使用场景描述之前使用arcpy python2.7写了一篇进行遥感数据波段组合计算的博文,使用场景是将所有单波段数据放在一个文件夹中(如下图),文件名称前缀一样,后缀波段编号用来区分不同的波段, ...

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

1.问题描述

2.实现思路

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

最新文章

热门文章