标签：ME 语言 Area median s1 矢量 medinc 地理 Income

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)

标签：ME,语言,Area,median,s1,矢量,medinc,地理,Income
来源： https://www.cnblogs.com/cqy-wt1124/p/15123310.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。