footprint.sh

#!/bin/bash
#
# ***************************************************************************
# Name                 : Footprint
# Description          : Create footprint of image with GeoJSON format
#
# Arguments: 
# $1: Image
#
# Dependencies         : gdal 1.10.1(gdal_calc.py, gdal_sieve.py, gdal_edit.py and gdal_polygonize.py, ogr2ogr)
#
# ***************************************************************************
# begin                : 2015-03-02 (yyyy-mm-dd)
# copyright            : (C) 2015 by Luiz Motta
# email                : motta dot luiz at gmail.com
# ***************************************************************************
#
# Revisions
#
# 2016-03-31:
# - Create multipolygon, added a field of total of parts 
#
# ***************************************************************************
#
# Example:
#   footprint.sh LC8_229-066_20140724_LGN00_r6g5b4.tif
#
# ***************************************************************************
# *                                                                         *
# *   This program is free software; you can redistribute it and/or modify  *
# *   it under the terms of the GNU General Public License as published by  *
# *   the Free Software Foundation; either version 2 of the License, or     *
# *   (at your option) any later version.                                   *
# *                                                                         *
# ***************************************************************************
#
msg_error(){
  local name_script=$(basename $0)
  echo "Usage: $name_script <image> <outsize>" >&2
  echo "<image> is the image for calculate footprint" >&2
  exit 1
}
#
totalargs=1
#
if [ $# -ne $totalargs ] ; then
  msg_error
  exit 1
fi
#
in_img=$1
#
if [ ! -f "$in_img" ]; then
  echo "The file '$in_img' not exist" >&2
  exit 1
fi
#
dir_img=$(dirname $in_img)
basename_img=$(basename $in_img)
name_img=${basename_img%.*}
#
footprint_geojson=$dir_img"/"$name_img".geojson"
zero_one_img=$dir_img"/"$name_img"_0x1.tif"
sieve_img=$dir_img"/"$name_img"_0x1_sieve.tif"
#
if [ -f $footprint_geojson ]; then
 rm $footprint_geojson
fi
# Processing
printf "'$name_img'...1"
gdal_calc.py -A $in_img --A_band 1 --type Byte --calc "A>0" --outfile $zero_one_img >/dev/null
#
code=$?
if [ "$code" != 0 ]
then
  exit $code
fi
#
printf ".2"
gdal_sieve.py -q -st 100 -8 $zero_one_img -nomask $sieve_img
code=$?
if [ "$code" != 0 ]
then
  exit $code
fi
#
gdal_edit.py -a_nodata 0 $sieve_img
code=$?
if [ "$code" != 0 ]
then
  exit $code
fi
#
printf ".3"
gdal_polygonize.py $sieve_img -q -b 1 -f "GeoJSON" $footprint_geojson
code=$?
if [ "$code" != 0 ]
then
  exit $code
fi
# 
tmp_file=$footprint_geojson".tmp"
ogr2ogr -dialect SQLITE -sql "SELECT DN, ST_Union( geometry ), COUNT( 1 ) as total_part FROM OGRGeoJSON WHERE DN = 1" -f "GeoJSON" $tmp_file $footprint_geojson
mv -f $tmp_file $footprint_geojson
#
ssed="s|\"DN\": 1,|\"path\": \"$dir_img\", \"image\": \"$basename_img\",|"
sed -i "$ssed" "$footprint_geojson"
printf "\rCreated: %-10s\n" $footprint_geojson
# Cleanup
rm $zero_one_img $sieve_img
#
exit 0