1

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(...)

これらすべての形状（線とポリゴン）を1つのプロットを使用してプロットし、すべてのレイヤー、複数のレイヤーを1つのプロットにオーバーレイしたいと思います。それは可能ですか？

4

2 に答える 2

3

ggplot2 では、複数のジオメトリを積み重ねることができます。これは次のようになります。

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

これは、すべてのデータセットが同じ投影法にあり、すべてが x 座標などに列名 x を使用することを前提としています。ggplot2 は、sp オブジェクトではなく、data.frames で動作することに注意してください。関数 fortify を使用して、SpatialPolygons オブジェクトを data.frame に変換します。

于 2012-04-05T14:49:06.123 に答える
3

次のようなものを使用できます。

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

関数add=T内で使用するだけです。plot()

于 2012-04-05T10:58:48.297 に答える