#!/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 27.011003494263 57.461157244289 -gcp 14763.0 0.0 27.156143188477 57.469281532634 -gcp 14763.0 9141.0 27.167859077454 57.421779658876 -gcp 0.0 9141.0 27.017869949341 57.413205665425 LVVA_F7404_US1_GV1767_0002.jpg LVVA_F7404_US1_GV1767_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"1200.0 5854.5\n1147.5 5791.5\n1089.0 5688.0\n1087.5 5607.0\n1107.0 5445.0\n1056.0 5362.5\n940.5 5235.0\n645.0 4935.0\n693.0 4558.5\n849.0 4612.5\n942.0 4623.0\n1020.75 4611.0\n1080.0 4657.5\n1174.5 4610.25\n1218.75 4495.5\n1291.5 4452.75\n1314.75 4429.5\n1310.25 4379.25\n1341.75 4348.5\n1365.75 4317.0\n1209.75 4219.5\n1047.75 4106.25\n873.0 3925.5\n714.0 3755.25\n707.625 3729.375\n838.5 3663.0\n901.5 3677.625\n942.0 3712.125\n1011.0 3727.875\n1071.75 3712.5\n1102.5 3632.25\n1149.375 3571.125\n1193.25 3523.5\n1253.625 3495.375\n1369.5 3463.5\n1486.875 3479.25\n1586.25 3500.25\n1642.5 3538.875\n1681.5 3595.5\n1897.5 3594.0\n1941.75 3493.5\n1956.0 3337.5\n1950.0 3148.5\n1988.25 3060.75\n2048.25 2985.75\n2140.5 2913.75\n2241.75 2903.25\n2367.0 2922.0\n2372.25 2905.5\n2465.25 2914.5\n2536.5 2903.25\n2646.0 2958.75\n2696.25 3013.5\n2708.25 3072.75\n2781.75 3136.5\n2826.375 3203.625\n2829.0 3322.125\n2814.0 3454.875\n2845.5 3592.5\n2891.25 3637.875\n2943.0 3689.625\n3027.75 3679.5\n3103.125 3637.125\n3156.75 3613.5\n3251.25 3618.75\n3361.5 3614.25\n3392.25 3654.375\n3413.25 3700.5\n3392.25 3757.5\n3307.5 3873.75\n3165.0 4023.0\n3157.5 4126.5\n3243.0 4194.0\n3333.75 4213.5\n3430.5 4195.5\n3518.25 4146.75\n3615.75 4099.5\n3666.0 4053.75\n3720.0 4023.75\n3774.0 3999.0\n3807.375 3937.5\n3872.25 4068.75\n3935.25 4232.25\n3925.5 4350.0\n3924.75 4433.25\n3929.25 4494.0\n3947.25 4545.0\n4004.25 4631.25\n4011.0 4718.25\n4030.5 4791.0\n4071.75 4864.5\n4119.0 4923.0\n4251.0 5142.75\n4296.0 5185.5\n4624.5 5313.75\n4717.5 5287.5\n4815.75 5262.0\n5011.5 5223.0\n5088.0 5237.25\n5145.0 5281.5\n5185.5 5502.0\n5249.25 5697.0\n5314.5 5970.75\n5320.5 6121.5\n5290.5 6226.5\n5224.5 6276.0\n5139.0 6367.5\n5053.5 6420.75\n4965.0 6457.5\n4792.5 7189.5\n4660.5 7194.375\n4488.0 7231.125\n4356.75 7263.0\n4288.5 7302.75\n4260.0 7348.5\n4213.5 7363.5\n4177.5 7317.0\n4168.5 7260.0\n4146.0 7235.25\n4123.5 7224.75\n4112.625 7342.5\n4080.0 7411.125\n3986.625 7455.375\n3918.0 7457.25\n3849.375 7427.25\n3706.5 7299.75\n3628.5 7084.5\n3575.25 7017.75\n3409.5 6940.5\n3202.5 6854.25\n3144.0 6849.0\n3072.0 6852.75\n2991.75 6852.0\n2891.25 6822.0\n2811.0 6762.0\n2765.25 6694.5\n2756.25 6660.75\n2760.75 6609.0\n2810.25 6470.25\n2740.5 6414.0\n2665.5 6433.5\n2605.5 6463.5\n2543.25 6466.5\n2475.75 6445.5\n2423.25 6402.0\n2400.75 6350.25\n2361.75 6315.0\n2306.25 6282.75\n2224.5 6250.5\n2161.5 6143.25\n2116.5 5987.25\n2089.5 5852.25\n1771.5 5950.5\n1486.5 5964.0\n1228.5 5988.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV1767_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV1767_0002.vrt LVVA_F7404_US1_GV1767_0002.cut.vrt

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