#!/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 29153.903225806487 11360.764112903242 27.43638 57.41957 -gcp 30060.71975806456 17352.568548387117 27.466659 57.342943 -gcp 20377.28225806455 33381.82862903227 27.261691 57.134477 -gcp 9592.034274193573 41764.675403225825 27.021462 57.020322 -gcp 25772.74596774197 41187.2661290323 27.39844354 57.03741414 LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0002.JPG LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "2064.0 42808.2639952858\n2004.0 36494.21331761933\n1764.0 29760.10606953447\n1656.0 22656.0\n1608.0 1571.8267530936946\n2754.0 1163.8621096051857\n9066.0 1326.0601060695344\n12312.0 1367.968179139658\n12624.0 1410.0070713022983\n14430.0 1422.0424278137891\n16578.0 1331.8939304655273\n17130.0 1380.03889216264\n19884.0 1416.045963464938\n29466.0 1368.0424278137891\n29520.0 1554.091926929876\n29796.0 1433.550972304066\n30030.0 2087.918680023571\n30087.0 3206.991160872127\n29896.5 3255.0\n29850.0 3324.5409546258106\n29994.0 3437.996464348851\n29784.0 3850.525633470831\n30216.0 9258.452563347084\n30336.0 18194.74366529169\n30342.0 31448.96287566293\n30393.0 42768.104301708896\n30339.0 43077.05833824396\n30228.0 43188.0\n2226.0 43176.07071302298" | \
  gdaltransform -tps -output_xy LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0002.vrt LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0002.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_0002.cut.vrt LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0002.xyz