#!/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.799582284178 57.581903977432 -gcp 14922.0 0.0 26.939334869385 57.590239270552 -gcp 14922.0 13275.0 26.957273483276 57.521261882715 -gcp 0.0 13275.0 26.815309524536 57.512273239056 LVVA_F7404_US1_GV1783_0002.jpg LVVA_F7404_US1_GV1783_0002.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"4854.0 4905.0\n9933.0 4992.0\n9975.0 1080.0\n11370.0 1102.5\n11545.5 1291.5\n11635.5 2005.5\n11697.0 2077.5\n11709.0 2134.5\n11703.0 2218.5\n11764.5 2328.0\n11920.5 2409.0\n12003.0 2533.5\n12364.5 2607.0\n12387.0 2778.0\n12360.0 2994.0\n12594.0 3136.5\n12646.5 3112.5\n12661.5 3037.5\n12772.5 3007.5\n12906.0 3001.5\n13036.5 2950.5\n13191.0 2737.5\n13501.5 2884.5\n13438.5 2975.25\n13520.25 3060.0\n13458.0 3102.75\n13339.5 3169.5\n13291.5 3237.75\n13204.5 3341.25\n13171.5 3537.0\n13197.0 3819.0\n13614.75 3914.25\n13744.5 3924.75\n13838.25 3980.25\n14027.25 4159.5\n14089.5 4267.5\n14065.5 4545.75\n14033.25 4795.5\n14004.0 4948.5\n13950.75 5072.25\n13869.75 5183.25\n13708.5 5345.25\n13560.0 5439.0\n13500.75 5960.25\n13584.0 6019.5\n13654.5 6206.25\n13922.25 6336.0\n13848.75 6459.75\n13788.75 6723.75\n13776.75 6804.75\n13789.5 6861.75\n13767.0 6949.5\n13724.25 7064.25\n13690.5 7265.25\n13596.0 7374.0\n13476.75 7322.25\n13377.0 7342.5\n13349.25 7423.5\n13398.0 7521.75\n13401.0 7608.0\n13389.75 7655.25\n13410.0 7731.0\n13422.0 7800.75\n13379.25 7887.75\n13356.0 7951.5\n13377.75 8014.5\n13428.0 8049.0\n13384.5 8073.75\n13372.5 8110.5\n13419.0 8223.75\n13329.75 8253.75\n13224.75 8278.5\n13117.5 8298.0\n13053.75 8385.75\n12838.5 8496.0\n12786.75 8459.25\n12735.75 8419.5\n12684.0 8399.25\n12639.75 8409.75\n12541.5 8475.0\n12511.5 8511.75\n12498.0 8564.25\n12525.75 8726.25\n12512.25 8845.5\n12445.5 8988.75\n12393.0 9058.5\n12327.0 9093.75\n12261.75 9089.25\n12256.5 9148.5\n12264.75 9211.5\n12320.25 9338.25\n12216.75 9367.5\n12181.5 9405.75\n12169.5 9449.25\n12227.25 9580.5\n12023.25 9809.25\n11982.75 9815.25\n11955.0 9795.0\n11935.5 9765.75\n11916.75 9710.25\n11856.0 9694.5\n11817.0 9669.75\n11767.5 9679.5\n11707.5 9669.75\n11652.75 9717.75\n11623.5 9781.5\n11601.0 9870.0\n11581.5 9972.75\n11559.0 10069.5\n11550.0 10123.5\n11572.5 10227.0\n11611.5 10314.0\n11560.5 10389.0\n11421.0 10519.5\n11322.0 10602.0\n11244.0 10659.0\n11247.0 10722.0\n11296.5 10824.0\n11295.0 10950.0\n11179.5 11011.5\n11071.5 11025.0\n10894.5 11028.0\n10761.0 11052.0\n10699.5 11158.5\n10629.0 11182.5\n10563.0 11140.5\n10581.0 11086.5\n9664.5 11611.5\n9634.5 11916.0\n9360.0 11994.0\n9165.0 11890.5\n9226.5 11479.5\n9133.5 11179.5\n9201.0 11124.0\n9138.0 10915.5\n9072.0 10909.5\n9066.0 10960.5\n8581.5 10875.0\n7722.0 9774.0\n6372.0 9990.0\n5442.0 8748.0\n4852.5 7068.0\" | \
  gdaltransform -tps -output_xy LVVA_F7404_US1_GV1783_0002.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  LVVA_F7404_US1_GV1783_0002.vrt LVVA_F7404_US1_GV1783_0002.cut.vrt

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