Download - Alunecari Mtii Suhard
Determinarea probabilității apariției alunecăriilor de teren folosind analiza spațială GIS
Studiu de caz: Munții Suhard
Îndrumător: Dr. Ștefan BILAȘCO
Masterand: Szilárd KOSZ
Date de intrare
• Modelul Digital de Elevație
• Limita unității de relief (Munții Suhard)
• Rețeaua Hidrografică
Adâncimea fragmentării 1
• Data Management Tools -> Feature Class -> Create Fishnet ( cu laturi de 1 km pentru, formula de calcul al adâncimii fragmentării se referă la kmp )
• Spatial Analyst Tools -> Zonal -> Zonal Statistics cu precizarea fishnetului și a modelului digital de elevație
• Clasificăm valoriile•
• Spatial Analyst Tools -> Zonal -> Zonal Histogram, cu o coloană care exprimă procentual prezența intervalului în raport cu totalul suprafeței, care pe urmă este intodus în formula de calcul:
(x * y) / 100 = z x = vpmax - vpmin
Adâncimea fragmentării 2
X Y Z bazaValoare
ProbabilitateVAL *
1000000.1 9.4139 0.0094139 0 0.0094139 941.390.2 12.12972 0.02425944 0.1 0.12425944 12425.944
0.19 39.7593 0.07554267 0.31 0.38554267 38554.2670.29 34.18009 0.099122261 0.51 0.609122261 60912.22610.2 4.516984 0.009033968 0.8 0.809033968 80903.3968
Calculul valorii de probabilitate
• Spatial Analyst Tools -> Reclass -> Reclassify (asocierea valorilor vechi cu valorile noi,valoarea de probabilitate)
• Spatial Analyst Tools -> Map Algebra -> Raster Calculator unde eliminăm înmulțirea cu 100000 introdusă în Excel.
• Datele se pot verifica prin încadrarea valorilor între 0 și 1.
Densitatea fragmentării • Analysis Tools -> Overlay -> Identity ( introducem
apele și caroiajul)
• Data Management Tools -> Generalization -> Dissolve ( separam apele după locul ocupat în caroiaj și calculăm lungimea și densitatea lor)
• Cu Join adăugăm aceste valori la caroiaj și îl transformăm în raster.
• Urmărim pașii menționați anterior cu precizarea că la reclasificare la no data la valorile noi să avem 0.
Densitatea fragmentării 2
• Se reduce dimensiunea pixelilor de la 1000 m la rezoluția modelului digital de elevație (25m)
Expozițiunea versanțiilor• Spatial Analyst Tools -> Surface -> Aspect
(MDE)
• Reclasificarea a fost efectuată de două ori. (E,NW)
Wetness Index – Flow direction, Flow
accumulation
• Ln((“Flow_accumulation”*400)/Tan(“Panta”))
• Spatial Analyst Tools -> Hydrology -> Flow Direction
• Spatial Analyst Tools -> Hydrology -> Flow Accumulation
• Spatial Analyst Tools -> Map Algebra -> Calculate Raster (cu formula sus menționată)
• La no data new values introducem 0. (Reclasificare)
Rezultatul
• Radical ((Panta * Expozițiunea versanțiilor)/7 * (Hipsometria + Adâncimea Fragmentării + Densitatea Fragmentării + Wetness Index + Stream Power Index))
• Spatial Analyst Tools -> Map Algebra -> Raster Calculator
• Reclasificare după intervalele date.