matlab基础上不变矩处理和直方图均衡化.docVIP

1、求图像不变矩效果图： 程序：I0=imread(1.jpg);if isrgb(I0)==1I1=rgb2gray(I0)endI2=imrotate(I1,-4,bilinear); %逆时针旋转4度I3=fliplr(I1); %镜像subplot(141),imshow(I1);subplot(142),imshow(I2);subplot(143),imshow(I3);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%计算I1七个不变矩%（1）计算I1的而质心I2I3%（2）计算各阶中心矩I2I3%（3）计算不变中心矩的值I2I3 以下的需要替换%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%I1=double(I1);[M,N]=size(I1);[x,y]=meshgrid(1:M,1:N);x=x(:);y=y(:);I1=I1(:);m00=sum(I1);m10=sum(x.* I1);m01=sum(y.* I1);xcenter=m10/m00;ycenter=m01/m00;% 计算归一化的各阶中心矩，具有平移不变性和伸缩不变性cm00=m00;cm02=sum((y-ycenter).^2.*I1)/(m00^2);cm11=sum((x-xcenter).*(y-ycenter).*I1)/(m00^2);cm12=sum((x-xcenter).*(y-ycenter).^2.*I1)/(m00^2.5);cm21=sum((x-xcenter).^2.*(y-ycenter).*I1)/(m00^2.5);cm03=sum((y-ycenter).^3.*I1)/(m00^2.5);cm30=sum((x-xcenter).^3.*I1)/(m00^2.5);cm20=sum((x-xcenter).^2.*I1)/(m00^2);%计算七个不变矩im(1)=abs(log(cm20+cm02))im(2)=abs(log((cm20-cm02)^2+4*cm11^2))im(3)=abs(log((cm30-3*cm12)^2+(3*cm21-cm03)^2))im(4)=abs(log((cm30+cm12)^2+(cm21+cm03)^2))im(5)=abs(log((cm30-3*cm12)*(cm30+cm12)*((cm30+cm12)^2-3*(cm21+cm03)^2)+(3*cm21-cm03)*(cm21+cm03)*(3*(cm30+cm12)^2-(cm21+cm03)^2)))im(6)=abs(log((cm20-cm02)*((cm30+cm12)^2-(cm21+cm03)^2)+4*cm11*(cm30-cm12)*(cm21+cm03)))im(7)=abs(log((3*cm21-cm02)*(cm21+cm30)*((cm30+cm12)^2-3*(cm21+cm03)^2)-(cm30-3*cm12)*(cm21+cm03)*(3*(cm30+cm12)^2-(cm21+cm03)^2))) 2、直方图均衡化效果图 程序： clear all%一，图像的预处理，读入彩色图像将其灰度化PS=imread(1.jpg);???????????????? %读入JPG彩色图像文件imshow(PS)????????????????????????????????? %显示出来 title(输入的彩色JPG图像)imwrite(rgb2gray(PS),PicSampleGray.bmp); %将彩色图片灰度化并保存PS=rgb2gray(PS);??????????????????????????? %灰度化后的数据存入数组 %二，绘制直方图[m,n]=size(PS);????????????????????????????? %测量图像尺寸参数GP=zeros(1,256);??????????????????????????? %预创建存放灰度出现概率的向量for k=0:255???? GP(k+1)=length(find(PS==k))/(m*n);????? %计算每级灰度出现的概率，将其存入GP中相应位置endfigure,bar(0:255,GP,g)??????????????????? %绘制直方图title(原图像直方图)