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
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkStreamingImageFilterTest.cxx,v $
Language: C++
Date: $Date: 2003-09-10 14:30:08 $
Version: $Revision: 1.9 $
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif
#include <iostream>
#include "itkImage.h"
#include "itkImageRegionIterator.h"
#include "itkShrinkImageFilter.h"
#include "itkStreamingImageFilter.h"
int itkStreamingImageFilterTest(int, char* [] )
{
// typedefs to simplify the syntax
typedef itk::Image<short, 2> ShortImage;
// Test the creation of an image with native type
ShortImage::Pointer if2 = ShortImage::New();
// fill in an image
ShortImage::IndexType index = {{0, 0}};
ShortImage::SizeType size = {{80, 125}};
ShortImage::RegionType region;
region.SetSize( size );
region.SetIndex( index );
if2->SetLargestPossibleRegion( region );
if2->SetBufferedRegion( region );
if2->Allocate();
itk::ImageRegionIterator<ShortImage> iterator(if2, region);
short i=0;
short scalar;
for (; !iterator.IsAtEnd(); ++iterator, ++i)
{
scalar = i;
iterator.Set( scalar );
}
// Create a filter
itk::ShrinkImageFilter< ShortImage, ShortImage >::Pointer shrink;
shrink = itk::ShrinkImageFilter< ShortImage, ShortImage >::New();
shrink->SetInput( if2 );
unsigned int factors[2] = { 2, 3 };
shrink->SetShrinkFactors(factors);
// shrink->DebugOn();
itk::StreamingImageFilter<ShortImage, ShortImage>::Pointer streamer;
streamer = itk::StreamingImageFilter<ShortImage, ShortImage>::New();
streamer->SetInput( shrink->GetOutput() );
streamer->SetNumberOfStreamDivisions( 4 );
streamer->Update();
std::cout << "Input spacing: " << if2->GetSpacing()[0] << ", "
<< if2->GetSpacing()[1] << std::endl;
std::cout << "Output spacing: " << streamer->GetOutput()->GetSpacing()[0]
<< ", "
<< streamer->GetOutput()->GetSpacing()[1] << std::endl;
// Test itkGetConstReferenceMacro and itkGetObjectMacro
const unsigned int & value = streamer->GetNumberOfStreamDivisions();
std::cout << "streamer->GetNumberOfStreamDivisions(): " << value << std::endl;
streamer->GetRegionSplitter();
//SplitterType * streamer->GetRegionSplitter();
//
// The rest of this code determines whether the shrink code produced
// the image we expected.
//
ShortImage::RegionType requestedRegion;
requestedRegion = streamer->GetOutput()->GetRequestedRegion();
itk::ImageRegionIterator<ShortImage>
iterator2(streamer->GetOutput(), requestedRegion);
bool passed = true;
for (; !iterator2.IsAtEnd(); ++iterator2)
{
short trueValue = (short) (shrink->GetShrinkFactors()[0] * iterator2.GetIndex()[0])
+ (region.GetSize()[0]
* shrink->GetShrinkFactors()[1] * iterator2.GetIndex()[1]);
if ( iterator2.Get() != trueValue )
{
passed = false;
std::cout << "Pixel " << iterator2.GetIndex() << " is incorrect" << std::endl;
}
}
if (passed)
{
std::cout << "ImageStreamingFilter test passed." << std::endl;
return EXIT_SUCCESS;
}
else
{
std::cout << "ImageStreaming Filter test failed." << std::endl;
return EXIT_FAILURE;
}
}
|
