#!/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 1245.7610116481417 486.1105563008101 11.99743929 63.99966641 -gcp 15100.929353802583 480.00246736336396 23.99750414 63.99987608 -gcp 16076.226442836196 10931.407984247082 23.99781135 59.99981764 -gcp 285.08416398713825 10932.44662687396 11.99775448 59.99961549 -gcp 8176.025300925005 807.126900479692 17.99745779 63.99977191 -gcp 8176.025300925005 11295.34475221948 17.9977707 59.99971721 -gcp 2394.0 589.606334841629 12.9992812 63.95726087 -gcp 3549.0 666.8665158371041 14.00272034 63.92559388 -gcp 4707.0 727.7707390648568 15.00594435 63.90061393 -gcp 5863.5 771.2926093514328 16.00569935 63.8827234 -gcp 7018.5 801.2843137254902 17.00273532 63.87034404 -gcp 9333.0 801.265460030166 18.9988457 63.86998252 -gcp 10488.0 772.7511312217193 19.99559252 63.88139043 -gcp 11647.5 727.8144796380091 20.99755023 63.89948066 -gcp 12802.5 660.2805429864254 21.99794508 63.92677857 -gcp 13952.93852770872 580.7707268504879 22.99691139 63.95897703 -gcp 15348.75 3093.8423831070886 24.0008542 62.95028628 -gcp 15588.75 5703.103318250377 23.99800341 61.93574287 -gcp 15832.5 8318.409502262444 23.99782923 60.95144845 -gcp 14763.75 11040.149321266968 22.99452065 59.9609035 -gcp 13449.0 11127.943438914028 21.99295461 59.92952394 -gcp 12133.5 11206.67043740573 20.99310399 59.90141322 -gcp 10815.0 11253.906862745098 19.99347713 59.88456235 -gcp 9499.5 11286.0 18.99790499 59.87385705 -gcp 6859.5 11285.40460030166 16.99940829 59.87334982 -gcp 5546.25 11255.41101055807 16.00473504 59.88406038 -gcp 4225.087023356181 11205.152876965712 15.00257527 59.90201229 -gcp 2908.5 11135.429487179488 14.00170928 59.9269346 -gcp 1594.5 11044.5 12.99805847 59.95994857 -gcp 521.25 8322.851055806937 11.99468652 60.95050362 -gcp 761.25 5710.253016591252 11.99430006 61.93423794 -gcp 1002.75 3096.820135746606 11.99507916 62.95086899 -gcp 8176.5 6057.281297134238 17.99863137 61.80113041 -gcp 11890.5 5973.302413273002 20.99495272 61.83284644 -gcp 4467.0 5967.281297134238 15.00495206 61.83573025 -gcp 5784.0 3396.3665158371036 16.00537918 62.83210067 -gcp 10570.5 3396.1855203619903 19.9955528 62.83160731 -gcp 10735.5 8634.997360482654 19.99567884 60.83468686 -gcp 5624.625 8637.202300150828 16.00428741 60.83408502 -gcp 3390.0 3282.4758672699845 14.00330347 62.87726985 -gcp 3069.0 8516.294871794871 14.00213621 60.87882845 -gcp 13288.5 8517.83031674208 21.99463196 60.87782178 -gcp 12966.0 3282.5460030165905 21.9978575 62.87609563 vlasenko/north/p-33_34.jpg p-33_34.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "1245.7610116481417 486.1105563008101\n15100.929353802583 480.00246736336396\n16076.226442836196 10931.407984247082\n285.08416398713825 10932.44662687396" | \
  gdaltransform -tps -output_xy p-33_34.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  p-33_34.vrt p-33_34.cut.vrt

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