r-与外部数据合并的空间聚合物，然后在属性上选择多边形

Question

我正在尝试在R中绘制我的研究区域。因此可以说，总共有1,000个OA。我只想查看其中的500个OA，其中我具有外部数据集。土地是我的空间多边形数据框架，其中包含1,000次OA。它具有OA名称" GeoCode"的一列。但是，它具有其坐标，因为您运行：

coordinates(land)

它产生所有坐标。

houseprice是一个包含房价的数据框架和带有地理位置的专栏。

远面的代码：

library("sp")
library("GISTools")
library("maptools")
library("spatial")
library("ggplot2")
library("rgdal")
land@data = merge(land@data,houseprice,by.x="geocode", 
by.y="geocode",sort=FALSE)
View(land)

这表明我只有500个条目(这是我想要的(。

plot(land)

这最终绘制了整个1,000个OA，而不仅仅是我想要的500个。

我只为500个OA，只需绘制数据才能做什么？感谢您的帮助

Answer 1

您可以直接在空间数据上使用merge，因为有一种方法。然后，您仅选择添加列中没有NA的行(在下面的示例中称为my_col(。
在这里，您仅修改了属性表，而不是几何

# Merge with external dataset
land2 <- merge(land, houseprice,by.x="geocode", by.y="geocode",sort=FALSE, all.x = TRUE)
# Select only rows of interest
land500 <- land2[which(!is.na(land2$my_col)),]
plot(land500)