Le portail francophone de la géomatique

Bonjour à tous,
je vois que mon souci n'attire pas les foules, soit le sujet est inintéressant, soit je me suis mal exprimé et personne n'a compris c'est quoi le problème.
je reprends:
je souhait extraire l'information (indice de végétation) contenu dans des pixels de manière rapide et efficace, j'ai les coordonnées géographique des ces pixels. pour l'instant j'utilise une méthode manuel avec le "corsor locator"

âme charitable, je ne demande pas le gite et le couvert, juste une solution à mon souci de télédétection.
merci

...ou c'est juste le mois d'août

Bonjour
Je pense que cette tache est faisable sur ENVI
Mais d'aprés mes connaissances, tu peut utiliser un logiciel SIG ( Arcgis or Qgis) qui supporte chargement des images satellites en un nombre de bandes et tu peux ressortir les valeurs des pixels pour lequel tu a les points.
Avant afin de s'assurer que les points se supporposent avec l'image sat que tu as; tu essaye de cré un fichier shp avec ces coordonnées et de le transformer sous format EVF ouvrable sous ENVI.
L'import de ce fichier te permet de s'assurer que le fichier est exact et il respond à la question de

Sous un logiciel SIG (Arcgis) tu ouvre l'image dans la bande que tu veux et le fichier shape et tu exploite la fonction extract values to points

J'espère que cette gymnastique de taches  va t'aidé

Bonjour
Il y a un outil sous Qgis qui permet de faire cette tache et qui est point sampling tools
https://pvanb.wordpress.com/2010/02/15/ … s-in-qgis/
CDT

Dernière modification par hichamfstgm (Thu 03 September 2015 10:45)

Bonjour

Quand tu essaye d'importer le fichier shp , ne spécifié pas le format du fichier cad le format ".EVF".
CDT

Bonjour babyshoone2,

Une fois que tu as créé ton shape, tu peux l'ouvrir avec ENVI. Comme l'a expliqué hichamfstgm, tu ne spécifies pas le type de fichier que tu souhaites ouvrir. Je ne sais pas quelle version d'ENVI tu utilises, mais soit il t'ouvrira directement le shape (c'est ce que fait ma version 5.1), soit il t'aidera à le convertir en .evf (c'est ce que fait la version 4.7 si mes souvenirs sont bons).

En ce qui concerne le Point Sampling Tool, il faut que tes 2 couches (Raster et vecteur) aient le même système de coordonnées. Mais attention, il ne s'agit pas simplement de changer le SRC de la couche vecteur dans Qgis, il faut que le shape soit créé d'ores-et déjà avec la bonne projection.
J'ai fait un test, j'avais également des valeurs nulles dans ma nouvelle table shape. J'ai changé le SRC de la couche vecteur pour qu'elle soit la même que mon image, résultat idem, valeurs nulles. J'ai recréé le chape avec la bonne projection et ça fonctionnait.

Bon courage,

Bénédicte

Bonjour les amies,

c'est formidable, merci, je vais de suit essayer.
pour info je suis parti faire un tour dans un forum anglophone, ils mon proposé d'utiliser IDL, et mon envoyer des "cript" (c'est comme ça qu'on dit ?) enfin, j'ai rien compris je ne m' y connais pas en script. je vous envoie leur répons si quelqu'un peut nous expliquer:

on ENVI & IDL is easy run the script roi_to_csv.pro.... And is great for many  layers stack!!!!



pro rois_to_csv
  ;Select file name to save
filename=dialog_pickfile(Title='Enter File Name')
if filename eq '' then return
Openw, lun, filename, /Get_lun
  envi_select, title='Input Filename', fid=fid, pos=pos
  if (fid eq -1) then return
  envi_file_query,fid,bname=bname,nb=nb
  bname=["ROI_Name",bname]
  printf,FORMAT='(1000(A,:,","))',lun, bname
  ;Check ROIs associated with the image
  roi_ids = envi_get_roi_ids(fid=fid, roi_names=roi_names,/Short_name)
  if (roi_ids[0] eq -1) then begin
    print, 'No regions associated with the selected file'
    return
  endif
  ; Compound widget for ROI selection
  base = widget_auto_base(title='ROI Selection')
  wm   = widget_multi(base, list=roi_names, uvalue='list', /auto)
  result = auto_wid_mng(base)
  if (result.accept eq 0) then return
  ptr = where(result.list eq 1, count)
  ; ROI data extraction
  for i=0l, count-1 do begin
     ;Initialize temp_data to hold the ROI's data
      data = envi_get_roi_data(roi_ids[ptr[i]], fid=fid, $
      pos=[0])
      temp_data = fltarr(n_elements(pos),n_elements(data))
      ;Extract the data for each ROI
       for j=0l, n_elements(pos)-1 do begin
               data = envi_get_roi_data(roi_ids[ptr[i]], fid=fid, $
               pos=pos[j])
               temp_data[j,*] = data
               help, temp_data
           endfor
       ;Prints the ROI data to filename
           for k = 0L, n_elements(temp_data[0,*])-1 do begin
               printf, FORMAT='(A,:,",",1000(F12.2,:,","))',lun, roi_names[i], temp_data[*,k]
           endfor
       ;Returns to the next ROI
  endfor
   Free_lun, lun
   print,'Done!'
end

Rebonjour,

ouf ça marche, je sais pas comment mais ça marche, au fait j'ai fait le contraire de ce que Bénédicte a suggéré, j'ai changé le CSR du shape et non créer un nouveau shape avec le nouveau csr pour que ça marche.

Maintenant j'aimerais que l'on revienne au script de nos amis les "Zenglishe" sous IDL si vous le voulez bien, si quelqu'un s'y connait il serait formidable de nous orienter quelque peu.

merci à tous