#!/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 0.0 0.0 25.237516656411 57.440276385982 -gcp 7134.0 0.0 25.337419045833 57.449584947732 -gcp 7134.0 7050.0 25.35452507785 57.392285987206 -gcp 0.0 7050.0 25.254622688428 57.382977425456 LVVA_F7404_US1_GV242_0002.jpg LVVA_F7404_US1_GV242_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"1401.0 5950.5\n1491.0 5866.5\n1533.0 5752.5\n1569.0 5706.0\n1566.0 5581.5\n1575.0 5502.0\n1528.5 5457.0\n1521.0 5376.0\n1480.5 5323.5\n1468.5 5289.0\n1426.5 5253.0\n1429.5 5188.5\n1387.5 5133.0\n1320.0 5112.0\n1270.5 5014.5\n1266.0 4963.5\n1302.0 4884.0\n1296.0 4825.5\n1314.0 4786.5\n1339.5 4665.0\n1306.5 4452.0\n1201.5 4443.0\n1201.5 4293.0\n1066.5 4162.5\n868.5 4021.5\n898.5 3945.0\n828.0 3889.5\n900.0 3783.0\n894.0 3709.5\n831.0 3456.0\n721.5 3261.0\n720.0 3049.5\n642.0 2994.0\n645.0 2841.0\n679.5 2734.5\n777.0 2662.5\n891.0 2593.5\n1000.5 2364.0\n1075.5 2293.5\n1147.5 2184.0\n1249.5 2116.5\n1326.0 2047.5\n1380.0 1981.5\n1450.5 1954.5\n1486.5 1885.5\n1543.5 1839.0\n1647.0 1792.5\n1728.0 1734.0\n1884.0 1615.5\n1918.5 1558.5\n1977.0 1527.0\n2046.0 1366.5\n1987.5 1282.5\n2014.5 1168.5\n2118.0 1090.5\n2230.5 1114.5\n2446.5 1069.5\n2631.0 1012.5\n2799.0 952.5\n3073.5 909.0\n3186.0 916.5\n3252.0 814.5\n3484.5 670.5\n3564.0 627.0\n3627.0 567.0\n3735.0 567.0\n3865.5 490.5\n4042.5 403.5\n4225.5 367.5\n4521.0 229.5\n6645.0 240.0\n6645.0 4968.0\n3157.5 4945.5\n3138.0 6201.0\n1293.0 6207.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV242_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV242_0002.vrt LVVA_F7404_US1_GV242_0002.cut.vrt

# Generate tiles
gdal2tiles.py -r $GDAL_SAMPLING_TILE --xyz -z $GDAL_TILE_ZOOMS -x --processes=$GDAL_TILE_PROCESSES LVVA_F7404_US1_GV242_0002.cut.vrt LVVA_F7404_US1_GV242_0002.xyz