#!/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 1800.3638466622606 2263.7823529411767 25.55 56.22222222222222 -gcp 14962.846992729676 2349.3991596638652 26 56.22222222222222 -gcp 14941.927296761402 14297.39075630252 26 56 -gcp 1711.6368142762724 14209.875630252101 25.55 56 Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Subate_6828_A1_L25_65_02.JPG Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Subate_6828_A1_L25_65_02.vrt

# Add cutline to VRT
gdalwarp -r $GDAL_SAMPLING_WARP -tps -dstalpha -cutline_srs EPSG:4326 \
  -cutline "POLYGON(($(echo -e \"1259.8730998017181 11185.653781512605\n1355.881031064111 6707.993277310924\n1439.7620621282222 1391.5092436974783\n2177.958360872439 737.9899159663867\n15846.582947785855 875.1731092436974\n16176.055518836747 1565.966386554622\n16176.0 3404.9848739495797\n16338.045604758758 3614.764705882353\n16362.111037673496 4664.988235294118\n16212.023793787179 5964.1949579831935\n16212.031725049572 14288.584873949581\n14941.927296761402 14297.39075630252\n1094.5915399867813 14207.99831932773\" | \
  gdaltransform -tps -output_xy Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Subate_6828_A1_L25_65_02.vrt | \
  awk 'NR==1{first=$0} {printf "%s %s,", $1,$2} END{print " " first}'))))" \
  Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Subate_6828_A1_L25_65_02.vrt Latvijas_dalas_Lietuvas_dalas_un_Igaunijas_top.kartes_Subate_6828_A1_L25_65_02.cut.vrt

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