1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100
|
/*=========================================================================
*
* Copyright Insight Software Consortium
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#include "itkMIRegistrationFunction.h"
#include "itkGaussianImageSource.h"
#include "itkTestingMacros.h"
#include <iostream>
/**
* This test uses two 2D-Gaussians (standard deviation RegionSize/2)
* One is shifted by 5 pixels from the other.
*
* This test exercises the MIRegistrationFunction
*
*/
int itkMIRegistrationFunctionTest( int, char* [] )
{
const unsigned int ImageDimension = 2;
typedef double PixelType;
typedef itk::Vector< PixelType, ImageDimension > DeformationPixelType;
typedef itk::Image< PixelType,ImageDimension > MovingImageType;
typedef itk::Image< PixelType,ImageDimension > FixedImageType;
typedef itk::Image< DeformationPixelType, ImageDimension >
DisplacementFieldType;
// Declare Gaussian image sources
typedef itk::GaussianImageSource< MovingImageType > MovingImageSourceType;
typedef itk::GaussianImageSource< FixedImageType > FixedImageSourceType;
// Create the two images
FixedImageType::SizeValueType fixedImageSize[] = { 100, 100 };
MovingImageType::SizeValueType movingImageSize[] = { 100, 100 };
FixedImageType::SpacingValueType fixedImageSpacing[] = { 1.0f, 1.0f };
MovingImageType::SpacingValueType movingImageSpacing[] = { 1.0f, 1.0f };
FixedImageType::PointValueType fixedImageOrigin[] = { 0.0f, 0.0f };
MovingImageType::PointValueType movingImageOrigin[] = { 0.0f, 0.0f };
MovingImageSourceType::Pointer movingImageSource = MovingImageSourceType::New();
FixedImageSourceType::Pointer fixedImageSource = FixedImageSourceType::New();
fixedImageSource->SetSize( fixedImageSize );
fixedImageSource->SetOrigin( fixedImageOrigin );
fixedImageSource->SetSpacing( fixedImageSpacing );
fixedImageSource->SetNormalized( false );
fixedImageSource->SetScale( 250.0f );
movingImageSource->SetSize( movingImageSize );
movingImageSource->SetOrigin( movingImageOrigin );
movingImageSource->SetSpacing( movingImageSpacing );
movingImageSource->SetNormalized( false );
movingImageSource->SetScale( 250.0f );
movingImageSource->Update();
fixedImageSource->Update();
MovingImageType::Pointer movingImage = movingImageSource->GetOutput();
FixedImageType::Pointer fixedImage = fixedImageSource->GetOutput();
// Set up the metric
typedef itk::MIRegistrationFunction< FixedImageType,
MovingImageType,
DisplacementFieldType >
MetricFunctionType;
MetricFunctionType::Pointer metricFunction = MetricFunctionType::New();
EXERCISE_BASIC_OBJECT_METHODS( metricFunction, MIRegistrationFunction,
PDEDeformableRegistrationFunction );
// Plug the images into the metric function
metricFunction->SetFixedImage( fixedImage );
metricFunction->SetMovingImage( movingImage );
std::cout << "Test finished." << std::endl;
return EXIT_SUCCESS;
}
|