10

还不如让猫从袋子里出来，说我已经尝试了大约一天来弄清楚当我要求plot一个类的对象时发生了什么SpatialPolygons*。我的直觉是寻找适用plot.SpatialPolygonsDataFrame于我见过和遇到的任何其他示例的方法，但这失败了。

快速示例：

美国各州的 .shp 文件在此处。

library(maptools)
us.states<-readShapePoly("cb_2014_us_state_5m.shp")
plot(us.states)

一个示例图。 它显示了放置在一个看不见的地球上的所有美国州和地区的多边形。 很难阅读，但情节不是这个问题的主要内容。

有点笨拙，但无论如何——关键是一个简单的调用plot和某个深处的某个地方R知道调用哪个方法，以便为我们提供所有形状的表示。

在搜索了一下之后，我尝试了各种 SO Q&As 建议的其他两种方法，首先是：

> methods(plot)
 [1] plot.aareg*                         plot.acf*                          
 [3] plot,ANY,ANY-method                 plot.cox.zph*                      
 [5] plot.data.frame*                    plot.decomposed.ts*                
 [7] plot.default                        plot.dendrogram*                   
 [9] plot.density*                       plot.ecdf                          
[11] plot.factor*                        plot.formula*                      
[13] plot.function                       plot.hclust*                       
[15] plot.histogram*                     plot.HoltWinters*                  
[17] plot.isoreg*                        plot.lm*                           
[19] plot.medpolish*                     plot.mlm*                          
[21] plot.ppr*                           plot.prcomp*                       
[23] plot.princomp*                      plot.profile.nls*                  
[25] plot.raster*                        plot.shingle*                      
[27] plot,SpatialGrid,missing-method     plot,SpatialLines,missing-method   
[29] plot,Spatial,missing-method         plot,SpatialPixels,missing-method  
[31] plot,SpatialPoints,missing-method   plot,SpatialPolygons,missing-method
[33] plot.spec*                          plot.spline*                       
[35] plot.stepfun                        plot.stl*                          
[37] plot.survfit*                       plot.table*                        
[39] plot.times*                         plot.trellis*                      
[41] plot.ts                             plot.tskernel*                     
[43] plot.TukeyHSD*                      plot.xyVector*                     
[45] plot.zoo*

答案似乎在插槽 27-32 中；然而，每一个后面都有一个令人烦恼的逗号！没有提到这是可能的?methods并且?plot,SpatialPolygons,missing-method是一个错误。快速搜索missing-method没有任何用处，并且在?methods.

好的; 继续。怎么样getS3method？getMethod? getMethods? getAllMethods??

> getS3method("plot","SpatialPolygonsDataFrame")
Error in getS3method("plot", "SpatialPolygonsDataFrame") : 
  S3 method 'plot.SpatialPolygonsDataFrame' not found

> getMethod("plot","SpatialPolygonsDataFrame")
Error in getMethod("plot", "SpatialPolygonsDataFrame") : 
  no method found for function 'plot' and signature SpatialPolygonsDataFrame

后两者已被弃用并且也不返回任何内容。

显然，这些功能只是我第一直觉的替代品。

那么现在怎么办？我如何知道plot在传递 a 时调用了哪个方法SpatialPolygonsDataFrame？是否有任何通用方法可以取代我上面使用的方法来解决这个问题？

编辑：

我偶然发现了这个（= ?`Spatial-Polygons-class` ），上面写着：

空间多边形的plot方法采用以下参数：

但它仍然没有确切地说该方法是什么。

4

2 回答 2

6

部分答案：方法是sp:::plot.SpatialPolygons。我不知道为什么这些逗号会出现在 ; 的结果中methods(plot)。plot,SpatialPolygons,missing-method对我来说毫无意义（我也看到了）。

我通过导致错误找到了答案。我没有你的形状文件，所以我从 中获取了示例help("SpatialPolygonsDataFrame-class")，我从中得到了?SpatialPolygons"。示例如下：

Sr1 = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))
Sr2 = Polygon(cbind(c(5,4,2,5),c(2,3,2,2)))
Sr3 = Polygon(cbind(c(4,4,5,10,4),c(5,3,2,5,5)))
Sr4 = Polygon(cbind(c(5,6,6,5,5),c(4,4,3,3,4)), hole = TRUE)
Srs1 = Polygons(list(Sr1), "s1")
Srs2 = Polygons(list(Sr2), "s2")
Srs3 = Polygons(list(Sr3, Sr4), "s3/4")
SpP = SpatialPolygons(list(Srs1,Srs2,Srs3), 1:3)
plot(SpP, col = 1:3, pbg="white")

然后我更改了绘图调用以给出错误：

plot(SpP, col = kasjdhfkjasdfhkas, pbg = "white")
# Error in plot.SpatialPolygons(x, ...) : object 'kasjdhfkjasdfhkas' not found

（除非您kasjdhfkjasdfhkas在工作区中将其定义为向量，否则应该可以使用）。果然，

sp:::plot.SpatialPolygons

显示了该方法的代码。

于 2015-09-04T01:51:24.040 回答
5

回答 Gregor 提出的几个问题。该methods函数以前不显示 S4 方法的结果，但现在可以了。在 NEWS 文档中进行搜索，我发现这是在 3.2.0 版中添加的：

methods() reports S4 in addition to S3 methods; output is simplified when the class 
argument is used. .S3methods() and methods::.S4methods() report S3 and S4 methods separately.

SpatialPolygon 是 S4 对象，因此具有调度的插槽和 S4 方法。您可以通过以下方式显示 S4 方法：

showMethods(f='plot', classes='SpatialPolygons', includeDefs=TRUE)
#    ---- result ---
Function: plot (package graphics)
x="SpatialPolygons", y="missing"
function (x, y, ...) 
plot.SpatialPolygons(x, ...)

这告诉您该类有一个 S4 函数。您可以执行sp:::SpatialPolygons或getAnywhere(plot.SpatialPolygons)查看功能代码。（通常调用 usingincludeDefs=TRUE会显示 R 代码，但此处不显示。）methods(plot)加载包“sp”时的结果告诉您，工作区中注册了 8 种不同的绘图 S4 方法。逗号后面的项目是用于函数调度的“签名”。“缺失方法”是指未指定第二个参数的情况，其中在missing函数体内执行的函数将返回TRUE：

  [1] plot,ANY,ANY-method                 plot,color,ANY-method              
  [3] plot,Spatial,missing-method         plot,SpatialGrid,missing-method    
  [5] plot,SpatialLines,missing-method    plot,SpatialPixels,missing-method  
  [7] plot,SpatialPoints,missing-method   plot,SpatialPolygons,missing-method

这将向您展示plot当传递给它的对象与任何 S4 方法不匹配时，然后用于分派 -call 的内容：

showMethods(f='plot', classes='ANY', includeDefs=TRUE)
Function: plot (package graphics)
x="ANY", y="ANY"
function (x, y, ...) 
UseMethod("plot")


x="color", y="ANY"
function (x, y, ...) 
{
    .local <- function (x, y, pch = 20, cex = 3) 
    pairs(coords(x), col = hex(x, fix = TRUE), pch = pch, cex = cex)
    .local(x, y, ...)
}

基本上，R 解释器首先检查是否有任何 S4 方法是合适的，如果失败，将开始遍历 S3 方法，直到它到达plot.default.

> getMethod(f='plot', signature=c(x='SpatialPolygons', y='missing'))
Method Definition:

function (x, y, ...) 
plot.SpatialPolygons(x, ...)
<environment: namespace:sp>

Signatures:
        x                 y        
target  "SpatialPolygons" "missing"
defined "SpatialPolygons" "missing"

顺便说一句，我在以下回复中看到了这一点?getMethods：

## Deprecated in 2010 and defunct in 2015 for \code{table = FALSE}:
getMethods(f, where, table = FALSE)
于 2015-09-04T03:19:52.317 回答