poly.clipper(UScensus2000)
poly.clipper()所属R语言包:UScensus2000
Selects the block, block group, or tracts which are contained within a given CDP.
选择块,块组,或大片都包含在一个给定的CDP。
译者:生物统计家园网 机器人LoveR
描述----------Description----------
poly.clipper allows the user to pull out all the blocks, block groups, or tracts that fall within the boundary of a given CDP and where the above grouping falls outside the boundaries poly.clipper will output the intersection of the two polygons and estimate the demographic variables based on the area of the intersection.
poly.clipper允许用户拉出所有的块,块组,或323属于一个给定的CDP的边界,其中上述分组不属于边界poly.clipper将输出交叉的两条多边形和估计人口统计变量的基础上的交点的区域。
用法----------Usage----------
poly.clipper(name, state, statefips = FALSE, level = c("tract", "blk", "blkgrp"), bb.epsilon = 0.006, sp.object = NULL, proj = NULL)
参数----------Arguments----------
参数:name
a character string, this string must be the name of CDP for a given state (if statefips=TRUE then this must be a CDP fips code).
一个字符串,这个字符串必须是名称的CDP对于一个给定的状态(如果statefips = TRUE,那么这必须是一个CDP FIPS代码)。
参数:state
a character string, can either be the full name of a state (e.g. "oregon"), the abbreviation (e.g. "or"), or the FIPS code (e.g. "41")– note that if you are using the FIPS code you have to change statefips to TRUE. This variable is insensitive to case.
一个字符串,可以是一个国家的全名(例如,“俄勒冈”),缩写(例如“或”),或的FIPS代码(例如“41”) - 请注意,如果您使用的是FIPS代码,你必须改变statefips到TRUE。这个变量是大小写不敏感的。
参数:statefips
logical, by default statefips=FALSE, change to TRUE when providing state with a FIPS code.
逻辑,默认情况下,statefips=FALSE,TRUE提供state用FIPS代码更改。
参数:level
character string, takes in one of three values: "tract", "blk", or "blkgrp". This defines the geographic level of data for the county.
字符的字符串,需要在以下三个值之一:“道”,“大牌”,或“blkgrp”。这定义县的GEO数据。
参数:bb.epsilon
numeric, by default bb.epsilon = 0.006, this value controls the number of Blocks/Block Groups/Tracts which are considered for clipping; if the function is not selecting some Blocks/Block Groups/Tracts then consider using a large value for bb.epsilon.
数字,默认情况下,bb.epsilon = 0.006,这个值控制块/块组/大港被认为是剪辑的数量;如果不选择一些功能块/块组/大港然后再考虑使用较大的值bb.epsilon。
参数:sp.object
SpatialPolygonsDataFrame, default NULL, allows the user to provide an sp object in which to perform this operation; primarily for use with demographics.add.
SpatialPolygonsDataFrame,默认NULL,允许用户提供了一个SP的对象,在其中执行此操作,主要用于与demographics.add。
参数:proj
CRS-class, takes a CRS object (e.g. CRS("+proj=utm +zone=10 +datum=NAD83") ); This is simply a wrapper for the spTransform function in rgdal . WARNING requires rgdal package.
CRS-级,需要CRS对象(如CRS(“+ PROJ = UTM +区域= 10 +基准= NAD83”)),这是一个简单的spTransform在rgdal函数的包装。警告需要rgdal包。
值----------Value----------
An object of class SpatialPolygonsDataFrame.
对象的类SpatialPolygonsDataFrame。
警告----------Warning ----------
You must have the packages UScensus2000blkgrp and UScensus2000blk installed to use levels "blkgrp" and "blk" respectively.
您必须拥有的包UScensus2000blkgrp和UScensus2000blk安装使用水平“blkgrp”和“BLK”。
(作者)----------Author(s)----------
Zack W. Almquist <a href="mailto:almquist@uci.edu">almquist@uci.edu</a>
参考文献----------References----------
Bureau, 2001. <br> http://www.census.gov/prod/cen2000/doc/sf1.pdf
参见----------See Also----------
city <br>
city参考
实例----------Examples----------
##Perform the clipping[#执行的剪裁]
portland<-poly.clipper(name="Portland",state="OR",level="tract")
######Plot Portland as a Choropleth map[#####绘制波特兰作为认识等值区域图]
############################################[###########################################]
## 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=unique(ceiling(quantile(dem))),dig.lab = 6)
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功能的图着色处理]
############################################[###########################################]
colors.use<-color.map(portland,portland$pop2000)
plot(portland,col=colors.use$colors)
#text(coordinates(alabama.blk),alabama.blk@data$name,cex=.3)[的文本(坐标(alabama.blk),alabama.blk的数据元的名称,CEX = 0.3)]
title(main="Census Tracts of\n Portland Oregon, 2000", sub="Quantiles (equal frequency)")
legend("bottomright",legend=colors.use$dem.cut,fill=colors.use$table.colors,bty="o",title="Population Count",bg="white")
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|