#!/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 1287.0 279.16407551360356 -179.99411787 67.99990061 -gcp 13116.75 270.9328151027207 -167.99700394 67.99968371 -gcp 14110.235237937282 10797.825478050963 -167.99743929 63.99963258 -gcp 270.0 10812.952803997778 -179.99748631 63.999843 -gcp 7200.75 560.2833148250971 -173.99701509 67.99979269 -gcp 7189.875 11128.269711271516 -173.99744883 63.99973828 -gcp 3261.0 441.2081447963801 -177.98227924 67.93293005 -gcp 5230.5 532.7541478129713 -175.9875385 67.89489937 -gcp 9180.0 528.4570135746607 -171.99858734 67.89544469 -gcp 11148.0 438.4117647058823 -170.00486996 67.93166457 -gcp 13366.5 2904.22850678733 -167.99592878 66.93984604 -gcp 13611.0 5536.846907993967 -167.99978757 65.92057158 -gcp 13857.75 8165.429487179487 -168.00098824 64.94190638 -gcp 11806.5 10989.381598793363 -170.00597321 63.93343374 -gcp 9498.0 11092.795625942685 -172.0055316 63.89820484 -gcp 4886.22541286895 11095.783559577678 -175.98828 63.89900914 -gcp 2579.3302782603996 10993.11105174092 -177.9858807 63.93581966 -gcp 520.5 8175.0 -179.99751972 64.94226414 -gcp 774.75 5544.9807692307695 -179.99953507 65.92205648 -gcp 1032.375 2914.962104072398 -179.99520757 66.93969664 -gcp 7191.0 5841.914404223227 -174.00241395 65.80732664 -gcp 7194.75 3204.792986425339 -174.00126305 66.82353149 -gcp 9336.0 5811.61387631976 -172.00515265 65.81800903 -gcp 7191.0 8485.060331825036 -174.00011436 64.82805468 -gcp 5054.25 5812.674962292609 -175.99181532 65.81914337 vlasenko/north/q-01_02.jpg q-01_02.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "1287.0 279.16407551360356\n13116.75 270.9328151027207\n14110.235237937282 10797.825478050963\n270.0 10812.952803997778" | \
  gdaltransform -tps -output_xy q-01_02.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  q-01_02.vrt q-01_02.cut.vrt

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