Làm cách nào tôi có thể thực hiện các lựa chọn không gian hoặc bộ lọc trong R? Tôi có thể làm điều này trong PostGIS bằng cách sử dụng (ví dụ):
SELECT * FROM mypoint, mypoly
WHERE ST_Within(mypoint.geom, mypoly.geom);
và nó tương tự như một miếng bánh trong hầu hết các phần mềm GIS khác. Vì vậy, tôi mong đợi một cái gì đó trong R trông giống như:
mypoint[within(mypoint, mypoly),]
nhưng tôi không thể tìm thấy bất kỳ khách hàng tiềm năng ở bất cứ đâu. Tôi đang tìm kiếm thứ gì đó sử dụng thuật ngữ vị ngữ không gian chung được sử dụng trong PostGIS / GEOS / JTS / Shapely (bên trong, chứa, giao nhau, chồng chéo, chạm, v.v.; DWithin sẽ là phần thưởng). Có chức năng hoặc thư viện này tồn tại cho R? Nếu không, có cái gì đó tương tự?
Lưu ý: Trong R, tôi sử dụng sp
lớp cho các đối tượng của mình ( SpatialPointsDataFrame
, v.v.).
point.in.polygon
là một khởi đầu tốt cho các đa giác cơ bản, mặc dù hơi thô (ví dụ):point.in.polygon(pnt@coords[,1], pnt@coords[,2], lake@polygons[[1]]@Polygons[[1]]@coords[,1], lake@polygons[[1]]@Polygons[[1]]@coords[,2])
Tôi thực sự mong đợirgeos
rằng bạn đã đề cập và có tất cả các vị từ nhị phân và hơn thế nữa.