Over a million developers have joined DZone.
{{announcement.body}}
{{announcement.title}}

Spatial Overlays with R - Retrieving Polygon Attributes for a Set of Points

DZone's Guide to

Spatial Overlays with R - Retrieving Polygon Attributes for a Set of Points

· Big Data Zone ·
Free Resource

Hortonworks Sandbox for HDP and HDF is your chance to get started on learning, developing, testing and trying out new features. Each download comes preconfigured with interactive tutorials, sample data and developments from the Apache community.

A short tutorial for spatial overlays using R-GIS...

library(sp)
library(dismo)

# spatial data
alt <- getData('alt', country = "AT")
gadm <- getData('GADM', country = "AT", level = 2)

# view
plot(alt)
plot(gadm, add=T)

# some addresses
pts <- geocode(c("Aldrans, Grubenweg", "Wien, Stephansdom", "Salzburg, Mozartplatz"))

# make it spatial
coords <- SpatialPoints(pts[, c("longitude", "latitude")])
spdf_pts <- SpatialPointsDataFrame(coords, pts)

# assign CRS/projection
proj4string(spdf_pts) <- CRS(" +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")

# check data
str(spdf_pts)

# plot it on top
plot(spdf_pts, cex = 2, col = 2, add = T)

# do an intersection (points in polygon)
# yielding the polygon's attribute data
over(spdf_pts, gadm)

Hortonworks Community Connection (HCC) is an online collaboration destination for developers, DevOps, customers and partners to get answers to questions, collaborate on technical articles and share code examples from GitHub.  Join the discussion.

Topics:

Published at DZone with permission of

Opinions expressed by DZone contributors are their own.

{{ parent.title || parent.header.title}}

{{ parent.tldr }}

{{ parent.urlSource.name }}