Re: [GRASS-user] stream network extraction
Hi Markus, That is what was missing in stream-network analyzes! Great to have this! I didnt have a deeper look yet, but I will test and probably enjoy the script (in fast C-code). Thank you for that! Achim Markus Metz schrieb: Hi all, in grass-addons is a new module for clean stream network extraction from flow accumulation. It can extract streams from accumulation maps generated with any method (D8, FD8, D-Inf, DEMON, etc) and any software, granted that exactly the same elevation map is provided that was used to calculate flow accumulation. In case of r.terraflow, that would be the flooded (filled) output elevation map, not the elevation map used as input for r.terraflow. For e.g. TAS flow accumulation as input, it must be the elevation map used within TAS to calculate flow accumulation. r.stream.extract output is available as raster and vector with unique stream IDs. The vector is a complete network with stream segments and points for stream heads, confluences and outlets. The raster output can be used as input for the r.stream.* modules of Jarek, also in grass-addons. Regards, Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] stream network extraction
Markus Metz pisze: Hi all, in grass-addons is a new module for clean stream network extraction from flow accumulation. It can extract streams from accumulation maps generated with any method (D8, FD8, D-Inf, DEMON, etc) and any software, granted that exactly the same elevation map is provided that was used to calculate flow accumulation. In case of r.terraflow, that would be the flooded (filled) output elevation map, not the elevation map used as input for r.terraflow. For e.g. TAS flow accumulation as input, it must be the elevation map used within TAS to calculate flow accumulation. r.stream.extract output is available as raster and vector with unique stream IDs. The vector is a complete network with stream segments and points for stream heads, confluences and outlets. The raster output can be used as input for the r.stream.* modules of Jarek, also in grass-addons. Regards, Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user Hi Markus, Hi all. rest of r.stream* modules including new r.stream.stats for Hortonian statistics are now available by grass addons page and by grass add-ons svn repository. I actually move tutorial on r.stream* to grass wiki page. Jarek ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] reclass
Hi all, I have a problem reclassing a vector of polygons: I have areas with unique cats and a column value as integer. I want to reclass from column value and -here is the problem: keep the value information to identify the areas. In the manual it says: No table is created for reclassed layer if column option is used and column type is integer (there is nothing which could be written to the table). But couldn't be there a table with value-integer as cat? So how can I identify areas with the value? Thanks, Achim PS: workaround would be: v.to.rast with old values update centros of reclassed map from that raster map - but it seems a bit in-elegant.. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] reclass
PS: workaround would be: v.to.rast with old values update centros of reclassed map from that raster map - but it seems a bit in-elegant.. ...and which is a bit complicated as centroids cannot be used in v.what.rast... ...and conversion into points produces more than one point for some centroids (where one cat has more than two separate areas), so v.what.rast produces NULLs here... -so the workaround isn't one- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] stream network extraction
Jarosław Jasiewicz wrote: rest of r.stream* modules including new r.stream.stats for Hortonian statistics are now available by grass addons page and by grass add-ons svn repository. I actually move tutorial on r.stream* to grass wiki page. Great, IMHO the tutorials on your website are very helpful and show nicely all the various results the new r.stream.* modules can produce. Maybe a link to the wiki page can be added to the manuals? Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] stream network extraction
Markus Metz pisze: Jarosław Jasiewicz wrote: rest of r.stream* modules including new r.stream.stats for Hortonian statistics are now available by grass addons page and by grass add-ons svn repository. I actually move tutorial on r.stream* to grass wiki page. Great, IMHO the tutorials on your website are very helpful and show nicely all the various results the new r.stream.* modules can produce. Maybe a link to the wiki page can be added to the manuals? Markus M Will be, Jarek ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] reclass - multipolygon?
Well, I really stuck here: I havent found anything about playing with the common geom multipolygon in GRASS. I think that is what I need and what I thought I could handle with v.reclass. But it does not work. -Is there something in GRASS to create multipolygons? -Do they have one centroid or multiple with the same cat (as I thought) Any suggestions are very welcome achim achim schrieb: PS: workaround would be: v.to.rast with old values update centros of reclassed map from that raster map - but it seems a bit in-elegant.. ...and which is a bit complicated as centroids cannot be used in v.what.rast... ...and conversion into points produces more than one point for some centroids (where one cat has more than two separate areas), so v.what.rast produces NULLs here... -so the workaround isn't one- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Strange r.mapcalc results
Hi, I am trying to find out if 2 raster layers are exactly identical. I am not an experienced mapcalc user, so I am wondering, where the 3200204 NULL values come from, in Test1. The result in Test2 is what I expected. My obviously wrong understanding was that I wouldn't need the special operator ||| in Test1, as NULL values are already treated specially by isnull() Any hint from the experts? Hermann Test1 $ r.mapcalc 'result = isnull(mapA) isnull(mapT) || mapA == mapT ? 1 : 2' $ r.stats -c result 1 799796 * 3200204 Test2 $ r.mapcalc 'result = isnull(mapA) isnull(mapT) ||| mapA == mapT ? 1 : 2' $ r.stats -c result 1 400 Test3 $ r.mapcalc 'result = isnull(mapA) isnull(mapT) ? 1 : mapA == mapT ? 2 : 3' $ r.stats -c result 1 3200204 2 799796 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Problem finding r.mask in Windows script
Hi - I have a script that I am running on the Windows version of GRASS 6.4.0 installed using OSGEO4W. The script seems to run okay but the I am getting a message that I can't figure out. Early in the script I have the command: r.mask -r maskcats=* I also tried: r.mask -r MASK In both cases when the script runs I get the following message: C:/OSGeo4W/apps/grass/grass-6.4.0svn/scripts/r.shrub.cloud.resample: r.mask: command not found The r.mask command works fine when I run it from the command window or using the menu. Can anyone shed light on why the message is popping up? Ned ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] CYGWIN
Does anyone know when the new release of Cygwin (for Windows) will come out? Thanks. Lubica __ Attention: This e-mail and any files transmitted with it from Hatch Mott MacDonald are confidential and intended solely for use of the individual or entity to whom they are addressed. If you have received this e-mail in error please immediately notify the sender. _ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problem finding r.mask in Windows script
Hi Ned, On Wed, Sep 23, 2009 at 7:50 PM, Ned Horning horn...@amnh.org wrote: Hi - I have a script that I am running on the Windows version of GRASS 6.4.0 installed using OSGEO4W. The script seems to run okay but the I am getting a message that I can't figure out. Early in the script I have the command: r.mask -r maskcats=* I also tried: r.mask -r MASK In both cases when the script runs I get the following message: C:/OSGeo4W/apps/grass/grass-6.4.0svn/scripts/r.shrub.cloud.resample: r.mask: command not found The r.mask command works fine when I run it from the command window or using the menu. Can anyone shed light on why the message is popping up? just a curiosity: are you sure to run this script in the GRASS terminal? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Problem finding r.mask in Windows script
Markus, I ran the script in the lower pane of the Output - GIS.m window. I expect this is a stupid question but in the Windows version where is the GRASS terminal? Ned Markus Neteler wrote: just a curiosity: are you sure to run this script in the GRASS terminal? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Strange r.mapcalc results
Hermann wrote: I am trying to find out if 2 raster layers are exactly identical. I usually do g.region rast=map1,map2 r.univar map1 r.univar map2 and make sure the output is the same. another way is r.mapcalc diffmap = map1 - map2 r.univar diffmap see also wish #618 https://trac.osgeo.org/grass/ticket/618 Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user