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)