如何解决如何将空间点数据合并到 R 中的多边形?
我正在尝试在 SpatialPointsDataFrame
和 SpatialpolygonsDataFrame
之间执行空间连接。
我有两组数据,名为 California_shp
和 Houses
。他们的信息如下:
dput(head(California_shp))
structure(list(NAME = 1:6),row.names = c(NA,6L),class = "data.frame")
class(California_shp)
[1] "SpatialpolygonsDataFrame"
attr(,"package")
[1] "sp"
> dput(head(Houses))
structure(list(houseValue = c(452600L,358500L,352100L,341300L,342200L,269700L),income = c(8.3252,8.3014,7.2574,5.6431,3.8462,4.0368),houseAge = c(41L,21L,52L,52L),rooms = c(880L,7099L,1467L,1274L,1627L,919L),bedrooms = c(129L,1106L,190L,235L,280L,213L),population = c(322L,2401L,496L,558L,565L,413L),households = c(126L,1138L,177L,219L,259L,193L),latitude = c(37.88,37.86,37.85,37.85),longitude = c(-122.23,-122.22,-122.24,-122.25,-122.25)),class = "data.frame")
class(Houses)
[1] "SpatialPointsDataFrame"
attr(,"package")
[1] "sp"
如何使用point-in-polygon
操作和mean
函数将Houses
对象合并到California_shp
?
我已经在 R 中编写并运行了以下代码:
library(spatialEco)
library(sp)
head(California_shp@data)
NAME
1 Alameda
2 Alpine
3 Amador
4 Butte
5 Calaveras
6 Colusa
head(Houses@data)
houseValue income houseAge rooms bedrooms population households latitude longitude
1 452600 8.3252 41 880 129 322 126 37.88 -122.23
2 358500 8.3014 21 7099 1106 2401 1138 37.86 -122.22
3 352100 7.2574 52 1467 190 496 177 37.85 -122.24
4 341300 5.6431 52 1274 235 558 219 37.85 -122.25
5 342200 3.8462 52 1627 280 565 259 37.85 -122.25
6 269700 4.0368 52 919 213 413 193 37.85 -122.25
pts.poly <- point.in.poly(Houses,California_shp)
pts.poly
class : SpatialPointsDataFrame
features : 20640
extent : -124.35,-114.31,32.54,41.95 (xmin,xmax,ymin,ymax)
crs : +proj=longlat +datum=NAD83 +no_defs
variables : 10
names : houseValue,income,houseAge,rooms,bedrooms,population,households,latitude,longitude,NAME
min values : 14999,0.4999,1,2,3,-124.35,1
max values : 500001,15.0001,52,39320,6445,35682,6082,41.95,58
不幸的是,它仍然对我没有帮助。我需要改进它们。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。