#!/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 5701.5 2085.244302390217 22.337216666666666 55.75 -gcp 20482.5 2439.30683713174 23.00388333333333 55.75 -gcp 20071.5 32816.4334352418 23.00388333333333 55 -gcp 4991.578652322826 32473.63243620037 22.337216666666666 55 Lietuv.un_Latv.dienvidu_dalas_top.kartes_Varni-Kelme-Baisogala_6828_A1_L35_2_03.JPG Lietuv.un_Latv.dienvidu_dalas_top.kartes_Varni-Kelme-Baisogala_6828_A1_L35_2_03.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "14236.5 2317.5224875069407\n5700.75 2083.5008328706276\n5421.75 12473.18711826763\n5195.25 23085.747084952804\n4989.0 32471.259578012214\n13599.0 32679.00749583565\n22216.5 32860.49750138812\n22743.0 2466.0049972237643" | \
  gdaltransform -tps -output_xy Lietuv.un_Latv.dienvidu_dalas_top.kartes_Varni-Kelme-Baisogala_6828_A1_L35_2_03.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  Lietuv.un_Latv.dienvidu_dalas_top.kartes_Varni-Kelme-Baisogala_6828_A1_L35_2_03.vrt Lietuv.un_Latv.dienvidu_dalas_top.kartes_Varni-Kelme-Baisogala_6828_A1_L35_2_03.cut.vrt

# Generate tiles
gdal2tiles.py -r $GDAL_SAMPLING_TILE --xyz -z $GDAL_TILE_ZOOMS -x --processes=$GDAL_TILE_PROCESSES Lietuv.un_Latv.dienvidu_dalas_top.kartes_Varni-Kelme-Baisogala_6828_A1_L35_2_03.cut.vrt Lietuv.un_Latv.dienvidu_dalas_top.kartes_Varni-Kelme-Baisogala_6828_A1_L35_2_03.xyz