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

已于 2022-11-02 23:00:51 修改
阅读量3.7k 收藏 26
点赞数 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）

 

qq_41638841
关注
  • 2
    点赞
  • 踩
  • 26
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
通过ENVI新版获取指定经纬度的影像波段反射率
stardadada的博客
05-02 2090
通过ENVI新版获取指定经纬度的影像波段反射率
哨兵二号介绍 哥白尼：哨兵2号（Sentinel-2）
11-08
来自于哨兵二号应用分析介绍，传感器和平台分析课程资源学习应用。
IDL学习——哨兵2 L1C数据辐射定标
weixin_42176976的博客
11-08 1988
对哨兵2L1C数据进行辐射定标，主要为大气表观反射率数据转换为辐射亮度数据。
关于哨兵2号数据辐射定标
weixin_44389786的博客
02-18 1598
最近在做sentinel2的TOA对比，下载到了L1C级别的数据但是DN值不是TOA值 查找了相关文档找到了一个计算公式 然后就去找相关参数 RADIO_ADD_OFFSET 和 QUANTIFICATION_VALUE 发现在下载数据的文件夹中 去找原数据的文档发现只有QUANTIFICATION_VALUE 所以会不会是直接用DN/QUANTIFICATION_VALUE，可能是这个数据中没有给定RADIO_ADD_OFFSET，所以就先这么做了，如果有大神知道可以指导一下。 ......
Sentinel-2数据详解（L1C，L2A做长时间序列必看）--新版本处理基线
最新发布
m0_59054747的博客
03-28 3308
Sentinel-2数据使用注意事项，TOA和SR的转换方法，批量预处理脚本
辐射定标、大气校正
热门推荐
GeoSuper的博客
10-30 3万+
辐射校正 指在光学遥感数据获取过程中，产生的一切与辐射有关的误差的校正（包括辐射定标和大气校正）。 三者关系如图：大气校正的准备过程为辐射定标 表观反射率：表观反射率就是指大气层顶的反射率,辐射定标的结果之一，大气层顶表观反射率，简称表观反射率，又称视反射率。英文表示为：apparent reflectance 辐射定标：将记录的原始DN值转换为大气外层表面反射率（或称为辐射亮度值）。...
关于哨兵2-LIC级数据辐射定标几个参数的获取
tyut_qc的博客
03-19 7539
Sentinel-2L1C级数据是经正射校正和亚像元级几何精校正后的大气表观反射率（TOA）产品，在我们需要使用辐射率（Radiance）数据时需要一些转换，所以该数据的辐射定标不在是简单的乘以一个Gain，加上一个Offset那么简单。而是利用公式关系反算辐射率数据。其详细计算公式参见(ENVI官博)http://blog.sina.com.cn/s/blog_764b1e9d0102y...
用ENVI软件对Sentinel-2（哨兵2）数据进行大气校正流程整理，亲测可用
02-12
用ENVI软件对Sentinel-2（哨兵2）数据进行大气校正流程整理，亲测可用用ENVI软件对Sentinel-2（哨兵2）数据进行大气校正流程整理，亲测可用用ENVI软件对Sentinel-2（哨兵2）数据进行大气校正流程整理，亲测可用
Sentinel-2数据辐射定标及大气校正
Corollary的博客
12-30 3791
使用工具： 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的路
Sentinel-2 大气校正 地表反射率大于1 scale 0.0001
wsdxxwr的博客
05-08 5549
Sentinel-2 大气校正之后地表反射率大于1情况说明) Sentinel-2 DN 正常取值范围在0-10000之间，除非使用SNAP（显示的卫地表反射率），使用其他的软件如ENVI，显示的数值就是DN值，取值范围0-10000 或者超过。至于为什么超过了10000，具体原因就不多说了。在有云的地方，可以查看，相应的像元DN值有超过10000的现象。请记住超过10000的部分就不用要了，以免...
Sentinel-2哨兵二号数据下载及处理教程.pdf
06-10
Sentinel-2哨兵二号数据下载及处理教程.pdf Sentinel-2 是高分辨率多光谱成像卫星，一颗卫星的重访周期为10天，两颗互补，重访周期为5天。分为2A和2B两颗卫星。2A于2015年6月23日01:52 UTC以“织女星”运载火箭发射...
Sentinel-2 MSI辐射表征和与Landsat-8 OLI的交叉校准
06-04
采集了撒哈拉沙漠中空间均匀区域分辨率为10、20和30 m的多光谱图像，用于直接比较MSI和OLI大气顶层（TOA）反射率。 本文介绍了Sentinel-2 MSI和Landsat 8 OLI传感器的8个对应光谱带的初始辐射交叉校准。 以经过良好...
新版哨兵数据（Sentinel-2 MSI Sentinel-3 OLCI）批量下载代码
10-16
之前以往用来下载哨兵数据的哥白尼开放...本代码能从新版哥白尼数据空间生态系统批量下载哨兵数据（Sentinel-2 MSI Sentinel-3 OLCI），卫星宫中号——“海研人” 后台回复“哨兵”直接领，别在CSDN上下。拉吉玩意。
Sentinel-2遥感卫星全球矢量数据（轨道、瓦片与覆盖范围）
03-05
该数据为本人自欧洲航天局（ESA）官网下载的Sentinel-2遥感卫星相关矢量数据（shp格式），主要包括Sentinel-2卫星的轨道覆盖范围，全球Grid瓦片与编号数据、可供下载的影像覆盖区域等。
sentinel2-cloud-detector:适用于Python的Sentinel-2图像的Sentinel Hub Cloud Detector
05-07
Sentinel Hub的Sentinel-2影像云探测器 注意：s2cloudless掩码现在可以作为Sentinel Hub中的预计算层使用。 查看 和。 s2cloudless Python软件包在Sentinel-2影像中提供了自动的云检测功能。 该分类基于Sentinel ...
snap软件中哨兵2A数据预处理及六种常用植被指数的计算
weixin_47649329的博客
04-19 1万+
1.使用snap软件对哨兵二号数据进行预处理。 2.使用envi处理哨兵二号数据，进行波段组合。 3.学习不同植被指数的计算公式，进行六种植被指数（归一化植被指数ndvi、增强型植被指数EVI、差值植被指数DVI、比值植被指数RVI的计算、土壤调整型植被指数MSAVI、土壤调整型植被指数(SAVI)）的计算。 二、实验数据： 天水地区的哨兵2A数据：S2A_MSIL2A_20201031T034911_N0214_R104_T48SVC_20201031T071812 三、实验步骤： 1、(1)
哨兵L1C数据从下载到数据处理
dayuhaitang1的博客
12-28 8830
因为最近要使用哨兵数据，所以花费了好几天在网上学习怎样获取哨兵数据，终于学会了怎样下载和处理哨兵数据，在这里就特意的做了一篇笔记来记录一下整个过程。 一、哨兵数据的下载 1、欧空局网址：https://scihub.copernicus.eu/dhus/#/home 二、Sen2Cor软件下载 ...
哨兵2介绍
rs_gis的博客
09-22 7356
1、哨兵2产品介绍 哨兵2共两颗卫星：2A和2B 发射时间： 2A：2015年6月23日 2B：2017年3月7日 周期： 单星10天。两组卫星同时运行，周期为5天 产品级别： Level-0： 原始数据。 Level-1A：包含元信息的几何粗校正产品。 Level-1B：辐射率产品，嵌入经GCP优化的几何模型但未进行相应的几何校正。 Level-1C：经正射校正和亚像元级几何精校正后的大气表观反射率产品。L1C级多光谱数据（MSI），Sentinel-2 L1C是经过几何精校正的大气表观反射率产品，是经过
以sentinel-2影像为例,说明时间分辨率、光谱分辨率、辐射分辨率、空间分辨率
05-30
Sentinel-2是一颗欧洲空间局（ESA）发射的卫星，搭载了多光谱成像仪（MSI）传感器，可以提供高分辨率、多光谱遥感影像。以下是Sentinel-2影像的各种分辨率的解释： - 时间分辨率：Sentinel-2的时间分辨率为5天。也就是说，每隔5天，Sentinel-2卫星会在地球上同一地点拍摄一张影像，这样可以捕捉到地球表面的变化和演变情况。 - 光谱分辨率：Sentinel-2的MSI传感器可以提供13个波段的数据。其中，4个波段的分辨率为10米，6个波段的分辨率为20米，3个波段的分辨率为60米。这些波段可以提供不同的信息，例如植被、水体、土地利用等等。 - 辐射分辨率：Sentinel-2的MSI传感器的辐射分辨率是12位，也就是说，每个像素点可以提供0到4095的数字化值，代表了接收到的电磁辐射强度。 - 空间分辨率：Sentinel-2的空间分辨率是10米、20米和60米。这些分辨率代表了卫星在地球表面上每个像素点所覆盖的实际面积大小。空间分辨率越高，代表着像素点所覆盖的面积越小，能够提供更加详细的地面信息。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值