R/GIS：使用其他来源的数据覆盖GADM数据？

Question

我有Global Administrative Areas Database的行政边界数据。用这些数据绘制细分是非常容易的。例如，如果我想overplot包含在文件mydata在佛罗里达县地理编码数据，我下载了美国的“2级”，然后使用此代码：

load("/home/anindya/Desktop/DELETE/USA_adm2.RData") 
temp = as.data.frame(gadm) *Code for eyeballing structure of data 
florida = gadm[344:410,] 
plot(florida); points(mydata$longitude, mydata$latitude) 

但你如何覆盖其他这方面的信息，比如说地形和水路？我对Natural Earth和Global Lakes and Wetlands Database的数据特别感兴趣。

后续编辑：非常感谢保罗解决问题。我已经在here at gis.stackexchange上转换了这个问题的一个变体，其中R.K.给出了很好的答案。请仔细阅读。

Answer 1

如果我正确理解了你，你会问如何叠加不同的地理层，就像使用GUI GIS一样。

这是直截了当的，如果数据是在一个Spatial*对象（例如SpatialPoints, SpatialLines, SpatialPolygons等从sp封装）或作为来自raster封装的Raster*对象。如果给出参数add=TRUE，则这两个包中的plot()方法都可以处理此叠加。

要将基于矢量的地理数据导入Spatial*对象，可以使用maptools包中的函数（例如readShapePoly()读取ESRI多边形shape文件）。可以使用raster()加载各种文件规格的光栅，如果安装了rgdal软件包，则可以使用其他格式。

以下是使用拼凑的地理数据叠加地理数据类型的示例。首先创建一个UTM网格一些示例数据：

library(sp) 
library(raster) 

## Create a RasterLayer object and fill with random values 
baseRaster <- raster(nrow=100, ncol=100, 
    xmn=0, xmx=100, ymn=0, ymx=100, 
    crs= "+proj=utm") 
baseRaster[] <- runif(ncell(baseRaster)) 

## Create a second raster to superimpose 
## It should contain NA values where it is transparent 
supRaster <- baseRaster 
supRaster[] <- NA 
supRaster[cellFromCol(supRaster, 48:52)] <- 1 


## Create SpatialPoints object to superimpose on these 
loc <- SpatialPoints(cbind(seq(10, 90, by=10), seq(10, 90, by=10))) 

现在做图：

## Plot base raster 
plot(baseRaster) 

## Superimpose second raster in a different colour 
## Turn off legend 
plot(supRaster, add=TRUE, col="blue", legend=FALSE) 

## Superimpose points and make them big and colourful 
plot(loc, add=TRUE, pch=20, cex=3, col="red") 

下面是结果：