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
|
/*=========================================================================
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 "vtkGDCMImageReader.h"
#include "vtkGDCMImageWriter.h"
#include "vtkImageData.h"
#include "vtkStringArray.h"
#include "vtkVersion.h"
#include "gdcmTesting.h"
#include "gdcmFilenameGenerator.h"
int main(int argc, char *argv[])
{
std::string filename;
if( argc <= 1 )
{
const char *directory = gdcm::Testing::GetDataRoot();
if(!directory) return 1;
std::string file = std::string(directory) + "/US-PAL-8-10x-echo.dcm";
filename = file;
}
else
{
filename = argv[1];
}
std::cout << "file: " << filename << std::endl;
vtkGDCMImageReader *reader = vtkGDCMImageReader::New();
reader->SetFileName( filename.c_str() );
reader->Update();
//reader->GetOutput()->Print( std::cout );
int dims[3];
reader->GetOutput()->GetDimensions( dims );
std::ostringstream os;
os << "singleframe";
os << "%04d.dcm";
gdcm::FilenameGenerator fg;
fg.SetPattern( os.str().c_str() );
unsigned int nfiles = dims[2];
fg.SetNumberOfFilenames( nfiles );
bool b = fg.Generate();
if( !b )
{
std::cerr << "FilenameGenerator::Generate() failed" << std::endl;
return 1;
}
if( !fg.GetNumberOfFilenames() )
{
std::cerr << "FilenameGenerator::Generate() failed somehow..." << std::endl;
return 1;
}
// By default write them as Secondary Capture (for portability)
vtkGDCMImageWriter *writer = vtkGDCMImageWriter::New();
vtkStringArray *filenames = vtkStringArray::New();
for(unsigned int i = 0; i < fg.GetNumberOfFilenames(); ++i)
{
filenames->InsertNextValue( fg.GetFilename(i) );
}
assert( filenames->GetNumberOfValues() == (int)fg.GetNumberOfFilenames() );
writer->SetFileNames( filenames );
filenames->Delete();
writer->SetFileDimensionality( 2 );
#if (VTK_MAJOR_VERSION >= 6)
writer->SetInputConnection( reader->GetOutputPort() );
#else
writer->SetInput( reader->GetOutput() );
#endif
writer->SetImageFormat( reader->GetImageFormat() );
writer->Write();
reader->Delete();
writer->Delete();
return 0;
}
|
