%%%%%%%%%%%簡(jiǎn)單的車牌識(shí)別系統(tǒng)%%%%%%%%%%%
%%%%%%%???????鐘培勛????????%%%%%%%%%%?
%%%%%%%?????QQ：930109175????%%%%%%%%%%

clc;
clear?all;
close?all;
%%?step1?讀入圖片?灰度化處理并顯示原圖、灰度圖

%選擇圖片路徑
[filenamepathname]=...
uigetfile（{‘*.jpg‘;‘*.bmp‘;‘*gif‘}‘選擇圖片‘）;
%合成路徑+文件名
str0=[pathname?filename];
%讀圖
I=imread（str0）;
%imshow（I）;
%I=imread（‘./original?pic/car7.jpg‘）;???%imread函數(shù)讀取圖像文件
[yxz]=size（I）;
if?y>400
????rate=400/y;
????I=imresize（Irate）;
end
%imshow（I）;
%I=imread（‘./original?pic/car7.jpg‘）;???%imread函數(shù)讀取圖像文件
%將彩色圖像轉(zhuǎn)換為黑白并顯示
I1?=?rgb2gray（I）;????%rgb2gray轉(zhuǎn)換成灰度圖
figure（1）imshow（I）title（‘原始彩色圖像‘）;????%figure命令同時(shí)顯示兩幅圖像
whos?I;
figure（2）imshow（I1）title（‘原始黑白圖像‘）;

%Step2?圖像預(yù)處理???對(duì)原始黑白圖像進(jìn)行開操作得到圖像背景??
%I1為灰度圖
[mn]=size（I1）;????????????????????????????????????????%測(cè)量圖像尺寸參數(shù)
GreyHist=zeros（1256）;???????????????????????????????????????%預(yù)創(chuàng)建存放灰度出現(xiàn)概率的向量
for?k=0:255
????GreyHist（k+1）=length（find（I1==k））/（m*n）;?????????????????%計(jì)算每級(jí)灰度出現(xiàn)的概率，將其存入GreyHist中相應(yīng)位置
end
figure（3）
subplot（222）;
bar（0:255GreyHist‘g‘）????????????????????????????????????%繪制直方圖???
title（‘拉伸前灰度直方圖‘）
xlabel（‘灰度值‘）
ylabel（‘出現(xiàn)概率‘）
subplot（221）imshow（I1）title（‘拉伸前黑白圖像‘）;
%灰度拉伸
I1=double（I1）;
ma=double（max（max（I1）））;
mi=double（min（min（I1）））;
I1=（255/（ma-mi））*I1-（255*mi）/（ma-mi）;
I1=uint8（I1）;
%figure（4）?
subplot（223）;
imshow（I1）;
title（‘灰度拉伸后黑白圖像‘）;
for?k=0:255
????GreyHist（k+1）=length（find（I1==k））/（m*n）;?????????????????
end
subplot（224）;
bar（0:255GreyHist‘b‘）??????????????????????????????????????
title（‘拉伸后的灰度直方圖‘）
xlabel（‘灰度值‘）
ylabel（‘出現(xiàn)概率‘）


%突出目標(biāo)對(duì)象
SE=strel（‘disk‘16）;%半徑為r=15的圓的模板
I2=imopen（I1SE）;%開運(yùn)算?????用模板SE對(duì)灰度圖I1進(jìn)行腐蝕，再對(duì)腐蝕后的結(jié)果進(jìn)行膨脹，使外邊緣圓滑
figure（4）imshow（I2）;title（‘背景圖像‘）;%輸出背景圖像
%用原始圖像與背景圖像作減法，增強(qiáng)圖像
I3=imsubtract（I1I2）;%兩幅圖相減
figure（5）imshow（I3）;title（‘增強(qiáng)黑白圖像‘）;%輸出黑白圖像

%Step3?取得最佳閾值，將圖像二值化
fmax1=double（max（max（I3）））;%I3的最大值并輸出雙精度型
fmin1=double（min（min（I3）））;%I3的最小值并輸出雙精度型
T=（fmax1-（fmax1-fmin1）/3）/255;%獲得最佳閾值
bw22=im2bw（I3T）;%轉(zhuǎn)換圖像為二進(jìn)制圖像
bw2=double（bw22）;
figure（6）imshow（bw2）;title（‘圖像二值化‘）;%得到二值圖像

%%%%%%%%%%%%%%%%??車牌定位模塊??%%%%%%%%%%%%%%%%%%%%%%%%%%
%%????數(shù)學(xué)形態(tài)學(xué)處理進(jìn)行車牌粗定位
%%????采用彩色像素點(diǎn)統(tǒng)計(jì)，行列掃描的方式實(shí)現(xiàn)車牌精確定位
%%????即改進(jìn)的像素中線掃描法（像素統(tǒng)計(jì)法）

%Step4?車牌粗定位，對(duì)得到二值圖像進(jìn)行邊緣檢測(cè)和開閉操作進(jìn)行數(shù)字形態(tài)學(xué)處理

grd=edge（bw2‘canny‘）%用canny算子識(shí)別強(qiáng)度圖像中的邊界
figure（7）imshow（grd）;title（‘Canny算子圖像邊緣提取‘）;%輸出圖像邊緣
bg1=imclose（grdstrel（‘rectangle‘[519]））;%取矩形框的閉運(yùn)算
figure（8）imshow（bg1）;title（‘圖像閉運(yùn)算[519]‘）;%輸出閉運(yùn)算的圖像
bg3=imopen（bg1strel（‘rectangle‘[519]））;%取矩形框的開運(yùn)算
figure（9）imshow（bg3）;title（‘圖像開運(yùn)算[519]‘）;%輸出開運(yùn)算的圖像
bg2=imopen（bg3strel（‘rectangle‘[115]））;%取矩形框的開運(yùn)算
%bg2=bwareaopen（bg2）;%消除細(xì)小對(duì)象
figure（10）imshow（bg2）;title（‘圖像開運(yùn)算[115]‘）;%輸出開運(yùn)算的圖像
bg2=bwareaopen（bg25）;%消除細(xì)小對(duì)象
figure（11）imshow（bg2）;title（‘消除小對(duì)象‘）;

%Step5?像素中線掃描（顏色紋理范圍定義，行列掃描的方式）粗定位和經(jīng)驗(yàn)閾值分割車牌

%%%%%%%%%%%%%%%%??Y方向?%%%%%%%%%%%%%%%%
%進(jìn)一步確定y方向（水平方向）的車牌區(qū)域
[yxz]=size（bg2）;??%y方向?qū)?yīng)行，x方向?qū)?yīng)列，z方向?qū)?yīng)深度，z=1為二值圖像
myI=double（bg2）;??%數(shù)據(jù)類型轉(zhuǎn)換，每個(gè)方向范圍在0~1??0為黑

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-05-05?16:21??LPRS\
?????文件????????1940??2018-05-05?16:32??LPRS\1.車牌灰度圖像.jpg
?????文件????????2456??2018-05-05?16:32??LPRS\2.車牌二值圖像.jpg
?????文件????????2387??2018-05-05?16:32??LPRS\3.消除小對(duì)象二值車牌圖像.jpg
?????文件????????2679??2018-05-05?16:32??LPRS\Plate彩色圖.jpg
?????文件????????1922??2018-05-05?16:32??LPRS\Plate校正后圖像.jpg
?????文件????????1922??2018-05-05?16:32??LPRS\Plate灰度圖.jpg
?????目錄???????????0??2018-05-01?16:45??LPRS\ceshi\
?????文件??????189981??2018-05-01?16:03??LPRS\ceshi\car.jpg
?????文件??????104309??2018-05-01?16:03??LPRS\ceshi\car13.jpg
?????文件??????174779??2018-05-01?16:01??LPRS\ceshi\云AG289N.jpg
?????文件??????106363??2018-05-01?16:03??LPRS\ceshi\吉AGH827.jpg
?????文件???????37993??2017-04-12?20:05??LPRS\ceshi\川F1A666.jpg
?????文件?????3985567??2018-05-01?16:03??LPRS\ceshi\粵A6ZC93.jpg
?????文件?????4070931??2018-05-01?16:03??LPRS\ceshi\粵A6ZC935.jpg
?????文件?????2460690??2018-05-01?16:03??LPRS\ceshi\粵AC609Z31.jpg
?????文件??????322453??2018-05-01?16:01??LPRS\ceshi\粵SX186Q.jpg
?????文件??????115385??2018-05-01?16:03??LPRS\ceshi\遼AQ6493.jpg
?????文件??????113645??2018-05-01?16:03??LPRS\ceshi\遼B5705D.jpg
?????文件????????1940??2018-05-05?16:32??LPRS\lisence.jpg
?????文件???????22803??2018-05-05?16:30??LPRS\main.m
?????文件?????????170??2018-05-01?16:00??LPRS\rando_bianhuan.m
?????目錄???????????0??2018-05-01?16:45??LPRS\word\
?????文件?????????697??2018-05-05?15:28??LPRS\word\1.jpg
?????文件?????????722??2018-05-05?16:32??LPRS\word\2.jpg
?????文件?????????767??2018-05-05?16:32??LPRS\word\3.jpg
?????文件?????????735??2018-05-05?16:32??LPRS\word\4.jpg
?????文件?????????758??2018-05-05?16:32??LPRS\word\5.jpg
?????文件?????????754??2018-05-05?16:32??LPRS\word\6.jpg
?????文件?????????843??2018-05-05?16:32??LPRS\word\7.jpg
?????文件?????????762??2018-05-01?19:12??LPRS\word\8.jpg
............此處省略94個(gè)文件信息