次の質問に対する回答を複製しようとしていますが、R Terra area コマンドでエラーが発生しています。
平射極図法を使用して衛星データからグリッド セル領域を推定する
私が実行している正確なコードは次のとおりです。
url1 <- "https://polarwatch.noaa.gov/erddap/griddap/nsidcCDRiceSQnhmday.nc?seaice_conc_monthly_cdr[(2019-12-16T00:00:00Z):1:(2019-12-16T00:00:00Z)][(5837500.0):1:(-5337500.0)][(-3837500.0):1:(3737500.0)]"
f <- 'nsidcCDRiceSQnhmday_935c_47bd_a147.nc'
download.file(url1, destfile=f, mode="wb")
library(terra)
r <- rast(f)
crs(r) <- "epsg:3976"
res(r)
a <- area(r, correct=TRUE, mask=TRUE, sum=FALSE)
印刷 r 収量
class : SpatRaster
dimensions : 448, 304, 1 (nrow, ncol, nlyr)
resolution : 25000, 25000 (x, y)
extent : -3850000, 3750000, -5350000, 5850000 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=stere +lat_0=-90 +lat_ts=-70 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84
+units=m +no_defs
source : nsidcCDRiceSQnhmday_935c_47bd_a147.nc
varname : seaice_conc_monthly_cdr (NOAA/NSIDC Climate Data Record of Passive
Microwave Monthly Northern Hemisphere Sea Ice Concentration)
name : seaice_conc_monthly_cdr
unit : 1
time : 2019-12-16
area() コマンドの行に次のエラーが表示されます。
Error: [area] area(x, sum=FALSE) will be removed. Use "cellSize(x)"
デフォルトの R インスタンスでこれを試し、クリーンなインスタンスを最初からセットアップしました。両方のインストールで同じエラーが発生します。どちらのシステムも macOS です。area コマンドを cellSize コマンドに変更しようとすると、R が完全に終了して再起動します。
各セットアップの sessionInfo() は次のとおりです。
システム 2: MacOS Big Sur のクリーン インストール
R version 4.1.0 (2021-05-18)
Platform: x86_64-apple-darwin17.0 (64-bit)
Running under: macOS Big Sur 11.2.3
Matrix products: default
LAPACK: /Library/Frameworks/R.framework/Versions/4.1/Resources/lib/libRlapack.dylib
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] raster_3.4-12 sp_1.4-5 Rcpp_1.0.7 ncdf4_1.17 terra_1.3-13
loaded via a namespace (and not attached):
[1] lattice_0.20-44 codetools_0.2-18 digest_0.6.27 grid_4.1.0
[5] evaluate_0.14 rlang_0.4.11 rmarkdown_2.9 tools_4.1.0
[9] xfun_0.24 yaml_2.2.1 compiler_4.1.0 htmltools_0.5.1.1
[13] knitr_1.33
システム 1: macOS Catalina
R version 4.0.2 (2020-06-22)
Platform: x86_64-apple-darwin17.0 (64-bit)
Running under: macOS Catalina 10.15.7
Matrix products: default
BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/4.0/Resources/lib/libRlapack.dylib
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] terra_1.2-10 raster_3.4-5 sp_1.4-4 downloader_0.4 XML_3.99-0.5
[6] RCurl_1.98-1.2
loaded via a namespace (and not attached):
[1] Rcpp_1.0.5 rerddap_0.7.4 pillar_1.4.7 compiler_4.0.2
[5] bitops_1.0-6 tools_4.0.2 ncdf4_1.17 digest_0.6.27
[9] lattice_0.20-41 jsonlite_1.7.1 evaluate_0.14 lifecycle_0.2.0
[13] tibble_3.0.4 gtable_0.3.0 pkgconfig_2.0.3 rlang_0.4.8
[17] rstudioapi_0.13 crul_1.1.0 curl_4.3 yaml_2.2.1
[21] xfun_0.19 dplyr_1.0.2 httr_1.4.2 knitr_1.30
[25] xml2_1.3.2 generics_0.1.0 vctrs_0.3.5 rappdirs_0.3.3
[29] grid_4.0.2 tidyselect_1.1.0 glue_1.4.2 httpcode_0.3.0
[33] data.table_1.13.2 R6_2.5.0 rmarkdown_2.5 purrr_0.3.4
[37] ggplot2_3.3.2 hoardr_0.5.2 magrittr_1.5 codetools_0.2-16
[41] scales_1.1.1 ellipsis_0.3.1 htmltools_0.5.0 colorspace_2.0-0
[45] tinytex_0.27 munsell_0.5.0 crayon_1.3.4