IDL define (divide) valores de píxel por encima del rango de acuerdo con las condiciones (artificiales)

Método de píxel por píxel (problemático, que consume mucho espacio):

inData = read_tiff(inputfile,GEOTIFF=Geovar,INTERLEAVE=2)
iSize = size(inData,/DIMENSIONS)
iNs = iSize[0]
iNl = iSize[1]
outData = fltarr(iNs,iNl)
for ii = 0,iNs-1 do begin
  for ij = 0,iNl-1 do begin
    if outData[ii,ij]  le 0  then begin
      outData[ii,ij] = 1
    endif else begin
      outData[ii,ij] = 2
    endelse
  endfor
endfor

Método de índice:

bgIdx=WHERE(inData le 0,bgCnt)
IF bgCnt GT 0 THEN BEGIN
  inData[bgIdx]=!VALUES.F_NAN 
ENDIF