zuoye.zip_分割图像分类_图像形状分类_形状 分割_形状识别_贝叶斯 图像

format long; i=imread('6.png'); %读取图像 figure，imshow(i); a=rgb2gray(i); %将彩色图像转换成灰度图像 %bw=edge(a,'canny'); %边缘检测 %bw=edge(a,'prewitt'); %bw=edge(a,'sobel'); %bw=edge(a,'roberts'); bw=edge(a,'log'); imshow(bw); se = strel('disk',4); %创造一个平坦的圆盘型结构元素，其半径为 2 bw = imclose(bw，se); %关闭图像 ------------------------- bw = imfill(bw,'holes');%填充图像中的孔洞 imshow(bw); L = bwlabel(bw); %在二维图像中标记连接分量s = regionprops(L， 'centroid'); % 测量图像区域的质心 s= regionprops(L,'centroid'); dt= regionprops(L,'area'); %测量图像区域的面积 cv = regionprops(L， 'perimeter'); %测量图像区域的周长 dim = size(s) % 返回图像的质心的行数和列数，维数 BW_filled = imfill(bw,'holes'); boundaries = bwboundaries(BW_filled); % 寻找边缘不包括孔imshow(bw); figure;imshow(i); hold on; for k=1:dim(1) b= boundaries{k}; dim = size(b) for i=1:dim(1) khoangcach{k}(1,i) = sqrt ( ( b(i,2) - s(k).Centroid(1) )^2 + ( b(i,1) - s(k).Centroid(2) )^2 ) end a=max(khoangcach{k}); b=min(khoangcach{k}); c=dt(k).Area; dolech=a-b; vuong = c/(4*b^2) chunhat=c/(4*b*(a^2-b^2)^0.5); tamgiacdeu=(c*4^0.5)/((a+b)^2); elip =c/(a*b*pi); thoi= (c*( a^2 - b^2 )^0.5) / (2*a^2*b) if dolech < 10 text(s(k).Centroid(1)-20, s(k).Centroid(2), '圆形') elseif (vuong < 1.05 ) && (vuong > 0.95 ) text(s(k).Centroid(1)-20,s(k).Centroid(2), '正方形 ') elseif (elip < 1.05 ) && (elip > 0.95 ) text(s(k).Centroid(1)-20, s(k).Centroid(2), '椭圆形 ') elseif (thoi < 1.05 ) && (thoi > 0.95 ) text(s(k).Centroid(1)-20, s(k).Centroid(2), '菱形') elseif ((chunhat <1.05) && (chunhat >0.95)) text(s(k).Centroid(1)-20, s(k).Centroid(2), '长方形 ') elseif (tamgiacdeu < 1.05 ) && (tamgiacdeu > 0.95 ) text(s(k).Centroid(1)-20, s(k).Centroid(2), '三角形 ') end end