Le portail francophone de la géomatique

Bonjour,

par métadonnées, je pense que vous voulez dire sans perdre le géo-ref du raster. Je n'ai pas ENVI avec moi, mais tu peux découper le raster via l'outil layer stacking dans le menu tools ( si je me souviens bien). vous sélectionnez les 4 bandes SPOT et vous pouvez logiquement lui indiquer un shape pour la découpe.

En fait pour faire un découpage selon les limites définies par un shape sous ENVI, il faut d'abord ouvrir le shape dans ENVI (open vector file) pour qu'il le convertisse au format evf. Ensuite c'est effectivement basic tools> Resize data ou layer stacking, puis choisir l'image que vous voulez découper. Il faut également cliquer sur Spatial Subset et choisir ROI/EVF et là vous pouvez sélectionner votre shape pour la découpe.

Je ne connais pas bien l'outil feature extraction de ENVI,  vous ne pouvez pas faire une segmentation forcé par le shape, classez les zones d'études puis resegmentez les zones classées?

C'est possible d'automatiser cela avec ENVI-IDL mais ça suppose de programmer en IDL, sinon si le nombre d'îlots est raisonnable, vous pouvez sous un logiciel de  SIG comme QGIS  sélectionner chaque entité et l'enregistrer comme nouvelle entité, ainsi vous aurez un shape par ilôts, il suffit ensuite d'importer chacun des îlots et de faire le découpage comme mentionné précedement.

Il faut aussi ajouter que sous R ce genre de manip est également possible (package raster, gdal) avec la fonction crop.

Pour conserver les métadonnées il faut utiliser une variable qui va stocker ces informations avec la fonction ENVI_GET_MAP_INFO de ce type :

map_info=ENVI_GET_MAP_INFO(fid=fid,UNDEFINED=u)

Et le spécifier dans l'écriture du fichier de sortie

ENVI_WRITE_ENVI_FILE,file,map_info=map_info