#!/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 507.0 1338.0 24 58 -gcp 10920.0 1338.0 25.5 58 -gcp 11106.0 14448.0 25.5 57 -gcp 420.0 14496.0 24 57 33_Cesis_1940.jpg 33_Cesis_1940.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"10675.5 1022.9625208217657\n10362.0 841.5566352026652\n10084.5 811.6782343142697\n9921.0 538.5233203775681\n9792.0 367.49167129372563\n9597.0 139.5241532481956\n9330.0 152.998334258745\n9297.0 321.07495835646864\n9204.75 329.97688784008886\n9174.0 525.8328706274292\n9197.25 657.0266518600778\n9423.75 903.127012770683\n9417.75 1094.2937257079402\n9382.5 1117.5037479178234\n9417.75 1239.8149639089395\n9380.25 1250.9833425874515\n9399.75 1322.2724875069407\n9351.0 1323.0603831204887\n9326.25 1188.732509716824\n9261.0 1199.9958356468628\n9244.5 1174.4950027762352\n9285.0 1127.3416157690172\n9276.75 1020.8124652970571\n9191.25 913.510410882843\n9108.0 885.7533314825097\n9099.75 855.7549972237647\n9150.0 789.7458356468628\n9166.5 691.594947251527\n9156.75 518.2491671293725\n9088.5 513.7549972237646\n9074.25 483.7554136590783\n9136.5 422.33495280399785\n9070.5 217.58911715713487\n8788.5 17.24000555247096\n8638.5 217.56996113270395\n8572.5 213.00749583564686\n8524.5 179.998334258745\n8329.5 147.01415880066622\n8028.0 208.51249305941155\n8026.5 336.0841199333702\n7902.0 334.40338700721827\n7521.0 982.4192115491394\n7391.25 1060.4825097168239\n7097.25 1055.902137701277\n6640.5 1335.7554136590784\n6198.75 1336.51790671849\n6158.25 1282.473348139922\n5943.0 1335.0016657412548\n507.0 1341.0\n423.0 14496.0\n11109.0 14448.0\n10923.0 1338.0\" | \
  gdaltransform -tps -output_xy 33_Cesis_1940.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  33_Cesis_1940.vrt 33_Cesis_1940.cut.vrt

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