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
|
/*=========================================================================
Program: GDCM (Grassroots DICOM). A DICOM library
Copyright (c) 2006-2011 Mathieu Malaterre
All rights reserved.
See Copyright.txt or http://gdcm.sourceforge.net/Copyright.html 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.
=========================================================================*/
#include "gdcmReader.h"
#include "gdcmFileMetaInformation.h"
#include "gdcmFile.h"
#include "gdcmTesting.h"
#include "gdcmMediaStorage.h"
#include "gdcmSystem.h"
#include "gdcmDirectory.h"
static int TestReadUpToTag2(const char* filename, bool verbose = false)
{
if( verbose )
std::cout << "TestRead: " << filename << std::endl;
std::ifstream is( filename, std::ios::binary );
gdcm::Reader reader;
reader.SetStream( is );
// Let's read up to end
gdcm::Tag maxDataSetTag(0xffff,0xffff);
gdcm::Tag pixeldata (0x7fe0,0x0010);
std::set<gdcm::Tag> skiptags;
// ... but do not read Pixel Data attribute (to skip mem allocation)
skiptags.insert( pixeldata );
if ( !reader.ReadUpToTag( maxDataSetTag, skiptags) )
{
std::cerr << "TestReadError: Failed to read: " << filename << std::endl;
return 1;
}
std::streamoff outStreamOffset = is.tellg();
bool iseof = is.eof();
//commenting out for warning avoidance
//const gdcm::FileMetaInformation &h = reader.GetFile().GetHeader();
const gdcm::DataSet &ds = reader.GetFile().GetDataSet();
//std::cout << ds << std::endl;
if( ds.FindDataElement( pixeldata ) )
{
std::cerr << "Found Pixel Data for: " << filename << std::endl;
return 1;
}
if( !iseof )
{
std::cerr << "Lost in out of space: " << outStreamOffset << " -> " << filename << std::endl;
return 1;
}
if(verbose)
std::cout << "{ \"" << filename << "\"," << outStreamOffset << " }," << std::endl;
const std::streamoff refoffset = -1;
if( refoffset != outStreamOffset )
{
std::cerr << filename << ": " << outStreamOffset << " " << iseof << " should be " << refoffset << std::endl;
return 1;
}
is.close();
return 0;
}
static int TestReadUpToTag2Extra()
{
const char *extradataroot = gdcm::Testing::GetDataExtraRoot();
if( !extradataroot )
{
return 1;
}
if( !gdcm::System::FileIsDirectory(extradataroot) )
{
std::cerr << "No such directory: " << extradataroot << std::endl;
return 1;
}
gdcm::Directory d;
unsigned int nfiles = d.Load( extradataroot, true ); // no recursion
std::cout << "done retrieving file list. " << nfiles << " files found." << std::endl;
gdcm::Directory::FilenamesType const & fns = d.GetFilenames();
int r = 0;
for( gdcm::Directory::FilenamesType::const_iterator it = fns.begin();
it != fns.end(); ++it )
{
const char *filename = it->c_str();
r += TestReadUpToTag2( filename );
}
return r;
}
int TestReaderUpToTag2(int argc, char *argv[])
{
if( argc == 2 )
{
const char *filename = argv[1];
return TestReadUpToTag2(filename, true);
}
// else
gdcm::Trace::DebugOff();
gdcm::Trace::WarningOff();
gdcm::Trace::ErrorOff();
int r = 0, i = 0;
const char *filename;
const char * const *filenames = gdcm::Testing::GetFileNames();
while( (filename = filenames[i]) )
{
r += TestReadUpToTag2( filename );
++i;
}
// puposely discard gdcmDataExtra test, this is just an 'extra' test...
int b2 = TestReadUpToTag2Extra(); (void)b2;
return r;
}
|
