#!/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 26.908650398254 57.523842827432 -gcp 15060.0 0.0 27.056309358889 57.533350660808 -gcp 15060.0 10632.0 27.069794882424 57.476984707384 -gcp 0.0 10632.0 26.921917519995 57.467685748974 LVVA_F7404_US1_GV1770_0002.jpg LVVA_F7404_US1_GV1770_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"2296.5 4989.0\n2253.0 4921.5\n2259.0 4747.5\n2175.0 4761.0\n2133.0 4719.0\n2071.5 4600.5\n1938.0 4186.5\n1953.0 4080.0\n1921.5 3972.0\n1552.5 3798.0\n943.5 3594.0\n934.5 1719.0\n936.75 688.5\n1593.0 699.75\n2076.0 711.0\n2102.25 556.5\n2103.75 507.0\n2040.75 447.75\n2014.5 425.25\n2019.0 392.25\n2058.0 361.5\n2091.75 326.25\n2084.25 285.75\n2073.75 256.5\n2055.0 228.0\n2049.75 140.25\n2081.25 96.0\n2096.25 111.0\n2106.0 201.0\n2140.5 256.5\n2199.0 348.0\n2187.0 508.5\n2202.0 567.0\n2686.5 697.5\n3129.0 834.0\n3195.0 829.5\n3328.5 882.0\n3517.5 870.0\n3696.0 888.0\n3736.5 943.5\n3792.0 988.5\n3852.0 1027.5\n3936.0 1047.0\n3999.0 1084.5\n4057.5 1141.5\n4140.0 1272.0\n4291.5 1294.5\n4423.5 1353.0\n4558.5 1435.5\n4645.5 1452.0\n4780.5 1462.5\n4860.0 1491.0\n4867.5 1681.5\n4935.0 1758.0\n5028.0 1791.0\n5112.0 1686.0\n5212.5 1518.0\n5328.0 1332.0\n5640.0 1414.5\n5646.0 1519.5\n5683.5 1573.5\n5796.0 1629.0\n5857.5 1701.0\n6039.0 1819.5\n6319.5 2011.5\n6507.0 2062.5\n7534.5 2376.0\n9163.5 3693.0\n9022.5 3930.0\n9079.5 4168.5\n8941.5 4369.5\n8796.75 4476.75\n8488.5 4579.5\n8220.75 4524.0\n8029.5 4542.0\n7681.5 4611.0\n7614.0 4744.5\n7375.5 4614.0\n7011.0 4714.5\n6711.0 4737.0\n6552.0 4711.5\n6502.5 4759.5\n6184.5 4720.5\n6244.5 4935.0\n5850.0 5040.0\n5643.0 4932.0\n5430.0 4917.0\n4909.5 4888.5\n4812.0 4899.0\n4719.0 4921.5\n4248.0 5175.0\n4071.0 5164.5\n3918.0 5220.0\n3720.0 5154.0\n3598.5 5125.5\n3483.0 5154.0\n3391.5 5202.0\n3289.5 5197.5\n3168.0 5169.0\n3055.5 5113.5\n2907.0 5152.5\n2878.5 5092.5\n2805.0 5052.0\n2722.5 5001.0\n2644.5 4959.0\n2592.0 4929.0\n2401.5 4968.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV1770_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV1770_0002.vrt LVVA_F7404_US1_GV1770_0002.cut.vrt

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