#! /bin/bash -
set -e
dir="$(dirname "${1}")"
cd "$dir"
pwd


ogr() {
#  echo docker run --rm -v /home:/home osgeo/gdal:alpine-normal-latest ogr2ogr $*
#  docker run --rm -u 1000:1004 -v /home:/home osgeo/gdal:alpine-normal-latest echo ogr2ogr $*
#  docker run --rm -u 1000:1004 -v /home:/home osgeo/gdal:ubuntu-full-latest ogr2ogr $*
#  echo docker run --rm -u 1000:1004 -v /home:/home bjreppen/gdal-mbtiles:latest ogr2ogr $*
  echo + $(date +"%T") $*
  docker run --rm -u 1000:1004 --env OGR_GEOJSON_MAX_OBJ_SIZE=2000 -v /home:/home bjreppen/gdal-mbtiles:latest ogr2ogr $*
  echo - $(date +"%T") $*
#  docker run --rm -u 1000:1004 -v /home:/home klokantech/gdal ogr2ogr $*
}

gdaladdo() {
#  echo docker run --rm -u 1000:1004 -v /home:/home bjreppen/gdal-mbtiles:latest gdaladdo $*
  echo + $(date +"%T") $*
  docker run --rm -u 1000:1004 -v /home:/home bjreppen/gdal-mbtiles:latest gdaladdo $*
  echo - $(date +"%T") $*
}

gdal_translate() {
#  echo docker run --rm -u 1000:1004 --env OGR_GEOJSON_MAX_OBJ_SIZE=2000 -v /home:/home bjreppen/gdal-mbtiles:latest gdal_translate $*
  echo + $(date +"%T") $*
  docker run --rm -u 1000:1004 -v /home:/home bjreppen/gdal-mbtiles:latest gdal_translate $*
  echo - $(date +"%T") $*
}


#if [[ -e polygon_no_dissolve.4326.geojson ]]; then
if [[ ! -e polygon.4326.geojson ]]; then
  echo Dissolving..
  if [[ -e polygon.4326.geojson ]]; then
  	rm polygon.4326.geojson
  fi
  ogr2ogr $PWD/polygon.4326.geojson  $PWD/polygon_no_dissolve.4326.geojson -dialect sqlite -sql "SELECT ST_Union(geometry), kode FROM polygons GROUP BY kode"
  rm polygon_no_dissolve.4326.geojson
fi

if [[ ! -e polygon.4326.geojson ]]; then
  ogr -overwrite -f GeoJSON $PWD/polygon.4326.geojson -t_srs EPSG:4326 $PWD/polygon.32633.geojson
fi

if [[ ! -e polygon.32633.geojson ]]; then
  ogr -overwrite -f GeoJSON $PWD/polygon.32633.geojson -t_srs EPSG:32633 $PWD/polygon.4326.geojson
fi

#if [[ ! -e polygon.3857.mbtiles ]]; then
  echo build polygon.3857.mbtiles
#  gdal_translate $PWD/polygon.4326.geojson $PWD/polygon.3857.mbtiles 
#rm polygon.3857.mbtiles
  ogr -of MVT $PWD/polygon.3857.mbtiles $PWD/polygon.4326.geojson -dsco MAXZOOM=12 -dsco BUFFER=0
#  ogr $PWD/polygon.3857.mbtiles $PWD/polygon.4326.geojson
#  gdaladdo $PWD/polygon.3857.mbtiles
#fi

if [[ ! -e grid.32633.tif ]]; then
  echo make grid.32633.tif
  gdal_rasterize -ot Byte -burn 255 -tr 1000 1000 -te -119000 6420000 1158000 8000000 $PWD/polygon.32633.geojson $PWD/grid.32633.tif
fi

if [[ ! -e grid.32633.png ]]; then
#  convert grid.32633.tif grid.32633.png
  echo make grid.32633.png
  gdal_translate -of PNG $PWD/grid.32633.tif $PWD/grid.32633.png 
  echo make grid.liten.png
  convert $PWD/grid.32633.png -resize 50% $PWD/grid_liten.png
fi