jnd算法，最小可觉差

function Sxy=jnd（A）
%最小可觉差Jnd计算，及人眼刺激图计算。
%A=imread（‘D:\我的文档\MATLAB\新建文件夹\评价方法\CCD.bmp‘）;
%计算背景灰度图I
[m n]=size（A）;
%A=rgb2gray（A）;
A=double（A）;
for i=1:m
    for j=1:n
       flag=2        if flag
            for a=1:5
                for b=1:5
                    C（ab）=A（i-3+aj-3+b）;
                end
            end
            C=double（C）;
            B=[11111;12221;12021;12221;11111]; 
            II=（1/32）*C.*B;
            I（ij）=sum（sum（II））;
        else
            I（ij）=A（ij）;
        end
    end
end
I=double（I）;
%计算某像素点处，其背景灰度值所对应的jnd值
for i=1:m
    for j=1:n
        if 0<=I（ij）&I（ij）<=60
            E（ij）=log2（1+（A（ij）-I（ij））/I（ij））;
            E0（ij）=log2（1+（-0.009*I（ij）+0.575））;
            Sxy（