Generating a polygon around points
On 16/08/14 02:51, Bob O'Hara wrote:
I've been struggling for half a day on what should be a simple problem... I have a data frame of lat/long coordinates that describe a region, and I want to draw a polygon around them so I can use that as a boundary (to be thrown at INLA, but those details aren't important). The coordinates are almost on a regular grid: there is some variation in latitude (because we're on a globe). If the coordinates were on a regular grid, I could use as.owin() to create a mask, and go on from there (I have code that will work). But as.owin() doesn't like unevenly spaced points. Can anyone suggest a way to sort this out? Preferable without having to mess around transforming the coordinates.
You need to be able to tell owin() the vertices of the polygon that you want to use as your window. It is possible that alphahull, as someone suggested, might give what you want; I don't know, I am not familiar with this function. OTOH R is still lacking a mind_read() function so it probably would NOT give you *exactly* what you want. If you want to get "exactly the right polygon" you probably will have to select out the vertices of the polygon (in the appropriate order) in some "by hand" manner. It's possible that clickpoly() might be of help to you. cheers, Rolf Turner
Rolf Turner Technical Editor ANZJS