Re: [FreeGeoCZ] Kriging interpolace veliciny v GPS souradnicích

[Martin Landa <landa.martin AT gmail.com>]

Zdravim,

2009/12/20 Jan Tošovský <j.tosovsky AT tiscali.cz>:

> Jde o velicinu merenou za soucasneho zaznamenani pozice. Vysledkem by byl
> patrne CSV soubor se tremi hodnotami na radek - dve GPS souradnice a
> namerena hodnota.
>
> Tim ze je mereni zatizene chybou i vzhledem k charakteru veliciny se
> osvedcilo vysledna data prohnat Kriging interpolaci s Nugget efektem. To se
> da s hodnotami merenymi v pevnem rastru (ne GPS) s vyuzitim sw typu Surfer.
> Pri hledani vhodneho nastroje na nepravidelne usporadana data jsem se dostal
> az ke Grassu, konkretne k wiki http://grass.osgeo.org/wiki/Kriging, kde se
> hovori o nejakem rozsireni pro nastroj R. Puvodne jsem vubec nehledal GIS
> nastroj, ale po letmem seznameni vidim, ze je to patrne presne to, co
> hledam.

pravda, GRASS nema vlastni implementaci pro kriging, nicmene nabizi
nekolik (add-ons) nastroju, ktere to resi pres Rko [1,2]. V
soucasnosti je asi nejucelenejsi v.krige [3]. v.krige byl vyvijen v
ramci Google SoC 2009 a pred nekolika dny byl zaclenen do oficialni
distribuce GRASSu 7.0 [4] a GRASS 6.5 [5]. Jde o Python skript, ktery
by mel fungovat i v GRASS 6.4 (krome GUI).

> Rad bych se zeptal, zda je mozne ma data nejak Grassu podstrcit (aby vedel,
> ze jde o GPS souradnice a namerenou hodnotu) a nasledne provest tu
> interpolaci? Bonusem by bylo nasledne vygenerovani iso kontur a export do

v.in.ascii [6].

> nejakeho vektoroveho formatu, treba SVG, kde by x,y souradnice jednotlivych

v.out.svg [7]

> bodu iso konturove polyline byly GPS souradnice. Pokud Grass nesvede

r.contour [8]

> vektorovy vystup iso kontur, potreboval bych nejak vyexportovat surova
> interpolovana data pro dalsi zpracovani.
>
> Cele bych to rad pouzil na Windows 7 64-bit, coz je priznavam trochu
> nestandardni kombinace.

To bude mirne problem, v soucasnosti se generuji binarky pouze GRASS
6.4 [9]. Proces kompilace je ale pomerne detailne popsan [10].

Martin

[1] http://grass.osgeo.org/wiki/GRASS_AddOns#v.autokrige
[2] 
http://grass.osgeo.org/wiki/GRASS_AddOns#v.surf.krige_.5Bdeprecated:_use_v.autokrige_instead.5D
[3] http://grass.osgeo.org/wiki/V.krige_GSoC_2009
[4] http://trac.osgeo.org/grass/browser/grass/trunk/scripts/v.krige
[5] 
http://trac.osgeo.org/grass/browser/grass/branches/develbranch_6/scripts/v.krige
[6] http://grass.osgeo.org/grass64/manuals/html64_user/v.in.ascii.html
[7] http://grass.osgeo.org/grass64/manuals/html64_user/v.out.svg.html
[8] http://grass.osgeo.org/grass64/manuals/html64_user/r.contour.html
[9] http://grass.osgeo.org/grass64/binary/mswindows/
[10] http://trac.osgeo.org/grass/wiki/CompileOnWindows

-- 
Martin Landa <landa.martin gmail.com> * http://gama.fsv.cvut.cz/~landa