ITK  4.9.0
Insight Segmentation and Registration Toolkit
WikiExamples/ImageProcessing/Upsampling.cxx
#include "itkImage.h"
#include <iostream>
int main( int argc, char * argv[] )
{
if( argc != 5 )
{
std::cerr << "Usage: "
<< std::endl
<< argv[0]
<< " inputImageFile outputImageFile nNewWidth nNewHeight"
<< std::endl;
return EXIT_FAILURE;
}
// Typedef's for pixel, image, reader and writer types
typedef unsigned char T_InputPixel;
// Doesn't work for RGB pixels
//typedef unsigned char T_OutputPixel;
//typedef itk::CovariantVector<unsigned char, 3> T_InputPixel;
//typedef itk::CovariantVector<unsigned char, 3> T_OutputPixel;
typedef itk::Image<T_InputPixel, 2> T_Image;
typedef unsigned char T_WritePixel;
typedef itk::Image<T_WritePixel, 2> T_WriteImage;
// Typedefs for the different (numerous!) elements of the "resampling"
// Identity transform.
// We don't want any transform on our image except rescaling which is not
// specified by a transform but by the input/output spacing as we will see
// later.
// So no transform will be specified.
T_Transform;
// If ITK resampler determines there is something to interpolate which is
// usually the case when upscaling (!) then we must specify the interpolation
// algorithm. In our case, we want bicubic interpolation. One way to implement
// it is with a third order b-spline. So the type is specified here and the
// order will be specified with a method call later on.
T_Interpolator;
// The resampler type itself.
T_ResampleFilter;
// Prepare the reader and update it right away to know the sizes beforehand.
T_Reader::Pointer pReader = T_Reader::New();
pReader->SetFileName( argv[1] );
pReader->Update();
// Prepare the resampler.
// Instantiate the transform and specify it should be the id transform.
T_Transform::Pointer _pTransform = T_Transform::New();
_pTransform->SetIdentity();
// Instantiate the b-spline interpolator and set it as the third order
// for bicubic.
T_Interpolator::Pointer _pInterpolator = T_Interpolator::New();
_pInterpolator->SetSplineOrder(3);
// Instantiate the resampler. Wire in the transform and the interpolator.
T_ResampleFilter::Pointer _pResizeFilter = T_ResampleFilter::New();
_pResizeFilter->SetTransform(_pTransform);
_pResizeFilter->SetInterpolator(_pInterpolator);
// Set the output origin. You may shift the original image "inside" the
// new image size by specifying something else than 0.0, 0.0 here.
const double vfOutputOrigin[2] = { 0.0, 0.0 };
_pResizeFilter->SetOutputOrigin(vfOutputOrigin);
// Compute and set the output spacing
// Compute the output spacing from input spacing and old and new sizes.
//
// The computation must be so that the following holds:
//
// new width old x spacing
// ---------- = ---------------
// old width new x spacing
//
//
// new height old y spacing
// ------------ = ---------------
// old height new y spacing
//
// So either we specify new height and width and compute new spacings (as
// we do here) or we specify new spacing and compute new height and width
// and computations that follows need to be modified a little (as it is
// done at step 2 there:
// http://itk.org/Wiki/ITK/Examples/DICOM/ResampleDICOM)
//
unsigned int nNewWidth = atoi(argv[3]);
unsigned int nNewHeight = atoi(argv[4]);
// Fetch original image size.
const T_Image::RegionType& inputRegion =
pReader->GetOutput()->GetLargestPossibleRegion();
const T_Image::SizeType& vnInputSize = inputRegion.GetSize();
unsigned int nOldWidth = vnInputSize[0];
unsigned int nOldHeight = vnInputSize[1];
// Fetch original image spacing.
const T_Image::SpacingType& vfInputSpacing =
pReader->GetOutput()->GetSpacing();
// Will be {1.0, 1.0} in the usual
// case.
double vfOutputSpacing[2];
vfOutputSpacing[0] = vfInputSpacing[0] * (double) nOldWidth / (double) nNewWidth;
vfOutputSpacing[1] = vfInputSpacing[1] * (double) nOldHeight / (double) nNewHeight;
// Set the output spacing. If you comment out the following line, the original
// image will be simply put in the upper left corner of the new image without
// any scaling.
_pResizeFilter->SetOutputSpacing(vfOutputSpacing);
// Set the output size as specified on the command line.
itk::Size<2> vnOutputSize = { {nNewWidth, nNewHeight} };
_pResizeFilter->SetSize(vnOutputSize);
// Specify the input.
_pResizeFilter->SetInput(pReader->GetOutput());
// Write the result
T_Writer::Pointer pWriter = T_Writer::New();
pWriter->SetFileName(argv[2]);
pWriter->SetInput(_pResizeFilter->GetOutput());
pWriter->Update();
return EXIT_SUCCESS;
}