search for: polygon

I have an XML file that has within it the coordinates of some polygons that I would like to extract and use in R. The polygons are nested rather deeply. For example, I found by trial and error that I can extract the coordinates of one of them using functions from the XML package: doc <- xmlInternalTreeParse('doc.kml') docroot <- xmlRoot(doc)...

Hola He hecho el típico mapa con gadm y spplot: library(sp) ES <- local(get(load(url("http://biogeo.ucdavis.edu/data/gadm2/R/ESP_adm2.RData ")))) acerco las canarias: A <- which(ES en data$NAME_1 == "Islas Canarias") L <- length(ES en polygons[[A]]@Polygons) L <- length(ES en polygons[[A]]@Polygons) for (i in 1:L){ ES en polygons[[A]]@Polygons[[i]]@coords <- cbind(ES en polygons [[A]]@Polygons[[i]]@coords[,1]+6,ES en polygons [[A]]@Polygons[[i]]@coords[,2]+5) } El tema es q quiero añadirle un recuadro a las islas usando spplot...

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 funct...

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 massachus...

I am trying to use the angle option in polygon to create polygons filled with horizontal and vertical lines. The polygons I am crating are irregular and it the angle function appears to set the angle of the shading perpendicular to the polygon sides rather than perpendicular to the axes. Is there any way to set the angle relative to the axes...

Dear all, is there any function to calculate the center of a polygon mass in R? Actually I need to find the best location within polygons to place labels. Thanks for any hint Jens Oehlschl?gel -- COMPUTERBILD 15/03: Premium-e-mail-Dienste im Test\ --------...{{dropped}}

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 cou...

Dear all, I am currently working with the spatstat package, using windows and pixel images. First: My aim is to transform a shapefile (see attached) into a pixel image. My idea is to start transforming the shapefile into a Spatial Polygon file: x <- readShapeSpatial("200001441.shp") y <- as(x, "SpatialPolygons") z <- as.owin(y) Given z, I want to identify each polygon with a single constant value. This is like adding marks to the SpatialPolygons file. Then I want to convert these polygons into an image...

...pico mapa con gadm y spplot: > > > library(sp) > ES <- > local(get(load(url("http://biogeo.ucdavis.edu/data/gadm2/R/ESP_adm2.RData > ")))) > > acerco las canarias: > > A <- which(ES en data$NAME_1 == "Islas Canarias") > L <- length(ES en polygons[[A]]@Polygons) > > L <- length(ES en polygons[[A]]@Polygons) > for (i in 1:L){ > ES en polygons[[A]]@Polygons[[i]]@coords <- cbind(ES en polygons > [[A]]@Polygons[[i]]@coords[,1]+6,ES en polygons > [[A]]@Polygons[[i]]@coords[,2]+5) > } > > El tema es q quiero aña...

Background: I have an algorithm which produces a large number of small polygons (of the spatial kind) which I would like to use within R using objects from sp. I can't predict the exact number of polygons a-priori, the polygons will be grouped into regions, and each region will be filled sequentially, so an appropriate C++ 'framework' (for the point of illustratio...

...e 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), col=c("red", "blue"), # border=c("green", "yellow"), border=c(1,10), lwd=c(1:10)) To create two triangles, with different line widths. But the polygon command doesn't seem to recycle the &q...

...e 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), col=c("red", "blue"), # border=c("green", "yellow"), border=c(1,10), lwd=c(1:10)) To create two triangles, with different line widths. But the polygon command doesn't seem to recycle the &q...

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 +...

Hi there, I want to integrate a function over an irregular polygon. Is there any function which can implement this easily? Otherwise, I am thinking of divide the polygon into very small rectangles and use "adapt" to approximate it. Do you have any suggestions to get the fine division? Any advice is appreciated. Haiyong

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 belo...

R2.7.1, WinXP Hi, I have a polygon inside a circle as follows: radius <- 3 x <- seq(-radius,radius,length=2000) y <- sqrt(radius^2-x^2) xx <- c(x,-x) yy <- c(y,-y) plot(xx,yy, xlim=c(-radius,radius),ylim=c(-radius,radius), type="l", ylab="", xlab="", axes=F) radius <- 2.7 x1 <- s...

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-----...

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...

Dear all, I had to place points at random, one in each of larger number of polygons (actually in objects of class 'SpatialPolygonsDataFrame' , see sp library), and tried first to do it using spsample (from sp). Surprisingly, every 5-15 trials, the output was a NULL value. The doc says that ' this may occur when trying to hit a small and awkwardly shaped polygon...

Am just feeling my way into the grid library, and cannot figure out how to define the plot limits. 3/5 of the example polygons below plot in the default 0-1 range viewport. But when I try to redefine the viewport the polygons plot in the same places. I also get the same result without employing push/pop. (As you can see from the scale I'm trying to introduce, I want to plot map polygons.) Grateful for an idea wher...