#!/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.85703071265 57.577831856425 -gcp 14913.0 0.0 26.999894035047 57.584194531587 -gcp 14913.0 9015.0 27.008406236017 57.536964238778 -gcp 0.0 9015.0 26.86554291362 57.530601563615 LVVA_F7404_US1_GV1786_0002.jpg LVVA_F7404_US1_GV1786_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"1629.0 4726.5\n1681.5 3700.5\n1882.5 2023.5\n1999.5 1863.0\n2086.5 1797.0\n2265.0 1800.0\n2361.0 1713.0\n2499.0 1746.0\n2625.0 1747.5\n2640.0 1785.0\n2778.0 1729.5\n2809.5 1531.5\n2769.0 1506.0\n2734.5 1476.0\n2764.5 1411.5\n2833.5 1444.5\n2883.0 1431.0\n2914.5 1444.5\n2958.0 1438.5\n3015.0 1311.0\n3082.5 1260.0\n3340.5 1249.5\n3604.5 1165.5\n3711.0 1045.5\n3786.0 1188.0\n3912.0 1089.0\n4026.0 1038.0\n3963.0 967.5\n3954.0 877.5\n3972.0 796.5\n4050.0 729.0\n4072.5 646.5\n4125.0 612.0\n4221.0 567.0\n4240.5 508.5\n4249.5 471.0\n4398.0 519.0\n4396.5 548.25\n4662.75 540.75\n4718.25 411.0\n4786.5 440.25\n4856.25 411.75\n4923.75 525.75\n4938.75 636.0\n5006.25 747.0\n5007.75 901.5\n5040.0 1042.5\n5076.0 1035.0\n5135.25 957.0\n5181.0 919.5\n5250.75 921.0\n5412.75 1026.0\n5508.0 1143.0\n5533.5 1223.25\n5552.25 1305.75\n5529.0 1671.0\n5564.25 1770.75\n5478.75 1840.5\n5397.0 1921.5\n5328.75 1966.5\n5253.75 1992.0\n5235.0 2039.25\n5288.25 2121.75\n5486.25 2253.0\n5556.75 2275.5\n5615.25 2329.5\n5709.0 2386.5\n5613.0 2511.75\n5565.0 2528.25\n5517.0 2573.25\n5454.0 2733.0\n5418.75 2840.25\n5459.25 2961.0\n5490.75 3013.5\n5484.0 3094.5\n5522.25 3120.0\n5556.0 3159.75\n5526.75 3255.0\n5542.5 3300.0\n5564.25 3340.5\n5640.0 3362.25\n5715.75 3409.5\n5745.0 3489.75\n5779.5 3529.5\n5817.75 3547.5\n5837.25 3583.5\n5824.5 3632.25\n5700.75 3666.75\n5525.25 3747.0\n5472.0 3737.25\n5450.25 3689.25\n5359.5 3939.75\n5268.0 4256.25\n5088.75 4365.0\n5010.0 4324.5\n4690.5 4684.5\n4681.5 4728.75\n4625.25 4788.0\n4524.75 4871.25\n4422.75 5305.5\n4365.75 5610.0\n4440.0 5642.25\n4458.0 5700.75\n4452.75 5723.25\n4542.0 5754.75\n4481.25 5797.5\n4403.25 5810.25\n5016.0 6873.0\n5184.0 6846.0\n5199.0 7420.5\n5239.5 7453.5\n5137.5 7518.0\n5081.25 7620.0\n5056.5 7851.0\n5076.75 8034.0\n4947.75 8229.0\n4908.0 8287.5\n4889.25 8326.5\n4720.5 7771.5\n4684.5 7710.0\n4664.25 7620.0\n4663.5 7563.75\n4673.25 7458.0\n4655.25 7391.25\n4629.75 7350.0\n4608.0 7359.0\n4584.0 7378.5\n4565.25 7404.75\n4532.25 7429.5\n4438.5 7476.75\n4415.25 7506.0\n4416.0 7617.0\n4375.5 7616.25\n4055.25 7620.0\n4005.0 7559.25\n3953.25 7509.0\n3850.5 7461.75\n3767.25 7403.25\n3735.75 7359.75\n3722.25 7287.0\n3732.75 7223.25\n3754.5 7158.0\n3743.25 7114.5\n3681.0 7065.75\n3624.75 7051.5\n3558.0 7056.75\n3481.5 7033.5\n3411.75 7096.5\n3280.5 7103.25\n3213.75 7120.5\n3179.25 7143.0\n2873.25 6459.0\n2832.0 6395.25\n2819.25 6312.75\n2781.0 6276.0\n2697.75 6257.25\n2686.5 6217.5\n2610.75 6184.5\n2551.5 6188.25\n2496.0 6201.0\n2473.5 6194.25\n2384.25 6021.0\n2429.25 5979.0\n2295.75 5874.75\n2119.5 5805.0\n2109.0 5835.0\n1894.5 5803.5\n1732.5 5760.0\n1571.25 5760.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV1786_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV1786_0002.vrt LVVA_F7404_US1_GV1786_0002.cut.vrt

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