Check if point is in spatial object which consists of multiple polygons/holes

by

You can do this simply with gContains(...) in the rgeos package.

gContains(sp1,sp2)

returns a logical depending on whether sp2 is contained within sp1. The only nuance is that sp2 has to be a SpatialPoints object, and it has to have the same projection as sp1. To do that, you would do something like this:

point <- data.frame(lon=10.2, lat=51.7)
sp2   <- SpatialPoints(point,proj4string=CRS(proj4string(sp1)))
gContains(sp1,sp2)

Here is a working example based on the answer to your previous question.

library(rgdal)   # for readOGR(...)
library(rgeos)   # for gContains(...)
library(ggplot2)

setwd("< directory with all your files >")
map <- readOGR(dsn=".", layer="vg250_gem", p4s="+init=epsg:25832")
map <- spTransform(map, CRS("+proj=longlat +datum=WGS84"))

nPolys <- sapply(map@polygons, function(x)length(x@Polygons))
region <- map[which(nPolys==max(nPolys)),]

region.df <- fortify(region)
points <- data.frame(long=c(10.5,10.51,10.15,10.4), 
                     lat =c(51.85,51.72,51.81,51.7),
                     id  =c("A","B","C","D"), stringsAsFactors=F)

ggplot(region.df, aes(x=long,y=lat,group=group))+
  geom_polygon(fill="lightgreen")+
  geom_path(colour="grey50")+
  geom_point(data=points,aes(x=long,y=lat,group=NULL, color=id), size=4)+
  coord_fixed()

Here, point A is in the main polygon, point B is in a lake (hole), point C is on an island, and point D is completely outside the region. So this code checks all of the points using gContains(...)

sapply(1:4,function(i)
  list(id=points[i,]$id,
       gContains(region,SpatialPoints(points[i,1:2],proj4string=CRS(proj4string(region))))))
#    [,1] [,2]  [,3] [,4] 
# id "A"  "B"   "C"  "D"  
#    TRUE FALSE TRUE FALSE

More Related Contents:

  1. How to efficiently calculate distance between pair of coordinates using data.table :=
  2. How to fill geom_polygon with different colors above and below y = 0 (or any other value)?
  3. Function to calculate geospatial distance between two points (lat,long) using R [duplicate]
  4. Different legends and fill colours for facetted ggplot?
  5. How to Convert data frame to spatial coordinates
  6. How to fill with different colors between two lines? (originally: fill geom_polygon with different colors above and below y = 0 (or any other value)?)
  7. Simple way to subset SpatialPolygonsDataFrame (i.e. delete polygons) by attribute in R
  8. Buffer (geo)spatial points in R with gbuffer
  9. R Plot Filled Longitude-Latitude Grid Cells on Map
  10. sp::over() for point in polygon analysis
  11. SpatialPolygons – Creating a set of polygons in R from coordinates
  12. How to plot a heat map on a spatial map
  13. Add raster to ggmap base map: set alpha (transparency) and fill color to inset_raster() in ggplot2
  14. Crop for SpatialPolygonsDataFrame
  15. Several questions on ggplot2 [closed]
  16. Error in if/while (condition) {: missing Value where TRUE/FALSE needed
  17. Interpreting “condition has length > 1” warning from `if` function
  18. How to change legend title in ggplot
  19. Importing Excel files into R, xlsx or xls
  20. Converting latitude and longitude points to UTM
  21. Solution. How to install_github when there is a proxy
  22. plot.new has not been called yet
  23. Count number of columns by a condition (>) for each row
  24. Remove all punctuation except apostrophes in R
  25. How to upgrade R in linux?
  26. no visible global function definition for ‘median’
  27. How to fix ‘tar: Failed to set default locale’ error?
  28. Duplicate list names in R
  29. Relocating Alaska and Hawaii on thematic map of the USA with ggplot2
  30. optimized rolling functions on irregular time series with time-based window

Leave a Comment