#!/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.013050079346 56.776291018255 -gcp 8805.0 0.0 26.090126037598 56.754510994511 -gcp 8805.0 13050.0 26.040000915527 56.695410863885 -gcp 0.0 13050.0 25.96489906311 56.717884566761 LVVA_F7404_US1_GV2032_0002.jpg LVVA_F7404_US1_GV2032_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"1332.0 10062.0\n1243.5 9898.5\n1120.5 9625.5\n1039.5 9391.5\n1279.5 9304.5\n1704.0 9294.0\n1749.0 9223.5\n1764.0 9133.5\n1762.5 8989.5\n1840.5 8856.0\n1857.0 8788.5\n1831.5 8686.5\n1807.5 8605.5\n1806.0 8287.5\n1734.0 8239.5\n1738.5 8130.0\n1806.0 8062.5\n1915.5 7963.5\n2023.5 7939.5\n2410.5 7692.0\n2436.0 7602.0\n2451.0 7456.5\n2419.5 7287.0\n2388.0 7161.0\n2341.5 7036.5\n2319.0 6720.0\n2337.0 6621.0\n2454.0 6585.0\n2493.0 6504.0\n2494.5 6433.5\n2457.0 6385.5\n2385.0 6367.5\n1689.0 6043.5\n1296.0 5721.0\n907.5 5325.0\n924.0 4425.0\n916.5 3970.5\n937.5 3565.5\n1062.0 3379.5\n1288.5 3271.5\n1518.0 3162.0\n1996.5 3115.5\n1996.5 2883.0\n1968.0 2776.5\n1980.0 2679.0\n1929.0 2635.5\n1900.5 2428.5\n1854.0 2380.5\n1854.375 2365.875\n1876.125 2354.25\n1883.625 2318.625\n1858.5 2268.0\n1876.875 2250.75\n1918.125 2265.0\n1917.375 2215.5\n1842.75 2128.125\n1815.75 2133.375\n1809.375 2111.625\n1840.125 1975.125\n1759.5 1975.875\n1743.0 1941.0\n1741.5 1910.25\n1758.0 1893.0\n1744.5 1865.25\n1710.75 1849.5\n1700.25 1792.5\n1705.5 1749.75\n1742.25 1733.25\n1746.0 1710.0\n1710.0 1686.0\n1713.75 1650.0\n1677.0 1623.75\n1674.75 1575.75\n1650.0 1522.5\n1670.25 1484.25\n1663.5 1449.0\n1608.75 1351.5\n1563.0 1364.25\n1533.75 1186.5\n1550.25 1104.0\n1604.25 1006.5\n1689.75 958.5\n1717.5 839.25\n1866.0 699.0\n2805.0 1609.5\n2977.5 1710.0\n3021.75 1668.0\n3095.25 1668.0\n3159.0 1689.0\n3202.5 1758.0\n3198.75 1865.25\n3370.5 1999.5\n3376.5 2088.0\n3465.0 2285.25\n3519.75 2307.75\n3570.75 2406.75\n3585.0 2537.25\n3618.0 2600.25\n3609.0 2676.75\n3632.25 2734.5\n3605.25 2794.5\n3608.25 2836.5\n3588.75 2884.5\n3585.75 2940.0\n3560.25 2973.0\n3566.25 3029.25\n3558.0 3080.25\n3520.5 3121.5\n3520.5 3160.5\n3547.5 3199.5\n3552.0 3244.5\n3565.5 3266.25\n3572.25 3316.5\n3563.25 3372.0\n3580.5 3398.25\n3595.5 3417.75\n3602.25 3441.0\n3599.25 3490.5\n3579.75 3528.0\n3588.75 3555.0\n3607.5 3565.5\n3624.75 3604.5\n3614.25 3660.75\n3621.75 3697.5\n3653.25 3738.75\n3666.75 3798.75\n3693.75 3828.0\n3728.25 3852.75\n3723.0 3911.25\n3729.0 3960.75\n3759.75 3978.75\n3755.25 4029.0\n3774.0 4064.25\n3810.75 4097.25\n3807.75 4145.25\n3839.25 4217.25\n3896.25 4383.0\n3899.25 4437.75\n3927.0 4514.25\n3990.0 4587.0\n4209.0 4710.0\n4258.5 4650.0\n4293.0 4581.75\n4381.5 4422.75\n4546.5 4273.5\n4515.0 4209.0\n4371.0 4170.0\n4330.5 4110.0\n4327.5 3996.0\n4358.25 3925.5\n4419.0 3889.5\n4659.0 3949.5\n4728.0 3910.5\n4887.0 3927.0\n4930.5 3991.5\n4990.5 4033.5\n5029.5 4092.0\n6117.0 3913.5\n6765.0 3972.0\n6879.0 4279.5\n6957.0 4453.5\n6975.0 4584.0\n6853.5 5029.5\n6922.5 5422.5\n6876.0 5785.5\n6859.5 5950.5\n6870.0 6057.0\n7110.0 6454.5\n7116.0 6586.5\n7330.5 6714.0\n7393.5 6685.5\n7426.5 6726.0\n7437.0 6772.5\n7492.5 6820.5\n7552.5 6826.5\n7626.0 6799.5\n7710.0 6862.5\n7830.0 6979.5\n7855.5 7095.0\n7944.0 7303.5\n8002.5 7474.5\n8082.0 7561.5\n8236.5 7666.5\n8221.5 7921.5\n8223.0 8182.5\n8179.5 8283.0\n8196.0 8368.5\n8172.0 8436.0\n8199.0 8509.5\n8185.5 8592.0\n8220.0 8671.5\n8197.5 8779.5\n8176.5 8937.0\n8265.0 9277.5\n8137.5 9277.5\n7893.0 9174.0\n7702.5 9147.0\n7251.0 8761.5\n6643.5 9222.0\n5266.5 9501.0\n4537.5 9658.5\n4588.5 9858.0\n4679.25 10038.0\n4791.0 10162.5\n5207.25 10569.75\n5278.5 11317.5\n5178.0 11296.5\n4843.5 11292.0\n4626.0 11313.0\n3954.0 11218.5\n2994.0 11110.5\n2802.0 11115.0\n2563.5 11184.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV2032_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV2032_0002.vrt LVVA_F7404_US1_GV2032_0002.cut.vrt

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