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

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

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

i'm really sorry for that. let me state the problem once again with the codes. graphics.off() library(ncdf) ncdata-open.ncdf(air.2m.1950-1999.nc) lon-get.var.ncdf(ncdata,lon) lat-get.var.ncdf(ncdata,lat) time-get.var.ncdf(ncdata,time) dat-get.var.ncdf(ncdata,air) Error: cannot allocate

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

On Mon, 6 Apr 2009, Don MacQueen wrote: 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. I never

Re: [R-sig-Geo] Spatial Interpolation of Regularly Gridded Data

Hi Greg, Yes, there are many possibilities for downscaling grids in R, so you are at the right place. :) 1. If you only wish to downscale climatic grids (e.g. using splines), then probably the most efficient (fastest; can handle large grids) way is to use the downscaling method in SAGA:

[R-sig-Geo] convert spatstat owin-object to splancs

Dear list, I have the following problem: my dataset contains some events (points) and a studyarea, I used some spatstat functions and now I want to use the bandwidth estimation mse2d() in the splancs package. My problem is to convert the spatstat owin-object (studyarea) to use it in the

[R-sig-Geo] Europe zip code shapefiles

Hi all, I'm looking for zip code shapefiles for Italy and France. Has anyone come across these by any chance? thank you! [[alternative HTML version deleted]] ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch

Re: [R-sig-Geo] problem to module MERGE GRID in RSAGA

-Original Message- From: r-sig-geo-boun...@stat.math.ethz.ch [mailto:r-sig-geo-boun...@stat.math.ethz.ch] On Behalf Of gianni lavaredo Sent: Tuesday, April 07, 2009 4:44 PM To: r-sig-geo@stat.math.ethz.ch Subject: [R-sig-Geo] problem to module MERGE GRID in RSAGA Dear R Users

[R-sig-Geo] suggestion for projection

Dear list, I am in search of a good projection with which I could run habitat utilization models on a continental scale. The projection I am looking for should preserve distances between points as good as possible and the map units should be meters. I think from what I have explored so far

Re: [R-sig-Geo] Spatial Interpolation of Regularly Gridded Data

Thanks for your reply Roger. I little while back I got kriging working in the gstat package with my data, but found that for my purposes it was too computationally expensive (read slow). Therefore I am looking for something a little quicker. I now appreciate I cannot use interp due to it

Re: [R-sig-Geo] suggestion for projection

I'm no expert here but I poked around spatialrefernce.org a bit and found some stuff that might get you going in the right direction. North American Equidistant Conic http://spatialreference.org/ref/esri/102010/ South American Equidistant Conic http://spatialreference.org/ref/esri/102032/ Any

Re: [R-sig-Geo] Spatial Interpolation of Regularly Gridded Data

Dear Greg, A while back I did a similar task with the function interp.surface() from the fileds package. It briefly does bilinear interpolation from a regular grid to arbitrary point. It is computationally quite inexpensive since the method of interpolation is quite simple. Hope I