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
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
|
/*
*
* Copyright (c) 2008-2012, OFFIS e.V. and ICSMED AG, Oldenburg, Germany
* Copyright (C) 2013-2024, J. Riesmeier, Oldenburg, Germany
* All rights reserved. See COPYRIGHT file for details.
*
* Test read and write methods of DICOM RT classes
*
* Generated manually
* File created on 2008-12-08
* Last modified on 2024-01-14 by Riesmeier
*
*/
#include "dcmtk/config/osconfig.h" /* make sure OS specific configuration is included first */
#include "dcmtk/dcmrt/drtdose.h"
#include "dcmtk/dcmrt/drtimage.h"
#include "dcmtk/dcmrt/drtplan.h"
#include "dcmtk/dcmrt/drtstrct.h"
#include "dcmtk/dcmrt/drttreat.h"
#include "dcmtk/dcmdata/dcfilefo.h"
#define TEST_COPY_CONSTRUCTOR
#define TEST_ASSIGNMENT_OPERATOR
#define RTDOSE_FILENAME "/share/dicom/images/dkfz-rt/Phantom/RD1.3.6.1.4.1.2452.6.2742071499.1301053063.852537260.1034784.dcm"
#define RTIMAGE_FILENAME "/share/dicom/images/dkfz-rt/Phantom/RI1.3.6.1.4.1.2452.6.2440400960.1328333180.1418898622.324303.dcm"
#define RTPLAN_FILENAME "/share/dicom/images/dkfz-rt/Phantom/RP1.3.6.1.4.1.2452.6.998701572.1323706291.2356805512.4049491.dcm"
#define RTSTRUCT_FILENAME "/share/dicom/images/dkfz-rt/Phantom/RS1.3.6.1.4.1.2452.6.2951369053.1166013784.715527822.1679960.dcm"
#define RTTREAT_FILENAME "/share/dicom/images/dkfz-rt/rt_set_20090707/RT_RECORD/SER01/IMG004"
int main()
{
OFCondition status;
DcmFileFormat fileformat;
OFLog::configure(OFLogger::INFO_LOG_LEVEL);
DCMRT_INFO("DCMRT Test Program" << OFendl << "------------------" << OFendl);
/* test RT dose */
DRTDoseIOD *rtdose = new DRTDoseIOD();
if (rtdose != NULL)
{
status = fileformat.loadFile(RTDOSE_FILENAME);
if (status.good())
status = rtdose->read(*fileformat.getDataset());
DCMRT_INFO("Read RT Dose: " << status.text());
if (status.good())
{
fileformat.clear();
status = rtdose->write(*fileformat.getDataset());
DCMRT_INFO("Save RT Dose: " << status.text());
if (status.good())
status = fileformat.saveFile("rtdose.dcm");
DCMRT_INFO("Write RT Dose: " << status.text());
#ifdef TEST_COPY_CONSTRUCTOR
DRTDoseIOD *rtdose2 = new DRTDoseIOD(*rtdose);
if (rtdose2 != NULL)
{
fileformat.clear();
status = rtdose2->write(*fileformat.getDataset());
DCMRT_INFO("Write copy of RT Dose: " << status.text());
if (status.good())
status = fileformat.saveFile("rtdose2.dcm");
DCMRT_INFO("Save copy of RT Dose: " << status.text());
delete rtdose2;
rtdose2 = NULL;
} else
DCMRT_ERROR("Cannot allocate memory for copy of RT Dose");
#endif
} else
DCMRT_INFO("Skipping further tests for RT Dose");
delete rtdose;
rtdose = NULL;
} else
DCMRT_ERROR("Cannot allocate memory for RT Dose");
DCMRT_INFO("" /* empty line */);
/* test RT image */
DRTImageIOD *rtimage = new DRTImageIOD();
if (rtimage != NULL)
{
status = fileformat.loadFile(RTIMAGE_FILENAME);
if (status.good())
status = rtimage->read(*fileformat.getDataset());
DCMRT_INFO("Read RT Image: " << status.text());
if (status.good())
{
fileformat.clear();
status = rtimage->write(*fileformat.getDataset());
DCMRT_INFO("Write RT Image: " << status.text());
if (status.good())
status = fileformat.saveFile("rtimage.dcm");
DCMRT_INFO("Save RT Image: " << status.text());
} else
DCMRT_INFO("Skipping further tests for RT Image");
delete rtimage;
rtimage = NULL;
} else
DCMRT_ERROR("Cannot allocate memory for RT Image");
DCMRT_INFO("" /* empty line */);
/* test RT plan */
DRTPlanIOD *rtplan = new DRTPlanIOD();
if (rtplan != NULL)
{
status = fileformat.loadFile(RTPLAN_FILENAME);
if (status.good())
status = rtplan->read(*fileformat.getDataset());
DCMRT_INFO("Read RT Plan: " << status.text());
if (status.good())
{
fileformat.clear();
status = rtplan->write(*fileformat.getDataset());
DCMRT_INFO("Write RT Plan: " << status.text());
if (status.good())
status = fileformat.saveFile("rtplan.dcm");
DCMRT_INFO("Save RT Plan: " << status.text());
#ifdef TEST_ASSIGNMENT_OPERATOR
DRTPlanIOD *rtplan2 = new DRTPlanIOD();
if (rtplan2 != NULL)
{
*rtplan2 = *rtplan;
fileformat.clear();
status = rtplan2->write(*fileformat.getDataset());
DCMRT_INFO("Write copy of RT Plan: " << status.text());
if (status.good())
status = fileformat.saveFile("rtplan2.dcm");
DCMRT_INFO("Save copy of RT Plan: " << status.text());
delete rtplan2;
rtplan2= NULL;
} else
DCMRT_ERROR("Cannot allocate memory for copy of RT Plan");
#endif
} else
DCMRT_INFO("Skipping further tests for RT Plan");
delete rtplan;
rtplan = NULL;
} else
DCMRT_ERROR("Cannot allocate memory for RT Plan");
DCMRT_INFO("" /* empty line */);
/* test RT structure set */
DRTStructureSetIOD *rtstruct = new DRTStructureSetIOD();
if (rtstruct != NULL)
{
status = fileformat.loadFile(RTSTRUCT_FILENAME);
if (status.good())
status = rtstruct->read(*fileformat.getDataset());
DCMRT_INFO("Read RT Structure Set: " << status.text());
if (status.good())
{
fileformat.clear();
status = rtstruct->write(*fileformat.getDataset());
DCMRT_INFO("Write RT Structure Set: " << status.text());
if (status.good())
status = fileformat.saveFile("rtstruct.dcm");
DCMRT_INFO("Save RT Structure Set: " << status.text());
} else
DCMRT_ERROR("Skipping further tests for RT Structure Set");
delete rtstruct;
rtstruct = NULL;
} else
DCMRT_ERROR("Cannot allocate memory for RT Structure Set");
DCMRT_INFO("" /* empty line */);
/* test RT treatment summary record */
DRTTreatmentSummaryRecordIOD *rttreat = new DRTTreatmentSummaryRecordIOD();
if (rttreat != NULL)
{
status = fileformat.loadFile(RTTREAT_FILENAME);
if (status.good())
status = rttreat->read(*fileformat.getDataset());
DCMRT_INFO("Read RT Treatment Summary Record: " << status.text());
if (status.good())
{
fileformat.clear();
rttreat->getRequestAttributesSequence().clear();
DCMRT_INFO("Fix invalid sequence elements: " << status.text());
status = rttreat->write(*fileformat.getDataset());
DCMRT_INFO("Write RT Treatment Summary Record: " << status.text());
if (status.good())
status = fileformat.saveFile("rttreat.dcm");
DCMRT_INFO("Save RT Treatment Summary Record: " << status.text());
} else
DCMRT_ERROR("Skipping further tests for RT Treatment Summary Record");
delete rttreat;
rttreat = NULL;
} else
DCMRT_ERROR("Cannot allocate memory for RT Treatment Summary Record");
DCMRT_INFO("" /* empty line */);
return 0;
}
|
