|
|
| (4 intermediate revisions by 2 users not shown) |
| Line 1: |
Line 1: |
| ==Upsampling.cxx== | | {{warning|1=The media wiki content on this page is no longer maintained. The examples presented on the https://itk.org/Wiki/* pages likely require ITK version 4.13 or earlier releases. In many cases, the examples on this page no longer conform to the best practices for modern ITK versions. |
| <source lang="cpp">
| | }} |
| /** | |
| \file Upsampling.cxx
| |
| \date 15 april 2011 | |
| \author Francis Girard
| |
|
| |
| A simple example showing how to use the ResampleImageFilter and the
| |
| BSplineInterpolateImageFunction to (up)scale an image using bicubic
| |
| interpolation.
| |
| */
| |
|
| |
|
| #include "itkImage.h"
| | [https://itk.org/ITKExamples[ITK Sphinx Examples]] |
| #include "itkImageFileReader.h"
| |
| #include "itkImageFileWriter.h"
| |
| #include "itkIdentityTransform.h"
| |
| #include "itkBSplineInterpolateImageFunction.h"
| |
| #include "itkResampleImageFilter.h"
| |
| | |
| #include <iostream>
| |
| | |
| #if defined(_MSC_VER)
| |
| #pragma warning ( disable : 4786 )
| |
| #endif
| |
| | |
| #ifdef __BORLANDC__
| |
| #define ITK_LEAN_AND_MEAN
| |
| #endif
| |
| | |
| | |
| 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;
| |
| }
| |
|
| |
| //
| |
| // -- 1 --
| |
| // Usual typedef's for pixel, image, reader and writer types
| |
| //
| |
| typedef unsigned char T_InputPixel;
| |
| typedef unsigned char T_OutputPixel;
| |
|
| |
| typedef itk::Image<T_InputPixel, 2> T_Image;
| |
| typedef itk::ImageFileReader<T_Image> T_Reader;
| |
|
| |
| typedef unsigned char T_WritePixel;
| |
| typedef itk::Image<T_WritePixel, 2> T_WriteImage;
| |
| typedef itk::ImageFileWriter<T_WriteImage> T_Writer;
| |
|
| |
| //
| |
| // -- 2 --
| |
| // 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.
| |
| typedef itk::IdentityTransform<double, 2>
| |
| 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.
| |
| typedef itk::BSplineInterpolateImageFunction<T_Image, double, double>
| |
| T_Interpolator;
| |
|
| |
| // The resampler type itself.
| |
| typedef itk::ResampleImageFilter<T_Image, T_Image>
| |
| T_ResampleFilter;
| |
|
| |
| //
| |
| // -- 3 --
| |
| // 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();
| |
|
| |
| //
| |
| // -- 4 --
| |
| // Prepare the resampler.
| |
| //
| |
|
| |
| // 4.1 Instantiate the transform and specify it should be the id transform.
| |
| T_Transform::Pointer _pTransform = T_Transform::New();
| |
| _pTransform->SetIdentity();
| |
|
| |
| // 4.2 Instantiate the b-spline interpolator and set it as the third order
| |
| // for bicubic.
| |
| T_Interpolator::Pointer _pInterpolator = T_Interpolator::New();
| |
| _pInterpolator->SetSplineOrder(3);
| |
|
| |
| // 4.3 Instantiate the resampler. Wire in the transform and the interpolator.
| |
| //
| |
| T_ResampleFilter::Pointer _pResizeFilter = T_ResampleFilter::New();
| |
| _pResizeFilter->SetTransform(_pTransform);
| |
| _pResizeFilter->SetInterpolator(_pInterpolator);
| |
|
| |
| // 4.4 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);
| |
|
| |
| // 4.5 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);
| |
|
| |
| // 4.6 Set the output size as specified on the command line.
| |
| //
| |
| itk::Size<2> vnOutputSize = {nNewWidth, nNewHeight};
| |
| _pResizeFilter->SetSize(vnOutputSize);
| |
|
| |
| // 4.7 Well, some input to our elaborate filter.
| |
| //
| |
| _pResizeFilter->SetInput(pReader->GetOutput());
| |
|
| |
| // -- 5 --
| |
| // Done ! Prepare the writer, do the wiring and "Update()" as usual.
| |
| T_Writer::Pointer pWriter = T_Writer::New();
| |
| pWriter->SetFileName(argv[2]);
| |
| pWriter->SetInput(_pResizeFilter->GetOutput());
| |
| pWriter->Update();
| |
|
| |
| // Go home.
| |
| return 0;
| |
| }
| |
| </source>
| |
| | |
| {{ITKCMakeLists|Upsampling|}}
| |
