从NETCDF提取数据

Question

我已经从这里下载了一月的海面温度https://oceancolor.gsfc.nasa.gov/l3/

并将其导入R。

我知道如何使用范围(Ymax，Ymin，Xmax，Xmin(进行裁剪，但是我不知道如何仅对一个站(53.9s，174,1W(或最接近该坐标的站点做。有什么方法可以为一个站裁剪数据吗？

val＆lt; - 提取物(174.1,53.9(.LOCAL(x，y，...(中的错误：扩展不重叠

SST_Jan <- brick("~https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A20021822018212.L3m_MC_SST_sst_9km.nc", stopIfNotEqualSpaced = FALSE, varname = "sst")
print(SST_Jan)
val<-extract(174.1, 53.9)
SST_Jan_station <- extract(SST_Jan, val)

我希望能够在12个月内在该特定位置绘制SST的更改

谢谢，

Answer 1

extract函数与数字向量无法使用。

您可以将坐标放在matrix-

中

pnt = matrix(c(174.1, 53.9), ncol = 2)
pnt
##       [,1] [,2]
## [1,] 174.1 53.9

，然后extract将起作用 -

extract(SST_Jan, pnt)
##      layer
## [1,]  8.24