Sentinel-2（哨兵-2）L1C数据辐亮度（辐射定标）和TOA反射率的获取说明

已于 2022-11-02 23:00:51 修改
阅读量4.9k 收藏 32
点赞数 2
分类专栏： 笔记 文章标签： sentinel python
于 2022-11-02 23:00:25 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41638841/article/details/127660752
版权

         哨兵-2 L1C级数据是我们常用的一种数据，但是如和从L1C数据中获取需要的大气顶表观反射率(TOA)以及大气顶的辐亮度（W/(sr\cdot m^{2}\cdot \mu m)）呢?

         获得反射率TOA一般就是用DN值除以10000（一万），这个值是固定的，在MTD_MSIL1C.xml文件里对应QUANTIFICATION_VALUE，可以看到这里是10000。

        反射率获取比较简单。辐亮度获取就麻烦了，找了半天没有明确公式。后来还是在envi网站上看到了，是描述envi哨兵2辐射定标插件。他这里给出了公式。

 

其中，\rho _{\lambda }：大气表观反射率（TOA）；L_{\lambda }:辐亮度，W/(sr\cdot m^{2}\cdot \mu m)；d：日地距离，注意这里似乎不是指绝对的米或千米，二是特指天文单位，太阳地球间的平均距离；ESUN_{\lambda}:太阳辐照度，(W/m^{2}\cdot \mu m);\theta:太阳高度角。

这个公式理论上是对的，参考辐射传输方程问题不大。接下来说说这几个值怎么获取。

  d:日地距离，打开MTD_MSIL1C.xml文件，在Product_Info节点中有个名称为U的节点，该节点数值即为日地距离。如0.980958599408787

辐射率已经说了，原始DN值除以10000.

太阳辐照度：MTD_MSIL1C.xml文件，Solar_Irradiance_List节点存了所有波段的辐照度。

太阳高度角：不在MTD_MSIL1C.xml文件里，在GRANULE/L1C_T51RTQ_A015101_20180514T024123下包含另一个名为MTD_TL的xml文件，打开这个xml文件，找到节点Mean_Sun_Angle，在此节点下包含两个角度，即ZENITH_ANGLE(天顶角)，AZIMUTH_ANGLE(方位角)。我们知道太阳高度角与天顶角互余，因此天顶角的余角即为太阳高度角。如xml中读取ZENITH_ANGLE值为19.3760450927414,那么太阳高度角为70.6239549072586,单位为度（这一段引自链接：https://blog.csdn.net/u013471015/article/details/88663209）

 

GEE|哨兵L1C和L2A求出来的NDVI到底差距多大
个人学习笔记
08-19 1953
这个问题我之前问过我的导师，早些时候没有SR数据，大家都用TOA数据，后来有了SR数据，很多研究就开始采用SR数据了，理论上SR更准确，具体来说大气校正前后Blue波段反射率有较大差异，其他波段也有差异，但是这种差异相对比较小，所以计算NDVI之类的指数，影响不大，但如果计算有Blue波段的指数，那就有影响了。对于植被的话，差别还是有的，毕竟NDVI的范围只有[-1,1]，差个十万八千里也不可能，但是经过这一对比实验，应该不用再纠结到底要不要选用大气校正后的数据了。还是用具体实验看看吧。
用ENVI软件对Sentinel-2（哨兵2）数据进行大气校正流程整理，亲测可用
02-12
用ENVI软件对Sentinel-2（哨兵2）数据进行大气校正流程整理，亲测可用用ENVI软件对Sentinel-2（哨兵2）数据进行大气校正流程整理，亲测可用用ENVI软件对Sentinel-2（哨兵2）数据进行大气校正流程整理，亲测可用
哨兵2号（Sentinel-2）介绍、下载、预处理及批处理
热门推荐
GeoSuper的博客
04-17 40万+
哨兵2号是高分辨率多光谱成像卫星，携带一枚多光谱成像仪(MSI)，用于陆地监测，可提供植被、土壤和水覆盖、内陆水路及海岸区域等图像，还可用于紧急救援服务。分为2A和2B两颗卫星。 第一颗卫星哨兵2号A于2015年6月23日01:52 UTC以“织女星”运载火箭发射升空。 6月29日，在轨运行4天的哨兵-2A卫星，传回了第一景数据，幅宽290km，卫星第一次扫描的范围是从瑞典开始，经过中欧和地中.........
一次讲清楚关于Sentinel-2 的处理级别信息，不在为预处理烦恼
qq_43805185的博客
02-19 864
大白话一次性讲清楚下载的Sentinel-2号的不同处理级别之间的差异，不要为预处理烦恼
关于哨兵2号数据辐射定标
weixin_44389786的博客
02-18 1874
最近在做sentinel2的TOA对比，下载到了L1C级别的数据但是DN值不是TOA值 查找了相关文档找到了一个计算公式 然后就去找相关参数 RADIO_ADD_OFFSET 和 QUANTIFICATION_VALUE 发现在下载数据的文件夹中 去找原数据的文档发现只有QUANTIFICATION_VALUE 所以会不会是直接用DN/QUANTIFICATION_VALUE，可能是这个数据中没有给定RADIO_ADD_OFFSET，所以就先这么做了，如果有大神知道可以指导一下。 ......
利用GEE提取哨兵-2卫星影像的多点反射率数据
qq_36686882的博客
07-17 724
在Google Earth Engine (GEE)平台上实现了对潍坊市2020年9月期间的哨兵-2卫星影像进行提取、处理和多点反射率数据的导出。主要功能包括定义多点坐标，去除云层干扰，提取指定波段的反射率数据并导出为CSV文件。
IDL学习——哨兵2 L1C数据辐射定标
weixin_42176976的博客
11-08 2457
对哨兵2L1C数据进行辐射定标，主要为大气表观反射率数据转换为辐射亮度数据。
利用Sen2cor对哨兵2号L1C级多光谱数据进行辐射定标和大气校正
05-29 1481
然后，在CMD中，切换到你数据所在的文件夹，如：cd G:\s2-gucheng（不熟悉CMD的，记得先使用G:切到G盘）；此时，你就可以在CMD中，根据Sen2cor语法要求（参见cmd给出的帮助信息），输入处理命令了，语法结构基本上是：L2A_Process+数据相对路径+可选参数，如：L2A_Process S2A_MSIL1C_20171005T030551_N0205_R075_T50SLJ_20171005T031513.SAFE --resolution=10 --refresh；
Sentinel-2数据辐射定标及大气校正
Corollary的博客
12-30 4750
使用工具： 1、SNAP平台 及 Sen2cor插件 处理过程： 1、下载SNAP平台和Sen2cor插件，建议安装最新版本 2、SNAP平台直接安装； 3、Sen2cor解压之后（文件夹名为：Sen2Cor-2.4.0-win64）将文件夹放置自己的用户文件夹下（比如：C:\Users\Administrator\AppData\Local） 双击运行其中的L2A_Process.bat文件。 4、打开windows CMD(Win键 + r键打开运行，输入CMD) 5、在CMD中定义到Sen2cor的路
使用Sen2cor对Sentinel-2辐射定标和大气校正
所有专家都曾始于新手
10-04 1252
其中地址要定位到.SAFE文件夹下，后面–resolution 10表示选择输出的分辨率，还有30m和60m可以选。对于一个文件夹下的多个文件，还有批处理的方法。用了17分钟，处理好了一景。先定位到插件所在文件夹下。
Sentinel-2 大气校正 地表反射率大于1 scale 0.0001
wsdxxwr的博客
05-08 6463
Sentinel-2 大气校正之后地表反射率大于1情况说明) Sentinel-2 DN 正常取值范围在0-10000之间，除非使用SNAP（显示的卫地表反射率），使用其他的软件如ENVI，显示的数值就是DN值，取值范围0-10000 或者超过。至于为什么超过了10000，具体原因就不多说了。在有云的地方，可以查看，相应的像元DN值有超过10000的现象。请记住超过10000的部分就不用要了，以免...
哨兵二号介绍 哥白尼：哨兵2号（Sentinel-2）
11-08
来自于哨兵二号应用分析介绍，传感器和平台分析课程资源学习应用。
GEE数据集——Harmonized Landsat Sentinel-2 (HLS) 卫星sentinel-2哨兵-2（HLS）
此星光明博客
12-26 1764
统一大地遥感卫星哨兵-2（HLS）项目通过虚拟卫星传感器群提供一致的地表反射率（SR）和大气层顶部亮度（TOA）数据。陆地成像仪（OLI）安装在美国宇航局/美国地质调查局的联合陆地卫星 8 号和陆地卫星 9 号上，而多光谱仪（MSI）则安装在欧洲的哥白尼哨兵-2A 号和哨兵-2B 号卫星上。通过综合测量，可以每 2 到 3 天以 30 米的空间分辨率对陆地进行全球观测。
GEE数据集：HLSL30（Landsat 8/9）和 HLSS30（Sentinel-2 A/B的比较（数量和大小）
最新发布
此星光明博客
04-29 897
统一大地遥感卫星哨兵-2（HLS）项目通过虚拟卫星传感器群提供一致的地表反射率（SR）和大气层顶部亮度（TOA）数据。陆地成像仪（OLI）安装在美国宇航局/美国地质调查局的联合陆地卫星 8 号和陆地卫星 9 号上，而多光谱仪（MSI）则安装在欧洲的哥白尼哨兵-2A 号和哨兵-2B 号卫星上。通过综合测量，可以每 2 到 3 天以 30 米的空间分辨率对陆地进行全球观测。
gee教程——sentinel-2数据的两种数据预（去云）处理方式QA60和SLC
此星光明博客
04-18 1764
哨兵二号卫星是欧洲空间局（ESA）发射的一颗地球观测卫星，主要用于监测地球表面的变化。在数据预处理阶段，去云是非常重要的一个步骤，因为云覆盖是影响遥感数据质量和可用性的主要因素之一。在这里，我们将具体介绍哨兵二号数据预处理阶段中的两种去云方式：QA60去云方式和SLC去云方式。QA60去云方式是一种基于光学遥感数据的去云方法。它是根据云覆盖的不同，将图像分为不同的云类别进行分类。在这种方法中，QA60是指一种云掩膜数据集，其中包含了60种云类别。
GF6 WFV卫星视角影响、与Sentinel2协同及精度分析
副组长的博客
05-03 4775
GF6 WFV卫星视角影响、与Sentinel2协同及精度分析 李国春 摘要：Sentinel2和GF6 WFV都是共享免费的重要数据源，定量化处理精度是人们广泛关心的问题。本文用同地区同日时间相差7分钟的两个Sentinel2和GF6 WFV数据集进行处理分析比较，认为使用RSD平台经本文流程处理后，两种数据相同（近）波段表现一致，可以协同交互混合使用。注意到并提出了GF6 WFV大视角扫描...
GEE 教程——分析特定区域的 Sentinel-2 图像数据，计算并导出中位数反射率
此星光明博客
11-24 169
这段代码用于在 Google Earth Engine 中处理 Sentinel-2 卫星图像，主要目标是计算特定区域的中位数反射率，并导出结果。以下是代码的主要功能和步骤：1. **设置地图和区域**：- 将地图设置为混合视图。- 定义感兴趣区域（AOI），可以选择区域 2 或区域 3。2. **设置时间参数**：- 定义年份（2018年）和起始、结束的日历天数（第91天到第258天）。- 设置云覆盖阈值（70%）。
S2A哨兵数据的波段合成、镶嵌、TOA（大气表观反射率）和裁剪的操作
yuanYW7556的博客
05-19 8014
声明：数据来自于网络。 如果看此博客看不懂具体的操作，可以看我发布的视频。 此操作的视频教程，我已发到今日头条上，可以关注我的今日头条进行学习。或者直接搜索我的今日头条号：遥感科学家小圆。 这是链接： https://www.ixigua.com/i6828473082307411975/ 这是两期的哨兵数据，都是2017年的，已经做过layer stacking处理。 波段合成的具体意义，可以观看我的视频教程学习，这里不详细介绍为什么做波段合成了。 可以看到共有6个波段，有第十一和第十二波段。 对两
哨兵2号数据处理
03-29
### Sentinel-2 数据处理方法及工具 #### 1. 数据获取 为了获得高质量的 Sentinel-2 卫星数据，推荐通过欧洲航天局（ESA）官方网站或其他授权平台进行下载[^1]。由于 USGS 已经停止提供哨兵系列数据的支持，因此需要转向 ESA 或其他合作机构来完成这一过程。 对于长时间序列分析的需求，可以参考知乎专栏文章中的详细介绍[^2]，其中提到 L1C 和 L2A 是两种主要的产品级别。L1C 提供的是经过几何校正但未大气校正的数据；而 L2A 则进一步包含了大气校正后的表面反射率产品，适合用于更精确的地表监测研究。 #### 2. 数据预处理 在实际应用之前通常要经历几个重要的预处理阶段： - **裁剪与拼接**：如果感兴趣区域跨越多个影像，则需先将这些图像拼接到一起并按边界条件裁切出目标范围内的部分。 - **辐射定标**：此操作会把原始DN值转换成物理单位如TOA反射率等，以便于后续定量计算使用。 - **大气校正**：采用专门软件比如 Sen2Cor 来执行该步骤，它能够有效去除由大气散射等因素引起的影响从而提高最终成果的质量[^3]。 #### 3. 地理编码转化 当涉及到特定地理位置时可能需要用到地理信息系统(GIS)，此时就需要知道如何正确表达位置信息。例如，在某些情况下可能会遇到名为“Military Grid Reference System (MGRS)” 的坐标体系描述方式——即所谓的 '10/S/DG' 参数形式表示某一点所在的具体格网编号[^4] 。可以通过在线服务快速实现从常规纬度经度到 MGRS 编码之间的相互转变工作。 #### 4. 使用MATLAB 处理 Sentinel-2 数据 针对科研人员而言，Matlab 软件提供了强大的矩阵运算能力和丰富的函数库支持复杂算法开发需求。下面给出一段简单的脚本示例展示读取 tiff 文件以及显示彩图的过程: ```matlab % 加载TIF文件 filename = 'path_to_your_file.tif'; img = imread(filename); % 显示彩色图片 figure; imshow(img); title('Sentinel-2 Image'); ``` 以上代码片段仅作为入门指南的一部分内容呈现出来，并不代表完整的解决方案。更多高级功能还需要深入学习官方文档或者查阅相关资料文献才能掌握全面的技术细节。 ---
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值