Bonjour
Je dispose de 5 points de mesure des PM2.5 et 15 points de mesure des PM10, dont 4 points sont en commun :
NOM X Y pm10 pm25
HVH 437.9310 2503.137 44.7 36.2
HRI 437.8720 2501.458 60.5 NA
AIL 500.8050 2547.541 52.2 NA
JUS 509.7810 2494.600 28.2 19.4
POS 518.3270 2479.166 20.4 NA
HON 448.0180 2492.879 40.7 NA
PQV 507.1650 2493.060 27.5 NA
EVT 513.0430 2447.685 22.2 NA
LIS 446.0212 2463.086 23.4 NA
CHV 401.0697 2469.523 49.4 36.7
IFS 403.7104 2464.852 50.4 NA
CHD 313.5191 2522.717 65.2 NA
ALEN 435.3121 2385.290 24.4 NA
STLO 350.1485 2463.035 42.4 NA
MERA 393.9536 2408.769 21.6 14.9
MPB 485.6040 2501.260 NA 16.8
Je dispose aussi d'une sortie de modèle déterministe sur une grille de 3x3 km.
Je m'intéresse aux PM2.5...
J'essaye de faire un co-krigeage co-localisé avec un db qui contient l'ensemble des points avec 2 mesures par points (NA pour l'une des deux sauf 4 points) .
data <- db.create(pm, flag.grid=F, ndim=2, autoname=F)
data <- db.locate(data, "pm25", "z", 1)
data <- db.locate(data, "pm10", "z", 2)
et avec un db qui contient la sortie de modèle que je considère comme une dérive externe :
derive <- db.create(z=modele[,"Z"], nx=c(nx,ny), x0=c(x[1],y[1]), dx=c(3,3), flag.grid=T, autoname=T)
derive <- db.locate(derive, "z", "f", 1)
Voici data :
Data Base Characteristics
=========================
Data Base Summary
-----------------
File is organized as a set of isolated points
Space dimension = 2
Number of fields = 5
Maximum Number of attributes = 5
Total number of samples = 16
Variables
---------
Field = 1 - Name = rank - Locator = NA
Field = 2 - Name = X - Locator = x1
Field = 3 - Name = Y - Locator = x2
Field = 4 - Name = pm10 - Locator = z2
Field = 5 - Name = pm25 - Locator = z1
et voici derive :
Data Base Characteristics
=========================
Data Base Summary
-----------------
File is organized as a regular grid
Space dimension = 2
Number of fields = 4
Maximum Number of attributes = 4
Total number of samples = 6958
Grid characteristics:
Origin : 277.500 2354.500
Mesh : 3.000 3.000
Number : 98 71
Angles : 0.000 0.000
Variables
---------
Field = 1 - Name = rank - Locator = rank
Field = 2 - Name = x1 - Locator = x1
Field = 3 - Name = x2 - Locator = x2
Field = 4 - Name = z - Locator = f1
Ensuite je fais une variographie très succincte et automatique :
vario.exp <- vario.calc(data)
vario.theo = model.auto(vario.exp, struct=c("Cubic", "Spherical"))
et enfin, le krigeage :
carte.colcok <- kriging(dbin=data, dbout=derive, model=vario.theo, neigh=neigh.uniq, uc=c(1,"f1"), radix="", rank.colcok=c(4,4))
Ma question : Que faut-il indiquer dans le paramètre rank.colcok ?
Merci d'avance
Michel