Sentinel-6数据下载与读取

最新推荐文章于 2025-03-23 10:52:53 发布

原创

最新推荐文章于 2025-03-23 10:52:53 发布 · 1.7k 阅读

6 ·

CC 4.0 BY-SA版权

文章标签：

#sentinel #python

本文介绍了如何从欧洲气象卫星中心下载哨兵6号卫星的雷达高度计数据，以及使用Python读取处理这些数据的过程，包括数据文件结构和处理缺失值的方法。

文章目录

前言
一、下载数据
二、读取数据（以Python为例）

前言

哨兵6号卫星(Sentinel-6)携带了先进的雷达高度计，用于监测全球和区域的海平面变化和海平面上升速度。它建立在Jason系列、CryoSat-2和哥白尼哨兵3号卫星的基础之上，采用由ESA开发的Poseidon-4双频(Ku波段和C波段)雷达高度计，提高了Jason系列雷达高度计的精度和稳定性的同时继承了CryoSat和哥白尼Sentinel-3的高分辨率合成孔径雷达模式，其设计允许合成孔径雷达和低分辨率模式同时提供数据。

一、下载数据

哨兵6数据并不在欧空局数据中心，而是在欧洲气象卫星中心网站（https://archive.eumetsat.int/）。
注册完成后进入data center
在这里插入图片描述每个产品有两个测量数据文件可用（标准和简化），每个都有不同数量的变量。标准数据文件包括Ku波段和C波段的1 Hz和20 Hz测量值以及1 Hz和一些20 Hz的地球物理校正值。简化的数据文件仅包含Ku和C波段的1Hz测量值以及1Hz的地球物理校正。

HR与LR的区别：HR数据产品仅包含Ku波段测量值
NTC/NRT/STC：产品于< 60 days/< 3 hours/< 36小时内处理结果

选取时间、研究区地理位置后可以在DETAILS中查看数据的轨道号。
在这里插入图片描述提交订单后，即可查看订单状态并受到邮件通知

数据准备好后，会在邮箱收到下载链接。

[在这里插入图片描述](https://img-blog.csdnimg.cn/614ee03f0f2e487ebe0ab82154ed933e.png

二、读取数据（以Python为例）

哨兵6数据的存储方式与普通nc文件相比有一点不同，它将数据分成了20hz和1hz两个group，使用
matlab读取哨兵6比较麻烦，在此以python为例进行读取。

from pathlib import Path
import netCDF4 as nc
import numpy as np
import csv

filepath = Path(r'F:/S6_2022/243')
filelist = list(filepath.glob('*.nc'))
for file in filelist:
        file_info = nc.Dataset(file,mode = 'r')
        sig0_oceanku = file_info.groups['data_20'].groups['ku'].variables['sig0_ocean'][()]
        latku = file_info.groups['data_20'].groups['ku'].variables['latitude'][()]
        altku = file_info.groups