ICESat-2 NSIDC Cloud Access

is2-nsidc-cloud.py

import subprocess
import s3fs
import h5py
import numpy as np

# Get Credentials

# See https://nsidc.org/support/how/how-do-i-programmatically-request-data-services#mac
inquiry = subprocess.run(['curl', "-b", "/home/ec2-user/.urs_cookies", \
                        "-c", "/home/ec2-user/.urs_cookies", \
                        "-L", "-n", "https://data.nsidc.earthdatacloud.nasa.gov/s3credentials"],\
                       stdout=subprocess.PIPE).stdout.decode('utf-8')
credentials = eval(inquiry)

# File created using search.earthdata.nasa.gov
filefile = 'earthdata-output-ronne.txt'
with open(filefile) as f:
    file_list = f.readlines()
print(file_list[0])

# Load the file
s3 = s3fs.S3FileSystem(key=credentials['accessKeyId'],
                       secret=credentials['secretAccessKey'],
                       token=credentials['sessionToken'])

s3url = 's3://nsidc-cumulus-prod-protected/ATLAS/ATL06/004/2019/06/26/ATL06_20190626143632_13640310_004_01.h5'

%time f = h5py.File(s3.open(s3url,'rb'),'r')
%time heights = np.array(f['/gt1l/land_ice_segments/h_li'])