Re: [R-sig-Geo] Creating a maximum raster layer from RasterStack

2013-10-16 Thread Agustin Lobo
Robert, I keep on because this is an important point for the analysis of time series of image: Actually I get stackApply() faster than max() if the nb of layers increases: Example with 3 larger layers: r1 - raster(matrix((1:900),3000,3000)) r2 -

Re: [R-sig-Geo] Package GWmodel Mixed GWR estimation

2013-10-16 Thread binbin lu
Hi Roosbeh, Dear members, I am using gwr.mixed and gwr.basic functions in package GWmodel. In the function specifications it has mentioned that argument: fixed.var, which is independent variables appeared in the formula to be treated as global I am wondering how should I specify (maybe

[R-sig-Geo] Long time to load shapefiles

2013-10-16 Thread philippe massicotte
Hi everyone. I'm trying to open some shapefiles, but the process is still not complete after 25 minutes. Looking at the shapefiles in ArcGIS, I find that these files have roughly 500K entries. I was wondering if this processing time to open the shapefiles is normal. Here is my configuration:

Re: [R-sig-Geo] Selecting cells in a raster

2013-10-16 Thread Agustin Lobo
Manuel, What do you mean by it didn't work? You should have got an object named forest, post the same info as for cob_r Agus On Tue, Oct 15, 2013 at 6:49 PM, Manuel Spínola mspinol...@gmail.com wrote: Dear list members, I have a raster with land use data (with values from 1 to 9). I want to

Re: [R-sig-Geo] Selecting cells in a raster

2013-10-16 Thread Manuel Spínola
Thank you Agus, Yes, it works. I did not recognize that the raster has 2 levels, 0 and 1, and the 1s were forest. I put 2 colors and now I can display the forest in the raster. Best, Manuel 2013/10/16 Agustin Lobo alobolis...@gmail.com Manuel, What do you mean by it didn't work? You

Re: [R-sig-Geo] Long time to load shapefiles

2013-10-16 Thread GD
Hi Phil, I have also found readOGR to be very slow when reading a large shapefile (in my case, a shapefile with 85000 polygons and ~ 200 columns in the attribute table). In my case, I would repeatedly re-read the shapefile each time I was working on the script. To speed things up, I found it

Re: [R-sig-Geo] Creating a maximum raster layer from RasterStack

2013-10-16 Thread Robert J. Hijmans
Hi Agus, Yes, you are right, in that case stackApply is faster (but note that you had a mistake in the indices argument that further increased the difference but is another reason why I would not use stackAppy unless I really needed it). I did not think it would be possible that a more