ITK-均值滤波

作者：翟天保Steven
版权声明：著作权归作者所有，商业转载请联系作者获得授权，非商业转载请注明出处

均值滤波原理

均值滤波是一种常用的图像平滑技术，用于减少图像中的噪声。其基本原理是通过计算图像中每个像素及其邻域像素的平均值来替换该像素的值，从而达到平滑图像的效果。

概括为以下几个步骤：

选择滤波窗口：定义一个窗口（通常是矩形或方形），该窗口覆盖目标像素及其周围的像素。窗口的大小由用户指定，常见的大小有 3×3、5×5等。
移动窗口：将窗口依次移动到图像中的每一个像素位置，逐个处理每个像素。
计算平均值：对于窗口覆盖的每一个位置，计算窗口内所有像素的平均值。
替换像素值：用计算得到的平均值替换当前窗口中心像素的值。

均值滤波可以有效地平滑图像，减少高频噪声，使图像看起来更加平滑。但它也有一些缺点，例如会模糊图像中的边缘，因为它同样平滑了边缘像素和非边缘像素。对于边缘检测等需要保留图像细节的任务，均值滤波可能不是最优选择。

环境准备

参见：Windows下用CMake编译ITK及配置测试_itk配置-CSDN博客

功能解析

1.引入必要的头文件：

#include <itkImage.h>
#include <itkImageFileReader.h>
#include <itkImageFileWriter.h>
#include <itkMeanImageFilter.h>
#include <itkRescaleIntensityImageFilter.h>
#include <itkJPEGImageIOFactory.h>
#include <itkCastImageFilter.h>

2.初始化图像类型和读写器：

typedef itk::Image<unsigned char, 2> CharImageType;
typedef itk::Image<float, 2> FloatImageType;
typedef itk::ImageFileReader<CharImageType> ReaderType;
typedef itk::ImageFileWriter<CharImageType> WriterType;
itk::JPEGImageIOFactory::RegisterOneFactory();
ReaderType::Pointer reader = ReaderType::New();
WriterType::Pointer writer = WriterType::New();

3.设置文件名：

// 设置
reader->SetFileName("test.jpg");            // 要读取的文件名
writer->SetFileName("mean_output.jpg");    // 写入的文件名

4.类型转换：

// 类型转换
typedef itk::CastImageFilter<CharImageType, FloatImageType> CastFilterType;
CastFilterType::Pointer castfilter = CastFilterType::New();
castfilter->SetInput(reader->GetOutput());

5.创建均值滤波器，并配置参数：

// 创建均值滤波器
typedef itk::MeanImageFilter<FloatImageType, FloatImageType> FilterType;
FilterType::Pointer filter = FilterType::New();
filter->SetInput(castfilter->GetOutput());
// 设置参数
FilterType::InputSizeType radius;
radius.Fill(3); // 设置均值滤波器的半径
filter->SetRadius(radius);

6.类型转换：

// 将浮点图像转换为unsigned char类型
typedef itk::RescaleIntensityImageFilter<FloatImageType, CharImageType> RescaleFilterType;
RescaleFilterType::Pointer castfilter2 = RescaleFilterType::New();
castfilter2->SetInput(filter->GetOutput());

7.连接过滤器输出到写入器并执行写入操作：

writer->SetInput(castfilter2->GetOutput());
try
{writer->Update();
}
catch (itk::ExceptionObject &error)
{std::cerr << "Error: " << error << std::endl;return EXIT_FAILURE;
}

完整代码

#include <itkImage.h>
#include <itkImageFileReader.h>
#include <itkImageFileWriter.h>
#include <itkMeanImageFilter.h>
#include <itkRescaleIntensityImageFilter.h>
#include <itkJPEGImageIOFactory.h>
#include <itkCastImageFilter.h>int main()
{// 初始化typedef itk::Image<unsigned char, 2> CharImageType;typedef itk::Image<float, 2> FloatImageType;typedef itk::ImageFileReader<CharImageType> ReaderType;typedef itk::ImageFileWriter<CharImageType> WriterType;itk::JPEGImageIOFactory::RegisterOneFactory();ReaderType::Pointer reader = ReaderType::New();WriterType::Pointer writer = WriterType::New();// 设置reader->SetFileName("test.jpg");            // 要读取的文件名writer->SetFileName("mean_output.jpg");    // 写入的文件名// 类型转换typedef itk::CastImageFilter<CharImageType, FloatImageType> CastFilterType;CastFilterType::Pointer castfilter = CastFilterType::New();castfilter->SetInput(reader->GetOutput());// 创建均值滤波器typedef itk::MeanImageFilter<FloatImageType, FloatImageType> FilterType;FilterType::Pointer filter = FilterType::New();filter->SetInput(castfilter->GetOutput());// 设置参数FilterType::InputSizeType radius;radius.Fill(3); // 设置均值滤波器的半径filter->SetRadius(radius);// 将浮点图像转换为unsigned char类型typedef itk::RescaleIntensityImageFilter<FloatImageType, CharImageType> RescaleFilterType;RescaleFilterType::Pointer castfilter2 = RescaleFilterType::New();castfilter2->SetInput(filter->GetOutput());// 输出writer->SetInput(castfilter2->GetOutput());try{writer->Update();}catch (itk::ExceptionObject &error){std::cerr << "Error: " << error << std::endl;return EXIT_FAILURE;}std::cout << "Mean filtering completed successfully." << std::endl;return EXIT_SUCCESS;
}