#!/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 2895.979838727791 10820.271169338313 27.44127358 57.43275547 -gcp 16092.148185500557 8902.988407243687 27.76107534 57.46415447 -gcp 21416.5303801001 40154.18149879336 27.91471957 57.0664005 -gcp 9479.864844031661 40432.38394658849 27.63430738 57.05800367 -gcp 3452.5645161290395 15342.99495967742 27.46607796 57.37532195 -gcp 3604.0516861281276 25116.13274861867 27.48135688 57.25133534 -gcp 11761.578629032256 30944.877520161328 27.68262821 57.1789205 -gcp 25518.315905691074 18785.312372606244 27.99730287 57.34228224 LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0003.JPG LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0003.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "924.0 32435.989393046555\n939.0 22206.0\n990.0 17076.00353565115\n1080.0 11934.0\n1152.0 10644.72626318712\n1464.0 1963.8822876179897\n1458.0 1697.916712937257\n2706.0 1727.9200444197668\n4470.0 1841.9900055524708\n5742.0 1848.0133259300387\n7038.0 1781.9966685174904\n11160.0 1703.586896168795\n14262.0 1752.0199888950583\n17760.0 1668.00666296502\n19920.0 1606.8872848417545\n26280.0 1450.1476957245973\n28536.0 1428.033314825097\n30870.0 1157.9233759022766\n30888.0 4056.0399777901166\n30960.0 6924.059966685175\n31152.0 17690.118822876182\n31224.0 23996.73514714048\n31302.0 31205.720155469186\n31374.0 34104.06329816769\n31404.0 36816.026651860084\n31548.0 42689.84089569829\n31488.0 42929.911608721275\n1680.0 42965.91868002357\n870.0 42707.96464348851" | \
  gdaltransform -tps -output_xy LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0003.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0003.vrt LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0003.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_0003.cut.vrt LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0003.xyz