Re: [R-sig-Geo] Help with GIS

2010-10-28 Thread Francisco Javier Santos Alamillos
Thank you for your help! I have finished my task. Best regards. 2010/10/28 Robert J. Hijmans r.hijm...@gmail.com Hi Mike, Thanks for that suggestion, It had not occurred to me to overload as.matrix. Done that now (raster version 1.6-7) such that as.matrix(projras) works. Best, Robert

[R-sig-Geo] !SPAM: prediction error variance

2010-10-28 Thread Pinar Aslantas Bostan
Dear all, I performed gwr() and obtained predictions for 31223 data points. Now I want to calculate prediction error variance at each estimation point. Is there a function to calculate them or do you have any idea about calculating these values? Thanks! Pinar

Re: [R-sig-Geo] !SPAM: prediction error variance

2010-10-28 Thread Roger Bivand
On Thu, 28 Oct 2010, Pinar Aslantas Bostan wrote: Dear all, I performed gwr() and obtained predictions for 31223 data points. Now I want to calculate prediction error variance at each estimation point. Is there a function to calculate them or do you have any idea about calculating these

[R-sig-Geo] Problem nblag_cumul

2010-10-28 Thread Hauser, Christoph
Hi all! As a newcomer to both r and spatial econometrics I am trying to become familiar with both worlds. I have a question which is hopefully not too trivial: I am trying to create a second order queen matrix for European regions. As indicated in the documentation on spdep nblag and

Re: [R-sig-Geo] Why I can't install RGDAL after R 2.12.0 release?

2010-10-28 Thread steven mosher
Thanks for the heads up I was going to switch to 2.12. Guess I will wait for resolution to your issue since I use Kyngchaos 2010/10/27 Ralf Schäfer sena...@ecotoxicology.de Hello, don`t know if this issue is already fixed given that Roger posted an answer. However, I just installed from

[R-sig-Geo] list to points, reproject, interpolate

2010-10-28 Thread Tom Kurkowski
I have a list of points as (X,Y,value) corresponding to longitude, latitude, value. The points come from files that are assigned crs=+proj=longlat +datum=WGS84 My goal is to convert this list into some type of geographically referenced vector format, and then to reproject it to NAD83 Alaska

Re: [R-sig-Geo] Problem nblag_cumul

2010-10-28 Thread Roger Bivand
On Thu, 28 Oct 2010, Hauser, Christoph wrote: Hi all! As a newcomer to both r and spatial econometrics I am trying to become familiar with both worlds. I have a question which is hopefully not too trivial: I am trying to create a second order queen matrix for European regions. As indicated

[R-sig-Geo] fix color scale

2010-10-28 Thread Peter Larson
Hello! I have a problem. I am using IDW to interpolate a daily series of geospatial observations. Thus, I want to produce a large number of sequential maps. I want them to all represent the same color scale. Is there any way to fix the color scale so that it is the same for all the plots?

Re: [R-sig-Geo] list to points, reproject, interpolate

2010-10-28 Thread Robert J. Hijmans
Dear Tom, Here are some hints: Do you really have a list object? Or do you have a matrix or data.frame? If you have a matrix 'xyz', create a data.frame like: xyv - data.frame(xyv) library(sp) #define crs crs - CRS(+proj=longlat +datum=WGS84) # create sp object pts -

Re: [R-sig-Geo] Why I can't install RGDAL after R 2.12.0 release?

2010-10-28 Thread Ralf Schäfer
Sorry, if my post lead to misinterpretation - I basically only wanted to give the code that allows for installing rgdal under 2.12 on Mac OS X. Somebody else posted this question a couple of days ago and the code below resolved the issue for me. So you don`t need to wait necessarily unless you

Re: [R-sig-Geo] Why I can't install RGDAL after R 2.12.0 release?

2010-10-28 Thread Peter Larson
Hello all, I have a grid of points and a spatial polygon shape file. I want to clip the grid so that I only have the point which lie within the shapefile. I tried the following code: ## create a grid onto which we will interpolate: ## first get the range in data x.range - c(37,50) y.range -

Re: [R-sig-Geo] Clipping problem

2010-10-28 Thread Kamran Safi
Hi Pete, you get the indices of the points inside the polygons and the associated polygon ID value. something like clip- grd[!is.na(overlay(grd, ccShapet)),] should do. Hope this helps, although I did not check the line in particular (just from the top of my head) Kamran Am

[R-sig-Geo] Changing the coordinates of a SpatialGridDataFrame

2010-10-28 Thread David.Clifford
Hi there, I have a SpatialGridDataFrame object and I wish to round the coordinates associated with the grid to three decimal places in order to match points in this object with points on another spatial object where the coordinates are recorded to three decimal places. Some information about

Re: [R-sig-Geo] Changing the coordinates of a SpatialGridDataFrame

2010-10-28 Thread Michael Sumner
Probably the best way to do this is to modify the GridTopology and recreate the grid. The coordinates of a SpatialGridDataFrame are calculated from the 6 values of the grid topology, only SpatialPixelsDataFrame stores every cell coordinate (possibly sparsely). It's best to recreate the object,