代码如下:
GSST<-matrix(1:992160,nrow = 1440,ncol = 689)
lonR<-360/1440
latR<-180/689
for (i in 2:1439){
for (j in 2:688) {
GSST[i,j]<-sqrt(((sst[i+1,j]-sst[i-1,j])/2)^2+((sst[i,j+1]-sst[i,j-1])/2)^2)
lon1[i]<-((i-1)*lonR+lon1[1]+(lonR/2))
lat1[j]<-((j-1)*latR+lat1[1]+(latR/2))
}
}