# Data fetched via HTTP from ['http://its-live-data.jpl.nasa.gov.s3.amazonaws.com/velocity_mosaic/landsat/v00.0/static/GRE_G0120_0000.nc']

# Data processed using the following steps:

gdalwarp -t_srs EPSG:3413 -r bilinear NETCDF:{input_dir}/GRE_G0120_0000.nc:ice {output_dir}/warped.tif

gdalwarp -cutline {assets_dir}/greenland_rectangle.geojson -crop_to_cutline -co COMPRESS=DEFLATE {input_dir}/warped.tif {output_dir}/velocity_mosaic_ice_mask.tif

gdal_translate -co TILED=YES -co COMPRESS=DEFLATE -co PREDICTOR=2 {input_dir}/velocity_mosaic_ice_mask.tif {output_dir}/compressed.tif

cp {input_dir}/compressed.tif {output_dir}/velocity_mosaic_ice_mask.tif && gdaladdo -r average {output_dir}/velocity_mosaic_ice_mask.tif 2 4 8 16