[R-sig-Geo] how to handle large datafile in R?

2009-04-06 Thread Anwesha Chakrabarti
hi all, i want to read a .netcdf file in R. the file contains one variable that has three dimensions lon, lat and time. and the data consists of a 192x94 matrix of the variable. whenever i'm trying to open the file i'm getting the following error msg "can not allocate vector of size 1.2gb"! so i gu

[R-sig-Geo] cross 2 SGDF with different extents

2009-04-06 Thread milton ruser
Dear all, I have two SpatialGridDataFrame objects with different extent: > str(reser) Formal class 'SpatialGridDataFrame' [package "sp"] with 6 slots ..@ data :'data.frame': 1826316 obs. of 1 variable: .. ..$ band1: int [1:1826316] 0 0 0 0 0 0 0 0 0 0 ... ..@ grid

Re: [R-sig-Geo] R interface to SEXTANTE?

2009-04-06 Thread Alexander Brenning
Hi Nick, SEXTANTE and SAGA are very closely related, many of their modules are identical. So at least the gain in functionality that can be achieved by interfacing SEXTANTE would probably be small, but maybe you know more about this. A SEXTANTE-R interface would be a bit different from SAGA-R

Re: [R-sig-Geo] raster package and "?focal" FUN

2009-04-06 Thread Alexander Brenning
Hi, I haven't used the focal function in raster, which should be able to do this, but in the RSAGA package there is a function called focal.function that does the same thing on ASCII grids. In your case this should work: focal.function("dem",fun="max",radius=3,mw.to.vector=TRUE,mw.na.rm=TRUE)

[R-sig-Geo] raster package and "?focal" FUN

2009-04-06 Thread gianni lavaredo
I am seeing RASTER package and "focal". I have a raster ascii grid format with a Z value. I am trying to do a mobile window with a specific radius (es: 10 m) to find the pixel with the max value of Z inside the window. I need to do this operation for each pixel and save a a new grid with only the

[R-sig-Geo] Plotting depends on search() order?

2009-04-06 Thread Don MacQueen
I am finding that plotting a SpatialPolygonsDataFrame object fails when the base graphics packages are ahead of sp in the search() path. This occurs when I put require(maptools) [ or require(rgdal), or require(sp) ] in my .Rprofile file. In the examples below, the object 'srdf' was create