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

2009-04-07 Thread Greg King
Thanks Miha, but it doesn't look like this function supports spatial data frames and great circle distances. On Tue, Apr 7, 2009 at 10:54 PM, Miha Staut wrote: > > Dear Greg, > > A while back I did a similar task with the function interp.surface() from > the fileds package. It briefly does bilin

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

2009-04-07 Thread Miha Staut
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 understood

[R-sig-Geo] georectify an image with tiepoints

2009-04-07 Thread Josh London
Hello I trying to work out method for georectifying QuickSCAT sea-ice GIFs (http://nereids.jpl.nasa.gov/cgi-bin/ice.cgi?show=scat). The daily images are consistent in their layout and include overlayed graticules. In ArcMap, I create tiepoints at the graticule intersections and geo-rectify from t

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

2009-04-07 Thread Alex Mandel
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 cha

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

2009-04-07 Thread Greg King
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 treatin

[R-sig-Geo] suggestion for projection

2009-04-07 Thread Kamran Safi
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 fr

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

2009-04-07 Thread Tomislav Hengl
> -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

[R-sig-Geo] Europe zip code shapefiles

2009-04-07 Thread Alina Sheyman
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 https://stat.ethz.ch/mailman

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

2009-04-07 Thread gianni lavaredo
Dear R Users I have a problem to understand the loop merge in RSAGA I have this file (not order name) DEM2644147 DEM2644148 DEM2644149 DEM2744150 DEM2744155 DEM2744168 I did a file list in CSV: path<- ("C:/mydata/") temp <- paste(path,"input.csv", sep="") input <- read.csv(temp, header = TRUE

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

2009-04-07 Thread Robert Hijmans
Anwesha, Assuming you want the data timestep by timestep (but you can adjust it for all timesteps, pixel by pixel) you can do something like the below. I copied it straight from ?get.var.ncdf library(ncdf) ncdata<-open.ncdf("air.2m.1950-1999.nc") lon<-get.var.ncdf(ncdata,"lon") lat<-get.var.ncd

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

2009-04-07 Thread Tomislav Hengl
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

2009-04-07 Thread Marco Helbich
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 splanc

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

2009-04-07 Thread Roger Bivand
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 u

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

2009-04-07 Thread Anwesha Chakrabarti
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:

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

2009-04-07 Thread Roger Bivand
On Mon, 6 Apr 2009, Greg King wrote: Hi, I'm fairly new to R, but finding the experience a good one. However I am a little overwhelmed by the number of packages and not sure I am necessarily using the most appropriate ones. Have you read the Spatial Task View on your nearest CRAN? There you

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

2009-04-07 Thread Roger Bivand
On Tue, 7 Apr 2009, Anwesha Chakrabarti wrote: 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 ms