From 0346a8cb8fab91d11efc05de24f2916189e15e2d Mon Sep 17 00:00:00 2001
From: Martin Spott <martin@sphere.telascience.org>
Date: Mon, 13 Dec 2010 11:35:33 -0800
Subject: [PATCH] Intermediate step while trying to clean CORINE data ....

---
 gisscripts/grass01import.sh  |  3 ++-
 gisscripts/grass03patch.sh   |  3 ++-
 gisscripts/grass04clean.sh   | 10 ++++++----
 gisscripts/grass05extract.sh |  8 ++++----
 4 files changed, 14 insertions(+), 10 deletions(-)

diff --git a/gisscripts/grass01import.sh b/gisscripts/grass01import.sh
index c80a6262..2325f40b 100755
--- a/gisscripts/grass01import.sh
+++ b/gisscripts/grass01import.sh
@@ -44,7 +44,8 @@ case ${MODE} in
 	shp)
 	    # Local, pre-intersected Shapefiles
 	    DSN=${HOME}/live/corine
-	    PREFIX=nl_
+#	    PREFIX=nl_
+	    PREFIX=clc00_
 	    SELECTION=${DSN}/${PREFIX}c[0-9][0-9][0-9].shp
 	    SNAP=1
 	    SPAT=""
diff --git a/gisscripts/grass03patch.sh b/gisscripts/grass03patch.sh
index 0f0ba247..914802ee 100755
--- a/gisscripts/grass03patch.sh
+++ b/gisscripts/grass03patch.sh
@@ -23,7 +23,8 @@
 # $MODE-values for the "case" clause, like 'grass03patch.sh_shp'
 MODE=`basename ${0} | cut -f 2 -d \_`
 #
-PATCHMAP=clc00_nl
+#PATCHMAP=clc00_nl
+PATCHMAP=clc00
 
 case ${MODE} in
 	shp)
diff --git a/gisscripts/grass04clean.sh b/gisscripts/grass04clean.sh
index 2d378d75..625151b2 100755
--- a/gisscripts/grass04clean.sh
+++ b/gisscripts/grass04clean.sh
@@ -26,7 +26,8 @@ COMMAND=`basename ${0} | cut -f 1 -d \_`
 MODE=`basename ${0} | cut -f 2 -d \_`
 RUN=`basename ${0} | cut -f 3 -d \_`
 
-PATCHMAP=clc00_nl
+#PATCHMAP=clc00_nl
+PATCHMAP=clc00
 
 case ${MODE} in
 	shp)
@@ -42,13 +43,14 @@ MIN_AREA=10
 case ${RUN} in
 	first)
 	    v.clean input=${PATCHMAP} output=${PATCHMAP}_clean tool=snap,bpol,rmdupl,rmline,break,rmdupl,rmdangle,rmarea thresh=${SNAP},0,0,0,0,0,0,${MIN_AREA} type=boundary --verbose
+	    v.dissolve input=${PATCHMAP}_clean output=${PATCHMAP}_dissolved
 	;;
 	second)
-	    v.clean input=${PATCHMAP}_clean output=${PATCHMAP}_clobber tool=break,rmdupl,rmdangle type=boundary --verbose
+	    v.clean input=${PATCHMAP}_dissolved output=${PATCHMAP}_clobber tool=break,rmdupl,rmdangle type=boundary --verbose
 	;;
 	loop)
-	    g.remove vect=${PATCHMAP}_clean
-	    g.rename vect=${PATCHMAP}_clobber,${PATCHMAP}_clean
+	    g.remove vect=${PATCHMAP}_dissolved
+	    g.rename vect=${PATCHMAP}_clobber,${PATCHMAP}_dissolved
 	    ${COMMAND}_${MODE}_second
 	;;
 esac
diff --git a/gisscripts/grass05extract.sh b/gisscripts/grass05extract.sh
index 0a715706..e03f8370 100755
--- a/gisscripts/grass05extract.sh
+++ b/gisscripts/grass05extract.sh
@@ -31,8 +31,9 @@ cd ${RUNDIR}
 #
 
 MAPPINGFILE=${BASEDIR}/CORINEtoCStest.txt
-CLEANMAP=clc00_nl_clobber
-EXTRACTPREFIX=clc00_nl_
+#EXTRACTPREFIX=clc00_nl_
+EXTRACTPREFIX=clc00_
+CLEANMAP=${EXTRACTPREFIX}clobber
 
 case ${MODE} in
 	shp)
@@ -49,8 +50,7 @@ SELECTION=`v.category input=${CLEANMAP} type=centroid option=print | sort -n | u
 for CATEGORY in ${SELECTION}; do
     LAYER=`grep \^${CATEGORY} ${MAPPINGFILE} | awk '{print $2}' | sed -e "s/^cs_/${EXTRACTPREFIX}/g"`
     v.extract list=${CATEGORY} input=${CLEANMAP} output=${EXTRACTPREFIX}${CATEGORY} type=area
-    v.dissolve input=${EXTRACTPREFIX}${CATEGORY} output=${EXTRACTPREFIX}${CATEGORY}_dissolved
-    v.out.ogr input=${EXTRACTPREFIX}${CATEGORY}_dissolved dsn=${DUMPDIR}/${LAYER}.shp
+    v.out.ogr input=${EXTRACTPREFIX}${CATEGORY} dsn=${DUMPDIR}/${LAYER}.shp
 done
 
 # EOF