#!/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 4313.25 1242.8153803442533 83.99685377 84 -gcp 14145.75 1318.7294558578567 108 84.00102482 -gcp 17268.75 16896.88450860633 108.000783 80.00095927 -gcp 922.5 16784.45197112715 84 80 -gcp 5130.0 1410.2581898945032 85.99721391 84.0009948 -gcp 5940.0 1546.748195446974 87.99757745 84.00100425 -gcp 6757.5 1653.2765130483065 89.99794394 84.00101238 -gcp 7578.0 1738.7973348139924 91.99831294 84.00101917 -gcp 8397.375 1792.578498056635 93.99868399 84.00102462 -gcp 9225.75 1811.3961687951137 95.99905665 84.00102872 -gcp 10047.75 1802.4053303720157 97.99943046 84.00103147 -gcp 10874.25 1763.4069961132702 99.99980496 84.00103286 -gcp 11693.25 1689.8895058300943 102.0001797 84.00103288 -gcp 12517.5 1592.3928373126043 104.00055422 84.00103155 -gcp 13327.5 1469.3878400888395 106.00092807 84.00102886 -gcp 14923.125 5206.629511382565 108.00111568 83.00100889 -gcp 15702.0 9099.477859522487 108.00097696 82.00099265 -gcp 16488.75 12995.47654081066 108.00086916 81.00097611 -gcp 15917.25 17142.51346474181 106.00055865 80.00096327 -gcp 14562.0 17346.23903387007 104.00033361 80.00096594 -gcp 13194.75 17496.996946141033 102.00010817 80.00096726 -gcp 11832.0 17610.234453081623 99.9998826 80.00096723 -gcp 10461.000048113201 17671.00777346173 97.99965716 80.00096585 -gcp 9088.5 17680.596404775126 95.99943215 80.00096313 -gcp 7717.125 17649.10577456968 93.99920783 80.00095907 -gcp 6350.625 17565.470363686843 91.99898447 80.00095368 -gcp 4985.625 17442.45599666852 89.99876235 80.00094695 -gcp 3624.75 17265.17573570239 87.99854174 80.0009389 -gcp 2268.75 17042.717448639647 85.99832291 80.00092954 -gcp 1771.125 12894.887423653527 83.99789774 81.0009356 -gcp 2621.25 9007.527588839534 83.99763701 82.00095204 -gcp 3465.375 5121.049764019988 83.99730148 83.00096819 -gcp 9124.5 13707.948639644643 95.99936967 81.00097999 -gcp 9185.25 5777.139228206553 95.99919089 83.00101279 -gcp 9153.375 9739.082940033315 95.99929149 82.00099654 -gcp 12065.25 5643.684827873404 102.00015413 83.00101694 -gcp 12446.25 9593.782516657413 102.00013496 82.00100068 -gcp 12815.4375 13538.241567184898 102.00012007 81.00098412 -gcp 5428.875 13490.24541921155 89.99862618 81.00096376 -gcp 5872.875 9542.883189894503 89.99845579 82.00098027 -gcp 6317.0625 5595.605809272626 89.99823652 83.00099648 vlasenko/north/u-45_46_47_48.jpg u-45_46_47_48.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "4313.25 1242.8153803442533\n14145.75 1318.7294558578567\n17268.75 16896.88450860633\n922.5 16784.45197112715" | \
  gdaltransform -tps -output_xy u-45_46_47_48.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  u-45_46_47_48.vrt u-45_46_47_48.cut.vrt

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