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

Sören Gebbert soerengebbert at googlemail.com
Wed Sep 21 16:02:23 EDT 2011


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 at 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 at 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 at 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 at 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 at 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 at 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 at 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.
>


More information about the grass-user mailing list