|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
先看处理结果! U* Y- @$ g; n- y% Y3 Z9 s
% p0 i L6 ?, K$ f7 ?
( Q) m8 D+ x. ~6 Q8 |8 y/ z2 \" g0 z6 o6 _+ X' l% l2 u
简单说一下实现思路:% ~8 } Q2 q \2 }( w
. p3 f0 P2 Y/ l3 o9 `' ]) I* |
读取图片,转灰度,计算灰度直方图,估算阈值(这里的阈值计算很重要,经过阈值算法,选取一个最恰当的阈值),之后二值化。显示图像即可。后面阈值的计算补更。
7 F' |# R Q, h: _, w4 ^3 t! G; Q9 {: j. [+ C
- %%车牌
- clear all
- clc
- PS=imread('chepai.png');
- subplot(1,2,1);
- imshow(PS)
- title('原图')
- p=rgb2gray(PS);
- subplot(1,2,2)
- imshow(p)
- title('原灰度图')
- [m,n]=size(p);
- GP=zeros(1,256);
- for k=0:255
- GP(k+1)=length(find(p==k))/(m*n); %计算每级灰度出现的概率,将其存入GP,这里的find函数寻找
- %与0-255相等的个数,然后通过length函数返回一个数,再除以
- %总的像素点个数,就是第0个灰度数据出现的次数,放入第一个
- %里面
- end
- figure
- subplot(1,2,1);
- bar(0:255,GP,'g')
- xlabel('灰度值')
- ylabel(' 出现概率')
- title('原灰度直方图')
- max_index=[];
- for i=3:length(GP)-2
- if((GP(i)>=GP(i+1))&(GP(i)>=GP(i-1)))&((GP(i+1)>=GP(i+2))&(GP(i-1)>=GP(i-2)))
- max_index(end+1)=i-1;%数组的扩展,有n个元素,那么第n+1个用来存放数据i-1(因为是从0开始的,所以-1)
- end
- end
- possible=GP(max_index);
- [max_value,index]=max(possible);%前者存值,后者放位置
- TT=max_index(index)-2;%取出数据
- [m,n]=size(p);
- R=zeros(m,n);
- for i=1:m
- for j=1:n
- if p(i,j)<TT
- R(i,j)=256;
- else R(i,j)=0;
- end
- end
- end
- subplot(1,2,2);
- imshow(R);
- title('二值图');
c, U5 R" {5 q& w" _
' z1 g8 l' x9 [ r% r+ Q# V) W4 @* K7 _% f7 B5 v- g
+ M; D0 K. t+ X
* Z, f; Q+ S1 B# f0 t2 J, d" Z/ {' r
' l; P# R* L2 X+ L4 l- d7 e; b1 C) D3 ?
5 L7 c7 J) o7 g# Z7 w6 ~+ r8 H5 g* a
5 J8 x& q& r( D0 R b" |% [ m' {9 ~) B/ J# }9 ~
- y( ^4 j; B8 g3 f3 @- F0 L
|
|