我正在尝试使用R文件将历史温度数据读取到.nc中。我无法从所需的坐标中提取温度数据。
我的代码：
sst <- brick(x='~/Desktop/HadISST_sst.nc', values=T, varname='sst')
nc_hadisst = nc_open("HadISST_sst.nc")
lat = ncvar_get(nc_hadisst,"latitude")
lon = ncvar_get(nc_hadisst,"longitude")
time = ncvar_get(nc_hadisst,"time")
sst = aperm(sst,c(2,1,3))
sst[which(sst< (-99), arr.ind=T)] = NA
mesh <- meshgrid(lon,lat)
LON <- mesh$X
LAT <- mesh$Y
LATvec <- as.vector(LAT)
LONvec <- as.vector(LON)
latlon <- cbind(LATvec,LONvec)
findPoint <- cbind(22.302, 38.99)
nearestPoint = apply(gDistance(SpatialPoints(latlon), SpatialPoints(findPoint), byid=TRUE), 1, which.min)
latPoint <- c(LATvec[nearestPoint])
lonPoint <- c(LONvec[nearestPoint])
coords <- which(LAT==latPoint & LON==lonPoint, arr.ind = T)
运行最后一行代码后，我不断收到此错误。
Error in arr.ind && !is.null(d <- dim(x)) : invalid 'x' type in 'x && y'
发生了什么事！我不明白这怎么了...