1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
| #longitude et latitude de la zone d'étude
long<-c(-4.30,-4.29,-4.28,-4.27,-4.26,-4.25,-4.24,-4.23,-4.22,-4.21,-4.20,-4.19,-4.18,-4.17,-4.16,-4.15,-4.14,-4.13,-4.12,-4.11,-4.10,-4.09,-4.08,-4.07,-4.06,-4.05,-4.04,-4.03,-4.02,-4.01,-4.00,-3.99,-3.98,-3.97,-3.96,-3.95,-3.94,-3.93,-3.92,-3.91,-3.90,-3.89,-3.88,-3.87,-3.86,-3.85,-3.84,-3.83,-3.82,-3.81,-3.80,-3.79,-3.78,-3.77,-3.76,-3.75,-3.74,-3.73,-3.72,-3.71,-3.70)
lat<-c(5.20,5.21,5.22,5.23,5.24,5.25,5.26,5.27,5.28,5.29,5.30,5.31,5.32,5.33,5.34,5.35,5.36,5.37,5.38,5.39,5.40,5.41,5.42,5.43,5.44,5.45,5.46,5.47,5.48,5.49,5.50,5.51,5.52,5.53,5.54,5.55,5.56,5.57,5.58,5.59,5.60,5.61,5.62,5.63,5.64,5.65,5.66,5.67,5.68,5.69,5.70,5.71,5.72,5.73,5.74,5.75,5.76,5.77,5.78,5.79,5.80)
#chargement des données pluies
load("c:/Users/HP/Documents/sorties_IDW2/2018-05-03_040000.Rdata")
a<-R_interp_1h
#création d'un data.frame
year <- 2018
month <- 5
day <- 3
hours <- 4
mesure1 <- as.vector(a)
mesure1[mesure1==0]=NaN
lonlat <- as.matrix(expand.grid(long,lat))
organiser <- data.frame(year,month,day,hours,lonlat,mesure1)
str(organiser)
names(organiser)
#renommer les variables
names(organiser)<-c("year","month","day","hours","long","lat","mesure1")
dossier <- dir("c:/Users/HP/Documents/sorties_IDW2",full.names = T)
n <- length(dossier) |
