#!/bin/bash

GDAL_TILE_PROCESSES=16
GDAL_TILE_ZOOMS=8-14
GDAL_SAMPLING_WARP=cubic
GDAL_SAMPLING_TILE=antialias

# Create virtual dataset with coordinates
gdal_translate -of VRT -a_srs EPSG:4326 -outsize 300% 300% -gcp 0.0 0.0 26.66072845 56.97695909 -gcp 21381.0 0.0 26.82578087 56.99059171 -gcp 21381.0 21498.0 26.84805393 56.89548042 -gcp 0.0 21498.0 26.67746544 56.88833084 GX825/GX_825_SK_frame_104_copy.jpg GX_825_SK_frame_104_copy.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "1077.0 20994.0\n594.0 20475.0\n0.0 19590.0\n0.0 11652.0\n77.25 11574.75\n276.0 11481.0\n285.0 10090.5\n105.75 9996.75\n0.0 9915.75\n0.0 3798.0\n46.5 2296.5\n529.5 2349.0\n1140.0 2364.0\n1377.0 2337.0\n1941.0 1968.0\n2115.0 1647.0\n2223.0 1311.0\n2265.0 228.0\n3054.0 225.0\n9486.0 294.0\n9741.0 312.0\n9939.0 630.0\n11364.0 634.5\n11476.5 429.0\n11601.0 331.5\n15514.5 357.0\n19512.0 400.5\n19503.0 2307.0\n21244.5 2329.5\n21141.0 10044.0\n21016.5 10186.5\n20851.5 10294.5\n20826.0 11007.0\n20835.0 11702.25\n21004.5 11812.5\n21111.0 11976.0\n21090.0 15861.0\n21108.0 17491.5\n21009.0 20775.0\n21018.0 21493.5\n11448.0 21495.0\n11157.0 21165.0\n9732.0 21150.0\n9501.0 21495.0\n4701.0 21426.0\n2352.0 21414.0\n1566.0 21396.0" | \
  gdaltransform -tps -output_xy GX_825_SK_frame_104_copy.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  GX_825_SK_frame_104_copy.vrt GX_825_SK_frame_104_copy.cut.vrt

# Generate tiles
gdal2tiles.py -r $GDAL_SAMPLING_TILE --xyz -z $GDAL_TILE_ZOOMS -x --processes=$GDAL_TILE_PROCESSES GX_825_SK_frame_104_copy.cut.vrt GX_825_SK_frame_104_copy.xyz