Hi All,
I'm somewhat of a novice at doing spatial analysis in R, so forgive me if my
question is obvious, or misposted.
I have a number of points located within the United States. My goal is to
compute the voronoi tessellation, bounded by the borders of the US. I would
like to eventually format this as a dataframe of coordinates with a grouping
ID column so that I can plot the polygons in ggplot2.
Here's what I've got so far.
cities <- structure(list(CitState = c("SiouxFalls SD", "Rockford IL",
"Kenosha WI",
"Duluth MN", "SiouxCity IA", "CedarRapids IA", "NewYork NY"),
Latitude = c(43.54599, 42.265973, 42.577791, 46.779135, 42.489678,
41.976662, 40.756054), Longitude = c(-96.731291, -89.086667,
-87.822644, -92.108243, -96.404948, -91.673155, -73.986951
)), .Names = c("CitState", "Latitude", "Longitude"), row.names = c(1L,
2L, 3L, 4L, 5L, 6L, 267L), class = "data.frame")
library(maps)
library(spatstat)
us <- data.frame(map("usa", regions = "main",plot = F)[c("x","y")])
manhattan <- data.frame(map("usa", regions = "manhattan",plot =
F)[c("x","y")])
spatstat.options(checkpolygons = FALSE)
us.win <- owin(poly = list(list(x = us$x, y = us$y),list(x = manhattan$x,y =
manhattan$y)))
us.ppp <- ppp(cities$Longitude, cities$Latitude, window = us.win)
But, if I use dirichlet() on us.ppp, it gives me back a pixel array. If I
use deldir() on us.ppp, the boundaries of the tesselations don't seem to be
constrained by the US borders, and I can't figure out how to extract the
polygons from the object.
plot(dirichlet(us.ppp))
##
map("usa")
plot(deldir(us.ppp), wlines = "tess",add = T)
If I could extract from some object, a list of polygons, I could probably
get the rest done myself.
Thanks in advance!
-Joe
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
