#!/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 13616.860887096766 8564.58870967742 27.10009575 57.17854853 -gcp 30102.368951612887 3738.9616935483905 27.47832 57.273741 -gcp 27816.39616935485 39306.42641129036 27.46461868 56.79448795 -gcp 11736.021169354839 34246.04939516133 27.08456039 56.84808099 -gcp 3350.334677419354 21382.12802419355 26.87805176 57.0063029 -gcp 30252.415471582004 22298.86232889613 27.50873566 57.01184276 -gcp 23229.30544354838 28233.420362903238 27.34788895 56.93299488 LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0004.JPG LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0004.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "2124.0 37523.92004441977\n2196.0 34606.45419211549\n2172.0 32137.03942254303\n2088.0 26412.739589117155\n2082.0 23399.000555247083\n2052.0 21011.19378123265\n2208.0 19738.953914491947\n2844.0 18342.266518600776\n4422.0 17664.049972237648\n5376.0 17177.97001665741\n5808.0 16801.09272626319\n6810.0 15678.116601887841\n7908.0 14129.916712937256\n10254.0 11315.75347029428\n10794.0 10446.023320377568\n11286.0 9491.890061077178\n11628.0 8886.059966685174\n12510.0 7937.930038867296\n15162.0 5669.75680177679\n16596.0 4608.186563020545\n19146.0 2939.673514714048\n20712.0 1919.7068295391455\n23208.0 5.953359244865169\n30786.0 6.003331482508656\n30756.0 23364.979455857858\n30984.0 40632.01332593004\n30996.0 41538.3098278734\n29574.0 41579.98334258745\n28476.0 41556.023320377564\n27636.0 41525.976679622436\n26538.0 41535.0\n25200.0 41514.01332593004\n24648.0 41595.02165463631\n24075.0 41619.12326485286\n21540.0 41688.759578012214\n14712.0 41713.8189894503\n11556.0 41748.01998889506\n8604.0 41701.232648528596\n2871.0 41532.01998889506\n2454.0 41601.03831204886\n1998.0 41568.03997779012" | \
  gdaltransform -tps -output_xy LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0004.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0004.vrt LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0004.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_0004.cut.vrt LVVA-006828-0002-000362-Celu kartes Sosejas un zemesceli_0004.xyz