similar to: attempt at making a polygon class failed

Objective: Determine if a set of x and y-latitude points are inside of a polygon using R. Lets say I have 9 polygons. Where I have labeled the polygons to be checked from 1-9. The problem I?m running into is running the point.in.polygon to check if those points are in one of several polygons, as my code overwrites the result. To accomplish this I am using the R-built in function called

I am a new user of the R spatstat package and am having problems creating a polygonal observation window with owin(). Code follows: library("maps") library ("sp")` library("spatstat") mass.map <- map("state", "massachusetts:main", fill=T) # This returns a data frame includding x and y components that form a polygon of massachusetts mainland`

Hello Friends!!! I would want draw a circular histogram, and I would like draw a polygon with a circular side. This is easy if I use the functions polygon and arc, but I want that the polygon with a circular side have background colour. The polygon created with function polygon can have background colour, but the surface created with function arc can?t have background colour. How I could create a

Dear list: I encountered some problems using the function point.in.polygon() of the sp package, when trying to determine whether some points lye inside, outside, on the border or on a vertice of a polygon. I have a list of point I know should lye right on the border of a polygon, but some of them are not classified as such by point.in.polygon() (see the example code below). To make a long story

Tena koe Simon plot(1:10, 1:10, type='n') polygon(c(2,3,6,8), c(2,5,5,3), density=20, angle=90) polygon(c(2,3,6,8), 5+c(2,5,5,3), density=20, angle=0) I don't understand your problem. Perhaps if you "provide[d] commented, minimal, self-contained, reproducible code" it would help. HTH .... Peter Alspach -----Original Message----- From: R-help [mailto:r-help-bounces at

The angle is not based on the polygon edges, but it can seem that way if you do not use the asp=1 argument in your plot. Try this example, > plot(1:10, 1:10, type='n') > x <- c(1, 3, 5, 3) > y <- c(3, 5, 3, 1) > polygon(x, y, angle=0, density=10) > polygon(x, y + 4.5, angle=45, density=10) > polygon(x + 4.5, y + 4.5, angle=90, density=10) > polygon(x + 4.5, y,

Thanks for the question. Here is a sample of the code for my plot: Top = c(34, 39, 42, 45, 46, 41, 41, 40, 43, 38, 33, 33) Bottom = c(24, 29, 32, 36, 32, 34, 32,41, 40, 39, 29, 24) plot(1,1, col = "white", xlim = c(1.3,11.7), ylim = c(0,80), axis = FALSE, xaxt = "n") axis(1, at = c(1:12)) polygon(c(c(1:12),c(12:1)), c(top, bottom), col =

I have a complex 2D polygon with thousands of vertices, and I'd like to be able to identify points from a large set contained within the polygon, and was wondering if there might be an efficient way of doing this? Any advice would be useful! Here is a small example of what I mean: # make polygon v1<-c(0,1,1,2,1,3,6,7) v2<-c(1,3,3,5,6,7,8,9) plot(v1, v2, type = "n" )

Dear R users; Is there any way to intersect a filled contour image and a polygon? My problem is that I want to create a kriging map and the boundaries of my map are given by the coordinates of the polygon. So far I can superompose the polygon in the filled.contour image but I don't know how to get rid of the contour image outside of the polygon boundaries. Any hint will be appreciated PM

Hi, I'm trying to use grid.polygon() to plot several polygons at once, with a view to putting coloured polygons beneath a curve. I'm struggling just to get the grid.polygon to plot anything # PLOT SOME POINTS x <- 1:100 y <- 1:100*0.5 + 3 plot(x, y, pch = ".") # PLOT 2 POLYGONS library(grid) grid.polygon(x = c(20, 20, 40, 40, 40, 40, 60, 60), y =

Dear forum, I would like to use the kernel2d or spkernel2d in the Splancs-package, but it does not recognize my polygon data. "Error in kernel2d(as.points(ptsbin), polygonprov, h0 = 2, nx = 100, : " is the error message. Invalid poly argument The data are defined as follows: polgonprov<-list(x=polyprov$X, y=polyprov$Y) with X and Y coordinates in the Lambert1972

HI, Sorry Carl, I received your message in my spam folder. Sarah proposed me a good example of code. Thank you Momadou  ________________________________ De : Carl Witthoft [via R] <ml-node+s789695n4188375h96@n4.nabble.com> Envoyé le : Mardi 13 Décembre 2011 3h34 Objet : Re: Polygon Please read the posting guide and provide a (small) reproducible example of your data. The statement

hi. below is a small test case (hopefully minimal, though i'm still a bit confused about initializers). i would have guessed (and maybe i still would have been right) that one could re-use the name of a generic function for functions with different numbers of arguments. in the case below, class A's bB() queries the status of a single A object, so bB(A) (where here "A" is an

Filling polygons with lines is a throwback to the time when the height of quality graphics was the mechanical pen plotter (a device that used a pen in a mechanical arm to draw the plot on a piece of paper). Computing and printing technology has advanced quite a bit from that day, so you may want to reconsider why you want polygons filled with lines instead of just a solid color (and I consider

Hi all, I would like to color the area between two time-series. I tried it by using the polygon() function but I keeps drawing lines between beginning and end points. Is there another more appropriate function or how could I close the polygon at the end en the beginning of the time series (e.g., drawing a straight line)? The following doesn't plot a polygon between the two

Preamble: Sorry for being dense. Now that that's done, here are my questions. I want to put a polygon on a plot of a time series. I'm going to add lines from a smooth.spline interpolation and other annotation to it. But here's the general idea: # Start code n <- 121 dat <- rnorm(n) plot(dat, type="n") polygon(c(1:n,n:1), c(dat,c(rep(0,n))), col =

Hello dear members of R-help and R-core mailing list, I am not sure if this request is a "ticket" that should be filled somewhere outside the mailing list. If so, I apologize for not doing and would like to know where I should have filled it. And to the subject matter: I would like to use a command like this: plot(c(1,8), 1:2, type="n") polygon(1:7, c(2,1,2,NA,2,1,2),

Hello dear members of R-help and R-core mailing list, I am not sure if this request is a "ticket" that should be filled somewhere outside the mailing list. If so, I apologize for not doing and would like to know where I should have filled it. And to the subject matter: I would like to use a command like this: plot(c(1,8), 1:2, type="n") polygon(1:7, c(2,1,2,NA,2,1,2),

R-listers, Given a polygon and a circle defined by its center coordinates and a radius, I would like to calculate the area of overlap. I know that I can create a polygon from the circle and then use available packages to get the area of the intersection. However, because the polygon is of a fixed size and I will be doing this for circles of varying sizes, I'm concerned about

Hello, I have the contour of the shoreline of a lake. I have measured data points away from the lake shore signifying a particular depth. The data is not ordered as a polygon and hence trying to draw it creates a series of lines all over the place. I would like to turn the points into a ordered set that results in a smooth polygon. I have tried a "nearest point" approach. An example