math_doit split window [message #38315] |
Sat, 28 February 2004 14:20 |
jmartine
Messages: 1 Registered: February 2004
|
Junior Member |
|
|
Hi,
I�m new in IDL and this is my first post.
I�m trying to perform a split window algorithm on a set of avhrr
images, and I'm doing with the envi math_doit routine. the problem is
that the algorithm
needs inputs from differents files and I don't know how put bands from
differents files in the expresion
exp = 'b4+(2+(0.28*W))*(b4-b5)-(0.4-(0.48*W))+(53-(4*W))*(1-e)+(14 9-26W)De
pro split_window
envi, /restore_base_save_files
envi_batch_init, log_file='batch.txt'
;
envi_open_file, '1995029_cal', r_fid=fid
if (fid eq -1) then begin
envi_batch_exit
return
endif
;
envi_file_query, fid, ns=ns, nl=nl, nb=nb
t_fid = [fid,fid]
dims = [-1, 0, ns-1, 0, nl-1]
pos = [0,1]
exp = 'b4+(2+(0.28*W))*(b4-b5)-(0.4-(0.48*W))+(53-(4*W))*(1-e)+(14 9-26W)De
out_name = '1995029_ts'
;
envi_doit, 'math_doit', $
fid=t_fid, pos=pos, dims=dims, $
exp=exp, out_name=out_name, $
r_fid=r_fid
;
; Exit ENVI
;
envi_batch_exit
end
for example this are the three files for an image
1995029_cal
1995029_emi
1995029_deltae
the first file is the avhrr calibrated image (5 bands), the two last
images has only one band, corresponding to the "e" and "De" parameters
of the expresion.
W is a constant.
Any help is appreciated
Thanks in advance
Javier Mart�nez
Instituto del Medio Ambiente
Universidad de la Frontera
|
|
|