# Data fetched via HTTP from ['https://daacdata.apps.nsidc.org/pub/DATASETS/nsidc0611_seaice_age_v4/data/iceage_nh_12.5km_20130101_20131231_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 37 NETCDF:{input_dir}/iceage_nh_12.5km_20130101_20131231_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