I have a set of points that form a polygon, except that they are in
the wrong order. For example, source into R the object "tmpsub" given
below. Then do
plot(tmpsub)
polygon(tmpsub)
You will see that although the points define a polygon, they are in
the wrong order.
I need an algorithm to p
Freddy,
You can also try the getpoly() function from the splancs package as a more
flexible way to do it.
library(sp)
library(splancs)
plot(meuse$x,meuse$y,type='p',pch='+')
bound <- getpoly()
bound <- rbind(bound,bound[1,]) # Closes the polygon coordinates
bound.bis.poly <- Polygon(bound.bis) #
Try the chull function to create a border.
On Thu, Oct 1, 2009 at 12:39 AM, Freddy López wrote:
> Hello folks,
>
> I'm working with geostatistics data and I'm using mainly geoR library.
> Now, in some examples, is useful to use the 'borders' to delimitate
> our kriging process, e.g., in the krig
Hi Freddy,
from where do you want to extract the information to which spatial
extent you want to interpolate? A shapefile with field borders? You
could also digitalize any polygon with locator(). Anyhow the borders
are not a priori in your data.
regards,
Tom
Am 01.10.2009 um 20:20 schrie
Dear list,
I just was wondering about sofware to perform 4D interpolations.
I've never tried and I'm probably not going to try it in the short time or but
I was asked to give (informal) advice about it and also I was curious.
I have heard GRASS gis can do that and also I've heard about EONFUSION
Hi Freddy,
to add borders to an object of class geodata just add them as an
element to the list (class geodata basically is a list):
obj$borders <- border.polygon
border.polygon here must match coordinates as described in ?polygon,
if I remember correctly.
regards,
Tom
_