过冷水：浅析数模大赛实用的Matlab六大图像处理技巧原程序

clc;clearA=imread('1.jpg');imshow(A)imwrite(A,'2.jpg','quality',0.2)B=imread('2.jpg');subplot(1,2,1);imshow(A);title('原始图像');subplot(1,2,2);imshow(B);title('压缩图像');

clc;clearA=imread('3.jpg');B=rgb2gray(A);C=imbinarize(B);subplot(1,3,1);imshow(A);title('原始图像');subplot(1,3,2);imshow(B);title('灰度图像');subplot(1,3,3);imshow(C);title('二值化');

%图像边缘检测clc;clearA=imread('4.jpg');B=rgb2gray(A);C=imbinarize(B);D=edge(C);subplot(1,2,1);imshow(C);title('原始图像');subplot(1,2,2);imshow(D);title('边缘检测');

clc;clearA=imread('4.jpg');B=rgb2gray(A);C=imbinarize(B);D=edge(C)[col,row]=find(D == 1);S=polyarea(col,row)

%检查小圆的个数clc;clearimage = imread('5.jpg');[height,width,channel] =size(image);%确定输入突变的长宽rc =sqrt(width^2 height^2)/2; %搜索圆的最大半径n=30;  %划分区间个数dr = rc/n;%确定半径区间%检测输入图片上的圆gray =rgb2gray(image);  %灰度化bw = imbinarize(gray);   %二值化bw_inverse =imcomplement(bw); %二值图像反转，即0变1,1变0[centers,radius] =imfindcircles(bw_inverse,[5,20]); %检测图像上半径在5~20之间的圆，也可以是其它，看实际情况需要particle_num =length(radius(:)); %统计圆的个数%显示图像figure;imshow(bw_inverse);hold on;viscircles(centers,radius,'EdgeColor','b');%计算每个单位圆内的原子个数[row,col] = size(centers);

%图像亮度处理warning offf=imread('6.jpg');subplot(2,2,1);imshow(f);g2=imadjust(f,[0;1],[1;0]);subplot(2,2,2);imshow(g2);g3=imadjust(f,[0;1],[0.5;0]);subplot(2,2,3);imshow(g3);g4=imadjust(f,[0;0.3],[1;0.4]);subplot(2,2,4);imshow(g4);

%%线性空间滤波f=imread('7.jpg');h1=fspecial('laplacian',0);g1=f-imfilter(f,h1);h2=[1 1 1; 1 -8 1; 1 1 1];g2=f-imfilter(f,h2);subplot(1,3,1),imshow(f);title('原图');subplot(1,3,2),imshow(g1);title('线性加噪');subplot(1,3,3),imshow(g2);title('除噪');

clearclcf=imread('9.jpg'); %读原图像f=rgb2gray(f);f1=imnoise(f,'salt & pepper',0.09); %加椒盐噪声g=medfilt2(f1); %进行中值滤波subplot(1,3,1),imshow(f),title('原图像')subplot(1,3,2),imshow(f1),title('被椒盐噪声污染的图像')subplot(1,3,3),imshow(g),title('中值滤波图像')

clearclcf=imread('8.jpg');imshow(f)f1=imnoise(f(:,:,1),'salt & pepper',0.09); %加椒盐噪声f2=imnoise(f(:,:,2),'salt & pepper',0.09); %加椒盐噪声f3=imnoise(f(:,:,3),'salt & pepper',0.09); %加椒盐噪声f_noise(:,:,1)=f1;f_noise(:,:,2)=f2;f_noise(:,:,3)=f3;imshow(f_noise)g1=medfilt2(f1);%中值滤波g2=medfilt2(f2);%中值滤波g3=medfilt2(f3);%中值滤波g_1(:,:,1)=g1;g_1(:,:,2)=g2;g_1(:,:,3)=g3;subplot(1,3,1),imshow(f),title('原图像');subplot(1,3,2),imshow(f_noise),title('添加椒盐噪声的图像');subplot(1,3,3),imshow(g_1),title('处理以后的图像');

clearA=imread('10.jpg');subplot(2,2,1);imshow(A);title('原图')[M,N]=size(A);%原始图像A的尺寸一MxNu1=4;u2=4;x1(1)=0.2;x2(1)=0.7;sumA=sum(sum(A)); k=mod(sumA,256)*1.0/255;x1(1)=(x1(1) 0.7804)/2;x2(1)=(x2(1) 0.7809)/2; y1(1)=(1/3.1415926)*asin(sqrt(x1(1)));y2(1)=(1/3.1415926)*asin(sqrt(x2(1))); for i=1:1:M*N-1 %产生2个Logistic混沌序列      x1(i 1)=u1*x1(i)*(1-x1(i));      x2(i 1)=u2*x2(i)*(1-x2(i));endfor i=1:1:M*N  %改造2个Logistic,得到两个y序列     y1(i)=(1/3.1415926)*asin(sqrt(x1(i)));        y2(i)=(1/3.1415926)*asin(sqrt(x2(i)));end n=1; for i=1:1:M   %由yl和y2序列对原图像进行值替代加密        for j=1:1:N               if mod(n,1)==0                     k(n)=mod(floor(y1(n)*10^15),256);              else           k(n)=mod(floor(y2(n)*10^15),256);              end        A1(i,j)=bitxor(A(i,j),k(n));  %得到加密像素             n=n 1;        endendsubplot(2,2,3);imshow(A1);title('混沌加密图像')  %输出得到的加密图像n=1; for i=1:1:M        for j=1:1:N                 if mod(n,1)==0                        k(n)=mod(floor(y1(n)*10^15),256);                else           k(n)=mod(floor(y2(n)*10^15),256);                end        A2(i,j)=bitxor(A1(i,j),k(n));  %得到加密像素             n=n 1;        endendsubplot(2,2,4);imshow(A2);title('解密图像')  %输出得到的解密图像

%彩色图像加密img=imread('11.jpg');%img代表原图像矩阵subplot(1,3,1);imshow(img);title('原图') img(1,1,1)img=double(img);sizel=size(img)img_1=img(:);%img_1代表拉直后的图像矩阵password=[8,12,14,16,3,7,9];vx=password[x1,x2]=size(password);%密码的长和宽k=1;while x2<1000    p2=password*3    password=[password p2]    [x1,x2]=size(password)    k=k 1end[m,n]=size(password)img_new=zeros(size(img));%新图像的矩阵c1=zeros(size(img));j=1;for x=1:sizel(1)    for y=1:sizel(2)        for z=1:sizel(3)            img_new(x,y,z)=img(x,y,z) password(j);            img_new(x,y,z)=mod(img_new(x,y,z),256);            j=j 1;            if j==n                j=1;            end        end    endendsubplot(1,3,2);imshow(uint8(img_new));title('加密图') %%图像解密img=double(img_new);j=1;for x=1:sizel(1)    for y=1:sizel(2)        for z=1:sizel(3)            img_new(x,y,z)=img(x,y,z)-password(j);            img_new(x,y,z)=mod(img_new(x,y,z),256);            j=j 1;            if j==n                j=1;            end        end    endendsubplot(1,3,3);imshow(uint8(img_new));title('解密图')