[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