0

我正在尝试在 R 中绘制我的研究区域。到目前为止，我已经下载了人口普查 OA 边界数据，并将属性数据合并到我希望查看的 OA 中。
因此，让我们从总共 1,000 个 OA 中说。我只想查看其中的 500 个 OA，我有外部数据集。Land 是我的空间多边形数据框，其中包含 1,000 个 OA。它有一列是 OA 名称“地理代码”。但是，它有它的坐标，因为你运行的是：

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，我必须做什么？感谢帮助

4

1 回答 1

0

您可以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)
于 2017-04-19T12:55:28.787 回答