merge_dom_cells.sh
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/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