Integration of advanced GIS technologies for decision support

Introduction

Data

Remote sensing data

Free satelite images: Landsat, Sentinel browsers:

The Copernicus Open Access Hub (previously known as Sentinels Scientific Data Hub) provides complete, free and open access to Sentinel-1, Sentinel-2, Sentinel-3 and Sentinel-5P user products

Landsat, SRTM and much more data provider:

GIS data

Urban Atlas

CORINE

GEOPORTAL PL

Open Street Map OSM

Data preparation

Coordinate System (CS) of the features

XYZ Tiles in QGIS - base maps are available

NO projection, only coordinate on the elipsoid! Problem!!!

CS of the project

QGIS is possible to change CS on-the-flight

CS EPSG codes

rys1

Change CS of the feature

We would like to avoid geographicala CS. In QGIS you can transform the data coordinates during export !!!!

right klick on the feature in legend-Export-Save feature as

you can:

if you clip raster you can change spatial resolution (pixel size), but take care QGIS cheats, check in feature properties and correct coordinates of map extent

Changing CS:

Vector analysis

Vector analysis - case study - site location

conditions:

  1. We are looking "barriers=gates"
  2. 10 m around residential streets

How we realize it?

  1. Query by attribute: SQL (osmpoints) SELECT * FROM osmpoints WHERE BARRIER='GATE'
  2. Query by attribute: SQL (osmlines) SELECT * osmlines WHERE HIGHWAY='RESIDENTIAL'
  3. Buffering Vector-Geoprocessing-Buffer (selected feature only), 10, round, dissolved result = condition1
  4. Overlaying 1 by 3 Vector-Geoprocessing-Intersection osm_points selected overlay on canditon1=condition2

Homework1:

  1. Install: QGIS, SAGA
  2. Download OSM and SRTM
  3. Perform vector analysis like above

Raster analysis

An example

Attribute querry:

  1. Show me the area below 205m

  2. Show me the area between 250m and 300m

old value1 old value2 new value

0 250 0

250 300 1

300 10000 0


Raster terrain analysis

Slope and aspect calculation

Slope Algorithms

Horn, B.K.P., 1981. Hill shading and the reflectance map. Proceedings of the IEEE 69:14-47

ESRI tutorial

How Slope works

How Hillshade works

How Aspect works

Raster analysis - case study - site location

We are looking for the steep area with south exposure

conditions:

  1. slope >15 degree
  2. aspect [135;225]

How we realize it?

  1. Slope calculation
  2. Query by attribute: slope>15
  3. Aspect calculation
  4. Query by attribute: aspect>=135 or <225
  5. Overlay, raster calculator, intesection

Common area, location that meets both conditions

multiplication of both raster layers (only 1 and 1 gives satisfactory result)

0=0*0

0=0*1

0=1*0

1=1*1


Digital Terrain Model

nmt

Slope

slope

Steep slopes

slope1

Aspect

aspect

South facing slopes

aspect1

Steep south facing slopes

wynik1

Activated Zevenbergen&Thorne formula for smooth landscapes

wynik1_ver2

Steep south facing slopes on Open Street Map

wynik1_osm