s1 <- readRDS("Income_schooling.rds")

s2 <- readRDS("dist.rds")

l1 <- readRDS("highway.rds")
#函数bbox获取空间数据边界，sp包
b1<-bbox(s1)
b2<-bbox(s2)
b3<-pmax(b1,b2)
b3[,1]<-pmin(b1[,1],b2[,1])

#aggregate()函数，融合边界,合为一体，raster包
ME<-aggregate(s1)

s1$med<-s1$Income>median(s1$Income)
#以某种属性进行融合
ME.inc<-aggregate(s1,by="med")
#group_by()对s1@data进行以med属性分组，然后summarize()统计各组的中位数
#.$medinc将统计结果赋值给
ME.inc$Income <- s1@data %>% group_by(med) %>%
  summarize(medinc = median(Income)) %>% .$medinc
#计算面积
ME.inc$Area<-gArea(ME.inc,byid=TRUE)/1000000

#intersect()截取相交部分,raster
clp1<-intersect(s1,s2)
#按照输入数据顺序，来决定输出数据类型，clp2为线空间数据
clp2<-intersect(l1,s1)
#clp3为面空间数据
clp3<-intersect(s1,l1)


#union()相交融合，raster包
un1<-union(s1,s2)