#!/bin/bash

baseDir="."
mindepth=1
maxdepth=100000
help=0

formats='.gpkg$|.shp$|.fgb$'

while [ "$#" -gt 0 ]; do
  case "$1" in
    --dir=*) baseDir="$(echo -n ${1#*=} |sed 's/[^\/\-\_0-9a-z\.]//gi')"; shift 1;;
    --formats=*) formats="$(echo -n ${1#*=} |sed 's/[^0-9a-z\,]//gi' | sed 's/\,/\$\|\./g')"; shift 1;;
    --mindepth=*) mindepth="$(echo -n ${1#*=} |sed 's/[^0-9]//g')"; shift 1;;
    --maxdepth=*) maxdepth="$(echo -n ${1#*=} |sed 's/[^0-9]//g')"; shift 1;;
    --help) help=1; shift 1;;
    *) shift 1;;
  esac
done
err="Usage:"
[[ ${#formats} -gt 0 ]] && formats=".${formats}\$"
[[ ${#formats} -le 1 ]] && err="'${formats}' --formats is empty" && help=1

if [ $help -eq 1 ] ; then
echo "$err"
cat<<EOF
$0
   --dir       : directory to search, default current directory
   --formats   : extensions to search for, default .gpkg, .shp and .fgb 
                 any vector format extension GDAL understands is valid.
   --mindepth  : min directory depth, default current directory
   --maxdepth  : max directory depth, default no limit 
   --help      : this output
Example: $0 --dir="/export/data" --formats="gpkg" 

EOF
   exit
fi


echo "path,layer,epsg,wkt_geometry"
for path in $(echo -e "$(find $baseDir -mindepth $mindepth -maxdepth $maxdepth |grep -Ei "$formats")" |tr "\n" " ") ; do
   cols="$(gdal vector info $path --of=json 2> /dev/null| jq '.layers[]|[.name,.geometryFields[].coordinateSystem.projjson.id.code,.geometryFields[].extent]|tostring' |sed 's/[^a-z0-9\._,-]//gi'|tr -d '\\'|grep ','|grep -v ',null' |tr "\n" " ")"
   for layer in $cols ; do
      IFS=',' read p lay epsg xmin ymin xmax ymax <<< "$path,$layer"
      printf -v csv '"%s","%s",%s,"POLYGON((%.6f %.6f,%.6f %.6f,%.6f %.6f,%.6f %.6f,%.6f %.6f))"' "$p" "$lay" "$epsg" "$xmin" "$ymin" "$xmin" "$ymax" "$xmax" "$ymax" "$xmax" "$ymin" "$xmin" "$ymin" || continue
      echo "$csv"
   done
done