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

2009-04-23 Thread Clint Bowman
Greg, I'm quite certain that kriging would not be appropriate with model output--you won't have any errors to produce the necessary statistics. Are you looking at smooth terrain (e.g., over the ocean) where topographic influences are negligible and you may want to interpolate the pressure field a

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

2009-04-23 Thread Dylan Beaudette
On Thursday 23 April 2009, Greg King wrote: > This reply is a little late, but after much searching, I cannot find an R > package/function that fits my needs. > > I am trying to perform either bilinear or bicubic interpolation using > great circle distances (ie using a spherical model of the earth

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

2009-04-23 Thread Greg King
This reply is a little late, but after much searching, I cannot find an R package/function that fits my needs. I am trying to perform either bilinear or bicubic interpolation using great circle distances (ie using a spherical model of the earth). Here is what I have tried: gstat: Can only perfo

Re: [R-sig-Geo] creating polygons

2009-04-23 Thread Roger Bivand
On Thu, 23 Apr 2009, Alina Sheyman wrote: I went back to the help page, but it still didn't clarify it for me. Yes, the two files are different length, but both have a zip variable that they can match on and then keep only those records that match. Or do I need to additionally create an index fo

Re: [R-sig-Geo] creating polygons

2009-04-23 Thread Alina Sheyman
I went back to the help page, but it still didn't clarify it for me. Yes, the two files are different length, but both have a zip variable that they can match on and then keep only those records that match. Or do I need to additionally create an index for the match to occur? I am understanding it t

Re: [R-sig-Geo] paste(" /" ") and paste(" /' ") for grass command

2009-04-23 Thread Greg Snow
It is not the paste function (as has been mentioned before), but the print function that is doing different things. The default behavior for printing character strings is to put "" around each string, if there is a " inside the string, then the print function will put the \ in the representatio

Re: [R-sig-Geo] Problems with Raster package?

2009-04-23 Thread Robert Hijmans
Hi Ned, that is the downside of using a package under development (albeit *quite* stable now). "x <- setRaster(x)" has been replaced by the (generic) "x <- raster(x)" "x <- setDatatype(x, value)" has been replaced by "dataType(x) <- value" As for the versions, you are not doing anything wrong. It

Re: [R-sig-Geo] spdep: stsls; tot.solve=1e-12

2009-04-23 Thread Roger Bivand
On Wed, 22 Apr 2009, Adam Okulicz-Kozaryn wrote: Hello, i am trying to run stsls from spdep, and I am getting error: Error in solve.default(ZpZp) : system is computationally singular: reciprocal condition number = 0 I would get the same error in lagsarlm, but I can fix it by setting I don'

Re: [R-sig-Geo] Find incircle of a polygon

2009-04-23 Thread Roger Bivand
On Wed, 22 Apr 2009, Etienne Bellemare Racine wrote: Thanks Roger it helped me a lot. Now I would like to iteratively find the incircle in a polygon until a certain radii (i.e. as if I was stuffing balls in my polygon, from the biggest to the smallest). So I have to subtract that circle to my

[R-sig-Geo] Help to OVERLAY (in sp) a "SpatialPointsDataFrame" with a 'SpatialPolygonsDataFrame"

2009-04-23 Thread gianni lavaredo
Dear R User , Sorry for this email, but i have a problem with overlay in sp. I have to shape file one is a points shape file (clip_point, with "pH" value) the second is a Polygons shapefile (plots, with ID number of plot) clip_point <- readOGR(".", "clip_point") proj4string(clip_point) <- CRS("+

[R-sig-Geo] Problems with Raster package?

2009-04-23 Thread Ned Horning
I just updated my R installation to 2.9 and then updated my packages. It looks like setRaster and setDatatype are not included in the Raster package that I downloaded (8.9-13). I noticed on the R-Forge site for the Raster package that version 8.9-14 is the most current version but when I run "i

Re: [R-sig-Geo] spplot [SEC=UNCLASSIFIED]

2009-04-23 Thread Edzer Pebesma
jin...@ga.gov.au wrote: > Thank you, Edzer. > My second question bas been solved by using: > proj4string(pred.mud.ked) <- CRS("+proj=longlat +datum=WGS84") > > For the first question, I have checked ?levelplot, but found no clue to do > what I intended. Thanks for any further help. > library(

Re: [R-sig-Geo] error using predict.sarlm

2009-04-23 Thread Roger Bivand
On Wed, 22 Apr 2009, evans...@umn.edu wrote: Hi Roger, Thanks for the response. Unfortunately, I'm still getting the same error running predict.sarlm ("Error in if (nrow(newdata) != length(listw$neighbours)) stop("mismatch between newdata and spatial weights") : argument is of length zero").

Re: [R-sig-Geo] spplot [SEC=UNCLASSIFIED]

2009-04-23 Thread Jin.Li
Thank you, Edzer. My second question bas been solved by using: proj4string(pred.mud.ked) <- CRS("+proj=longlat +datum=WGS84") For the first question, I have checked ?levelplot, but found no clue to do what I intended. Thanks for any further help. Cheers, Jin -Original Message- From: Edz