1

我正在尝试绘制一个包含 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

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

提前致谢！

4

2 回答 2

3

我不确定循环是否是这种情况下的最佳实践，但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]]))
于 2013-07-03T03:04:12.783 回答
3

在这种情况下，您不想使用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，它只是在 a 的每个元素上运行一个函数list。

lapply(list.of.regions,plot)

但是，如果您想变得更漂亮，并且想在网格中绘制所有区域的图，则可以使用latticeor ggplot。尝试寻找ggmap一些很好的例子。

于 2013-07-03T03:08:43.290 回答