# Data fetched via HTTP from ['ftp://sidads.colorado.edu/DATASETS/NOAA/G02135/north/monthly/geotiff/03_Mar/N_201803_concentration_v3.0.tif']

# Data processed using the following steps:

gdal_calc.py --calc 'A / 10.0' -A {input_dir}/*.tif --outfile={output_dir}/downscaled.tif

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