Le portail francophone de la géomatique

bonjour;
j'ai écrit un petit programme en idl "fenêtre glissante sur une image afin de regrouper les zones urbaines et autres" j'arrive a lire l'image sur laquelle je vais appliquer ce programme mais j'ai eu un problème pour faire tourner la boucle pour la création de masque.
le programme est le suivant :

pro FG,homog,moyenne,compt
;--------
; fonction : détection des zones urbaines
;
;
;-----------------------------------
;lecture du fichier
;-----------------------------------
dir_in='C:\Master\merguellil\SPOT\correction\coor_radio090414\corr_finale\'
dir_out='C:\Master\merguellil\SPOT\correction\coor_radio090414\corr_finale\'
pattern='*homog*.tif'
filename=file_search(dir_in,pattern,/fold_case,count=count_file)
basename=file_basename(filename)
envi_open_file,filename_homog
;-----------------------------------------------------------
;fenêtre glissante
;------------------------------------------------------------

siz_a=size(homog)
nl=5
nc=5


for i=5,nl-5 do begin
  for j=5,nc-5 do begin
  compt=0
  for p=i-4,i+4 do begin
     for q=j-4,j+4 do begin
           moyenne=moyenne+hom(i,j)
           compt=compt+1
      endfor
  endfor
       moyenne=moyenne/compt
  if moyenne>0.2 then begin
     hom(i,j)=1
   endif else begin
    hom(i,j)=0
     endelse
   endfor
  endfor
end

est que quelqu'un peut m'aider à le corriger??
merci