哨兵二号（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
    评论
专栏目录
gee教程——sentinel-2数据的两种数据预（去云）处理方式QA60和SLC
qq_31988139的博客
04-18 710
哨兵二号卫星是欧洲空间局（ESA）发射的一颗地球观测卫星，主要用于监测地球表面的变化。在数据预处理阶段，去云是非常重要的一个步骤，因为云覆盖是影响遥感数据质量和可用性的主要因素之一。在这里，我们将具体介绍哨兵二号数据预处理阶段中的两种去云方式：QA60去云方式和SLC去云方式。QA60去云方式是一种基于光学遥感数据的去云方法。它是根据云覆盖的不同，将图像分为不同的云类别进行分类。在这种方法中，QA60是指一种云掩膜数据集，其中包含了60种云类别。
哨兵2号数据处理
年年的博客
01-17 4080
本文主要是针对像我一样，刚学遥感专业的同学们，对哨兵2号数据处理的一些心得，交流学习一下。
哨兵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】哨兵系列数据下载（哨兵2号Sentinel-2波段介绍、下载）
Remote Sensing迷途小书童
06-16 7622
以哨兵二号为示例，分享如何下载欧空局的哨兵系列数据
ENVI5.6加哨兵2火灾面积提取
最新发布
qq_65772557的博客
05-21 492
该方法提供三种方法去检测，每种方法出图也都不一样，可以都试试看效果，如果不懂可以选择左下角的蓝色问号，官方网站会介绍每一种方法，很详细的。3、绘制好以后，File→save as→ENVI ，之后出现Data Selection窗口，源文件包括11号和21号，所以裁剪区域同样要裁剪两次不同日期的底图，波段选择4，mask选择上一步绘制的感兴趣区域。每一种方法提取出来的都不一样，如果对结果不满意可以叠加不同方法所提取的结果，选择最适合的方式。数据来源，欧氏空间局需要注册登录，有点麻烦可以搜索教程多试几次。
【RS】Sentinel-2哨兵二号L1C数据预处理教程（Sen2Cor大气校正、SNAP重采样、ENVI波段组合）
Remote Sensing迷途小书童
03-06 5554
刚分享过新版欧空局的数据下载教程，正好把哨兵2号预处理的教程也分享一下，主要就是使用官方插件Sen2or对L1C级数据进行大气校正，使用SNAP软件对L2A级数据进行重采样并导出ENVI可以打开的格式，最后使用ENVI对重采样后的数据进行波段组合，以便于后期的定量分析。
GIS遥感专栏03：Sentinel-2影像获取及处理流程
weixin_55429615的博客
06-03 6900
哨兵2号（Sentinel-2）是欧洲空间局哥白尼计划下的一个地球观测任务，该任务主要对地球表面进行观测以提供相关遥测服务，例如森林监测、土地覆盖变化侦测、天然灾害管理。该计划是由2颗相同的卫星哨兵2号A（Sentinel-2A） 与B（Sentinel-2B）组成的卫星群。为了能在短时间内重复拍摄同一区域，和高任务可用性，哨兵2号任务包含了两颗同时运作的相同卫星哨兵2号A（Sentinel-2A） 与B（Sentinel-2B）。
基于ENVI5.6对哨兵2数据进行去云操作
Pyyyson的博客
03-11 1685
基于ENVI5.6对哨兵2数据进行去云操作
中科星图——Sentinel-2_MSI_L2A数据集
qq_31988139的博客
01-23 2278
哨兵2号（Sentinel-2）卫星是高分辨率多光谱成像卫星，携带一枚多光谱成像仪（MSI），分为2A和2B两颗卫星，其中一颗卫星的重访周期为10天，两颗互补，重访周期为5天，主要用于包括陆地植被、土壤以及水资源、内河水道和沿海区在内的全球陆地观测。Sentinel-2_MSI_L2A数据集中的数据为L2A级产品数据，是经过大气校正的地表反射率产品，包括12个UINT16光谱带，还包含了场景分类产品、气溶胶光学厚度（AOT）、水蒸气及部分云和降雪概率产品，CLM波段是具有云掩码信息的位掩码频段。
哨兵影像命名规则
u013622978的博客
06-27 1692
A：S1卫星的C波段合成孔径雷达（C-band Synthetic Aperture Radar，SAR）。- B：S1卫星的X波段合成孔径雷达（X-band Synthetic Aperture Radar，SAR）。- S1：哨兵-1，合成孔径雷达（Synthetic Aperture Radar，SAR）卫星。- S3：哨兵-3，海洋与陆地观测卫星。- S2：哨兵-2，多光谱成像卫星。- S5：哨兵-5，高光谱成像卫星。- S6：哨兵-6，地球辐射计卫星。- S7：哨兵-7，环境安全卫星。
Sentinel2 哨兵2数据下载方法（USGS）-史上最全讲解
xiatiandexia123的博客
12-23 4万+
本文讲述了从USGS官网下载Sentinel2的方法。主要包括USGS官网的注册和下载流程，以及选定研究区中所出现的问题、解决办法。
哨兵2号数据波段说明
09-22
哨兵2号是最近发射升空的卫星，相比于landSat等卫星分辨率高，最高可达10米，文件是详细的哨兵各个波段的说明，对于遥感方面想利用哨兵数据做研究的同志用处很大
哨兵数据处理教程.zip
05-12
哨兵二号数据下载后必然离不开数据预处理，一方面2017年以后的哨兵二号数据在ENVI5.3打不开，这个教程就包括如何让哨兵二号能在ENVI5.3中打开，以及哨兵二号处理软件的下载和预处理流程
哨兵二号预处理-sentinel-2A
04-23
哨兵二号（Sentinel-2）是欧洲航天局（ESA） Copernicus 计划的一部分，主要用于地球观测，提供高分辨率的多光谱图像，用于土地覆盖分类、农业、森林监测、灾害管理等多个领域。预处理是数据分析的重要步骤，它涉及...
Sentinel-2哨兵二号数据下载及处理教程.pdf
06-10
Sentinel-2哨兵二号卫星是欧洲航天局（ESA） Copernicus计划的一部分，用于地球观测，提供高分辨率的多光谱成像数据。它由两颗卫星组成，即Sentinel-2A和Sentinel-2B，它们的重访周期设计得非常短，以确保全球覆盖的...
遥感【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-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。 ...
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、付费专栏及课程。

余额充值