Kriging

Z Wikipedie, otevřené encyklopedie
Skočit na navigaci Skočit na vyhledávání

Kriging či počeštěně krigování je ve statistice (původně v geostatistice) metoda interpolace, kde jsou interpolované hodnoty modelovány gaussovským procesem podle apriorních kovariancí. Za vhodných předpokladů dává kriging nejlepší lineární nestrannou předpověď střední hodnoty. Interpolační metody založené na jiných kritériích, jako je například hladkost, nemusejí přinést nejpravděpodobnější střední hodnoty. Tato metoda se běžně používá v oblasti prostorové analýzy a počítačových experimentů. Technika je také známá jako Kolmogorova-Wienerova predikce.

Teoretický základ pro metodu byl vyvinut francouzským matematikem Georgem Matheronem na základě diplomové práce Daniela G. Kriga. Průkopníka plotruváženého průměru zlatých nalezišť ve Westminském na útesovém komplexu v jižní Africe. Krige se snažil odhadnout nejpravděpodobnější rozložení zlata na základě vzorků z několika vrtů.

Hlavní zásady[editovat | editovat zdroj]

Související pojmy a techniky[editovat | editovat zdroj]

Jednoduchá myšlenka je předpoklad funkční hodnoty daného bodu vypočítaného váženého průměru známých funkčních hodnot sousedních bodů. Metoda je matematicky úzce související s regresní analýzou. Obě teorie odvozují nejlepší lineární nestranný odhad, založený na předpokladu kovariance. Využitím gauss-Markovy teorie dokážeme nezávislost odhadů a omylů a využívá velmi podobné vzorce. Nicméně jsou použitelné rozdílné rámce: kriging je dělaný pro odhad jednu realizaci z náhodných polí, zatímco regresní model je založen na základě pozorování vícerozměrné datové sady.

Odhad krigingu může být také viděn jako křivka v reprodukujícím se jádra Hilbertova prostoru s reprodukujícím jádrem dané kovarianční funkce. Na rozdíl od klasického krigingu je provedena interpretace: zatímco křivka je motivována minimální normou interpolací založenou na Hilbertově prostorové struktuře, kriging je motivován očekávanou kvadratickou chybou na základě stochastického modelu.

Kriging s polynomiálním trendem ploch je matematicky identický jako zobecněný polygon nejmenších čtverců křivky.

Kriging může být známý jako forma Bayesianova závěru. Kriging začíná s základním rozšířením přes funkce. Tento základ bere formu gausova procesu. N vzorků z funkce která má normální rozdělení a kovariance mezi dvěma vzorky je kovariance funkce (nebo jádra) gausova procesu hodnoceny prostorovým rozmístěním dvou bodů. Sada hodnot je pak pozorována, každá hodnota souvisí s prostorovým rozmístěním. Takže nová hodnota může být předpovězena v nějakém novém prostorovém rozmístění, kombinací základní gausovskou funkcí s gausovskou věrohodnostní funkcí pro každou pozorovanou hodnotu. Výsledná distribuce pozadí je také Gaussian s průměrnou hodnotou a kovariancí, která může být snadno vypočítána z pozorovaných hodnot, jejich rozdílnosti a jádra matic odvozených od předchozích.


Geostatistický odhad[editovat | editovat zdroj]

V geostatistickém modelu jsou některá data interpretována jako výsledky náhodného procesu. Skutečnost, že tyto modely obsahují nejistotu v jeho konceptualizaci, neznamená, že tento jev – les, zavodnění nebo minerální usazenina – mají výsledky z náhodného procesu ale pouze umožňují vytvořit metodologické základy pro prostorové závěry z množství nevypozorovaných lokalit a kvantifikovat nejistoty spojené s odhadem. Stochastický proces je v rámci tohoto modelu, způsob jak přistupovat k souboru údajů shromážděných ze vzorků. První krok v geostatistických modulacích je vytvoření náhodného procesu, který nejlépe popisuje sadu experimentálních pozorovaných dat. Hodnota prostorově umístěná v (obecné označení datové sady s geografickým souřadnicovým systémem) je interpretovaná jako realizace náhodné proměnné . V prostoru , kde datová sada vzorků je rozptýlená, existuje realizace těchto náhodných veličin, korelujících mezi sebou. Datová sada náhodných veličin, představuje náhodnou funkci pouze u těch u kterých je známá pouze jedna realizace - soubor experimentálních dat. S pouze jednou realizací každé náhodné proměnné, je teoreticky možné určit všechny statistické parametry jednotlivých proměnných nebo funkcí.

Navrhované řešení geostatistických formalismů spočívá v předpokladu různých stupňů stacionarity v náhodné funkci, aby bylo možné se domnívat závěru některých statistických hodnot.

Například, jestliže skupina vědců předpokládá odpovídající, je vhodné na základě homogenity vzorků v oblasti , kde je distribuována proměnná, hypotetizovat, že první část je stacionární (všechny náhodné proměnné mají stejný průměr), z čehož vyplývá, že průměr může být odhadnut aritmetickým průměrem vybraných hodnot. Soudě hypotézu jako je tato, tak v případě potřeby dostatečně ověřte homogenitu vzorových hodnot, jestli jsou reprezentativní.

Hypotéza stacionarity týkající se druhé části je definována následujícím způsobem: korelace mezi dvěma náhodnými proměnnými závisí pouze na prostorové vzdálenosti, která je odděluje a je závislá na jeho poloze:



kde


Tato hypotéza umožňuje odvozovat kroky - variogram a kovariogram - na základě vzorků:

kde


Lineární odhad[editovat | editovat zdroj]

Prostorový závěr nebo odhad kvantity na neurčeném místě je vypočítán z lineární kombinace pozorovaných hodnot a vah :


Váhy jsou určeny ke shrnutí dvou extrémně důležitých procedur v závěru prostorového procesu:

  • odrážet strukturální "blízkost" vzorků odhadovaného umístění,
  • ve stejnou dobu by měli mít degradující efekt, aby se předešlo zkreslení případnými shluky

Když vypočítáváme váhy , tak jsou v geostatistice dva závěry: nestranný a minimální odchylka odhadu.

Jestliže mračno reálných hodnotje vykresleno proti odhadovaným hodnotám tak kritérium pro souhrnnou nestrannost, vnitřní stacionarita nebo široký smysl stacionarity oboru, znamená, že průměr odhadů musí být roven průměru skutečných hodnot.

Druhé kritérium říká, že průměrná kvadratická odchylka musí být minimální, což znamená, že když se mračno odhadovaných hodnot oproti mraku skutečných hodnot více rozptýlí, odhad je více nepřesný.


Metody[editovat | editovat zdroj]

V závislosti na stochastických vlastnostech náhodného pole a různým stupněm předpokládané stacionarity, můžou vzniknout různé metody pro výpočet hmotnosti. Klasické metody jsou:

  • Běžný kriging předpokládá stacionaritu prvního momentu všech náhodných veličin: , kde je neznámé.
  • Jednoduchý kriging předpokládá známý stacionární průměr: , kde je známé.
  • Univerzální kriging počítá s obecným polynomickým modelem trendu, např. lineární model trendu .
  • IRFk-kriging předpokládá, že je polynomická neznámá v .
  • Indikační kriging používá indikační funkce namísto samotného procesu, aby bylo možné odhadnout pravděpodobnosti přechodů.
  • Více-indikační kriging je verze indikačního krigingu pracující s rodinou ukazatelů. Nicméně, V-IK přestal v posledních letech vyhovovat jako interpolační technika. K tomu došlo hlavně díky problémům spojeným s provozem a validací modelu. Podmíněná simulace se v tomto oboru rychle stává uznávanou náhradní technikou.
  • Disjunktivní kriging je nelineární zobecnění krigingu.
  • Logaritmicko-normální kriging interpoluje pozitivní data pomocí logaritmů.


Běžný kriging[editovat | editovat zdroj]

Neznámá hodnota je interpretována jako náhodná proměnná nacházející se v , jakož i hodnoty sousedních vzorků . Odhad je také interpretován jako náhodná proměnná nacházející se v , výsledek lineární kombinace proměnných.

Aby bylo možné vyvodit krigovací systém pro předpoklad modelu, následující chyba při odhadování v je deklarována:

Obě kritéria kvality uvedená dříve mohou nyní být vyjádřena z hlediska střední hodnoty a rozptylu nové náhodné proměnné :

Objektivnost

Vzhledem k tomu, že náhodná funkce je stacionární , uvádí se následující omezení:

Aby se zajistilo, že model je objektivní, součet vah musí být jedna.

Minimální variance: minimalizuje

Dva odhady můžou mít , ale disperze kolem jejich střední hodnoty určí rozdíl mezi kvalitou odhadů.

* navštivte covariance matrix pro detailní vysvětlení (v angličtině)

* kde literály představují stand for .

Jednou definovaná kovariance modelu nebo variogram, nebo , platné ve všech oblastech analýzy , pak můžeme definovat vztah pro odhad rozptylu libovolného odhadu v závislosti na kovarianci mezi vzorky:

Některé závěry mohou být z těchto výrazů stanoveny. Rozptyl odhadu:

  • není kvantifikovatelný na kterýkoliv lineární odhad, jakmile se předpokládá stacionarita průměru a prostorových kovariancí nebo variogramy.
  • roste, když se kovariance mezi vzorky a místem pro odhad snižuje. To znamená, že když vzorky jsou daleko od , nejhorší je odhad.
  • roste zároveň s "prioritním" rozptylem proměnné . Pokud je proměnná méně rozptýlená, rozptyl je nižší v každém bodě prostoru .
  • nezávisí na hodnotách vzorků. To znamená, že stejné prostorové uspořádání (se stejnými geometrickými vztahy mezi vzorky a body odhadu) vždy reprodukuje stejný odhad rozptylu v jakékoli části areálu . Tímto způsobem rozptyl neměří nejistotu odhadu produkovanou místní proměnnou.


Systém Kriging[editovat | editovat zdroj]


Řešení tohoto optimalizačního problému má za následek vznik systému kriging:

další parametr je Lagrangeův násobič použitý při minimalizaci chyby krigingu pro dosažení nestrannosti.

Jednoduchý kriging[editovat | editovat zdroj]

Jednoduchý kriging je matematicky nejjednodušší, ale nejméně obecný. To předpokládá, že očekávaná hodnota náhodného pole je známá, a opírá se o kovarianční funkci. Nicméně, ve většině aplikací nejsou očekávání, ani kovariance známy předem.

Praktické předpoklady pro uplatnění jednoduchého krigingu jsou:

  • Široká stacionarita pole.
  • Očekávání je všude nulové: .
  • Známá kovarianční funkce


Systém Kriging[editovat | editovat zdroj]

Krigovací váhy jednoduchého krigingu nemají nestrannou podmínku a jsou dány systémem rovnic jednoduchého krigingu:

Toto je analogické s lineární regresí na jiném .

Odhad[editovat | editovat zdroj]

Interpolace pomocí jednoduchého krigingu je dána:

Krigovací chyba ja dána:

což vede ke generalizované verzi nejmenších čtverců Gauss-Markovova teorému:

Vlastnosti[editovat | editovat zdroj]

  • Odhad krigingu je nestranný::
  • Odhad krigingu využívá skutečně pozorované hodnoty: (za předpokladu, že nevznikla chyba ve měření, při kterém vznikly)
  • Odhad krigingu je nejlepší lineární nestranný odhad pro v případě, že předpoklad trvá. Avšak:
    • Stejně jako u jiných metod: Pokud předpoklady netrvají, krigování může být špatné.
    • Mohou existovat lepší nelineární a/nebo zkreslené metody.
    • Vlastnosti nejsou garantovány, pokud je použit nesprávný variogram. Nicméně obvykle je dosaženo ještě "dobré" interpolace.
    • Nejlepší nemusí být nutně dobré: např. V případě neprostorové závislosti je krigovací interpolace jen tak dobrá jako aritmetický průměr.
  • Kriging používá jako měřítko přesnosti. Nicméně toto opatření se opírá o správnost variogramu.


Aplikace[editovat | editovat zdroj]

Ačkoli kriging byl původně vyvinut pro použití v geostatistice, je to obecná metoda statistické interpolace, která může být použita v jakékoliv disciplíně na vzorku dat z náhodných polí, které splňují příslušné matematické předpoklady.

K dnešnímu dni se kriging používá v celé řadě oborů, včetně následujících:

  • vědy o životním prostředí
  • hydrogeologie
  • těžba
  • přírodní zdroje
  • dálkový průzkum
  • hodnocení realit

a mnoho dalších.


Softwary, které využívají kriging[editovat | editovat zdroj]

  • R packages
  1. BACCO - Bayesian analysis of computer code software
  2. tgp - Treed Gaussian processes
  3. DiceDesign, DiceEval, DiceKriging, DiceOptim - metamodeling packages of the Dice Consortium.
  • Matlab/GNU Octave
  1. mGstat - Geostistics toolbox for Matlab.
  2. DACE - Design and Analysis of Computer Experiments. A matlab kriging toolbox.
  3. GPML - Gaussian Processes for Machine Learning.
  4. STK - Small (Matlab/GNU Octave) Toolbox for Kriging for design and analysis of computer experiments.
  5. scalaGAUSS - Matlab kriging toolbox with a focus on large datasets
  • Scilab
  1. DACE-Scilab - Scilab port of the DACE kriging matlab toolbox
  2. krigeage - Kriging toolbox for Scilab
  3. KRISP - Kriging based regression and optimization package for Scilab
  • Python
  1. scikit-learn - machine learning in Python

Reference[editovat | editovat zdroj]

V tomto článku byl použit překlad textu z článku Kriging na anglické Wikipedii.