#!/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.814880371094 57.548076358567 -gcp 14574.0 0.0 26.954441070557 57.555720858802 -gcp 14574.0 9090.0 26.965856552124 57.510152518178 -gcp 0.0 9090.0 26.826210021973 57.502406229479 LVVA_F7404_US1_GV1789_0002.jpg LVVA_F7404_US1_GV1789_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"2445.75 3961.5\n2496.75 3854.25\n2463.0 3777.0\n2469.75 3744.75\n2547.0 3637.5\n2445.75 3523.5\n2427.0 3465.0\n2414.25 3302.25\n2324.25 3084.75\n2248.5 2841.75\n2248.5 2748.0\n2284.5 2700.75\n2355.75 2652.75\n2413.5 2657.25\n2478.0 2616.75\n2502.0 2559.75\n2514.0 2496.0\n2646.0 2389.5\n2723.25 2305.5\n2675.25 2131.5\n2698.5 2076.0\n2768.25 2032.5\n2793.75 1976.25\n2789.25 1937.25\n2761.5 1887.0\n2848.5 1824.75\n2859.75 1878.0\n2996.25 1911.75\n3110.25 1928.25\n3220.5 1909.5\n3309.0 1670.25\n3361.5 1647.75\n3376.5 1608.0\n3419.25 1542.75\n3426.75 1457.25\n3409.5 1413.0\n3378.0 1404.0\n3372.75 1349.25\n3401.25 1050.0\n3544.5 1058.25\n3584.25 1041.0\n3635.25 1039.5\n3697.5 1076.25\n3745.5 1134.0\n3867.0 1212.75\n3995.25 1125.0\n4098.75 1036.5\n4137.75 965.25\n4247.25 840.0\n4359.0 744.0\n4431.0 794.25\n4473.75 847.5\n4509.75 872.25\n4535.25 845.25\n4570.5 840.75\n4596.75 847.5\n4674.0 888.0\n4789.5 906.0\n4896.75 947.25\n4980.0 944.25\n5037.0 959.25\n5081.25 984.75\n5134.5 1066.5\n5121.75 1140.0\n5007.0 1264.5\n4939.5 1290.0\n4872.0 1403.25\n4809.0 1543.5\n4888.5 1595.25\n4850.25 1658.25\n4844.25 1737.375\n4803.375 1768.875\n4771.125 1814.25\n4779.0 1913.625\n4757.25 1974.75\n4749.375 2019.375\n4762.125 2071.5\n4754.25 2120.625\n4783.875 2164.125\n4824.75 2180.25\n4892.625 2191.125\n4962.0 2256.0\n4977.375 2301.0\n4959.75 2397.0\n5042.625 2439.375\n5125.5 2490.0\n5466.0 3435.0\n6993.0 3206.25\n7131.0 3541.5\n7336.5 3746.25\n7522.5 3969.0\n7547.25 4095.0\n7499.25 4118.25\n7528.5 4201.5\n7660.5 4236.0\n7662.75 4338.75\n7915.5 4383.75\n8181.0 4467.0\n8274.75 4774.5\n8313.0 4835.25\n8325.75 4924.5\n8247.0 5295.75\n7284.0 5553.0\n6832.5 7142.25\n6552.75 7152.75\n6334.5 7139.25\n6216.75 7106.25\n6098.25 7050.0\n5313.75 7611.0\n5286.0 7692.75\n5184.0 7761.0\n5137.5 7728.75\n4951.5 7630.5\n4779.75 7639.5\n4716.75 7636.5\n4647.75 7602.0\n4587.0 7551.0\n4490.25 7427.25\n4434.75 7311.75\n4422.75 6993.0\n4392.0 6827.25\n4305.0 6604.5\n4251.75 6468.75\n4235.25 6382.5\n4323.0 6160.5\n4425.75 6000.75\n4668.75 5243.25\n4628.25 4851.0\n4507.5 4845.75\n4412.25 4902.75\n4312.5 4950.75\n4218.0 4944.0\n4119.0 4992.75\n4098.0 5016.75\n4092.0 5057.25\n4067.25 5068.5\n4041.0 5060.25\n3999.0 5058.75\n3974.25 5034.75\n3885.0 5042.25\n3736.5 4964.25\n3711.75 4964.25\n3686.25 4965.75\n3645.0 5004.0\n3583.5 5005.5\n3547.5 4983.75\n3511.5 4887.0\n3413.25 4834.5\n3381.0 4839.75\n3339.75 4868.25\n3262.5 4909.5\n3190.5 4905.75\n3120.0 4862.25\n3111.0 4816.5\n3026.25 4751.25\n2883.75 4620.75\n2856.75 4566.0\n2811.0 4501.5\n2672.25 4344.75\n2666.25 4294.5\n2652.0 4270.5\n2607.75 4244.25\n2535.0 4211.25\n2401.5 4080.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV1789_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV1789_0002.vrt LVVA_F7404_US1_GV1789_0002.cut.vrt

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