[GRASS-user] r.terraflow: water.cc:309: failed assertion `winp->getPriority() > prevWin.getPriority()'

[Xiaoyan Li]

I tried to use r.terraflow with SRTM data loaded on http://65.123.203.154:8099/grassweb/srtm.tif
but failed. Test with other SRTM data, same problem.

I used another data loaded on http://65.123.203.154:8099/grassweb/s1e00.tif
run r.terraflow successfully.

Is there anyone running into the same problem? Is it a bug in r.terraflow?


Platform: Mac OSx
GRASS Version is 6.1

GRASS 6.1.cvs (newloc):~ > g.region rast=s1
GRASS 6.1.cvs (newloc):~ > r.terraflow elev=s1 filled=s1filled direction=direction_s1 swatershed=sink_s1 accumulation=accum_s1 tci=tci_s1
r.terraflow December 2003
region size is 252 x 266
cell s1 header compatible with region header
elevation stored as FLOAT (4B) WARNING: raster s1 is of type CELL_TYPE --you should use r.terraflow.short

STREAM temporary files in /var/tmp  (THESE INTERMEDIATE STREAMS WILL NOT BE DELETED IN CASE OF ABNORMAL TERMINATION OF THE PROGRAM. TO SAVE SPACE PLEASE DELETE THESE FILES MANUALLY!)
file stats.out exists - renaming.
MFD flow direction
D8CUT=999999986991104.000000
memory size: 300.00M (0) bytes
Memory manager registering memory in MM_WARN_ON_MEMORY_EXCEEDED mode.
reading data from s1 to stream /var/tmp/STREAM_EOjs7e:   99%
total elements=67032, nodata elements=0
largest temporary files: 
                 FILL: 3.58M (0) [67032 elements, 56B each]
                 FLOW: 5.11M (0) [67032 elements, 80B each]
Will need at least 10.23M (0) space available in /var/tmp
------------------------------
COMPUTING FLOW DIRECTIONS
classifying nodata (inner & boundary)
finding nodata
relabeling nodata
merging relabeled grid
----------
assigning preliminary directions
----------
finding flat areas (plateaus and depressions)
plateau.cc:starting memscan
plateau.cc:memscan done
removing duplicate plateaus
relabeling plateaus
generating plateau statistics
----------
assigning directions on plateaus
--------------
generating watersheds and watershed graph
creating windows
warning: using slower scan
generateWatersheds
starting generate watersheds main loop
done with generate watersheds
creating windows
----------
flooding depressions
available memory: 299MB (314469302B)
warning: watershed 1 (R=1) not done
------------------------------
REASSIGNING DIRECTIONS
----------
finding flat areas (plateaus and depressions)
plateau.cc:starting memscan
plateau.cc:memscan done
removing duplicate plateaus
relabeling plateaus
generating plateau statistics
----------
assigning directions on plateaus
creating flowStream: [AMI_STREAM /var/tmp/flowStream 0]
compute flow directions done.
writing stream /var/tmp/STREAM_HsvfdD to cell file direction_s1:   99%
writing stream /var/tmp/STREAM_dI0qA4 to cell file s1filled:   99%
writing stream /var/tmp/STREAM_VtqrQd to cell file sink_s1:   99%
------------------------------
COMPUTING FLOW ACCUMULATION
creating sweep stream from fill output stream
sweep stream size: 5.11MB (67032 items, item size=80 B
 sorting sweep stream (5.11MB) in priority order
sorting sweep stream
sweeping: sweep:initialize flow data structure
FLOW_DATASTRUCTURE: adaptive pqueue
  99%
sorting sweep output stream
writing stream /var/tmp/STREAM_ohxZij to cell files accum_s1, tci_s1:   99%
r.terraflow done


----------------------------------------------
GRASS 6.1.cvs (newloc):~ > g.region rast=srtm
GRASS 6.1.cvs (newloc):~ > r.terraflow elev=srtm filled=srtmfilled direction=direction_srtm swatershed=sink_srtm accumulation=accum_srtm tci=tci_srtm
r.terraflow December 2003
region size is 420 x 638
cell srtm header compatible with region header
elevation stored as FLOAT (4B) 
STREAM temporary files in /var/tmp  (THESE INTERMEDIATE STREAMS WILL NOT BE DELETED IN CASE OF ABNORMAL TERMINATION OF THE PROGRAM. TO SAVE SPACE PLEASE DELETE THESE FILES MANUALLY!)
file stats.out exists - renaming.
MFD flow direction
D8CUT=999999986991104.000000
memory size: 300.00M (0) bytes
Memory manager registering memory in MM_WARN_ON_MEMORY_EXCEEDED mode.
reading data from srtm to stream /var/tmp/STREAM_On2FWU:   99%
total elements=267960, nodata elements=0
largest temporary files: 
                 FILL: 14.31M (0) [267960 elements, 56B each]
                 FLOW: 20.44M (0) [267960 elements, 80B each]
Will need at least 40.89M (0) space available in /var/tmp
------------------------------
COMPUTING FLOW DIRECTIONS
classifying nodata (inner & boundary)
finding nodata
relabeling nodata
merging relabeled grid
----------
assigning preliminary directions
----------
finding flat areas (plateaus and depressions)
plateau.cc:starting memscan
plateau.cc:memscan done
removing duplicate plateaus
relabeling plateaus
generating plateau statistics
----------
assigning directions on plateaus
--------------
generating watersheds and watershed graph
creating windows
warning: using slower scan
generateWatersheds
starting generate watersheds main loop
water.cc:309: failed assertion `winp->getPriority() > prevWin.getPriority()'
Abort trap



Xiaoyan Li