Page 1 of 1

[SOLVED] Problème de définition de voisinage

PostPosted: Tue Jul 05, 2016 9:56 am
by bez
Deux problèmes sont documentés ici :
- la différence de résultat quand on utilise neigh.input() et neigh.init()
- la non conformité des résultats avec les paramètres du voisinage

Ci-dessous deux façons de définir le même voisinage.
neigh.1 <- neigh.init(nmini=1,nmaxi=2,dmax=0.5, flag.aniso=T,coeffs=c(1,0.2))
neigh.2 <- neigh.input()
# avec les réponses suivantes 2 1 2 n y n 0.5 2.5
l'édition des objets neigh.1 et neigh.2 confirme que ce sont les mêmes, ie mêmes paramètres affichés.
Cependant le krigeage ne donne pas les mêmes résultats.
On procède au krigeage d'un point par un autre pour tester le voisinage.
On s'attend à ce que le résultat soit une ellipse égale au voisinage avec la valeur connue au point central.
Cependant :
plot(kriging(db.create(x1=0,x2=0,z1=2),db.create(flag.grid = T,x0=c(-1,-1),dx=c(0.01,0.01),nx=c(200,200)),
model=model.create(3),neigh=neigh.1))
plot(kriging(db.create(x1=0,x2=0,z1=2),db.create(flag.grid = T,x0=c(-1,-1),dx=c(0.01,0.01),nx=c(200,200)),
model=model.create(3),neigh=neigh.2))
ne donnent pas la même chose.
Avec neigh.1 le résultat ne présente pas l'anisotropie souhaitée et présente une rotation.
Avec neigh.2 le résultat présente l'anisotropie souhaitée mais avec une rotation

La question est : avec quel voisinage le krigeage est-il réalisé ?

Re: Problème de définition de voisinage

PostPosted: Wed Jul 06, 2016 7:45 pm
by Didier Renard
Nice point.

The problem was linked to the neigh.input() interactive procedure.
This function was asking the minimum and maximum range ... but through the internal conversion used for storage purpose, these values were corrupted.
The neigh.input function has been fixed.

Re: Problème de définition de voisinage

PostPosted: Sun Jul 17, 2016 6:51 pm
by Didier Renard
The new version of neigh.init(), neigh.input(), melem.init() and model.create() are availabel.... only for test
on the private download site (my personal cloud).
When you have downloaded and reasonably tested the new version, please update the title of this post (adding [SOLVED] to its title) so that I can publish this new version.

If the last modifications (dated August 2nd) are satisfactory, please modify the title of this post by adding "[SOLVED]" at its beginning