Vježbe 8 Analiza rastera

Transcription

1 Vježbe 8 Analiza rastera Sadržaj današnjih vježbi: a) Primjena filtera b) Analiza nagiba padine c) Reklasifikacija rastera d) Clip e) Cost distance analiza f) Buffer analiza g) Ekstrakcija rasterskih vrijednosti 2 1

2 Što na kraju morate predati: 3 1. Pokrenite novi projekt u QGIS-u i dodajte podatke: AsterDEM_Vj8.tif PotencijalneLokacije.shp Naselja.shp *U sljedećim koracima cilj će vam biti izraditi procjenu potencijalnih lokacija za izgradnju odlagališta otpada (PotencialneLokacije.shp). Glavni uvjeti za izbor najbolje lokacije su udaljenost od naselja i cijena izgradnje i korištenja prilazne ceste do odlagališta. UVJETI ODABIRA LOKACIJE: 1) Odlagalište NE smije biti smješteno bliže od 1.5 km do bilo kojeg naselja (Naselja.shp). 2) Prilazna cesta se NE smije graditi na nagibima terena iznad 10. 3) Cijena izgradnje i vožnje kamiona po prilaznoj cesti bi trebao biti što manji. Trebala bi izbjegavati strmije dijelove terena! Glavni utjecaj na cijenu izgradnje pojedinog segmenta ceste u ovom primjeru biti će nagib padine. Taj utjecaj procjeniti ćete primjenom cost distance analize na temelju digitalnog modela reljefa (AsterDEM_Vj8.tif). 2

3 Novigradsko more Obrovac Zrmanja Karinsko more *Artefakte možemo ukloniti prethodnom primjenom filtera za zaglađivanje površine (Low pass filter) na izvorne DEM podatke 2. U Windows-ima pokrenite Notepad (Start>Programs>Accessories) i upišite tekst točno kako je prikazano na slici: Spremite tekst pod imenom: LowPassFilter3x3.txt (u isti direktorij sa ostalim podacima ) Sadržaj filter datoteke: Oblik matrice 3x3 TAB! Vrijednost sa kojom se vrijednosti matrice dijele 6 3

4 3. U Processing Toolbox-u pod GRASS commands>raster pronađite i duplim klikom pokrenite alat r.mfilter. Unesite sljedeće postavke: Input layer: AsterDEM_Vj8 Filter file: (*kliknite browse) LowPassFilter3x3 Number of times to repeat the filter: 15 Apply filter only to zero data values: Isključiti Filtered: Save to file...>asterdem_filtriran (*Browse: spremite u isti direktorij sa ostalim podacima kao.tif) Isključite Open output file after running algorithm Ostale postavke ostavite kakve jesu i kliknite Run 4. Dodajte podatke u projekt: AsterDEM_filtriran.tif 7 Prije filtriranja (AsterDEM_VJ8): Poslije filtriranja (AsterDEM_filtriran): Rezolucija rastera je ista! 4

5 *Ponovite izradu karte nagiba ovaj puta na filtriranim podacima 5. Pod Raster>Terrain Analysis pokrenite alat Slope. Unesite sljedeće postavke: Elevation layer: AsterDEM_filtriran Output layer: AsterDEM_slope Output format: GeoTIFF Z factor: 1.0 Uključite Add result to project 9 Rezultat - karta nagiba padina: (svijetle nijanse=veći nagib) Ako prikažemo puni raspon podataka raster, vidimo da se nagib padine na promatranom području kreće od 0 do oko 43 5

6 *Da nismo prethodno primijenili low pass filter rezultat bi izgledao ovako: Rezultirajuća karta nagiba sadržavala bi artefakte (vodoravne i okomite linije ) koji su rezultat artefakata u izvornim podacima digitalnog modela reljefa (bilo zbog metode prikupljanja ili kasnije obrade podataka) *Izvorni DEM: *U sljedećem koraku ćete reklasificirati kartu nagiba tako da raspon vrijednosti od 0 do 43 svedete u zadane klase, te svakoj klasi prododate određenu cijenu (Cost) koja će odražavati faktor porasta cijene izgradnje i vožnje po terenu zbog povećanog nagiba (veći nagib - > veća cijena). Prethodno je potrebno izradtiti reklasifikacijsku tablicu. 6. U Windows-ima pokrenite Notepad (Start>Programs>Accessories) i upišite tekst točno kako je prikazano na slici: Čelije u rezultirajućem rasteru će imati jednu od 7 mogućih vrijednosti (1, 2, 4, 8, 16, 32, NULL) Spremite tekst pod imenom: ReklasifikacijskaTablica (u isti direktorij sa ostalim podacima ) Cijena raste eksponencijalno sa svakom sljedećom klasom nagiba Nagib veći od 10 dobiti će NULL vrijednost i neće ulaziti u daljnju analizu 12 6

7 7. Pod GRASS commands>raster pronađite i pokrenite alat r.reclass. Unesite sljedeće postavke: Input raster layer: AsterDEM_slope File containing reclass rules: ReklasifikacijskaTablica.txt Output raster layer: Save to file...>nagib_reklasificiran (*Browse: spremite u isti direktorij sa ostalim podacima kao.tif) Isključite Open output file after running algorithm Ostale postavke ostavite kakve jesu i kliknite Run 8. Dodajte podatke u projekt: Nagib_Reklasificiran.tif 13 7

8 *Stil rastera Nagib_Reklasificiran može se podesiti tako da se svaka od 6 (7) klasa prikaže zasebnom bojom (ne utječe na tijek daljnje analizu) *Iz reklasificirane karte nagiba želimo ukloniti područja mora i rijeka koja fizički ne mogu biti lokacije odlagališta niti želimo da niski (pogodan) nagib tih područja utječe na analizu 9. Dodajte podatke u projekt: Poligon_BezVodenihPovrsina.shp 10. Pod menijem Raster>Extraction pokrenite alat Clipper i podesite sljedeće postavke: Input file: Nagib_Reklasificiran Output file: Nagib_Reklasificiran_Clip (*Browse: spremite u isti direktorij sa ostalim podacima kao.tif) Uključite No data value: 0 Clipping mode: Mask layer Mask layer: Poligon_BezVodenihPovrsina Uključite Load into canvas when finished 16 8

9 Nagib_Reklasificiran_Clip.tif (Stil podešen za prikaz pojedinih klasa): Područje mora i rijeke Zrmanje su izrezani iz rastera 11. Dodajte podatke u projekt: Ceste.shp 12. Pod SAGA>Vector point tools pronađite i pokrenite alat Convert lines to points. Unesite sljedeće postavke: Lines: Ceste Insert Additional Points: Uključite Insert distance: 100 Points: Ceste_Tocke (*Browse: spremite u isti direktorij sa ostalim podacima kao.shp) Uključite Open output file after running algorithm Kliknite Run 18 9

10 13. Budući da prethodni alat ne preuzima koordinatni sustav, novostvoreni shapefile je potrebno ponovno spremiti sa odgovarajućom projekcijom: Desni klik na sloj>spremi kao... Spremi kao: Ceste_TockeHTRS (*Browse: spremite u isti direktorij sa ostalim podacima kao.shp) CRS: Selected CRS (Traži> HTRS96/Croatia TM) Uključite Add saved file to map Pod GRASS commands>raster pronađite i pokrenite alat r.cost.points (ili r.cost.full). Unesite sljedeće postavke: Unit cost layer: Nagib_Reklasificiran_Clip Start points: Ceste_TockeHTRS Use the Knight s move : Yes Keep null values in output raster layer: Yes Cumulative cost: Save to file...>procjenacijene (*Browse: spremite u isti direktorij sa ostalim podacima kao.tif) Isključite Open output file after running algorithm Ostale postavke ostavite kakve jesu i kliknite Run 15. Dodajte podatke u projekt: ProcjenaCijene.tif 20 10

11 ProcjenaCijene.tif (Stil podešen za prikaz): 16. Pod menijem Vector>Alati Geoprocesiranja pokrenite alat Zaštita/e (fixed distance buffer) i podesite postavke: Unesi vektorski sloj: Naselja Udaljenost međuprostora (buffer): 1500 Izlazna Shape datoteka: Naselja_Buffer (*Traži: spremite u isti direktorij sa ostalim podacima kao.shp) Uključite Dodaj rezultat na prikaz karte 22 11

12 Pogodne lokacije ali su pre blizu naselja 17. Pod menijem Vector>Alati Geoprocesiranja pokrenite alat Razlika i podesite postavke: Unesi vektorski sloj: PotencijalneLokacije Sloj razlike: Naselja_Buffer Izlazna Shape datoteka: PotencijalneLokacije_Izvan1500m (*Traži: spremite u isti direktorij sa ostalim podacima kao.shp) Uključite Dodaj rezultat na prikaz karte 24 12

13 18. Pod SAGA>Vector to raster pronađite i pokrenite alat Add raster values to points. Unesite sljedeće postavke: Points: PotencijalneLokacije_Izvan1500m Grids: ProcjenaCijene Interpolation: Nearest Neighbor Result: PotencijalneLokacije_Cijena (*Browse: spremite u isti direktorij sa ostalim podacima kao.shp) Uključite Open output file after running algorithm Kliknite Run Podesite Stil sloja potencijalnih lokacija (PotencijalneLokacije_Cijena) kategorizirano prema vrijednosti kolone PROCJENACI isto (ili slično) kao na slici 26 13

14 20. Podesite prikaz prema slici dolje (ProcjenaCijene.tif, Naselja.shp, PotencijalneLokacije_Cijena.shp) spremite ga kao.jpg sliku pomoću opcije Project>Spremi kao sliku... Sliku nazovite RPUG_Vjezbe8_Ime_Prezime.jpg: 27 Sliku (.jpg) pošaljite kao attachment na: 28 14