直方图规格化.doc

% 对图像进行规格化处理（匹配）filename,pathname=uigetfile;img=imread(strcat(pathname,filename);m,n,k=size(img);figure(1);subplot(221);imshow(img);title(Original Bmp);pr=zeros(256,1);pg=zeros(256,1);pb=zeros(256,1);% 对R值的规格化for i=1:m for j=1:n temp=img(i,j,1); pr(temp+1)=pr(temp+1)+1; end;end;pr=pr/(m*n);pr_sum=zeros(256,1);for t=1:256 for i=1:t pr_sum(t)=pr(i)+pr_sum(t); end;end;% 对G值的规格化for i=1:m for j=1:n temp=img(i,j,2); pg(temp+1)=pg(temp+1)+1; end;end;pg=pg/(m*n);pg_sum=zeros(256,1);for t=1:256 for i=1:t pg_sum(t)=pg(i)+pg_sum(t); end;end;% 对B值的规格化for i=1:m for j=1:n temp=img(i,j,3); pb(temp+1)=pb(temp+1)+1; end;end;pb=pb/(m*n);pb_sum=zeros(256,1);for t=1:256 for i=1:t pb_sum(t)=pb(i)+pb_sum(t); end;end;% 读取目标图像Objfilename,Objpathname=uigetfile;Objimg=imread(strcat(Objpathname,Objfilename);om,on,ok=size(Objimg);Opr=zeros(256,1);Opg=zeros(256,1);Opb=zeros(256,1);subplot(222);imshow(Objimg);title(Objection Bmp);% 对R值的规格化for i=1:om for j=1:on temp=Objimg(i,j,1); Opr(temp+1)=Opr(temp+1)+1; end;end;Opr=Opr/(om*on);Opr_sum=zeros(256,1);for t=1:256 for i=1:t Opr_sum(t)=Opr(i)+Opr_sum(t); end;end;% 对G值的规格化for i=1:om for j=1:on temp=Objimg(i,j,2); Opg(temp+1)=Opg(temp+1)+1; end;end;Opg=Opg/(om*on);Opg_sum=zeros(256,1);for t=1:256 for i=1:t Opg_sum(t)=Opg(i)+Opg_sum(t); end;end;% 对B值的规格化for i=1:om for j=1:on temp=Objimg(i,j,3); Opb(temp+1)=Opb(temp+1)+1; end;end;Opb=Opb/(om*on);Opb_sum=zeros(256,1);for t=1:256 for i=1:t Opb_sum(t)=Opb(i)+Opb_sum(t); end;end;%tempGray=0;for x=1:m for y=1:n temp=1.0; t=img(x,y,1); for i=t+1:255 if(abs(pr_sum(t+1)-Opr_sum(i)temp) temp=abs(pr_sum(t+1)-Opr_sum(i); tempGray=i; end; end; img(x,y,1)=tempGray; end;end;%tempGray=0;for x=1:m for y=1:n temp=1.0; t=img(x,y,2); for i=t+1:255 if(abs(pg_sum(t+1)-Opg_sum(i)temp) temp=abs(pg_sum(t+1)-Opg_sum(i); tempGray=i; end; end; img(x,y,2)=tempGray; end;end;%tempGray=0;for x=1:m for y=1:n temp=1.0; t=img(x,y,3); for i=t+1:255 if(abs(pb_sum(t+1)-Opb_sum(i)temp) temp=abs(pb_sum(t+1)-O