[R-sig-Geo] Put polygon points in correct order

2009-10-01 Thread Don MacQueen
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

Re: [R-sig-Geo] Creating borders.-

2009-10-01 Thread Pierre Roudier
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) #

Re: [R-sig-Geo] Creating borders.-

2009-10-01 Thread Henrique Dallazuanna
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

Re: [R-sig-Geo] Creating borders.-

2009-10-01 Thread Tom Gottfried
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

[R-sig-Geo] curious about 4D interpolations

2009-10-01 Thread Pilar Tugores Ferra
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

Re: [R-sig-Geo] Creating borders.-

2009-10-01 Thread Tom Gottfried
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 _