我想把圆形度小于0.5的区域留下，其他区域去除要怎样改程序，小白求教

clear;close all;
%%
%读取源图像
I = imread('999.png');
figure;imshow(I);
%%
%灰度化、取反
h = rgb2gray(I);
figure;imshow(h);%灰度图像
h = imcomplement(h);%取反
figure;imshow(h);
%%
%中值滤波、二值化
h = medfilt2(h,[4,4]);
bw = im2bw(h,graythresh(h));
%%
%消除噪点
se = strel('disk',2);
bw = imclose(bw,se);
figure;imshow(bw);
%%
%填补闭合图形，填充色为白色
bw  = imfill(bw,'holes');
%%
%边界寻找
[B,L] = bwboundaries(bw,'noholes');
% 为每个闭合图形设置颜色显示
figure;imshow(label2rgb(L, @jet, [.5 .5 .5]))
hold on
for k = 1:length(B)
  boundary = B{k};
  plot(boundary(:,2), boundary(:,1), 'w', 'LineWidth', 2)
end
%%
%计算面积
stats = regionprops(L,'Area','Centroid');
threshold = 0.64;
% 循环处理每个边界，length(B)是闭合图形的个数,即检测到的陶粒对象个数
for k = 1:length(B) % 获取边界坐标'
  boundary = B{k}; % 计算周长
  delta_sq = diff(boundary).^2;
  perimeter = sum(sqrt(sum(delta_sq,2))); 
 % 对标记为K的对象获取面积
  area = stats(k).Area; 
 % 圆度计算公式4*PI*A/P^2
  metric = 4*pi*area/perimeter^2;  
% 结果显示
  metric_string = sprintf('%2.2f',metric);
if metric < threshold
    centroid = stats(k).Centroid;
    plot(centroid(1),centroid(2),'ko');
  end
%设置显示字体
  text(boundary(1,2)-35,boundary(1,1)+13,metric_string,'Color','y',...
       'FontSize',14,'FontWeight','bold');
end

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答

本次提问扣除的有问必答次数，已经为您补发到账户，我们后续会持续优化，扩大我们的服务范围，为您带来更好地服务。