python遥感图像处理应用篇(二十四):python绘制遥感图像各波段热力图(相关系数矩阵)(代码片段)

空中旋转篮球 空中旋转篮球     2022-12-11     288

关键词:

给多光谱遥感图像各个波段绘制热力图,首先需要计算波段之间的相关系数矩阵,而计算遥感图像波段相关系数矩阵有不同的方法,常用的我们可以采用遥感图像处理软件计算,比如ENVI软件就可以计算相关系数矩阵,使用工具箱中的Statistics工具即可进行多种统计运算。

我们这里不使用遥感图像软件计算,直接使用Python计算。要实现整个目标,我们首先需要提取遥感图像波段中各个像元的信息。这里我们可以将遥感图像读取并保存为csv文件,实现这一步我们可以参照之前的一篇文章,链接如下:

Python遥感图像处理应用篇(十四):GDAL 读取多光谱数据为二维数组并存入csv文件_空中旋转篮球的博客-CSDN博客

这里不做赘述了,直接跳过这一部,我们使用pandas直接读取csv文件计算相关系数。

Python使用到的库文件:

import seaborn as sns
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

pandas读取文件,并对文件进行截取,选取目标波段:

# 输入数据
data = pd.read_csv(filepath)
data = data.iloc[:,1:10]  #2-10列,根据自己需求截取

计算相关系数矩阵&#

python遥感图像处理应用篇(二十二):python+gdal批量等距离裁剪影像-续

...数量来进行影像等距离裁剪的博客,链接如下:Python遥感图像处理应用篇(二十二):Python+GDAL批量等距离裁剪影像_空中旋转篮球的博客-CSDN博客_python遥感图像处理本篇对上一篇博客代码进行少量修改,实现遥感图像按照指定大小... 查看详情

python遥感图像处理应用篇(二十二):python+gdal批量等距离裁剪影像-续

...数量来进行影像等距离裁剪的博客,链接如下:Python遥感图像处理应用篇(二十二):Python+GDAL批量等距离裁剪影像_空中旋转篮球的博客-CSDN博客_python遥感图像处理本篇对上一篇博客代码进行少量修改,实现遥感图像按照指定大小... 查看详情

python遥感图像处理应用篇(二十八):python绘制遥感图像分类结果混淆矩阵和计算分类精度

1.使用数据IndiansPines高光谱数据,使用SVM分类方法(选取10%样本量)计算得到的结果。参考数据:分类数据: 查看详情

python遥感图像处理应用篇(二十八):python绘制遥感图像分类结果混淆矩阵和计算分类精度

1.使用数据IndiansPines高光谱数据,使用SVM分类方法(选取10%样本量)计算得到的结果。参考数据:分类数据: 查看详情

python遥感图像处理应用篇(二十八):python绘制遥感图像分类结果混淆矩阵和计算分类精度

1.使用数据IndiansPines高光谱数据,使用SVM分类方法(选取10%样本量)计算得到的结果。参考数据:分类数据: 查看详情

python遥感图像处理应用篇(二十五):python+gdal波段组合(代码片段)

...默认排序,读取的时候也是按照排列顺序读取。Python遥感图像处理基础篇(一):arcpy遥感图像波段组合批量处理_空中旋转篮球的博客-CSDN博客_python处理遥感 2.实现思 查看详情

python遥感图像处理应用篇(二十五):python+gdal波段组合(代码片段)

...默认排序,读取的时候也是按照排列顺序读取。Python遥感图像处理基础篇(一):arcpy遥感图像波段组合批量处理_空中旋转篮球的博客-CSDN博客_python处理遥感 2.实现思 查看详情

python遥感图像处理应用篇(二十九):遥感图像拉伸处理

1.实验目的将原始遥感图像各个波段拉伸到0-255范围显示。运行环境:windows10pycharmpython3.7.7GDAL-3.2.3-cp37-cp37m-win_amd64.whl包等。2.使用公式x_nor=((data-d2)/(d98-d2))*255x_nor为新的像元值,data为原始像元值,d2为2%像元值,d98为98%像元值,255... 查看详情

python遥感图像处理应用篇(二十二):python+gdal批量等距离裁剪影像

...如果需要按照指定大小进行影像裁剪可以参考:Python遥感图像处理应用篇(二十二):Python+GDAL批量等距离裁剪影像-续_空中旋转篮球的博客-CSDN博客2.实现代码运行环境:运 查看详情

python遥感图像处理应用篇(二十七):python绘制遥感图像各波段热力图(相关系数矩阵)(续)

续-https://soderayer.blog.csdn.net/article/details/125757807上一篇中使用csv文件计算的相关系数热力图,本篇我们直接使用遥感图像来计算图像波段之间的相关系数。方法一:已有软件ENVI计算实际上,目前已有的软件,如ENVI就可以直接计算... 查看详情

python遥感图像处理应用篇(二十六):python+gdalsentinel-2数据波段组合(代码片段)

...有tif文件放置在一个文件夹种实现的波段组合。Python遥感图像处理应用篇(二十五):Python+GDAL波段组合文件路径很多,我们需要先读取到需要组合波段的完整路径。E:\\S2AB_2020_hubei\\S2A_MSIL2A_20200817T0255 查看详情

python遥感图像处理应用篇(二十六):python+gdalsentinel-2数据波段组合(代码片段)

...有tif文件放置在一个文件夹种实现的波段组合。Python遥感图像处理应用篇(二十五):Python+GDAL波段组合文件路径很多,我们需要先读取到需要组合波段的完整路径。E:\\S2AB_2020_hubei\\S2A_MSIL2A_20200817T0255 查看详情

python遥感图像处理应用篇(二十一):python+gdal批量计算遥感图像ndvi指数(代码片段)

1.NDVI指数计算公式之前也写过NDVI批量计算的实现方法,采用的是Acrpypython2.7实现的,这里我们采用GDAL实现NDVI指数的批量计算。使用数据,多波段遥感图像数据。Landsat05Colection2Level2数据。NDVI=(NIR-RED)/(NIR+RED) 就算之前需要搞清... 查看详情

python遥感图像处理应用篇(二十一):python+gdal批量计算遥感图像ndvi指数(代码片段)

1.NDVI指数计算公式之前也写过NDVI批量计算的实现方法,采用的是Acrpypython2.7实现的,这里我们采用GDAL实现NDVI指数的批量计算。使用数据,多波段遥感图像数据。Landsat05Colection2Level2数据。NDVI=(NIR-RED)/(NIR+RED) 就算之前需要搞清... 查看详情

python遥感图像处理应用篇(二十三):python+gdal批量拼接图像(代码片段)

遥感影像拼接分多种不同的情况,比如比较常用的是两幅影像拼接或者多个图像合并为一个图像。GDAL中实现影像拼接的方式也有多种,比如常用可以通过numpy读取影像数组,并计算各个图像的numpy数组范围,合并多个numpy数组为... 查看详情

python遥感图像处理应用篇(二十三):python+gdal批量拼接图像(代码片段)

遥感影像拼接分多种不同的情况,比如比较常用的是两幅影像拼接或者多个图像合并为一个图像。GDAL中实现影像拼接的方式也有多种,比如常用可以通过numpy读取影像数组,并计算各个图像的numpy数组范围,合并多个numpy数组为... 查看详情

python遥感图像处理应用篇(二十):python+gdal批量提取多波段图像为单波段图像(代码片段)

1.实现目标批量将多波段图像提取为各个单波段图像运行环境:windows10 pycharm  python3.7.7  GDAL-3.2.3-cp37-cp37m-win_amd64.whl包等。2.实现代码"""此段代码将遥感图像背景值去除处理"""importnumpyasnpfromosgeoimportgdal,gdalconstimpo 查看详情

python遥感图像处理应用篇(二十):python+gdal批量提取多波段图像为单波段图像(代码片段)

1.实现目标批量将多波段图像提取为各个单波段图像运行环境:windows10 pycharm  python3.7.7  GDAL-3.2.3-cp37-cp37m-win_amd64.whl包等。2.实现代码"""此段代码将遥感图像背景值去除处理"""importnumpyasnpfromosgeoimportgdal,gdalconstimpo 查看详情