[GRASS-user] Re: r.terraflow: water.cc:309: failed assertion

[Xiaoyan Li]

These links should work:
http://laits.gmu.edu/testpages/s1e00.tif
http://laits.gmu.edu/testpages/stream.tif


Thanks!

Xiaoyan Li
> ------------------------------
> Message: 6
> Date: Fri, 26 Oct 2007 14:07:17 -0400
> From: Xiaoyan Li <xlia at gmu.edu>
> Subject: [GRASS-user] r.terraflow: water.cc:309: failed assertion
> 	`winp->getPriority() > prevWin.getPriority()'
> To: grassuser at grass.itc.it
> Cc: hmitaso at ncsu.edu
> Message-ID: <f520c95626635.4721f495 at gmu.edu>
> Content-Type: text/plain; charset=us-ascii
> 
> 
> I tried to use r.terraflow with SRTM data loaded on 
> http://65.123.203.154:8099/grassweb/srtm.tifbut failed. Test with 
> other SRTM data, same problem.
> 
> I used another data loaded on 
> http://65.123.203.154:8099/grassweb/s1e00.tifrun 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_s1r.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
>