Question

Im利用R绘制地图文档。

setwd("c:\apa\")
# Load libraries
library("maptools")
library("RColorBrewer")
library("classInt")
library("R2wd")
library("Hmisc")
library("plyr")
library("ggplot2")
library("rgdal")

file.path <- "c:\apa\"
file.name <- "apa_planalto_central"
apa <-readShapePoly(paste(file.path,file.name,sep=""))

file.name <- "AUTORIZACAO_09889652011"
auth <-readShapePoly(paste(file.path,file.name,sep=""))

file.name <- "CLIP_HIDROGRAFIA_APA"
hidro <-readShapeLines(paste(file.path,file.name,sep=""))

file.name <- "ferrovias"
ferrovias <-readShapeLines(paste(file.path,file.name,sep=""))

file.name <- "lagos_df"
lagos <-readShapePoly(paste(file.path,file.name,sep=""))

file.name <- "limite_df"
limites <-readShapeLines(paste(file.path,file.name,sep=""))

file.name <- "rodovias"
rodovias <-readShapeLines(paste(file.path,file.name,sep=""))


plot(apa)
plot(...)

我想用一块地划定所有这些形状(线和聚合物),把所有层层层推向一片土地上,多层层。这是可能的吗?

Answer 1

你们可以使用以下一些东西:

plot(apa)
plot(auth, add=T)

Just use the add=T within the plot ( function.

Answer 2

In ggplot2 you can stack multiple geometries. This would something like:

ggplot(aes(x = x, y = y), data = pointset1) +
  geom_point() +
  geom_polygon(aes(fill = z), data = polyset1) +
  etc

假设所有数据集都在同一预测中,所有数据集都使用X栏的名称,用于X坐标等。注:ggplot2使用数据、镜子、而不是 sp。利用起义功能将空间工程物体转化为数据。

友情链接