D - SCHNEIDER, MARIO JORGE.pdf - Universidade Federal do ...

; Ajuste das faixas dos dados
ptr = where(hue lt 0, count)
if(count gt 0) then hue(ptr) = hue(ptr) + 360.
ptr = where(sat lt .000001, count)
if (count gt 0)then hue(ptr) = 400.
Trecho de código para conversão HSV para RGB
;Constantes do usuário
k2 = sqrt(2) / 2.0
k6 = sqrt(6) / 6.0
k3 = sqrt(3) / 3.0
k7 = sqrt(6) / 3 .0
hue=float(envi_get _tile(tile_id 1 ,i,ys-ys,ye=ye))
sat=float(envi_get_tile(tile_id2,i))
val=float(envi_get_tile(tile_id3,i))
; Calcula as coordenada RGB
hue = temporary(hue) * (ípi/180.)
bl = sat * cos(hue)
xl = sat * sin(hue)
red = float(k3 * vai - k2 * xl - k6 * bl)
grn = float(k2 * xl + k3 * vai - k6 * bl )
blu = float(k7 * bl + k3 * vai)
; Clip to 0 and 255 and convert to byte
red = fix((red > 0) < 2048)
grn - fix((gm > 0) < 2048)
blu = fix((blu > 0) < 2048)
Trecho de código para fusão HSI
sharpen_doit, fid=[fid, fid, fid], pos=pos, f_fid=f_fid, method=0, $
f_pos=f_pos, f_dims=f_dims, interp=0, out_name=out_name, $
r_fid=r_fid
;end
envi_stats_doit, fid=r_fid, pos=pos, dims=dims, dmin=dmin, $
dmax=dmax, mean=mean, stdv=stdv, cov=cov, evec=evec, eval=eval, $
comp_flag=4, report_flag=5, to_screen=l, rep_name=rep_out_name