如何解决分水岭的过度分割

user5622144 发表于 Dev

用户名

我的输入图片是

输出图像是

预期的输出是这样的

可以看到有些类似椭圆的结构已与矩形合并。此外，我无法分离每个标签来获得椭圆。使用的算法是分水岭

clear; close all;
I = imread('Sub.png');
I = rgb2gray(I);
figure; imshow(I)
I2 = imtophat(I, strel('square', 45));
figure; imshow(I2)
% Alpha=.047;
% h = fspecial('motion', 10, 5);
% w=gausswin(I2,Alpha)  % you'll have to play with N and alpha
% I2 = imfilter(I2,h,'same','symmetric'); % something like these options
level = .047;
BW = im2bw(I2,level);
D = -bwdist(~BW,'chessboard');
D(~BW) = -Inf;

L = watershed(D);
imshow(label2rgb(L,'jet','w'))

最终开放代码：

ImageSource=imread('cameraman.tif');
ImTmp=ImageSource
ImResidue = zeros(size(ImageSource));
ImIndicator= zeros(size(ImageSource));
ImValues= zeros(size(ImageSource));
For size= 1 : N
se = strel('square',N);
ImOp  = imopen(ImageSource,se);
ImDiff=imabsdiff(ImOp,ImTmp)
if ImResidue < ImDiff then
               ImResidue = ImDiff
               ImIndicator = size
               ImValues = ImOp
end
ImTmp=ImOp
end

FiReTiTi

如果您想要准确的内容，则必须使用带有标记的分水岭，但这将变得更加棘手。默认情况下，基本分水岭在路段上，因为它使用每个局部最小值作为标记。

因此，您必须对图像进行一些预处理，以增加要分割的对象之间的间隔，然后使用标记来引导分水岭。

[根据您的编辑编辑]如果您只想要椎骨之间的小结构，那么我建议进行小的侵蚀，以增加它们之间的间隙，然后进行最终开口。当椎骨需要较大的椎骨时，所需的结构将以较小的半径消失。

不要忘记在渐变图像上使用标记。

[编辑2，初步结果]我对您的问题感到好奇，因此尝试了一下。我没有去到椎骨之间的小区域（要分割的区域），而是尝试先分割椎骨（想要在它们之间的区域）。

这是我所做的：

小开口（正方形1阶，正方形更快，更细，因为您可以使患者在图像中很好地定向，否则可以取盘或六角形），以增加椎骨及其附近区域之间的间隙。
区域开口（表面23，但这并不重要），以使不同区域平坦，因此消除峰。
区域闭合（表面23，但这并不重要），以便再次使不同区域变平，因此填充孔。此时查看图像结果。瞧，一切看起来都比较流畅，但是不同的边界/边缘仍然完好无损。
终极开放（UO）。
根据最终开放结果（残留，值和指标）的阈值。参见椎骨近似。我没有值了（我已经删除了代码），但是您可以查看一下UO结果并找到它们。但是，如果您愿意，我仍然有UO结果。
打开（磁盘顺序7）以擦除所有伪像和误报（椎骨很大）
与5相同的操作，以近似要分割的小图案。查看结果
为了使外部标记（在椎骨之间）具有在步骤7上产生的较小腐蚀结果。我用步骤6的结果的大扩张边界（磁盘顺序11）完成了外部标记。
这是我得到的标记。
用计算的标记进行分水岭，这是初步结果