TP1

Initialisation

## Warning: le package 'fields' a été compilé avec la version R 4.3.3

## Le chargement a nécessité le package : spam

## Warning: le package 'spam' a été compilé avec la version R 4.3.3

## Spam version 2.10-0 (2023-10-23) is loaded.
## Type 'help( Spam)' or 'demo( spam)' for a short introduction 
## and overview of this package.
## Help for individual functions is also obtained by adding the
## suffix '.spam' to the function name, e.g. 'help( chol.spam)'.

## 
## Attachement du package : 'spam'

## Les objets suivants sont masqués depuis 'package:base':
## 
##     backsolve, forwardsolve

## Le chargement a nécessité le package : viridisLite

## 
## Try help(fields) to get started.

Source

source("plot.cor.R")
source("cartosurf2.R")

Chargement des données

load("TEMPERATURE5X5.Rdata")
load("INDEX.Rdata")

Extraire la période commune

TEMP.c=TEMP$data[TEMP$data[,1]>=1950 & TEMP$data[,1]<=1990,]
INDEX.ENSO=INDEX$enso[INDEX$enso[,1]>=1950 & INDEX$enso[,1]<=1990,]
INDEX.NAO=INDEX$nao[INDEX$nao[,1]>=1950 & INDEX$nao[,1]<=1990,]

Extraire les températures d’hiver

id.H=which(TEMP.c[,2]<=3)
TEMP.H=TEMP.c[id.H,]

Extraire les températures d’été

id.E=which(TEMP.c[,2]>=6 & TEMP.c[,2]<=8)
TEMP.E=TEMP.c[id.E,]

Extraire les indices d’hiver et d’été

NAO.H=INDEX.NAO[,2:4]
NAO.E=INDEX.NAO[,7:9]
ENSO.H=INDEX.ENSO[,2:4]
ENSO.E=INDEX.ENSO[,7:9]

Création de matrice vides

TEMP.MOY.H=matrix(NA,nrow=41,ncol=2594)
TEMP.MOY.E=matrix(NA,nrow=41,ncol=2594)

Remplissage des cellules par les moyennes saisonnières

k = 1
i=1950
for(i in 1950:1990){
TEMP.MOY.H[k,]=apply(TEMP.H[which(TEMP.H[,1]==i),],2,mean)
  k=k+1
}
k = 1
i=1950
for(i in 1950:1990){
TEMP.MOY.E[k,]=apply(TEMP.E[which(TEMP.E[,1]==i),],2,mean)
  k=k+1
}

Calcul et stockage pour les indices

NAO.H.MOY=apply(NAO.H,1,mean)
NAO.E.MOY=apply(NAO.E,1,mean)
ENSO.H.MOY=apply(ENSO.H,1,mean)
ENSO.E.MOY=apply(ENSO.E,1,mean)

Calcul des corrélations indices-T

cor.ENSO.T.H=cor(ENSO.H.MOY,TEMP.MOY.H [,3:2594], use="pairwise.complete.obs")
cor.ENSO.T.E=cor(ENSO.E.MOY,TEMP.MOY.E [,3:2594], use="pairwise.complete.obs")
cor.NAO.T.H=cor(NAO.H.MOY,TEMP.MOY.H [,3:2594], use="pairwise.complete.obs")
cor.NAO.T.E=cor(NAO.E.MOY,TEMP.MOY.E [,3:2594], use="pairwise.complete.obs")

Géoréférencer

latitudes=TEMP$lat_long[1,]
longitudes=TEMP$lat_long[2,]

Graphique

plot.cor(longitudes,latitudes,cor.ENSO.T.H,title="Corrélation entre l'indice ENSO et les températures hivernales moyennes")

## [1] "well done"

plot.cor(longitudes,latitudes,cor.ENSO.T.E,title="Corrélation entre l'indice ENSO et les températures estivales moyennes")

## [1] "well done"

plot.cor(longitudes,latitudes,cor.NAO.T.H,title="Corrélation entre l'indice NAO et les températures hivernales moyennes")

## [1] "well done"

plot.cor(longitudes,latitudes,cor.NAO.T.E,title="Corrélation entre l'indice NAO et les températures estivales moyennes")

## [1] "well done"