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
|
/*=========================================================================
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 "gdcmImageHelper.h"
#include "gdcmMediaStorage.h"
#include "gdcmDataSet.h"
#include "gdcmAttribute.h"
#include "gdcmDirectionCosines.h"
#include "gdcmImage.h"
#include "gdcmFile.h"
int TestImageHelper3(int, char *[])
{
std::vector<double> impos;
impos.resize(3);
impos[0] = 1;
impos[1] = 2;
impos[2] = 3;
std::vector<double> spacing;
spacing.resize(3);
spacing[0] = 0.6;
spacing[1] = 0.5;
spacing[2] = 0.4;
std::vector<double> dircos;
dircos.resize(6);
dircos[0] = 0;
dircos[1] = 1;
dircos[2] = 0;
dircos[3] = 1;
dircos[4] = 0;
dircos[5] = 0;
gdcm::Image img;
img.SetNumberOfDimensions(3);
img.SetOrigin( impos.data() );
img.SetSpacing( spacing.data() );
img.SetDirectionCosines( dircos.data() );
// Try SC
gdcm::File file;
{
gdcm::ImageHelper::SetSecondaryCaptureImagePlaneModule( false );
gdcm::MediaStorage ms( gdcm::MediaStorage::SecondaryCaptureImageStorage );
gdcm::DataSet&ds = file.GetDataSet();
gdcm::DataElement de( gdcm::Tag(0x0008, 0x0016) );
const char* msstr = gdcm::MediaStorage::GetMSString(ms);
de.SetByteValue( msstr, (uint32_t)strlen(msstr) );
de.SetVR( gdcm::Attribute<0x0008, 0x0016>::GetVR() );
ds.Insert( de );
gdcm::ImageHelper::SetDirectionCosinesValue( ds, dircos );
gdcm::ImageHelper::SetOriginValue( ds, img );
gdcm::ImageHelper::SetSpacingValue( ds, spacing );
std::cout << ds << std::endl;
gdcm::Tag nominal(0x0018,0x2010 );
if( !ds.FindDataElement( nominal ) ) return 1;
gdcm::Tag iop(0x0020,0x0037);
if( ds.FindDataElement( iop ) ) return 1;
std::vector<double> ret = gdcm::ImageHelper::GetSpacingValue( file );
// warning;: only two dim:
if( ret[0] != spacing[0] || ret[1] != spacing[1] ) {
std::cerr << ret[0] << "," << ret[1] << std::endl;
return 1;
}
}
gdcm::File cp2330file;
{
// New CP 2330 behavior
gdcm::ImageHelper::SetSecondaryCaptureImagePlaneModule( true );
gdcm::MediaStorage ms( gdcm::MediaStorage::SecondaryCaptureImageStorage );
gdcm::DataSet&ds = cp2330file.GetDataSet();
gdcm::DataElement de( gdcm::Tag(0x0008, 0x0016) );
const char* msstr = gdcm::MediaStorage::GetMSString(ms);
de.SetByteValue( msstr, (uint32_t)strlen(msstr) );
de.SetVR( gdcm::Attribute<0x0008, 0x0016>::GetVR() );
ds.Insert( de );
// New CP 2330 behavior
gdcm::ImageHelper::SetDirectionCosinesValue( ds, dircos );
gdcm::ImageHelper::SetOriginValue( ds, img );
gdcm::ImageHelper::SetSpacingValue( ds, spacing );
std::cout << ds << std::endl;
// previous call added the attribute
gdcm::Tag pixelspacing(0x0028,0x0030);
if( !ds.FindDataElement( pixelspacing ) ) return 1;
gdcm::Tag iop(0x0020,0x0037);
if( !ds.FindDataElement( iop ) ) return 1;
std::vector<double> ret = gdcm::ImageHelper::GetSpacingValue( cp2330file );
if( ret != spacing ) {
std::cerr << ret[0] << "," << ret[1] << "," << ret[2] << std::endl;
return 1;
}
// make sure legacy still works:
ret = gdcm::ImageHelper::GetSpacingValue( file );
if( ret[0] != spacing[0] || ret[1] != spacing[1] || ret[2] != 1 ) {
std::cerr << ret[0] << "," << ret[1] << std::endl;
return 1;
}
}
std::cout << "success" << std::endl;
return 0;
}
|
