Question

我想使用spplot + sp.lines（点阵）代替plot + segments。你知道一个简单的方法来实现这个，例如R: Calculating the shortest distance between two point layers

library(dismo)  
require(rgdal)
require(FNN)

laurus <- gbif("Laurus", "nobilis")
locs <- subset(laurus, !is.na(lat) & !is.na(lon),
               select = c("country", "lat", "lon"))

locs.uk  <- subset(locs, locs$country=="United Kingdom")
locs.ire <- subset(locs, locs$country=="Ireland")

uk_coord <- SpatialPoints(locs.uk[,c("lon","lat")])
ire_coord <- SpatialPoints(locs.ire[,c("lon","lat")])
crs.geo<-CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")  
proj4string(uk_coord) <- crs.geo 
proj4string(ire_coord) <- crs.geo 

uk_coord  <- spTransform(uk_coord, CRS("+init=epsg:27700"))
ire_coord <- spTransform(ire_coord, CRS("+init=epsg:27700"))


g = get.knnx(coordinates(uk_coord), coordinates(ire_coord),k=1)

可视化此

plot(uk_coord, col=2, xlim=c(-1e5,6e5))
plot(ire_coord, add=TRUE)
segments(coordinates(ire_coord)[,1], 
         coordinates(ire_coord)[,2], 
         coordinates(uk_coord[g$nn.index[,1]])[,1], 
         coordinates(uk_coord[g$nn.index[,1]])[,2])

可能会转换为类似

的内容

ire <- list("sp.points", ire_coord)

spplot(uk_coord, sp.layout=list(ire))

但是有一种简单的方法可以将segments转换为SpatialLines，即list("sp.lines", Lines(...))

Answer 1

从panel.segments() - 包中尝试lattice：

library("lattice")
spplot(rbind(uk_coord, ire_coord), auto.key=FALSE,
       panel=function(...) {
         panel.xyplot(...)
         panel.segments(coordinates(ire_coord)[,1], 
                        coordinates(ire_coord)[,2],
                        coordinates(uk_coord[g$nn.index[,1]])[,1],
                        coordinates(uk_coord[g$nn.index[,1]])[,2])
       })

Answer 2

了解面板功能比依赖sp.layout中的spplot更强大 - 因此直接使用lattice或grid功能。 sp.layout的解决方案可能如下所示：

spplot(uk_coord, auto.key=FALSE, col.regions = 'black', 
  sp.layout = list(ire, 
    list("panel.segments", 
        coordinates(ire_coord)[,1],
        coordinates(ire_coord)[,2],
        coordinates(uk_coord[g$nn.index[,1]])[,1],
        coordinates(uk_coord[g$nn.index[,1]])[,2])), 
    xlim = c(-140000,700000))

请注意，它不仅限于sp.lines等功能;在即将发布的sp 1.1-0中，函数名称的引号也可以省略。

spplot尝试默认情况下用颜色绘制要素的属性，这在这里没有意义，所以你基本上想要的是具有受控纵横比（xyplot）的asp="iso"。 / p>

SpatialLines而不是segment（spplot）

2 个答案: