Python遥感图像处理应用篇(十四)：GDAL 读取多光谱数据为二维数组并存入csv文件

已于 2022-10-11 09:22:59 修改
阅读量2.3k 收藏 21
点赞数 5
分类专栏： Python遥感图像处理 文章标签： 图像处理 GDAL
于 2022-06-13 12:07:06 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/soderayer/article/details/125256957
版权
本文介绍了一种使用Python的GDAL库直接读取Landsat 05多光谱数据并将其转换为二维数组，然后存储到CSV文件的方法。详细讲述了如何通过GDAL读取数据，创建CSV文件并写入像元ID，以及如何将各个波段数据写入CSV文件。同时，还探讨了代码优化和批量处理的实现。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

1.使用数据

之前写过一遍文章是采用Arcgis读取多光谱数据为csv文件的，相对来说有一些麻烦，本文有个更简单一点的方法，使用python直接读取。

本文使用Landsat05数据，采用GEE下载的多光谱数据，该数据下载下来包含19个波段。这些波段信息可以查看数据该数据详细说明。1-7是其对应的多光谱波段，第六波段是地表温度数据（ST),其他的是(SR) 地表反射率产品数据。数据显示如下（321波段）:

2.gdal读取数据

2.1 使用的库

from osgeo import gdal
import pandas as pd

gdal库的安装可以参考：python中如何导入gdal包？_空中旋转篮球的博客-CSDN博客_gda

了解本专栏 订阅专栏 解锁全文 超级会员免费看
Python遥感图像处理应用篇(十六)：GDAL 将归一化处理csv数据转化为遥感影像
soderayer的博客
06-16 1543
将上一篇文章中得到的计算结果作为转换数据。链接如下：Python遥感图像处理应用篇(十五)：GDAL 读取多光谱数据做归一化处理_空中旋转篮球的博客-CSDN博客基本思路，先读取csv数据，现读取每一列数据为一维数组，然后根据原始影像数据的大小（row，col）将一维数组塑形，即使用reshape函数变换为（row，col）形状对的二维数组。使用gdal将二维数组转换为tif图像。完整代码如下： 修改一下以上代码，可以将数据转化为和原始数据相同波段数的多波段遥感影像。可以使用for循环一个通道一个通道的
Python遥感图像处理应用篇(二十四)：Python绘制遥感图像各波段热力图（相关系数矩阵）
soderayer的博客
07-13 2284
给多光谱遥感图像各个波段绘制热力图，首先需要计算波段之间的相关系数矩阵，而计算遥感图像波段相关系数矩阵有不同的方法，常用的我们可以采用遥感图像处理软件计算，比如ENVI软件就可以计算相关系数矩阵，使用工具箱中的Statistics工具即可进行多种统计运算。我们这里不使用遥感图像软件计算，直接使用Python计算。要实现整个目标，我们首先需要提取遥感图像波段中各个像元的信息。这里我们可以将遥感图像读取并保存为csv文件，实现这一步我们可以参照之前的一篇文章，链接如下：Python遥感图像处理应用篇(十四)：G
Python空间数据处理1： GDAL读写遥感图像
vonuo
07-08 3万+
GDAL是空间数据处理的开源包，支持多种数据格式的读写。遥感图像是一种带大地坐标的栅格数据，遥感图像的栅格模型包含两部分的内容：栅格矩阵：由正方形或者矩形栅格点组成，每个栅格点所对应的数值为该点的像元值，在遥感图像中用于表示地物属性值；遥感图像有单波段与多波段，波段表示地物属性的种类，每个波段表示地物一种属性。大地坐标：空间数据参考表示地图的投影信息；仿射矩阵能将行列坐标映射到面坐标上。GDAL读写
解锁多光谱遥感 “密码”！从数据获取到 Python ，一文打通地质、农业、土壤应用全链路
最新发布
2401_89571885的博客
04-21 274
解锁多光谱遥感 “密码”！从数据获取到 Python ，一文打通地质、农业、土壤应用全链路
Python读取多种格式的高光谱数据
m0_64037204的博客
11-11 1040
python读取遥感影像数据
Python读取16位的多光谱遥感图像
欧阳罢笔的博客
09-15 3728
# 参考《python地理空间分析指南》 from osgeo import gdal_array import operator from functools import reduce import matplotlib.pyplot as plt import cv2 # 打开文件查看 src = 'phase2/20170218_UH_CASI_S4_NAD83.pix' # 50波段遥感数据 # 来自DFC2018 Houston # 下载https://hyperspectral.ee.uh
基于python多光谱遥感数据处理、图像分类、定量评估及机器学习方法
weixin_xiao5kou4chang6kai4的博客
12-28 2245
通过对光谱、图像等数据处理，掌握岩矿、土壤、植被等地物的光谱特征和图像特征，结合ENVI等专业软件、Python开发工具平台，开展多光谱数据预处理、图像分类、定量评估、机器学习等方法的实践和开发，提高运用多光谱遥感技术解决实际问题能力。
“Python+”集成技术高光谱遥感数据处理与机器学习深度应用
aishangyanxiu的博客
04-17 883
结合Python编程工具，专注于解决高光谱数据读取、数据预处理、高光谱数据机器学习等技术难题，通过复现高光谱数据处理和分析过程，并解析代码，使学员掌握python高光谱数据处理技巧。通过矿物识别、农业应用、木材含水量提取、土壤有机碳评估等案例，提供可借鉴的高光谱应用技术方案，结合Python科学计算、可视化、数据处理和机器学习库，深入讲解应用开发。深入探讨了高光谱遥感数据处理技术，涵盖了基本概念、成像原理、数据处理和分析方法，以及运用机器学习和深度学习模型提取和应用高光谱信息的技术。
GDAL读取遥感影像转化为数组，并储存为文件
qq_44407908的博客
10-15 1989
基于Python的GDAL库，将遥感影像的TIFF文件，转换为数组形式，并储存为文件。
python将tif（多光谱，高光谱）遥感图像转为mat格式数据
Mr_zhang66的博客
06-23 4803
import scipy.io as sio import skimage.io 读入tif图 imgpath = r'E:\高光谱数据\014140790010_01\014140790010_01_P001_MUL\19AUG26035324-M2AS_R1C1-014140790010_01_P001.tif' imggt = skimage.io.imread(imgpath) 转为mat sio.savemat(r"E:\高光谱数据\014140790010_01\014140790010_0
python多光谱遥感数据处理、图像分类、定量评估及机器学习方法应用
PhyliciaFelicia的博客
02-22 2582
python、多光谱遥感数据处理、图像分类、定量评估、机器学习、遥感、近红外、热红外、高光谱、作物估产、长势评估、ENVI、矿物识别、土壤质量评价、图像处理、岩矿、土壤、植被
Python遥感图像处理应用篇(十七)：GDAL 将归一化处理csv数据转化为多波段遥感影像
soderayer的博客
06-19 1025
本文所使用数据和上一篇数据相同，上一篇链接：Python遥感图像处理应用篇(十六)：GDAL 将归一化处理csv数据转化为遥感影像_空中旋转篮球的博客-CSDN博客上一篇文章将csv文件转化为单波段图像，这一篇将csv文件合成为一个多波段图像。注意一个前提是需要知道原图像的基本信息，如数据行列号，投影信息，数据类型信息等。实现写入多光谱数据的时候，核心问题在于使用循环写入各个通道数据。这里直接写入的是二维数组数据，看别人的文章的时候，有的使用的三维数组，如格式（channel,row,col）,我在使用的时
GDAL读取多通道图像并False Color可视化
weixin_52200494的博客
03-11 590
首先，确保你已经安装了GDAL库以及Python的GDAL绑定（如GDAL或rasterio库）。"""开发时间：2024/2/26 21:56@author:xd"""else:plt.show()如有错误的地方，还请指正！有疑惑的同学可以在评论区留言，无论能否解决，笔者看到一定回复。
Python分类文件（大疆精灵4多光谱版PM4影象分类）
weixin_45452300的博客
11-03 1390
任务描述 大疆精灵4多光谱版有六个通道，每次拍摄都会生成六张图片。图片的顺序为 0. 可见光 1. 蓝光 2. 绿光 3. 红光 4. 红边 5. 近红外。 在用其他软件进行拼接时需要对六个通道的照片分别进行拼接，而不同通道的图片仅是最后一位数字不同，因此要将相同通道的图片放到一起，易于利用软件拼接。 基于Python的os库和shutil库进行文件转移操作，代码如下`import shutil import os path = r’E:\博士期间记录\照片\站内\高光谱\20210811’ files
ORG/GDAL利用csv文件生成shapefile点文件
liuph_的博客
08-20 3549
输入csv文件，确保前三列为ID/经纬度，后面的列数及列的类型不限制，生成对应的点文件，默认设置地理坐标系为WGS84，投影坐标为通用墨卡托49N。        #include #include "ogrsf_frmts.h" #include "gdal.h" #include "gdal_priv.h" #include "cpl_string.h" #include #inclu
基于python的gdal读取遥感影像
a_aharon的博客
07-21 8337
空间数据文件访问——遥感影像的读写
使用GDAL库随机选取指定像元值的栅格经纬度并保存在CSV文件
qq_43827597的博客
03-05 668
writer.writerow(['Index', 'Longitude', 'Latitude']) # 写入表头。raise Exception("值为1的像元数量不足5000个")raise Exception("无法打开栅格文件")print(f"随机点已保存到 {output_csv}")# 将经纬度坐标保存到CSV文件。# 获取值为1的像元位置。# 随机选择500个像元。# 转换为经纬度坐标。
python环境下的遥感图像dct水印插入实现
梦醒时分
05-05 2090
毕设在做在Python环境下的遥感图像水印算法 最近看了挺多博客 下面是自己整理的 有瑕疵的话请大神指正 刚开始学 希望大家多交流我的水印用的是灰度图像 但是我处理以后只有0和255两个颜色 238*238我的遥感图像是6通道8431*7646import gdalimport numpyfrom skimage import data,img_as_float,img_as_ubyteimpo...
Python中reshape的用法及多个二维数组合并为三维数组的方法
热门推荐
Bebr的博客
01-11 3万+
reshape(shape) ： 不改变数组元素，返回一个shape形状的数组，原数组不变。是对每行元素进行处理 resize(shape) ： 与.reshape()功能一致，但修改原数组 In [1]: a = np.arange(20) #原数组不变 In [2]: a.reshape([4,5]) Out[2]: array([[ 0, 1, 2, 3, 4], ...
评论 10
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

空中旋转篮球

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值