Fatal编程技术网
  • r/
  • R在空间多边形中查找空间多边形

R在空间多边形中查找空间多边形

r

R在空间多边形中查找空间多边形,r,R,我有一个美国住宅区多边形，它是该住宅区的简单轮廓。这是一个空间多边形数据帧 CRS("+proj=longlat +datum=WGS84")) CRS("+proj=longlat +datum=WGS84") 接下来，我有一个州（得克萨斯州）人口普查区SpatialPolygonsDataFrame CRS("+proj=longlat +datum=WGS84")) CRS("+proj=longlat +datum=WGS84") 我的目标是找到适合美国住宅区的人口普查区组 我的

我有一个美国住宅区多边形，它是该住宅区的简单轮廓。这是一个
空间多边形数据帧

CRS("+proj=longlat +datum=WGS84"))

CRS("+proj=longlat +datum=WGS84")
接下来，我有一个州（得克萨斯州）人口普查区
SpatialPolygonsDataFrame

CRS("+proj=longlat +datum=WGS84"))

CRS("+proj=longlat +datum=WGS84")
我的目标是找到适合美国住宅区的人口普查区组

我的方法似乎有严重缺陷

mysp <- raster::aggregate(rbind(SUBBLOCK_ALL.shp, CD32.shp))

mysp这里是一种使用北卡罗来纳州的
sf
包和附带的
nc
数据集的方法。这显然与你的问题不同，但前提应该是相同的。我们可以创建一个覆盖北卡罗来纳州一些县的大正方形“区”，然后使用
过滤器
和
st_in
仅保留边界内的那些县。您可以使用其他谓词，例如
st_intersects
，具体取决于您希望对位于边界上的谓词执行的操作。这也适用于更多的地区，因为
sf
geometry谓词实际上比较了两个参数中的所有几何体


库（tidyverse）
图书馆（sf）
#>链接至GEOS 3.6.1、GDAL 2.2.3、项目4.9.3
nc%
读取\u sf（）%>%
st_变换（3857）
地区%
st_证监会（crs=3857）
含%
过滤器（st_在（、地区、稀疏=假））
绘图（nc$几何图形）
绘图（地区，添加=真）
绘图（包含$geometry，add=TRUE，col=“red”）





由（v0.2.0）于2018年5月11日创建。
欢迎来到StackOverflow，请花点时间回顾一下教程：，如何创建一个最小的、完整的、可验证的示例：，以及如何提出好的问题-如果您在来SO寻求帮助之前对问题进行一些思考，您将获得更好的反馈和有用的答案，通过在代码中包含一个特定的问题，您已经尝试自己解决这个问题。请格式化您的问题。。。