[GRASS-dev] wxGUI GRASS 7: how to drape a map over a DEM?

hi,

I fail to drape a map over a DEM in the wxGUI of GRASS 7... (Linux).

Any pointers how to achieve that? The view is just empty...

thanks
Markus

On Fri, May 24, 2013 at 1:03 PM, Markus Neteler <neteler@osgeo.org> wrote:

hi,

I fail to drape a map over a DEM in the wxGUI of GRASS 7... (Linux).

Any pointers how to achieve that? The view is just empty...

So the problem is that wxNviz is not working? Please read NOTE on man page
[1] and follow the instructions. Otherwise I have no idea what's wrong.

Anna

[1] http://grass.osgeo.org/grass70/manuals/wxGUI.nviz.html

thanks
Markus
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

On Sat, May 25, 2013 at 2:17 PM, Anna Petrášová <kratochanna@gmail.com> wrote:

On Fri, May 24, 2013 at 1:03 PM, Markus Neteler <neteler@osgeo.org> wrote:

hi,

I fail to drape a map over a DEM in the wxGUI of GRASS 7... (Linux).

Any pointers how to achieve that? The view is just empty...

So the problem is that wxNviz is not working? Please read NOTE on man page
[1] and follow the instructions.

Oh, thanks, I overlooked this hint. Added to
http://grasswiki.osgeo.org/wiki/WxNviz#Troubleshooting

Will try again,
Markus

Otherwise I have no idea what's wrong.

Anna

[1] http://grass.osgeo.org/grass70/manuals/wxGUI.nviz.html

thanks
Markus
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

On Sat, May 25, 2013 at 5:47 PM, Markus Neteler <neteler@osgeo.org> wrote:

On Sat, May 25, 2013 at 2:17 PM, Anna Petrášová <kratochanna@gmail.com> wrote:

On Fri, May 24, 2013 at 1:03 PM, Markus Neteler <neteler@osgeo.org> wrote:

hi,

I fail to drape a map over a DEM in the wxGUI of GRASS 7... (Linux).

Any pointers how to achieve that? The view is just empty...

So the problem is that wxNviz is not working? Please read NOTE on man page
[1] and follow the instructions.

Oh, thanks, I overlooked this hint. Added to
http://grasswiki.osgeo.org/wiki/WxNviz#Troubleshooting

Will try again,

... no change... nothing is drawn.

The NC exercise is:

# set computational region to map
g.region rast=elev_lid792_1m -p

r.watershed elevation=elev_lid792_1m accumulation=flowacc

- load both elev_lid792_1m and flowacc into Layer Manager
- 3D view: wxNVIZ

... (for me) it does not show the perspective view.

Markus

On Sat, May 25, 2013 at 10:25 PM, Markus Neteler <neteler@osgeo.org> wrote:

On Sat, May 25, 2013 at 5:47 PM, Markus Neteler <neteler@osgeo.org> wrote:
> On Sat, May 25, 2013 at 2:17 PM, Anna Petrášová <kratochanna@gmail.com>
wrote:
>>
>>
>>
>> On Fri, May 24, 2013 at 1:03 PM, Markus Neteler <neteler@osgeo.org>
wrote:
>>>
>>> hi,
>>>
>>> I fail to drape a map over a DEM in the wxGUI of GRASS 7... (Linux).
>>>
>>> Any pointers how to achieve that? The view is just empty...
>>
>>
>> So the problem is that wxNviz is not working? Please read NOTE on man
page
>> [1] and follow the instructions.
>
> Oh, thanks, I overlooked this hint. Added to
> http://grasswiki.osgeo.org/wiki/WxNviz#Troubleshooting
>
> Will try again,

... no change... nothing is drawn.

The NC exercise is:

# set computational region to map
g.region rast=elev_lid792_1m -p

r.watershed elevation=elev_lid792_1m accumulation=flowacc

- load both elev_lid792_1m and flowacc into Layer Manager
- 3D view: wxNVIZ

... (for me) it does not show the perspective view.

ok, I get it. You shouldn't load flowacc as a layer because than it is
displayed as another surface which we don't want in this case. Or when it
is loaded (when 2D) just uncheck it and adjust the view to see the DEM.
Then, go to the Data -> Surface -> Surface attributes and change color to
flowacc. This should be the same concept as for the old nviz or not? BTW,
it's in the video tutorial:)

The problem why you don't see anything when you load both maps is that
flowacc has large extreme values, so the initial view is confused. Anyway
you should be able see something when you adjust the view to some normal
height.

Anna

On Sun, May 26, 2013 at 10:32 AM, Anna Petrášová <kratochanna@gmail.com> wrote:
...

ok, I get it. You shouldn't load flowacc as a layer because than it is
displayed as another surface which we don't want in this case. Or when it is
loaded (when 2D) just uncheck it and adjust the view to see the DEM.

OK - since I come in the exercise from the 2D view, this may do the trick.
So, in essence, all maps active in the Layer Manager will be treated
as elevation layers in wxNVIZ (?).

Then,
go to the Data -> Surface -> Surface attributes and change color to flowacc.
This should be the same concept as for the old nviz or not? BTW, it's in
the video tutorial:)

The old NVIZ I rigorously started from CMD line :slight_smile:

nviz elev=elevation colo=flowacc

The problem why you don't see anything when you load both maps is that
flowacc has large extreme values, so the initial view is confused.

Indeed, I am surprised (extra unrelated issue) that flowacc has this
range:

r.univar flowacc
100%
total null and non-null cells: 525000
total null cells: 0

Of the non-null cells:
----------------------
n: 525000
minimum: -253898
maximum: 15278.2
range: 269176
mean: -275.776
mean of absolute values: 311.452
standard deviation: 4199.2
variance: 1.76333e+07
variation coefficient: -1522.69 %
sum: -144782397.558392

Anyway
you should be able see something when you adjust the view to some normal
height.

Thanks, Anna. I'll get to speed with it for sure.

Markus

On Sun, May 26, 2013 at 10:50 AM, Markus Neteler <neteler@osgeo.org> wrote:

Indeed, I am surprised (extra unrelated issue) that flowacc has this
range:

r.univar flowacc
100%
total null and non-null cells: 525000
total null cells: 0

Of the non-null cells:
----------------------
n: 525000
minimum: -253898
maximum: 15278.2

Why is this a surprise?

Note that "The absolute value of each cell in this output map layer is
the amount of overland flow that traverses the cell. This value will
be the number of upland cells plus one [...]. Negative numbers
indicate that those cells possibly have surface runoff from outside of
the current geographic region. Thus, any cells with negative values
cannot have their surface runoff and sedimentation yields calculated
accurately."

Markus M

range: 269176
mean: -275.776
mean of absolute values: 311.452
standard deviation: 4199.2
variance: 1.76333e+07
variation coefficient: -1522.69 %
sum: -144782397.558392

Anyway
you should be able see something when you adjust the view to some normal
height.

Thanks, Anna. I'll get to speed with it for sure.

Markus
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

On Sun, May 26, 2013 at 10:50 AM, Markus Neteler <neteler@osgeo.org> wrote:

On Sun, May 26, 2013 at 10:32 AM, Anna Petrášová <kratochanna@gmail.com>
wrote:
...
> ok, I get it. You shouldn't load flowacc as a layer because than it is
> displayed as another surface which we don't want in this case. Or when
it is
> loaded (when 2D) just uncheck it and adjust the view to see the DEM.

OK - since I come in the exercise from the 2D view, this may do the trick.
So, in essence, all maps active in the Layer Manager will be treated
as elevation layers in wxNVIZ (?).

yes, loaded raster map = surface

> Then,
> go to the Data -> Surface -> Surface attributes and change color to
flowacc.
> This should be the same concept as for the old nviz or not? BTW, it's in
> the video tutorial:)

The old NVIZ I rigorously started from CMD line :slight_smile:

nviz elev=elevation colo=flowacc

> The problem why you don't see anything when you load both maps is that
> flowacc has large extreme values, so the initial view is confused.

Indeed, I am surprised (extra unrelated issue) that flowacc has this
range:

r.univar flowacc
100%
total null and non-null cells: 525000
total null cells: 0

Of the non-null cells:
----------------------
n: 525000
minimum: -253898
maximum: 15278.2
range: 269176
mean: -275.776
mean of absolute values: 311.452
standard deviation: 4199.2
variance: 1.76333e+07
variation coefficient: -1522.69 %
sum: -144782397.558392

> Anyway
> you should be able see something when you adjust the view to some normal
> height.

Thanks, Anna. I'll get to speed with it for sure.

you're welcome

Markus

On Sun, May 26, 2013 at 10:32 AM, Anna Petrášová <kratochanna@gmail.com> wrote:

The problem why you don't see anything when you load both maps is that
flowacc has large extreme values, so the initial view is confused. Anyway
you should be able see something when you adjust the view to some normal
height.

That means I can not change the raster map to be used as surface?
Because I see that two different surfaces are used after updating,
resetting, and changing height to a reasonable value, even though only
one surface map is given. BTW, nviz in 6.x accepts multiple raster
maps as surfaces.

Markus M

On Sun, May 26, 2013 at 11:08 AM, Markus Metz <markus.metz.giswork@gmail.com

wrote:

On Sun, May 26, 2013 at 10:32 AM, Anna Petrášová <kratochanna@gmail.com>
wrote:
>
> The problem why you don't see anything when you load both maps is that
> flowacc has large extreme values, so the initial view is confused. Anyway
> you should be able see something when you adjust the view to some normal
> height.

That means I can not change the raster map to be used as surface?
Because I see that two different surfaces are used after updating,
resetting, and changing height to a reasonable value, even though only
one surface map is given. BTW, nviz in 6.x accepts multiple raster
maps as surfaces.

Sorry, I don't really understand what you mean. You can display more

surfaces, you can unload (uncheck or remove) surface and add new.

Anna

On Sun, May 26, 2013 at 11:21 AM, Anna Petrášová <kratochanna@gmail.com> wrote:

On Sun, May 26, 2013 at 11:08 AM, Markus Metz
<markus.metz.giswork@gmail.com> wrote:

On Sun, May 26, 2013 at 10:32 AM, Anna Petrášová <kratochanna@gmail.com>
wrote:
>
> The problem why you don't see anything when you load both maps is that
> flowacc has large extreme values, so the initial view is confused.
> Anyway
> you should be able see something when you adjust the view to some normal
> height.

That means I can not change the raster map to be used as surface?
Because I see that two different surfaces are used after updating,
resetting, and changing height to a reasonable value, even though only
one surface map is given. BTW, nviz in 6.x accepts multiple raster
maps as surfaces.

Sorry, I don't really understand what you mean. You can display more
surfaces, you can unload (uncheck or remove) surface and add new.

OK, got it. All active raster maps in the layer manager are loaded as surfaces.

Markus M

I have updated the videos for wxnviz including working with multiple surfaces, color draping and adding a constant surface
and cutting planes here:
http://courses.ncsu.edu/mea582/common/GIS_anal_grass/GIS_Anal_grvisual.html
you can see there wxnviz and old nviz side-by-side.

I will re-do the videos once we have a final release of 6.4.3

The only thing that you need to be aware of is that in GRASS6 change in region is ignored by wxnviz - it always
uses the region which was set when you first switched to it. You need to quit wxgui, change region and restart wxgui.
In GRASS7 region can be changed - one more reason to have GRASS7 pre-release.

Anna please correct me if I am wrong in anything above.

I hope this helps, Helena

Helena Mitasova
Associate Professor
Department of Marine, Earth, and Atmospheric Sciences
2800 Faucette Drive, Rm. 1125 Jordan Hall
North Carolina State University
Raleigh, NC 27695-8208
hmitaso@ncsu.edu

"All electronic mail messages in connection with State business which are sent to or received by this account are subject to the NC Public Records Law and may be disclosed to third parties.”

On May 26, 2013, at 7:16 AM, Markus Metz wrote:

On Sun, May 26, 2013 at 11:21 AM, Anna Petrášová <kratochanna@gmail.com> wrote:

On Sun, May 26, 2013 at 11:08 AM, Markus Metz
<markus.metz.giswork@gmail.com> wrote:

On Sun, May 26, 2013 at 10:32 AM, Anna Petrášová <kratochanna@gmail.com>
wrote:

The problem why you don't see anything when you load both maps is that
flowacc has large extreme values, so the initial view is confused.
Anyway
you should be able see something when you adjust the view to some normal
height.

That means I can not change the raster map to be used as surface?
Because I see that two different surfaces are used after updating,
resetting, and changing height to a reasonable value, even though only
one surface map is given. BTW, nviz in 6.x accepts multiple raster
maps as surfaces.

Sorry, I don't really understand what you mean. You can display more
surfaces, you can unload (uncheck or remove) surface and add new.

OK, got it. All active raster maps in the layer manager are loaded as surfaces.

Markus M
_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

On Sun, May 26, 2013 at 3:47 PM, Helena Mitasova <hmitaso@ncsu.edu> wrote:

I have updated the videos for wxnviz including working with multiple
surfaces, color draping and adding a constant surface
and cutting planes here:
http://courses.ncsu.edu/mea582/common/GIS_anal_grass/GIS_Anal_grvisual.html
you can see there wxnviz and old nviz side-by-side.

I will re-do the videos once we have a final release of 6.4.3

The only thing that you need to be aware of is that in GRASS6 change in
region is ignored by wxnviz - it always
uses the region which was set when you first switched to it. You need to
quit wxgui, change region and restart wxgui.
In GRASS7 region can be changed - one more reason to have GRASS7
pre-release.

Anna please correct me if I am wrong in anything above.

Agreed. As far as I remember it was not easy to backport the relevant part
of code in core library, at least for me.

Anna

I hope this helps, Helena

Helena Mitasova
Associate Professor
Department of Marine, Earth, and Atmospheric Sciences
2800 Faucette Drive, Rm. 1125 Jordan Hall
North Carolina State University
Raleigh, NC 27695-8208
hmitaso@ncsu.edu

"All electronic mail messages in connection with State business which are
sent to or received by this account are subject to the NC Public Records
Law and may be disclosed to third parties.”

On May 26, 2013, at 7:16 AM, Markus Metz wrote:

> On Sun, May 26, 2013 at 11:21 AM, Anna Petrášová <kratochanna@gmail.com>
wrote:
>>
>>
>>
>> On Sun, May 26, 2013 at 11:08 AM, Markus Metz
>> <markus.metz.giswork@gmail.com> wrote:
>>>
>>> On Sun, May 26, 2013 at 10:32 AM, Anna Petrášová <
kratochanna@gmail.com>
>>> wrote:
>>>>
>>>> The problem why you don't see anything when you load both maps is that
>>>> flowacc has large extreme values, so the initial view is confused.
>>>> Anyway
>>>> you should be able see something when you adjust the view to some
normal
>>>> height.
>>>
>>> That means I can not change the raster map to be used as surface?
>>> Because I see that two different surfaces are used after updating,
>>> resetting, and changing height to a reasonable value, even though only
>>> one surface map is given. BTW, nviz in 6.x accepts multiple raster
>>> maps as surfaces.
>>>
>> Sorry, I don't really understand what you mean. You can display more
>> surfaces, you can unload (uncheck or remove) surface and add new.
>
> OK, got it. All active raster maps in the layer manager are loaded as
surfaces.
>
> Markus M
> _______________________________________________
> grass-dev mailing list
> grass-dev@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-dev