ITK系列27_ 边缘保留平滑滤波对脑部MHA文件进行三维滤波

女爷i 2023-01-22 06:54 38阅读 0赞

### **实例27** **边缘保留平滑滤波对脑部MHA****文件进行三维滤波** ###

#include "itkImage.h"
    #include "itkCastImageFilter.h"
    #include "itkCurvatureFlowImageFilter.h"
    #include "itkImageFileReader.h"
    #include "itkImageFileWriter.h"
    
    int main(int argc, char* argv[])
    {
        //定义浮点像素类型 InternalPixelType
        typedef   float           InternalPixelType;
        const     unsigned int    Dimension = 3;
        //定义浮点、维数2的图像类型 InternalImageType
        typedef itk::Image< InternalPixelType, Dimension >  InternalImageType;
        //定义无符号字符 输出像素类型 OutputPixelType
        typedef unsigned char                            OutputPixelType;
        //定义无符号字符、维数2 输出图像类型 OutputImageType
        typedef itk::Image< OutputPixelType, Dimension > OutputImageType;
        //定义滤波输入图像类型InternalImageType   滤波输出图像类型OutputImageType
        typedef itk::CastImageFilter< InternalImageType, OutputImageType >
            CastingFilterType;
        //实例化滤波器对象caster
        CastingFilterType::Pointer caster = CastingFilterType::New();
        //定义图像读取类型ReaderType
        typedef  itk::ImageFileReader< InternalImageType > ReaderType;
        //定义图像写类型WriterType
        typedef  itk::ImageFileWriter<  OutputImageType  > WriterType;
    
        ReaderType::Pointer reader = ReaderType::New();
        WriterType::Pointer writer = WriterType::New();
    
        reader->SetFileName("BrainProtonDensity3Slices.mha");
        writer->SetFileName("lvbo_CastingFilter.mha");
    
        //实例化滤波对象smoothing
        typedef itk::CurvatureFlowImageFilter< InternalImageType, InternalImageType >
            CurvatureFlowImageFilterType;
        CurvatureFlowImageFilterType::Pointer smoothing =
            CurvatureFlowImageFilterType::New();
    
        现在连接管道
        smoothing->SetInput(reader->GetOutput());
        //输出滤波后的效果图
        caster->SetInput(smoothing->GetOutput());
        writer->SetInput(caster->GetOutput());
    
        /*IsolatedConnectedImageFilter 期望用户指定一个门限和两个种子。在这个例子中，我们
        从命令行得到它们*/
        smoothing->SetNumberOfIterations(12);
        smoothing->SetTimeStep(0.125);
    
        try
        {
            writer->Update();
        }
        catch (itk::ExceptionObject & excep)
        {
            std::cerr << "Exception caught !" << std::endl;
            std::cerr << excep << std::endl;
        }
    
        return EXIT_SUCCESS;
    }

输入三维图像（BrainProtonDensity3Slices.mha）：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 1][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 2][]

**切片1 切片2 切片3**

输出三维滤波图像（lvbo\_CastingFilter.mha）：

![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 3][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 4][] ![watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 5][]

**切片1 切片2 切片3**

**ITK系列目录：**

**[1 ITK图像数据表达之图像][1 ITK]**

**[2 ITK图像处理之图像滤波][2 ITK]**

**[3 ITK图像处理之图像分割][3 ITK]**

**注：例程配套素材见系列目录**

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70]: https://img-blog.csdnimg.cn/20201203105054261.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 1]: https://img-blog.csdnimg.cn/2020090817513321.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 2]: /images/20221020/2bff7e9b235342bc8c71219ad65ccb71.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 3]: https://img-blog.csdnimg.cn/2020120310490365.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 4]: https://img-blog.csdnimg.cn/20201203104925202.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz,size_16,color_FFFFFF,t_70
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzMyODA5MDkz_size_16_color_FFFFFF_t_70 5]: /images/20221020/15376fef920f44d3a9f3e1805d111365.png
[1 ITK]: https://blog.csdn.net/qq_32809093/article/details/108307568
[2 ITK]: https://blog.csdn.net/qq_32809093/article/details/108882952
[3 ITK]: https://blog.csdn.net/qq_32809093/article/details/108327734