Question

我有一个数据框，其中包含有关犯罪（变量x）以及犯罪发生地的纬度和经度的信息。我有一个与圣保罗市各区的形状文件。我需要合并这两个数据，这样我才能得出每个地区的犯罪总数。有没有办法做到这一点？我使用

将数据框转换为空间数据框

df.sp <- SpatialPointsDataFrame(cbind(df$longitude,df$latitude ), df)

但是我不知道如何实现这种合并以获得我所需要的。在df上，我有1万多个obs，例如：

    latitude  longitude  n_homdol
1    -23.6     -46.6     1
2    -23.6     -46.6     1
3    -23.6     -46.6     1
4    -23.6     -46.6     1
5    -23.6     -46.6     1
6    -23.6     -46.6     1

形状文件如下：

                       geometry      NOME_DIST
1 POLYGON ((352436.9 7394174,... JOSE BONIFACIO
2 POLYGON ((320696.6 7383620,...    JD SAO LUIS
3 POLYGON ((349461.3 7397765,...    ARTUR ALVIM
4 POLYGON ((320731.1 7400615,...        JAGUARA
5 POLYGON ((338651 7392203, 3...  VILA PRUDENTE
6 POLYGON ((320606.2 7394439,...        JAGUARE

我需要按地区划分n_homdol的总和。我正在尝试合并两个数据框，但是没有成功。

Answer 1

如果您愿意从sp切换到sf软件包，那么您将有一种简单的方法，可以使用类似dplyr的语法进行空间连接：{{ 1}}。

它会像这样工作（我的电脑上没有R，所以可能会有一些“笔滑”）

st_join

您可以执行以下操作：


library(sf)
library(dplyr)

#Instead of data.frame of class "sp", create "simple features"-data.frame 

sf_df = st_as_sf(df, coords = c("longitude", "latitude"), crs = 4326)

#You'll have to convert your shapefile to sf, too. 
#Depending what class it is you can use "st_as_sf()"


#Then join the shapefile with sf_df via the "st_contains" which merges two rows 
#if a point from sf_df falls within a polygon from the shapefile.

shape_df <- st_join(shapefile, sf_df , join = st_contains)

如果您想坚持使用shape_df %>% group_by(NOME_DIST) %>% summarise(crime = sum(n_homdol))，建议您在评论中的Dave2e链接中查看答案。

合并具有经度和纬度的数据框以及一个形状文件

1 个答案: