#!/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 6801.352415026834 3627.151750972763 24.08246014753693 57.46669675333489 -gcp 31594.58810375671 2594.390077821012 24.665336583296195 57.47864986689594 -gcp 25828.298747763867 38430.56698165648 24.529163944060155 57.02672955146567 -gcp 3099.0661896243287 40552.02056698166 24 57 mapster/Estland-Lettland_100K_33_SW_Valmiera_(Wolmar)_ca_1941_incomplete_nnv8xdV_BN_Sygn._ZZK_S-4_152_A_0001.jpg Estland-Lettland_100K_33_SW_Valmiera_(Wolmar)_ca_1941_incomplete_nnv8xdV_BN_Sygn._ZZK_S-4_152_A_0001.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "3231.0143112701253 20820.03001667593\n3315.0143112701253 977.969983324069\n35177.96601073345 1125.0917176209005\n35234.83363148479 28036.589216231238\n33066.91949910554 28595.659811006117\n29441.406082289803 29754.146748193438\n27918.060822898035 30456.04002223457\n27611.95706618962 30928.669260700386\n27606.04293381038 31474.842690383546\n27738.110912343473 31703.706503613117\n27624.00357781753 32225.933296275707\n27624.035778175312 32772.09005002779\n27498.0 32988.206781545305\n27389.91771019678 33755.95664257921\n27486.03935599284 34397.6798221234\n27600.078711985687 34967.8165647582\n27582.032200357782 35610.01667593107\n27489.019677996424 36360.09338521401\n27605.735241502683 36839.726514730406\n28608.286225402502 38304.1600889383\n28822.597495527727 39819.08171206226\n28500.314847942755 39782.78488048916\n27779.982110912344 40103.97331851028\n27120.400715563504 40523.89994441357\n3099.0661896243287 40552.02056698166" | \
  gdaltransform -tps -output_xy Estland-Lettland_100K_33_SW_Valmiera_(Wolmar)_ca_1941_incomplete_nnv8xdV_BN_Sygn._ZZK_S-4_152_A_0001.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  Estland-Lettland_100K_33_SW_Valmiera_(Wolmar)_ca_1941_incomplete_nnv8xdV_BN_Sygn._ZZK_S-4_152_A_0001.vrt Estland-Lettland_100K_33_SW_Valmiera_(Wolmar)_ca_1941_incomplete_nnv8xdV_BN_Sygn._ZZK_S-4_152_A_0001.cut.vrt

# Generate tiles
gdal2tiles.py -r $GDAL_SAMPLING_TILE --xyz -z $GDAL_TILE_ZOOMS -x --processes=$GDAL_TILE_PROCESSES Estland-Lettland_100K_33_SW_Valmiera_(Wolmar)_ca_1941_incomplete_nnv8xdV_BN_Sygn._ZZK_S-4_152_A_0001.cut.vrt Estland-Lettland_100K_33_SW_Valmiera_(Wolmar)_ca_1941_incomplete_nnv8xdV_BN_Sygn._ZZK_S-4_152_A_0001.xyz