#!/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 1225.5793408360128 468.06996113270407 23.99750414 63.99987608 -gcp 15105.198231511255 483.9444752915047 35.99767807 64.00007586 -gcp 16076.0627236906 10980.369783719289 35.99796387 60.00001024 -gcp 243.1183426163255 10982.437811074731 23.99781135 59.99981764 -gcp 8172.375 797.6718489727928 29.99757783 63.99997777 -gcp 8158.875 11363.000971682399 29.99787597 59.99991567 -gcp 2385.0 570.0 25.00655152 63.95832808 -gcp 3540.0 657.9404223227752 26.00578311 63.92398134 -gcp 4696.5 729.933257918552 27.00504825 63.89536845 -gcp 5858.25 766.6885369532429 28.00627869 63.88103122 -gcp 7007.25 794.4404223227754 28.99574993 63.87033284 -gcp 9326.03493267508 799.2003473767015 30.99196819 63.86947659 -gcp 10492.28493267508 775.1882810116941 31.99705815 63.87973536 -gcp 11644.28493267508 723.4122629121466 32.99192655 63.90124172 -gcp 12800.25 666.8989441930619 33.99216103 63.92469112 -gcp 13956.0 581.3559577677224 34.99597091 63.95994869 -gcp 15344.25 3108.308446455505 35.99467079 62.95020925 -gcp 15589.5 5733.0 35.98382371 61.93554409 -gcp 15829.5 8350.587481146304 35.99529636 60.95315201 -gcp 14763.0 11085.120663650074 34.99489272 59.96275071 -gcp 13447.5 11181.168929110105 33.993161 59.92864548 -gcp 12124.5 11262.209653092006 32.98884914 59.89991992 -gcp 10812.0 11316.37330316742 31.9953089 59.88076362 -gcp 9487.5 11340.38310708899 30.99462475 59.87231502 -gcp 6845.900871986471 11347.883861236802 29.00098904 59.86978801 -gcp 5522.900871986471 11317.868778280543 28.00231747 59.88049071 -gcp 4204.242721313416 11265.360680182965 27.00574294 59.89916933 -gcp 2887.242721313416 11188.822218644504 26.00838067 59.92638183 -gcp 1565.742721313416 11095.7754614801 25.00517402 59.95947679 -gcp 487.48050037434814 8352.218855374791 23.99664503 60.95129238 -gcp 739.5 5727.116515837104 24.00196124 61.93307588 -gcp 980.25 3106.56334841629 23.99859862 62.94603972 -gcp 8163.75 6082.039969834088 29.99345108 61.79983249 -gcp 11877.75 5987.5475113122175 32.98512831 61.83637326 -gcp 12959.25 3289.677978883861 33.99038685 62.87811449 -gcp 13282.5 8556.526018099546 33.99025662 60.87717225 -gcp 10728.0 8677.995098039217 31.99260434 60.83240141 -gcp 10566.0 3402.330316742082 31.99161069 62.83314424 -gcp 4452.75 5993.6674208144805 27.00718178 61.83261682 -gcp 3376.5 3289.6161387631973 26.00618603 62.8749898 -gcp 5776.5 3402.1726998491704 28.00686269 62.83166 -gcp 5605.5 8673.0 28.00570281 60.83204305 -gcp 3049.5 8555.380467571644 26.00703207 60.87684868 vlasenko/north/p-35_36.jpg p-35_36.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "1225.5793408360128 468.06996113270407\n15105.198231511255 483.9444752915047\n16076.0627236906 10980.369783719289\n243.1183426163255 10982.437811074731" | \
  gdaltransform -tps -output_xy p-35_36.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  p-35_36.vrt p-35_36.cut.vrt

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