IDL相关系数计算

pro correlation
    
    ;初始化工作
    COMPILE_OPT idl2
    ENVI, /restore_base_save_files
    envi_batch_init
    firstfile = dialog_pickfile(title = '请输入第一幅多波段影像')
    envi_open_file, firstfile, r_fid=ffid
    if (ffid eq -1) then begin
        envi_batch_exit
        return
    endif
    envi_file_query, ffid, dims=fdims, nb=fnb
    secondfile = dialog_pickfile(title = '请输入第二幅多波段影像')
    envi_open_file, secondfile, r_fid=sfid
    if (sfid eq -1) then begin
        envi_batch_exit
        return
    endif
    envi_file_query, sfid, dims=sdims, nb=snb  ;初始化工作结束
    
    if fnb lt snb then begin
        nb = fnb
    endif else begin
        nb = snb
    endelse
    
    fpos = lindgen(fnb)
    spos = lindgen(snb)          ;定义用于保存波段号的数组
    outdata = make_array(nb)   ;定义用于保存各波段间相关系数的数组