1.icesat-2 alt03数据下载方法
第一步:直接复制网址,点击进入https://nsidc.org/data/atl03/versions/5
第二步,登录;第三步,选取时间、空间范围
第四步,下拉,点击“订单文件”;第五步,点击确定
第六步,账号绑定邮箱会收到通知,然后直接进入https://nsidc.org/order-history,下拉,找到链接,点击下载
以上即是icesat-2 alt03数据下载方法。
2.icesat-2 ATL03数据字典
详情见:https://nsidc.org/sites/default/files/icesat2_atl03_data_dict_v003_0.pdf
可以看到字段名称、数据类型、名称描述、单位、数据详细描述等情况
关于icesat-2强弱波束判断的问题,目前网上的介绍较少,在这里引用一下孙伟老师等的方法[1]:
ATLAS地面轨迹沿卫星飞行方向从左至右依次被定义为gt1L、gt1R、gt2L、gt2R、gt3L、gt3R。ATLAS探测器每年将被重新定向两次,当ATLAS探测器定向为前向(与卫星飞行方向一致)时,gt1L、gt2L和gt3L对应每组波束左边的弱波束,gt1R、gt2R和gt3R对应每组波束的强波束;当ATLAS探测器定向为后向(与卫星飞行方向相反)时,gt1L、gt2L和gt3L对应每组波束左边的强波束,gt1R、gt2R和gt3R对应每组波束的弱波束。(如下图所示)
也给出本人的办法:ancillary_data/atlas_engineering/transmit/tx_pulse_energy,即读取h5文件·这里的内容进行波束强度判断,即可。
参考文献:
[1]孙伟,金建文,李国元,么嘉棋.激光测高卫星ICESat-2监测太湖水位精度评价[J].测绘科学,2021,46(11):6-11.DOI:10.16251/j.cnki.1009-2307.2021.11.002.
3.icesat-2 ATL03光子数据读取
这里给出python读取方法:
import pandas as pd
import h5py
import numpy as np
fname = r'ATL03_XXXXXXXXX_05070606_005_01.h5'#你的文件路径
group = ['/gt1l', '/gt1r', '/gt2l', '/gt2r', '/gt3l', '/gt3r']
for k, g in enumerate(group):#循环六个波束
with h5py.File(fname, 'r') as fi:
lat = fi[g + '/heights/lat_ph'][:]#读取光子纬度
time= fi[g +'/heights/delta_time'][:]#读取光子相对沿轨时间
height=fi[g +'/heights/h_ph'][:]#读取光子高度