rast() 不会从 NetCDF 文件中提取空间信息：[rast] 未知范围

Question

我一直在尝试使用

rast

中的

terra

函数加载 NetCDF，但我收到以下警告消息：

> SM_rast <- rast("Soil_moisture_v3.nc", "soil_moisture") 
Warning message:
[rast] unknown extent

对应的SpatRaster没有crs，分辨率错误，坐标错误：

> SM_rast
class       : SpatRaster 
dimensions  : 2160, 4320, 172  (nrow, ncol, nlyr)
resolution  : 1, 1  (x, y)
extent      : 0, 4320, 0, 2160  (xmin, xmax, ymin, ymax)
coord. ref. :  
source      : Soil_moisture_v3.nc 
names       : Soil_~_v3_1, Soil_~_v3_2, Soil_~_v3_3, Soil_~_v3_4, Soil_~_v3_5, Soil_~_v3_6, ... 

您可以在这里下载NetCDF文件。

按照this question的答案，我尝试将

drivers="NETCDF"

添加到

rast()

函数，但它甚至无法读取文件：

> SM_rast <- rast("Soil_moisture_v3.nc", "soil_moisture", drivers="NETCDF") 
Error: [rast] cannot open this file as a SpatRaster: Soil_moisture_v3.nc
In addition: Warning message:
`Soil_moisture_v3.nc' not recognized as a supported file format. (GDAL error 4) 

将 SpatRaster 转换为数据帧时，所有包含 NA 的单元格都丢失并且坐标列错误：

> df_rast <- terra::as.data.frame(SM_rast[[1]], xy = TRUE)
> summary(df_rast)
       x                y          Soil_moisture_v3_1
 Min.   :   2.5   Min.   : 414.5   Min.   :0.02000   
 1st Qu.:1525.5   1st Qu.: 937.5   1st Qu.:0.07621   
 Median :2403.5   Median :1240.5   Median :0.17518   
 Mean   :2340.6   Mean   :1197.1   Mean   :0.20558   
 3rd Qu.:2841.5   3rd Qu.:1440.5   3rd Qu.:0.30225   
 Max.   :4317.5   Max.   :1960.5   Max.   :0.91196  

如果我使用

brick()

包中的

raster

做同样的事情，一切正常：

> SM_brick <- brick("data-raw/TROPOMI/Soil_moisture_v3.nc", varname = "soil_moisture")
> SM_brick
class      : RasterBrick 
dimensions : 2160, 4320, 9331200, 172  (nrow, ncol, ncell, nlayers)
resolution : 0.0833333, 0.0833333  (x, y)
extent     : -180, 179.9998, -90, 89.99993  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : Soil_moisture_v3.nc 
names      : X2018.02.02, X2018.02.10, X2018.02.18, X2018.02.26, X2018.03.06, X2018.03.14, X2018.03.22, X2018.03.30, X2018.04.07, X2018.04.15, X2018.04.23, X2018.05.01, X2018.05.09, X2018.05.17, X2018.05.25, ... 
Date       : 2018-02-02, 2021-10-24 (min, max)
varname    : soil_moisture 

> df_brick <- terra::as.data.frame(SM_brick[[1]], xy = TRUE)
> summary(df_brick)
       x                    y              X2018.02.02     
 Min.   :-179.95833   Min.   :-89.95834   Min.   :0        
 1st Qu.: -89.97920   1st Qu.:-44.97919   1st Qu.:0        
 Median :  -0.00008   Median : -0.00004   Median :0        
 Mean   :  -0.00008   Mean   : -0.00004   Mean   :0        
 3rd Qu.:  89.97905   3rd Qu.: 44.97911   3rd Qu.:0        
 Max.   : 179.95818   Max.   : 89.95826   Max.   :1        
                                          NA's   :8181031 
 

我仍然需要使用

terra

，因为它更快，并且具有

raster

中没有的更多有用功能。

Answer 1

我觉得一切都很好：

library(terra)
# terra 1.7.20
rast("Soil_moisture_v3.nc")
#class       : SpatRaster 
#dimensions  : 2160, 4320, 172  (nrow, ncol, nlyr)
#resolution  : 0.0833333, 0.0833333  (x, y)
#extent      : -180, 179.9998, -90, 89.99993  (xmin, xmax, ymin, ymax)
#coord. ref. : lon/lat WGS 84 (EPSG:4326) 
#source      : Soil_moisture_v3.nc 
#varname     : soil_moisture 
#names       : soil_~ure_1, soil_~ure_2, soil_~ure_3, soil_~ure_4, soil_~ure_5, soil_~ure_6, ... 
#unit        :       m3/m3,       m3/m3,       m3/m3,       m3/m3,       m3/m3,       m3/m3, ... 
#time (days) : 2018-02-02 to 2021-10-24 

gdal()
#[1] "3.6.2"

g <- gdal(drivers=T)
g[g$name=="netCDF", ]
#     name   type        can   vsi                  long.name
#88 netCDF raster read/write FALSE Network Common Data Format

如果你没有看到这个，那可能是因为你有另一个版本的 terra 和/或 GDAL，和/或因为你的操作系统以及你如何安装“terra”。

从你的 zip 文件的内容我可以看出你使用的是 Mac。该问题可能与 CRAN 上二进制版本附带的 GDAL 有关。您可以尝试使用 R-Universe 安装

install.packages('terra', repos='https://rspatial.r-universe.dev')

或者，如果这不起作用，请从源代码安装。请参阅说明。

还有这个解决方法：

library(raster)
# add zero to force the values out of the file.
# that takes a while
b <- brick("Soil_moisture_v3.nc") + 0
r <- rast(b)