Rの空間点データの周りにバッファを作成し、バッファ内のポイント数をカウントします。

Question

Rを使用してガソリンスタンドポイントの空間密度を分析しようとしています。ガソリンスタンドの周りにバッファ（円）バッファ内のガソリンスタンドの数。バッファー距離で遊んで、面白い何かを見るには合理的なバッファーが何であるかを見ていく必要があります。これらは私が作業しているファイルです：https://dl.dropboxusercontent.com/u/45095175/sbc_gas.shp; https://dl.dropboxusercontent.com/u/45095175/sbc_gas.shx; https://dl.dropboxusercontent.com/u/45095175/sbc_gas.dbf

# Install packages 
x <- c("ggmap", "rgdal", "rgeos", "maptools", "ks") 
lapply(x, library, character.only = TRUE) 
all <- readShapePoints("sbc_gas.shp") 
all.df <- as(all, "data.frame") 
locs <- subset(all.df, select = c("OBJECTID", "Latitude", "Longitude")) 
head(locs) # a simple data frame with coordinates 
coordinates(locs) <- c("Longitude", "Latitude") # set spatial coordinates 
plot(locs) 

すべてのヘルプは大歓迎!!

Answer 1

私の質問への回答が見つかりました：fivekm <- cbind(coordinates(locs), X=rowSums(distm (coordinates(locs)[,1:2], fun = distHaversine)/1000 <= 5)) # number of points within 5 km

Answer 2

.shpファイルだけでは十分ではないため、提供されたデータを使用することはできません。このデータをロードするには、少なくとも.shxファイルと.dbfファイルも提供する必要があります。

しかし、動作するはずのパッケージはgeosphereです。これにはという関数が含まれています。それを使って各ガソリンスタンドから他のすべてのステーションまでの距離を取得できます。距離マトリックスから、指定された距離内のすべてのステーションを選択することができます。