DGtal  1.3.beta
Functions
testITKImage.cpp File Reference
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/kernel/SpaceND.h"
#include "DGtal/kernel/domains/HyperRectDomain.h"
#include "DGtal/images/ImageContainerByITKImage.h"
#include <boost/foreach.hpp>
#include <itkExtractImageFilter.h>
Include dependency graph for testITKImage.cpp:

Go to the source code of this file.

Functions

bool testITKImage ()
 
bool testITKMethod ()
 
bool testITKImageWithMetadata ()
 
bool testITKImageWithShiftDomain ()
 
int main (int argc, char **argv)
 

Detailed Description

This program is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR a PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

Author
David Coeurjolly (david.nosp@m..coe.nosp@m.urjol.nosp@m.ly@l.nosp@m.iris..nosp@m.cnrs.nosp@m..fr ) Laboratoire d'InfoRmatique en Image et Systèmes d'information - LIRIS (CNRS, UMR 5205), CNRS, France
Date
2010/12/09

Functions for testing class ITKImage.

This file is part of the DGtal library.

Definition in file testITKImage.cpp.

Function Documentation

◆ main()

int main ( int  argc,
char **  argv 
)

Definition at line 411 of file testITKImage.cpp.

412 {
413  trace.beginBlock ( "Testing class ITKImage" );
414  trace.info() << "Args:";
415  for ( int i = 0; i < argc; ++i )
416  trace.info() << " " << argv[ i ];
417  trace.info() << endl;
418 
419  bool res = testITKImage() && testITKMethod() &&
421  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
422  trace.endBlock();
423  return res ? 0 : 1;
424 }

References DGtal::Trace::beginBlock(), DGtal::Trace::emphase(), DGtal::Trace::endBlock(), DGtal::Trace::info(), testITKImage(), testITKImageWithMetadata(), testITKImageWithShiftDomain(), testITKMethod(), and DGtal::trace.

◆ testITKImage()

bool testITKImage ( )

Example of a test. To be completed.

Definition at line 67 of file testITKImage.cpp.

68 {
69  unsigned int nbok = 0;
70  unsigned int nb = 0;
71 
72  trace.beginBlock ( "ITK Image init..." );
73 
74  typedef DGtal::int32_t Integer;
75  typedef SpaceND<3, Integer > Space3Type;
77  typedef Domain::Point Point;
78 
80 
81  const Integer t[ ] = { 1, 1, 1};
82  const Integer t2[ ] = { 5, 5, 5};
83  const Integer t3[ ] = { 2, 2, 2};
84  Point lowerBound ( t );
85  Point upperBound ( t2 );
86  Point c ( t3 );
87  Integer val;
88 
89  Image myImage ( Domain(lowerBound, upperBound) );
90 
91  trace.info() << myImage << std::endl;
92  trace.info() << "getvalue= " << myImage(c) << endl;
93  trace.info() << "set value 23 " << endl;
94  myImage.setValue( c, 23);
95 
96  val = myImage(c);
97 
98  if (val == 23)
99  nbok++;
100  trace.info() << "getvalue= " << val << endl;
101  nb++;
102 
103  //Iterator test
104  trace.info() << "Simple Iterator=";
105  for (Image::ConstIterator it = myImage.begin(), itend = myImage.end();
106  it != itend;
107  ++it)
108  trace.warning() << myImage(it) << " ";
109  trace.info() << endl;
110 
111  //We rewrite the image
112  int nbVal = 0;
113  for (Image::Iterator it = myImage.begin(), itend = myImage.end();
114  it != itend;
115  ++it)
116  myImage.setValue(it, nbVal++);
117 
118  trace.info() << "Set Iterator=";
119  for (Image::ConstIterator it = myImage.begin(), itend = myImage.end();
120  it != itend;
121  ++it)
122  trace.warning() << myImage(it) << " ";
123  trace.info() << endl;
124 
125  auto & container = myImage.container();
126  (void)container;
127 
128 
129  trace.info() << "(" << nbok << "/" << nb << ") "
130  << "true == true" << std::endl;
131  trace.endBlock();
132 
133  return nbok == nb;
134 }

References DGtal::Trace::beginBlock(), DGtal::Trace::endBlock(), DGtal::Trace::info(), DGtal::Image< TImageContainer >::setValue(), DGtal::trace, and DGtal::Trace::warning().

Referenced by main().

◆ testITKImageWithMetadata()

bool testITKImageWithMetadata ( )

Definition at line 216 of file testITKImage.cpp.

217 {
218  trace.beginBlock ( "ITK Image With Metadata init..." );
219 
220  unsigned int nbok = 0;
221  unsigned int nb = 0;
222 
223  typedef DGtal::int32_t Integer;
224  typedef SpaceND<3, Integer > Space3Type;
226  typedef Domain::Point Point;
227 
229 
230  // REMEMBER: Origin in ITK is the physical position of the index {0,0,...}
231  // even when that zero-index is not included in the LargestPossibleRegion of the image.
232  const Integer t[ ] = { 1, 1, 1};
233  const Integer t2[ ] = { 5, 5, 5};
234  const Integer t3[ ] = { 2, 2, 2};
235  Point lowerBound ( t );
236  Point upperBound ( t2 );
237  Point c ( t3 );
238  Integer val;
239 
240  Image myImage ( Domain(lowerBound, upperBound) );
241  // Fill the image
242  int nbVal = 0;
243  for (Image::Iterator it = myImage.begin(), itend = myImage.end();
244  it != itend; ++it) {
245  myImage.setValue(it, nbVal++);
246  }
247  // Change the default metadata (physical properties) of an itk_image.
248  auto itk_image = myImage.getITKImagePointer();
249  Image::ITKImage::PointType origin;
250  origin.Fill(10.0);
251  Image::ITKImage::SpacingType spacing;
252  spacing.Fill(2.0);
253  Image::ITKImage::DirectionType direction;
254  direction.SetIdentity();
255  itk_image->SetOrigin(origin);
256  itk_image->SetSpacing(spacing);
257  itk_image->SetDirection(direction);
258 
259  // Check that the value of index points is not affected by a change of metadata.
260  val = myImage.operator()(lowerBound);
261  nbok += (val == 0); nb++;
262  trace.info() << "Index: " << lowerBound << ". Value: " << val << ". Expected: " << 0 << std::endl;
263  val = myImage.operator()(upperBound);
264  nbok += (val == 124); nb++;
265  trace.info() << "Index: " << upperBound << ". Value: " << val << ". Expected: " << 124 << std::endl;
266  val = myImage.operator()(c);
267  nbok += (val == 31); nb++;
268  trace.info() << "Index: " << c << ". Value: " << val << ". Expected: " << 31 << std::endl;
269 
270  Image::PhysicalPoint physical_point;
271  Image::PhysicalPoint expected_physical_point;
272  Image::Point index_point;
273 
274  // when shiftDomain is zero, index points (ITK) and domain points (DGtal) are equal
275  index_point = myImage.getIndexFromDomainPoint(lowerBound);
276  nbok += (index_point == lowerBound); nb++;
277  physical_point = myImage.getPhysicalPointFromDomainPoint(index_point);
278  expected_physical_point = myImage.getLowerBoundAsPhysicalPoint();
279  nbok += (physical_point[0] == 12.0); nb++;
280  nbok += (physical_point == expected_physical_point); nb++;
281  trace.info() << "Index: " << index_point << ". PhysicalPoint: " << physical_point << ". Expected: " << expected_physical_point << std::endl;
282 
283  index_point = myImage.getIndexFromDomainPoint(upperBound);
284  nbok += (index_point == upperBound); nb++;
285  physical_point = myImage.getPhysicalPointFromDomainPoint(index_point);
286  expected_physical_point = myImage.getUpperBoundAsPhysicalPoint();
287  nbok += (physical_point[0] == 20.0); nb++;
288  nbok += (physical_point == expected_physical_point); nb++;
289  trace.info() << "Index: " << index_point << ". PhysicalPoint: " << physical_point << ". Expected: " << expected_physical_point << std::endl;
290 
291  auto index_back = myImage.getDomainPointFromPhysicalPoint(physical_point);
292  nbok += (index_back == upperBound); nb++;
293  trace.info() << "PhysicalPoint: " << physical_point << ". Index (back): " << index_back << ". Expected: " << upperBound << std::endl;
294 
295  trace.endBlock();
296 
297  return nbok == 10 && nb == 10;
298 }

References DGtal::Trace::beginBlock(), DGtal::Trace::endBlock(), DGtal::Trace::info(), DGtal::Image< TImageContainer >::setValue(), and DGtal::trace.

Referenced by main().

◆ testITKImageWithShiftDomain()

bool testITKImageWithShiftDomain ( )

Definition at line 300 of file testITKImage.cpp.

301 {
302  trace.beginBlock ( "ITK Image With ShiftDomain init..." );
303 
304  unsigned int nbok = 0;
305  unsigned int nb = 0;
306 
307  typedef DGtal::int32_t Integer;
308  typedef SpaceND<3, Integer > Space3Type;
310  typedef Domain::Point Point;
311 
313 
314  const Integer t[ ] = { 1, 1, 1};
315  const Integer t2[ ] = { 5, 5, 5};
316  const Integer t3[ ] = { 2, 2, 2};
317  Point lowerBound ( t );
318  Point upperBound ( t2 );
319  Point c ( t3 );
320  Integer val;
321 
322  Image myImage ( Domain(lowerBound, upperBound) );
323 
324  // Fill the image
325  int nbVal = 0;
326  for (Image::Iterator it = myImage.begin(), itend = myImage.end();
327  it != itend; ++it) {
328  myImage.setValue(it, nbVal++);
329  }
330 
331  // Change the default metadata (physical properties) of an itk_image.
332  auto itk_image = myImage.getITKImagePointer();
333  Image::ITKImage::PointType origin;
334  origin.Fill(10.0);
335  Image::ITKImage::SpacingType spacing;
336  spacing.Fill(2.0);
337  Image::ITKImage::DirectionType direction;
338  direction.SetIdentity();
339  itk_image->SetOrigin(origin);
340  itk_image->SetSpacing(spacing);
341  itk_image->SetDirection(direction);
342 
343 
344  // Apply a domainShift
345  const Integer sd[ ] = { -20, -20, -20};
346  Point domainShift(sd);
347  myImage.updateDomain(domainShift);
348  Point new_lowerBound = myImage.domain().lowerBound();
349  Point new_upperBound = myImage.domain().upperBound();
350  nbok += ( new_lowerBound == lowerBound + domainShift); nb++;
351  trace.info() << "lowerBound: " << new_lowerBound << ". Expected: " << lowerBound + domainShift << std::endl;
352  nbok += ( new_upperBound == upperBound + domainShift); nb++;
353  trace.info() << "upperBound: " << new_upperBound << ". Expected: " << upperBound + domainShift << std::endl;
354 
355  // Check that the shifted domain points to the correct indices of the image.
356  val = myImage.operator()(new_lowerBound);
357  // It should have the same value than lowerBound had before applying the domainShift
358  nbok += (val == 0); nb++;
359  trace.info() << "Index: " << new_lowerBound << ". Value: " << val << ". Expected: " << 0 << std::endl;
360  val = myImage.operator()(new_upperBound);
361  nbok += (val == 124); nb++;
362  trace.info() << "Index: " << new_upperBound << ". Value: " << val << ". Expected: " << 124 << std::endl;
363  val = myImage.operator()(myImage.getDomainPointFromIndex(c));
364  nbok += (val == 31); nb++;
365  trace.info() << "Index: " << c << ". Value: " << val << ". Expected: " << 31 << std::endl;
366 
367  Image::PhysicalPoint physical_point;
368  Image::PhysicalPoint expected_physical_point;
369  Image::Point index_point;
370  Image::Point domain_point;
371 
372  index_point = lowerBound;
373  domain_point = new_lowerBound;
374  physical_point = myImage.getPhysicalPointFromDomainPoint(domain_point);
375  expected_physical_point = Image::PhysicalPoint(12.0, 12.0, 12.0);
376  nbok += (physical_point == expected_physical_point); nb++;
377  trace.info() << "Domain: " << domain_point <<
378  ". Index: " << index_point <<
379  ". PhysicalPoint: " << physical_point <<
380  ". Expected: " << expected_physical_point << std::endl;
381 
382  index_point = myImage.getIndexFromDomainPoint(new_lowerBound);
383  nbok += ( index_point == lowerBound ); nb++;
384  trace.info() << "index_point: " << index_point << ". Expected: " << lowerBound << std::endl;
385  physical_point = myImage.getPhysicalPointFromDomainPoint(new_lowerBound);
386  expected_physical_point = myImage.getLowerBoundAsPhysicalPoint();
387  nbok += (physical_point[0] == 12.0); nb++;
388  nbok += (physical_point == expected_physical_point); nb++;
389  trace.info() << "Domain: " << new_lowerBound <<
390  ". Index: " << index_point <<
391  ". PhysicalPoint: " << physical_point <<
392  ". Expected: " << expected_physical_point << std::endl;
393 
394  index_point = upperBound;
395  domain_point = new_upperBound;
396  physical_point = myImage.getPhysicalPointFromDomainPoint(domain_point);
397  expected_physical_point = myImage.getUpperBoundAsPhysicalPoint();
398  nbok += (physical_point[0] == 20.0); nb++;
399  nbok += (physical_point == expected_physical_point); nb++;
400  trace.info() << "Domain: " << new_lowerBound <<
401  ". Index: " << index_point <<
402  ". PhysicalPoint: " << physical_point <<
403  ". Expected: " << expected_physical_point << std::endl;
404 
405  return nbok == 11 && nb == 11;
406 }

References DGtal::Trace::beginBlock(), DGtal::Image< TImageContainer >::domain(), DGtal::Trace::info(), DGtal::Image< TImageContainer >::setValue(), and DGtal::trace.

Referenced by main().

◆ testITKMethod()

bool testITKMethod ( )

Definition at line 136 of file testITKImage.cpp.

137 {
138  unsigned int nbok = 0;
139  unsigned int nb = 0;
140 
141  trace.beginBlock ( "Test the use of a pure ITK method..." );
142 
143  typedef DGtal::int32_t Integer;
146  typedef Domain::Point Point;
147 
148 
150 
151  Point lowerBound ( 0, 0 );
152  Point upperBound ( 9, 9);
153  Domain domain(lowerBound, upperBound);
154 
155  Image myImage(domain);
156  trace.info() << myImage << std::endl;
157 
158  //We fill the image
159  Integer nbVal = 0;
160  for (Image::Iterator it = myImage.begin(), itend = myImage.end();
161  it != itend;
162  ++it)
163  myImage.setValue(it, nbVal++);
164 
165  trace.info() << "Input image=";
166  for (Image::ConstIterator it = myImage.begin(), itend = myImage.end();
167  it != itend;
168  ++it)
169  trace.warning() << myImage(it) << " ";
170  trace.info() << endl;
171 
172 
173  // We define a cropFilter
174  typedef itk::ExtractImageFilter< Image::ITKImage, Image::ITKImage > CropFilter;
175 
176  // Crop filter region
177  Image::ITKImage::SizeType size;
178  size[0] = 5;
179  size[1] = 5;
180 
181  Image::ITKImage::IndexType index;
182  index[0] = 2;
183  index[1] = 2;
184 
185  Image::ITKImage::RegionType regionToExtract(index,size);
186 
187  // Crop filter process
188  CropFilter::Pointer cropFilter = CropFilter::New();
189  cropFilter->SetInput( myImage.getITKImagePointer() );
190  cropFilter->SetExtractionRegion( regionToExtract );
191  cropFilter->Update();
192 
193  // Pointer to the filter output
194  Image::ITKImagePointer handleOut = cropFilter->GetOutput();
195  Image myImageOut ( handleOut );
196 
197 
198  trace.info() << "Output image=";
199 
200  for (Image::ConstIterator it = myImageOut.begin(), itend = myImageOut.end();
201  it != itend;
202  ++it)
203  {
204  nbok += (it.Value() == (it.GetIndex()[1]*10 + it.GetIndex()[0]));
205  nb++;
206  trace.warning() << it.Value() << "(" << (it.GetIndex()[1]*10 + it.GetIndex()[0]) << ")" << " ";
207  }
208  trace.info() << endl;
209 
210  trace.info() << "(" << nbok << "/" << nb << ") " << "true == true" << std::endl;
211  trace.endBlock();
212 
213  return nbok == 25 && nb == 25;
214 }

References DGtal::Trace::beginBlock(), domain, DGtal::Trace::endBlock(), index(), DGtal::Trace::info(), DGtal::Image< TImageContainer >::setValue(), DGtal::trace, and DGtal::Trace::warning().

Referenced by main().

DGtal::HyperRectDomain
Aim: Parallelepidec region of a digital space, model of a 'CDomain'.
Definition: HyperRectDomain.h:99
DGtal::Trace::endBlock
double endBlock()
DGtal::ImageContainerByITKImage
Aim: implements a model of CImageContainer using a ITK Image.
Definition: ImageContainerByITKImage.h:93
DGtal::int32_t
boost::int32_t int32_t
signed 32-bit integer.
Definition: BasicTypes.h:72
DGtal::Trace::emphase
std::ostream & emphase()
DGtal::trace
Trace trace
Definition: Common.h:154
DGtal::ImageContainerBySTLVector< Domain, Value >::ConstIterator
std::vector< Value >::const_iterator ConstIterator
Definition: ImageContainerBySTLVector.h:265
DGtal::Trace::beginBlock
void beginBlock(const std::string &keyword="")
DGtal::SpaceND
Definition: SpaceND.h:95
DGtal::Trace::info
std::ostream & info()
Image
ImageContainerBySTLVector< Domain, Value > Image
Definition: testSimpleRandomAccessRangeFromPoint.cpp:45
Domain
HyperRectDomain< Space > Domain
Definition: testSimpleRandomAccessRangeFromPoint.cpp:44
index
unsigned int index(DGtal::uint32_t n, unsigned int b)
Definition: testBits.cpp:44
DGtal::ImageContainerBySTLVector< Domain, Value >::Iterator
std::vector< Value >::iterator Iterator
Definition: ImageContainerBySTLVector.h:264
Integer
Point::Coordinate Integer
Definition: examplePlaneProbingParallelepipedEstimator.cpp:44
DGtal::Image
Aim: implements association bewteen points lying in a digital domain and values.
Definition: Image.h:69
domain
Domain domain
Definition: testProjection.cpp:88
DGtal::PointVector< dim, Integer >
testITKMethod
bool testITKMethod()
Definition: testITKImage.cpp:136
Space2Type
SpaceND< 2, int > Space2Type
Definition: testFreemanChain.cpp:51
testITKImageWithShiftDomain
bool testITKImageWithShiftDomain()
Definition: testITKImage.cpp:300
DGtal::Trace::warning
std::ostream & warning()
testITKImage
bool testITKImage()
Definition: testITKImage.cpp:67
Point
MyPointD Point
Definition: testClone2.cpp:383
testITKImageWithMetadata
bool testITKImageWithMetadata()
Definition: testITKImage.cpp:216