#!/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 1959.2518175809648 2719.698319327731 24.5 57.5 -gcp 16080.584930601453 2706.3008403361346 25 57.5 -gcp 16155.499669530735 16111.846218487393 25 57.25 -gcp 1935.358889623265 16102.864705882352 24.5 57.25 Latvijas_dalas_Lietuvas_dalas_un_Igau_T82JyKV.kartes_Straupe_6828_A1_L25_44_02.JPG Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Straupe_6828_A1_L25_44_02.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"1959.2518175809648 2719.698319327731\n15935.90482485129 2709.001680672269\n15791.950429610046 2825.9764705882353\n15801.001982815598 3026.963025210084\n15662.869134170522 3477.1109243697483\n15774.033707865168 3641.9882352941177\n15785.857237276932 3851.673949579832\n15782.984137475214 4341.021848739496\n15830.93654990086 4754.862184873949\n15981.019828155982 4952.719327731093\n15908.92267019167 5009.695798319328\n15993.001982815598 5189.781512605042\n15779.833443489755 5363.963025210084\n16068.023793787179 5670.0\n16064.996034368804 5816.946218487395\n15890.88301387971 5933.971428571429\n16070.994051553205 6153.031932773109\n16112.996034368804 10116.396638655462\n16155.499669530735 16111.846218487393\n1935.358889623265 16102.864705882352\" | \
  gdaltransform -tps -output_xy Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Straupe_6828_A1_L25_44_02.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Straupe_6828_A1_L25_44_02.vrt Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Straupe_6828_A1_L25_44_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_Straupe_6828_A1_L25_44_02.cut.vrt Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Straupe_6828_A1_L25_44_02.xyz