OIF计算代码matlab

multBand = imread（‘C:\1111.tif‘）;
% 组合数
n = 3;
% 波段数
band = size（multBand）;
combin = combntns（1:6n）;
sizeCombin = size（combin）;
OIF = zeros（sizeCombin（1）sizeCombin（2）+1）;
OIF（:1:end-1） = combin;
for row = 1:sizeCombin（1）
    bands = multBand（::combin（row:））;
    STDS = 0;
    for b = 1:n
        STDS = STDS + std2（bands（::b））;
    end
    CORRS = 0;
    for b1=