空間ポリゴンを外部データと結合し、属性のポリゴンを選択します

Question

私は研究領域をRにプロットしようとしています。これまでは国勢調査の境界データをダウンロードしましたが、私は見たいOAに属性データをマージしました。
したがって、合計1,000のOAから言えます。私は外部データセットを持っているこれらのOAの500だけを見たいと思います。 Landは、1,000sのOAを含む私の空間ポリゴンデータフレームです。これには、OA名の「ジオコード」という1つの列があります。ただし、あなたが実行されているので、それはその座標を持っています：

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の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)