当 lat 和 lon 在 R 中是双倍时，从 NetCDF 文件中提取“tos”数据

Question

这是我第一次使用堆栈溢出，我的编码技术非常糟糕，我正在使用一个历史tos的NetCDF文件。我想提取特定lat和lon的tos数据。我有一个三维数组中的tos数据，lon和lat分别在一个二维矩阵中。问题是，我选择的lon和lat的行列组合与tos数组的行-列组合不一致。下面是我目前掌握的代码

# Open a NetCDF file
library('ncdf4')
library('raster') 
library('rgdal') 
library('ggplot2')
hist_acc <- "tos_Omon_ACCESS-ESM1-5_historical_r1i1p1f1_gn_185001-201412.nc"
hist_acc1 <- nc_open(hist_acc)
# Get Latitude, Longitude and time from the Netcdf File
lon_1 <- ncvar_get(hist_acc1, "longitude")
lat_1 <- ncvar_get(hist_acc1, "latitude")
tt <- ncvar_get(hist_acc1, "time")
#Extract Temperature of the surface data and substitue fillvalue with NA and close the NetCDF file to leave more space
tos_array <- ncvar_get(hist_acc1, "tos")
fillvalue<- ncatt_get(hist_acc1, "tos", "_FillValue")
tos_array[tos_array == fillvalue$value] <- NA
nc_close(hist_acc1)
# Extract row and column numbers with specific longitudes
wnf_lon <- which(lon_1 > 7 & lon_1 < 8, arr.ind=TRUE) 
rows_wnf_lon <- wnf_lon[,1] 
col_wnf_lon <- wnf_lon[, 2] 
#Extract row and column numbers with specific latitudes
wnf_lat <- which(lat_1 > 61 & lat_1 < 62, arr.ind=TRUE)
rows_wnf_lat <- wnf_lat[,1] 
col_wnf_lat <- wnf_lat[,2] 
# Extract 1 matrix from array (e g day 1)
day1 <- tos_array[, ,1]

我被困在这里，因为我的纬度和纬度矩阵的行和列数与 tos day1 数组的行和列号不对应。

如果你不明白我的问题，请告诉我，我会尽量说得更具体。如果你在R中有NetCDF文件操作的资源，请提前告诉我，谢谢你

Answer 1

访问此类数据的更简单方法是这样的

library(terra)
f <- "tos_Omon_ACCESS-ESM1-5_historical_r1i1p1f1_gn_185001-201412.nc"
r <- rast(f)
r
plot(r, 1)

如果您不熟悉R中的空间（光栅）数据操作，请查看此处