Nikos Alexandris wrote: > > > > r.cross fails on G7 while it works fine in G6.
> > > > G7: > > > > g.version -r > > > > GRASS 7.0.svn (2013) > > > > libgis Revision: 52468 > > > I tested this in another machine -- the same G7 revision. r.cross justs > > > "waits...." -- no action at all. Glynn Clements wrote: > > FWIW, it works for me; sort of. > > However: I note that it generates a distinct category for each cell > > which is null in all inputs. If the number of cells is large and the > > inputs contain many null cells, this could result in memory exhaustion > > and/or extreme computation times. > > I believe that the problem is due to an incomplete conversion from the > > 4.x G_get_map_row() function (which maps null to zero) to the 7.0 > > Rast_get_c_row() function (which leaves null as-is, i.e. -2^31). I.e. > > it now reads nulls as nulls (instead of zeros) but doesn't actually > > make any attempt to deal with the nulls which it gets. [..] > ...r.cross might "result in memory exhaustion and/or extreme computation > times" in G7, in a Location defined by > +proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +no_defs +a=6378137 > +rf=298.257222101 +towgs84=0,0,0,0,0,0,0 +to_meter=1 > Crossing maps composed of > - 263672500 and 263672500 > > Total cells in computational region are > 65914160 (100m res) Oh... now I see it: more important I guess is the number of NULL cells of the maps to be crossed! FWIW, they were/are (in the failing case): 61149574 and 63837480 Thanks again, Nikos
signature.asc
Description: This is a digitally signed message part.
_______________________________________________ grass-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-dev
