foreach 的神秘错误

Question

在我第一次尝试 foreach 时，我遇到了这个错误：

   Error in function (handle)  :
  no function to return from, jumping to top level

我为每个工作人员得到一个错误。

我正在使用 doMC/多核后端。

起初我认为这是由我的 foreach 调用中的函数试图发送到控制台作为副作用的文本进度条引起的，但在关闭它后我仍然收到错误。尽管存在错误，评估仍在继续，所有仓鼠似乎都在跑步。里面的函数恰好是 reclass() 库（栅格）中的 calc()。谁能猜出这可能意味着什么？我想我会在等待输出是否有效时询问。

---

更新

这是代码的核心：

    foreach( cover= names( classes), .packages= "raster") %dopar% {
      class <- classes[[ cover]]
      calc( mlct$pri,
           function( pri) {
             ifelse( is.na( pri), NA,
                    ifelse( pri ==class, 1, 0))
           },
           filename= paste(
             mlctName,
             paste( cover, ".tif", sep=""),
             sep="_"),
           overwrite= TRUE, ...)
    }

如果不熟悉library( raster)，它提供了一种巧妙的机制，可以与地理空间栅格数据集进行交互，而不必将它们放入内存中。这些函数逐块读取和写入数据。在这种情况下，工作人员都从相同的源数据读取，但写入单独的输出文件。我想知道这个文件 I/O 是否与上面报告的（看似良性的）错误有关。

---

现在它已经吐出了我期望的所有新文件，但这出现在控制台上：

GDAL Error 1: TIFFFetchDirectory:Sanity check on directory count failed, this is probably not a valid IFD offset

尽管核心很安静，但解释器提示符没有重新出现。我猜想，文件 I/O 不知何故出了问题的另一个线索。有什么想法吗？

Answer 1

我想我已经有了答案。我在 GDAL-dev 列表上从 Even Rouault 那里得到了一些帮助。请参阅 GDAL-dev 档案中的此主题 （也可以通过 Nabble 获得）。如果他出现了，他应该得到这个答案的荣誉，而不是我

。问题的关键是我必须将光栅数据集的文件名传递给 foreach 工作人员，并在 foreach 闭包中重新实例化 raster() 对象。每个工作线程都需要有自己的 rgdal 句柄来读取输入数据。

新代码如下所示：

  fracDoparFun <- function( priFilename, ...) {
    foreach( cover= names( classes), .packages= "raster") %dopar% {
      class <- classes[[ cover]]
      frac <-
        calc( raster( priFilename),
             function( pri) {
               ifelse( is.na( pri), NA,
                      ifelse( pri ==class, 1, 0))
             },
             filename= paste(
               mlctName,
               paste( cover, ".tif", sep=""),
               sep="_"),
             overwrite= TRUE, ...)
      return( filename( frac))
    }
  }
mlct$fracs <-
  brick( stack( fracDoparFun( filename( mlct$pri), ...)),
        filename= fracsBrickFile,
        overwrite= TRUE,
        ...)

现在 fracDoparFun() 返回一个文件名列表，这些文件名是工作人员生成的输出，然后由 brick() 将其收集到单个多频段结果中 函数。

感谢大家的关注。我希望我不会过早地关闭它，但我已经对其进行了测试，并且我想捕捉我所学到的东西。