#!/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 630.9599447616446 726.7916421614527 20.833854675292972 56.333288974155664 -gcp 5447.52473526557 711.0804587645223 21.499729156494144 56.333574460188444 -gcp 5479.378522147543 5062.613910979498 21.498184204101562 56.000300824577394 -gcp 647.8074533288913 5081.578978856264 20.83162307739258 56.000684786976336 2_izdevums/04_Jurmalciems.jpg 04_Jurmalciems.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "230.625 265.875\n5853.75 247.875\n5862.0 5523.0\n3076.5 5509.5\n286.5 5547.0" | \
  gdaltransform -tps -output_xy 04_Jurmalciems.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  04_Jurmalciems.vrt 04_Jurmalciems.cut.vrt

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