Thanks for quick response and sorry for delayed response from my side.<br>This is the second time , I have noticed that v.to.rast changes the region settings equal to display resolution in GRASS70.<br>As, it cannot be reproduced at test side.<br>
I will send more data the next time. <br>At present, I have created another mapset for working on v.to.rast and my base rasters and vectors are kept securely in "PERMANENT" mapset. I haven't faced the problem since.<br>
Thanks again. <br><br><div class="gmail_quote">On Tue, Apr 10, 2012 at 9:30 PM, <span dir="ltr"><<a href="mailto:grass-user-request@lists.osgeo.org">grass-user-request@lists.osgeo.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Send grass-user mailing list submissions to<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:grass-user-request@lists.osgeo.org">grass-user-request@lists.osgeo.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:grass-user-owner@lists.osgeo.org">grass-user-owner@lists.osgeo.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of grass-user digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Re: r.stream.extract Syntax Problem (Markus Metz)<br>
2. Re: v.to.rast in grass70 (Markus Metz)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 10 Apr 2012 09:14:55 +0200<br>
From: Markus Metz <<a href="mailto:markus.metz.giswork@googlemail.com">markus.metz.giswork@googlemail.com</a>><br>
Subject: Re: [GRASS-user] r.stream.extract Syntax Problem<br>
To: Rich Shepard <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>><br>
Cc: <a href="mailto:grass-users@lists.osgeo.org">grass-users@lists.osgeo.org</a><br>
Message-ID:<br>
<CAG+h=FFXNBRd8E_K88NeQPDGZXLDa71sxyDLJzez+a=<a href="mailto:tyGjR8g@mail.gmail.com">tyGjR8g@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
[Answering to your latest post in a new thread, but I prefer to keep<br>
it all in one thread.]<br>
<br>
<br>
On Tue, Apr 10, 2012 at 1:54 AM, Rich Shepard <<a href="mailto:rshepard@appl-ecosys.com">rshepard@appl-ecosys.com</a>> wrote:<br>
> On Mon, 9 Apr 2012, Markus Metz wrote:<br>
><br>
>> The computational region should be aligned to the DEM in order to get<br>
>> accurate results. Aligning to the resolution of the DEM is not enough,<br>
>> i.e. you should use g.region align=DEM instead of g.region res=XXX -a.<br>
>> r.stream.extract should work if only a subregion is to be analysed by<br>
>> setting the computational region accordingly and/or creating a MASK e.g.<br>
>> for a specific basin.<br>
><br>
><br>
> Markus, madi, et al.:<br>
><br>
> The region does need to be set with a mask (in my case, and based on the<br>
> v.to.rast output of the individual subbasins).<br>
<br>
v.to.rast? Basins as created by r.watershed, r.stream.basins and<br>
r.water.outlet are already raster maps.<br>
><br>
> The last advice in the example on the manual page is to compare the<br>
> weighted and unweighted extracted stream maps and "see which looks better."<br>
> As a newcomer to this module I have no experience on what "better" looks<br>
> like.<br>
><br>
> Attached are two .pdf files. One shows the weighted output streams, the<br>
> other the unweighted output streams. Both have the vector streams overlaid<br>
> on them. The weighted output has a solid yellow background and I don't know<br>
> what that represents.<br>
<br>
In this case there is indeed no answer to "what looks better" because<br>
both look terrible, to be honest. Something terrible must have<br>
happened to the DEM, or there are indeed completely flat areas in the<br>
west. If the aim is to extract streams that resemble the digitized<br>
vector lines, the DEM you used seems to be not suitable. You could try<br>
to run r.stream.extract with the DEM as the only input and a larger<br>
threshold. The region should be aligned to the DEM first with g.region<br>
-p align=dem. If that produces reasonable results, you can try to<br>
fine-tune by adjusting the threshold and supplying a (weighed)<br>
accumulation map as additional input.<br>
<br>
Markus M<br>
<br>
PS: Please do not change the thread all the time. With regard to this<br>
subject, you have started now the 7th or so thread, I lost count,<br>
impossible for anyone to follow.<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 10 Apr 2012 10:01:11 +0200<br>
From: Markus Metz <<a href="mailto:markus.metz.giswork@googlemail.com">markus.metz.giswork@googlemail.com</a>><br>
Subject: Re: [GRASS-user] v.to.rast in grass70<br>
To: Paulo van Breugel <<a href="mailto:p.vanbreugel@gmail.com">p.vanbreugel@gmail.com</a>><br>
Cc: <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
Message-ID:<br>
<CAG+h=<a href="mailto:FFKW0U9XPbw%2B4L93K1gci2tMbU8o-B4umSXudUAUc18Uw@mail.gmail.com">FFKW0U9XPbw+4L93K1gci2tMbU8o-B4umSXudUAUc18Uw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
On Mon, Apr 9, 2012 at 12:08 PM, Paulo van Breugel<br>
<<a href="mailto:p.vanbreugel@gmail.com">p.vanbreugel@gmail.com</a>> wrote:<br>
> I can confirm there are issues in how GRASS 7 handles region settings. In my<br>
> case, when using v.to.rast, the region settings are reset to the default,<br>
> i.e., the raster is converted using the default region settings rather then<br>
> the current. I am now at revision 51177. In a previous revision (51117) I<br>
> had a problem with region settings when using r.patch, see Ticket #1636<br>
> (<a href="http://trac.osgeo.org/grass/ticket/1636" target="_blank">http://trac.osgeo.org/grass/ticket/1636</a>). These issues seem to be<br>
> independend of the mapset (it happens in both PERMANENT and other mapsets).<br>
><br>
<br>
I can not reproduce this problem. I tested in the North Carolina<br>
dataset, mapset user1.<br>
<br>
I displayed boundary_county, zoomed to the default region. Then I<br>
zoomed to a small part of the boundary_county vector and set the<br>
region from display. Then I zoomed to a different part of the vector<br>
that does not overlap with the current region. Thus the current<br>
region, the default region, and the display extents were all different<br>
from each other. Then I ran v.to.rast on the boundary_county vector.<br>
v.to.rast correctly used the current region and nothing was changed<br>
after v.to.rast was completed. This is r51306.<br>
<br>
Markus M<br>
<br>
> Paulo<br>
><br>
><br>
><br>
><br>
> On 04/09/2012 08:44 AM, Pankaj Kr Sharma wrote:<br>
><br>
> Dear Grass users and developers,<br>
><br>
> I just noticed that in grass70 , the v.to.rast module changes computational<br>
> region / current region<br>
> resolution to the display extent and resolution after finishing.<br>
><br>
> My maps are in "PERMANENT" mapset. And, I have tried a number of options ,<br>
> but cannot reset the resolution of raster thus created to the required<br>
> resolution.<br>
><br>
> Thanks.<br>
><br>
><br>
> _______________________________________________<br>
> grass-user mailing list<br>
> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
><br>
><br>
><br>
> _______________________________________________<br>
> grass-user mailing list<br>
> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
><br>
<br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
<br>
<br>
End of grass-user Digest, Vol 72, Issue 25<br>
******************************************<br>
</blockquote></div><br>