#!/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 1620.0 35765.1545302946 23.837216666666666 56.75 -gcp 1060.5 2653.613396331295 25.337216666666666 56.75 -gcp 20724.0 2112.2743190661477 25.337216666666666 56.25 -gcp 21148.5 35871.29099499722 23.837216666666666 56.25 Latvijas_top.kartes_Rigas_aprinka_dienv.dala_un_Jaunjelgavas_6828_A1_L37_1_03.JPG Latvijas_top.kartes_Rigas_aprinka_dienv.dala_un_Jaunjelgavas_6828_A1_L37_1_03.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "11190.0 2345.9583564686286\n1599.0 2629.6524153248192\n1629.0 35762.91004997224\n11304.0 35831.963353692394\n21147.0 35870.99666851749\n21112.5 27511.542476402\n21003.0 18999.004997223765\n20901.0 10564.505830094391\n20722.5 2110.4441976679623" | \
  gdaltransform -tps -output_xy Latvijas_top.kartes_Rigas_aprinka_dienv.dala_un_Jaunjelgavas_6828_A1_L37_1_03.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  Latvijas_top.kartes_Rigas_aprinka_dienv.dala_un_Jaunjelgavas_6828_A1_L37_1_03.vrt Latvijas_top.kartes_Rigas_aprinka_dienv.dala_un_Jaunjelgavas_6828_A1_L37_1_03.cut.vrt

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