clear;
clc;
%%%%%%%%%1、預(yù)處理%%%%%%%%%%%
I=imread（‘川A99999.jpg‘）;
figure（1）;
imshow（I）;
title（‘原始圖像‘）

%?剪切后的圖像
JianQie=imcrop（I[5?7?130?29?]）;
figure（2）;
imshow（JianQie）title（‘剪切后的車牌圖像‘）?????
?????
%灰度圖像
I1=rgb2gray（JianQie）;
figure（3）
imshow（I1）;
title（‘灰度圖像‘）?
???
%同態(tài)濾波
J=I1;
J=double（J）;
f=fft2（J）;
g=fftshift（f）;
[MN]=size（f）;
d0=10;
r1=0.5;
rh=2;
c=4;
n1=floor（M/2）;
n2=floor（N/2）;
for?i=1:M
????for?j=1:N
????????d=sqrt（（i-n1）^2+（j-n2）^2）;
????????h=（rh-r1）*（1-exp（-c*（d.^2/d0.^2）））+r1;
????????g（ij）=h*g（ij）;
????end
end
g=ifftshift（g）;
g=uint8（real（ifft2（g）））;
figure（4）;
imshow（g）;
title（‘同態(tài)濾波‘）;
??????
%Otsu方法二值化圖像
?T=graythresh（g）;
?ChePaiErZhi=im2bw（gT）;
?figure（5）;
?imshow（ChePaiErZhi）;
?title（‘Otsu方法二值化圖像‘）;


%刪除小對(duì)象
ChePaiLvBo=bwareaopen（ChePaiErZhi10）;%10011002---100Car---12
figure（6）;
imshow（ChePaiLvBo）;
title（‘刪除小對(duì)象之后的圖像‘）

[yxz]=size（