天天看點

itk中的基本數學運算

數字圖像處理的本質就是數學運算,今天我們看下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();
           

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