sp包—aggregate函数

最新推荐文章于 2021-06-25 16:47:43 发布

HHHHYN

最新推荐文章于 2021-06-25 16:47:43 发布

阅读量294

点赞数

分类专栏： R语言函数

本文链接：https://blog.csdn.net/weixin_52550583/article/details/118164055

版权

R语言函数专栏收录该内容

19 篇文章 2 订阅

订阅专栏

描述

空间对象聚合

空间对象中主题信息的空间聚合

用法

## S3 method for class 'Spatial'
aggregate(x, by = list(ID = rep(1, length(x))),
	FUN, ..., dissolve = TRUE, areaWeighted = FALSE)

参数

x：对象派生自空间，具有属性

by：聚合谓词;如果by是一个Spatial对象，则x中的属性被聚合的几何形状;如果by是一个列表，则按属性聚合，参见aggregate.data.frame

FUN：聚集函数，如mean;详见细节

...：参数传递给函数FUN，除非minDimension被指定，它被传递给函数

dissolve：逻辑值;当基于属性进行聚合时，生成的几何图形应该被分解吗?注意，如果x有SpatialPointsDataFrame类，则返回一个SpatialMultiPointsDataFrame类的对象

areaWeighted：逻辑值;x的总和是否应该用它与每个特征的交点的面积来加权?见值

细节

FUN应该是一个函数，它的第一个参数是一个向量，并返回单个数字。典型的例子是均值和求和。计数特征是在对一个处处为1的属性变量求和时获得的。

值

属性值x的聚合要么通过对几何图形的over进行空间匹配，要么通过对属性值的聚合，使用聚合函数FUN。如果areaWeighted为TRUE，则忽略FUN，并计算数值变量的面积加权平均值，或者如果所有属性都是因子，则返回区域主导因子水平(区域模式)。这将计算x和by的gIntersection;请参见下面的例子。

如果缺少，则聚合所有特性。

例子

data("meuse")
coordinates(meuse) <- ~x+y
data("meuse.grid")
coordinates(meuse.grid) <- ~x+y
gridded(meuse.grid) <- TRUE
i = cut(meuse.grid$dist, c(0,.25,.5,.75,1), include.lowest = TRUE)
j = sample(1:2, 3103,replace=TRUE)
## Not run: 
if (require(rgeos)) {
	# aggregation by spatial object:
	ab = gUnaryUnion(as(meuse.grid, "SpatialPolygons"), meuse.grid$part.a)
	x = aggregate(meuse["zinc"], ab, mean)
	spplot(x)
	# aggregation of multiple variables
	x = aggregate(meuse[c("zinc", "copper")], ab, mean)
	spplot(x)
	# aggregation by attribute, then dissolve to polygon:
	x = aggregate(meuse.grid["dist"], list(i=i), mean)
	spplot(x["i"])
	x = aggregate(meuse.grid["dist"], list(i=i,j=j), mean)
	spplot(x["dist"], col.regions=bpy.colors())
	spplot(x["i"], col.regions=bpy.colors(4))
	spplot(x["j"], col.regions=bpy.colors())
}

## End(Not run)

x = aggregate(meuse.grid["dist"], list(i=i,j=j), mean, dissolve = FALSE)
spplot(x["j"], col.regions=bpy.colors())

if (require(gstat) && require(rgeos)) {
	x = idw(log(zinc)~1, meuse, meuse.grid, debug.level=0)[1]
	spplot(x[1],col.regions=bpy.colors())
	i = cut(x$var1.pred, seq(4, 7.5, by=.5), 
		include.lowest = TRUE)
	xa = aggregate(x["var1.pred"], list(i=i), mean)
	spplot(xa[1],col.regions=bpy.colors(8))
}

if (require(rgeos)) {
# Area-weighted example, using two partly overlapping grids:

  gt1 = SpatialGrid(GridTopology(c(0,0), c(1,1), c(4,4)))
  gt2 = SpatialGrid(GridTopology(c(-1.25,-1.25), c(1,1), c(4,4)))

  # convert both to polygons; give p1 attributes to aggregate
  p1 = SpatialPolygonsDataFrame(as(gt1, "SpatialPolygons"), 
		  data.frame(v = 1:16, w=5:20, x=factor(1:16)), match.ID = FALSE)
  p2 = as(gt2, "SpatialPolygons")

  # plot the scene:
  plot(p1, xlim = c(-2,4), ylim = c(-2,4))
  plot(p2, add = TRUE, border = 'red')
  i = gIntersection(p1, p2, byid = TRUE)
  plot(i, add=TRUE, density = 5, col = 'blue')
  # plot IDs p2:
  ids.p2 = sapply(p2@polygons, function(x) slot(x, name = "ID"))
  text(coordinates(p2), ids.p2)
  # plot IDs i:
  ids.i = sapply(i@polygons, function(x) slot(x, name = "ID"))
  text(coordinates(i), ids.i, cex = .8, col = 'blue')

  # compute & plot area-weighted average; will warn for the factor
  ret = aggregate(p1, p2, areaWeighted = TRUE)
  spplot(ret)

  # all-factor attributes: compute area-dominant factor level:
  ret = aggregate(p1["x"], p2, areaWeighted = TRUE) 
  spplot(ret)
}