#!/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 3070.5 2342.312117843246 24.837216666666666 57 -gcp 24469.5 2451.246803779878 25.837216666666666 57 -gcp 24658.5 32835.3902167871 25.837216666666666 56.25 -gcp 2766.0 32652.970122290164 24.837216666666666 56.25 Latvijas_top.kartes_Jaunjelgavas_un_Jekabpils_aprinkis_6828_A1_L37_2_03.JPG Latvijas_top.kartes_Jaunjelgavas_un_Jekabpils_aprinkis_6828_A1_L37_2_03.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "3069.0 2343.0\n2944.5 17498.25\n2765.25 32652.750416435312\n24661.5 32835.0\n24511.5 2447.9983324068926" | \
  gdaltransform -tps -output_xy Latvijas_top.kartes_Jaunjelgavas_un_Jekabpils_aprinkis_6828_A1_L37_2_03.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  Latvijas_top.kartes_Jaunjelgavas_un_Jekabpils_aprinkis_6828_A1_L37_2_03.vrt Latvijas_top.kartes_Jaunjelgavas_un_Jekabpils_aprinkis_6828_A1_L37_2_03.cut.vrt

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