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. Some attempts that _work_ on G7, are projected in +proj=utm +no_defs +zone=34 +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1 Crossing maps composed of - 182402720 and 165183040 cells or - 165183040 and 165183040 Total cells in computational region are - 458870596 (30m res) - 165234916 (50m res) But, 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) Hope this info is of some use. Best, Nikos _______________________________________________ grass-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-dev
