Re: [R-sig-Geo] Question about a spatial non-square cells error in R

2010-01-05 Thread Roger Bivand
On Tue, 5 Jan 2010, Mailhiot, Joshua wrote: Hi folks, Long-time lurker, first time poster here. I have a conundrum that perhaps you might have come across. I'm getting an error when I try to compile some of my R spatial code: Error in write.asciigrid(xx2_res["NPHS"], newHSmap_asc, attr = 1, na

Re: [R-sig-Geo] One-to-many table relationship to shapefile polygon in spplot / plot - possible?

2010-01-05 Thread Roger Bivand
On Wed, 6 Jan 2010, Noli Sicad wrote: Hi, I have problem transposing my database table which is to be merged / joined with the shapefile file (*.dbf) inorder to have one-to-one relationship with shapefile polygon. As we know that inorder for the polygon to be rendered in GIS, it has to have one

Re: [R-sig-Geo] Testing if a point-pattern is on a regular grid

2010-01-05 Thread Roger Bivand
On Wed, 6 Jan 2010, Pierre Roudier wrote: Dear list, I'm looking for an easy way to test if a point pattern (I got a data frame of x and y coordinates) is on a regular grid or not. I'm aware of the gridded() function in sp(), but this function is devoted to assign rather than to test, as far a

[R-sig-Geo] One-to-many table relationship to shapefile polygon in spplot / plot - possible?

2010-01-05 Thread Noli Sicad
Hi, I have problem transposing my database table which is to be merged / joined with the shapefile file (*.dbf) inorder to have one-to-one relationship with shapefile polygon. As we know that inorder for the polygon to be rendered in GIS, it has to have one-to-ope relationship. One to many table:

[R-sig-Geo] Testing if a point-pattern is on a regular grid

2010-01-05 Thread Pierre Roudier
Dear list, I'm looking for an easy way to test if a point pattern (I got a data frame of x and y coordinates) is on a regular grid or not. I'm aware of the gridded() function in sp(), but this function is devoted to assign rather than to test, as far as I know. I'm pretty sure there's an easy way

[R-sig-Geo] Question about a spatial non-square cells error in R

2010-01-05 Thread Mailhiot, Joshua
Hi folks, Long-time lurker, first time poster here. I have a conundrum that perhaps you might have come across. I'm getting an error when I try to compile some of my R spatial code: Error in write.asciigrid(xx2_res["NPHS"], newHSmap_asc, attr = 1, na.value = -) : Asciigrid does not support

[R-sig-Geo] Kriging concept question

2010-01-05 Thread Tobin Cara
Hello, I have recently read an interesting article about integrating Limited Area Models (LAMs) into kriging with external drift for temperature (Libertà et al. link below). www.wmo.int/pages/prog/www/IMOP/.../IOM.../P2(05)_Perini_Italy.doc As I understand, it seems that the authors generated

Re: [R-sig-Geo] mixed geographically weighted regression

2010-01-05 Thread Marco Helbich
Dear Roger, thank you for your quick response! If I understand it correctly, the hat matrix is calculated using all explanatory variables. In my case, however, I would need to restrict the column space to those covariates where I assume varying coefficients (as in eq. (3)), and for this purpo

Re: [R-sig-Geo] mixed geographically weighted regression

2010-01-05 Thread Roger Bivand
On Tue, 5 Jan 2010, Marco Helbich wrote: Dear Roger, thank you for your quick response! If I understand it correctly, the hat matrix is calculated using all explanatory variables. In my case, however, I would need to restrict the column space to those covariates where I assume varying coeffi

Re: [R-sig-Geo] mixed geographically weighted regression

2010-01-05 Thread Roger Bivand
On Tue, 5 Jan 2010, Marco Helbich wrote: Dear list, I am trying to fit a mixed geographically weighted regression model (with adaptive kernel) using the spgwr package, i.e. I want to hold some of the coefficients fixed at the global level. Thus, I have the following questions: 1. Which is th

[R-sig-Geo] mixed geographically weighted regression

2010-01-05 Thread Marco Helbich
Dear list, I am trying to fit a mixed geographically weighted regression model (with adaptive kernel) using the spgwr package, i.e. I want to hold some of the coefficients fixed at the global level. Thus, I have the following questions: 1. Which is the most efficient way to estimate such a mode

Re: [R-sig-Geo] prism maps

2010-01-05 Thread Hisaji ONO
Hello. How about this. http://translate.google.co.jp/translate?js=y&prev=_t&hl=ja&ie=UTF-8&layout=1&eotf=1&u=http%3A%2F%2Fwww.okada.jp.org%2FRWiki%2F%3FGoogleEarth%25A4%25C8R%23d46aa784&sl=ja&tl=en Original language is Japanese. Regards. --- Thomas Adams

Re: [R-sig-Geo] Spatial correlogram - Problems with creating neigbourhood list and weights

2010-01-05 Thread Roger Bivand
On Tue, 5 Jan 2010, Strubbe Diederik wrote: Dear all, I have a dataset with several lat-long coordinates (48) and corresponding residuals values (not a grid). Using the spdep package, I try to create a correlogram of these residuals according to several lags (i.e. distance classes) but I can

Re: [R-sig-Geo] raster / polygon overlay

2010-01-05 Thread laure velez
Thanks again for the fast reply, In fact this problem is only a tiny part of the project and I'll need to run this overlay operation over more than 600 polygons (watersheds) to get mean values (and standard error) of ~10 environmental variables (rasters that are not at the same resolution,

[R-sig-Geo] Spatial correlogram - Problems with creating neigbourhood list and weights

2010-01-05 Thread Strubbe Diederik
Dear all, I have a dataset with several lat-long coordinates (48) and corresponding residuals values (not a grid). Using the spdep package, I try to create a correlogram of these residuals according to several lags (i.e. distance classes) but I can’t seem to get the neighborhood and weights rig

Re: [R-sig-Geo] prism maps

2010-01-05 Thread Thomas Adams
Adam, Do you mean PRISM maps such as these: http://www.prism.oregonstate.edu/ ? If yes, then the question would be, do you want to duplicate the PRISM methodology or effectively get similar results? One could use universal kriging using gstat or methods proposed by Tom Adam A wrote: Hello

Re: [R-sig-Geo] prism maps

2010-01-05 Thread Roger Bivand
On Tue, 5 Jan 2010, Adam A wrote: Hello! Can I make prism maps in R library(fortunes) fortune("Yoda") But should you want to - they do seem to be much more open to misinterpretation than more conservative representations? Roger

Re: [R-sig-Geo] raster / polygon overlay

2010-01-05 Thread Roger Bivand
On Tue, 5 Jan 2010, laure velez wrote: Many thanks for the fast replies, @ Roger : polygons and raster are ok (see attached file). @ Johannes : I have already try this but the result was not satisfying, using a SpatialPointsDataFrame do not solve the problem as the points are not inside

Re: [R-sig-Geo] raster / polygon overlay

2010-01-05 Thread laure velez
Many thanks for the fast replies, @ Roger : polygons and raster are ok (see attached file). @ Johannes : I have already try this but the result was not satisfying, using a SpatialPointsDataFrame do not solve the problem as the points are not inside the polygon (see attached file). I gues

[R-sig-Geo] prism maps

2010-01-05 Thread Adam A
Hello! Can I make prism maps in R Kibicuj biało-czerwonym w Zakopanem Największa impreza sportowa tej zimy! Zobacz zdjęcia: http://klik.wp.pl/?adr=http%3A%2F%2Fcorto.www.wp.pl%2Fas%2Fskoki_galeria.html&sid=942

Re: [R-sig-Geo] (no subject)

2010-01-05 Thread Roger Bivand
On Tue, 5 Jan 2010, laure velez wrote: Dear list member, Although the answer to this might be straightforward I'm struggling with this issue for weeks now : I have a shapefile (only one polygon) and a raster (.tif). I would like to know the value of each pixel on the raster which overlap wi

[R-sig-Geo] (no subject)

2010-01-05 Thread laure velez
Dear list member, Although the answer to this might be straightforward I'm struggling with this issue for weeks now : I have a shapefile (only one polygon) and a raster (.tif). I would like to know the value of each pixel on the raster which overlap with the polygon in the shapefile. I ha