摘要：我有一个包含 NA 值的栅格数据集，并且想要计算它的变异函数，而忽略 NA。我怎样才能做到这一点？
我有一个使用readGDAL函数加载到 R 中的图像，存储为im. 为了使这种可重现，dput图像上的结果可在https://gist.github.com/2780792 获得。我正在尝试显示此数据的变异函数并且正在挣扎。我将经历到目前为止我已经尝试过的：
我尝试了该gstat软件包，但似乎无法获得有效的函数调用。我已经收集到，基本上我需要的是数据值本身 ( im@data$band1) 和坐标 ( coordinates(im))。我尝试了各种命令，例如：
> variogram(locations=coordinates(im), y = im@data$band1)
Error in is.list(object) : 'object' is missing
和
> variogram(coordinates(im), im@data$band1)
Error in variogram.default(coordinates(im), im@data$band1) :
argument object and locations should be lists
我在这里做错了什么？
由于这似乎不起作用，我尝试了geoR使用以下方法调用的软件包：
> variog(coords=coordinates(im), data=im@data$band1)
variog: computing omnidirectional variogram
Error in FUN(X[[1L]], ...) : NA/NaN/Inf in foreign function call (arg 4)
该错误似乎与其中包含 …
我已经通过安装了R和RStudio brew。对于项目，我需要R库中的一个函数，而该函数geoR又依赖于tcl-tk。我也通过brew：安装了它brew install tcl-tk。但是，尝试安装geoR仍然会导致以下错误-我没有得到。重新安装R / RStudio没有帮助。也没有重新启动Mac。哦，尝试geoR从CRAN下载的二进制文件进行安装也无济于事……最后，尝试tcl-tk从R网站安装软件包也不能解决此问题。
我也将其发布在的GitHub上brew-但也许其他人可以访问StackOverflow领域。
我的brew状态
brew config
HOMEBREW_VERSION: 1.5.2
ORIGIN: https://github.com/Homebrew/brew
HEAD: 60a30e966b7cece5bd4823dae3fb981ab85106ea
Last commit: 11 days ago
Core tap ORIGIN: https://github.com/Homebrew/homebrew-core
Core tap HEAD: 8c1c4f86a269cb23d9d92008abf1a37eabb297b6
Core tap last commit: 2 hours ago
HOMEBREW_PREFIX: /usr/local
CPU: quad-core 64-bit skylake
Homebrew Ruby: 2.3.3 => /System/Library/Frameworks/Ruby.framework/Versions/2.3/usr/bin/ruby
Clang: 9.0 build 900
Git: 2.16.1 => /usr/local/bin/git
Curl: 7.54.0 => /usr/bin/curl
Perl: /usr/bin/perl …