[GRASS-user] problems with r.to.rast3elev

Hello all,

I’m experiencing some issues with r.to.rast3elev , that I didn’t had before. I’m trying to reprocess some volume calculation
on dune fields in southern Brazil (see this: http://geomorphometry.org/Grohmann2011) so I can evolve what I presented in
the Geomorphometry Meeting into a full paper, but at this point I can’t calculate the volumes any more!

(BTW, Helena, I’m thinking about comparing the volumes from r.volume and r.to.rast3elev, since I think there will be differences…)

OS: Linux, Ubuntu 11.04
In GRASS 6.4.1, installed from Ubuntu repositories, I got this (running inside a python session):

grass.run_command(‘r.to.rast3elev’, input=clump, elevation=dem, output=dem3d, lower=1, overwrite=True)
Creating 3D raster map
ERROR: G3d_cache_hash_remove_name: name not in hashtable

and in GRASS 7.0svn updated today, I got a segmentation fault.

the region is set to match the ‘clump’ map (which is the same as the mask, to limit the calculations to the active dunes).

any help is appreciated

Carlos


Prof. Carlos Henrique Grohmann - Geologist D.Sc.
Institute of Geosciences - Univ. of São Paulo, Brazil

http://www.igc.usp.br/pessoais/guano
http://digitalelevation.wordpress.com/
http://lattes.cnpq.br/5846052449613692 (CV)

Twitter: @CarlosGrohmann
http://carlosgrohmann.tumblr.com/
Linux User #89721


Can’t stop the signal.

Hi Carlos,
can you please provide more information about the grass7 segfault?

Can you please print your region settings (g.region -p3)?

Can you please use gdb or valgrind to check at which line in the code
the segfault appears?

I will try to fix this issue in grass7, but i cant reproduce it on my system.

Best regards
Soeren

2011/9/20 Carlos Grohmann <carlos.grohmann@gmail.com>:

Hello all,

I'm experiencing some issues with r.to.rast3elev , that I didn't had before.
I'm trying to reprocess some volume calculation
on dune fields in southern Brazil (see this:
http://geomorphometry.org/Grohmann2011) so I can evolve what I presented in
the Geomorphometry Meeting into a full paper, but at this point I can't
calculate the volumes any more!

(BTW, Helena, I'm thinking about comparing the volumes from r.volume and
r.to.rast3elev, since I think there will be differences..)

OS: Linux, Ubuntu 11.04
In GRASS 6.4.1, installed from Ubuntu repositories, I got this (running
inside a python session):

grass.run_command('r.to.rast3elev', input=clump, elevation=dem,
output=dem3d, lower=1, overwrite=True)
Creating 3D raster map
ERROR: G3d_cache_hash_remove_name: name not in hashtable

and in GRASS 7.0svn updated today, I got a segmentation fault.

the region is set to match the 'clump' map (which is the same as the mask,
to limit the calculations to the active dunes).

any help is appreciated

Carlos

--
Prof. Carlos Henrique Grohmann - Geologist D.Sc.
Institute of Geosciences - Univ. of São Paulo, Brazil
---
http://www.igc.usp.br/pessoais/guano
http://digitalelevation.wordpress.com/
http://lattes.cnpq.br/5846052449613692 (CV)
---
Twitter: @CarlosGrohmann
http://carlosgrohmann.tumblr.com/
Linux User #89721
________________
Can’t stop the signal.

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

Thanks Soeren,

As a side note, it did run on my Mac, using Kygschaos package (6.4.1, I guess)…

here is the output if g.region -p3:

GRASS 7.0.svn (santa_catarina@dunas2):~ > g.region -p3
projection: 1 (UTM)
zone: -22
datum: wgs84
ellipsoid: wgs84
north: 6901400
south: 6899300
west: 731900
east: 732900
top: 60.00000000
bottom: 0.00000000
nsres: 0.5
nsres3: 0.5
ewres: 0.5
ewres3: 0.5
tbres: 0.5
rows: 4200
rows3: 4200
cols: 2000
cols3: 2000
depths: 120
cells: 8400000
cells3: 1008000000

and here of gdb:

GRASS 7.0.svn (santa_catarina@dunas2):~ > r.to.rast3elev --overwrite input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic output=garopaba_vol05m lower=1
Creating 3D raster map
Segmentation fault
GRASS 7.0.svn (santa_catarina@dunas2):~ > gdb which r.to.rast3elev
GNU gdb (Ubuntu/Linaro 7.2-1ubuntu11) 7.2
Copyright (C) 2010 Free Software Foundation, Inc.
License GPLv3+: GNU GPL version 3 or later <http://gnu.org/licenses/gpl.html>
This is free software: you are free to change and redistribute it.
There is NO WARRANTY, to the extent permitted by law. Type “show copying”
and “show warranty” for details.
This GDB was configured as “i686-linux-gnu”.
For bug reporting instructions, please see:
<http://www.gnu.org/software/gdb/bugs/>…
Reading symbols from /usr/local/grass-7.0.svn/bin/r.to.rast3elev…(no debugging symbols found)…done.
(gdb) run input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic output=garopaba_vol05m lower=1
Starting program: /usr/local/grass-7.0.svn/bin/r.to.rast3elev input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic output=garopaba_vol05m lower=1
Creating 3D raster map
0%
Program received signal SIGSEGV, Segmentation fault.
0xb7fb5ae4 in Rast_set_d_null_value () from /usr/local/grass-7.0.svn/lib/libgrass_raster.7.0.svn.so
(gdb) bt full
#0 0xb7fb5ae4 in Rast_set_d_null_value ()
from /usr/local/grass-7.0.svn/lib/libgrass_raster.7.0.svn.so
No symbol table info available.
#1 0xb7fd3611 in Rast3d_set_null_value ()
from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#2 0xb7fd6a92 in Rast3d_set_null_tile_type ()
from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#3 0xb7fd6e3e in Rast3d_read_tile () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#4 0xb7fc9288 in ?? () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#5 0xb7fc89fb in Rast3d_cache_elt_ptr () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#6 0xb7fd642d in Rast3d_get_tile_ptr () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#7 0xb7fd4761 in Rast3d_put_double () from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#8 0x08049880 in elev_raster_to_g3d ()
No symbol table info available.
#9 0x08049e6e in main ()
No symbol table info available.
(gdb) l
No symbol table is loaded. Use the “file” command.
(gdb) frame 2
#2 0xb7fd6a92 in Rast3d_set_null_tile_type ()
from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
(gdb) l
No symbol table is loaded. Use the “file” command.
(gdb)

best,

Carlos

2011/9/21 Sören Gebbert <soerengebbert@googlemail.com>

Hi Carlos,
can you please provide more information about the grass7 segfault?

Can you please print your region settings (g.region -p3)?

Can you please use gdb or valgrind to check at which line in the code
the segfault appears?

I will try to fix this issue in grass7, but i cant reproduce it on my system.

Best regards
Soeren

2011/9/20 Carlos Grohmann <carlos.grohmann@gmail.com>:

Hello all,

I’m experiencing some issues with r.to.rast3elev , that I didn’t had before.
I’m trying to reprocess some volume calculation
on dune fields in southern Brazil (see this:
http://geomorphometry.org/Grohmann2011) so I can evolve what I presented in
the Geomorphometry Meeting into a full paper, but at this point I can’t
calculate the volumes any more!

(BTW, Helena, I’m thinking about comparing the volumes from r.volume and
r.to.rast3elev, since I think there will be differences…)

OS: Linux, Ubuntu 11.04
In GRASS 6.4.1, installed from Ubuntu repositories, I got this (running
inside a python session):

grass.run_command(‘r.to.rast3elev’, input=clump, elevation=dem,
output=dem3d, lower=1, overwrite=True)
Creating 3D raster map
ERROR: G3d_cache_hash_remove_name: name not in hashtable

and in GRASS 7.0svn updated today, I got a segmentation fault.

the region is set to match the ‘clump’ map (which is the same as the mask,
to limit the calculations to the active dunes).

any help is appreciated

Carlos


Prof. Carlos Henrique Grohmann - Geologist D.Sc.
Institute of Geosciences - Univ. of São Paulo, Brazil

http://www.igc.usp.br/pessoais/guano
http://digitalelevation.wordpress.com/
http://lattes.cnpq.br/5846052449613692 (CV)

Twitter: @CarlosGrohmann
http://carlosgrohmann.tumblr.com/
Linux User #89721


Can’t stop the signal.


grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Prof. Carlos Henrique Grohmann - Geologist D.Sc.
Institute of Geosciences - Univ. of São Paulo, Brazil

http://www.igc.usp.br/pessoais/guano
http://digitalelevation.wordpress.com/
http://lattes.cnpq.br/5846052449613692 (CV)

Twitter: @CarlosGrohmann
http://carlosgrohmann.tumblr.com/
Linux User #89721


Can’t stop the signal.

Hi Carlos,
please try svn revision 48398 of grass7.

I think the reason was the wrong computation of the tile size, because
the chosen variable types are not large enough. This should be fixed
now.

I hope it works now.

Best regards
Soeren

2011/9/21 Carlos Grohmann <carlos.grohmann@gmail.com>:

Thanks Soeren,

As a side note, it did run on my Mac, using Kygschaos package (6.4.1, I
guess)...

here is the output if g.region -p3:

GRASS 7.0.svn (santa_catarina@dunas2):~ > g.region -p3
projection: 1 (UTM)
zone: -22
datum: wgs84
ellipsoid: wgs84
north: 6901400
south: 6899300
west: 731900
east: 732900
top: 60.00000000
bottom: 0.00000000
nsres: 0.5
nsres3: 0.5
ewres: 0.5
ewres3: 0.5
tbres: 0.5
rows: 4200
rows3: 4200
cols: 2000
cols3: 2000
depths: 120
cells: 8400000
cells3: 1008000000

and here of gdb:

GRASS 7.0.svn (santa_catarina@dunas2):~ > r.to.rast3elev --overwrite
input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic
output=garopaba_vol05m lower=1
Creating 3D raster map
Segmentation fault
GRASS 7.0.svn (santa_catarina@dunas2):~ > gdb `which r.to.rast3elev`
GNU gdb (Ubuntu/Linaro 7.2-1ubuntu11) 7.2
Copyright (C) 2010 Free Software Foundation, Inc.
License GPLv3+: GNU GPL version 3 or later
<http://gnu.org/licenses/gpl.html&gt;
This is free software: you are free to change and redistribute it.
There is NO WARRANTY, to the extent permitted by law. Type "show copying"
and "show warranty" for details.
This GDB was configured as "i686-linux-gnu".
For bug reporting instructions, please see:
<http://www.gnu.org/software/gdb/bugs/&gt;\.\.\.
Reading symbols from /usr/local/grass-7.0.svn/bin/r.to.rast3elev...(no
debugging symbols found)...done.
(gdb) run input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic
output=garopaba_vol05m lower=1
Starting program: /usr/local/grass-7.0.svn/bin/r.to.rast3elev
input=mask_volume_garopaba elevation=garopaba_SOLO_bicubic
output=garopaba_vol05m lower=1
Creating 3D raster map
0%
Program received signal SIGSEGV, Segmentation fault.
0xb7fb5ae4 in Rast_set_d_null_value () from
/usr/local/grass-7.0.svn/lib/libgrass_raster.7.0.svn.so
(gdb) bt full
#0 0xb7fb5ae4 in Rast_set_d_null_value ()
from /usr/local/grass-7.0.svn/lib/libgrass_raster.7.0.svn.so
No symbol table info available.
#1 0xb7fd3611 in Rast3d_set_null_value ()
from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#2 0xb7fd6a92 in Rast3d_set_null_tile_type ()
from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#3 0xb7fd6e3e in Rast3d_read_tile () from
/usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#4 0xb7fc9288 in ?? () from
/usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#5 0xb7fc89fb in Rast3d_cache_elt_ptr () from
/usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#6 0xb7fd642d in Rast3d_get_tile_ptr () from
/usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#7 0xb7fd4761 in Rast3d_put_double () from
/usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
No symbol table info available.
#8 0x08049880 in elev_raster_to_g3d ()
No symbol table info available.
#9 0x08049e6e in main ()
No symbol table info available.
(gdb) l
No symbol table is loaded. Use the "file" command.
(gdb) frame 2
#2 0xb7fd6a92 in Rast3d_set_null_tile_type ()
from /usr/local/grass-7.0.svn/lib/libgrass_g3d.7.0.svn.so
(gdb) l
No symbol table is loaded. Use the "file" command.
(gdb)

best,

Carlos

2011/9/21 Sören Gebbert <soerengebbert@googlemail.com>

Hi Carlos,
can you please provide more information about the grass7 segfault?

Can you please print your region settings (g.region -p3)?

Can you please use gdb or valgrind to check at which line in the code
the segfault appears?

I will try to fix this issue in grass7, but i cant reproduce it on my
system.

Best regards
Soeren

2011/9/20 Carlos Grohmann <carlos.grohmann@gmail.com>:
> Hello all,
>
> I'm experiencing some issues with r.to.rast3elev , that I didn't had
> before.
> I'm trying to reprocess some volume calculation
> on dune fields in southern Brazil (see this:
> http://geomorphometry.org/Grohmann2011) so I can evolve what I presented
> in
> the Geomorphometry Meeting into a full paper, but at this point I can't
> calculate the volumes any more!
>
> (BTW, Helena, I'm thinking about comparing the volumes from r.volume and
> r.to.rast3elev, since I think there will be differences..)
>
> OS: Linux, Ubuntu 11.04
> In GRASS 6.4.1, installed from Ubuntu repositories, I got this (running
> inside a python session):
>
> grass.run_command('r.to.rast3elev', input=clump, elevation=dem,
> output=dem3d, lower=1, overwrite=True)
> Creating 3D raster map
> ERROR: G3d_cache_hash_remove_name: name not in hashtable
>
> and in GRASS 7.0svn updated today, I got a segmentation fault.
>
> the region is set to match the 'clump' map (which is the same as the
> mask,
> to limit the calculations to the active dunes).
>
> any help is appreciated
>
> Carlos
>
>
>
>
> --
> Prof. Carlos Henrique Grohmann - Geologist D.Sc.
> Institute of Geosciences - Univ. of São Paulo, Brazil
> ---
> http://www.igc.usp.br/pessoais/guano
> http://digitalelevation.wordpress.com/
> http://lattes.cnpq.br/5846052449613692 (CV)
> ---
> Twitter: @CarlosGrohmann
> http://carlosgrohmann.tumblr.com/
> Linux User #89721
> ________________
> Can’t stop the signal.
>
> _______________________________________________
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>
>

--
Prof. Carlos Henrique Grohmann - Geologist D.Sc.
Institute of Geosciences - Univ. of São Paulo, Brazil
---
http://www.igc.usp.br/pessoais/guano
http://digitalelevation.wordpress.com/
http://lattes.cnpq.br/5846052449613692 (CV)
---
Twitter: @CarlosGrohmann
http://carlosgrohmann.tumblr.com/
Linux User #89721
________________
Can’t stop the signal.