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
118
119
120
121
122
123
124
125
126
127
128
129
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkNrrdDiffusionTensor3DImageReadTest.cxx,v $
Language: C++
Date: $Date: 2009-04-05 23:09:36 $
Version: $Revision: 1.3 $
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 <fstream>
#include "itkImageFileReader.h"
#include "itkImage.h"
#include "itkDiffusionTensor3D.h"
int itkNrrdDiffusionTensor3DImageReadTest( int ac, char* av[] )
{
if(ac < 1)
{
std::cerr << "Usage: " << av[0] << " Input\n";
return EXIT_FAILURE;
}
typedef itk::DiffusionTensor3D<float> PixelType;
typedef itk::Image<PixelType, 3> myImage;
itk::ImageFileReader<myImage>::Pointer reader
= itk::ImageFileReader<myImage>::New();
reader->SetFileName(av[1]);
try
{
reader->Update();
}
catch (itk::ExceptionObject & e)
{
std::cerr << "exception in file reader " << std::endl;
std::cerr << e << std::endl;
return EXIT_FAILURE;
}
myImage::Pointer image = reader->GetOutput();
// Pick off some pixels in the test image (the first, the last, and
// a few in between) and make sure that the values are very close to
// what we know to be correct, given that we know exactly what
// volume is being read. I/O errors will generate more than tiny
// differences near representational precision.
float err = 0;
myImage::IndexType coord;
PixelType sample;
coord[0] = 0;
coord[1] = 0;
coord[2] = 0;
sample = image->GetPixel(coord);
err += vcl_fabs(sample(0,0) - 4.0248222);
err += vcl_fabs(sample(0,1) - -0.2367909);
err += vcl_fabs(sample(0,2) - 0.23370844);
err += vcl_fabs(sample(1,1) - 1.2593846);
err += vcl_fabs(sample(1,2) - -0.042955428);
err += vcl_fabs(sample(2,2) - -0.15239859);
coord[0] = 4;
coord[1] = 5;
coord[2] = 6;
sample = image->GetPixel(coord);
err += vcl_fabs(sample(0,0) - 6.3746634);
err += vcl_fabs(sample(0,1) - -0.10792637);
err += vcl_fabs(sample(0,2) - 0.078167915);
err += vcl_fabs(sample(1,1) - 6.2988205);
err += vcl_fabs(sample(1,2) - 0.40197921);
err += vcl_fabs(sample(2,2) - 4.8993769);
coord[0] = 3;
coord[1] = 3;
coord[2] = 3;
sample = image->GetPixel(coord);
err += vcl_fabs(sample(0,0) - 4.9852095);
err += vcl_fabs(sample(0,1) - 0.51356757);
err += vcl_fabs(sample(0,2) - -1.1457335);
err += vcl_fabs(sample(1,1) - 0.6538533);
err += vcl_fabs(sample(1,2) - -0.19546235);
err += vcl_fabs(sample(2,2) - 1.3520061);
coord[0] = 3;
coord[1] = 0;
coord[2] = 5;
sample = image->GetPixel(coord);
err += vcl_fabs(sample(0,0) - 6.1530757);
err += vcl_fabs(sample(0,1) - 0.40738893);
err += vcl_fabs(sample(0,2) - 0.11354899);
err += vcl_fabs(sample(1,1) - 5.968379);
err += vcl_fabs(sample(1,2) - -0.22650599);
err += vcl_fabs(sample(2,2) - 5.0935121);
coord[0] = 0;
coord[1] = 3;
coord[2] = 6;
sample = image->GetPixel(coord);
err += vcl_fabs(sample(0,0) - 7.3227096);
err += vcl_fabs(sample(0,1) - -0.34219909);
err += vcl_fabs(sample(0,2) - -0.2011447);
err += vcl_fabs(sample(1,1) - 5.6443777);
err += vcl_fabs(sample(1,2) - 0.43205333);
err += vcl_fabs(sample(2,2) - 5.3755102);
double thresh = 0.00000041;
if (err > thresh)
{
std::cout << "failure because err == " << err
<< "> " << thresh << std::endl;
return EXIT_FAILURE;
}
else
{
return EXIT_SUCCESS;
}
}
|
