Classe RasterSet : chargement d'un jeu de données raster

[Dolite]

RasterSet

Classe : RasterSet

Fichier concerné : Raster.py

Le but est de charger des informations sur un jeu de données raster et de la manipuler. Il y aura deux constructeurs

Constructeur from_list

Prend en entrée un chemin (fichier ou objet) vers une liste et un SRS (string) : la liste en entrée contient des chemins (fichier ou objet) vers des images. La fonction va récupérer pour chaque image de la liste les informations suivantes (a priori fait par la classe Raster) :

• les caractéristiques des canaux : le nombres et leur format
• la bounding box
• la taille en pixel

Elle tiendra à jour la liste des caractéristiques distinctes rencontrées (si toutes les images ont les mêmes, un seul élément) et la bbox globale.

Constructeur from_descriptor

Juste à partir du descripteur qui est écrit par la méthode d'instance, recharge toutes les informations. Il faut ajouter un constructeur de Raster à partir des informations sans lire le fichier source (si plusieurs descripteur, respecter la logique mise en place dans Pyramid ou Layer : from_descriptor, from_path, from_parameters, from_other pour un clone...)

Attributs

• Liste de Raster
• la bbox globale (attention, en 4326 on veut les latitudes en premier, ce que GDAL retourne de toute manière)
• le tableau des caractéristiques des canaux
• le SRS

Méthode d'instance

• write_descriptor (chemin) : écrit un JSON qui consigne toutes les informations, globales et par image. Passer par une propriété serializable (voir Pyramid)
• optionnel : write_vrt : écrit le fichier XML en respectant le format VRT : uniquement si une seule caractéristique de canaux

<VRTDataset rasterXSize="310000" rasterYSize="450000">
  <GeoTransform>  9.8100000000000000e+05,  1.0000000000000028e-01,  0.0000000000000000e+00,  6.7550000000000000e+06,  0.0000000000000000e+00, -1.0000000000000028e-01</GeoTransform>
  <VRTRasterBand dataType="Byte" band="1">
    <SimpleSource>
      <SourceFilename relativeToVRT="1">BDP_14_1011_6720_010_L93.tif</SourceFilename>
      <SourceBand>1</SourceBand>
      <SourceProperties RasterXSize="10000" RasterYSize="10000" DataType="Byte" BlockXSize="10000" BlockYSize="16" />
      <SrcRect xOff="0" yOff="0" xSize="10000" ySize="10000" />
      <DstRect xOff="300000" yOff="350000" xSize="10000" ySize="10000" />
    </SimpleSource>
    <SimpleSource>
      <SourceFilename relativeToVRT="1">BDP_14_1011_6721_010_L93.tif</SourceFilename>
      <SourceBand>1</SourceBand>
      <SourceProperties RasterXSize="10000" RasterYSize="10000" DataType="Byte" BlockXSize="10000" BlockYSize="16" />
      <SrcRect xOff="0" yOff="0" xSize="10000" ySize="10000" />
      <DstRect xOff="300000" yOff="340000" xSize="10000" ySize="10000" />
    </SimpleSource>
  </VRTRasterBand>
</VRTDataset>

Enrichissement

• Documentation
• Tests unitaires

[XavDmz]

la bbox globale (attention, en 4326 on veut les latitudes en premier, ce que GDAL retourne de toute manière)

L'ordre des coordonnées dans la bbox globale en fonction du CRS risque d'être compliqué à gérer. Je n'ai pas encore trouvé comment avoir cette information grâce au module python osgeo.gdal.

Pour exemple, voici le retour de la fonction osgeo.gdal.Dataset.GetGeoTransform() selon la projection pour une même image :

• en EPSG:3857 (lon, lat) : (242800.0971, 49.98151770491799, 0.0, 6262874.0058, 0.0, -50.0396596747973)
• en EPSG:4326 (lat, lon) : (2.181110382, 0.00044899161311475426, 0.0, 48.929302475, 0.0, -0.0002958633284552825)

[XavDmz]

Suite du commentaire précédent.

Les fonctions osgeo.gdal.Dataset.GetProjection(), osgeo.gdal.Dataset.GetProjectionRef(), et osgeo.gdal.Dataset.GetSpatialRef() retournent toutes trois une chaîne de caractère contenant des informations similaires sur le SRS ou CRS de la donnée. Parmis ces informations se trouvent les descriptions des axes du référentiel, dans l'ordre des coordonnées. Ça peut être un moyen de connaître l'ordre attendu, mais ce n'est pas idéal.

Exemple de substring décrivant les axes en retour de GetProjection() :

• en EPSG:3857 (lon, lat) : AXIS["Easting",EAST],AXIS["Northing",NORTH]
• en EPSG:4326 (lat, lon) : AXIS["Latitude",NORTH],AXIS["Longitude",EAST]

Le hic, c'est que ces informations concernent le CRS isolé de l'image.
Comment utiliser les retours de GetGeoTransform() et GetProjection() pour savoir dans quel ordre exprimer la bbox ? Je ne sais pas.

[XavDmz]

J'ai une incertitude sur l'une des attributs de la classe, le tableau des caractéristiques des canaux.

Je vois ce tableau comme une liste dont chaque élement est un dictionnaire {'bands': int, 'format': ColorFormat}, mais je ne sais pas ce que doivent représenter les éléments dans le cas d'un set hétérogène.
Doivent-ils représenter :

• les caractéristiques de chaque image avec le même indice que dans la liste délements Raster ?
  • ça ferait doublon avec cette liste
• les différentes caractéristiques distinctes rencontrées dans le set ?
  • avec un mapping sous forme d'une troisième propriété qui liste les index des instances Raster concernées
  • sans aucun mapping, juste pour synthétiser les codages couleur utilisés dans le set

Dans tous ces cas, si le set est homogène la liste ne contient qu'un seul élément.
Peut-être que j'ai mal compris et qu'aucune de ces propositions ne convient, mais j'ai de toute façon besoin de précisions sur ce point.

EDIT:

Effectivement, en relisant:

Elle tiendra à jour la liste des caractéristiques distinctes rencontrées (si toutes les images ont les mêmes, un seul élément) et la bbox globale.

Donc le dernier cas dans mes propositions.