www.pudn.com > Brovery.rar > Brovery.m

clear; 
x0=imread('d:\data\象素层数据\tm1.bmp'); 
[a,b,c]=size(x0); 
x=double(x0)/255;  
y=imread('d:\data\象素层数据\spot1.bmp'); 
y1=double(y)/255; 
x1=zeros(a,b); 
x2=zeros(a,b); 
x3=zeros(a,b); 
for f=1:a 
    for e=1:b 
        xx(f,e)=x(f,e,1)+x(f,e,2)+x(f,e,3); 
        x1(f,e)=x(f,e,1)*y1(f,e)/xx(f,e); 
        x2(f,e)=x(f,e,2)*y1(f,e)/xx(f,e); 
        x3(f,e)=x(f,e,3)*y1(f,e)/xx(f,e); 
    end 
end 
for i=1:a 
    for j=1:b 
        p(i,j,1)=x1(i,j); 
        p(i,j,2)=x2(i,j); 
        p(i,j,3)=x3(i,j); 
    end 
end 
fu=uint8(round(p*255)); 
figure,imshow(fu); 
 
 
%%%%%%%%%%%%%%%%%%%%%%%% 以下为各主成分归一化后再融合 %%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%% 第一主成分 %%%%%%%%%%%%%%%%%% 
h=x1(1,1);%最大值 
for i=1:a 
    for j=1:b 
        if x1(i,j)>h 
            h=x1(i,j); 
        else 
        end 
    end 
end 
o=x1(1,1); % 最小值 
for i=1:a 
    for j=1:b 
        if x1(i,j)h 
            h=x2(i,j); 
        else 
        end 
    end 
end 
o=x2(1,1); % 最小值 
for i=1:a 
    for j=1:b 
        if x2(i,j)h 
            h=x3(i,j); 
        else 
        end 
    end 
end 
o=x3(1,1); % 最小值 
for i=1:a 
    for j=1:b 
        if x3(i,j)