#!/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 1817.625 252.81308995002775 59.99575681 80.00091887 -gcp 10041.375 235.20891171571347 83.99864051 80.00091888 -gcp 11690.625 7973.336410327596 83.99864051 76.00084917 -gcp 230.25 7997.362437534703 59.99575681 76.00084917 -gcp 5934.75 653.4136500754148 71.99686289 80.00082843 -gcp 3866.25 558.1798642533936 65.99628952 80.00076705 -gcp 7999.5 548.4400452488687 77.99747065 80.00087923 -gcp 8841.0 8404.02149321267 77.99818435 76.00081009 -gcp 5962.125 8559.842383107089 71.99774808 76.00076003 -gcp 3078.0 8425.218514328808 65.99733648 76.00069954 -gcp 2499.0 376.6057745696835 61.99592949 80.00072057 -gcp 3179.25 478.66240977234867 63.99610713 80.00074435 -gcp 4555.5 615.9232370905052 67.99647643 80.00078865 -gcp 5246.25 645.9398945030539 69.99666763 80.00080912 -gcp 6624.747561599345 643.6307606885064 73.99706198 80.00084657 -gcp 7314.747561599345 609.1116046640755 75.99726465 80.00086351 -gcp 8682.75 464.414908384231 79.99767973 80.00089371 -gcp 9363.75 360.8574403109384 81.99789164 80.00090693 -gcp 10450.5 2169.9215713492504 83.99827643 79.00090187 -gcp 10863.242140822353 4103.261576084013 83.99841819 78.00088458 -gcp 11273.625 6038.318087173792 83.99853798 77.00086701 -gcp 10746.0 8152.672820655191 81.99848655 76.00083739 -gcp 9797.25 8296.716546363132 79.99833443 76.00082436 -gcp 7883.25 8494.6457523598 75.99803647 76.0007946 -gcp 6921.75 8540.39117157135 73.99789099 76.00077791 -gcp 5002.163610876545 8544.863151543126 69.99760791 76.000741 -gcp 4041.413610876545 8505.0910804715 67.99747065 76.00072083 -gcp 2125.0828052842044 8312.344458551575 63.99720556 76.00067717 -gcp 1176.0 8167.596612992781 61.99707804 76.00065374 -gcp 631.875 6060.032898389783 59.99672436 77.00064621 -gcp 1025.25 4123.927887284842 59.99645597 78.00066293 -gcp 1420.125 2189.2554136590784 59.99613838 79.00067945 -gcp 5949.75 4607.5127706829535 71.99737981 78.00079473 -gcp 8420.25 4478.747015546918 77.99788742 78.00084519 -gcp 3475.125 4488.104108828428 65.99690091 78.00073376 -gcp 5954.25 6584.356399222654 71.99757824 77.00077751 -gcp 5941.5 2630.327248750694 71.997145 79.00081171 vlasenko/north/t-41_42_43_44.jpg t-41_42_43_44.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "1817.625 252.81308995002775\n10041.375 235.20891171571347\n11690.625 7973.336410327596\n230.25 7997.362437534703" | \
  gdaltransform -tps -output_xy t-41_42_43_44.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  t-41_42_43_44.vrt t-41_42_43_44.cut.vrt

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