[Gfoss] Help interpolazione spaziale con R gstat

Amedeo Fadini fame a libero.it
Gio 30 Ott 2014 14:12:57 CET


Ok, qualche progresso l'ho fatto... allego i dati con cui sto lavorando.

Il mio variogramma è questo:

> segnale.vgm
   np      dist     gamma dir.hor dir.ver   id
1   1  347.8193 0.6517344       0       0 var1
2   3  456.4563 0.1712414       0       0 var1
3   6  607.5393 0.2174966       0       0 var1
4   5  830.0280 0.3035470       0       0 var1
5  15  981.8709 0.7539914       0       0 var1
6  15 1129.0870 1.4327383       0       0 var1
7  21 1316.4489 1.2885020       0       0 var1
8   8 1485.5518 6.6526687       0       0 var1
9  16 1700.7135 2.8137764       0       0 var1
10 21 1825.2846 3.0281450       0       0 var1
11 18 2028.3255 3.6824488       0       0 var1
12 22 2189.8199 4.8961914       0       0 var1
13 15 2369.3774 2.2049531       0       0 var1
14 15 2567.6370 2.4152567       0       0 var1

e dopo qualche tentativo credo di averlo interpolato in maniera decente con

> segnale.fit <- fit.variogram(segnale.vgm, vgm(1, "Wav", 2000))

Le coordinate sono metriche e ho pensato che il range di 2Km,
trattandosi di segnale 3G è più che sufficiente.

Il primo argomento di vgm() è il nugget o il sill? Ho messo 1 solo
perché funziona..

Commenti?

Un grazie a tutti

amefad
-------------- parte successiva --------------
Un allegato non testuale è stato rimosso....
Nome:        amefad.RData
Tipo:        application/octet-stream
Dimensione:  1730 bytes
Descrizione: non disponibile
URL:         <http://lists.gfoss.it/pipermail/gfoss/attachments/20141030/5bcdbff2/attachment.obj>


Maggiori informazioni sulla lista Gfoss