Makina Blog

Le blog Makina-corpus

Draper des lignes sur un MNT avec PostGIS


Comment obtenir des géométries 3D en drapant des géométries sur un MNT avec PostGIS 2

Dans Geotrek les géométries des tronçons sont saisies en 2D, puis drapées en 3D, afin d'obtenir des informations relatives à l'altimétrie (pente, dénivelées, …)

Nous présentons ici les fonctions PostGIS utilisées pour y parvenir.


Charger le MNT

Si votre Modèle Numérique de Terrain (MNT) est compatible avec GDAL, vous pouvez alors charger le raster très facilement dans la base de données avec les commandes suivantes.

Reprojette vers le SRID spécifié, rogne sur l'étendue spécifiée, et écrit le résultat dans un fichier binaire:

gdalwarp -t_srs EPSG:32632 -te 289942 4845809 400671 4947295 dem_file.geotif output.bin

Découpe en tuiles de 100 pixels et convertit en instructions SQL:

raster2pgsql -c -C -I -M -t 100x100 output.bin mnt > output.sql

Charge le SQL dans la base :

psql -d yourdb < output.sql


Draper les géométries

Il y a pour le moins 3 stratégies pour draper les géométries vers la 3D.

En gardant la résolution des géométries

On obtient une altitude pour chaque point de la ligne.

Avantages: On garde le niveau de définition original (nombre de points)

Inconvénients: On perd potentiellement beaucoup d'information ("sauts")

WITH line AS -- Prendre une ligne arbitraire (SELECT 'SRID=32632;LINESTRING (348595 4889225,352577 4887465,354784 4883841)'::geometry AS geom), points2d AS -- Extraire ses points (SELECT (ST_DumpPoints(geom)).geom AS geom FROM line), cells AS -- Obtenir l'altitude pour chacun d'eux (SELECT p.geom AS geom, ST_Value(mnt.rast, 1, p.geom) AS val FROM mnt, points2d p WHERE ST_Intersects(mnt.rast, p.geom)), -- Instancier des points 3D points3d AS (SELECT ST_SetSRID(ST_MakePoint(ST_X(geom), ST_Y(geom), val), 32632) AS geom FROM cells) - Contruire la ligne 3D depuis les points 3D ELECT ST_MakeLine(geom) FROM points3d;

Selon la résolution du MNT

On obtient une altitude pour chaque cellule du raster.

Avantages: On profite au maximum des informations du MNT

Inconvénients: On augmente énormément la résolution des géométries, parfois pour rien.

WITH line AS -- Prendre une ligne arbitraire (SELECT 'SRID=32632;LINESTRING (348595 4889225,352577 4887465,354784 4883841)'::geometry AS geom), cells AS -- Obtenir l'altitude pour chaque cellule intersectée (SELECT ST_Centroid((ST_Intersection(mnt.rast, line.geom)).geom) AS geom, (ST_Intersection(mnt.rast, line.geom)).val AS val FROM mnt, line WHERE ST_Intersects(mnt.rast, line.geom)), -- Instancier des points 3D, dans l'ordre de la ligne points3d AS (SELECT ST_SetSRID(ST_MakePoint(ST_X(cells.geom), ST_Y(cells.geom), val), 32632) AS geom FROM cells, line ORDER BY ST_Distance(ST_StartPoint(line.geom), cells.geom)) - Construire la ligne 3D depuis les points 3D ELECT ST_MakeLine(geom) FROM points3d;

En échantillonnant

On obtient une altitude par pas de X unités (mètres).

Avantages: On controle la résolution des géométries obtenues

Inconvénients: Parfois difficile un bon compromis selon la longueur des géométries

WITH line AS -- Prendre une ligne arbitraire (SELECT 'SRID=32632;LINESTRING (348595 4889225,352577 4887465,354784 4883841)'::geometry AS geom), linemesure AS -- Ajouter une mesure pour découper selon un pas (ici, de 50) (SELECT ST_AddMeasure(line.geom, 0, ST_Length(line.geom)) as linem, generate_series(0, ST_Length(line.geom)::int, 50) as i FROM line), points2d AS (SELECT ST_GeometryN(ST_LocateAlong(linem, i), 1) AS geom FROM linemesure), cells AS -- Obtenir l'altitude pour chacun d'eux (SELECT p.geom AS geom, ST_Value(mnt.rast, 1, p.geom) AS val FROM mnt, points2d p WHERE ST_Intersects(mnt.rast, p.geom)), -- Instancier des points 3D points3d AS (SELECT ST_SetSRID(ST_MakePoint(ST_X(geom), ST_Y(geom), val), 32632) AS geom FROM cells) - Construire la ligne 3D ELECT ST_MakeLine(geom) FROM points3d;

En tant que fonction PostgreSQL

Il est possible d'enfermer tout ce code dans une fonction :

CREATE OR REPLACE FUNCTION drape(line geometry) RETURNS geometry AS $$ DECLARE line3d geometry; BEGIN WITH … … … … SELECT ST_MakeLine(geom) INTO geom3d FROM points3d; RETURN geom3d; END; $$ LANGUAGE plpgsql;

Et draper les géométries d'une table:

-- Ajouter une colonne ALTER TABLE yourtable ADD COLUMN geom_3d geometry(LineStringZ, 32632);

-- geom_3d avec la version drapée UPDATE yourtable SET geom_3d = drape(geom);

Formations associées

Formations Outils et bases de données

Formation PostgreSQL

A distance (foad) Du 10 au 12 décembre 2025

Voir la Formation PostgreSQL

Formations SIG / Cartographie

Formation PostGIS

Toulouse Du 22 au 24 octobre 2025

Voir la Formation PostGIS

Formations SIG / Cartographie

Formation QGIS

Toulouse Du 1 au 3 octobre 2025

Voir la Formation QGIS

Actualités en lien

Geotrek et OpenS­treet­Map : Mise en place d’une passe­relle pour une connais­sance du terri­toire enri­chie

08/09/2025

Dans l’uni­vers des logi­ciels open-source, les plus belles inno­va­tions naissent souvent de la rencontre entre des commu­nau­tés qui partagent les mêmes valeurs. Aujour­d’hui, nous célé­brons une avan­cée majeure pour Geotrek : la créa­tion d’une passe­relle avec OpenS­treet­Map (OSM), la plus grande base de données carto­gra­phique colla­bo­ra­tive au monde. Plus qu’une simple fonc­tion­na­lité, ce projet est le fruit d’un travail d’in­gé­nie­rie et de recherche appro­fondi.
Voir l'article
Image
Logo d'illustration pour la passerelle entre OSM et Geotrek

Instal­ler Geotrek : avec ou sans segmen­ta­tion dyna­mique ?

08/09/2025

Geotrek-admin propose deux modes de fonc­tion­ne­ment pour gérer les objets liés aux tronçons : avec ou sans segmen­ta­tion dyna­mique. Ce choix a un impact impor­tant sur la manière dont sont stockées et gérées les données, et sur les possi­bi­li­tés d’édi­tion, de cohé­rence topo­lo­gique et d’in­ter­opé­ra­bi­lité avec d’autres systèmes. Dans cet article, on vous explique ce qu’est la segmen­ta­tion dyna­mique ainsi que le réfé­ren­ce­ment linéaire, ses avan­tages, ses limites, et dans quels cas il est perti­nent (ou non) de les utili­ser.
Voir l'article
Image
Réseau de tronçons dans Geotrek

Geotrek-Mobile, zones sensibles désor­mais acces­sibles

17/07/2025

Une nouvelle fonc­tion­na­lité vient enri­chir l’ap­pli­ca­tion Geotrek-Mobile : l’af­fi­chage des zones sensibles direc­te­ment sur la carte, le long des itiné­raires.
Voir l'article
Image
Geotrek-Mobile contenus outdoor

Inscription à la newsletter

Nous vous avons convaincus