#!/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.865606307983 57.528082554385 -gcp 14817.0 0.0 27.00792889808 57.536791378882 -gcp 14817.0 9249.0 27.017612457275 57.488338823997 -gcp 0.0 9249.0 26.875289867178 57.4796299995 LVVA_F7404_US1_GV1790_0002.jpg LVVA_F7404_US1_GV1790_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"2173.5 7815.0\n703.5 7779.0\n844.5 3234.0\n735.0 3117.0\n591.0 3010.5\n856.5 2853.0\n1313.25 2556.0\n1368.75 2525.25\n1743.75 2595.75\n2145.75 2608.5\n2640.0 1065.0\n3118.5 906.0\n3573.75 792.75\n3771.0 886.5\n3862.5 771.0\n4002.0 867.0\n3984.0 1038.0\n3979.5 1177.5\n4258.5 1425.0\n4396.5 1840.5\n4603.5 2268.0\n5370.0 2316.0\n5364.0 4098.0\n5844.0 4185.0\n6129.0 4485.0\n6297.0 4947.0\n6417.0 5430.0\n6292.5 5444.25\n6267.75 5517.75\n6375.0 5699.25\n6393.75 5801.25\n6387.0 5959.5\n6440.25 5991.75\n6459.0 6051.0\n6504.0 6023.25\n6561.0 6012.0\n6594.75 6045.75\n6604.5 6121.5\n6586.5 6159.0\n6595.5 6215.25\n6591.75 6276.0\n6552.0 6294.75\n6511.5 6346.5\n6463.5 6387.0\n6424.5 6385.5\n6397.5 6406.5\n6388.5 6452.25\n6364.5 6479.25\n6309.0 6459.75\n6278.25 6406.5\n6206.25 6431.25\n6171.0 6432.0\n6144.0 6406.5\n6106.5 6357.0\n6057.75 6393.0\n5991.75 6404.25\n5968.5 6351.0\n5935.5 6319.5\n5901.0 6228.75\n5898.75 6144.75\n5872.5 6119.25\n5846.25 6084.75\n5834.25 6024.0\n5819.25 5999.25\n5783.25 5983.5\n5708.25 5957.25\n5607.0 5934.0\n5533.5 5959.5\n5482.5 6042.0\n5449.5 6136.5\n5442.75 6242.25\n5496.75 6321.0\n5561.25 6391.5\n5580.0 6497.25\n5634.0 6554.25\n5646.75 6646.5\n5634.0 6739.5\n5616.75 6813.75\n5577.0 6873.75\n5503.5 6891.0\n5416.5 6956.25\n5340.0 6966.0\n5289.0 7232.25\n5252.25 7262.25\n5193.75 7271.25\n5125.5 7276.5\n5097.0 7308.0\n5082.0 7350.0\n5058.0 7442.25\n5074.5 7484.25\n5112.75 7504.5\n5128.5 7540.5\n5062.5 7576.5\n4989.75 7593.0\n4932.75 7638.75\n4849.5 7711.5\n4794.0 7730.25\n4755.75 7789.5\n4716.0 7802.25\n4677.75 7802.25\n4641.0 7846.5\n4600.5 7890.75\n4624.5 7937.25\n4680.75 8157.0\n4665.0 8192.25\n4621.5 8222.25\n4602.0 8252.25\n4605.75 8292.75\n4583.25 8317.5\n4497.0 8311.5\n4445.25 8325.75\n4330.5 8380.5\n4298.25 8443.5\n3780.0 8416.5\n3664.5 8101.5\n3397.5 8083.5\n3076.5 8086.5\n2790.0 8152.5\n2454.0 8178.0\n2205.0 8122.5\n2088.0 8028.0\n2145.0 7899.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV1790_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV1790_0002.vrt LVVA_F7404_US1_GV1790_0002.cut.vrt

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