#!/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 798.4857898215465 766.8840336134455 24.5 56.5 -gcp 15284.426966292136 807.90756302521 25 56.5 -gcp 15342.38664904164 14260.658823529413 25 56.25 -gcp 753.4719101123595 14215.851260504201 24.5 56.25 Latvijas_dalas_Lietuvas_dalas_un__A0OvRdU.kartes_Skaistkalne_6828_A1_L25_48_02.JPG Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Skaistkalne_6828_A1_L25_48_02.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"1029.0178453403832 13344.015126050419\n1002.0654329147388 12413.741176470588\n1011.0158625247852 12147.028571428571\n972.0456047587575 11777.89075630252\n891.001982815598 11564.949579831933\n803.9920687376074 11532.025210084033\n798.027759418374 4402.681512605042\n837.025776602776 2279.828571428571\n826.5183410442828 768.0004201680671\n8055.0 840.028571428571\n9017.456708526108 845.9697478991593\n9098.191011235955 893.9949579831937\n11889.172504957038 932.9848739495797\n12659.387309980171 945.042016806723\n13199.488433575676 902.954621848739\n15291.003965631196 1023.1210084033614\n15342.38664904164 14260.658823529413\n999.1625908790481 14220.00168067227\" | \
  gdaltransform -tps -output_xy Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Skaistkalne_6828_A1_L25_48_02.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Skaistkalne_6828_A1_L25_48_02.vrt Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Skaistkalne_6828_A1_L25_48_02.cut.vrt

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