# Data fetched via HTTP from ['https://daacdata.apps.nsidc.org/pub/DATASETS/nsidc0611_seaice_age_v4/data/iceage_nh_12.5km_20100101_20101231_v4.1.nc']

# Data processed using the following steps:

gdal_translate -a_srs "+proj=laea +lat_0=90 +lon_0=0 +x_0=0 +y_0=0 +a=6371228 +b=6371228 +units=m +no_defs" -b 14 NETCDF:{input_dir}/iceage_nh_12.5km_20100101_20101231_v4.1.nc:age_of_sea_ice {output_dir}/age_of_sea_ice.tif

cp {input_dir}/age_of_sea_ice.tif {output_dir}/edited.tif && gdal_edit.py -a_ullr -4518421 4518421 4506579 -4506579 {output_dir}/edited.tif

gdalwarp -t_srs EPSG:3413 -r bilinear -tr 12500 12500 {input_dir}/edited.tif {output_dir}/warped.tif

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

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

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