Skip to content

Instantly share code, notes, and snippets.

@jpswinski
Last active July 7, 2022 15:53
Show Gist options
  • Save jpswinski/419c05ab44bdb16e246e484015e10d49 to your computer and use it in GitHub Desktop.
Save jpswinski/419c05ab44bdb16e246e484015e10d49 to your computer and use it in GitHub Desktop.
Submits a large YAPC request that takes a lot of time to process
# Imports
from sliderule import icesat2
# Settings
url = "icesat2sliderule.org"
asset = "nsidc-s3"
verbose = True
# Initialize Client
icesat2.init(url, verbose=verbose)
# Region of Interest
poly=[-60.491806, -64.28095,
-56.623029, -69.028066,
-60.603624, -71.175799,
-64.825318, -69.704525,
-66.571939, -67.96301,
-63.715481, -65.673347,
-60.491806, -64.28095]
# Calclate Area of Region
gdf = icesat2.toregion(poly)["gdf"].to_crs({'init': 'epsg:32633'})
print("Area: ", gdf.area / 10**6 )
# Parameters
res=20
cycle=1
params= { 'poly':poly,
'cnf':0,
'len':res,
'res':res,
'ats':res/2,
'cnt':10,
'cycle':cycle,
'maxi': 1,
'yapc':{"score":190, "knn":0, "win_h":3, "win_x":15},
'pass_invalid':False }
# Perform ATL06 Request
D6=icesat2.atl06p(params, asset='nsidc-s3', version="005", resources=["ATL03_20181220163836_12680112_005_01.h5"])
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment