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

Localizarea zonei de interes

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.

Adâncimea fragmentării 3

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)

Panta

• Spatial Analyst Tools -> Surface -> Slope (MDE)

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)

Wetness Index – Flow direction, Flow

accumulation

Stream Power Index• Ln((“Flow Accumulation”+0.001)* (“Panta”/100

+0.001))

Hipsometria• Reclasificare mde, inversarea unor intrări în Excel.

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.

Rezultatul

Validare și Concluzii