Re: [R-sig-Geo] 3rd+ rasters added using plot(r, add=TRUE) are misregistered

2013-12-05 Thread GD
Good catch. I have no solution -- but I get the same problem on an earlier version of R in a different environment (so it is not just because of your setup). My session info is below. Cheers, Gareth. sessionInfo() R version 3.0.1 (2013-05-16) Platform: x86_64-unknown-linux-gnu (64-bit)

Re: [R-sig-Geo] covariance matrix

2013-12-05 Thread Benedikt Gräler
Dear Erin, you might want to take a look at the gstat function variogramLine and set the parameter covariance to TRUE. The argument dist_vector takes as well a distance matrix and will then return a covariance matrix. HTH, best wishes, Ben Hodgess, Erin hodge...@uhd.edu wrote: Hi! Here is

Re: [R-sig-Geo] covariance matrix

2013-12-05 Thread Edzer Pebesma
Erin, http://ifgi.uni-muenster.de/~epebe_01/mstp/lec5.html http://ifgi.uni-muenster.de/~epebe_01/mstp/lec7.html contain some minimal, one-dimensional kriging functions that I used in my class recently. R markdown source files are on https://github.com/edzer/mstp On 12/05/2013 04:36 AM,

[R-sig-Geo] Strange lines when projecting a world map

2013-12-05 Thread Francesco Carotenuto
Dear list, as in the title, I found that when projecting a world map by mollweide projection and setting a  longitude value, some strange parallels are drawn. The real problem is that these lines “interact” with a raster when performing the “over” function. How can I manage it? Here I

Re: [R-sig-Geo] original coordinates of the redwood data?

2013-12-05 Thread Barry Rowlingson
On Wed, Dec 4, 2013 at 8:01 PM, Marcelino de la Cruz marcelino.delac...@upm.es wrote: .. The plot of the data published by Strauss (1975) was scanned and digitised by Sandra Pereira, University of Western Australia, 2002. So you can digitize yourself the original true coordinates from the

Re: [R-sig-Geo] covariance matrix

2013-12-05 Thread Hodgess, Erin
Thanks to all for the good help...it's great! Sincerely, Erin From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] on behalf of Edzer Pebesma [edzer.pebe...@uni-muenster.de] Sent: Thursday, December 05, 2013 2:56 AM To:

Re: [R-sig-Geo] Strange lines when projecting a world map

2013-12-05 Thread Roger Bivand
On Thu, 5 Dec 2013, Francesco Carotenuto wrote: Dear list, as in the title, I found that when projecting a world map by mollweide projection and setting a?? longitude value, some strange parallels are drawn. The real problem is that these lines ???interact??? with a raster when performing the

Re: [R-sig-Geo] Strange lines when projecting a world map

2013-12-05 Thread Francesco Carotenuto
Thank you for the answer. The object is the prediction location in an interpolation process. The function of overlay is involved in an automatic process and the longitude (and the kind of projection too) can change according to the dataset that is, in some cases, worldwide distributed.   Dr.

Re: [R-sig-Geo] 3rd+ rasters added using plot(r, add=TRUE) are misregistered

2013-12-05 Thread Michael Sumner
I just wonder if this is related to an existing resize problem when overplotting? I see this on the windows() graphics device, you need to try these two sets of plot commands and resize the window (especially stretch it in one direction) to see the problem at ## 1 and ## 2 v - list(x = seq(0, 1,

Re: [R-sig-Geo] Strange lines when projecting a world map

2013-12-05 Thread Roger Bivand
On Thu, 5 Dec 2013, Francesco Carotenuto wrote: Thank you for the answer. The object is the prediction location in an interpolation process. The function of overlay is involved in an automatic process and the longitude (and the kind of projection too) can change according to the dataset that

Re: [R-sig-Geo] 3rd+ rasters added using plot(r, add=TRUE) are misregistered

2013-12-05 Thread Josh O'Brien
Hi Michael, Thanks for those thoughts. The issue you raise (which has bit me before as well) appears to be unrelated. In particular, your resizing issue disappears if one sets useRaster=FALSE, whereas my problem persists whether or not rasterImage() is used to draw the raster. Try this to get a

Re: [R-sig-Geo] 3rd+ rasters added using plot(r, add=TRUE) are misregistered

2013-12-05 Thread Michael Sumner
Interesting and difficult! Thanks for pointing out the useRaster difference, I had mixed up TRUE/FALSE in my original report (offlist) on that one. That might help me find the solution, though sorry I have nothing for your issue. :) Cheers, Mike. On Fri, Dec 6, 2013 at 11:58 AM, Josh O'Brien

[R-sig-Geo] original coordinates of the redwood data?

2013-12-05 Thread Adrian Baddeley
Lee De Cola ldec...@comcast.net writes: does anyone know the original, true coordinates of the spatstat redwoodfull data? i can't find them in: Strauss, D. J. (1975). A Model for Clustering. Biometrika 62(2): 467-475. i think reporting spatial data in rescaled units is unscientific.