Question

我正在努力制作一个12个生态区域的图，这些区域来自shapefile。举个例子：

ER_10.2<-Level.2.ecoregs[Level.2.ecoregs$NA_L2CODE=="10.2",]
> ER_10.2
class       : SpatialPolygonsDataFrame 
nfeatures   : 26 
extent      : -1693158, 44930.55, -2591002, -691719.8  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=laea +lat_0=45 +lon_0=-100 +x_0=0 +y_0=0 +a=6370997 +b=6370997 +units=m +no_defs 
nvariables  : 8
names       : NA_L2CODE,    NA_L2NAME, NA_L1CODE,              NA_L1NAME,           NA_L2KEY,                   NA_L1KEY,  Shape_Leng,   Shape_Area 
min values  :      10.2, WARM DESERTS,        10, NORTH AMERICAN DESERTS, 10.2  WARM DESERTS, 10  NORTH AMERICAN DESERTS,    11613.69,      8382714 
max values  :      10.2, WARM DESERTS,        10, NORTH AMERICAN DESERTS, 10.2  WARM DESERTS, 10  NORTH AMERICAN DESERTS, 11456404.58, 510159399963

我需要在循环中执行此操作，因为我还包括其他分析。

Ecoregions.list <- c("ER_10.2", "ER_12.1", "ER_14.3","ER_13.2",
  "ER_09.6", "ER_09.5", "ER_14.1", "ER_13.3", "ER_09.4", "ER_08.3", "ER_13.1", "ER_11.1")

Ecoregions<-unique(as.character(Ecoregions.list))

for(i in 1:length(Ecoregions))
    {
    Ecoregions<-unique(as.character(Ecoregions.list))
    ER=as.name(Ecoregions[i])
    plot (ER)
    }

但是当我试图在图中读取它时，我总是得到这个错误：

Error in xy.coords(x, y, xlabel, ylabel, log) : invalid first argument

有人对如何解决此问题有任何建议吗？

提前致谢！

Answer 1

我不确定循环是否是这种情况下的最佳做法，但在您当前的代码中使用as.name会让您感到悲伤。

也不需要循环1:length(x)，for语句将迭代地遍历像Ecoregions这样的向量的每个元素，而不会被明确告知这样做。

使用get：

尝试使用此简化示例

ER_10.2 <- data.frame(v1=1:10)
ER_12.1 <- data.frame(v1=2:11)
Ecoregions.list <- c("ER_10.2", "ER_12.1")
Ecoregions <- unique(Ecoregions.list)

for(i in Ecoregions) {
    ER <- get(i)[[1]]
    # add the below line if you want 12 separate plots
    # dev.new() 
    plot(ER)
    # insert other code here which would presumably negate
    # your ability to use the 'lapply' function.
}

我还应该注意，在R中使用get等函数和文本字符串来表示数据名称通常不是一个好主意。通常情况下，如果您有许多相关的data.frames或其他对象，则可以将它们放在list中，如：

Ecoregions <- list(ER_10.2, ER_12.1)

然后，您可以使用lapply或sapply将功能应用于每个组件，如下所示：

lapply(Ecoregions,function(x) plot(x[[1]]))

Answer 2

在这种情况下，您不想使用as.name，但我认为您正在寻找get之类的内容。因此，只需将ER=as.name(Ecoregions[i])替换为ER=get(Ecoregions[i])即可。

但这并不是解决这种情况的好方法。您可以制作list个区域，而不是为每个区域创建新变量。这样，您遍历list元素而不是变量名称。

例如，代替ER_10.2<-Level.2.ecoregs[Level.2.ecoregs$NA_L2CODE=="10.2",]尝试执行以下操作：

# Split your regions by the L2 Code.
list.of.regions<-split(Level.2.ecoregs,Level.2.ecoregs$NA_L2CODE)

现在，如果您的数据发生了变化（例如添加了新的NA_L2CODE），您就不必更改代码，因为它将是您列表中的新元素。

现在，您可以遍历列表中的元素：

# Loop over each element of the list.
for (region in list.of.regions) 
    plot(region)

而且，如果你想得到想象，你可以使用lapply，它只在list的每个元素上运行一个函数。

lapply(list.of.regions,plot)

如果你想让自己变得更加漂亮，并且想要绘制网格中所有区域的情节，可以使用lattice或ggplot。尝试查找ggmap以获取一些不错的示例。

绘制shapefile循环错误

2 个答案: