Fixed zone map region bug, introduced with the -r

flag for map specific regions.
Modified the output of the -g flag in case of zones. Now the shell
style is consistent. Instead of "\nzone num cat\n\n" -> "zone=num;cat"
is printed to stdout. 


git-svn-id: https://svn.osgeo.org/grass/grass/trunk@66602 15284696-431f-4ddb-bdfa-cd5b030d7da7

raster/r.univar/r.univar_main.c

@@ -178,11 +178,11 @@ int main(int argc, char *argv[])
 	if(param.use_rast_region->answer) {
     	    mapset = G_find_raster2(*p, "");
 	    Rast_get_cellhd(*p, mapset, &region);
+	    /* Set the computational region */
+	    Rast_set_window(&region);
         } else {
     	    G_get_window(&region);
         }
-	/* Set the computational region */
-	Rast_set_window(&region);
 
 	fd = open_raster(*p);
 

raster/r.univar/stats.c

@@ -134,19 +134,22 @@ int print_stats(univar_stat * stats)
 	sprintf(sum_str, "%.15g", stats[z].sum);
 	G_trim_decimal(sum_str);
 
-	if (zone_info.n_zones) {
-	    int z_cat = z + zone_info.min;
-	    
-	    fprintf(stdout, "\nzone %d %s\n\n", z_cat, Rast_get_c_cat(&z_cat, &(zone_info.cats)));
-	}
 
 	if (!param.shell_style->answer) {
+	    if (zone_info.n_zones) {
+	    	int z_cat = z + zone_info.min;
+	    	fprintf(stdout, "\nzone %d %s\n\n", z_cat, Rast_get_c_cat(&z_cat, &(zone_info.cats)));
+	    }
 	    fprintf(stdout, "total null and non-null cells: %lu\n", stats[z].size);
 	    fprintf(stdout, "total null cells: %lu\n\n", stats[z].size - stats[z].n);
 	    fprintf(stdout, "Of the non-null cells:\n----------------------\n");
 	}
 
 	if (param.shell_style->answer) {
+	    if (zone_info.n_zones) {
+	    	int z_cat = z + zone_info.min;
+	    	fprintf(stdout, "zone=%d;%s\n", z_cat, Rast_get_c_cat(&z_cat, &(zone_info.cats)));
+	    }
 	    fprintf(stdout, "n=%lu\n", stats[z].n);
 	    fprintf(stdout, "null_cells=%lu\n", stats[z].size - stats[z].n);
 	    fprintf(stdout, "cells=%lu\n", stats->size);

raster/r.univar/testsuite/test_r_univar.py

@@ -18,8 +18,10 @@ class TestRasterUnivar(TestCase):
         cls.del_temp_region()
 
     def tearDown(self):
+    	pass
         self.runModule("g.remove", type="raster", name="map_a")
         self.runModule("g.remove", type="raster", name="map_b")
+        self.runModule("g.remove", type="raster", name="zone_map")
 
     def setUp(self):
         """Create input data
@@ -29,6 +31,8 @@ class TestRasterUnivar(TestCase):
                        overwrite=True)
         self.runModule("r.mapcalc", expression="map_b = 200 + row() + col()",
                        overwrite=True)
+        self.runModule("r.mapcalc", expression="zone_map = if(row() < 20, 1,2)",
+                       overwrite=True)
 
     def test_1(self):
         # Output of r.univar
@@ -135,6 +139,40 @@ class TestRasterUnivar(TestCase):
         self.assertModuleKeyValue(module="r.univar", map=["map_a","map_b"], flags="rg",
                                   reference=univar_string, precision=3, sep='=')
 
+
+    def test_1_zone(self):
+        """
+        multiple maps and zone
+        :return:
+        """
+
+        # Output of r.univar
+        univar_string="""zone=1;
+			n=1710
+			null_cells=0
+			cells=1710
+			min=102
+			max=209
+			range=107
+			mean=155.5
+			mean_of_abs=155.5
+			sum=265905
+			zone=2;
+			n=6390
+			null_cells=0
+			cells=1710
+			min=121
+			max=280
+			range=159
+			mean=200.5
+			mean_of_abs=200.5
+			sum=1281195"""
+
+        self.runModule("g.region", res=1)
+        self.assertModuleKeyValue(module="r.univar", map=["map_a"], 
+                                  zones="zone_map",flags="g",
+                                  reference=univar_string, precision=3, sep='=')
+
 class TestAccumulateFails(TestCase):
 
     def test_error_handling(self):