Le portail francophone de la géomatique

Bonjour Pascal,

nous utilisons Postgis Raster depuis un petit moment maintenant, mais pour des opérations basiques (ce sont nos besoins qui le sont).

Concrètement nous avons intégrer le MNT de la BD TOPO ainsi que des données de worldclim pour caractériser les stations d'espèces de papillons et de libellules (altitude, pente, exposition, météo). Les données d'observation arrivent régulièrement dans la base PostGIS et cela nous permet d'actualiser ces illustrations en temps réel, sans export ni traitement externe des données : http://si.cenlr.org/%C3%A9tiquettes/postgis-raster

On l'a aussi utilisé pour sélectionner au sein d'une couche d’occupation du sol tous les postes qui nous intéressaient entre deux altitudes.

Les performances sont très satisfaisantes. Je ne peux pas comparer avec les outils raster que tu site mais l’avantage indéniable est de pouvoir "tout faire en SQL".

Mon collègue affichait directement les rasters dans QGis.

Concernant la documentation, il y a la page de Pierre racine (http://trac.osgeo.org/postgis/wiki/WKTRaster), mais à compléter par la doc de PostGIS 2.0 (beaucoup d'évolutions sont apparues entre ces deux documents). PostGIS In Action y consacre son chapitre 13.

Enfin je n'arrive pas à croire que tu hésites!

Bonne journée,

Je vois que le sujet n'a pas franchement de succès !!
Alors j'ai décidé de faire quelques tests. Malgré la doc qui ne semble pas franchement fournie, j'ai pu faire quelques manipulations intéressantes et j'y vois un bon potentiel dans certaines chaines de traitement.

Les débuts ont été un peu fastidieux : je n'arrivai même pas à importer un TIF et à le visualiser dans QGIS ! Mais après un peu d'insistance et de recherche,  je pense désormais utiliser les fonctions Raster de Postgis dans mes prochaines analyses.

Côté performance, ca semble correct mais on est en dessous je pense de GRASS ou de Spatial Analyst d'ESRI. Il y a peut-être moyen d'optimiser mieux certaines requêtes mais je n'ai pas creusé.

Je pense qu'on peut dire Bravo aux développeurs de PostGIS Raster et on ne peut que les encourager à continuer !

Personnellement je suis très intéressé par l'utilisation de Postgis 2. mais je butte déjà dans l'import des données. l'import semble se réaliser, mais je n'ai qu'un tableau à 3 colonnes  avec le chemin des dalles. Quels sont les différents paramètres que vous avez utilisés avec le module raster2psql? ou dans la création de la base de données quel est le modèle que vous avez utilisé template_postgis20?

Bonjour,
j'ai importé des raster (232 tif de cadastre) dans une base postgis via raster2pgsql. Sans tuilage, le traitement a planté.
Donc j'ai spécifié une option de tuilage (chaque image a été divisée en 100), ainsi que l'option pour stocker les raster dans la base. Le traitement a duré (très) longtemps, mais ça a marché. marche.
Par contre, je constate que les requêtes sur les raster sont très lentes. Par exemple, un select qui ramène les images qui intersectent un polygone :

select r.rid, r.rast
from raster_lot9_1000 as r, zapm_feurs_en_forez_2 as p
where p.n_pm='74' and st_intersects(r.rast,p.geom)
;--715 lignes, 524643ms soit plus de 8 min !

Ensuite, je veux visualiser cette requête dans QGIS avec le dbManager, mais une erreur survient.
Je crée alors une table à partir de la requête précédente. Avec le dbManager, clic droit sur la table, ajouter au canevas, mais au bout de quelques minutes, QGIS plante.
Impossible de visualiser mes 715 raster dans QGIS !

Je cherche un retour d’expérience sur l'affichage de raster dans QGIS.

Alexis

Bonjour,

Je vais déplacer cette discussion sur le forum GeoBD plus à sa place.

Pour info, les fonctions Raster dans PostGIS n'ont pas été architecturé pour stocker des raster pour faire de l'analyse. Le stockage n'est pas performant et ne le sera pas (source : le designer du raster dans PostGIS).

Y.

Yves, tu voulais dire ceci ?

Dernière modification par MathieuB (Fri 28 February 2014 12:08)

Je confirme

Y.

Euh.. petite précision, la dernière phrase n'a rien à voir avec les rasters !
"PostGIS c'est vraiment de la bombe", c'est ma façon de dire que depuis que j'ai découvert postGIS, je trouve vraiment que c'est génial. Et c'est aussi une occasion de remercier tous les contributeurs de ce forum.
Il n'en reste pas moins que cette question de visualisation de raster est pour l'instant non résolue.
Alexis

Bonjour,
très bien. Je suis complètement largué sur ce sujet de tuilage et d'overview. Si vous aviez un exemple concret, cela m'aiderai grandement, car j'ai un peu de mal avec l'aide en ligne sur postgis raster, qui contient beaucoup de notions à comprendre, je trouve...
Ma problématique est classique:
J'ai constitué un jeu de données test : 1 seule image tiff de cadastre au 1:100, de taille 10000x10000 pixels. J'avoue que dans ce cas, je ne perçois pas l’intérêt du tuilage (car 1 tiff de cadastre couvre une petite surface) , mais si je ne spécifie pas d'option -t (ou si je spécifie -t  10000x10000 ce qui correspond à 1 seule image donc pas de tuilage), raster2pgsql plante dans le traitement sql (trop de pixels dans l'image d'origine ?).
J'en déduis donc que le tuilage est obligatoire (pouvez-vous me le confirmer ?).
Pour mon test, j'ai essayé une option de tuilage "raisonnable" : diviser en 4 mon image initiale ; à ce stade, je ne vois pas l’intérêt de découper 1 image de cadastre au 1:1000 en 1 grand nbre de morceaux; si mon raisonnement n'est pas bon, dites-le moi.

Au niveau des options, malgré la lecture de l'aide en ligne (difficile à digérer en partant de zéro sur le sujet), j'ai du mal. La plupart des exemples spécifient -I (Index gist) -C (create new table) -M (vacuum analyze), donc j'ai essayé avec ça. Au passage, Vacuum analyze, c'est vraiment utile ?
Ensuite, je n'ai pas compris l’intérêt (et les limites ?) de -Y (copy statements).
De même, beaucoup m'ont déconseillé de stocker les raster en base, mais je ne vois pas d'exemple avec -R (out-db) raster) ; est-ce que ça pose des pb pour la visualisation dans qgis ?
Bref, au final j'ai testé les options suivantes :  -I -C -M -l 2 -F. Ça crée bien la table et l'overview. Je peux afficher chacune via le db manager (au passage, les couleurs noir et blanc sont inversée).
Par contre, lorsque je crée une vue qui crée une géométrie polygone à la manière de l'index gist sur le raster (comme vous me l’avez suggéré),
cela crée des polygone mais je ne parviens pas à avoir le visuel du cadastre de mes raster. Donc au final, je n'arrive pas à exploiter votre conseil. Pouvez-vous me donner un exemple ?
Cordialement,
Alexis GAILLOT
Le traitement

Bonjour,

Effectivement, Postgis raster n'est pas fait pour le stockage en vue d'affichage mais uniquement pour du traitement sur les rasters (calcul de pentes par exemple).
Ca implique des choses au niveau organisation:

Un raster 10000x10000 non tuilé stocké dans postgis, ca correspond a une table avec une seule ligne contenant TOUT le raster, qui est enorme: Les bases de données n'aiment pas du tout cela (une seule enorme ligne): cela prend bcp de temps de lire cette ligne et d'en faire une image avant affichage.
Le tuilage est donc quasiment indispensable lorsqu'on stocke un raster dans postgis:

la table contenant le raster represente la couverture spatiale du raster. Chaque ligne de la table représente une partie du raster, ou tuile (comme chaque ligne d'une table postgis vectorielle des communes contient un polygone de commune):
Il faut mieux, pour la performance de traitement de postigs, avoir 10 000 lignes avec des petites tuiles.

Si votre raster sert uniquement a l'affichage, stockez-le sur le disque, dans un format compatible avec l'affichage (tiff avec overviews, ecw, raster découpé avec un index shape des tuiles, etc...)
Qgis se chargera tout seul de n'afficher que la partie du raster concernée par l'extension spatiale de la carte.

Nico

Bonsoir Alexis,

Pour de la visualition il est préférable soit de partager un espace disque avec tous les postes nécessitant un accès à ces raster. Ceux-ci seraient accessibles au format TIFF ou ECW (ou MrSID) car ces formats gèrent la récupération d'une partie des données. Les formats PNG, GIF et JPEG, quant à eux, ont l'indécences de ne pas donner cette possibilité, par conséquent meme pour une petite zone tu devras lire l'ensemble du fichier

Autre possibilité : tu met en place un flux WMS ce qui évite le partage d'un espace disque. Les flux WMS passe par des ports standards du web (le port 80 ou 8080) donc pas de soucis s'il y a des parefeux en place. L'avantage des flux WMS est qu'ils sont facilement mis en cache via le standard WMTS (tout deux gérer par QGIS).

Y.

Bonsoir Yves,
merci pour tes 2 solutions, sur lesquelles je reviens quelques instants :
a)partage d'un espace disque qui contient mes TIFF.
Ce format gère la récupération d'une partie des données. --> Comment ? Il faut utiliser des outils du style GDAL ?

b)Autre possibilité : tu met en place un flux WMS. L'avantage des flux WMS est qu'ils sont facilement mis en cache via le standard WMTS (tout deux gérer par QGIS). --> tu sous-entend l'installation d'un serveur carto ? quand tu dis QGIS, tu parles de QGIS Server ?
Sur ces sujets (flux WMS, serveur carto), je suis novice total ; il me manque quelques explications pour mettre en œuvre ton conseil.

Cordialement, Alexis GAILLOT