#!/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.924707341551 57.548957199051 -gcp 14817.0 0.0 27.071509787486 57.557955147499 -gcp 14817.0 9048.0 27.082207359453 57.510334688873 -gcp 0.0 9048.0 26.935404913518 57.501336740425 LVVA_F7404_US1_GV1772_0002.jpg LVVA_F7404_US1_GV1772_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"2364.75 4083.75\n2367.0 3728.25\n2356.5 3636.0\n2309.25 3534.75\n2189.25 3332.25\n2007.75 3288.75\n1804.5 3261.0\n1569.0 3246.75\n1476.75 3250.5\n1380.75 3276.75\n1427.25 3194.25\n1755.0 3182.25\n2133.0 2975.25\n2263.5 2909.25\n2419.5 2876.25\n2532.75 2814.75\n2637.0 2868.0\n2718.75 2944.5\n2763.0 3021.0\n2799.0 2988.0\n2844.75 2977.5\n2890.5 3018.75\n2925.75 3081.75\n2910.75 3147.0\n2921.25 3216.0\n2970.75 3299.25\n3112.5 3267.75\n3258.0 3260.25\n3376.5 3204.0\n3495.75 3175.5\n3498.75 2799.75\n3481.5 2745.75\n3399.75 2604.75\n3370.5 2481.0\n3408.75 2348.25\n3472.5 2310.75\n3534.75 2283.0\n3660.0 2282.25\n3708.0 2127.0\n3822.0 2116.5\n3851.25 2087.25\n3916.5 1878.75\n3918.75 1750.5\n3965.25 1635.75\n4027.5 1561.5\n4149.0 1588.5\n4299.75 1483.5\n4359.75 1474.5\n4413.0 1469.25\n4456.5 1307.25\n4493.25 1236.0\n4552.5 1175.25\n4594.5 984.75\n4581.75 838.5\n4560.0 765.0\n4521.75 723.75\n4482.0 742.5\n4427.25 650.25\n4360.5 554.25\n4531.5 537.0\n4677.75 478.5\n4907.25 564.75\n5061.0 565.5\n5087.25 642.0\n5189.25 705.75\n5256.75 784.5\n5334.0 765.75\n5400.0 804.0\n5512.5 1099.5\n5583.0 1417.5\n5625.0 1410.0\n5596.5 1503.0\n5608.5 1612.5\n5673.0 1854.0\n5749.5 1894.5\n5826.0 1944.0\n5898.0 2070.75\n5942.25 2263.5\n5993.25 2348.25\n6006.0 2484.0\n6036.75 2577.0\n6084.75 2656.5\n6288.75 2732.25\n6539.25 2774.25\n7317.0 3228.0\n7386.0 3522.0\n7360.5 3765.0\n7348.5 4002.0\n7458.75 4132.5\n7520.25 4111.5\n7590.0 4119.0\n7545.75 4282.5\n7513.5 4289.25\n7563.0 4522.5\n7653.75 4670.25\n7747.5 4812.0\n8240.25 5790.0\n8224.5 5865.0\n8177.25 5892.0\n8120.25 5914.5\n7644.0 5964.0\n7355.25 5949.0\n7099.5 5918.25\n7101.75 6255.0\n7086.75 6342.75\n7132.5 6356.25\n7323.75 6513.75\n7510.5 6633.75\n7592.25 6748.5\n7667.25 6880.5\n7582.5 6981.0\n7405.5 7087.5\n7310.25 7137.0\n7254.0 7203.0\n7197.0 7246.5\n7122.75 7274.25\n7054.5 7328.25\n7019.25 7413.75\n7047.75 7509.75\n7191.0 7642.5\n7246.5 7713.75\n7252.5 7755.0\n7251.75 7796.25\n7224.0 7869.0\n6921.0 8187.0\n5308.5 6892.5\n4233.0 6579.0\n4113.0 6568.5\n3397.5 6051.0\n3406.5 5940.0\n3126.0 5865.0\n3219.0 5599.5\n3244.5 5187.0\n3138.0 5173.5\n3079.5 5145.0\n3031.5 5089.5\n2983.5 4870.5\n2806.5 4848.0\n2632.5 4813.5\n2758.5 4487.25\n2663.25 4385.25\n2595.0 4242.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV1772_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV1772_0002.vrt LVVA_F7404_US1_GV1772_0002.cut.vrt

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