[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