#!/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 2364.0 2436.0 21.326666666666668 56.5 -gcp 22059.0 2508.0 21.776666666666667 56.5 -gcp 22140.0 22764.0 21.776666666666667 56.25 -gcp 2310.0 22677.0 21.326666666666668 56.25 lvva/Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas__Priekules_apkart__6828__A1_L22_33_02.JPG Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas__Priekules_apkart__6828__A1_L22_33_02.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "2364.0 2436.0\n22059.0 2508.0\n22140.0 22764.0\n7857.0 22713.04497501388\n7686.0 22933.43503609106\n7573.5 22981.400888395336\n7266.0 23074.439200444198\n7141.5 23170.466685174903\n7003.5 23326.483342587453\n6883.5 23541.084952803998\n6789.0 23591.91837867851\n6589.5 23608.413381454746\n6471.0 23534.868406440866\n6369.0 23456.8450860633\n6244.5 23473.375902276515\n6195.0 23554.416712937258\n6117.0 23602.420877290395\n6000.0 23648.933370349805\n5916.0 23771.988339811218\n5788.5 23915.986674069958\n5607.0 24010.43503609106\n5440.5 24140.981676846197\n5106.0 24346.504164353137\n4513.5 24639.00499722376\n4393.5 24670.51749028318\n3963.0 24689.915047196002\n3768.0 24772.43253747918\n3585.0 24811.43087173792\n3451.5 24922.481676846197\n3361.5 24940.46085508051\n3255.0 25099.39755691283\n3067.5 25300.49250416435\n2955.0 25342.50333148251\n2709.0 25389.011660188786\n2497.5 25273.41837867851\n2344.5 25416.024153248196\n2290.5 25455.001665741256\n2310.0 22677.0" | \
  gdaltransform -tps -output_xy Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas__Priekules_apkart__6828__A1_L22_33_02.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas__Priekules_apkart__6828__A1_L22_33_02.vrt Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas__Priekules_apkart__6828__A1_L22_33_02.cut.vrt

# Generate tiles
gdal2tiles.py -r $GDAL_SAMPLING_TILE --xyz -z $GDAL_TILE_ZOOMS -x --processes=$GDAL_TILE_PROCESSES Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas__Priekules_apkart__6828__A1_L22_33_02.cut.vrt Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas__Priekules_apkart__6828__A1_L22_33_02.xyz