Tue, Feb 7 2006
17:54:42
|
|
Request created by guest
|
|
Subject: r.carve: output identical to input and segfault when using points=
Platform: GNU/Linux/x86
grass binary for platform: Compiled from Sources
GRASS Version: CVS 2006.02.07
$ r.carve rast=nnbathy_D_no vect=cieki_atlas_patched output=nnbathy_D_no_carved
points=nnbathy_D_no_carved_pkt width=10.000000 depth=0.1 -n
Reading raster file...
100%
Processing lines...Segmentation fault
More info needed to sort it out?
But without option points= it seems to work, despite the warning:
$ r.carve rast=nnbathy_D_no vect=cieki_atlas_patched output=nnbathy_D_no_carved
width=10.000000 depth=0.1 -n Reading raster file...
100%
Processing lines...WARNING: Vect runs out of region and re-enters - this case
is not yet
implemented!
Moreover, the r.carve output is _identical_ to input. Substracting the r.carve
output DEM from original DEM gives 0.0. Looks like r.carve is not doing anything...
too bad.
Has anydoby succeeded using it? Maybe there is something special about my case.
I will keep on trying tomorrow with different input and any hints are wellocme.
Maciek |
|
Tue, Feb 7 2006
18:07:02
|
|
Mail sent by msieczka
|
|
I have just tried with the same input DEM but very simple, single vector line
indicating input streamline. I made sure it doesn't exceed the input DEM
extent this time.
$ r.carve rast=nnbathy_D_no vect=carveit output=nnbathy_D_no_carved
points=nnbathy_D_no_carved_pkt width=10.000000 depth=0.1 -n Reading raster file...
100%
Processing lines...
Writing raster file... done.
Building topology ...
0 primitives registered
0 areas built
0 isles built
Attaching islands:
Attaching centroids: Topology was built.
Number of nodes : 0
Number of primitives: 0
Number of points : 0
Number of lines : 0
Number of boundaries: 0
Number of centroids : 0
Number of areas : 0
Number of isles : 0
r.carve complete. r.carve -n rast="nnbathy_D_no" vect="carveit"
output="nnbathy_D_no_carved" points="nnbathy_D_no_carved_pkt" width=10.000000
depth=0.1
No segfault this time.
Still, the output DEM is _identical_ to input DEM. And the ouput points= has
0
geometry.
No good. Any chances it will be fixed in close future?
Maciek
|
|
Sun, Feb 19 2006
17:35:56
|
|
Mail sent by mneteler
|
|
For now, r.carve is disabled again. Please contact Brad
for details.
Markus |
|
Mon, May 15 2006
15:24:25
|
|
Mail sent by mneteler
|
|
Hi,
does the problem persist?
Markus |
|
Mon, May 15 2006
18:08:12
|
|
Mail sent by msieczka
|
|
Markus wrote:
> does the problem persist?
Markus,
I'm in touch with Brad, who works on the issue, as time allows it to him.
The segfault is fixed, but there are still problems with the vector output:
GRASS 6.1.cvs (caves_utm33):~ > r.carve -n rast=hipso_natneigh_validvect=cieki03
output=hipso_natneigh_valid_carved
points=hipso_natneigh_valid_carve_points depth=0.5
Reading raster file... 100%
Processing lines... 98%
Writing raster file... done.
Building topology ... Registering lines:
r.carve: index.c:157: RTreeInsertRect: Assertion `r->boundary[i] <=
r->boundary[3 +i]' failed. Aborted [Raster MASK present]
Brad said he will take care of it when he can.
Best,
Maciek
|
|
Wed, Jul 26 2006
18:19:52
|
|
User changed to tutey@o2.pl by msieczka
|
|