#!/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 3333.2963709677415 3602.655241935482 27.48109818 57.25108896 -gcp 26480.782258064508 4219.820564516129 28.03367615 57.25323642 -gcp 24470.262096774153 32615.10483870965 28.00243378 56.88416434 -gcp 1076.667338709618 38968.5 27.46427536 56.7949581 -gcp 8060.2935872565295 22436.153224116242 27.61323452 57.00870633 -gcp 24857.148786584305 20835.417409502283 28.00775528 57.03788122 LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0005.JPG LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0005.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "1533.0 20472.011660188786\n2028.0 11.92670738478622\n31062.0 11.993337034980566\n30960.0 2903.9000555247085\n31008.0 8952.093281510271\n30972.0 13067.920044419767\n30996.0 16044.999444752917\n31032.0 19044.12659633537\n31020.0 21660.03997779012\n31068.0 22913.18711826763\n31110.0 23729.79011660189\n31080.0 24485.98001110494\n31110.0 25308.273181565797\n31140.0 26981.103831204888\n31224.0 29489.983342587453\n31248.0 31920.10660744031\n31380.0 36023.900055524704\n31560.0 40526.032204330935\n26454.0 40764.02332037757\n26433.0 41010.00333148251\n22704.0 41220.219877845644\n16824.0 41427.33148250972\n10620.0 41700.07995558024\n7308.0 41695.219322598554\n5316.0 41550.00666296502\n3552.0 41382.26985008329\n2106.0 41322.19655746808\n858.0 41196.146585230425\n912.0 40889.99666851749" | \
  gdaltransform -tps -output_xy LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0005.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0005.vrt LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0005.cut.vrt

# Generate tiles
gdal2tiles.py -r $GDAL_SAMPLING_TILE --xyz -z $GDAL_TILE_ZOOMS -x --processes=$GDAL_TILE_PROCESSES LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0005.cut.vrt LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0005.xyz