NDVI パラメータを R で表示し、続いて各領域でこのパラメータの値を抽出したいと思います。
私のファイルは、次のリンクからダウンロードした NetCDF (.nc) 形式です: https://www.ncei.noaa.gov/data/avhrr-land-normalized-difference-vegetation-index/access/1981/
そこで、次のコードの 2 つのパッケージ「raster」と「ncdf4」を使用して、最初のファイルを「AVHRR-Land_v005_AVH13C1_NOAA-07_19810624_c20170610041337.nc」として R にインポートしました。
library(ncdf4)
library(raster)
setwd("path/.../data")
r=raster("AVHRR-Land_v005_AVH13C1_NOAA-07_19810624_c20170610041337.nc",varname="NDVI")
しかし、私のプログラムを適用すると、このエラーが発生します:
Error in CRS(as.character(projection(crs))) :
no arguments in initialization list
以下のコードから変数の名前 (「NDVI」) を選択します。問題は座標系の形式にあると思います!
library(ncdf4)
nc <- nc_open("AVHRR-Land_v005_AVH13C1_NOAA-07_19810624_c20170610041337.nc")
attributes(nc$var)$names
この問題の解決策はありますか!
前もって感謝します
私の R バージョンと私のシステム PC:
R version 3.6.2 (2019-12-12)
Copyright (C) 2019 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
編集
Robert Hijmans さん、回答ありがとうございます (以下の回答)。
パッケージを更新しました (一部のパッケージを更新する際のエラー メッセージ)
update.packages(ask=F)
Warning: package 'BH' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'callr' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'cli' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'DT' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'fansi' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'farver' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'gh' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'mime' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'pillar' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'plyr' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'prettyunits' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'stringi' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'vctrs' in library '/usr/local/lib/R/site-library' will not be updated
Warning: package 'foreign' in library '/usr/lib/R/library' will not be updated
コードを適用しましたが、それでも同じエラーが発生します。
r <- raster(f, varname="NDVI")
Error in CRS(as.character(projection(crs))) :
no arguments in initialization list.
traceback()
No traceback available .
sessionInfo()
R version 3.6.2 (2019-12-12)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 18.04.3 LTS
Matrix products: default
BLAS: /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.7.1
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.7.1
locale:
[1] LC_CTYPE=en_GB.UTF-8 LC_NUMERIC=C LC_TIME=en_GB.UTF-8
[4] LC_COLLATE=en_GB.UTF-8 LC_MONETARY=en_GB.UTF-8 LC_MESSAGES=en_GB.UTF-8
[7] LC_PAPER=en_GB.UTF-8 LC_NAME=C LC_ADDRESS=C
[10] LC_TELEPHONE=C LC_MEASUREMENT=en_GB.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] ncdf4_1.17 raster_3.0-7 sp_1.3-2
loaded via a namespace (and not attached):
[1] compiler_3.6.2 rgdal_1.4-8 tools_3.6.2 yaml_2.2.0 Rcpp_1.0.3 codetools_0.2-16
[7] grid_3.6.2 lattice_0.20-38