#!/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.753339767456 57.577448263938 -gcp 14391.0 0.0 26.887922286987 57.589595226662 -gcp 14391.0 9174.0 26.900753974915 57.542549214785 -gcp 0.0 9174.0 26.766128540039 57.529464967466 LVVA_F7404_US1_GV1796_0002.jpg LVVA_F7404_US1_GV1796_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"8037.0 7252.5\n8164.5 7164.0\n8164.5 7071.0\n8239.5 6849.0\n8257.5 6666.0\n8305.5 6385.5\n8409.0 6217.5\n8500.5 6042.0\n8499.0 5943.0\n8466.0 5788.5\n8635.5 5716.5\n8721.0 5650.5\n8803.5 5544.0\n8905.5 5463.0\n8973.0 5362.5\n9064.5 5311.5\n9157.5 5184.0\n9373.5 5194.5\n9510.0 5127.0\n9634.5 5082.0\n9741.0 5077.5\n9880.5 4936.5\n9982.5 4897.5\n10057.5 4807.5\n10351.5 4374.0\n10324.5 3924.0\n10375.5 3768.0\n10380.0 3604.5\n9942.0 3192.0\n9729.0 3126.0\n9129.0 3000.0\n8919.0 2743.5\n8766.0 2611.5\n8710.5 2527.5\n8376.0 2667.0\n8308.5 2461.5\n8286.0 2310.0\n8269.5 2238.0\n8280.0 2179.5\n8283.0 2121.0\n8250.0 2070.0\n8344.5 1957.5\n8391.0 1789.5\n8625.0 1704.0\n8802.0 1714.5\n8892.0 1678.5\n8820.0 1545.0\n8947.5 1467.0\n9066.0 1557.0\n9166.5 1714.5\n9316.5 1780.5\n9450.0 1810.5\n9564.0 1854.0\n9787.5 1866.0\n9979.5 2040.0\n10146.0 2256.0\n10696.5 2064.0\n10651.5 1968.0\n10663.5 1854.0\n10786.5 1759.5\n10857.0 1639.5\n11122.5 1635.0\n11214.0 1618.5\n11293.5 1554.0\n11301.0 1414.5\n11365.5 1299.0\n11472.0 1215.0\n11506.5 1081.5\n11679.0 1126.5\n11715.0 1198.5\n11658.0 1219.5\n11629.5 1284.0\n11805.0 1314.0\n12076.5 1249.5\n12192.0 1389.0\n12298.5 1822.5\n12538.5 2194.5\n12858.0 2020.5\n12889.5 1939.5\n12942.0 1995.0\n13016.25 1947.75\n13047.0 1842.0\n13107.75 1771.5\n13155.0 1683.0\n13330.5 1711.5\n13479.75 1639.5\n13570.5 1519.5\n13626.0 1521.0\n13644.75 1626.0\n13784.25 1689.0\n13844.25 1872.0\n13875.0 1939.5\n13920.75 1949.25\n13946.25 2078.25\n13758.75 2223.0\n13719.75 2493.0\n13332.75 2622.75\n13269.0 2708.25\n13257.0 2814.75\n13249.5 2922.75\n13300.5 2964.0\n13366.5 3096.0\n13425.0 3213.0\n13430.25 3282.0\n13380.75 3486.75\n13330.5 3582.0\n13305.75 3621.0\n13284.75 3702.0\n13204.5 3672.0\n13012.5 3670.5\n12819.75 3711.0\n12729.75 3766.5\n12663.0 3864.75\n12600.75 3975.0\n12366.0 5443.5\n12277.5 5445.75\n12082.5 6220.5\n12267.0 6310.5\n12144.0 7416.0\n11946.0 7470.0\n11829.0 7551.0\n11700.0 7677.0\n11506.5 7774.5\n11074.5 7741.5\n10996.5 7704.0\n10947.0 7809.0\n10921.5 7843.5\n10840.5 7821.0\n10831.5 7858.5\n10905.0 7882.5\n11007.0 7873.5\n10908.0 7995.0\n10573.5 7845.0\n10564.5 7783.5\n10588.5 7753.5\n10677.0 7792.5\n10738.5 7786.5\n10837.5 7621.5\n10446.0 7428.0\n10110.0 7387.5\n9699.0 7345.5\n9663.0 7456.5\n9864.0 7825.5\n9804.0 7968.0\n9598.5 8076.0\n9238.5 7906.5\n8901.0 7437.0\n8964.0 7330.5\n9042.0 7291.5\n9070.5 7147.5\n9040.5 7062.0\n8946.0 7020.0\n8926.5 6945.0\n8809.5 6823.5\n8655.0 6766.5\n8614.5 6867.0\n8541.0 6961.5\n8595.0 7125.0\n8587.5 7195.5\n8568.0 7282.5\n8526.0 7342.5\n8326.5 7395.0\n8326.5 7462.5\n8292.0 7522.5\n8260.5 7498.5\n8227.5 7450.5\n8167.5 7465.5\n8106.0 7425.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV1796_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV1796_0002.vrt LVVA_F7404_US1_GV1796_0002.cut.vrt

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