아래 소스코드는 MODIS 영상의 4, 6번 밴드를 이용해서 NDSI라는 눈지수를 계산하는 소스 코드 입니다
궁금한점이 b4 = envi_get_data(fid=fid1, dims=dims1, pos=pos1[1]) 부분이 4번 밴드를 불러들이는 것 같은데 pos=pos[1]이 pos=pos[3] 이 되어야 하는게 아닌가 생각이되는데 이부분이 확실치 않아서 문의 드립니다.
pro snow07_ndsi_gui, event
compile_opt IDL2
; Open corrected image
data1 = envi_pickfile(filter=’*.img’, title=’Select input corrected image File’)
if (data1 eq ”) then return
print, ‘The selected filename is: ‘, data1
envi_open_file, data1, r_fid=fid1
envi_file_query, fid1, dims=dims1, nb=nb1, nl=nl1, ns=ns1
pos1 = lindgen(nb1)
; Band 4
b4 = envi_get_data(fid=fid1, dims=dims1, pos=pos1[1])
; Band 6
b6 = envi_get_data(fid=fid1, dims=dims1, pos=pos1[5])
; Output file name
out_name = envi_pickfile(filter=’*.img’, /output, title=’Output File Name’)
if (out_name eq ”) then return
print, ‘Output file name is: ‘, out_name
; Compute NDSI
ndsi = (b4-b6)/(b4+b6)
; Get Info
inherit1 = envi_set_inheritance(fid1, dims1, pos1, /full)
bnames = [‘NDSI’]
; Save data
ENVI_WRITE_ENVI_FILE, ndsi, $
out_name = out_name, $
bnames = bnames, $
inherit = inherit1
end