帮我写一段通过MATLAB数学形态学算法识别出一张只有黑白两色且黑色为底色，白色为纯白色小圆形粒子的图像中纯白色小圆形粒子的个数但是分纯白色小圆形粒子存在部分重叠，需要代码能够识别两个或多个部分重叠的白色圆形粒子并将它们单独区分开来加入计算结果的代码

以下是识别圆形粒子数量并区分部分重叠的MATLAB代码：

% 读取图像
img = imread('circles.png');

% 转换为二值图像
bw = im2bw(img, 0.5);

% 对二值图像进行开操作，去除噪声
se = strel('disk', 5);
bw = imopen(bw, se);

% 对二值图像进行闭操作，填充内部空洞
se = strel('disk', 15);
bw = imclose(bw, se);

% 对二值图像进行形态学重建，得到圆形粒子的掩模
bw = imfill(bw, 'holes');
se = strel('disk', 50);
bw = imdilate(bw, se);
bw = imreconstruct(imcomplement(bw), imcomplement(bw) & imcomplement(bw));

% 通过圆形Hough变换检测圆形粒子
circles = imfindcircles(bw, [5 30]);

% 绘制圆形粒子检测结果
imshow(img);
viscircles(circles(:, 1:2), circles(:, 3), 'EdgeColor', 'g');

% 统计圆形粒子数量
num_circles = size(circles, 1);

% 区分部分重叠的圆形粒子
for i = 1:num_circles
    % 提取当前圆形粒子的掩模
    [x, y] = meshgrid(1:size(bw, 2), 1:size(bw, 1));
    mask = hypot(x - circles(i, 1), y - circles(i, 2)) <= circles(i, 3);
    
    % 检测当前圆形粒子与其他圆形粒子的重叠部分
    overlap_masks = {};
    for j = 1:num_circles
        if i == j
            continue;
        end
        [x, y] = meshgrid(1:size(bw, 2), 1:size(bw, 1));
        other_mask = hypot(x - circles(j, 1), y - circles(j, 2)) <= circles(j, 3);
        overlap_mask = mask & other_mask;
        if any(overlap_mask(:))
            overlap_masks{end+1} = overlap_mask;
        end
    end
    
    % 区分重叠的圆形粒子
    if ~isempty(overlap_masks)
        for j = 1:length(overlap_masks)
            % 提取当前重叠部分的掩模
            overlap_mask = overlap_masks{j};
            
            % 计算当前重叠部分的圆形粒子的圆心和半径
            stats = regionprops(overlap_mask, 'Centroid', 'MajorAxisLength', 'MinorAxisLength');
            x = stats.Centroid(1);
            y = stats.Centroid(2);
            r = mean([stats.MajorAxisLength stats.MinorAxisLength]) / 2;
            
            % 在圆形粒子检测结果中添加当前重叠部分的圆形粒子
            circles(end+1, :) = [x y r];
            
            % 在原图像中绘制当前重叠部分的圆形粒子
            hold on;
            viscircles([x y], r, 'EdgeColor', 'r');
        end
    end
end

% 输出圆形粒子数量
num_circles = size(circles, 1);
disp(['Number of circles: ' num2str(num_circles)]);

注意，这段代码假设输入图像中只有黑色和白色两种颜色，并且黑色为底色，白色为纯白色。如果输入图像中存在其他颜色或灰度值，请先进行颜色或灰度值的处理。