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

Any question regarding the Interpolation method using Kriging

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

Postby bez » Tue Jul 05, 2016 9:56 am

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é ?
bez
 
Posts: 21
Joined: Tue Oct 23, 2012 10:31 am

Re: Problème de définition de voisinage

Postby Didier Renard » Wed Jul 06, 2016 7:45 pm

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.
Didier Renard
 
Posts: 337
Joined: Thu Sep 20, 2012 4:22 pm

Re: Problème de définition de voisinage

Postby Didier Renard » Sun Jul 17, 2016 6:51 pm

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
Didier Renard
 
Posts: 337
Joined: Thu Sep 20, 2012 4:22 pm


Return to Kriging

Who is online

Users browsing this forum: No registered users and 3 guests

cron