areaPoly(UScensus2000)
areaPoly()所属R语言包:UScensus2000
Area of the polygons in SpatialPolygonsDataFrame
多边形SpatialPolygonsDataFrame区
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Calculates the area of each polygon in SpatialPolygonsDataFrame.
计算面积中每个多边形SpatialPolygonsDataFrame。
用法----------Usage----------
areaPoly(sp.object=NULL)
参数----------Arguments----------
参数:sp.object
SpatialPolygonsDataFrame, must be a SpatialPolygonsDataFrame object.
SpatialPolygonsDataFrame,必须一个SpatialPolygonsDataFrame对象。
值----------Value----------
a numeric vector.
一个数值向量。
(作者)----------Author(s)----------
Zack W. Almquist
参考文献----------References----------
实例----------Examples----------
############################################[###########################################]
## Helper function for handling coloring of the map[#Helper功能的图着色处理]
############################################[###########################################]
color.map<- function(x,dem,y=NULL){
l.poly<-length(x@polygons)
dem.num<- cut(dem,breaks=ceiling(quantile(dem)),dig.lab = 10)
dem.num[which(is.na(dem.num)==TRUE)]<-levels(dem.num)[1]
l.uc<-length(table(dem.num))
if(is.null(y)){
##commented out, but creates different color schemes [#注释掉,而造成不同的配色方案]
## using runif, may take a couple times to get a good color scheme.[#的runif,可能需要几次得到一个很好的配色方案。]
##col.heat<-rgb( runif(l.uc,0,1), runif(l.uc,0,1) , runif(l.uc,0,1) )[col.heat <-RGB(runif(l.uc,0,1),runif(l.uc,0,1),runif(l.uc,0,1))]
col.heat<-heat.colors(16)[c(14,8,4,1)] ##fixed set of four colors[#固定的四种颜色]
}else{
col.heat<-y
}
dem.col<-cbind(col.heat,names(table(dem.num)))
colors.dem<-vector(length=l.poly)
for(i in 1:l.uc){
colors.dem[which(dem.num==dem.col[i,2])]<-dem.col[i,1]
}
out<-list(colors=colors.dem,dem.cut=dem.col[,2],table.colors=dem.col[,1])
return(out)
}
############################################[###########################################]
## Helper function for handling coloring of the map[#Helper功能的图着色处理]
############################################[###########################################]
data(oregon.tract)
## Build density using areaPoly()[#建立密度,使用areaPoly()]
den00<-oregon.tract$pop2000/areaPoly(oregon.tract)
oregon.tract$den00<-den00
colors.use<-color.map(oregon.tract,oregon.tract$den00)
plot(oregon.tract,col=colors.use$colors,border="transparent")
title(main="Census Tracts \n of Oregon, 2000", sub="Quantiles (equal frequency)")
legend("bottomright",legend=colors.use$dem.cut,fill=colors.use$table.colors,bty="o",title="Population Density",bg="white")
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|