《itk实用demo》-canny算子

最新推荐文章于 2021-12-04 20:07:14 发布

原创最新推荐文章于 2021-12-04 20:07:14 发布 · 438 阅读

0 ·

CC 4.0 BY-SA版权

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

94 篇文章

订阅专栏

本文介绍了一种基于Canny算子的边缘检测算法实现过程，包括中值滤波、顶帽变换等预处理步骤，并详细展示了如何使用ITK库进行图像处理，最终通过Canny算子提取图像边缘。

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

canny算子

void Segment()
{
    ImageType::Pointer sliceimage = image_data;
    itk::TimeProbe clock;
    clock.Start();
    //-------------------------------------------
    //---------canny
    float variance = 1.0;
    float maximumError = 0.10;
    float upperThreshold = 0.5;

    typedef itk::MedianImageFilter<ImageType, ImageType > FilterType;
    FilterType::Pointer medianFilter = FilterType::New();
    FilterType::InputSizeType radius;
    radius.Fill(3);
    medianFilter->SetRadius(radius);
    medianFilter->SetInput( sliceimage);
    //-------------------------------------------
    typedef itk::BinaryBallStructuringElement<PixelType, ImageDimension > SRType;
    SRType kernel;
    kernel.SetRadius( 22 );
    kernel.CreateStructuringElement();

    //typedef itk::WhiteTopHatImageFilter<ImageType,ImageType,SRType > TopHatFilterType;
    typedef itk::BlackTopHatImageFilter<ImageType,ImageType,SRType > TopHatFilterType;
    TopHatFilterType::Pointer topHatFilter = TopHatFilterType::New();
    topHatFilter->SetInput( medianFilter->GetOutput() );
    topHatFilter->SetSafeBorder( false);
    topHatFilter->SetKernel( kernel );

    //------------------------------------------------------------
    typedef double        RealPixelType;  //  Operations
    typedef itk::Image< RealPixelType, ImageDimension> RealImageType;
    typedef itk::CastImageFilter< ImageType, RealImageType >    CastToRealFilterType;    
    typedef itk::CannyEdgeDetectionImageFilter< RealImageType, RealImageType > CannyFilterType;
    typedef itk::RescaleIntensityImageFilter< RealImageType, ImageType > RescaleFilterType;

    CastToRealFilterType::Pointer toReal      = CastToRealFilterType::New();
    CannyFilterType::Pointer      cannyFilter = CannyFilterType::New();
    RescaleFilterType::Pointer    rescale     = RescaleFilterType::New();

    toReal->SetInput( topHatFilter->GetOutput() );
    cannyFilter->SetInput( toReal->GetOutput() );
    rescale->SetInput( cannyFilter->GetOutput() );

    cannyFilter->SetVariance( variance );
    cannyFilter->SetUpperThreshold( upperThreshold );
    cannyFilter->SetMaximumError( maximumError );
    //------------------------------------------------------------
    rescale->Update();
    image_canny = rescale->GetOutput();

    coronary_data = image_canny;

    clock.Stop();
    std::cout << "Algorithm Total time: " << clock.GetTotal() << std::endl;
}