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: DICOM for VTK
Copyright (c) 2012-2024 David Gobbi
All rights reserved.
See Copyright.txt or http://dgobbi.github.io/bsd3.txt 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 "vtkDICOMMetaData.h"
#include "vtkDICOMFileSorter.h"
#include "vtkDICOMReader.h"
#include "vtkDICOMWriter.h"
#include "vtkDICOMMRGenerator.h"
#include "vtkSmartPointer.h"
#include "vtkStringArray.h"
#include <sstream>
#include <string.h>
#include <stdlib.h>
int main(int argc, char *argv[])
{
int rval = 0;
const char *exename = argv[0];
// remove path portion of exename
const char *cp = exename + strlen(exename);
while (cp != exename && cp[-1] != '\\' && cp[-1] != '/') { --cp; }
exename = cp;
vtkSmartPointer<vtkDICOMFileSorter> sorter =
vtkSmartPointer<vtkDICOMFileSorter>::New();
vtkSmartPointer<vtkStringArray> files =
vtkSmartPointer<vtkStringArray>::New();
for (int i = 1; i < argc; i++)
{
files->InsertNextValue(argv[i]);
}
sorter->SetInputFileNames(files);
sorter->Update();
int m = sorter->GetNumberOfStudies();
for (int j = 0; j < m; j++)
{
cout << "Study" << j << ":\n";
int k = sorter->GetFirstSeriesForStudy(j);
int kl = sorter->GetLastSeriesForStudy(j);
for (; k <= kl; k++)
{
cout << " Series " << k << ":\n";
vtkStringArray *a = sorter->GetFileNamesForSeries(k);
vtkSmartPointer<vtkDICOMReader> reader =
vtkSmartPointer<vtkDICOMReader>::New();
reader->SetFileNames(a);
//reader->SetMemoryRowOrderToFileNative();
reader->UpdateInformation();
vtkDICOMMetaData *meta = reader->GetMetaData();
vtkSmartPointer<vtkDICOMMRGenerator> generator =
vtkSmartPointer<vtkDICOMMRGenerator>::New();
vtkSmartPointer<vtkDICOMWriter> writer =
vtkSmartPointer<vtkDICOMWriter>::New();
writer->SetInputConnection(reader->GetOutputPort());
writer->SetGenerator(generator);
writer->SetPatientMatrix(reader->GetPatientMatrix());
writer->SetMetaData(meta);
//writer->SetMemoryRowOrderToFileNative();
writer->SetFilePrefix("/tmp");
writer->SetFilePattern("%s/IM-0001-%04.4d.dcm");
writer->Write();
}
}
return rval;
}
|
