第六章 图像分割
实例:灰度图像的阈值分割
clc;
clear;
%读入图像
A = imread('shape.jpg');
figure(1),subplot(2,2,1),imshow(A),title('原图像');
%获取灰度直方图,根据直方图人工选取阈值
h = imhist(A);
%h = imhist(A) / numel(A);
h1 = h(1:10:256);
horz = 1:10:256;
figure(2),bar(horz, h1),axis([0 255 0 15000]),title('直方图');
%人工设定阈值进行图像分割
thr1 = 100;
[m, n] = size(A);
for i = 1:m
for j = 1:n
if A(i,j) < thr1
A1(i,j) = 0;
else
A1(i,j) = 1;
end
end
end
figure(1),subplot(2,2,2),imshow(A1),title('人工设定阈值');
%使用迭代法选取阈值进行图像分割
%计算阈值初始值
thr2 = 0.5 * (double(min(A(:))) + double(max(A(:))));
%设定循环终止条件
done = false;
%循环计算阈值和分割图像
while ~done
%获取二值分割结果
A2 = A >= th