Software Pachetul GLOBEC Kriging – EasyKrig2.1, 1 mai 2000

Original: http://globec.whoi.edu/software/kriging/easy_krig/easy_krigV21.html

Copyright (c) 1998 proprietatea Dezhang Chu și Institutul Oceanografic Woods Hole. �Toate drepturile rezervate.

Software-ul kriging descris în acest document a fost elaborat de Dezhang Chu cu finantare de la National Science Foundation prin Serviciul proiect ne GLOBEC Georges Bank a programului și de birou de management al datelor. Acesta a fost inspirat de un set de instrumente MATLAB dezvoltat de Yves Gratton și Caroline Lafleur (INRS-Ocanologie, Rimouski, Qc, Canada), și Jeff Runge (Institut Maurice-Lamontagne). Acest software poate fi reprodus pentru necomercial

scopuri numai.

Acest program este distribuit în speranța că va fi util, dar FĂRĂ NICI O GARANȚIE. Contactați Dr. Chu la [email protected] cu accesorii sau sugestii pentru schimbări.

Cuprins:

1. Introducere

1.1 �Informații generale

1.1.1 Sisteme de operare

1.1.2 descrieri scurte

1.2 �Noțiuni de bază

1.2.1 Jos-încărcare a programului

1.2.2 de inițiere rapidă

2. Etapele de prelucrare a datelor

2.1 �DatePregătire

2.2 �Semi-variogramelor

2.3 �Kriging

2.4 �Vizualizare

2.5 �Salvarea rezultatelor

3. Exemple

4. Referințe

1. Introducere

1.1 Informații generale

1.1.1 Sisteme de operare

Software-ul a fost inițial dezvoltat în Matlab 5.2 pe un PC și destinat a fi calculator și / sau sistem de operare independent de platforma. Programul a fost testat pe mai multe masini si efectueaza bine. [Versiunea curentă a software-ului easyKrig nu este compatibil cu Matlab 6. O versiune care este compatibilă va fi disponibilă în curând. 10/23/03] Următorul tabel arată masini testate / Sistem de operare:

Computer ������������������������������OS

————————————————– ——————————————–

Windows 95 compatibil��������������������� PC / NT, Windows 98 *

PC ���������������������Linux (Red Hat 5.2)

SUN ���������������������UNIX (Solaris)

SGI ���������������������UNIX (IRIX)

Macintosh ���������������������OS 8.6

————————————————– ——————————————–

* Windows98 nevoie pentru a obține patch-uri de la Mathworks de a stabili interfața grafică (GUI) bug-uri legate, vezi instrucțiunile din Mathworks site-ul web:

http://www.mathworks.com/

1.1.2 Scurtă descriere

Programul utilizează o interfață grafică cu utilizatorul (GUI) pentru a simplifica modul de funcționare. Este nevoie de Matlab 5.2 sau mai mare cu sau fără optimizare set de instrumente (a se vedea secțiunea 2.2) și este format din cinci componente, sau etape de prelucrare: (1) pregătire de date, (2) calcul variogramelor, (3) kriging, (4) vizualizare și (5 ) rezultate de economisire. Acesta permite utilizatorului să prelucreze datele anizotrope, selectați o listă de modele variogramelor, și o gamă de metode kriging, care sunt, de asemenea, caracteristici comune ale alte pachete software existente. Unul dintre avantajele majore ale acestui pachet de programe este că programul minimizează cerințele utilizatorilor de a “ghici” parametrii inițiali și automat generează parametrii impliciți necesare. În plus, deoarece foloseste o interfață grafică, modificările din setările inițiale ale parametrilor pot fi efectuate cu ușurință. O altă caracteristică a acestui pachet de programe este faptul că are o bibliotecă built-in on-line de ajutor care permite utilizatorului să înțeleagă utilizarea butoanelor de parametri și de operare cu doar un simplu click.

Deși acest pachet software lipsesc câteva abilități, cum ar fi 3D si Co-kriging, acesta oferă un instrument convenabil pentru aplicații geostatistice și ar trebui, de asemenea, oamenii de stiinta de la a ajuta alte domenii.

1.2 Noțiuni de bază

1.2.1 Jos-încărcare a programului

�Puteți

(a) a descărca fișierul comprimat de la site-ul web GLOBEC;

ftp: //globec.whoi.ed/pub/software/kriging/easy_krig/V2.1/

Windows 95/98 / NT:��������������������������� easy_krig21.zip

LINUX: �������������������������������������������������easy_krig21.zip

UNIX(LINUX): ������������������������������������easy_krig21.tar.Z

Macintosh: ������������������������������������������� easy_krig21mac.zip **

și apoi executați

�������������������������������dezarhivați pentru Windows (95/98 / NT), Linux și

decomprimare + untar pentru masina UNIX-based (inclusiv Linux)

unstuff pentru Macintosh (sau software echivalent)

a decomprima fișierul. Dupã decomprimare fișierul, sunteți gata pentru a începe și a rula programul.

sau

(b) obține fișierele necomprimate originale (și sub-directoare) direct:

ftp://globec.whoi.edu/pub/software/kriging/easy_krig/V2.1/easy_krig2.1/

(folosind modul �binar pentru UNIX, automată pentru Macintosh)

** Necesită software Stuffit Expander pentru unstuff fișierul comprimat. Dacă nu-l au, puteți să-l gratuit de pe Internet: ��������������� http://www.aladdinsys.com/downloads/index.html������

1.2.2 De inițiere rapidă

Începeți MATLAB și du-te la directorul home easy_krig2.1 desemnat (pentru utilizatorii Macintosh, ar putea fi necesar pentru a crește de memorie alocată cel puțin 48MB). Trebuie doar să tastați “startkrig” în fereastra de comandă MATLAB, o fereastră se va deschide. Această fereastră este fereastra de bază, numit fereastra Navigator. Bara de meniu în această fereastră conține multe opțiuni, puteți alege. Acum sunteți gata să merg mai departe.

Notă: Puteți adăuga directorul home kriging și executați programul de la alte directoare. Cu toate acestea, trebuie să vă asigurați că nu există funcții de propria dvs. cu același nume ca și cele utilizate pentru easy_krig2.1. Având în vedere că programul vă va permite să încărcați și salva fișiere folosind un browser de fișiere, este recomandat să rulați programul din directorul easy_krig2.1 acasă.

Programul oferă multe butoane de ajutor on-line; utilizarea de aproape toate funcțiile selectabile și parametrii pot fi găsite făcând clic pe aceste butoane de ajutor (mici butoane fiecare cu un semn de întrebare “?” pe el). Este destul de auto-explicative și ușor de utilizat.



2. Etapele de prelucrare a datelor

Există mai multe de prelucrare a datelor etape (sarcini), care pot fi selectate din bara de meniu din partea de sus a ferestrei principale Navigator. Prin selectarea sau click pe oricare dintre sarcinile, o fereastră care corespunde sarcinii va pop-up. În fiecare fereastră sarcină, butoane mici marcat cu un semn de întrebare “?” furniza ajutorul on-line. Făcând clic pe butonul, o fereastră mică va deschide și furnizează informații despre utilizarea butonului parametru sau operație. Pe orice fereastră, faceți clic pe butonul “Navigator”, sau selectând “Navigator” din bara de meniu și apoi “Open” din meniul sub trage în jos va activa fereastra de navigare.

2.1 �DatePregătire

Selectarea “încărca datele” din bara de meniu de orice fereastră și apoi “Open” din sub-meniu conduce la fereastra de pregătire de date. Faceți clic pe butonul “Load” pentru a încărca fisierul de date brute. Formatul fișierului de date ar trebui să fie un fișier ASCII 3-coloană. Utilizator (i) poate oferi propria lor m-file pentru a încărca datele. Aceasta functie permite utilizatorului pentru a genera datele de intrare direct din datele brute în loc de a genera un fișier ASCII. Formatul cerut pentru intrare și de ieșire parametrii pot fi găsite făcând clic pe butonul mic patrat marcat cu un ‘?’ marca de lângă butonul “Browse” (dezactivată în mod implicit) asociate cu opțiunea “program extern”. Se poate seta ceilalți parametri înainte de a încărca fișierul. După ce fișierul este încărcat, în cazul în care utilizatorul dorește să schimbe alte setări, cum ar fi filtru de date, faceți clic pe “Aplică” după ce au fost stabilite parametrii. Diferența dintre “Load” și “Aplică” este că acesta din urmă nu va re-încărca fișierul de date. De exemplu, pentru a schimba “xlabel” sau “ylabel” după încărcate datele, este nevoie de a re-încărca datele.

2.2 Semi-variogramelor

Selectarea “variogramelor” din bara de meniu și apoi “Open” din sub-meniu conduce la o fereastră pop-up variogramelor. În primul rând, faceți clic pe butonul “Compute” pentru a genera un semi-variogramelor pe bază de date (sau correlogram). Apoi, caută un semi-variogramelor bazate pe modelul (sau correlogram) pentru a se potrivi variogramelor calculat. Relația dintre y semi-variogramelor (h) și correlogram C (h) este:

C(h) = 1 – γ(h),

unde h este distanța lag. Un variogramă (correlogram) modelul adecvat poate fi selectat în funcție de forma de variogramelor (correlogram) calculat din datele. Explicațiile parametrilor asociați cu modelul pot fi găsite folosind butonul on-line help (“?” Buton). Faceți clic pe butonul “Aplică”, iar programul va calcula pe baza modelului variogramelor / correlogram utilizând setările curente. Setările parametrilor pot fi modificate fie prin schimbarea poziției cursorului sau prin introducerea numărului direct în fereastra de editare de text. Dacă utilizați cursorul, curba teoretică se va schimba în mod automat, în timp ce în cazul în care intră în numărul în fereastra de editare de text, utilizatorul trebuie să faceți clic pe butonul “Aplică” pentru a re-calcula curba.

��������������� Dacă este instalat setul de instrumente de optimizare (recomandat), utilizatorul poate folosi mai mici pătrate funcție potrivit pentru a se potrivi datelor, făcând clic pe butonul “LSQ potrivesc”. În cazul în care caseta de instrumente de optimizare nu este instalat, programul va dezactiva automat mai mici pătrate funcția fit.

Rețineți că, deoarece se potrivesc LSQ are un anumit interval pentru fiecare parametru, diferite setări inițiale (parametri) va produce rezultate diferite.

În unele cazuri, atunci când datele sunt anizotrope, o (correlogram) modelul anizotrop variogramelor poate fi activat. Făcând clic pe butonul “2D-variogramelor permite”, setările corespunzătoare sunt activate. Dacă nu știți sau nu sunteți sigur de a raportului de orientare și aspect al caracteristicii anizotrope de date, puteți utiliza setările implicite și apoi ajustați pentru a obține rezultate satisfăcătoare.

2.3 Kriging

Selectarea “Kriging” din bara de meniu și apoi “Open” din sub-meniu deschide fereastra kriging. Parametrii variogramelor (sau correlogram) necesare funcționării kriging sunt trecute automat din variogramelor fereastra la fereastra kriging, dar poate fi schimbat de la fereastra kriging. Faceți clic pe butonul “Refresh” pentru a obține cele mai recente parametrii semi-variogramelor și faceți clic pe butonul “Krig” pentru a porni Kriging. Există două treceri de kriging; prima este folosită în producerea harta kriging iar al doilea este utilizat pentru a efectua o validare încrucișată a modelului kriging.

2.4 Vizualizare (Display)

Alegerea “Afișaj” din bara de meniu și apoi “Open” din sub-meniu conduce la fereastra de vizualizare. Există trei cifre în această fereastră: harta kriging (stânga sus), harta kriging-varianță (dreapta sus) și eco-validare (stânga jos). Primele două hărți sunt frecvent utilizate în literatura de specialitate kriging în timp ce ultima este o abordare unică a acestui program.

Faceți clic pe butonul “Arată Plot” pentru a afișa rezultatele kriging actuale și asupra “Load” pentru a încărca rezultatele salvate anterior (.mat fișier). Puteți încărca, de asemenea, fișierul selectând “File ‘din bara de meniu și selectând” Load Data “.

În cazul în care datele au fost transformate în fereastra de pregătire de date, puteți alege fie transformat (implicit) sau rezultate ale ONU-a transformat.

Prin alegerea opțiunii unice complot, o fereastră care conține numai cifra selectat va apărea.

2.5 Salvarea rezultatelor

Selectând “Fișier” din bara de meniu și apoi alegeți “Save date” pentru a salva de ieșire ca fișier .mat cu un nume de fișier definite de utilizator și la locația selectată de utilizator. Programul salvează parametrii și structuri de date (para si date) în fișierul. Pentru informații cu privire la structurile, utilizatorul trebuie să Selectează ‘File’ din bara de meniu și apoi “Structura variabilă” pentru a obține explicații pentru toate variabilele (rețineți că multe variabile nu sunt importante pentru utilizatori).

3. Exemple

Există unele fișiere de date eșantion incluse în pachet. Toate fișierele sunt în format ASCII, care este cerut de program. În următoarele două exemple, directorul krig desemnat se presupune a fi “easy_krig2.1”.

Exemplul 1: O imagine aeriană a zooplancton Abundentei date

(1) Start MATLAB; în fereastra de comandă MATLAB, schimba directorul la directorul home al EasyKrig2.1. Tip “startkrig” pentru a lansa programul.

(2) Reglați dimensiunea ferestrei și poziția și selectați opțiunea “Salvați Window poziție” din bara de meniu și apoi pe “Salvați poziție” din sub-meniu. Faceți clic pe butonul “Instrucțiuni” pentru a încărca acest document în editorul MATLAB.

(3) Selectați “Load date” din bara de meniu și apoi “Open” din sub-meniu pentru a aduce fereastra de pregătire de date.

(a) Selectați ‘latitudine’ din lista Deblocarea de fișier de date Col. 1 (implicit).

(b) Selectati ‘longitudine’ din lista Deblocarea de fișier de date Col. 2 (implicit).

(c) Lasă alți parametri neschimbate.

(d) Faceți clic pe butonul “Load” pentru a aduce o fereastră browser de fișiere, selectați directorul “date”, apoi selectați “zooplankton.dat” fișier în directorul “date” pentru a încărca fișierul selectat în memorie.

(4) ����Selectați “variogramelor” din bara de meniu și apoi “Open” din sub-meniu pentru a aduce fereastra variogramelor.

(a) Faceți clic pe butonul “Compute” pentru a calcula semi-variogramelor (correlogram). Aceasta parcele semi-variogramelor (correlogram) ca cercuri deschise discrete.

(b) Selectati “General exponențiale-Bessel” (implicit) de la modelul variogramelor listbox din dreapta sus. Faceți clic pe butonul “LSQ Fit”. Dacă nu aveți de optimizare set de instrumente, faceți clic pe butonul “Aplică”, o teoretic semi-variogramelor pe bază de model va fi suprapusă pe original semi-variogramelor pe bază de date. Cele rezultate Parametrii semi-variogramelor de a se potrivi mai mici pătrate (dacă aveți optimizarea set de instrumente) sau din setările implicite (dacă nu aveți optimizarea set de instrumente) sunt prezentate în editarea de text ferestrele în cadru “Modelul Parametri”.

(c) să faceți clic pe cursoarele asociate cu parametrii variogramelor, curbele teoretice se va modifica în consecință pentru a reflecta modificările de parametri. Parametrii de asemenea pot fi schimbate prin introducerea valorilor direct în ferestrele de editare text și faceți clic pe butonul “Aplică”.

(d) Tip 0 în fereastra de editare de text “Nugget”, 1,0 în fereastra de editare de text “prag”, 0,14 în fereastra de editare de text “lungime”, 1,5 in fereastra “putere”, și 0 în “SCL gaura” (scala gaura) fereastră de editare text. Apoi, faceți clic pe butonul “Aplică”. O curbă teoretică revizuit al semi-variogramelor pe bază de model va fi reprezentată grafic.

(5) ����Selectează ‘Kriging “din bara de meniu și apoi” Open “din sub-meniu pentru a aduce fereastra kriging.

(a) Pur și simplu faceți clic pe butonul “Krig” pentru a începe kriging prin acceptarea setările implicite. O fereastră de stare apare și afișează progresul de prelucrare. În cazul în care se face Kriging, faceți clic pe butonul “Quit” pentru a închide fereastra de stare.

(6) ����Selectati ‘Display’ din bara de meniu și apoi “Open” din sub-meniu pentru a aduce fereastra de vizualizare.��

(a) Pur și simplu faceți clic pe butonul “Arată Plot” pentru a afișa rezultatele kriging. Rețineți că pozițiile stație de pe harta kriging sunt pentru a reflecta nivelul de acordul rezultatelor kriged cu datele originale culori. Cifra cross-validare (stânga jos) utilizează metoda dublu-kriging (default).

(b) Click pe “Track de linie” din bara de meniu și apoi selectați dintre “Color”, “Black / White și” Niciuna “pentru a observa schimbarea de culoare a locațiile datele originale. Selectarea “valoare” sau “eroare”, și specificarea “Color” și “Fontsize” permite utilizatorului să afișeze valorile observate sau diferența dintre valorile observate și cele din kriging.

(c) tipul 1,0 în “Afișare prag” fereastra de text și faceți clic pe butonul ‘Arată Plot “. Există regiuni goale pe harta kriging, în cazul în care variația kriging normalizat (normalizat prin varianța) depășește specificat “Pragul de afișare” 1.0. Schimba acest număr la o altă valoare pentru a vedea cum harta kriging modificările. Maxim “prag de afișare” este definit ca “Eroare relativă” în fereastra kriging (valoare implicită este de 1,5).

(d) Selectează ‘Q1’ din caseta listă cross-validare și așteptați până când calculul se face. O linie roșie verticală la circa 0,07 între cele două linii negre la aproximativ 0,5 și -0.5 indică parametrii modelului sunt bune bazate pe criteriul Q1.

(e) Selectează ‘Q2’ de pe cruce-de validare rezultatele casetă listă într-un alt complot. O linie roșie verticală se află în regiune accepta la aproximativ 0,78, ceea ce înseamnă că parametrii modelului sunt susceptibile de “bun”.

(7) ����Selectează ‘File’ din bara de meniu și apoi selectați “Salvare date” pentru a aduce fereastra de browser de fișiere.

(a) Selectează ‘ieșire’ director, precum și orice nume doriți (implicit numele de fișier este cal83sal.mat) și faceți clic pe butonul “Salvați”.

(b) Click pe “Quit” din bara de meniu și apoi “Quit Easykrig” din sub-meniu pentru a iesi din program. Reporniți programul din nou prin tastarea “startkrig” în fereastra de comandă MATLAB. Faceți clic pe butonul “Afișare” și apoi “Open” din sub-meniu pentru a aduce fereastra de vizualizare și faceți clic pe butonul “Load” pentru a încărca fișierul de ieșire care tocmai a fost salvat cu ajutorul browser-ul de fișiere.

(c) Faceți clic pe butonul “Arată Plot” pentru a afișa rezultatele kriging, care ar trebui să fie la fel ca înainte.

(8) ��Selectează ‘Navigator “din bara de meniu și apoi” Open “din sub-meniu pentru a aduce fereastra de navigare.

(a) ��Selectează ‘Despre’ din bara de meniu și apoi “Structura variabilă” din sub-meniu pentru a deschide fișierul de descriere structură de date.

(b) ��Urmați instrucțiunile și să încerce să încarce și complot variabilele de ieșire tine.

Exemplul 2: o secțiune verticală a Salinitate date

(1) Start MATLAB; în fereastra de comandă MATLAB, schimba directorul la directorul home al EasyKrig2.1. Tip “startkrig” pentru a lansa programul.

(2) Selectați “Load date” din bara de meniu și apoi “Open” din sub-meniu pentru a aduce fereastra de pregătire de date.

(a) Selectați ‘X’ din lista Deblocarea de fișier de date Col. 1.

(b) Selectează ‘Adâncimea’ din lista Deblocarea de fișier de date Col. 2.

(c) Lasă alți parametri neschimbate.

(d) Faceți clic pe butonul “Load” pentru a aduce o fereastră browser de fișiere, selectați directorul “date”, apoi selectați “salinity.dat” fișier în directorul “date” pentru a încărca fișierul selectat în memorie.

(e) Selectează ‘inversă axa Y “pentru a inversa axa y.

(3) Selectati ‘variogramelor “din bara de meniu și apoi” Open “din sub-meniu pentru a aduce fereastra variogramelor.

(a) Faceți clic pe butonul “Compute” pentru a calcula semi-variogramelor (correlogram). Aceasta parcele semi-variogramelor (correlogram) ca cercuri deschise discrete.

(b) Selectati “liniar” din modelul variogramelor listbox din dreapta sus. Faceți clic pe butonul “LSQ Fit” [Dacă nu aveți de optimizare set de instrumente, treceți la pasul (c).] Un teoretic semi-variogramelor bazat pe modelul (o linie dreaptă, în acest caz) va fi suprapusă pe de date originale pe bază de semi-variogramelor.

(c) tipul 0 în fereastra de parametri de editare “Nugget”, și introduceți în fereastra de 1,7 parametrul de editare “prag”, și apoi faceți clic pe butonul “Aplică”. O curbă teoretică revizuit al semi-variogramelor pe bază de model va fi reprezentată grafic.

(4) ��Selectați “Kriging” din bara de meniu și apoi “Open” din sub-meniu pentru a aduce fereastra kriging.��������

(a) Pur și simplu faceți clic pe butonul “Krig” pentru a începe kriging prin acceptarea setările implicite.

(5) ��Selectati ‘Display’ din bara de meniu și apoi “Open” din sub-meniu pentru a aduce fereastra de vizualizare.��

(a) Pur și simplu faceți clic pe butonul “Arată Plot” pentru a afișa rezultatele kriging. Rețineți că pista navei pe harta kriging este de a reflecta nivelul de acordul rezultatelor kriged cu datele originale culori. Cifra cross-validare (stânga jos) utilizează metoda dublu-kriging (default).

(b) Selectati “inversă axa Y” pentru a inversa axa y.

(c) Selectează ‘Q1’ din caseta listă cross-validare și așteptați până când calculul se face. O linie roșie verticală între cele două linii negre la aproximativ 0,1 și -0.1 indică parametrii modelului sunt bune bazate pe criteriul Q1.

(d) Selectează ‘Q2’ de pe cruce-de validare rezultatele casetă listă într-un alt complot. O linie roșie verticală se află dincolo de regiunea accepta la aproximativ 0,22, ceea ce înseamnă că parametrii modelului sunt susceptibile de nu “bun” și sunt necesare modificări suplimentare ale parametrilor.

(6) ��Selectează ‘variogramelor “din bara de meniu și apoi” Open “din sub-meniu pentru a aduce înapoi fereastra variogramelor.

(a) Selectați “permite 2D-variogramelor”, pentru a permite 2D variogramelor care permite 2D anizotropie.

(b) Faceți clic pe butonul “calcul” în “permite 2-D variogramelor” caseta pentru a calcula și afișa variogramelor 2D (imagine).

(c) Având în vedere că imaginea variogramelor prezintă un model alungit având un raport de aspect de 1: 4 și un unghi de orientare cu privire la 5 grade, intrăm 0.25 în “Raportul de aspect” fereastra de editare și 5 în fereastra de editare “Angle”, și apoi click pe butonul “calcul” din nou. A 2-D semi-variogramelor recalculate prezinta practic un model izotrop.

(d) Faceți clic pe butonul “Compute” privind dreptul din nou pentru a re-calcula 1D semi-variogramelor din datele folosind parametrii anizotrope salvate automat obținute în (c). Acest pas este necesar pentru a obține corect semi-variogramelor.

(e) Reglați controlul de afișare, făcând clic pe butoanele culisante (“Lag” și “valoare”) de sub fereastra grafica.

(f) Selectați modelul variogramelor ‘General exponențiale-Bessel “(modelul 13).

(g) Introduceți 0,6 în fereastra de intervalul editare și faceți clic pe butonul “LSQ potrivesc”. Dacă nu aveți de optimizare set de instrumente, setați “prag” pentru 0,46 “lungime” la 0.34, “putere”, a 1,76 și “SCL gaura” la 0.

(h) Schimbarea “Nugget” la aproximativ 0.0165 prin introducerea numărului în fereastra de editare și apoi faceți clic pe “Apply”. Reglați scara de afișare, făcând clic pe butoanele culisante de sub fereastra grafic 1-D semi-variogramelor, lag de scară orizontală și valoare pentru scară verticală.

(7) ��Selectează ‘Kriging “din bara de meniu și apoi” Open “din sub-meniu pentru a aduce înapoi fereastra kriging și faceți clic pe butonul” Krig “pentru a re-calcula harta kriging.

(8)�� Selectează ‘display’ din bara de meniu pentru a aduce înapoi fereastra de vizualizare.

(a) Click pe “Arată Plot” pentru a actualiza hărțile.

(b) La fel ca în (5), selectați “Q1” și “Q2” din nou, am descoperit că atât Q1 și Q2 sunt acum în regiunile acceptabile. Rețineți că, prin utilizarea unui model anizotrope, unele îmbunătățiri pot fi realizate, dar harta kriging sine nu se schimbă mult.

(9) ��Selectează ‘File’ din bara de meniu și apoi selectați “Salvare date” pentru a aduce fereastra de browser de fișiere.

(a) Selectează ‘ieșire’ director, precum și orice nume doriți (implicit numele de fișier este cal83sal.mat) și faceți clic pe butonul “Salvați”.

(b) Click pe “Quit” din bara de meniu și apoi “Quit Easykrig” din sub-meniu pentru a iesi din program. Reporniți programul din nou prin tastarea “startkrig” în fereastra de comandă MATLAB. Faceți clic pe butonul “Afișare” și apoi “Open” din sub-meniu pentru a aduce fereastra de vizualizare și faceți clic pe butonul “Load” pentru a încărca fișierul de ieșire care tocmai a fost salvat cu ajutorul browser-ul de fișiere.

(c) Faceți clic pe butonul “Arată Plot” pentru a afișa rezultatele kriging, care ar trebui să fie la fel ca înainte.

Exemplul 3: �proces discontinuu de Densitate Potential la diferite Adâncime

(1) Start MATLAB; în fereastra de comandă MATLAB, schimba directorul la directorul home al EasyKrig. Tip “startkrig” pentru a lansa programul.

(2) ��Selectați “Load date” din bara de meniu și apoi “Open” din sub-meniu pentru a aduce fereastra de pregătire de date.

�(a)�� Selectați ‘latitudine’ din lista Deblocarea de fișier de date Col. 1 (implicit).

(b) Selectati ‘longitudine’ din lista Deblocarea de fișier de date Col. 2 (implicit).

(c) ��Selectează ‘program extern “, browser-ul program extern este activată și apare o fereastră de nume de fișier.

�������(d) Faceți clic pe broswer dosar, pentru a schimba directorul extern_program și selectați “remove_nan.m”, apoi faceți clic pe “Open” permite programul pentru a accesa programul (Dacă acest program extern nu este ales, locațiile în care valorile datelor sunt de NaN va fi afisat pe harta). Acest program oferă un exemplu pentru a scrie propria rutină pentru a obține formatul datelor necesare pentru EasyKirg dintr-un fișier de date arbitrar.

���������������(e) Lasă alți parametri neschimbate.

(f) ���Faceți clic pe butonul “Load” pentru a aduce o fereastră browser de fișiere, selectați “de date / pdensity” director, și apoi selectați “layer12.dat” fișier în directorul pentru a încărca fișierul selectat în memorie.

(3) ��Selectati ‘variogramelor “din bara de meniu și apoi” Open “din sub-meniu pentru a aduce fereastra variogramelor.

(a) Faceți clic pe butonul “Compute” pentru a calcula semi-variogramelor (correlogram). Aceasta parcele variogramelor semi (correlogram) ca cercuri deschise discrete.

(b) Selectati “General exponențiale-Bessel” (implicit) de la modelul variogramelor listbox din dreapta sus. Faceți clic pe butonul “LSQ Fit”. Dacă nu aveți de optimizare set de instrumente, faceți clic pe butonul “Aplică”, o teoretic semi-variogramelor pe bază de model va fi suprapusă pe original semi-variogramelor pe bază de date. Parametrii de a se potrivi mai mici pătrate (dacă aveți optimizarea set de instrumente) sau din setările implicite (dacă nu aveți optimizarea set de instrumente) sunt prezentate în editarea de text ferestrele în cadru “Modelul Parametri”.

���������������(c) ��Puteți modifica modelul, după cum doriți sau acceptați parametrii de a trece la pasul următor.

(4) ��Selectați “Kriging” din bara de meniu și apoi “Open” din sub-meniu pentru a aduce fereastra kriging.��

(a) să stabilească coordonatele ca:

�����������������������������������������������minimă orizontală a -22

��������������������������������orizontală maximă la -12

�����������������������������������������������rezolutie orizontala de 0,4

�����������������������������������������������minim vertical la -25

�����������������������������������������������verticală maximă la -20

�����������������������������������������������Rezoluție verticală la 0,2

(butonul radio) (b) Selectează ‘Procesare lot “, pentru a permite procesarea lotului.

(c) Faceți clic pe butonul “Browse” de lângă “Load-list Filename dosar” pentru a iniția browserul de fișiere. Schimba directorul la “date” și selectați “batch_file.dat” pentru a încărca fișierul lista filename în care este specificată o listă de nume de fișiere de date.

(d) Faceți clic pe butonul “Browse” de lângă “Salvați fișierul jurnal lot” de a iniția browserul de fișiere. Schimba directorul de “ieșire”, sau oriunde doriți. Introduceți un nume de fișier pentru fișierul jurnal, în cazul în care este salvat de informație de intrare / ieșire.

(e) Faceți clic pe butonul “lot Krig” pentru a începe lot kriging. Fișierele de ieșire vor fi salvate în același director ca și fișierul jurnal procesul de lot.

(5) ��Selectati ‘Display’ din bara de meniu și apoi “Open” din sub-meniu pentru a aduce fereastra de vizualizare.

(a) Faceți clic pe butonul “Arată Plot” pentru a afișa rezultatele kriging din ultimul set de date: “layer14”.

(b) Verificați de ieșire pentru celelalte straturi. Pentru a verifica rezultatele de prelucrare, aveți nevoie pentru a încărca fișierele de ieșire. �Faceți clic pe butonul “Load”, apoi schimba directorul în care este salvat fișierul jurnal. Încărcați unul dintre fișierele de ieșire (“layer10.dat” – “layer14.dat”) și afișa prin click pe butonul “Arată Plot”.

(c) În cazul în care nu sunteți mulțumit de rezultatele kriging, poti sa te duci înapoi la fereastra “variogramelor” și sau fereastra / “Kriging” pentru a modifica parametrii implicate, și re-kriging până când sunteți mulțumit. De exemplu, prin setarea pepită = 0,05, pervazul = 1,16, lungimea = 0.253, putere = 1,5, iar gaura SCL = 0 în fereastra variogramelor și efectuarea din nou kriging, veti obtine un rezultat mai bun, adică, Q1 și Q2 rezultat va intra în regiunile acceptabile.

4. Referințe.�

Deutsch, C. V și AG Journel, 1992. GSLIB: Geostatistical Software Library
și Ghidul utilizatorului. Oxford University Press, Oxford, 340 p. (Deutsch, C. V and A. G. Journel, 1992. GSLIB: Geostatistical Software Library and User’s Guide. Oxford University Press, Oxford, 340 p.)
Journel, AG & CJ Huijbregts, 1992 Utilaje Geostatistics. Academic Press, New
York, 600 p. (
Journel, A.G. & C.J. Huijbregts, 1992. Mining Geostatistics. Academic Press, New York, 600 p.)
Kitanidi, P.K. 1997. Introducere în Geostatistics. Aplicații în hidrogeologie. Cambridge University Press. 249 p. (Kitanidi, P.K. 1997. Introduction to Geostatistics. Applications in hydrogeology. Cambridge University Press. 249 p.)
Marcotte, D. 1991. Cokriging cu MATLAB. Calculatoare & Geosciences. 17 (9): 1265 – 1,280. (Kitanidi, P.K. 1997. Introduction to Geostatistics. Applications in hydrogeology. Cambridge University Press. 249 p.)

Comments are closed.