ITK  4.13.0
Insight Segmentation and Registration Toolkit
WikiExamples/ImageProcessing/ScaleTransform.cxx
#include "itkImage.h"
typedef itk::Image<unsigned char, 2> ImageType;
static void CreateImage(ImageType::Pointer image);
int main(int, char *[])
{
ImageType::Pointer image = ImageType::New();
CreateImage(image);
typedef itk::ImageFileWriter<ImageType> WriterType;
WriterType::Pointer inputWriter = WriterType::New();
inputWriter->SetFileName("input.png");
inputWriter->SetInput(image);
inputWriter->Update();
// typedef itk::ScaleTransform<float, 2> TransformType; // If you want to use float here, you must use:
// typedef itk::ResampleImageFilter<ImageType, ImageType, float> ResampleImageFilterType; later.
typedef itk::ScaleTransform<double, 2> TransformType;
TransformType::Pointer scaleTransform = TransformType::New();
scale[0] = 1.5; // newWidth/oldWidth
scale[1] = 1.5;
scaleTransform->SetScale(scale);
center[0] = image->GetLargestPossibleRegion().GetSize()[0]/2;
center[1] = image->GetLargestPossibleRegion().GetSize()[1]/2;
scaleTransform->SetCenter(center);
typedef itk::ResampleImageFilter<ImageType, ImageType> ResampleImageFilterType;
ResampleImageFilterType::Pointer resampleFilter = ResampleImageFilterType::New();
resampleFilter->SetTransform(scaleTransform);
resampleFilter->SetInput(image);
resampleFilter->SetSize(image->GetLargestPossibleRegion().GetSize());
resampleFilter->Update();
WriterType::Pointer outputWriter = WriterType::New();
outputWriter->SetFileName("output.png");
outputWriter->SetInput(resampleFilter->GetOutput());
outputWriter->Update();
return EXIT_SUCCESS;
}
void CreateImage(ImageType::Pointer image)
{
start.Fill(0);
size.Fill(101);
ImageType::RegionType region(start, size);
image->SetRegions(region);
image->Allocate();
image->FillBuffer(0);
// Make a white square
for(unsigned int r = 40; r < 60; r++)
{
for(unsigned int c = 40; c < 60; c++)
{
pixelIndex[0] = r;
pixelIndex[1] = c;
image->SetPixel(pixelIndex, 255);
}
}
itk::ImageRegionIterator<ImageType> imageIterator(image, image->GetLargestPossibleRegion());
// Draw a white border
while(!imageIterator.IsAtEnd())
{
if(imageIterator.GetIndex()[0] == 0 || imageIterator.GetIndex()[0] == static_cast<int>(image->GetLargestPossibleRegion().GetSize()[0]) - 1 ||
imageIterator.GetIndex()[1] == 0 || imageIterator.GetIndex()[1] == static_cast<int>(image->GetLargestPossibleRegion().GetSize()[1]) - 1)
{
imageIterator.Set(255);
}
++imageIterator;
}
}