#!/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 2658.0 2634.8340874811465 21.337216666666666 56.75 -gcp 19029.0 2617.713423831071 24.337216666666666 56.75 -gcp 19176.0 12429.207390648566 24.337216666666666 55.75 -gcp 2572.5 12539.073906485672 21.337216666666666 55.75 NLL241_Latvijas_satiksmes_ceļu_karte_1928-k_001_ktl1-3-124_2.jpg NLL241_Latvijas_satiksmes_ceļu_karte_1928-k_001_ktl1-3-124_2.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "624.0 13923.0\n708.0 8919.042232277527\n381.0 8219.971342383107\n417.0 4769.954751131221\n780.0 3752.933634992459\n810.0 1503.0\n20767.5 1448.9984917043741\n20880.0 13758.0" | \
  gdaltransform -tps -output_xy NLL241_Latvijas_satiksmes_ceļu_karte_1928-k_001_ktl1-3-124_2.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  NLL241_Latvijas_satiksmes_ceļu_karte_1928-k_001_ktl1-3-124_2.vrt NLL241_Latvijas_satiksmes_ceļu_karte_1928-k_001_ktl1-3-124_2.cut.vrt

# Generate tiles
gdal2tiles.py -r $GDAL_SAMPLING_TILE --xyz -z $GDAL_TILE_ZOOMS -x --processes=$GDAL_TILE_PROCESSES NLL241_Latvijas_satiksmes_ceļu_karte_1928-k_001_ktl1-3-124_2.cut.vrt NLL241_Latvijas_satiksmes_ceļu_karte_1928-k_001_ktl1-3-124_2.xyz