[R-sig-Geo] Dilation and Window Definition in Spatstat

2013-11-02 Thread Adrian Baddeley
On 11/01/13 04:40, Ted Rosenbaum wrote: > I am using spatstat version 1.34-0 and the pop command is not correctly > recognizing the interior of a window after I use the dilation command on a > polygon (i.e., non rectanglular window). This was a bug in 'is.subset.owin' which affected the functions

Re: [R-sig-Geo] Dilation and Window Definition in Spatstat

2013-10-31 Thread Adrian Baddeley
In spatstat 1.34-0 there appears to be a bug in union.owin which also affects dilation.owin. A fix will be distributed shortly. Until this is fixed, the fastest way to do what you want is to convert the window to a pixel mask (using as.mask) before applying the dilation. Adrian Baddeley Pro

Re: [R-sig-Geo] Dilation and Window Definition in Spatstat

2013-10-31 Thread Rolf Turner
It is hard to tell without reproducible code, but I suspect (since you use the word "dilate") that you have used the function dilation() and are not getting the window that you think you are getting. Try: W <- owin(poly=as.owin(demopat)$bdry[[1]]) WD <- dilation(W,100) set.seed(4

[R-sig-Geo] Dilation and Window Definition in Spatstat

2013-10-31 Thread Ted Rosenbaum
I am using spatstat version 1.34-0 and the pop command is not correctly recognizing the interior of a window after I use the dilation command on a polygon (i.e., non rectanglular window). I have window that is a convex hull of the lower 48 states. I have a set of points that are county centroids.