person3_102411.m

有些对人脸图象进行学习的初步程序,还在进一步优化过程中,与大家分享.

clear;
%%function L = downsample(H,n)
s1='w_34_';
% s3='.jpg';
% s4='a.jpg';
s3='.bmp';
s4='b.bmp';

H=cell(51,1); 
L=cell(51,1);
n=4;    %downsample by n
%  XL = downsample(H,n);
% for k=1:1
for k=1:51
    s2 = num2str(k);
    str1=strcat(s1,s2,s3);      %str1=w_34_k.bmp          %input string
    str2=strcat(s1,s2,s4);      %str2=w_34_ka.bmp         %output string
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
         h=imread(str1);%错误
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    figure,imshow(h);
    H{k}=h;      
%     for i=1:(256/n)
%        for j=1:(256/n)
             for i=1:(100/n)
                for j=1:(100/n)
          l(i,j)= h(n*i,n*j);       %将采样
       end
    end
    figure,imshow(l);
    L{k}=l;
    [w,h]=size(l);
    imwrite(l,str2);
    %figure,imshow(l);
end
%%patch processing
W=cell(1024,51);
R=cell(1024,51);
for k=1:51
    %for n=1:256
        n=1;
        for a=1:32 %low-resolution images
            for b=1:32
                 for i=1:2 
                     for j=1:2
               % w(i,j)=L{k,1}((n-1)*4+i,(n-1)*4+j);
                           W{n,k}(i,j)=L{k,1}((a-1)*2+i,(b-1)*2+j);
                           %n=n+1;
                     end
                 end
                 n=n+1;
            end
        end
        %W{n,k}=w(i,j);
        n=1;
         for a=1:32  %high-resolution images
            for b=1:32
                 for i=1:8
                       for j=1:8
                          %r(i,j)=H{k,1}((n-1)*16+i,(n-1)*16+j);
                          R{n,k}(i,j)=H{k,1}((a-1)*8+i,(b-1)*8+j);
                          %n=n+1;
                       end
                 end
                 n=n+1;
            end
         end
        %R{n,k}=r(i,j);
    %end
end
%input low-resolution image
T=cell(1024,1);
%for n=1:256
n=1;
for a=1:32
        for b=1:32
            for i=1:2
                for j=1:2
                %test(i,j)=L{1,1}((n-1)*4+i,(n-1)*4+j);
                    T{n}(i,j)=L{1,1}((a-1)*2+i,(b-1)*2+j);
                    %n=n+1;
                end
            end
            n=n+1;
        end
end
%T{n}=test(i,j);
%%match and restruction
n=1;
for a=1:32
    for b=1:32
       % for n=1:256
            for k=2:51
            %q=abs(T{n}-W{n,k})
            %a{k}=sum(sum(abs(T{n}-W{n,k})))
            q=abs(double(W{n,k})-double(T{n}))
            p(k)=sum(sum(q))
            %             if (k>1)&(a(k)>a(k-1))
%                   b=k-1;
%                 else
%                   b=k;
%             end
            end  
           [C,I]=min(p(2:51));
           for i=1:2
               for j=1:2
                    testL((a-1)*2+i,(b-1)*2+j)=W{n,I+1}(i,j);
               end
           end
           for i=1:8
               for j=1:8
                    testH((a-1)*8+i,(b-1)*8+j)=R{n,I+1}(i,j);
               end
           end
           n=n+1;
    end
end
imwrite(testL,'1-1024-3-123L.jpg');
imwrite(testH,'1-1024-3-123H.jpg');
figure,imshow(testL);
figure,imshow(testH);