layer stacking question [message #56568] |
Sun, 28 October 2007 03:02 |
woopik
Messages: 18 Registered: September 2007
|
Junior Member |
|
|
hi
i make layer stacking programaticly with idl
i have an problem , when i link two images
one - i make with idl, second - i make
with envi GUI interface and i verify cursor location
value , data in corresponding pixels are not the same
why ???
i use teh same .tif files for input in the same order
they have the same output size , map info ,
(the .tif are georeferenced )
long integer data type , same resampling
any suggestion ?
; ************************
tb_fid = lonarr(nblay)
pos = lonarr(nblay)
dims = lonarr(5,nblay)
for i=0,nblay-1 do begin
envi_open_data_file, ncd[i],r_fid=fid
tb_fid[i]=fid
envi_file_query, tb_fid[i], ns=t_ns, nl=t_nl, nb=t_nb
dims[0,i] = [-1,0,t_ns-1,0,t_nl-1]
pos[i] = 0
endfor
; ************************
; ************************
out_proj = envi_get_projection(fid=tb_fid[0], pixel_size=out_ps)
out_dt = 3
; ************************
envi_doit, 'envi_layer_stacking_doit', fid=tb_fid, pos=pos, dims=dims,
$
out_dt=out_dt, out_name=out_name, interp=1, out_ps=out_ps, $
out_proj=out_proj, r_fid=r_fid
|
|
|