Вычислить корреляцию между температурой поверхности моря и количеством осадков, получена ошибка в R - PullRequest
Новая
       70

Вычислить корреляцию между температурой поверхности моря и количеством осадков, получена ошибка в R

0 голосов
/ 18 июня 2020

Я ищу пространственную корреляцию между температурой поверхности моря и количеством осадков. Осадки, которые я использую, представляют собой осадки для наблюдения за одной точкой, и я хочу выяснить, какие области выпадения осадков сильно коррелируют с температурой поверхности моря

Я ищу корреляции между количеством осадков и температурой поверхности моря, используя следующий код: 

# load libraries ----------------------------------------------------------

library(chron)
library(RColorBrewer)
library(lattice)
library(ncdf4)
library(maptools)
library(fields)
library(readxl)
library(raster)
library(sp)
library(rgdal)

# import data -------------------------------------------------------------

# set directory
alamat_file <- "E:/project/paper/Prediksi CH dengan SST/data/sst/"
alamat_simpan <- "E:/project/paper/Prediksi CH dengan SST/proses/"
setwd(alamat_file)

sst.nc <- nc_open("simultan.nc")
rainfall.xlsx <- read_excel("ch_das_palu_82_18.xlsx")
data(wrld_simpl)


# processing --------------------------------------------------------------

# get variable
rainfall <- matrix(as.numeric(rainfall.xlsx$ch))
sst <- ncvar_get(sst.nc, varid = "sst")
sst <- sst-273.15

lon_sst <- ncvar_get(sst.nc, varid = "longitude")
lat_sst <- ncvar_get(sst.nc, varid = "latitude")

# calculate correlation between sst and rainfall
corr_pcp_sst<-matrix(0, dim(sst)[1], dim(sst)[2]) # created a container

for (i in 1:dim(sst)[1]){
  for (j in 1:dim(sst)[2]){
    corr_pcp_sst[i,j]<-cor(sst[i,j,],rainfall)
  }
}

mypalette<- brewer.pal(11, "RdYlBu")
image.plot(lon_sst,rev(lat_sst), corr_pcp_sst, col = mypalette,
           xlab ="Longitude", ylab="Latitude",
           main ="Correlation between Rainfall and SST")

plot(wrld_simpl,add=TRUE)

, но я получил пустую корреляцию вроде этой

enter image description here

как решить эту проблему, чтобы появилось число корреляции? спасибо

Я использую следующие данные:

данные об осадках:

здесь

данные о температуре поверхности моря:

здесь

...