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
|
/*=========================================================================
Program: Visualization Toolkit
Module: TestMetaIO.cxx
Copyright (c) Ken Martin, Will Schroeder, Bill Lorensen
All rights reserved.
See Copyright.txt or http://www.kitware.com/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 notice for more information.
=========================================================================*/
// .NAME Test of vtkMetaIO / MetaImage
// .SECTION Description
//
#include "vtkDataObject.h"
#include "vtkImageData.h"
#include "vtkImageMathematics.h"
#include "vtkMetaImageReader.h"
#include "vtkMetaImageWriter.h"
#include "vtkObjectFactory.h"
#include "vtkOutputWindow.h"
#include <cmath>
int TestMetaIO(int argc, char* argv[])
{
vtkOutputWindow::GetInstance()->PromptUserOn();
if (argc <= 1)
{
cout << "Usage: " << argv[0] << " <meta image file>" << endl;
return 1;
}
vtkMetaImageReader* reader = vtkMetaImageReader::New();
reader->SetFileName(argv[1]);
reader->Update();
cout << "10, 10, 10 : (1) : " << reader->GetOutput()->GetScalarComponentAsFloat(10, 10, 10, 0)
<< endl;
cout << "24, 37, 10 : (168) : " << reader->GetOutput()->GetScalarComponentAsFloat(24, 37, 10, 0)
<< endl;
vtkMetaImageWriter* writer = vtkMetaImageWriter::New();
writer->SetFileName("TestMetaIO.mha");
writer->SetInputConnection(reader->GetOutputPort());
writer->Write();
reader->Delete();
writer->Delete();
vtkMetaImageReader* readerStd = vtkMetaImageReader::New();
readerStd->SetFileName(argv[1]);
readerStd->Update();
vtkMetaImageReader* readerNew = vtkMetaImageReader::New();
readerNew->SetFileName("TestMetaIO.mha");
readerNew->Update();
double error = 0;
int* ext = readerStd->GetOutput()->GetExtent();
for (int z = ext[4]; z <= ext[5]; z += 2)
{
for (int y = ext[2]; y <= ext[3]; y++)
{
for (int x = ext[0]; x <= ext[1]; x++)
{
error += fabs(readerStd->GetOutput()->GetScalarComponentAsFloat(x, y, z, 0) -
readerNew->GetOutput()->GetScalarComponentAsFloat(x, y, z, 0));
}
}
}
readerStd->Delete();
readerNew->Delete();
if (error > 1)
{
cerr << "Error: Image difference on read/write = " << error << endl;
return 1;
}
cout << "Success! Error = " << error << endl;
return 0;
}
|
