自商图像SQI的MATLAB程序

对光照变化下的图像进行自商图像（SQI）处理的一种方法，增加了单亮区和单暗区，便于对光照进行更好的处理。

clc;
clear;
X=imread（‘D:\人脸识别论文\程序\程序\第四章（自商图像）\自商图像（可以用）\002.bmp‘）;
figure
subplot（221）;
imshow（X）

[mn]=size（X）;
allmean=mean（mean（X））;
sum1=0;sum2=0;
for i=1:m
    for j=1:n
        if X（ij）            sum1=sum1+X（ij）;
            sum2=sum2+1;
        end
    end
end
threshold=sum2/sum1;
A=zeros（mn）;
for i=1:m
    for j=1:n
        if X（ij）>threshold
            A（ij）=1;
        else
            A（ij）=0;
        end
    end
end
subplot（222）;
imshow（A）
%对图像进行二值化

%利用模板确定图像中比较暗的部分
M=damb（mn）;
Afan=zeros（mn）;
for i=1:m
    for j=1:n
        if A（ij）==0 && M（ij）==1
            Afan（ij）=1;
        else
            Afan（ij）=0;
        end
    end
end
%断开连接处理
se1=strel（‘square‘8）;
Afan=imerode（Afanse1）;
Afan=imdilate（Afanse1）; 
subplot（223）;
imshow（Afan）

[Afannum]=bwlabeln（Afan4）;
area1=0.01*m*n
area2=0.07*m*n
B=zeros（mn）;
for i=1:num
    l=length（find（Afan==i））;
    if l>=area1 && l<=area2
        B（find（Afan==i））=1;
    end
end
%确定出需要调整的暗区

X（find（B==1））=50;
subplot（224）;
imshow（X）
imwrite（X‘D:\abcdefj.bmp‘‘bmp‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2014-11-12 10:33  SQI\
     文件        1214  2014-10-05 17:01  SQI\alonedark.m
     文件        1416  2014-10-05 17:20  SQI\alonelight.m
     文件         320  2014-10-05 11:01  SQI\damb.m
     文件         950  2014-10-23 21:45  SQI\iterative.m
     文件        2099  2014-10-19 15:47  SQI\SQI.m