PostGIS : Gestion de données raster

Se rendre dans pgAdmin et afficher le « contenu » de la vue « raster_columns » qui doit contenir une seule ligne, représentant le raster importé.

Cette vue est l'équivalent de la vue « geometry_columns » pour les rasters.

Pour l'instant, la plupart des colonnes sont vides concernant le raster importé, ce qui empêche son exploitation avec des outils externes tels que gdalinfo ou QGIS. De plus, l'identifiant de référence spatiale (SRID) n'a pas été précisé lors de l'import et n'est donc pas renseigné, ni dans les contraintes, ni dans le raster lui-même (colonne « rast » de la table mnt_25m_asc).

La fonction UpdateRasterSRID met à jour le SRID. Le système de projection du MNT est le Lambert93, qui porte le numéro 2154 d'après la table spatial_ref_sys.

Dans la fenêtre SQL de pgAdmin, exécuter :

À présent, il faut ajouter toutes les contraintes pour rendre le Raster exploitable par les outils externes. La fonction « AddRasterConstraints » fait cela automatiquement. Exécuter :

Retourner dans la vue « raster_columns » et rafraîchissons-la (F5). On obtient les valeurs suivantes :

• r_table_catalog = « formation » : c'est nom de la base de données

• r_table_schema = « public » : c'est le nom du schéma dans lequel la table se trouve

• r_table_name = « mnt_25m » : nom de la table

• r_raster_column = « raster » : nom de la colonne portant le raster

• srid = 2154 : identifiant de référence satiale du raster

• scale_x = 25 et scale_y = -25 : taille du pixel en coordonnées de projection géographique (largeur et hauteur, respectivement). La hauteur est négative car ici comme pour de nombreux rasters, la coordonnée Y en pixel évolue positivement de haut en bas quand la coordonnée Y géographique évolue positivement de l'origine vers le nord

• blocksize_x = 401 et blocksize_y = 401 : taille de chaque tuile, sujet que nous aborderons plus loin via la question du tuilage

• same_alignment = TRUE

• regular_blocking = FALSE

• num_bands = 1 : il n'y a qu'une seule bande car il s'agit d'un MNT simple

• pixel_types : comme chaque pixel est codé (2 bits, 4 bits, 8 bits, nombre à virgule flottante, etc., pour chaque bande

• nodata_values : valeurs interprétées comme « pas de données » (no data) pour chaque bande

• out_db : indique si le raster est stocké en interne ou en fichier externe, pour chaque bande (qui ont toutes la même valeur, TRUE ou FALSE, pour un raster stocké)

• extent : géométrie polygone représentant l'enveloppe géographique du raster