itk中的基本數學運算

數字圖像處理的本質就是數學運算，今天我們看下itk中的有哪些基本的數學運算。

//1、加，減，乘（注意，沒有除哦）
#include "itkAddImageFilter.h"//加
#include "itkSubtractImageFilter.h"//減
#include "itkMultiplyImageFilter.h"//乘

typedef itk::AddImageFilter <ImageType, ImageType >	mathFilterType;
mathFilterType::Pointer mathFilter	= mathFilterType::New ();
mathFilter->SetInput1(input1);
mathFilter->SetInput2(input2);
mathFilter->Update();

typedef itk::SubtractImageFilter <ImageType, ImageType >	mathFilterType;
mathFilterType::Pointer mathFilter	= mathFilterType::New ();
mathFilter->SetInput1(input1);
mathFilter->SetInput2(input2);
mathFilter->Update();

typedef itk::MultiplyImageFilter <ImageType, ImageType >	mathFilterType;
mathFilterType::Pointer mathFilter	= mathFilterType::New ();
mathFilter->SetInput1(input1);
mathFilter->SetInput2(input2);
mathFilter->Update();

//2、與，或，非
#include "itkAndImageFilter.h"//與
#include "itkOrImageFilter.h"//或
#include "itkBinaryNotImageFilter.h"//非

typedef itk::AndImageFilter <ImageType, ImageType >	mathFilterType;
mathFilterType::Pointer mathFilter	= mathFilterType::New ();
mathFilter->SetInput1(input1);
mathFilter->SetInput2(input2);
mathFilter->Update();

typedef itk::OrImageFilter <ImageType, ImageType >	mathFilterType;
mathFilterType::Pointer mathFilter	= mathFilterType::New ();
mathFilter->SetInput1(input1);
mathFilter->SetInput2(input2);
mathFilter->Update();

typedef itk::BinaryNotImageFilter <ImageType>       mathFilterType; 
mathFilterType::Pointer mathFilter = mathFilterType::New();
mathFilter ->SetInput(input1);
mathFilter ->Update();

//3、膨脹，腐蝕，開，閉
#include "itkBinaryBallStructuringElement.h"//基本球形
#include "itkBinaryMorphologicalOpeningImageFilter.h"//開運算
#include "itkBinaryMorphologicalClosingImageFilter.h"//閉運算
#include "itkGrayscaleErodeImageFilter.h"//灰階膨脹
#include "itkGrayscaleDilateImageFilter.h"//灰階膨脹
#include "itkBinaryDilateImageFilter.h"//二值膨脹
#include "itkBinaryErodeImageFilter.h"//二值膨脹

typedef itk::BinaryBallStructuringElement<ImageType::PixelType, ImageType::ImageDimension>
StructuringElementType;
StructuringElementType structuringElement;
structuringElement.SetRadius(rad);
structuringElement.CreateStructuringElement();

typedef itk::BinaryMorphologicalOpeningImageFilter <ImageType, ImageType, StructuringElementType> BinaryMorphologicalOpeningImageFilterType;
BinaryMorphologicalOpeningImageFilterType::Pointer mathOperationFilter
			= BinaryMorphologicalOpeningImageFilterType::New();
mathOperationFilter->SetInput(*input);
mathOperationFilter->SetKernel(structuringElement);
mathOperationFilter->Update();

typedef itk::BinaryMorphologicalClosingImageFilter <ImageType, ImageType, StructuringElementType>	BinaryMorphologicalClosingImageFilterType;
BinaryMorphologicalClosingImageFilterType::Pointer mathOperationFilter
			= BinaryMorphologicalClosingImageFilterType::New();
mathOperationFilter->SetInput(*input);
mathOperationFilter->SetKernel(structuringElement);
mathOperationFilter->Update();

typedef itk::GrayscaleErodeImageFilter <ImageType, ImageType, StructuringElementType>
			GrayscaleErodeImageFilterType;
GrayscaleErodeImageFilterType::Pointer mathOperationFilter = GrayscaleErodeImageFilterType::New();
mathOperationFilter->SetInput(*input);
mathOperationFilter->SetKernel(structuringElement);
mathOperationFilter->Update();

typedef itk::GrayscaleDilateImageFilter <ImageType, ImageType, StructuringElementType>
			GrayscaleDilateImageFilterType;
GrayscaleDilateImageFilterType::Pointer mathOperationFilter = GrayscaleDilateImageFilterType::New();
mathOperationFilter->SetInput(*input);
mathOperationFilter->SetKernel(structuringElement);
mathOperationFilter->Update();

typedef itk::BinaryErodeImageFilter  <ImageType, ImageType, StructuringElementType>
			BinaryErodeImageFilter ;
BinaryErodeImageFilter ::Pointer mathOperationFilter = BinaryErodeImageFilter ::New();
mathOperationFilter->SetInput(*input);
mathOperationFilter->SetKernel(structuringElement);
mathOperationFilter->Update();

typedef itk::BinaryDilateImageFilter  <ImageType, ImageType, StructuringElementType>
			BinaryDilateImageFilter ;
BinaryDilateImageFilter ::Pointer mathOperationFilter = BinaryDilateImageFilter ::New();
mathOperationFilter->SetInput(*input);
mathOperationFilter->SetKernel(structuringElement);
mathOperationFilter->Update();

typedef itk::BinaryErodeImageFilter  <ImageType, ImageType, StructuringElementType>
			BinaryErodeImageFilter ;
BinaryErodeImageFilter ::Pointer mathOperationFilter = BinaryErodeImageFilter ::New();
mathOperationFilter->SetInput(*input);
mathOperationFilter->SetKernel(structuringElement);
mathOperationFilter->SetBackgroundValue(0);
mathOperationFilter->SetForegroundValue(1);
mathOperationFilter->Update();

typedef itk::BinaryDilateImageFilter  <ImageType, ImageType, StructuringElementType>
			BinaryDilateImageFilter ;
BinaryDilateImageFilter ::Pointer mathOperationFilter = BinaryDilateImageFilter ::New();
mathOperationFilter->SetInput(*input);
mathOperationFilter->SetKernel(structuringElement);
mathOperationFilter->SetBackgroundValue(0);
mathOperationFilter->SetForegroundValue(1);
mathOperationFilter->Update();
           

Σ( ° △ °|||)︴，代碼寫的多了，文采果然越來越爛，除了注釋，居然不知道寫些什麼。這樣下去不是辦法，我會盡力從最簡單的寫起。