数字图像处理分块程序.doc

彩色图像的直方图均衡化处理 %彩色图像的直方图均衡化护理 sourcepic=imread(C:\MATLAB7\work\qq2015.jpg); figure,imshow(sourcepic) graypic=sourcepic(:,:,1); %计算各灰度出现的概率 gp=zeros(1,256); for i=1:256 gp(i)=length(find(graypic==(i-1)))/(m*n); end newgp=zeros(1,256); s1=zeros(1,256);s2=zeros(1,256); tmp=0; for i=1:256 tmp=tmp+gp(i); s1(i)=tmp; s2(i)=round(s1(i)*256); end for i=1:256 newgp(i)=sum(gp(find(s2==i))); end newgraypic=graypic; for i=1:256 newgraypic(find(graypic==(i-1)))=s2(i); end nr=newgraypic; graypic=sourcepic(:,:,2); gp=zeros(1,256); for i=1:256 gp(i)=length(find(graypic==(i-1)))/(m*n); end newgp=zeros(1,256); s1=zeros(1,256);s2=zeros(1,256); tmp=0; for i=1:256 tmp=tmp+gp(i); s1(i)=tmp; s2(i)=round(s1(i)*256); end for i=1:256 newgp(i)=sum(gp(find(s2==i))); end newgraypic=graypic; for i=1:256 newgraypic(find(graypic==(i-1)))=s2(i); end ng=newgraypic; graypic=sourcepic(:,:,3); gp=zeros(1,256); for i=1:256 gp(i)=length(find(graypic==(i-1)))/(m*n); end newgp=zeros(1,256); s1=zeros(1,256);s2=zeros(1,256); tmp=0; for i=1:256 tmp=tmp+gp(i); s1(i)=tmp; s2(i)=round(s1(i)*256); end for i=1:256 newgp(i)=sum(gp(find(s2==i))); end newgraypic=graypic; for i=1:256 newgraypic(find(graypic==(i-1)))=s2(i); end nb=newgraypic; res=cat(3,nr,ng,nb); figure,imshow(res,[]); RGB转化为HSI和HSI转化为RGB function hsi = rgb2hsi(rgb) %抽取图像分量 rgb = im2double(rgb); r = rgb(:,:,1); g = rgb(:,:,2); b = rgb(:,:,3); %执行转换方程 num = 0.5*((r-g) + (r-b)); den = sqrt((r-g).^2 + (r-b).*(g-b)); theta = acos(num./(den+eps));%防止除数为0 H = theta; H(bg) = 2*pi - H(bg); H=H/(2*pi); num=min(min(r,g),b); den=r+g+b; den(den==0)=eps;%防止除数为0 S=1-3.*num./den; H(S==0)=0; I=(r+g+b)/3; %将3个分量联合成为一个HSI图像 hsi=cat(3,H,S,I); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function rgb=hsi2rgb(hsi) H=hsi(:,:,1)*2*pi; S=hsi(:,:,2); I=hsi(:,:,3); R=zeros(size(hsi,1),size(hsi,2)); G=zeros(size(hsi,1),size(hsi,2)); B=zeros(size(hsi,1),size(hsi,2)); idx=find((0=H)(H2*pi/3)); B(idx)=I(idx).*(1-S(idx))