Géoréférencer et projeter
Contrairement aux géométries, les rasters sont des grilles de points et n'ont pas de coordonnées géographiques explicites pour chacun des points : c'est l'ordonnancement régulier des points qui permet de faire l'équivalence entre un point du raster et son emplacement géographique, grâce aux informations de géoréférencement.
Le géoréférencement définit les coordonnées géographiques de l'origine des pixels ainsi que la dimension géographique du pixel (identique pour tous les pixels). Un paramètre de transvection X et Y est également possible (« skewX » et « skewY » en Anglais), mais nous ne l'aborderons pas ici.
La fonction « ST_GeoReference » permet de récupérer les informations de géoréférencement du raster, au format World File [http://fr.wikipedia.org/wiki/World_file] selon 2 variantes : GDAL ou ESRI.
Essayons :
SELECT ST_GeoReference( ST_Union( rast ), 'GDAL' ) FROM mnt_25m_asc;
Le format de géoréférencement GDAL se structure en 6 lignes, ainsi :
largeur du pixel en unités géographiques
paramètre de transvection Y
paramètre de transvection X
hauteur du pixel en unités géographiques (nombre négatif, le plus souvent)
coordonnée X du coin supérieur gauche
coordonnée Y du coin supérieur gauche
Pour davantage d'informations, nous nous référerons à la documentation : http://postgis.net/docs/RT_ST_GeoReference.html
Taille géographique du pixel
Taille géographique du pixel
Les fonctions ST_PixelWidth et ST_PixelHeight calculent respectivement la hauteur et la largeur du pixel en coordonnées géographiques (selon le référentiel spatial – SRID – du raster).
Exemples :
SELECT ST_PixelWidth( ST_Union( rast ) ) FROM mnt_25m_asc;
SELECT ST_PixelHeight( ST_Union( rast ) ) FROM mnt_25m_asc;
Projection
Pour obtenir l'identifiant de référence spatiale, et par là la projection utilisée, c'est la fonction « ST_SRID » comme pour les vecteurs :
SELECT ST_SRID( ST_Union( rast ) ) FROM mnt_25m_asc;
