Dear all, I'm running r.stream.order stream_rast=stream direction=dir elevation=elv accumulation=flow strahler=order_strahler memory=15000 --o on a quite large stream raster 25200 x 25200 pixels.
gdalinfo -mm Computed Min/Max=11001.000,52594532.000 I'm giving 15G as memory availability so all the raster can be in the RAM. The computation proceeds with errors, nonetheless the output presents some streams that are labeled as -1 rather than having the stream order number. If I crop the stream raster input to 2400 x 6000 pixels, gdalinfo -mm Computed Min/Max=7632381.000,52594532.000 then output stream order results are mapped correctly. Seems to me that in case of the large raster r.stream.order hits some data type limit that produces the output error. I have tried to look into the r.stream.order code but I did not see where the error can be. Any suggestions? Thanks Giuseppe -- Giuseppe Amatulli, Ph.D. Research scientist at School of the Environment Yale University New Haven, CT, USA - 06511 Tweeter: @BigDataEcology Teaching: http://spatial-ecology.net Work: https://environment.yale.edu/profile/giuseppe-amatulli/
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user