Как прочитать файл NetCDF "NDVI NOAA AVHRR" как растр в R? - PullRequest
Новая
       125

Как прочитать файл NetCDF "NDVI NOAA AVHRR" как растр в R?

1 голос
/ 23 января 2020

Я хотел бы отобразить параметр NDVI с помощью R и затем извлечь значения этого параметра в каждом регионе.

Мой файл в формате NetCDF (.n c), загруженный по следующей ссылке: https://www.ncei.noaa.gov/data/avhrr-land-normalized-difference-vegetation-index/access/1981/

Итак, я импортировал первый файл как «AVHRR-Land_v005_AVH13C1_NOAA-07_19810624_c20170610041337.n c» под R, используя два пакета «растр» и «ncdf4» из следующих код: 

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 и моя система P C: 

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)

РЕДАКТИРОВАТЬ

Спасибо, Роберт Хейманс, за ваш ответ (ответ ниже).

Я обновил пакеты (сообщения об ошибках при обновлении некоторых пакетов) 

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

1 Ответ

0 голосов
/ 24 января 2020

Скачать файл 

url <- "https://www.ncei.noaa.gov/data/avhrr-land-normalized-difference-vegetation-index/access/1981/"
f <- "AVHRR-Land_v005_AVH13C1_NOAA-07_19810624_c20170610041337.nc"
if (!file.exists(f)) download.file(file.path(url, f), f, mode="wb")

С R 3.6.2 и обновленными пакетами я получаю: 

library(raster)
#Loading required package: sp
r <- raster(f, varname="NDVI")
#Loading required namespace: ncdf4
r
#class      : RasterLayer 
#dimensions : 3600, 7200, 25920000  (nrow, ncol, ncell)
#resolution : 0.05, 0.05  (x, y)
#extent     : -180, 180, -89.99999, 90  (xmin, xmax, ymin, ymax)
#crs        : +init=EPSG:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
#source     : AVHRR-Land_v005_AVH13C1_NOAA-07_19810624_c20170610041337.nc 
#names      : NOAA.Climate.Data.Record.of.Normalized.Difference.Vegetation.Index 
#z-value    : 1981-06-24 
#zvar       : NDVI

Можете ли вы запустить update.packages(ask=F) и попробовать еще раз? Если по-прежнему происходит сбой, вы можете показать traceback() после возникновения ошибки, а также свой sessionInfo()?

...