extract_coordinate_dims.sh

#!/bin/bash

FILE_IN=$1

file_name=$( basename $FILE_IN | sed 's/.tif//' )

function gdal_extent_gdalwarp_te() {
    if [ -z "$1" ]; then
        echo "Missing arguments. Syntax:"
        echo "  gdal_extent <input_raster>"
        return
    fi
    EXTENT=$(gdalinfo "$1" |\
        grep "Lower Left\|Upper Right" |\
        sed "s/Lower Left  //g;s/Upper Right //g;s/).*//g" |\
        #tr "\n" " " |\
        sed 's/ *$/ /g' |\
        tr -d "\n[(]'"|\
        sed 's/,/ /g')
    echo -n "$EXTENT"
}

#echo $FILE_IN

ExtractedExtent=$(gdal_extent_gdalwarp_te "${FILE_IN}")
#echo $ExtractedExtent

set -- $ExtractedExtent
#echo "flip extent (ulx uly lrx lry) = "$1 $4 $3 $2

ulx=$1
uly=$4
lrx=$3
lry=$2

echo "$ulx"
echo "$uly"
echo "$lrx"
echo "$lry"