使用GDAL读取Sentinel数据

最新推荐文章于 2023-03-07 22:23:11 发布
VIP文章 最新推荐文章于 2023-03-07 22:23:11 发布
阅读量4.1k 收藏 26
点赞数 1
分类专栏： GIS Python
本文为博主原创文章，转载请注明原文出处！
本文链接：https://blog.csdn.net/theonegis/article/details/82194108
版权

使用GDAL读取Sentinel数据

GDAL 2.1已经原生支持对于Sentinel数据的读取，我这里使用Sentinel-2光学卫星数据给出使用GDAL工具对其进行读取的方法。

这里我们要大概知道Sentinel数据的组织。下载下来的Sentinel数据是一个ZIP压缩包，里面包含了JPEG2000格式的影像数据以及一些XML格式的元数据文件。

GDAL将Sentinel数据看做一个数据集（概念上类似HDF格式的数据集），里面包含了很多子数据文件。所以，对于Sentinel数据的读取就和对于HDF数据的读取是相同的啦。

对于HDF或者NetCDF格式数据的读取参考我的博文：读取HDF或者NetCDF格式的栅格数据

使用GDAL命令行读取Sentinel数据的元数据信息

直接使用gdalinfo [文件名]可以查看Sentinel文件的元信息，如下图所示：

使用GDAL命令行读取Sentinel数据

从上面的图中我们可以看到所有的Subdatasets的文件全名，这样我们可以继续使用gdalinfo [子数据集全路径]的方式查看具体的子数据集的元数据信息

下图显示的数据子集中包含四个波段的

最低0.47元/天 解锁文章
卡尔曼和玻尔兹曼谁曼
关注
  • 1
    点赞
  • 踩
  • 26
    收藏
    觉得还不错? 一键收藏
  • 15
    评论
专栏目录
GDAL读取矢量数据的点坐标集.txt
10-10
基于C++的GDAL/OGR读取矢量数据（.shp）的点线面的坐标集，并且将坐标写入txt文档中。
GDAL\OGR读取数据示例 C#版本
01-09
GDAL\OGR读取数据示例 C#版本 有读取栅格和矢量，并显示数据里面的基本信息，可以用来作为入门例子 GDAL\OGR读取数据示例 C#版本 有读取栅格和矢量，并显示数据里面的基本信息，可以用来作为入门例子
c++GDAL读取显示shp数据
02-15
使用C++并利用GDAL库实现对GIS中shp数据的读取和显示。该版本没有放缩功能，有简陋版放缩功能的版本，如有需要可以下载后评论可获取。
VS2013 C++ 配置GDAL 提供对Sentienl-2 的驱动支持
遥感基础算法爱好者
07-06 939
1. gdal下载下载路径是：https://trac.osgeo.org/gdal/wiki/DownloadSource2. 下载Jasper（支持JPEG2000 格式，即Sentinel-2的数据格式）下载路径是：http://download.osgeo.org/gdal/jasper-1.900.1.uuid.tar.gz编译jasper，注意勾掉jiv字工程，因为它不是必须，而且...
GDAL学习1:初识GDAL以及读取处理保存影像信息
君子之学必日新，日新者日进也。不日新者必日退，未有不进而不退者
04-22 4831
GDAL学习1:初识GDAL以及读取影像信息前言一、GDAL是什么？栅格数据组织二、使用C++ GDAL3读取处理保存栅格数据1.引入库2.读取影像及信息函数3.主函数三、使用python GDAL读取处理保存栅格数据1.引入包2.读取影像及信息总结 前言 随着遥感、GIS的不断发展，GDAL这个库也越来越重要，很多人都开启了学习GDAL之旅，本文就介绍了GDAL的基础知识和利用GDAL读取和保存影像信息的c++和python代码。 一、GDAL是什么？ GDAL(Geospatial Data Abs
SentinelHASP读狗工具
04-30
SentinelHASP读狗工具，Sentinel HASP SRM读狗工具，目前最新版本
用sentinel-2光学影像通过GAN生成Sentinel-1影像
个人学习笔记
03-07 237
参考的代码是：https://github.com/ZirvedaAytimur/GAN_Electro-Optic_to_SAR_Images。记录在学习用sentinel-2光学影像通过GAN生成Sentinel-1影像时遇到的问题。
python+gdal+excel构建矢量图
x19979的博客
04-25 604
本文介绍一个从excel中读取数据，从而构建点状矢量文件的方法（当然构建面状，线状也是类似的，如有需要可以留言），使用的python+gdal库+xlrd。当然本程序不一定是最简化的，各位读者可以自行优化。本文并不是附上所有的程序而是针对重要的代码段进行讲解。 first： import xlrd ###读取excel文件，比较简单，读取出每个点的 ##信息比如经纬度，还有其他的属性。 file...
加密狗型号查看读取工具
02-28
usb通用加密狗型号查看工具，读狗插件，部分 加密狗/加密锁 VID识别加密狗具体型号USB\\VID_096E&PID;_0001 飞天诚信(ftsafe) 飞天4有驱型 加密锁 ROCKEY 4 Rockey4 R4HID\\VID_096E&PID;_0006 飞天诚信(ftsafe) 飞天4无驱型 加密锁 rockey 4 nd ROCKEY4ND R4ndUSB\\VID_096E&PID;_0201 飞天诚信(ftsafe) 飞天2无驱型 加密锁 rockey 2 ROCKEY2 R2USB\\VID_096E&PID;_0304 坚石诚信(Rockey) ET199 Token ET199USB\\VID_096E&PID;_0303 坚石诚信(Rockey) ET99 Token ET99USB\\VID_096E&PID;_0201 飞天诚信(FTSafe) Rockey 2 加密狗,无驱型USB\\VID_096E&PID;_0202 飞天诚信(FTSafe) Rockey 3 智能狗,无驱型USB\\VID_096E&PID;_0205 飞天诚信(FTSafe) Rockey 1 智能狗,无驱型USB\\VID_096E&PID;_0207 飞天诚信(FTSafe) Rockey 1 智慧型,无驱型USB\\VID_096E&PID;_0209 坚石诚信(FTSafe) 　　 ARM 时钟智能狗USB\\VID_096E&PID;_0403 飞天诚信(ftsafe) 飞天6无驱型 加密锁 rockey6 smart R6USB\\VID_83D3&PID;_3773 域天(域之天) 域之天 域天 专业型/经济易用型 加密锁 USB\\VID_3689&PID;_8762 域天(域之天) 域之天 域天 密码型 加密锁 USB\\VID_1199&PID;_8282 域天(域之天) 域之天 域天 简单型 加密锁 USB\\VID_1BC0&PID;_1001 深思洛克(SenseLock) 灵锐 I(Living) 加密狗,无驱型USB\\VID_1BC0&PID;_8101 深思洛克(SenseLock) 灵锐 定制加密狗,无驱型USB\\VID_04B4&PID;_4A58 深思洛克(SenseLock) 深思 3(sense3)加密狗,无驱型USB\\VID_04B4&PID;_4A59 深思洛克(SenseLock) 深思 3(sense3)加密狗,有驱型USB\\VID_1BC0&PID;_8013 深思洛克(SenseLock) 精锐 E(Elite-E)智能狗,无驱型USB\\VID_0471&PID;_485D 深思洛克(senseLock) SenseIV2.x 深思洛克 sense4 senseIV 深思4 有驱USB\\VID_0471&PID;_485E 深思洛克(senseLock) SenseIV2.x 深思洛克 sense4 senseIV 深思4 无驱USB\\VID_1312&PID;_2012 域天(域之天)域之天 域天 32位智能型 易用型智能锁 加密锁 无驱USB\\VID_1312&PID;_1912 域天(域之天)域之天 域天 32位智能型(有驱型)USB\\ 　　VID_1312&PID;_2012 域天(域之天)域天(域之天)简单智能型USB\\VID_3689&PID;_3689 EncryptPE 订制型(域天)EPE OEM加密狗 域天密码型 EncryptPE订制型 USB\\VID_5042&PID;_FFFF 安全芯 schip 安全芯UKF-002USB\\VID_08E2&PID;_0002 赛孚耐(SafeNet) safenet microdog UMI 彩虹加密狗 微狗USB\\VID_08E2&PID;_0004 赛孚耐(SafeNet) safenet softdog UDA 彩虹加密狗 软件狗USB\\VID_08E2&PID;_0005 赛孚耐(SafeNet) 精灵狗UGA(GeniusDog),有驱型USB\\VID_08E2&PID;_0006 赛孚耐(SafeNet) safenet GrandDog 宏狗UGRAUSB\\VID_08E2&PID;_0008 赛孚耐(SafeNet) safenet SuperDog 超级狗USB\\VID_04B9&PID;_0300.DeviceDesc="SafeNet USB SuperPro/UltraPro" USB\\VID_04B9&PID;_0301.DeviceDesc="SafeNet USB Security Security Token" USB\\VID_04B9&PID;_0302.DeviceDesc="SafeNet USB Security Security Token" USB\\VID_04B9&PID;_0303.DeviceDesc="SafeNet USB Security Security Token& 　　quot; USB\\\\VID_04B9&PID;_0304.DeviceDesc="SafeNet USB Security Security Token" USB\\VID_04B9&PID;_8000.DeviceDesc="SafeNet Sentinel Hardware Key" USB\\VID_04B9&PID;_8001.DeviceDesc="SafeNet Sentinel Hardware Key RTC" USB\\VID_04B9&PID;_8002.DeviceDesc="SafeNet Sentinel Hardware Key Mega-Memory" USB\\VID_04B9&PID;_8003.DeviceDesc="SafeNet Sentinel Hardware Key Mega-Memory RTC" USB\\VID_04B9&PID;_8004.DeviceDesc="SafeNet Sentinel Dual Hardware Key" USB\\VID_04B9&PID;_8005.DeviceDesc="SafeNet Sentinel Dual Hardware Key RTC"USB\\VID_0925&PID;_A508 北京磐石诚信NT系列加密狗USB\\VID_7F5E&PID;_DE5D 广州飞盾NT系列加密锁USB\\VID_5042&PID;_???? 江波龙(NetCom) 安全芯(SChip) 加密狗,有驱型USB\\VID_5149&PID;_???? 江波龙(NetCom) 易锁(SimpleKey)加密狗,无驱型 　　USB\\VID_0925&PID;_7816 世纪龙脉(Century Longmai) NOX 时钟狗,无驱型USB\\VID_0925&PID;_9627 世纪龙脉(Century Longmai) NOX+ 单机狗,无驱型USB\\VID_0925&PID;_7A13 世纪龙脉(Century Longmai) DAM2+网络狗,无驱型USB\\VID_0925&PID;_5906 世纪龙脉(Century Longmai) NOX 定制狗,无驱型USB\\VID_0925&PID;_6B05 世纪龙脉(Century Longmai) NOX 加密狗,无驱型USB\\VID_0925&PID;_7A21 世纪龙脉(Century Longmai) NOX2 加密狗,无驱型USB\\VID_0925&PID;_7122 世纪龙脉(Century Longmai) DAM2+加密狗,无驱型USB\\VID_0925&PID;_5624 世纪龙脉(Century Longmai) DAM2 加密狗,无驱型USB\\VID_0925&PID;_8705 世纪龙脉(Century Longmai) NOX5 加密狗 智能狗,无驱型USB\\VID_0925&PID;_7504 世纪龙脉(Century Longmai) NOX5 主锁USB\\VID_0529&PID;_0001 阿拉丁(Aladdin) HASP HL Time 有驱型USB\\VID_0529&PID;_0001 阿拉丁(Aladdin) HandLock 有驱型USB\\VID_0830&PID;_2002 和升达 SK-USB180 无驱型USB\\VID_5149&PID;_EACE 杭州唯赛 易锁(Simple Key) 无驱型USB\\VID_ 　　6A75&PID;_9801 上海十条电子 JSP-C USB Key 有驱型USB\\VID_10F7&PID;_8469 MegaLock 无驱型USB\\VID_064F&PID;_0BD8 WIBU WIBU RU+加密狗,有驱型
python引入gdal包读取dem数据
08-07
利用python引入GDAL包，读取栅格tif数据信息
Sentinel2-L2A数据的下载、读取及基础操作
sievert的博客
06-18 1万+
哨兵2号是高分辨率多光谱成像卫星，携带一枚多光谱成像仪(MSI)，用于陆地监测，可提供植被、土壤和水覆盖、内陆水路及海岸区域等图像，还可用于紧急救援服务。分为2A和2B两颗卫星。 第一颗卫星哨兵2号A于2015年6月23日01:52 UTC以“织女星”运载火箭发射升空。 第二颗卫星哨兵2号B于2017年3月07日北京时间9时49分 UTC以“织女星”运载火箭发射升空。 欧洲航天局介绍说，“哨兵-2B”卫星与2015年6月发射的“哨兵-2A”卫星为同一组，携带高分辨率多光谱成像装置，主要用于监测土地环境，可提
gdal学习笔记2-数据读写
qq_42031694的博客
10-30 2124
今天常识读写sentinel-2的数据。 先写上gdal支持的数据类型 GDT_Unknown : 未知数据类型 GDT_Byte : 8bit正整型 (C++中对应unsigned char) GDT_UInt16 : 16bit正整型 (C++中对应 unsigned short) GDT_Int16 : 16bit整型 (C++中对应 short 或 short int) GDT_UInt3...
Google Earth Engine（GEE）——sentinel2数据介绍
qq_31988139的博客
11-25 5835
之前有粉丝在公众号问我，对于sentinel2是否需要进行大气校正的问题，这里就简单写一下关于这个数据的问题： 现在回答是因为所有光学影像都是需要去云的，因为遇到没有云的影像概率太小了，况且数据预处理阶段的去云是必不可少的一个关键环节，一般在GEE平台需要的处理就是去云处理，而很多modis数据和Landsat8数据，以及sentinel数据都是通过QA波段的位掩技术实现的。所以一般无论做反演还是监测都需要进行这一步修正。 现在先给大家推荐一个官方的文档介绍的链接： Sentinel-2 User H
04-GDAL读取Sentinel 2数据
suntongxue100的博客
03-01 684
Sentinel2 数据在GDAL中读取类似于科学数据集的获取方式，如图： 如常用的蓝、绿、红、近红外常用波段，就存储于第一个数据集list[0][0],list[0][1]为描述信息。 GDAL中读取数据示例如下： from osgeo import gdal import numpy as np import os os.environ['CPL_ZIP_ENCODING'] = 'UTF-8' filename = "G:\mypython3\gdal\data\S2A_MSIL1C_2017
python gdal 无法打开 Sentinel-2 L2A 数据的问题
遥感基础算法爱好者
04-22 999
在Sentinel-2 L2A 地表反射率产品过程中发现，GDAL不支持 读取 L2A 产品。这是由于GDAL版本较低的原因。可以安装2.4.0+版本。在gdal安装过程中发现，采用pip install gdal 无法成功安装gdal的情况，这个时候可以使用whl格式文件安装。链接地址如下：https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal，在这里...
Numpy中找出array中最大值所对应的行和列
热门推荐
给永远比拿愉快
03-07 7万+
Python特别灵活，肯定方法不止一种，这里介绍一种我觉得比较简单的方法。 如下图，使用x == np.max(x) 获得一个掩模矩阵，然后使用where方法即可返回最大值对应的行和列。 where返回一个长度为2的元组，第一个元素保存的是行号，第二个元素保存的是列号。
Numpy中如何给矩阵增加一行或一列
给永远比拿愉快
03-06 7万+
使用Python的numpy的array结构，如何给矩阵增加一行或者一列呢？ 下面提供一种方法，当然numpy还提供了很多API函数可供选择。
Python中一个for循环循环多个变量
给永远比拿愉快
10-25 5万+
首先，熟悉一个函数zip，如下是使用help(zip)对zip的解释。 Help on built-in function zip in module __builtin__: zip(...)     zip(seq1 [, seq2 [...]]) -> [(seq1[0], seq2[0] ...), (...)]     Return a list of tuple
Python中如何优雅地使用switch语句
给永远比拿愉快
03-08 3万+
文章目录Python中如何优雅地使用switch语句案例一（简单情况）案例二（带条件判断） 版权声明：本文为博主原创文章，转载请注明原文出处！ 写作时间：2019-03-07 13:49:45 Python中如何优雅地使用switch语句 我们知道Python中没有类似C++或者Java中的switch...case语句，我们可以使用多个if...elif...else进行模拟，但是这样的写法让代...
gdal获取栅格数据范围
最新发布
05-14
使用GDAL库获取栅格数据范围可以使用以下代码： ```python from osgeo import gdal # 打开栅格数据集 ds = gdal.Open('path/to/raster.tif') # 获取栅格数据集的范围 xmin, xres, _, ymax, _, yres = ds.GetGeoTransform() xmax = xmin + (ds.RasterXSize * xres) ymin = ymax + (ds.RasterYSize * yres) ``` 其中，`GetGeoTransform()`方法返回一个元组，包含了数据集左上角的x坐标、地理坐标系下的x和y方向上的像素分辨率、数据集左上角的y坐标、地理坐标系下的y方向上的像素分辨率和旋转角度。通过这些参数可以计算出栅格数据集的范围。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 15
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

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

余额充值