#!/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.780681530981344 57.6845473403928 -gcp 14340.0 0.0 27.121539863616064 57.68558727675949 -gcp 14340.0 16359.0 27.1217403870652 57.477790069363756 -gcp 0.0 16359.0 26.782363888112833 57.47756705770168 225-Apukalns_1923.jpg 225-Apukalns_1923.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"7192.5 1441.8725013881176\n1936.125 1420.1256246529706\n1898.625 8004.750624652971\n1885.125 14587.123958911718\n7165.5 14594.62312604109\n12447.75 14608.87437534703\n12454.5 8042.626041088284\n12457.125 1472.8107301499167\" | \
  gdaltransform -tps -output_xy 225-Apukalns_1923.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  225-Apukalns_1923.vrt 225-Apukalns_1923.cut.vrt

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