itk Casting and Intensity Mapping

1 Linear Mappings

主要的类：

itk::CastImageFilter;

itk::RescaleIntensityImageFilter;

itk::ShiftScaleImageFilter;

itk::NormalizeImageFilter.

以上四个类的变换规则分别是：

    (1) itk::CastImageFilter:

只对输入图像的每一个像素的类型进行变换，而不对灰度值作其他的算术运算。应用itk::CastImageFilter就相当于调用C风格的类型转换：

outputPixel = static_cast<OutputPixelType>(inputPixel);

(2) itk::RescaleIntensityImageFilter:

把输入图像灰度值的最小最大值分别映射到用户指定的灰度值最小最大值，所采用的映射公式为：

 

输出图像灰度值最小最大值的设置方法：SetOutputMinimum(); SetOutputMaximum().

(3) itk::ShiftScaleImageFilter:

也是对输入图像的灰度值进行线性变换，但是这种线性变换是以用户指定的偏移量和缩放比例应用所得最终结果的，采用的映射公式为：

 

偏移量和缩放比例的设置方法：SetShift(); SetScale().

(4) itk::NormalizeImageFilter:

通过类里面的一些参数设置，使得输出图像的灰度级统计信息为平均数是0，方差是1。采用的映射公式为：

 

2 Non Linear Mappings

主要的类：itk::SigmoidImageFilter

 

类itk::SigmoidImageFilter对输入图像的灰度值进行变换的规则如上式所示，用户调用类itk::SigmoidImageFilter里的方法SetOutputMaximum(), SetOutputMinimum()设置输出图像的灰度值最大值最小值，也就是式中的Max和Min。I表示输入图像的灰度值，I’表示输出图像的灰度值，α表示输入图像的灰度范围，β表示输入图像灰度范围的中心位置(这两个值有点类似窗宽和窗位)。它们的值的设置由方法：SetAlpha(), SetBeta()来完成。一般认为图像的灰度范围为[-3α,3α]之间，α和β的值可正可负。

也就是：对于输入图像的灰度值，如果其像素灰度值低于β-3α或者高于β+3α时，则这些灰度值将被映射到用户设置的Minimum和Maximum的值，如果输入图像的灰度值是在[β-3α, β+3α]之间时，则就会调用以上公式计算出输出图像的灰度值。