該樓層疑似違規已被系統折疊?隱藏此樓查看此樓

I=imread('C:\Users\Administrator\Desktop\123.jpg');

figure(1)

subplot(3,2,1),imshow(I), title('原始圖像');

I1=rgb2gray(I);

subplot(3,2,2),imshow(I1),title('灰度圖像');

I2=edge(I1,'roberts',0.09,'both');

subplot(3,2,3),imshow(I2),title('邊緣檢測后圖像');

se=[1;1;1];

I3=imerode(I2,se);

subplot(3,2,4),imshow(I3),title('腐蝕后邊緣圖像');

se=strel('rectangle',[25,25]);

I4=imclose(I3,se);

subplot(3,2,5),imshow(I4),title('填充后圖像');

I5=bwareaopen(I4,2000);

subplot(3,2,6),imshow(I5),title('形態濾波后圖像');

[y,x,z]=size(I5);

I6=double(I5);

Y1=zeros(y,1);

for i=1:y

for j=1:x

if(I6(i,j,1)==1)

Y1(i,1)= Y1(i,1)+1;

end

end

end

[temp MaxY]=max(Y1);

figure(2)

subplot(3,2,1),plot(0:y-1,Y1),title('行方向像素點灰度值累計和'),xlabel('行值'),ylabel('像素');

PY1=MaxY;

while((Y1(PY1,1)>=50)&&(PY1>1))

PY1=PY1-1;

End

PY2=MaxY;

while((Y1(PY2,1)>=50)&&(PY2

PY2=PY2+1;

End

IY=I(PY1:PY2,:,:);

X1=zeros(1,x);

for j=1:x

for i=PY1:PY2

if(I6(i,j,1)==1)

X1(1,j)= X1(1,j)+1;

End

End

End

subplot(3,2,2),plot(0:x-1,X1),title('列方向像素點灰度值累計和'),xlabel('列值'),ylabel('像數');

PX1=1;

while((X1(1,PX1)<3)&&(PX1

PX1=PX1+1;

End

PX2=x;

while((X1(1,PX2)<3)&&(PX2>PX1))

PX2=PX2-1;

End

PX1=PX1-1;

PX2=PX2+1;

dw=I(PY1:PY2,PX1:PX2,:);

subplot(3,2,3),imshow(dw),title('定位剪切后的彩色車牌圖像')

if isrgb(dw)

I1 = rgb2gray(dw);

Else

I1=dw;

End

g_max=double(max(max(I1)));

g_min=double(min(min(I1)));

T=round(g_max-(g_max-g_min)/3);

[m,n]=size(I1);

%h=graythresh(I1);

I1=im2bw(I1,T/256);

subplot(3,2,4);

imshow(I1),title('二值化車牌圖像');

I2=bwareaopen(I1,20);

subplot(3,2,5);

imshow(I2),title('形態學濾波后的二值化圖像');