2012-04-05 4 views
1

マップファイルをプロットするためにRを使用しています。Rを使用して1つのフレームに複数のレイヤーを配置する方法?

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

Iは、単一プロット上の全ての層、複数の層を重ね、あるプロットを使用して、すべてのこれらの形状(ラインとポリゴン)をプロットしたいです。それは可能ですか?

答えて

3

次のようなものを使用することができます

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

ちょうど plot()関数内 add=Tを使用しています。

2

はggplot2では、複数のジオメトリを積み重ねることができます。

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

これは、すべてのデータセットはすべて、そのggplot2がdata.frames、ないSPのオブジェクトで動作しますx座標などのために列名xを使用し、同じ投影であることを前提としています。これは、何かがしたいと思います。関数fortifyを使用して、SpatialPolygonsオブジェクトをdata.frameに変換します。

関連する問題