《itk实用demo》-删除canny图像中比较小的部分

最新推荐文章于 2024-01-22 20:40:37 发布

winston600

最新推荐文章于 2024-01-22 20:40:37 发布

阅读量516

点赞数

CC 4.0 BY-SA版权

分类专栏： ITK医学图像处理

本文链接：https://blog.youkuaiyun.com/rabbitbride/article/details/82422278

ITK医学图像处理专栏收录该内容

94 篇文章

订阅专栏

本文介绍了一种使用ITK库进行Canny边缘检测后的图像处理方法，通过调整像素强度、二值化转换、区域标记和去除小面积区域等步骤，实现了对Canny检测结果的有效后处理。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

删除canny图像中比较小的部分

    typedef itk::RescaleIntensityImageFilter<ImageType2D, ImageType2D> RescaleIntensityImageFilterType;
    RescaleIntensityImageFilterType::Pointer rescaleFliter = RescaleIntensityImageFilterType::New();
    rescaleFliter->SetInput(image);
    rescaleFliter->SetOutputMinimum(0);
    rescaleFliter->SetOutputMaximum(5);
    //------------------------------------------------------------
    typedef itk::BinaryImageToLabelMapFilter<ImageType2D> BinaryImageToLabelMapFilterType;
    BinaryImageToLabelMapFilterType::Pointer binaryImageToLabelMapFilter = BinaryImageToLabelMapFilterType::New();
    binaryImageToLabelMapFilter->SetInput(rescaleFliter->GetOutput());
    binaryImageToLabelMapFilter->SetFullyConnected(true);
    binaryImageToLabelMapFilter->SetInputForegroundValue(5);
    binaryImageToLabelMapFilter->Update();

    std::cout << "There are " << binaryImageToLabelMapFilter->GetOutput()->GetNumberOfLabelObjects() << " objects." << std::endl;

    std::vector<unsigned long> labelsToRemove;

    for(unsigned int i = 0; i < binaryImageToLabelMapFilter->GetOutput()->GetNumberOfLabelObjects(); i++)
    {
        // Get the ith region
        BinaryImageToLabelMapFilterType::OutputImageType::LabelObjectType* labelObject = binaryImageToLabelMapFilter->GetOutput()->GetNthLabelObject(i);
        if (ifDebug)
        {
            std::cout << "Region " << i << " has " << labelObject->Size() << " pixels." << std::endl;
        }        

        // Mark every other label to be removed
        //if(i%2 == 0)
        if(labelObject->Size() < 20)
        {
            labelsToRemove.push_back(labelObject->GetLabel());
        }
    }
    std::cout << "Removing " << labelsToRemove.size() << " objects." << std::endl;
    // Remove all regions that were marked for removal.
    for(unsigned int i = 0; i < labelsToRemove.size(); ++i)
    {
        binaryImageToLabelMapFilter->GetOutput()->RemoveLabel(labelsToRemove[i]);
    }

    std::cout << "There are " << binaryImageToLabelMapFilter->GetOutput()->GetNumberOfLabelObjects() 
        << " objects remaining." << std::endl;

    typedef itk::LabelMapToLabelImageFilter<BinaryImageToLabelMapFilterType::OutputImageType, ImageType2D > L2IType;
    L2IType::Pointer l2i = L2IType::New();
    l2i->SetInput( binaryImageToLabelMapFilter->GetOutput());
    l2i->Update();

    typedef itk::BinaryThresholdImageFilter<ImageType2D, ImageType2D>  ThresholdFilterType;
    ThresholdFilterType::Pointer thresholdfilter = ThresholdFilterType::New();
    thresholdfilter->SetInput( l2i->GetOutput() );
    thresholdfilter->SetLowerThreshold(1);
    thresholdfilter->SetUpperThreshold(255);
    thresholdfilter->SetInsideValue(255);
    thresholdfilter->SetOutsideValue(0);
    thresholdfilter->Update();
    //----------------------------------
    ImageType2D::Pointer image_temp =ImageType2D::New();
    image_temp->Graft(thresholdfilter->GetOutput());

    return image_temp;