图像模糊处理：
源程序：
clc;
clear all
I=imread('56.jpg');
A=[1 1 1 1 1;1 1 1 1 1;1 1 1 1 1;1 1 1 1 1;1 1 1 1 1];
I1=im2double(I);
[row,col]=size(I1);
[m,n]=size(A);
L=floor(m/2);
W=floor(n/2);
image=zeros(row+2*L,col+2*W);
for i=1:1:row
for j=1:1:col
image(i+L,j+W)=I1(i,j);
end
end
output=I1;
for i=1+L:1:row-L
for j=1+W:1:col-W
mask=image(i-L:i+L,j-W:j+W);
numlmatrix=mask.*A;
output(i-L,j-W)=sum(sum(numlmatrix))/(m*n); end
end
figure, subplot(121);imshow(I1);title('原图');
subplot(122);imshow(output);title('均值滤波后图像');
图像锐化处理：
源程序：
clc;
clear all
I=imread('56.jpg');
A=[-1 -1 -1;-1 8 -1;-1 -1 -1];
I1=im2double(I);
[row,col]=size(I1);
[m,n]=size(A);
L=floor(m/2);
W=floor(n/2);
image=zeros(row+2*L,col+2*W);
for i=1:1:row
for j=1:1:col
image(i+L,j+W)=I1(i,j);
end
end
output=I1;
for i=1+L:1:row-L
for j=1+W:1:col-W
mask=image(i-L:i+L,j-W:j+W);
numlmatrix=mask.*A;
output(i-L,j-W)=sum(sum(numlmatrix))/(m*n);
end
end
figure, subplot(131);imshow(I1);title('原图');
subplot(132);imshow(output);title('锐化滤波后图像');
ima=I1(:,:,2);
ima=double(ima);
bw1 = edge(ima,'sobel');
subplot(133);imshow(bw1);title('sobel锐化');
直方图均衡化：
源程序：
clear all
clc
I=imread('56.jpg');
PS=rgb2gray(I);
figure,subplot(121);imshow(PS);
title('输入的灰度图');
[m,n]=size(PS);
GP=zeros(1,256);
for k=0:255
GP(k+1)=length(find(PS==k))/(m*n);
end
subplot(122);plot([0:255],GP,'r')
title('原图直方图')
xlabel('灰度值')
ylabel('出现概率')
S1=zeros(1,256);
for i=1:256
for j=1:i
S1(i)=GP(j)+S1(i);
end
end
S2=round((S1*256)+0.5);
for i=1:256
GPeq(i)=sum(GP(find(S2==i)));
end
figure,subplot(122);plot([0:255],GPeq,'g') ;
title('均衡化后直方图');
xlabel('灰度值');
ylabel('出现概率');
PA=PS;
for i=0:255
PA(find(PS==i))=S2(i+1);
end
subplot(121);imshow(PA);
title('均衡化后图像');
边缘检测：
源程序：
clear all;
A = imread ('56.jpg');
I1=im2double(A);
I=rgb2gray(I1);
BW_sobel = edge(I,'sobel');
BW_prewitt = edge(I,'prewitt');
BW_roberts = edge(I,'roberts');
BW_laplace = edge(I,'log');
BW_canny = edge(I,'canny');
figure(1);
subplot(2,3,1),imshow(I),xlabel('?-í?');
subplot(2,3,2),imshow(BW_sobel),xlabel('sobel检测');
subplot(2,3,3),imshow(BW_prewitt),xlabel('prewitt检测');
subplot(2,3,4),imshow(BW_roberts),xlabel('robert检测');
subplot(2,3,5),imshow(BW_laplace),xlabel('laplace检测');
subplot(2,3,6),imshow(BW_canny),xlabel('canny检测');
I2 = imnoise(I,'salt & pepper',0.01);
BW_sobel1 = edge(I2,'sobel');
BW_prewitt1 = edge(I2,'prewitt');
BW_roberts1 = edge(I2,'roberts');
BW_laplace1 = edge(I2,'log');
BW_canny1 = edge(I2,'canny');
figure(2);
subplot(2,3,1),imshow(I2),xlabel('椒盐sobel检测');
subplot(2,3,2),imshow(BW_sobel1),xlabel('椒盐sobel检测');
subplot(2,3,3),imshow(BW_prewitt1),xlabel('椒盐prewitt检测');
subplot(2,3,4),imshow(BW_roberts1),xlabel('椒盐roberts检测');
subplot(2,3,5),imshow(BW_laplace1),xlabel('椒盐laplace检测');
subplot(2,3,6),imshow(BW_canny1),xlabel('椒盐canny检测');
图像分割_1：
源程序：
clear all;
clc;
I=imread('56.jpg');
I1=rgb2gray(I);
I2= edge(I1,'roberts');
se=[1;1;1];
I3=imerode(I2,se);
figure,subplot(131),imshow(I3),title('原图');
se=strel('rectangle',[25,25]);
I4=imclose(I3,se);
subplot(132),imshow(I4),title('腐蚀');
I5=bwareaopen(I4,2000);
subplot(133),imshow(I5),title('去除杂点');
[x,y]=size(I5);
I6=double(I5);
X1=zeros(x,1);
for i=1:x
for j=1:y
if(I6(i,j)==1)
X1(i,1)= X1(i,1)+1;
end
end
end
[tempx MaxX]=max(X1);
PX1=MaxX;
while ((X1(PX1,1)>=10)&&(PX1>1))
PX1=PX1-1;
end
PX2=MaxX;
while ((X1(PX2,1)>=10)&&(PX2<y))
PX2=PX2+1;
end
PX1=PX1-10;
PX2=PX2+10;
IX=I(PX1:PX2,:,:);
I7=rgb2gray(IX);
g_1=double(max(max(I7)));
g_2=double(min(min(I7)));
T=round(g_1-(g_1-g_2)/3);
I8=(double(I7)>=T);
imtool(I8);
[m,n]=size(I8);
Y1=zeros(n,1);
for j=1:n
for i=1:m
if(I8(i,j)==1)
Y1(j,1)= Y1(j,1)+1;
end
end
end
[tempy MaxY]=max(Y1);
PY1=MaxY;
while ((Y1(PY1,1)>1)&&(PY1<n-1))
PY1=PY1+1;
end
PY2=MaxY;
while ((Y1(PY2,1)>1)&&(PY2>1))
PY2=PY2-1;
end
PY1=PY1-10;
PY2=PY2+10;
dw=I(PX1:PX2,PY2:PY1,:);
figure,subplot(1,2,1),imshow(IX),title('行方向数字区域'); subplot(1,2,2),imshow(dw),title('数字区域截取');
imwrite(dw,'dw.jpg');
