当前位置：首页 > news >正文

landsat卫星遥感影像下载、处理教程

news 来源：原创 2025/5/23 5:18:37

1. landsat数据下载

USGS网址：
EarthExplorer

参考链接：

USGS下载遥感影像——以Landsat影像下载为例_usgs怎么下载遥感影像-CSDN博客

L1TP数据进行了几何校正和辐射校正，L2SP数据在此基础上，进一步处理后的数据，通常包括地表反射率和其他相关的地表特征信息。但是L1和L2的选择需要根据实际需求。

这里我们选择下载landsat8-9 L1数据

目前是可以直接在浏览器下载，无须发送至邮箱。

2. 安装ENVI软件-QGIS

参考链接：

ENVI5.6详细安装教程，亲测可用（附安装包和破解文件）-CSDN博客

3. 辐射定标和大气校正

3.1 修改头文件

需要先修改头文件，否则导入ENVI会失败。头文件指的是下载的landsat文件夹中的以”_MTL.txt”结尾的文件。

修改Landsat8头文件的参考链接：ENVI无法打开Landsat8的头文件问题和解决_envi出现cannotread tiff header-CSDN博客

Landsat9的头文件修改方式：在landsat8修改方式的基础上，将头文件中
SPACECRAFT_ID = "LANDSAT_9" 修改为 SPACECRAFT_ID = "LANDSAT_8"
参考链接：【ENVI】Landsat9影像数据打开方式（envi）_envi打开landsat9图像-CSDN博客

3.2 辐射定标和大气校正

辐射定标和大气校正的参考链接：
ENVI5.3.1使用Landsat 8影像进行辐射定标和大气校正实例操作_使用郑州地区的landsat-8遥感影像数据,完成郑州地区植被覆盖度计算: (1)辐射校正:-CSDN博客

辐射定标后的文件命名为：adiometric_result
大气校正后的文件命名为：FLAASH_result

4. 裁剪图像

在ENVI中裁剪图像分为两种方式：导入shp文件或者制作ROI。

4.1.导入shp文件

4.1.1 矢量边界Geojson获取
矢量边界Geojson获取方式：DataV.GeoAtlas地理小工具系列
使用阿里云DataV.GeoAtlas工具（范围选择/边界生成/层级生成）使用省/市/县区的矢量边界或者手动绘制矢量边界，下载geojson文件。
4.1.2 Geojson文件转为shp文件
使用mapshaper工具将geojson文件转为shp文件
网址：mapshaper

4.2.制作ROI

参考链接：ENVI5.3制作ROI并保存为shp全流程_roi转shp-CSDN博客
（只需要看免费内容的步骤即可）

4.3 使用shp或者ROI裁剪

在Toolbox中搜索Subset Data from ROIs工具，双击，选择要被裁剪的图像，
选择ROI或者shp文件，mask pixels outside of ROIs选择yes，下方的值设置为0，选择文件导出的位置和文件名。

4.4 xml转ROI

参考链接：

ENVIroi标定结果输出为shp文件步骤_enviroi转shp-CSDN博客

5. 计算各种遥感指数

5.1 NDVI

5.1.1 计算方式1
ENVI 5.3 计算NDVI的几种方法及结果的差异分析_ndvi计算-CSDN博客
第一种方式的公式为：(float(b5)-float(b4))/(float(b5)+float(b4))

5.2 WET

参考上方NDVI链接中的第一种计算方式（使用Band Math工具）
计算公式为：
(float(b2)*0.1509+float(b3)*0.1973+float(b4)*0.3279+float(b5)*0.3406+float(b6)*(-0.7112)+float(b7)*(-0.4572))/10000
5.3 NDBSI

参考上方NDVI链接中的第一种计算方式（使用Band Math工具）
计算公式为：
(((2*float(b6)/(float(b6)+float(b5))-(float(b5)/(float(b5)+float(b4))+float(b3)/(float(b3)+float(b6))))/(2*float(b6) /(float(b6)+float(b5))+(float(b5)/(float(b5)+float(b4))+float(b3)/(float(b3)+float(b6)))))+(((float(b6)+float(b4)) - (float(b5)+float(b2)))/((float(b6)+float(b4)) + (float(b5)+float(b2)))))/2