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. 여기에는 distGeo이라는 기능이 있습니다. 각 주유소에서 다른 모든 주유소까지의 거리를 구할 수 있습니다. 거리 매트릭스에서 지정된 거리 내의 모든 스테이션을 선택할 수 있어야합니다.