Slovenia CLSS LiDAR — Flood Susceptibility & Forest Health Map

Overview

This project turns a large raw LiDAR dataset from Slovenia’s official CLSS national scanning programme into a deployable flood-analysis product. I built both the offline analysis pipeline and the interactive delivery layer: processing region tiles, calibrating dataset-wide factor ranges, generating overlay rasters, ranking the highest-risk cells, and publishing the results as a static web map with no backend.

The portfolio value here is not only the visualization. The harder engineering work is in transforming roughly 15 GB of raw .laz data into lightweight web assets that remain spatially consistent, comparable across tiles, and cheap to host. The final system visualizes flood susceptibility, forest coverage, land classification, and ranked high-risk locations across a contiguous 9 × 9 km block over Ljubljana, with additional regional examples.

Official Dataset

The source data comes from Geodetska uprava RS through Slovenia’s CLSS programme, which makes this a project grounded in official government geospatial infrastructure rather than synthetic or toy inputs.

• Source: official Slovenia national airborne LiDAR programme
• Coordinate systems: EPSG:3794 and EPSG:4326 in the published map workflow
• Grid resolution: 2 m × 2 m × 2 m voxels in the derived analysis pipeline
• Raw format: .laz point-cloud tiles processed into deployable raster overlays

What The System Visualizes

• Flood susceptibility from a weighted composite of topographic wetness, canopy interception, NDVI health, plan curvature, and terrain roughness
• Forest coverage / NDVI health from per-cell vegetation scoring using 16-bit NIR and red channels
• Land classification separating ground, low vegetation, medium vegetation, high vegetation, and building returns
• High-risk locations from globally ranked candidate cells exported into interactive map markers

Visual Outputs

Flood Susceptibility

The flood layer is where the pipeline becomes more than a renderer. Each pixel score combines topographic wetness, vegetation condition, canopy interception, terrain roughness, and curvature into a normalized susceptibility surface. The map then exposes the highest-ranked candidate cells as clickable points so the result is not just pretty heatmapping, but actionable prioritization.

Forest Coverage (NDVI)

This view shows how the same processed data can be reinterpreted as a forest health / vegetation condition layer. Instead of flood susceptibility, the interface surfaces the gradient from stressed to healthy vegetation, which is useful both as a direct environmental signal and as an input into the flood model itself.

Land Classification

The land-classification layer is important because it proves the pipeline is not producing a single visualization in isolation. It extracts multiple semantically meaningful map products from the same official dataset, which makes the site feel like a lightweight geospatial analysis platform rather than a one-off demo.

What I Built

• A dataset-wide calibration stage so overlay scores remain comparable across all processed tiles
• A Python pipeline that transforms raw GKOT_*.laz files into deployable susceptibility, NDVI, and classification rasters
• A ranking/export layer that produces globally ranked high-risk cells and GeoJSON markers
• A static frontend that overlays the outputs interactively and publishes them through GitHub Pages

Technical Depth

The core technical problem was making flood-risk scores comparable across the entire dataset, not just visually plausible within a single tile. I solved that by introducing a calibration stage that derives dataset-wide normalization ranges once, stores them with a dataset fingerprint, and warns when the underlying tile set changes and recalibration is needed.

The canonical pipeline.py then processes all GKOT_*.laz files into deployable assets:

• per-tile susceptibility rasters
• per-tile NDVI rasters
• per-tile land-classification rasters
• a manifest consumed by the web app
• a globally ranked candidate list
• a GeoJSON file for top flood-risk points

That separation between raw source data and committed derived outputs matters. The heavyweight LiDAR tiles stay local and gitignored, while the smaller analysis products in web/data/ are versioned and deployed, which keeps the project reproducible without making the repository impractical.

Architecture

The system is intentionally simple at runtime:

• Python pipeline for acquisition, calibration, raster generation, candidate ranking, and GeoJSON export
• Static web frontend for rendering overlays and interaction
• GitHub Pages deployment from the web/ directory on every push

This is a strong example of choosing the right architecture for the problem. Instead of introducing unnecessary services, I kept the runtime fully static and moved complexity into the offline data-processing layer where it belongs.

Why It Matters

This project shows that I can work comfortably at the intersection of geospatial data engineering, analysis design, and product delivery. I did not just compute one metric on one file; I built a repeatable pipeline that scales across many official LiDAR tiles, preserves comparability at dataset scope, and ships the results in a form that is immediately inspectable by a user in the browser.

Tech Stack

• Data processing: Python, laspy, lazrs, numpy, scipy, pyproj, Pillow
• Geospatial outputs: PNG overlays, manifest-driven tile loading, GeoJSON risk points
• Frontend / deployment: Static web app, GitHub Pages