#!/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 688.0556061987236 684.7707390648568 20.833333333333332 57.333333333333336 -gcp 5368.188582497721 707.2801659125189 21.5 57.333333333333336 -gcp 5372.396308113035 5056.69419306184 21.5 57 -gcp 682.5512762078395 5044.655731523379 20.833333333333332 57 2_izdevums/01_Uzava.jpg 01_Uzava.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "233.25 221.25113122171948\n244.49683766690092 5126.247917823432\n3164.25158116655 5103.752082176568\n3168.001054111033 5023.498334258745\n5375.243148278285 5058.763742365352\n5366.0618411806045 235.1256246529706" | \
  gdaltransform -tps -output_xy 01_Uzava.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  01_Uzava.vrt 01_Uzava.cut.vrt

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