哨兵二号（Sentinel-2）卫星数据处理

已于 2024-04-23 10:58:42 修改
阅读量1.1w 收藏 122
点赞数 20
分类专栏： 科研 文章标签： arcpy
于 2022-04-11 15:42:56 首次发布
未经允许禁止转载
本文链接：https://blog.csdn.net/qq_42375636/article/details/122470730
版权

Github项目

特别提醒：本文所有操作所需时间都很长，因此在自己电脑上操作的朋友记得将睡眠模式关闭。

一、数据下载

1. 打开网址：

[USGS](EarthExplorer (usgs.gov)) ；

注：USGS 上的哨兵系列数据已经下架，建议访问欧空局获取数据，后续的数据处理流程是一样的。

2. 选择数据集：

1

3. 选择地区：

2

  • 点击最下方红框中的陕西省即可选择；
  • 该步骤如果搜索不到国内的省份，要么是省份的英文写错了，要么是需要科学上网。

4. 选择日期和云度：

3
4

5. 选择其他参数（可选）：

5

6. 生成结果：

6

7. 选择图像：

7

  • 对每一页都做相同操作，可以看到图中的纸箱子 icon 变成了绿色，表示选定；

8. 生成下载链接：

8

9. 安装下载器：

9

10
11

  • 要安装该软件，还需要安装 java 10 64bit ，安装方法请自行百度；
  • 配好 java 的环境变量后，就可以安装该软件了；

10. 提交下载链接：

12
13
14
15
16

11. 通过软件下载：

  • 打开 bda 软件；
  • 登录账号；
  • 在打开的 Open Order 窗口选择你的下载申请，若该窗口关闭，可在 File 菜单栏打开；
  • 选择你的目标 order ；
  • 点击 Begin Download 进行下载；
  • 所下载的文件在 bda 软件的安装目录下；
  • 一般晚上 1 点下载速度较快且连接稳定；
  • 如果下载失败，出现了 error ，不要点击 Redownload ，该操作会删除下载失败的文件。此时应先下载完其他文件，最后重启 bda ，Select 错误的部分点击 Begin Download 即可继续下载；
  • 当 error 发生后还有一种可行的方法就是不去管它，有时候 bda 软件会自动重新接续下载；
  • 重启 bda 不成功，可能是因为它在后台运行，需要先打开资源管理器关闭它的进程。

二、数据介绍

1. 文件介绍

以 S2AMSIL1C20190122T025021N0208_R132T50RNV_20190122T065329.zip 为例：

  1. S2A：表示 Sentinel-2A 卫星；
  2. MSL：表示多光谱数据；
  3. L1C：表示产品等级为 L1C 级别，该级别只经过了几何校正，L2A 级别产品还经过了辐射校正，但该级别需要用户自己处理；
  4. 20190122T025021：表示数据获取时间，即 2019 年 1 月 22 日凌晨 2 点 50 分 21 秒，该时间为格林威治时间，比北京时间提前 8 个小时；
  5. N0208_R132：处理基线编号与相对轨道编号；
  6. T50RNV：拼接域编号；
  7. 20190122T065329：估计是产品生成的时间；

2. 卫星介绍

S2AS2AS2BS2B
Band NumberBand nameCentral wavelength (nm)Bandwidth (nm)Central wavelength (nm)Bandwidth (nm)Spatial resolution (m)
1Coastal aerosol443.927442.34560
2Blue496.698492.19810
3Green560.0455594610
4Red664.5386653910
5Vegetation Red Edge703.919703.82020
6Vegetation Red Edge740.218739.11820
7Vegetation Red Edge782.528779.72820
8NIR835.114583313310
8ANarrow NIR864.8338643220
9Water vapour945.026943.22760
10SWIR – Cirrus1373.5751376.97660
11SWIR1613.71431610.414120
12SWIR2202.42422185.723820

三、数据显示与处理

1. 软件下载：

2. 软件安装：

  • SNAP 无脑下一步，需要 Python 环境的用户请注意软件支持的版本；
  • Sen2Cor 加压即可，如果想要在任意文件夹下使用该工具，需要将其添加入系统环境变量；

17

3. 大气校正：

  • 将所有数据文件进行解压，放入一个单独的文件夹 Sentinel-2 中；
  • 打开 cmd ；
  • 进入 Sen2Cor 文件夹（如果添加了环境变量可跳过该步骤）；
  • 输入命令：for /D %s in (D:\Documentation\Project\Sentinel-2\S2A_MSIL1C*) do L2A_process --resolution=10 %s ，其中 --resolution=10 表示处理分辨率为 10 的图像；D:\Documentation\Project\Sentinel-2\S2A_MSIL1C* 表示所下载的数据集的位置；
  • 如果需要处理 S2B 星的数据，仅需要将上述操作中的 S2A_MSIL1C* 改为 S2B_MSIL1C* 即可。
  • 等待系统处理即可；

18

4. 导出 TIFF ：

  • 创建一个 Python 环境；

  • 下载 gdal (注意版本须与 python 版本匹配);

  • 将其放入你的 python 环境中的 Scripts 文件夹中（推荐使用虚拟环境）；

  • 在 Terminal 中激活你的 python （通常是在同一目录下执行 activate.bat 文件）；

  • 在 Terminal 中输入：pip install GDAL-xxxxxx ，注意是你自己的 GDAL 版本；

  • 按需配置其他环境（例如 numpy 等）；

  • 使用下述代码进行批处理：

    from osgeo import gdal
import os
import numpy as np
from osgeo import gdal, osr, ogr
import glob
# os.environ['CPL_ZIP_ENCODING'] = 'UTF-8'

def S2tif(filename):
    # 打开栅格数据集
    print(filename)
    root_ds = gdal.Open(filename)
    # print(type(root_ds))
    # 返回结果是一个list，list中的每个元素是一个tuple，每个tuple中包含了对数据集的路径，元数据等的描述信息
    # tuple中的第一个元素描述的是数据子集的全路径
    ds_list = root_ds.GetSubDatasets()  # 获取子数据集。该数据以数据集形式存储且以子数据集形式组织
    visual_ds = gdal.Open(ds_list[0][0])  # 打开第1个数据子集的路径。ds_list有4个子集，内部前段是路径，后段是数据信息
    # print(visual_ds)
    print(f'打开数据为：{ds_list[0][1]}')
    # print(f'投影信息：{visual_ds.GetProjection()}')
    # print(f'栅格波段数：{visual_ds.RasterCount}')
    # print(f'栅格列数（宽度）：{visual_ds.RasterXSize}')
    # print(f'栅格行数（高度）：{visual_ds.RasterYSize}')
    visual_arr = visual_ds.ReadAsArray()  # 将数据集中的数据读取为ndarray

    # 创建.tif文件
    band_count = visual_ds.RasterCount  # 波段数
    xsize = visual_ds.RasterXSize
    ysize = visual_ds.RasterYSize
    out_tif_name = filename.split(".SAFE")[0] + ".tif"
    driver = gdal.GetDriverByName("GTiff")
    out_tif = driver.Create(out_tif_name, xsize, ysize, band_count, gdal.GDT_Float32)
    out_tif.SetProjection(visual_ds.GetProjection())  # 设置投影坐标
    out_tif.SetGeoTransform(visual_ds.GetGeoTransform())

    for index, band in enumerate(visual_arr):
        band = np.array([band])
        for i in range(len(band[:])):
            # 数据写出
            out_tif.GetRasterBand(index + 1).WriteArray(band[i])  # 将每个波段的数据写入内存，此时没有写入硬盘
    out_tif.FlushCache()  # 最终将数据写入硬盘
    out_tif = None  # 注意必须关闭tif文件

if __name__ == "__main__":
    from osgeo import gdal
    SAFE_Path = (r'D:\Documentation\Project\Grassland ecology\Sentinel2_pretreatment\data')
    data_list = glob.glob(SAFE_Path + "\\*.SAFE")

    for i in range(len(data_list)):
        data_path = data_list[i]
        filename = data_path + "\\MTD_MSIL2A.xml"
        S2tif(filename)
        print(data_path + "-----转tif成功")
    print("----转换结束----")

  • 代码中第 17 行是取第一个数据子集（band 2、3、4、8），要使用其他 band 须自行更改（第 18 行代码可以查看每个数据子集的信息）；

  • 由于读取子数据集时，gdal 工具会将波段进行重组，导致波段的顺序与原数据集（sentinel-2 原始数据集）的波段顺序不同，因此可以通过 gdalinfo [文件名或子数据集名] 来获取重组后波段的顺序结构，方法如下：

    • 首先，打开并进入你所下载的 gdal 文件夹下的 Lib\site-packages\osgeo 目录；

    • 该目录下检查是否有一个 gdalinfo.exe 程序；

    • 在该目录下打开 cmd 命令行；

    • 输入 gdalinfo xxx/xxx.xml 查看数据信息，下面给出一个具体的例子以供参考：

      $ gdalinfo S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml
# 上面这行是我的输入，其中数据文件夹为 S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE
# 由于 gdalinfo 不能直接读取 S2A 的文件，因此应该读取 xml 文件：/MTD_MSIL2A.xml
# 程序支持写入相对或绝对路径，此处我是将数据文件放入了该目录下进行读取
# 下面是文件信息的输出：
Driver: SENTINEL2/Sentinel 2
Files: S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml
Size is 512, 512
Metadata:
  AOT_QUANTIFICATION_VALUE=1000.0
  AOT_QUANTIFICATION_VALUE_UNIT=none
  AOT_RETRIEVAL_ACCURACY=0.0
  BOA_QUANTIFICATION_VALUE=10000
  BOA_QUANTIFICATION_VALUE_UNIT=none
  CLOUD_COVERAGE_ASSESSMENT=1.049327
  CLOUD_SHADOW_PERCENTAGE=0.440092
  DARK_FEATURES_PERCENTAGE=30.990644
  DATATAKE_1_DATATAKE_SENSING_START=2021-02-02T03:29:41.024Z
  DATATAKE_1_DATATAKE_TYPE=INS-NOBS
  DATATAKE_1_ID=GS2A_20210202T032941_029330_N99.99
  DATATAKE_1_SENSING_ORBIT_DIRECTION=DESCENDING
  DATATAKE_1_SENSING_ORBIT_NUMBER=18
  DATATAKE_1_SPACECRAFT_NAME=Sentinel-2A
  DEGRADED_ANC_DATA_PERCENTAGE=0.0
  DEGRADED_MSI_DATA_PERCENTAGE=0
  FOOTPRINT=POLYGON((108.8262199395114 34.3223655300343, 110.01937756706981 34.3377366079947, 110.03062458149692 33.347577738129736, 108.85113286341723 33.3327671075422, 108.8262199395114 34.3223655300343))
  FORMAT_CORRECTNESS=PASSED
  GENERAL_QUALITY=PASSED
  GENERATION_TIME=2021-12-10T13:10:28.715000Z
  GEOMETRIC_QUALITY=PASSED
  HIGH_PROBA_CLOUDS_PERCENTAGE=0.061825
  MEDIUM_PROBA_CLOUDS_PERCENTAGE=0.484693
  NODATA_PIXEL_PERCENTAGE=0.000043
  NOT_VEGETATED_PERCENTAGE=47.824116
  PREVIEW_GEO_INFO=Not applicable
  PREVIEW_IMAGE_URL=Not applicable
  PROCESSING_BASELINE=99.99
  PROCESSING_LEVEL=Level-2A
  PRODUCT_START_TIME=2021-02-02T03:29:41.024Z
  PRODUCT_STOP_TIME=2021-02-02T03:29:41.024Z
  PRODUCT_TYPE=S2MSI2A
  PRODUCT_URI=S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE
  RADIATIVE_TRANSFER_ACCURACY=0.0
  RADIOMETRIC_QUALITY=PASSED
  REFLECTANCE_CONVERSION_U=1.03080315559675
  SATURATED_DEFECTIVE_PIXEL_PERCENTAGE=0.000000
  SENSOR_QUALITY=PASSED
  SNOW_ICE_PERCENTAGE=0.002857
  SPECIAL_VALUE_NODATA=0
  SPECIAL_VALUE_SATURATED=65535
  THIN_CIRRUS_PERCENTAGE=0.502809
  UNCLASSIFIED_PERCENTAGE=10.570553
  VEGETATION_PERCENTAGE=6.625548
  WATER_PERCENTAGE=2.496857
  WATER_VAPOUR_RETRIEVAL_ACCURACY=0.0
  WVP_QUANTIFICATION_VALUE=1000.0
  WVP_QUANTIFICATION_VALUE_UNIT=cm
Subdatasets:
  SUBDATASET_1_NAME=SENTINEL2_L2A:S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml:10m:EPSG_32649
  SUBDATASET_1_DESC=Bands B2, B3, B4, B8 with 10m resolution, UTM 49N
 # 上述两行，第一行是子数据的位置信息，第二行是子数据的波段信息
    SUBDATASET_2_NAME=SENTINEL2_L2A:S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml:20m:EPSG_32649
  SUBDATASET_2_DESC=Bands B5, B6, B7, B8A, B11, B12, AOT, CLD, SCL, SNW, WVP with 20m resolution, UTM 49N
  SUBDATASET_3_NAME=SENTINEL2_L2A:S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml:60m:EPSG_32649
  SUBDATASET_3_DESC=Bands B1, B9, AOT, CLD, SCL, SNW, WVP with 60m resolution, UTM 49N
  SUBDATASET_4_NAME=SENTINEL2_L2A:S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml:TCI:EPSG_32649
  SUBDATASET_4_DESC=True color image, UTM 49N
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)

      $ gdalinfo SENTINEL2_L2A:S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml:10m:EPSG_32649
# 上面一行代码的作用是获取第一个子数据集的信息，子数据集的位置在上一节代码中有注释
ERROR 1: PROJ: proj_create_from_database: Cannot find proj.db
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
Driver: SENTINEL2/Sentinel 2
Files: S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/MTD_TL.xml
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R10m/T49SCT_20210202T032941_B04_10m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R10m/T49SCT_20210202T032941_B03_10m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R10m/T49SCT_20210202T032941_B02_10m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R10m/T49SCT_20210202T032941_B08_10m.jp2
# 上面的输出结果显示了波段重组后的信息
Size is 10980, 10980
Origin = (300000.000000000000000,3800040.000000000000000)
Pixel Size = (10.000000000000000,-10.000000000000000)
Metadata:
  AOT_QUANTIFICATION_VALUE=1000.0
  AOT_QUANTIFICATION_VALUE_UNIT=none
  AOT_RETRIEVAL_ACCURACY=0.0
  BOA_QUANTIFICATION_VALUE=10000
  BOA_QUANTIFICATION_VALUE_UNIT=none
  CLOUD_COVERAGE_ASSESSMENT=1.049327
  CLOUD_SHADOW_PERCENTAGE=0.440092
  DARK_FEATURES_PERCENTAGE=30.990644
  DATATAKE_1_DATATAKE_SENSING_START=2021-02-02T03:29:41.024Z
  DATATAKE_1_DATATAKE_TYPE=INS-NOBS
  DATATAKE_1_ID=GS2A_20210202T032941_029330_N99.99
  DATATAKE_1_SENSING_ORBIT_DIRECTION=DESCENDING
  DATATAKE_1_SENSING_ORBIT_NUMBER=18
  DATATAKE_1_SPACECRAFT_NAME=Sentinel-2A
  DEGRADED_ANC_DATA_PERCENTAGE=0.0
  DEGRADED_MSI_DATA_PERCENTAGE=0
  FORMAT_CORRECTNESS=PASSED
  GENERAL_QUALITY=PASSED
  GENERATION_TIME=2021-12-10T13:10:28.715000Z
  GEOMETRIC_QUALITY=PASSED
  HIGH_PROBA_CLOUDS_PERCENTAGE=0.061825
  MEDIUM_PROBA_CLOUDS_PERCENTAGE=0.484693
  NODATA_PIXEL_PERCENTAGE=0.000043
  NOT_VEGETATED_PERCENTAGE=47.824116
  PREVIEW_GEO_INFO=Not applicable
  PREVIEW_IMAGE_URL=Not applicable
  PROCESSING_BASELINE=99.99
  PROCESSING_LEVEL=Level-2A
  PRODUCT_START_TIME=2021-02-02T03:29:41.024Z
  PRODUCT_STOP_TIME=2021-02-02T03:29:41.024Z
  PRODUCT_TYPE=S2MSI2A
  PRODUCT_URI=S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE
  RADIATIVE_TRANSFER_ACCURACY=0.0
  RADIOMETRIC_QUALITY=PASSED
  REFLECTANCE_CONVERSION_U=1.03080315559675
  SATURATED_DEFECTIVE_PIXEL_PERCENTAGE=0.000000
  SENSOR_QUALITY=PASSED
  SNOW_ICE_PERCENTAGE=0.002857
  SPECIAL_VALUE_NODATA=0
  SPECIAL_VALUE_SATURATED=65535
  THIN_CIRRUS_PERCENTAGE=0.502809
  UNCLASSIFIED_PERCENTAGE=10.570553
  VEGETATION_PERCENTAGE=6.625548
  WATER_PERCENTAGE=2.496857
  WATER_VAPOUR_RETRIEVAL_ACCURACY=0.0
  WVP_QUANTIFICATION_VALUE=1000.0
  WVP_QUANTIFICATION_VALUE_UNIT=cm
Image Structure Metadata:
  COMPRESSION=JPEG2000
Corner Coordinates:
Upper Left  (  300000.000, 3800040.000)
Lower Left  (  300000.000, 3690240.000)
Upper Right (  409800.000, 3800040.000)
Lower Right (  409800.000, 3690240.000)
Center      (  354900.000, 3745140.000)
# 下面的代码详细介绍了波段重组后各波段的信息
Band 1 Block=128x128 Type=UInt16, ColorInterp=Red
  Description = B4, central wavelength 665 nm
  Overviews: 5490x5490, 2745x2745, 1373x1373, 687x687, 344x344
  Metadata:
    BANDNAME=B4
    BANDWIDTH=30
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=1512.06
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=665
    WAVELENGTH_UNIT=nm
Band 2 Block=128x128 Type=UInt16, ColorInterp=Green
  Description = B3, central wavelength 560 nm
  Overviews: 5490x5490, 2745x2745, 1373x1373, 687x687, 344x344
  Metadata:
    BANDNAME=B3
    BANDWIDTH=35
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=1823.24
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=560
    WAVELENGTH_UNIT=nm
Band 3 Block=128x128 Type=UInt16, ColorInterp=Blue
  Description = B2, central wavelength 490 nm
  Overviews: 5490x5490, 2745x2745, 1373x1373, 687x687, 344x344
  Metadata:
    BANDNAME=B2
    BANDWIDTH=65
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=1959.66
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=490
    WAVELENGTH_UNIT=nm
Band 4 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = B8, central wavelength 842 nm
  Overviews: 5490x5490, 2745x2745, 1373x1373, 687x687, 344x344
  Metadata:
    BANDNAME=B8
    BANDWIDTH=115
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=1041.63
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=842
    WAVELENGTH_UNIT=nm

      $ gdalinfo SENTINEL2_L2A:S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml:20m:EPSG_32649
ERROR 1: PROJ: proj_create_from_database: Cannot find proj.db
Driver: SENTINEL2/Sentinel 2
Files: S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/MTD_TL.xml
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R20m/T49SCT_20210202T032941_B05_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R20m/T49SCT_20210202T032941_B06_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R20m/T49SCT_20210202T032941_B07_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R20m/T49SCT_20210202T032941_B8A_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R20m/T49SCT_20210202T032941_B11_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R20m/T49SCT_20210202T032941_B12_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R20m/T49SCT_20210202T032941_AOT_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/QI_DATA/MSK_CLDPRB_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R20m/T49SCT_20210202T032941_SCL_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/QI_DATA/MSK_SNWPRB_20m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R20m/T49SCT_20210202T032941_WVP_20m.jp2
Size is 5490, 5490
Origin = (300000.000000000000000,3800040.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)
Metadata:
  AOT_QUANTIFICATION_VALUE=1000.0
  AOT_QUANTIFICATION_VALUE_UNIT=none
  AOT_RETRIEVAL_ACCURACY=0.0
  BOA_QUANTIFICATION_VALUE=10000
  BOA_QUANTIFICATION_VALUE_UNIT=none
  CLOUD_COVERAGE_ASSESSMENT=1.049327
  CLOUD_SHADOW_PERCENTAGE=0.440092
  DARK_FEATURES_PERCENTAGE=30.990644
  DATATAKE_1_DATATAKE_SENSING_START=2021-02-02T03:29:41.024Z
  DATATAKE_1_DATATAKE_TYPE=INS-NOBS
  DATATAKE_1_ID=GS2A_20210202T032941_029330_N99.99
  DATATAKE_1_SENSING_ORBIT_DIRECTION=DESCENDING
  DATATAKE_1_SENSING_ORBIT_NUMBER=18
  DATATAKE_1_SPACECRAFT_NAME=Sentinel-2A
  DEGRADED_ANC_DATA_PERCENTAGE=0.0
  DEGRADED_MSI_DATA_PERCENTAGE=0
  FORMAT_CORRECTNESS=PASSED
  GENERAL_QUALITY=PASSED
  GENERATION_TIME=2021-12-10T13:10:28.715000Z
  GEOMETRIC_QUALITY=PASSED
  HIGH_PROBA_CLOUDS_PERCENTAGE=0.061825
  MEDIUM_PROBA_CLOUDS_PERCENTAGE=0.484693
  NODATA_PIXEL_PERCENTAGE=0.000043
  NOT_VEGETATED_PERCENTAGE=47.824116
  PREVIEW_GEO_INFO=Not applicable
  PREVIEW_IMAGE_URL=Not applicable
  PROCESSING_BASELINE=99.99
  PROCESSING_LEVEL=Level-2A
  PRODUCT_START_TIME=2021-02-02T03:29:41.024Z
  PRODUCT_STOP_TIME=2021-02-02T03:29:41.024Z
  PRODUCT_TYPE=S2MSI2A
  PRODUCT_URI=S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE
  RADIATIVE_TRANSFER_ACCURACY=0.0
  RADIOMETRIC_QUALITY=PASSED
  REFLECTANCE_CONVERSION_U=1.03080315559675
  SATURATED_DEFECTIVE_PIXEL_PERCENTAGE=0.000000
  SENSOR_QUALITY=PASSED
  SNOW_ICE_PERCENTAGE=0.002857
  SPECIAL_VALUE_NODATA=0
  SPECIAL_VALUE_SATURATED=65535
  THIN_CIRRUS_PERCENTAGE=0.502809
  UNCLASSIFIED_PERCENTAGE=10.570553
  VEGETATION_PERCENTAGE=6.625548
  WATER_PERCENTAGE=2.496857
  WATER_VAPOUR_RETRIEVAL_ACCURACY=0.0
  WVP_QUANTIFICATION_VALUE=1000.0
  WVP_QUANTIFICATION_VALUE_UNIT=cm
Image Structure Metadata:
  COMPRESSION=JPEG2000
Corner Coordinates:
Upper Left  (  300000.000, 3800040.000)
Lower Left  (  300000.000, 3690240.000)
Upper Right (  409800.000, 3800040.000)
Lower Right (  409800.0ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
00, 3690240.000)
Center      (  354900.000, 3745140.000)
Band 1 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = B5, central wavelength 705 nm
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=B5
    BANDWIDTH=15
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=1424.64
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=705
    WAVELENGTH_UNIT=nm
Band 2 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = B6, central wavelength 740 nm
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=B6
    BANDWIDTH=15
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=1287.61
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=740
    WAVELENGTH_UNIT=nm
Band 3 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = B7, central wavelength 783 nm
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=B7
    BANDWIDTH=20
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=1162.08
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=783
    WAVELENGTH_UNIT=nm
Band 4 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = B8A, central wavelength 865 nm
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=B8A
    BANDWIDTH=20
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=955.32
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=865
    WAVELENGTH_UNIT=nm
Band 5 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = B11, central wavelength 1610 nm
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=B11
    BANDWIDTH=90
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=245.59
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=1610
    WAVELENGTH_UNIT=nm
Band 6 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = B12, central wavelength 2190 nm
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=B12
    BANDWIDTH=180
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=85.25
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=2190
    WAVELENGTH_UNIT=nm
Band 7 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = AOT, Aerosol Optical Thickness map (at 550nm)
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=AOT
Band 8 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = CLD, Raster mask values range from 0 for high confidence clear sky to 100 for high confidence cloudy
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=CLD
Band 9 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = SCL, Scene Classification
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Categories:
      0: NODATA
      1: SATURATED_DEFECTIVE
      2: DARK_FEATURE_SHADOW
      3: CLOUD_SHADOW
      4: VEGETATION
      5: NOT_VEGETATED
      6: WATER
      7: UNCLASSIFIED
      8: CLOUD_MEDIUM_PROBA
      9: CLOUD_HIGH_PROBA
     10: THIN_CIRRUS
     11: SNOW_ICE
  Metadata:
    BANDNAME=SCL
Band 10 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = SNW, Raster mask values range from 0 for high confidence NO snow/ice to 100 for high confidence snow/ice
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=SNW
Band 11 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = WVP, Scene-average Water Vapour map
  Overviews: 2745x2745, 1373x1373, 687x687, 344x344, 172x172
  Metadata:
    BANDNAME=WVP

      $ gdalinfo SENTINEL2_L2A:S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml:60m:EPSG_32649
ERROR 1: PROJ: proj_create_from_database: Cannot find proj.db
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
ERROR 1: PROJ: proj_create: no database context specified
Driver: SENTINEL2/Sentinel 2
Files: S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/MTD_MSIL2A.xml
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/MTD_TL.xml
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R60m/T49SCT_20210202T032941_B01_60m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R60m/T49SCT_20210202T032941_B09_60m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R60m/T49SCT_20210202T032941_AOT_60m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/QI_DATA/MSK_CLDPRB_60m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R60m/T49SCT_20210202T032941_SCL_60m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/QI_DATA/MSK_SNWPRB_60m.jp2
       S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE/GRANULE/L2A_T49SCT_A029330_20210202T032944/IMG_DATA/R60m/T49SCT_20210202T032941_WVP_60m.jp2
Size is 1830, 1830
Origin = (300000.000000000000000,3800040.000000000000000)
Pixel Size = (60.000000000000000,-60.000000000000000)
Metadata:
  AOT_QUANTIFICATION_VALUE=1000.0
  AOT_QUANTIFICATION_VALUE_UNIT=none
  AOT_RETRIEVAL_ACCURACY=0.0
  BOA_QUANTIFICATION_VALUE=10000
  BOA_QUANTIFICATION_VALUE_UNIT=none
  CLOUD_COVERAGE_ASSESSMENT=1.049327
  CLOUD_SHADOW_PERCENTAGE=0.440092
  DARK_FEATURES_PERCENTAGE=30.990644
  DATATAKE_1_DATATAKE_SENSING_START=2021-02-02T03:29:41.024Z
  DATATAKE_1_DATATAKE_TYPE=INS-NOBS
  DATATAKE_1_ID=GS2A_20210202T032941_029330_N99.99
  DATATAKE_1_SENSING_ORBIT_DIRECTION=DESCENDING
  DATATAKE_1_SENSING_ORBIT_NUMBER=18
  DATATAKE_1_SPACECRAFT_NAME=Sentinel-2A
  DEGRADED_ANC_DATA_PERCENTAGE=0.0
  DEGRADED_MSI_DATA_PERCENTAGE=0
  FORMAT_CORRECTNESS=PASSED
  GENERAL_QUALITY=PASSED
  GENERATION_TIME=2021-12-10T13:10:28.715000Z
  GEOMETRIC_QUALITY=PASSED
  HIGH_PROBA_CLOUDS_PERCENTAGE=0.061825
  MEDIUM_PROBA_CLOUDS_PERCENTAGE=0.484693
  NODATA_PIXEL_PERCENTAGE=0.000043
  NOT_VEGETATED_PERCENTAGE=47.824116
  PREVIEW_GEO_INFO=Not applicable
  PREVIEW_IMAGE_URL=Not applicable
  PROCESSING_BASELINE=99.99
  PROCESSING_LEVEL=Level-2A
  PRODUCT_START_TIME=2021-02-02T03:29:41.024Z
  PRODUCT_STOP_TIME=2021-02-02T03:29:41.024Z
  PRODUCT_TYPE=S2MSI2A
  PRODUCT_URI=S2A_MSIL2A_20210202T032941_N9999_R018_T49SCT_20211210T131028.SAFE
  RADIATIVE_TRANSFER_ACCURACY=0.0
  RADIOMETRIC_QUALITY=PASSED
  REFLECTANCE_CONVERSION_U=1.03080315559675
  SATURATED_DEFECTIVE_PIXEL_PERCENTAGE=0.000000
  SENSOR_QUALITY=PASSED
  SNOW_ICE_PERCENTAGE=0.002857
  SPECIAL_VALUE_NODATA=0
  SPECIAL_VALUE_SATURATED=65535
  THIN_CIRRUS_PERCENTAGE=0.502809
  UNCLASSIFIED_PERCENTAGE=10.570553
  VEGETATION_PERCENTAGE=6.625548
  WATER_PERCENTAGE=2.496857
  WATER_VAPOUR_RETRIEVAL_ACCURACY=0.0
  WVP_QUANTIFICATION_VALUE=1000.0
  WVP_QUANTIFICATION_VALUE_UNIT=cm
Image Structure Metadata:
  COMPRESSION=JPEG2000
Corner Coordinates:
Upper Left  (  300000.000, 3800040.000)
Lower Left  (  300000.000, 3690240.000)
Upper Right (  409800.000, 3800040.000)
Lower Right (  409800.000, 3690240.000)
Center      (  354900.000, 3745140.000)
Band 1 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = B1, central wavelength 443 nm
  Overviews: 915x915, 458x458, 229x229
  Metadata:
    BANDNAME=B1
    BANDWIDTH=20
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=1884.69
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=443
    WAVELENGTH_UNIT=nm
Band 2 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = B9, central wavelength 945 nm
  Overviews: 915x915, 458x458, 229x229
  Metadata:
    BANDNAME=B9
    BANDWIDTH=20
    BANDWIDTH_UNIT=nm
    SOLAR_IRRADIANCE=812.92
    SOLAR_IRRADIANCE_UNIT=W/m2/um
    WAVELENGTH=945
    WAVELENGTH_UNIT=nm
Band 3 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = AOT, Aerosol Optical Thickness map (at 550nm)
  Overviews: 915x915, 458x458, 229x229
  Metadata:
    BANDNAME=AOT
Band 4 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = CLD, Raster mask values range from 0 for high confidence clear sky to 100 for high confidence cloudy
  Overviews: 915x915, 458x458, 229x229
  Metadata:
    BANDNAME=CLD
Band 5 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = SCL, Scene Classification
  Overviews: 915x915, 458x458, 229x229
  Categories:
      0: NODATA
      1: SATURATED_DEFECTIVE
      2: DARK_FEATURE_SHADOW
      3: CLOUD_SHADOW
      4: VEGETATION
      5: NOT_VEGETATED
      6: WATER
      7: UNCLASSIFIED
      8: CLOUD_MEDIUM_PROBA
      9: CLOUD_HIGH_PROBA
     10: THIN_CIRRUS
     11: SNOW_ICE
  Metadata:
    BANDNAME=SCL
Band 6 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = SNW, Raster mask values range from 0 for high confidence NO snow/ice to 100 for high confidence snow/ice
  Overviews: 915x915, 458x458, 229x229
  Metadata:
    BANDNAME=SNW
Band 7 Block=128x128 Type=UInt16, ColorInterp=Undefined
  Description = WVP, Scene-average Water Vapour map
  Overviews: 915x915, 458x458, 229x229
  Metadata:
    BANDNAME=WVP

  • rgb 图像展示：
    19
    20

    • 通过改变代码中的子数据集以及 gdalinfo 工具查询的信息可以获取 Sentinel-2 数据集的不同波段。

如果本项目帮助到了你，不如打赏作者一杯奶茶吧~~~

BeZer0
关注
  • 20
    点赞
  • 踩
  • 122
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 14
    评论
专栏目录
Sentinel-2哨兵二号数据下载及处理教程.pdf
06-10
Sentinel-2哨兵二号数据下载及处理教程.pdf Sentinel-2 是高分辨率多光谱成像卫星，一颗卫星的重访周期为10天，两颗互补，重访周期为5天。分为2A和2B两颗卫星。2A于2015年6月23日01:52 UTC以“织女星”运载火箭发射升空。2B于2017年3月07日北京时间9时49分 UTC以“织女星”运载火箭发射升空。两者同时进入运行状态后，每5天可完成一次对地球赤道地区的完整成像，而对于纬度较高的地区，这一周期仅需3天。 Sentinel-2 卫星携带一枚多光谱仪器(MSI)，可覆盖13个光谱波段，地面分辨率分别有10m、20m和60m：
哨兵2号（Sentinel-2）介绍、下载、预处理及批处理
热门推荐
GeoSuper的博客
04-17 35万+
哨兵2号是高分辨率多光谱成像卫星，携带一枚多光谱成像仪(MSI)，用于陆地监测，可提供植被、土壤和水覆盖、内陆水路及海岸区域等图像，还可用于紧急救援服务。分为2A和2B两颗卫星。 第一颗卫星哨兵2号A于2015年6月23日01:52 UTC以“织女星”运载火箭发射升空。 6月29日，在轨运行4天的哨兵-2A卫星，传回了第一景数据，幅宽290km，卫星第一次扫描的范围是从瑞典开始，经过中欧和地中.........
【RS】Sentinel-2哨兵二号L1C数据预处理教程（Sen2Cor大气校正、SNAP重采样、ENVI波段组合）
Remote Sensing迷途小书童
03-06 5651
刚分享过新版欧空局的数据下载教程，正好把哨兵2号预处理的教程也分享一下，主要就是使用官方插件Sen2or对L1C级数据进行大气校正，使用SNAP软件对L2A级数据进行重采样并导出ENVI可以打开的格式，最后使用ENVI对重采样后的数据进行波段组合，以便于后期的定量分析。
Sentinel-2（哨兵2号）SNAP-ENVI预处理
最新发布
weixin_63175072的博客
06-02 880
这一步主要是为了将所有波段的分辨率统一，因为哨兵2号的波段空间分辨率有10m、20m和60m，不便于后续分析。进行处理，可存为ENVI或TIFF格式，即可用ENVI进行打开。右键点击文件，选择Open RGB window可以查看彩色图像。一、由SNAP打开下载好的Sentinel-2数据。选择resampling parameters，由SNAP打开解压后MTD开头的XML文件。确定待处理数据、命名、保存格式、路径；选择I/0 parameters，调节参数，点击run，等待完成。
【RS】哨兵系列数据下载（哨兵2号Sentinel-2波段介绍、下载）
Remote Sensing迷途小书童
06-16 7670
以哨兵二号为示例，分享如何下载欧空局的哨兵系列数据
基于ENVI5.6对哨兵2数据进行去云操作
Pyyyson的博客
03-11 1725
基于ENVI5.6对哨兵2数据进行去云操作
Sentinel-2 哨兵二号数据（Level-1C）下载及预处理教程
C_hqpxd的博客
05-15 3万+
欧空局（ESA）仅发布了哨兵2号（S2）的L1C级多光谱数据（MSI），L1C级数据是经过几何精校正的正射影像，并没有进行辐射定标和大气校正。同时，ESA还对S2 L2A级数据进行了定义，L2A级数据主要包含经过辐射定标和大气校正的大气底层反射率数据（Bottom-of-Atmosphere corrected reflectance），但这个L2A数据需要用户根据需求自行生产，为此，ESA发布了专门生产L2A级数据的插件Sen2cor。 ...
GIS遥感专栏03：Sentinel-2影像获取及处理流程
weixin_55429615的博客
06-03 6976
哨兵2号（Sentinel-2）是欧洲空间局哥白尼计划下的一个地球观测任务，该任务主要对地球表面进行观测以提供相关遥测服务，例如森林监测、土地覆盖变化侦测、天然灾害管理。该计划是由2颗相同的卫星哨兵2号A（Sentinel-2A） 与B（Sentinel-2B）组成的卫星群。为了能在短时间内重复拍摄同一区域，和高任务可用性，哨兵2号任务包含了两颗同时运作的相同卫星哨兵2号A（Sentinel-2A） 与B（Sentinel-2B）。
哨兵2号预处理
m0_60674379的博客
03-30 3567
打开ENVI，在Toolbox选择Raster Management，找到Layer Stacking（波段融合），将想要融合的波段加载进来，点击Reorder Files将波段进行排序，选择投影坐标和输出路径，选择重采样方式，最后确定。同时，将输出格式选择为ENVI，确定输出位置，点击Run。PS：如果显示图像offline（时间太久远就会被暂时下架），需要申请上线（添加购物车，点下载后会出现提示栏，绿色说明已提交申请，红色是提交失败，需要过会儿再试试），大概两三个小时就可以下载了。（这样说明处理完成）
哨兵2号数据介绍
clever101的专栏
12-19 1万+
1、数据介绍   哨兵2号是高分辨率多光谱成像卫星，携带一枚多光谱成像仪(MSI)，用于陆地监测，可提供植被、土壤和水覆盖、内陆水路及海岸区域等图像，还可用于紧急救援服务。分为2A和2B两颗卫星。   第一颗卫星哨兵2号A于2015年6月23日01:52 UTC以“织女星”运载火箭发射升空。   6月29日，在轨运行4天的哨兵-2A卫星，传回了第一景数据，幅宽290km，卫星第一次扫描的范围是从瑞典开始，经过中欧和地中海，到阿尔及利亚结束。   第二颗卫星哨兵2号B于2017年3月07日北京时间9时49分
关于哨兵sentinel 2所有知识都在这里
m0_57191404的博客
01-09 1万+
目录 Sentinel2介绍 处理工具 光谱特性 大气校正 文件结构及命名 全色锐化 数据下载 Sentinel2介绍 Sentinel 2是高分辨率多光谱成像卫星，搭载多光谱成像仪（MSI），由“双胞胎” Sentinel 2A 和Sentinel 2B 两颗卫星组成，它们的轨道彼此相差180º，覆盖13个工作波段，放置在地球的两端。图像采集时间分辨率是每颗卫星10 天，两颗就是5天。因此，你将在不超过五天的时间段内获得同一区域更新的卫片。卫星拍摄的每幅图像都显示了290公里..
哨兵二号预处理-sentinel-2A
04-23
关于对哨兵二号预处理。
遥感【SAR数据处理】【SNAP使用指导】【哨兵1号2号】如何预处理sentinel-1和2数据
06-10
SNAP是一款强大的遥感数据处理工具，它集成了多种处理模块，支持Sentinel系列卫星的数据处理。无论是Sentinel-1的SAR数据还是Sentinel-2的光学数据，SNAP都能提供完整的预处理流程，包括辐射定标、几何校正、滤波等...
Sentinel2-L2A数据的下载、读取及基础操作.pdf
06-10
Sentinel2卫星是欧洲航天局（ESA）推出的高分辨率多光谱成像卫星，主要用于陆地监测、灾害响应以及农业、环境等领域的应用。Sentinel2由两颗卫星组成，即Sentinel2A和Sentinel2B，它们分别在2015年和2017年发射升空...
Sentinel-5P卫星1级数据介绍
08-31
Sentinel-5P（哨兵5号前期）卫星是一项由欧洲航天局（ESA）主导的重要任务，旨在监测地球大气中的关键成分，以提供全球空气质量、气候影响及平流层臭氧层状况的信息。这颗卫星在低地球轨道上运行，其主要载荷是荷兰...
中科星图——Sentinel-2_MSI_L2A数据集
qq_31988139的博客
01-23 2298
哨兵2号（Sentinel-2）卫星是高分辨率多光谱成像卫星，携带一枚多光谱成像仪（MSI），分为2A和2B两颗卫星，其中一颗卫星的重访周期为10天，两颗互补，重访周期为5天，主要用于包括陆地植被、土壤以及水资源、内河水道和沿海区在内的全球陆地观测。Sentinel-2_MSI_L2A数据集中的数据为L2A级产品数据，是经过大气校正的地表反射率产品，包括12个UINT16光谱带，还包含了场景分类产品、气溶胶光学厚度（AOT）、水蒸气及部分云和降雪概率产品，CLM波段是具有云掩码信息的位掩码频段。
Sentinel 2（哨兵2）卫星遥感数据的波段和组合介绍
qq_45245354的博客
08-04 3万+
摘要: Sentinel 2（哨兵2） 的波段共有 13 个，每个波段的像素大小为 10、20 或 60 米不等。Sentinel 2 由 2 颗卫星组成： 第一个是 2015 年发行的 Sentinel 2A 卫星，另一个是 2017 年发行的 Sentinel ...
Sentinel-2（哨兵2号）数据下载及预处理
baidu_34785556的博客
12-07 2万+
1.数据介绍： 哨兵2号是高分辨率多光谱成像卫星，携带一枚多光谱成像仪(MSI)，用于陆地监测，可提供植被、土壤和水覆盖、内陆水路及海岸区域等图像，还可用于紧急救援服务。分为2A和2B两颗卫星。 第一颗卫星哨兵2号A于2015年6月23日01:52 UTC以“织女星”运载火箭发射升空。 6月29日，在轨运行4天的哨兵-2A卫星，传回了第一景数据，幅宽290km，卫星第一次扫描的范围是从瑞典开始，经过中欧和地中海，到阿尔及利亚结束。 第二颗卫星哨兵2号B于2017年3月07日北京时间9时49分 UTC
哨兵二号数据处理
Little_fine的博客
12-02 1626
对下载的哨兵二号数据进行处理最终获取TIFF遥感影像
sentinel-2数据处理流程
07-16
处理 Sentinel-2 数据的一般流程如下： 1. 数据获取： Copernicus Open Access Hub 或其他数据源下载 Sentinel-2 数据。数据以 SAFE (Sentinel Application For Earth) 格式提供，其中包含多个文件和文件夹。 2. 数据解压缩：解压缩 SAFE 文件，得到原始数据文件。 3. 辐射校正：校正原始数据以消除大气影响，常用的方法有 Sen2Cor 和 iCOR。 4. 波段合成：对不同波段的数据进行合成，生成多光谱影像。Sentinel-2 数据包含多个波段，包括可见光、近红外和短波红外波段。 5. 几何校正：对影像进行几何校正以消除地球表面的变形，常用的方法有地理纠正和影像匹配。 6. 云掩膜：根据云覆盖信息生成云掩膜，以便后续分析中排除云影响。 7. 图像增强：应用各种增强技术，如对比度增强、直方图均衡化等，以提高图像质量和可视化效果。 8. 特征提取：根据应用需求，提取感兴趣的地物特征，比如植被指数 (NDVI)、土壤水分指数 (SWI) 等。 9. 数据分析：根据特定的应用场景，对处理后的数据进行分析和解释，如土地覆盖分类、环境监测等。 10. 结果输出：根据需求将分析结果输出为图像、报告或其他形式的数据产品。 请注意，具体的数据处理流程可能因应用需求和软件工具而有所不同。上述流程仅提供了一般的处理步骤。
评论 14
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

BeZer0

打赏一杯奶茶支持一下作者吧~~

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

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

打赏作者

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

抵扣说明：

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

余额充值