Re: [R-sig-Geo] matching points by their coordinates

2009-01-19 Thread Roger Bivand
On Mon, 19 Jan 2009, Struve, Juliane wrote: Dear list members, I would like to look for attributes (e.g. distance from the shore) of spatial points in other data sets (e.g. dist_shoreDF) after matching the points in fishlocationDF and dist_shoreDF by their coordinates. The most direct way

Re: [R-sig-Geo] colorkey print order is incorrect??

2009-01-19 Thread Edzer Pebesma
Jim, your percentage variable is read as a factor; just check: > class(tx3_sp$CK_PCT) this is because R does not understand % symbols (unless you tell it). Convert it to numeric by > tx3_sp$CK_PCTN = as.numeric(unlist(strsplit(tx3_sp$CK_PCT,"%"))) and use spplot on the new variable. The as.

Re: [R-sig-Geo] colorkey print order is incorrect??

2009-01-19 Thread Roger Bivand
On Mon, 19 Jan 2009, Jim Burke wrote: Question. When I use the code below to print a list of values to the right of my graph, the colorkey list appears unusual. At the top the list of percentages tops out at 90.87% instead of the expected 100%. The 100% appears at the bottom below the lowest pe

Re: [R-sig-Geo] OCK in gstat and co-located [SEC=UNCLASSIFIED]

2009-01-19 Thread Edzer Pebesma
Hi Jin, you're not the first to run into this one; read the fit.lmc documentation. fit.lmc fits partial sill matrices that may have perfect correlation among variables. To avoid so, you could instead of jin...@ga.gov.au wrote: ... mud.ock.fit=fit.lmc(x, mud.ock) try: mud.ock.fit = fit.lmc

Re: [R-sig-Geo] OCK in gstat and co-located [SEC=UNCLASSIFIED]

2009-01-19 Thread Jin.Li
Dear all, One more error: > mud.pred<-predict(mud.ock.fit, newdata=data.file.pred) non-positive definite coefficient matrix in structure 2Error in predict.gstat(mud.ock.fit, newdata = data.file.pred) : gstat: value not allowed for: variograms do not satisfy a legal model Thanks for any suggest

[R-sig-Geo] OCK in gstat and co-located [SEC=UNCLASSIFIED]

2009-01-19 Thread Jin.Li
Dear there, I am running OCK for a number of datasets. The codes worked for the first five datasets and then I got an error message: Linear Model of Coregionalization found. Good. [using ordinary cokriging] Linear Model of Coregionalization found. Good. [using ordinary cokriging] Linear Model of

Re: [R-sig-Geo] contour lines on a plot

2009-01-19 Thread Virgilio Gomez Rubio
Hi, El lun, 19-01-2009 a las 22:19 +0100, José Manuel Blanco Moreno escribió: > Just a doubt... there seems that the contours are not shifted but > "mirrored" (flipped) along an horizontal line? Very well spotted!! The right code should be as follows (the change is in the definition of 'spvolcan

[R-sig-Geo] matching points by their coordinates

2009-01-19 Thread Struve, Juliane
Dear list members, I would like to look for attributes (e.g. distance from the shore) of spatial points in other data sets (e.g. dist_shoreDF) after matching the points in fishlocationDF and dist_shoreDF by their coordinates. fishlocationDF$dist=dist_shoreDF$dist[match(dist_shoreDF$coordinat

[R-sig-Geo] conversion from tessellation in spatstat package to sp classes ???

2009-01-19 Thread Volkan Kepoglu
Is there any way to convert from tessellation to sp classes like SpatialLines or SpatialPolygons? library(spatstat) library(maptools) x <- runifpoint(42) t <- dirichlet(x) plot(t) plot(x, add=TRUE) sp <- as(t, "SpatialPolygons") Error in .classEnv(thisClass) : unable to find an environment cont

Re: [R-sig-Geo] contour lines on a plot

2009-01-19 Thread José Manuel Blanco Moreno
Just a doubt... there seems that the contours are not shifted but "mirrored" (flipped) along an horizontal line? Virgilio Gomez Rubio escribió: Hi, El lun, 19-01-2009 a las 18:37 +0100, Maxime Pauwels escribió: Dear Virgilio, I think you're right. Unfortunately, I am not familiar enough

Re: [R-sig-Geo] contour lines on a plot

2009-01-19 Thread Virgilio Gomez Rubio
Hi, El lun, 19-01-2009 a las 18:37 +0100, Maxime Pauwels escribió: > Dear Virgilio, > > I think you're right. Unfortunately, I am not familiar enough with the > analysis of spatial data to avoid such conflict. > Do you know wath function I should use to generate contour before > conversion. Th

Re: [R-sig-Geo] contour lines on a plot

2009-01-19 Thread Maxime Pauwels
Dear Virgilio, I think you're right. Unfortunately, I am not familiar enough with the analysis of spatial data to avoid such conflict. Do you know wath function I should use to generate contour before conversion. I tried contour but it didn't work: contour<-ContourLines2SLDF(contour(p, nlevel

Re: [R-sig-Geo] contour lines on a plot

2009-01-19 Thread Virgilio Gomez Rubio
Dear Maxime, I have not been able to reproduce your example because I lack your data file (file.txt). However, I believe that the problem is at > spplot(p,zcol="K1.pred",col.regions=gray(0:100/100), cuts=40, > sp.layout=list(pts), pretty=T) > contour(p, add=T, nlevels=2) You are combining latti

[R-sig-Geo] contour lines on a plot

2009-01-19 Thread Maxime Pauwels
Dear users, First of all, I'd like to wish you a happy new year, with many nice results from R. I have a graphical problem when plotting kriging results. I first want to see results on a plot and, secondary, I'd like to add contour lines to the plot. The problem is that lines are at a wrong pl