ITK上的形态学分水岭标记过滤器

Morphological Watershed From Markers filter on ITK

本文关键字:过滤器 分水岭 形态学 ITK      更新时间:2023-10-16

我正在尝试使用ITK中的库创建一个用于图像分割的管道。但是,当我应用itkMorphologicalWatershedFromMarkersFilter时,结果是一个空白图像(只有1的二进制图像)。

有人知道如何正确应用此筛选器吗?

我的输入图像应该是图像的梯度,标记图像应该是在同一图像上应用分水岭滤波器的结果。

输入图像

标记图像

这是过滤器的声明和应用:

typedef itk::MorphologicalWatershedFromMarkersImageFilter < OutputImageType, OutputImageType >
      MorphologicalWatershedFromMarkersImageFilterType;
MorphologicalWatershedFromMarkersImageFilterType::Pointer CwatershedFilter
      = MorphologicalWatershedFromMarkersImageFilterType::New();

CwatershedFilter->SetInput1(reader1->GetOutput());
CwatershedFilter->SetMarkerImage(reader2->GetOutput());
CwatershedFilter->SetMarkWatershedLine(true);
try{
    CwatershedFilter->Update();
}
catch (itk::ExceptionObject & error)
{
    std::cerr << "Error: " << error << std::endl;
    getchar();
    return EXIT_FAILURE;
}

此外,这是该过滤器文档的链接,来自itk.org/

http://www.itk.org/Doxygen48/html/classitk_1_1MorphologicalWatershedFromMarkersImageFilter.html#a20e3b8de42219606ba759e822be0aaa2

非常感谢!!

虽然不是C++ITK,但有一个SimpleITK Notebooks演示了它的用法:

http://insightsoftwareconsortium.github.io/SimpleITK-Notebooks/32_Watersheds_Segmentation.html

标记图像只是二进制的,而不是标签图像。我的意思是,你的图像只有0和1(或255)。在链接的示例中,请注意以下内容:

min_img = sitk.RegionalMinima(feature_img, backgroundValue=0, foregroundValue=1.0, fullyConnected=False, flatIsMinima=True)
marker_img = sitk.ConnectedComponent(min_img, fullyConnected=False)

"min_img"是一个二进制图像,但随后使用"ConnectedComponent"图像过滤器对图像进行处理,该过滤器为每个"岛"提供一个唯一的编号。这是WatershedFromMarker过滤器的标记(或标签)图像的预期值。

我还将注意到,您的输入图像有一些边界线,您可能不希望将其作为输入。