tuxing.m

利用matlab对图形进行辨识

clear;clc;
raw = imread('作业2.bmp');
raw = padarray(raw, [1 1], 1, 'both');  
gray_bmp = rgb2gray(raw);
%注意matlab是对图像的白色部分进行处理的
labeled_bmp = ~bwlabel(gray_bmp);   
sep1 = ~bwselect(labeled_bmp, 100, 100);
sep2 = ~bwselect(labeled_bmp, 550, 100);
sep3 = ~bwselect(labeled_bmp, 150, 250);
sep4 = ~bwselect(labeled_bmp, 550, 250);
sep5 = ~bwselect(labeled_bmp, 100, 500);
sep6 = ~bwselect(labeled_bmp, 400, 500);
sep7 = ~bwselect(labeled_bmp, 730, 490);
sep8 = ~bwselect(labeled_bmp, 100, 750);
sep9 = ~bwselect(labeled_bmp, 500, 760);
sep10 = ~bwselect(labeled_bmp, 777, 777);

black = 0;
white = 1;

for i = 1:10
    switch i
        case 1
            sep_bmp = sep1;
        case 2
            sep_bmp = sep2;
        case 3
            sep_bmp = sep3;
        case 4
            sep_bmp = sep4;
        case 5
            sep_bmp = sep5;
        case 6
            sep_bmp = sep6;
        case 7
            sep_bmp = sep7;
        case 8
            sep_bmp = sep8;
        case 9
            sep_bmp = sep9;
        case 10
            sep_bmp = sep10;
    end
    euler(i) = bweuler(sep_bmp);    %计算欧拉数
end
hole = find(euler > 0);    
for i = 1:9
    switch i
        case 1
            sep_bmp = sep1;
        case 2
            sep_bmp = sep2;
        case 3
            sep_bmp = sep3;
        case 4
            sep_bmp = sep4;
        case 5
            sep_bmp = sep5;
        case 6
            sep_bmp = sep6;
        case 7
            sep_bmp = sep7;
        case 8
            sep_bmp = sep8;
        case 9
            sep_bmp = sep9;
    end
    [x, y] = find(sep_bmp == black);    
    top(i) = min(x);   %找到图形的上下左右四个点
    left(i) = min(y);
    bottom(i) = max(x);
    right(i) = max(y);
    area(i) = bwarea(~sep_bmp);     %计算出每个图形的面积
    area_matrix(i) = area(i) / ((bottom(i) - top(i)) * (right(i)-left(i)))*100;    
end

bow = find(area_matrix < 40);
triangle = find(area_matrix>40 & area_matrix<60);
round = find(area_matrix>60 & area_matrix<95);
matrix = find(area_matrix<105 & area_matrix>95);

%按面积大小区别大小矩形
[x, y1] = max([area(matrix(1)) area(matrix(2))]);   
[x, y2] = min([area(matrix(1)) area(matrix(2))]);
big_matrix = matrix(y1);
small_matrix = matrix(y2);

surface(1) = (right(triangle(1)) - left(triangle(1))) / (bottom(triangle(1)) - top(triangle(1)));
surface(2) = (right(triangle(2)) - left(triangle(2))) / (bottom(triangle(2)) - top(triangle(2)));
[x, y1] = max(surface); 
[x, y2] = min(surface);
hor_triangle = triangle(y1);
upright_triangle = triangle(y2);


sep2 = bwlabel(sep2);  
sep9 = bwlabel(sep9);
[x, y] = size(bow);
for i = 1:y
    
    switch bow(i)
        case 2
            sep_bmp = sep2;
            stat1 = regionprops(sep_bmp,'Orientation');
        case 9
            sep_bmp = sep9;
            stat2 = regionprops(sep_bmp,'Orientation');
    end
end
%大于零的是上弦，小于零的是下弦
[x, y1] = max([stat1.Orientation stat2.Orientation]);
[x, y2] = min([stat1.Orientation stat2.Orientation]);
up_bow = bow(y1);
down_bow = bow(y2);

for i = 1:3
    box(i) = (right(round(i)) - left(round(i))) / (bottom(round(i)) - top(round(i)));
end
[x, y] = max(box);
ellipse = round(y); %区别出椭圆来了
%看看图形最上面的像素和其最小矩形的上边重合的像素个数多少，多的是八边形
count_black1 = 0;    %重合的像素个数
count_black2 = 0;
round(1) = round(2);    %从round中去除椭圆，还剩下圆和八边形
round(2) = round(3);
for i = 1:844
    if sep5(top(5), i) == black;
        count_black1 = count_black1 + 1;
    elseif sep8(top(8), i) == black;
        count_black2 = count_black2 + 1;
    end
end
[x, y] = max([count_black1 count_black2]);
eight_border = round(y);
[x, y] = min([count_black1 count_black2]);
round_round = round(y);

num_bmp = [big_matrix small_matrix up_bow down_bow hole upright_triangle hor_triangle eight_border ellipse round_round];
for i = 1:10
    switch num_bmp(i)
        case 1
            bmp = sep1;
        case 2
            bmp = sep2;
        case 3
            bmp = sep3;
        case 4
            bmp = sep4;
        case 5
            bmp = sep5;
        case 6
            bmp = sep6;
        case 7
            bmp = sep7;
        case 8
            bmp = sep8;
        case 9
            bmp = sep9;
        case 10
            bmp = sep10;
    end
    subplot(2, 5, i);imshow(bmp);
end