Idrisi Andes DjgiííjJľJj j/JDíJyJy -terénu (Q-rj) DJVJT v GJS J dri si And *-\ ^* Ing. Martin KLIMANEK, Ph.D. Ing. Petr DOUDA 411 Ústav geoinformačních technologií Lesnická a dřevařská fakulta, Mendelova zemědělská a lesnická univerzita v Brně n^^% OJSÜ22?* ^OldO^. DjgjííiJľjj moíÁBbj í^ľáiJLJ (o) Tvorba DJVJT DMT 6-8 ŠLP „Masarykův les &» ifti yí< í 5 Data: Idrisi Andes i" Krtiny ^S-J' irA fcjfe*F ä íl. EL/ I 1 ml- výškopisu ZM ČR 24-41 1-16 1(1:10 000) S:\DMT\ (\\195.178.78.202\student\dmt\) Idrisi Andes Import dat do Idrisi Zkontrolovat: Start - Ovládací panely - Místní a jazykové nastavení - Desetinný oddělovač „ . " (ne „ , ") Spustit Idrisi Andes Idrisi Explorer - Projects Idrisi - File - Import - Software-Specific Formats - ESRI Formats - modul SHAPEIDR Struktura souboru (vektorové kolekce) Database Works ho. Vytvořit vektorovou vrstvu z atributu výšky Zjistit maximální a minimální souřadnice pro výpočet velikosti pixelu (X v Idrisi = columns, Y v Idrisi = rows; ve skutečnosti je ale souřadnice X souřadnicí Y v S-JTSK a naopak) • Min X: -593 665.4 m, Max X: -588 438.5 m • Min Y: -1 157 834.6 m, Max Y: -1 153 544.8 m Idrisi Andes INTERPOL INTERCON TIN Interpolation ► Kriging ► TH1E5SEN TREND TIN TIN5URF GENERALIZATION LINTOPNT TINPREP Interpolation ► Geostatistics ► Topographic Variables ► Feature Extraction ► 5LOPE ASPECT HILL SHADE CURVATURE FRACTAL Spatial Dependence Modeler Model Fitting Kriging and Simulation CONTOUR TOPOSHAPE PIT REMOVAL RUNOFF FLOW RU5LE WATERSHED SLOPELENGTH SEGMENT SEDIMENTATION DMT 6-8 Modul GENERALIZATION rJh líj';iJZ;.7JrJJ 1 ■ ŕi'.Vjr -.inú J1.1 -J ,:;.'Jj:.Jj!lJ^'.!rJVJJ File format type 1 Vector option f* |Uw !^.^^!Í^^| r Point generalisation Input vector Hb: Output vector Be: J J 1 . 1 . Generalisation option <* Point selection f* Low-pass fätering ** Tolerance band Selection rete 1 Tide: ÜK Ckw Help Idrisi Andes l£jlALJi;TJrJJ I - t-il VjS mjj'J /■.;.-l-a yjťJjyj'ulíiísjiJyjj File lt.irviai lj.ip& f*1 Vector generalisation c Raster generalisation Vector option C line generalisation (i jpcjirrt generalisation! Input vector He: Output vector We: Tolerance distance: Tide: Sä A A OK I Gore Help Idrisi Andes Modul TINPREP a LINTOPNT 1 f* Add points along tie <" Delete pants afang line Input vector line fte: 244116>e ^ OulpUl vector (ne lie: J244116_in*s J Tolerance diiünoe. (• Absol Je distanc f Relative distanc jelinilion: e Tolerance: |20 OuipUid«yfftenta(ion.. I OK Close Help Idrisi Andes t»-" * ■ ■ ■ ■ ■ ■ " ř i . " fr I tĽ4T + * ' t » . , '•■•■■ j I nov »n BI ti mm sorm SSáJJ JMErt as k 4mu i" ". wja ■ DMT 6-8 Modul INTERPOL Jj ITíj'j-'OL - jUr/iJĽ-í jjji-j>:Ajijujj Function <* InteJpolate digital elevation model <~ Calculate potential sof«* Input vector He: |2441l6_pdnts j Uutput mage: |24411 Gjnterpol ... Distanc« wí*jW exponent: B.O f U» vector port IP'sas betohtí T U« afofcute valued for height*: R? Use a au point search radus Idrisi Andes DMT: pixel 10 x 10 m H Output Be ipeciňcalkms H«i«nX: I-593665.4375 Mrtiwn Y: 1157834 625 MaďnumX: -588438.5625 M&ďmmY: -1153544.75 - Output ne specficatiom Dala Type Columns Row«: V Copy from enctiig 13a: [«ear 3 523 429 Output c&umenlslnn... |-------PK-------j| Close Help Modul RASTERVECTOR Idrisi Andes n aster/Vector - Vector/rl aster (* Vector to raster C Rastet Iq vector _ í~ Point to raster (• Line to laster C Polygon to raster Vector line file: Image file to be updated : |244116Jnes -I |2441ie_lines J Het Of Close i\ (•' Copy spabal parameters from another mage C Deine spatial parameters ndhidualy Output mage: | D:\dmft24411 6_.nes.rst "1 Image to copy parameters from; 124411G nterpo 1 ^J Output data type: Initial value: |ieal _l 1» Output documentation... 1 H* OK Close DMT 6-8 Modul INTERCON File type í* Image Input file name: Output file name: Current data type: C Vector |244116_lines |244116_lines2 Real Current file type: Conversion type (if applicable] C Truncation <* Rounding Binary Output data type Output file type <* Integer <~ ASCII C Real (* Binary C Byte £~ Packed binary OK Clo;e Help Idrisi Andes Input image: Background value: 24411S Iines2 A Height of the top left corner: Height of the top right corner: Height of the bottom left corner: Height of the bottom right corner: J350 ;335 29Ů |385 Output image: 24411S intercon OK J Close- Help Idrisi Andes Unconstrained (left) and constrained (right) Delaunay triangula- tions. Solid lines represent isolines. B/T Edges b Critical Points Re-Tri angulation c d a: Contours at the top of a hill; b: triangulation of highest contour, with BIT edges identified; c: placement of critical points on B/T edges; Tunnel Ectees HritiL^- l .lIí^ľ i in m 100 m 500m ÓOOm ĎOOm SOOm a b a: Contours at a stream; b: contours at a "saddle" feature. Contours are shown with solid lines, constrained triangle edges with dashed lines. B/T edges are shown in red. hO hO+Ah hO+2Ah hO hO+Ah hO+2Ah PO^...................h\............. r'p2 Jj*J hO hO+Ah hO+2Ah PO hO hO+Ah 'pí / V. +Ei...........i,n f-i v / v ^ hO+Ah hO f^ fe/w£l \ P2 4 PO/PÍ \^. f f .^ / hO š / PO \ hO+Ah ..P.U...ÍP2 \J hO+2Ah hO+Ah hO A H P2 ........................ ^JM PO.-* ir------------- ----------*s H PI/ #0 ^ÉP2 " H PÍ____P2 PO......-IK^----- ■■■■* Idrisi Andes Modul TIN Samce of points for iriangulaHon: ■f» Une vertices í" Poinls Type o\ iii-angulelisn: f" Cpnsdained (• Notvconsirained |2«11S_fc»S Input veetff lie: QHputTINIle: |244116_tírt !^ Perform Biidge .and Tunnel Edge removal Biidge arid kainel «dg» removal; <* Parabole f Optimized ünea( r- Ljieaf 1^ Cieate rasier surface Raster lite name: dmt_24411 SJnetinJ1 Acfd comet points Oulput documentation... 0ÍC CID» He* ^J J J rjii;'ji!f L--illrJ:!l It! - ■ J J T . '.t ľj'.Jli P ]j I -.S' Input TIM He; Output rasier surface image: B^ckgrrxnd value Jj 2M11S_tň |drrt_24411&_linet*i_P P" Output docufnenrjiion., 0Í Close Help 'J' J Output lie špecifickom: Copy fiorm e/islríg file: MinimimumX: Maximum X: Minimum Y: Maximum Y: Number d columns: Numbs: d iows: I -5936654375 -588433.5625 -11E7834.G25 ■1153544.75 523 |*29 totéž provést pro jiné varianty „B/T edge removal" a pro bodová data analýza rozdílů (modul OVERLAY, HISTO) Idrisi Andes Modul TREND Idrisi Andes ľľí) '!■' - !J-JjJ -■'.)!! : ■ ■ ■.;»-'...-!. ._J Input image: 244116jhes j| Output image: [244116jiend J Older ůf surface to be lilted: 1* á 1 I- Include rao data eels TiHe Que 1 UK l| l Hefp Idrisi Andes Modul THIESSEN nejprve relativně 5x redukovat počet lomových bodů vrstevnic (TI N PREP) vyjmout lomové body (LINTOPNT) a konvertovat formát dat (CONVERT) Input vector line file name: Output vector point file name: OK File type C Image Input file name: |244116_lines2 |244116_points2 Llü'.f- Melp & Vector 1244116_points2, vet J J J Output file name: |24411S_points2_con ^ Current data type: Real Current file type: Binary Output data type <* Integer P Real Output file type C ASCII vector export formet (.vxp) <* Binary I_________,___________________^_ Conversion type (if applicable] f~ Truncation ff Rounding OK ] Oca | Help ľ r 11 ^. ■ .-■ E11 - Tii i -.■.■ n ^'j1 ..-.■.■.i 11 i ■■■■! ■:!;■ i; r ■ .-.■! j'ISÜmJi Point file type: í Vector file <~ Image file Input vector file: Output image file: 244116_points2_con |244116_thiessen Output documentation.. J [:::::::::::ie::::iii ciose h^ j j Boundaries of study region: Minimum*: ■533635.4375 Maximum X: ■538433.5625 Minimum V: ■1157334.625 Maximum Y: ■1153544.75 Output file parameters: Data type: Columns: Rows: iri^-Hí-r w 1523 14 29 Analýzy DJVJT Idrisi Andes modul CONTOUR - extrakce vrstevnic modul PROFILE - řezy terénem modul HILLSHADE - vytvoření stínovaného reliéfu modul SURFACE - výpočet sklonu a expozice svahů modul CURVATURE - výpočet horizontálního a vertikálního zakřivení modul FRACTAL - výpočet fraktálních dimenzí modul SEGMENT - homogenizace (segmentace) rastru DMT modul VIEWSHED - analýza viditelná modul WATERSHED - výpočet povod modul FLOW - směry povrchového odtoku modul RUNOFF - akumulovaný odtok modul TOPOSHAPE - identifikace tvarů terénu modul RUSLE - revidovaná universální rovnice ztráty půdy, včetně možnosti samostatného výpočtu souvislé délky svahu v lokálním povodí (SLOPELENGTH) a ztráty a depozice půdy (SEDIMENTATION). Modul CONTOUR Idrisi Andes .-'jjrrtiUiť Inpu lasia image: OiMpui vecl« lie: Minmum contour value: Manmum corfou value: ConAou irrieival 17 Geneiatze contour dala T Define background aieas lhat stadd not be pioce;;«d J JJL-JJJ J-'. !-'J !'£. . IBlalta dn^2441lG .=] | content ^ 2CC1 480 |20 TUs: I OK II do» Help vytvořit kompozici s DMT a původními vrstevnicemi vizuálně porovnat DMT 6-8 Modul SURFACE .'Uí'r,.'. £ - _- LJ j f-.:1. í -.: j j -.: 1 .-■ _- J _ ■BIBIBI Calculate: C Slope : C Aspect f.: clrjsfíf idcNom of [o]iojydjihi|i' io.rhirei OiijxJiraoe: ftow)d«yt0nďÍQ ',=1 a "^1 'á Fí tertjnw Mw r Gatathn r MMlMA r MtflJM r Ms r. m r Adaptiv« Bo« r Male r SlandadOeviéíůn f L*pl*OBnEdöeCnh*n«meiH r HítflPa» r SobdEdgsDtlcda r UifltHnsd |v«mfcle u»knns| modul FOURIER modul FILTER Idrisi Andes Modul RUSLE Idrisi Andes J'.IĎLÉ -ľ.*/Lvi uij];-tfj. :'.\ j'j}\ \'jjj ■-"■]' J ■: í i -j J i J Jj Input Files DEM image: fv Use field image: Precipitation image (R factor): | _j Land-cover image (C factor): ...I Sou image [K factor]: _J Management image (P factor): | _J Control Specifications Slope threshold {%): ±1 Aspect threshold (deg): Smallest patch size: _ -H«' Maximum slope length: -H|F«,I jj C Rounded to shorter (* Rounded to longer I- Average soil factor within patches Default background value : |0 Output Files Prefix Specifications Patch output files prefix j J Field output files prefix : -1 Retrieve parameters Save parameters czirnzzii Close Help Suk, P. Možnosti kvantifikace erozního ohrožení půdy v GIS Idrisi. [S.I.], 2005. 55 s. LDF MZLU v Brně http://mapserver.mendelu.cz/skripta/GIS_IV/CVIC05.pdf DjgjííiJľjj moíÁBbj í^ľáiJLJ (u) Geostaiisticke násirojs Idrisi Andes Gstat v Idrisi Andes Gstat (http://www.gstat.org/) je integrován do 3 grafických uživatelských rozhraní Modul „Spatial Dependence Modeler" - analýza složek prostorové variability dat • experimentální stanovování šířky a počtu tříd (lag) pro výpočet variability jednotlivých párů • prostorová variabilita se chová isotropně nebo anisotropně -všesměrové pole se rozděluje na směr minimální variability a směr maximální variability (při dané úhlové toleranci, např. 22,5°). Spatial Dependence Modeler v^ivtzlíull u •Tsjíľi lis« i K\ lili semivariogramy pro směry minimální a maximální variability (soubory s příponou var), Display Type Method Layer Type Variable Directional -1 [Sernivariograrn (Moments Est.) - | [Vector A 1 : b050G_vrst_ Raw t| Zero Separation: Surface chart pal. Use method default -| h-scatterplot Variability Minimum ■ Maximum Direction: 133.8 Graph ] Save About Help Close Data Transformation Lags JW° Cutoff X ~T] [Regular T] _J |i[ 4 000- j 3 500 3 000 Jy 2 500-y 2 000 f 1 =00 1 000 y^ ^j***^ 500 0- 500 1 000 1 500 2 000 distance between pairs 2 500 Series (t 135° C Omni r 45* r (ľ Stats C On (T Off Redraw Gstat v Idrisi Andes Modul „Model Fitting" - stanovení a výpočet teoretického semivariogramu Idrisi Andes matematická funkce (11 variant), která nejlépe vystihuje zobrazený semi-variogram • vizuálně co nejlépe stanovit její parametry (dosah, práh a zbytkový rozptyl). 1 lze modifikovat i měřítko anisotropie a přidat do grafu dva další semivariogramy. • systém poté automaticky zpřesní tyto parametry pomocí jedné ze tří metod (dvě váhové metody nejmenších čtverců - WLS a metoda maximální pravděpodobnosti - REML). ;ii: Model Fř Sample variogrann (o fit model (o (.var |*| elev45maH Optional files [for display only) Sample variograrm 2 : |_^J □r a rMe Sample variograrm 3 Method:— WLS1 C WLS 2 C REML 5.003170 |Spherical 7] 1509.744 Structure 2 |_J J j\ \JJ~ jjj [7J Spherical t| 0 Structure 3 pTj ~3UD JJ UTZB Spherical »J 0 Structure 4 [T]~ | ~^j |7j~ | ~~j] d urzB \7 Change number of lags F^ T jLag stats; Iteration Limit: |50 ^j Fit Limit: Combined Structures 200 400 600 distance between pairs Angle 1 : [45 ' Angle 2: Independent Structures 200 400 £00 distance between pairs Gstat v Idrisi Andes Modul „Kriging and Simulation" - prostorová interpolace povrchu Idrisi Andes • výběr z variant typu výpočtu povrchu (DMT) • pro každý z teoretických semivariogramů vznikají až čtyři rastry. Dva obsahují hodnoty predikovaného povrchu (jeden pro souvislý povrch a jeden pro bodové hodnoty křížové validace) a zbývající dva obsahují hodnoty variancí (opět pro povrch a pro křížovou validaci). • kromě těchto rastrových výstupů je možné vytvořit i textový soubor se statistikou vstupních dat a predikovaného povrchu. Estimation and Statistical Options Simple Cokriging Ordinary Cokriging Gaussian Simulation Universal Kriging Indicator Simulation Local Linear Model Prediction Cross-validate on sample locations? "" Block average variogram values Size in x: Size in v : I I- Krige more than one stratum How many? FT Model Specification "Model source: — ß" .prd File: f edit existing or new variograms Input File(s) Input Data File: Data transformation: Local neighborhood Select options: I- minimum number of sample points ~~ maximum number of sample points | radius for sample selection I- force neighborhood selection |"~ use variogram distance for neighborhood selection "" set maximum observations per quadrant Mask File [Optional if Input Data File is raster] Output File(s) (Only first output filename required] I------------------_l Prediction File About Help Idrisi Andes ELEV2 - - 45° směr minimální variability (maximálni kontinuity) ORDINARY Kriging corr(Obs, Pred) : 0.9986 [using ord: in ar y kriging] observed predicted pred.-obs. pred.std. sscore minimum 260 260.3 -14 3.126 -1.2 1st q. 340 336.8 -0.8164 6.157 -0.1143 median 360 360.2 -0.005162 7.177 -0.0007764 3rd q. 380 381.6 0.6987 8.273 0.1067 maximum 420 420.8 13.44 16.08 1.686 n 2402 2402 237 8 2402 237 8 mean 357.8 357 .8 -0.01592 7.288 0.001087 std.dev. 37.78 37.36 2.048 1.668 0.2498 ELEV3 - -135° směr maximální variability (minimální kontinuity) ORDINARY Kriging | corr(Obs, Pred) : 0.9978 [using ord: in ar y kriging] observed predicted pred.-obs. pred.std. sscore minimum 260 260.8 -11.2 5.695 -1.246 1st q. 340 335.9 -1.397 7.14 -0.17 8 median 360 360.3 -0.02638 7.7 07 -0.003044 3rd q. 380 382.2 1.213 8.323 0.1631 maximum 420 421.2 14.7 5 13.17 1.67 2 n 2402 2402 2402 2402 2402 mean 357.8 357 .8 -0.04454 7.811 -0.002429 std.dev. 37.78 37.27 2.519 0.9451 0.3104 textový soubor se statistikou vstupních dat a predikovaného povrchu (příklad) I - Spatial Dependence Modeler 1) zobrazení dat Idrisi Andes jiijij Rainfall Data, West Africa o ♦ 33.00 * ♦ 45.19 • * 57.38 ♦ * * 69.56 o • O 81.75 íl • # 93.94 O o 0 O 106.13 O o * O 118-31 O 130.50 O O * o íl O 1^2.69 O 154.88 3"= O o * o • íl Q 167.06 o íl * o • o 0179-25 o O ° O 0 o 0 íl o ■íl íl o í (j| 191.44 ° • • o O ° ?=0 O o O O ° Q) o 9 % * o o O o ' 6 *v* O O 5 £f 203.63 ^■215.81 ■■ 228.00 o o° O n O ° o% o if oo^ " * O oo o o ° O 00 O • O ° o o O O o o o o D ° f? O o ° ° o o 0 ooo o o • o» o u n ° • 0 0° g O Oo °o ° 0 ° o O o ° O0 O • • • • U O o • * O O • * o ~ • • RAIN Display Parameters | Properties ] Visibility Autoscaling Options C None (direct) Number of Classes ^ Equal Intervals i Quantiles Standard Scores |25G _^J Display M in/Max Contrast Settings ., [33" Display Min /.......... DMay Max '................J (ŠäŠ" Symbol File: | quant _J Advanced Palette/Symbol Selection Apply Revert Save ■OK" a Close Help Idrisi Andes 2) Spuštění Spatial Dependence Modeler GIS Analysis -» Surface -» Geostatistics -» Spatial Dependence Modeler a) zvolíme RAIN jako vstupní proměnnou („Variable") b) zkontrolujeme nastavení „Display type" na „Surface" c) necháme vytvořit „Variogram surface graph" stiskem tlačítka „Graph" Idrisi Andes d) hledání prostorové závislosti ve srážkových datech c;::*"»I Dependence Modeler [Displai) Type ^ Method Layer Type Variable Data Transformation Lags iDirectional lemivariograrn (Moments Est.) Ťf |Vector ~^\ 1 : |rain ... | [Ňč na těchto datech je vidět jeden za základních axiomů geografie: „Data", která jsou si v prostoru blíže si budou podobnější než ta, která jsou od sebe dále. Idrisi Andes ^JjíIÍJ::! Ud\jdllÚ'-lltd jji'j'Jikj" .atterplot Series : Lag: |1 ^ Graph Lag i Threshold value: 0.0 Appli) Threshold ^J Zoom 100% Cancel Display Type Method Layer Type |Directional _»J |Sernivariogram (Moments Est.) _rj |Vector _^J 1 : r; Raw ~B Zero Separation : |Use method default -řj Surface chart pal. : |idris256 Variable Data Transformation Lags Cutoffs ■ I JNone _J |Regular T] _J J33.33 1 2364 | ' Se,ies ~ (• Omni Help 220- * 210- . . .* . .• 200- * . ** 190- * *v* • * 180-170- •» • * *• » «*f 911 * * * * H ** < 160-^ 150 •«VJtWf *: ** H | 140-: 130 O * * t * H ** b- 120 * 0 . _ * *** ** • * 110- •.....** • 100- h t *v* :* 90 * * • * 80 * * * * * * 70- * * * 60 * * * * 50- * * 40- 4 0 60 80 100 120 140 160 180 200 Trom" Valne 220 DMT 6-8 - Lag Specification" C" Automatic (* Manual w: Number of lags : |40 Lag width Cancel 20 mrecrionai t | pennivanogram Use method default wl thod Layer Type (Moments Est.) t| Vector ^1 1 : rain Raw Zero Separation : Surface chart pal. Jidris256 _l h-scatterplot Variability Minimum ■ Maximum Direction: 114. S 0 Lag: 9.4 Graph ; Save About Help Close Idrisi Andes Variable Data Transformation Lags | |None t| [Regular "3. 100 Omnidirectional override' Enter direction angle graphicallv with mouse or enter value ISO Direction angle: [ŠŠ Angular tolerance: É------- 82 první interval o šířce 20 km obsahuje pouze jeden datový pár, se vzájemnou vzdáleností 11.26 km —» lze předpokládat, že prvních několik intervalů bude méně spolehlivých (málo datových párů), proto se budeme snažit, aby v každém intervalu bylo alespoň 30 párů (získání reprezentativního průměru) Idrisi Andes jU;;!J:j! DiťUijjdijj':^ ii'sjúaížir Display Type Method Layer Type [Directional t Sernivariogrann (Moments Est.) -r i IVector ^J 1 : Irain Raw tJ Zero Separation : |Use method default t| h-scatterplot Surface chart pal. Jidris256 J Series Statistics Lag Statistics Lag Pairs Average Distance 1 11.2633 •mm 3 4 5 6 7 8 9 10 11 12 13 14 15 16 IT 18 32 34 48 47 51 58 49 60 73 81 65 56 58 76 .01 ■ 50.4571 70.2812 91.0488 109.253 130.182 150.556 170.327 189.231 208.881 229.481 250.365 271.276 290.411 311.371 331.192 Si Graph Save About Help Close Variable EJUS Data Transformation Lags Cutoffs _J |None _*J JRegular T] _J |l00 200 400 600 distance between pairs Redraw 180 I- Omnidirectional override Enter direction angle graphically with mouse or enter value Direction angle: |95 Angular tolerance: É------- Idrisi Andes j'.üy.ij!ijj; JjjÍťíj " Lag Specification-f~" Automatic JManuaí OK Number of lags : 20 Lag width: [40 Spatial Dependence Modeler Display Type Method Layer Type Variable ||Pirectional t\ [Sernivariogram [Moments Est.) t j I Vector "▼] 1 : Irain lJ1 Raw ~3 Zero Separation : |Use method default ^| Surface chart pal. : Data Transformation Lags Cutoff % _J [None _^J |F!egular T] _J 1100 _ Series-ď Omni r 35* r 35' r 5- *jjd • Stats — T On (T Off Redraw r, Direction: 62.5 0 Lag: 10.0 Maximum Graph j Save About Help Close 180 17 Omnidirectional override Enter direction angle graphically with mouse or enter value Direction angle: Angular tolerance: É------- necháme si zobrazit graf pro směr 95°, 5° a Omnidirectional Idrisi Andes uložíme „směrové" grafy jako experimentální semivariogramy: i) variogram 95° -> RAIN-major-95 Ü) variogram 5° -> RAIN-minor-5 iii) variogram omnidirectional -> RAIN-omni ■Save directional model to (* Variogram file (.var) C Clipboard C Printer Enter name of file to create (no extension): JDADATAMdrisiVgeostatistika OK Cancel nyní máme dostatek informací, abychom mohli odvodit prostorově souvislý model pro srážková data DMT 6-8 II - Model Fitting Modeling Zonal Anisotropy • jako „Sample Variogram" zvolir • pokusíme proložit křivku body t ii\'jU'ť\ FJÍÍJJJLJ Sample variogram to fit model to (.var) |*1 ]rain_major-95 Optional files (for display only) : Sample variogram 2 : |_fj I Sample variogram 3 : j j deg, J Ranges CK) Sills Anisotropy |Y) Ratios jNugget Structure 1 70 (Spherical t| Structure 2 20000 1 0.00001 [Exponential t| Structure 3 GO ELT 500 ±]W J±J UT3Ľ Power t| Structure 4 Q" 750 ET Q" Fit Model View Log Close Idrisi Andes ne RAIN-MAJOR-95W ak, aby je co nejlépe vystihovala - n x "Method:— <* íwlsj1 T WLS2 T REML I- Change number of lags I- Lag stats Actual Sills 1 0.3333 2 166.65 3 249.975 Iteration Limit: Fit Limit: |l.0e-S Save Model Append to... About Help Combined Structures i i i i i i i i 200 400 600 distance between pairs Angle 1 : |95 Angle 2 : J95 Angle 3 : Í95~ Independent Structures 200 400 G00 800 distance between pairs DMT 6-8 - Ordinary Kriging Estimation and Statistical Options [Ordinary Kriging "3 Kriging Options ll? Cross-validate on sample locations? I Block average variogram values Size in y: Size in x: Krige more than one stratum How many? 1 -^j I Model Specification Model source : r prd File: |rain-pred ... | ídit existing or new variograms Specify model(s) to use : Var. 1 Var. 2 ,5,0.00001) + 166.6667 Exp(60) + 250 Pow(2); (I 0 Nug(0) Covar. 0 Nug(0) Sph 1 Exp | Pow | Lin | Log Gau | Cir | Bes | Pen | Per | a) cross validation Idrisi Andes Input File[s] Input Data File: J? (• Vector Raster Data transformation: 31 Local neighborhood Select options: minimum number of sample points mm. E maximum number of sample points |~~ radius for sample selection |~~ force neighborhood selection |~~ use variogram distance for neighborhood selection set maximum observations per quadrant 30 Mask File (Optional if Input Data File is raster) rainmask H Output File(s) (Only first output filename required) RAIN-XL-PR E D Prediction File OK Output FiSe(s) (Only first output filename required) RAIN-XL-VAR Variance File Idrisi Andes b) interpolace souvislého povrchu J.JJL--JJJL-- mel Simulation Estimation and Statistical Options Ordinary Kriging ~B Kriging Options 11- Cross-validate on sample locations? I |~~ Block average variogram values Size in x; I Size in y : |~~ Krige more than one stratum How many ? [i -^-j Model Specification Model source: — .prd File: |rain-pred (* edit existing or new vasograms _l Specify model(s) to use : Var. 1 170.000000 Nug(0) + 0.3333333 Sph(20000,5,0 G Var. 2 0 Nug(O) Covar. |0Nug(0) Sph | Exp | Pow | Lin | Log Gau | Cir | Bes | Pen | Per | IHIBIB Input File(s) Input Data File: rain Data transformation: _l (* Vector f* Raster ~B Local neighborhood Select options : minimum number of sample points p1 maximum number of sample points max. radius for sample selection force neighborhood selection use variogram distance for neighborhood selection set maximum observations per quadrant 30 Mask File (Optional if Input Data File is raster) rainmask ... | Output File(s] (Only first output filename required) D rain-pred Prediction File OK Output File(s) (Only first output filename required) Variance File Idrisi Andes c) Porovnání interpolovaného povrchu se zdrojovými daty pomocí modulu OVERLAY O'/ĽíL First innage: Second image: Output image: rain-samples rain-pred ...I | D:\DATA\ldrisi\geostatistika\cover. ^J Overlay options: First + Second C~ First to the power of the Second First - Second P Minimum FirstK Second f Maximum First 1 Second (• First covers Second except where zero First - Second I First + Second Output documentation... ■OK'" Close Help 8204