merge_dom_cells.sh
5fbf2f59e661master

merge_dom_cells.sh
View Options

			#!/bin/bash

			## == Print inputs (for logging purposes only)
			echo
			echo Entered grass script
			echo
			echo Path to DOM cells: $1
			echo File with DOM cells to merge: $2
			echo New layer name: $3

			## == read DOM cells from the file
			MYPATH=$1

			while read DOM_CELL; do

			echo $DOM_CELL
			name="${DOM_CELL%.*}"

			r.in.gdal -o input=$MYPATH\/$name\.asc output=TMP_$name

			done < $2

			## == get a list of all DOM cells and set the region spanning all
			MAPS=`g.list type=raster sep=',' pattern=TMP_*`
			g.region -p raster=$MAPS

			# patch all cells together to one
			r.patch input=$MAPS output=$3 --overwrite

			# remove the individual parts
			g.remove -f type=raster pattern=TMP_*

			echo
			echo Finished executing grass script
			echo