# 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:

gdal_calc.py --calc="A*B" --outfile={output_dir}/masked_velocity_mosaic.tif -A NETCDF:{input_dir}/GRE_G0120_0000.nc:v -B NETCDF:{input_dir}/GRE_G0120_0000.nc:ice

gdalwarp -t_srs EPSG:3413 -r bilinear {input_dir}/masked_velocity_mosaic.tif {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.tif

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

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