#!/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 1272.90192926045 336.8757634647418 -179.99287044 71.99991692 -gcp 11025.871061093247 298.66615769017216 -167.99636862 71.99947212 -gcp 12117.933762057877 10759.680732926152 -167.99700394 67.99968371 -gcp 309.8724259274456 10814.225861820261 -179.99705895 67.99990066 -gcp 6150.0 556.539353137146 -173.99638208 71.99984779 -gcp 6209.8125 11077.745280399777 -173.99701509 67.99979269 -gcp 2899.875 460.59244863964466 -177.99641314 71.99992178 -gcp 4525.875 535.6066074403109 -175.99639542 71.99988487 -gcp 7788.75 520.6744863964465 -171.99637316 71.9998106 -gcp 9404.25 436.67323709050527 -169.99636865 71.99977334 -gcp 11308.5 2913.9698778456413 -167.99655306 70.99972287 -gcp 11579.25 5532.969461410327 -167.99671874 69.99970975 -gcp 11847.0 8145.172820655191 -167.9968683 68.99969669 -gcp 10144.5 10929.997362576347 -169.99700401 67.99972008 -gcp 8183.25 11038.001943364798 -171.99700773 67.99975642 -gcp 4239.0 11052.220710716269 -175.99702609 67.99982885 -gcp 2274.75 10971.926152137701 -177.99704072 67.99986485 -gcp 551.25 8199.818711826763 -179.9969258 68.99991511 -gcp 792.0 5582.390338700722 -179.99677899 69.99992957 -gcp 1032.0 2959.5966129927824 -179.99661634 70.99994403 -gcp 6184.5 5821.77224458634 -173.99673096 69.99982017 vlasenko/north/r-01_02.jpg r-01_02.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "1272.90192926045 336.8757634647418\n11025.871061093247 298.66615769017216\n12117.933762057877 10759.680732926152\n309.8724259274456 10814.225861820261" | \
  gdaltransform -tps -output_xy r-01_02.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  r-01_02.vrt r-01_02.cut.vrt

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