Question

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

并将其导入R。

我知道如何使用范围（ymax，ymin，xmax，xmin）进行裁剪，但是我无法弄清楚如何仅对一个测站（53.9S，174,1W）或与该坐标最接近的测站进行裁剪。有没有一种方法可以只为一个站裁剪数据？

val <-提取（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

从NetCDF提取数据

1 个答案: