Final Assignment for Python in GIS
- gdal version 1.11.3 (install:
pip install gdal
) - tested on Python 2.7.12
- Download or clone this repository and if necessary unzip all files.
- Download Landcover image (this script has been tested with and optimized for g100_clc12_V18_5.tif)
- Download DEM (this script has been tested with and optimized for
dgm1_5meter.img
. This file was provided within the course)
Note: Two of three of the provided shapefiles are missing a projection. Therefore the .prj
of the third shapefile needs to be duplicated twice and renamed according to the filenames of the other shapefiles.
- Open evaluater.py in a text editor
- Add paths to input and output files in lines 23-37 (description of allowed types and variable definitions is given within evaluater.py)
- Open IDLE
- Open
evaluater.py
- Click on
Run > Run Module
The script creates three outputs:
- reprojected shapefile with CRS of DEM and added field
DEM_elev
- reprojected shapefile with CRS of Landcover file and added field
Land_cov
- png file containing
- bar chart of mean elevation differences between GPS and DEM elevation, classified in landcover types
- time series of elevation differences between GPS and DEM elevation, classified in landcover types
Example .png outputs have been added to this repository