#!/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 2265.0 3507.0 20.876666666666665 56.75 -gcp 21636.0 3546.0 21.326666666666668 56.75 -gcp 21741.0 16941.0 21.326666666666668 56.583333333333336 -gcp 2271.0 16986.0 20.876666666666665 56.583333333333336 lvva/Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas_Liepajas_apkartne___6828__A1_L22_21_02.JPG Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas_Liepajas_apkartne___6828__A1_L22_21_02.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e "2265.0 3507.0\n21636.0 3546.0\n21657.0 3506.99666851749\n23085.0 3480.0016657412552\n23052.0 10290.043309272627\n23043.0 23612.983342587453\n10512.0 23617.838978345368\n8454.0 23622.033314825097\n8359.5 23569.482509716825\n8169.0 23578.53914491949\n8004.0 23646.119933370348\n7710.0 23640.076624097725\n7122.0 23447.973348139923\n6300.0 23454.033314825097\n6042.0 23345.97001665741\n5748.0 23351.880066629652\n5112.0 23363.966685174903\n4656.0 23484.02665186008\n3444.0 23507.953359244864\n2256.0 23568.079955580233\n2271.0 16986.0" | \
  gdaltransform -tps -output_xy Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas_Liepajas_apkartne___6828__A1_L22_21_02.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}')))" \
  Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas_Liepajas_apkartne___6828__A1_L22_21_02.vrt Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas_Liepajas_apkartne___6828__A1_L22_21_02.cut.vrt

# Generate tiles
gdal2tiles.py -r $GDAL_SAMPLING_TILE --xyz -z $GDAL_TILE_ZOOMS -x --processes=$GDAL_TILE_PROCESSES Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas_Liepajas_apkartne___6828__A1_L22_21_02.cut.vrt Kurzemes,Vidzemes,Kaunas_un_Vitebskas_gubernas_Liepajas_apkartne___6828__A1_L22_21_02.xyz