0

我正在尝试使用 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

太感谢了！！

4

1 回答 1

0

虽然不是 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 过滤器的标记（或标签）图像所期望的。

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

于 2016-04-01T13:36:45.573 回答