数据处理
数据处理
从官网下载的海洋卫星数据是一个经过两重压缩的压缩包：
H1D_OPER_CZI_L2A_20201102T053358_20201102T053453_02072_10.tar.gz
先把.gz格式包解压成.tar格式，压缩比最大接近20%（解压Python程序，注意tar只有打包功能没有压缩功能，gz才是压缩格式）。
HY-1A和HY-1B卫星均配置了海洋水色水温扫描仪、海岸带成像仪（coastal zone imager，CZI）；
HY-1C和HY-1D卫星均配置了水色水温扫描仪（Chinese ocean color and temperature scanner，COCTS）、海岸带成像仪（coastal z one imager，CZI）、紫外成像仪（u ltra-v iolet imager, UVI）、定标光谱仪（satellite calibration spectrometer，SCS）和 船 舶 自 动 识 别 系 统（a u t o m a t i c identification system，AIS）。
一景数据的文件为压缩包：
H1D_OPER_CZI_L2A_20201102T053358_20201102T053453_02072_10.tar.gz
需要对其进行解压（gz）和解包（tar）得到1个xml数据说明文件、4个波段的图片、1个hdf5格式文件：
H1D_OPER_CZI_L2A_20201102T053358_20201102T053453_02072_10/H1D_OPER_CZI_L2A_20201102T053358_20201102T053453_02072_10.h5
HY 1 C/D L2A产品有：
1 文件属性(SDS Global Attributes)
2 传感器波段特性(Sensor Band Parameters)
3 扫描行属性数据(Scan-Line Attributes)
4 图象定位数据(Navigation Data)
5 地球物理参数 (Geophysical Data)
使用import h5py导入hdf文件库，打开h5文件：
with h5py.File(h5_file, 'r') as f:
for key in f.keys():
print(key)
print(f[key].name)
print(f[key].values)
结果：
Geophysical Data
/Geophysical Data
Navigation Data
/Navigation Data
Scan Line Attributes
/Scan Line Attributes
Sensor Band parameters
/Sensor Band parameters
读取h5文件的key：01 |
GDAL读取tif文件：02 |
GDAL输出tif文件：03 |
输出真彩色tif文件：04
输出图片（带格式）：05 （logo/向北针图片）
输出图片（不带格式）：06
CZI 通道与应用对象
波段编号 波段/μm 应用对象
1 0.42~0.50 叶绿素、污染、冰、浅海地形
2 0.52~0.60 叶绿素、低浓度泥沙、污染、滩涂
3 0.61~0.69 中等浓度泥沙、植被、土壤
4 0.76~0.89 植被、高浓度泥沙、大气校正
对应Rrc460、Rrc560、Rrc650、Rrc825
b460 = f['Geophysical Data']['Rrc460'][:, :]
可输出其大小和形状
print(b460.size)
print(b460.shape)
165680000
(7600, 21800)
