[GRASS-user] spgrass6

2012-02-22 Thread Paolo Cavallini
Hi all. I'm working with spgrass6 R module, and I found that spplot() is horribly slow (several *minutes* for a 5k cell raster). Is this a common, known problem, or a local one? Thanks. -- Paolo Cavallini See: http://www.faunalia.it/pc ___

Re: [GRASS-user] spgrass6

2012-02-22 Thread Markus Neteler
On Wed, Feb 22, 2012 at 9:17 AM, Paolo Cavallini cavall...@faunalia.it wrote: Hi all. I'm working with spgrass6 R module, and I found that spplot() is horribly slow (several *minutes* for a 5k cell raster). Is this a common, known problem, or a local one? I am afraid that it is a generic R

Re: [GRASS-user] spgrass6

2012-02-22 Thread Paolo Cavallini
Il 22/02/2012 10:38, Markus Neteler ha scritto: I am afraid that it is a generic R problem. Other example: While GRASS made a PCA in 1-2 min on a set of three radar maps I had to kill R after 20min because I got impatient... OK, I see, thanks. Are R devs aware of this? Should we open tickets?

Re: [GRASS-user] spgrass6

2012-02-22 Thread Markus Neteler
On Wed, Feb 22, 2012 at 11:32 AM, Paolo Cavallini cavall...@faunalia.it wrote: OK, I see, thanks. Are R devs aware of this? I guess so but.. Should we open tickets? ... no idea, perhaps grass-st...@lists.osgeo.org would be the right place to discuss this issue. Markus

Re: [GRASS-user] The algorithm r.cost is not working

2012-02-22 Thread Elisa Solano Villareal
Thanks for your answer. The cost surface used as input for r.cost must have only positive values, negative values will cause infinite loops. Please make sure your input cost surface has only positive values. but we have seen our data and it doesn't have negative values, the minimun value is

Re: [GRASS-user] spgrass6

2012-02-22 Thread Benjamin Ducke
As opposed to GRASS, R has not been designed with computational and/or memory efficiency as a priority. It has its limits when dealing with large datasets (not only GIS datasets). Maybe your analysis would allow you to run your computations on a representative sample instead of the whole

Re: [GRASS-user] spgrass6

2012-02-22 Thread Paolo Cavallini
Il 22/02/2012 12:21, Benjamin Ducke ha scritto: As opposed to GRASS, R has not been designed with computational and/or memory efficiency as a priority. Oh! I was not fully aware of this. Maybe your analysis would allow you to run your computations on a representative sample instead of the

Re: [GRASS-user] spgrass6

2012-02-22 Thread Benjamin Ducke
Have you tried a plotting function that is more efficient for gridded R data? Perhaps plot() from the raster package? I think this is rather a question for the R mailing list, though. Ben On Wed, Feb 22, 2012, at 13:02, Paolo Cavallini wrote: Il 22/02/2012 12:21, Benjamin Ducke ha scritto:

[GRASS-user] r.mask broken in 6.4.2?

2012-02-22 Thread Huub Munstege
Hello fellow 'Grass-rooters' With 6.4.2RC3 and also the recently released 6.4.2 (stable version), I run into this problem after trying to create a mask (either through gui or terminal = r.mask): /opt/grass-6.4.2/scripts/r.mask: eval: line 87: unexpected EOF while looking for matching `''

Re: [GRASS-user] r.mask broken in 6.4.2?

2012-02-22 Thread Micha Silver
On 02/22/2012 03:34 PM, Huub Munstege wrote: Hello fellow 'Grass-rooters' With 6.4.2RC3 and also the recently released 6.4.2 (stable version), I run into this problem after trying to create a mask (either through

[GRASS-user] Re: spgrass6

2012-02-22 Thread Duccio Rocchini
Rely on hexbin() instead of plot() hb - hexbin(your_independent_var, your_dependent_var, xbins = 40) plot(hb) Duccio ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Re: spgrass6

2012-02-22 Thread Roger Bivand
Just 2k may mean 4 million rectangles. R display is vector, hard-copy, with some recent support for raster grids when the rectangles are in fact square. As has been said, the graphics engine is not designed for fast screen output, but for scientific statistical graphics. spplot uses lattice

Re: [GRASS-user] Problem with i.cluster

2012-02-22 Thread Markus Neteler
On Wed, Mar 30, 2011 at 6:28 PM, ndekoning nilsdekon...@hotmail.com wrote: Hello, I have been trying to do an unsupervised classification of my WorldView2 satellite image. When I am using i.cluster I receive the following message: ERROR: Illegal value(s) of sample intervals (100,100). These

Re: [GRASS-user] Re: spgrass6

2012-02-22 Thread Paolo Cavallini
Il 22/02/2012 15:54, Roger Bivand ha scritto: Using the improved raster graphics handling for square cells with image() rather than spplot() and useRaster=TRUE - equivalent to image.SpatialGridDataFrame() and useRasterImage=TRUE with the same matrix takes 1.2 seconds on x11/cairo. You didn't say

[GRASS-user] Re: spgrass6

2012-02-22 Thread Duccio Rocchini
Paolo: at a first glance I was persuaded you would perform exploratory analysis between two sets. The same concepts expressed by Roger Bivand about mapping lattice matrices apply to scatterplots. In case you would perform exploratory analysis on the data I suggest hexagon binning to make quick

Re: [GRASS-user] r.mask broken in 6.4.2?

2012-02-22 Thread Hamish
Huub wrote: With 6.4.2RC3 and also the recently released 6.4.2 (stable version), I run into this problem after trying to create a mask (either through gui or terminal = r.mask): /opt/grass-6.4.2/scripts/r.mask: eval: line 87: unexpected EOF while looking for matching `''

Re: [GRASS-user] The algorithm r.cost is not working

2012-02-22 Thread Markus Metz
On Wed, Feb 22, 2012 at 11:31 AM, Elisa Solano Villareal villar...@fbk.eu wrote: Thanks for your answer. The cost surface used as input for r.cost must have only positive values, negative values will cause infinite loops. Please make sure your input cost surface has only positive values.