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
|
/*=========================================================================
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.
=========================================================================*/
/*
* The image was a broken file where the Pixel Data element was 8 times too big
* Apparently multiplying the BitsAllocated to 4 and multiplying the number of
* frames by 2 would solve the problem
*
* This C++ code can be used to patch the header.
*/
#include "gdcmReader.h"
#include "gdcmImageReader.h"
#include "gdcmWriter.h"
#include "gdcmDataSet.h"
#include "gdcmAttribute.h"
int main(int argc, char *argv[])
{
if( argc < 3 )
{
return 1;
}
const char *f = argv[1];
const char *out = argv[2];
gdcm::Reader r;
r.SetFileName( f );
if( !r.Read() )
{
return 1;
}
gdcm::File &file = r.GetFile();
gdcm::DataSet& ds = file.GetDataSet();
// (0028,0100) US 16 # 2, 1 BitsAllocated
// (0028,0101) US 16 # 2, 1 BitsStored
// (0028,0102) US 15 # 2, 1 HighBit
//
{
gdcm::Attribute<0x28,0x100> at;
at.SetFromDataElement( ds.GetDataElement( at.GetTag() ) );
if( at.GetValue() != 8 )
{
return 1;
}
at.SetValue( 32 );
ds.Replace( at.GetAsDataElement() );
}
{
gdcm::Attribute<0x28,0x101> at;
at.SetFromDataElement( ds.GetDataElement( at.GetTag() ) );
if( at.GetValue() != 8 )
{
return 1;
}
at.SetValue( 32 );
ds.Replace( at.GetAsDataElement() );
}
{
gdcm::Attribute<0x28,0x102> at;
at.SetFromDataElement( ds.GetDataElement( at.GetTag() ) );
if( at.GetValue() != 7 )
{
return 1;
}
at.SetValue( 31 );
ds.Replace( at.GetAsDataElement() );
}
// (0028,0008) IS [56] # 2, 1 NumberOfFrames
{
gdcm::Attribute<0x28,0x8> at;
at.SetFromDataElement( ds.GetDataElement( at.GetTag() ) );
at.SetValue( at.GetValue() * 2 );
ds.Replace( at.GetAsDataElement() );
}
gdcm::Writer w;
w.SetFile( file );
w.SetCheckFileMetaInformation( false );
w.SetFileName( out );
if( !w.Write() )
{
return 1;
}
// Now let's see if we can read it as an image:
gdcm::ImageReader ir;
ir.SetFileName( out );
if(!ir.Read())
{
return 1;
}
gdcm::Image &image = ir.GetImage();
unsigned long len = image.GetBufferLength();
const gdcm::ByteValue *bv = ir.GetFile().GetDataSet().GetDataElement( gdcm::Tag(0x7fe0,0x0010) ).GetByteValue();
if( !bv || len != bv->GetLength() )
{
return 1;
}
std::cout << bv->GetLength() << " " << len << std::endl;
std::cout << "Sucess to rewrite image !" << std::endl;
image.Print( std::cout );
return 0;
}
|
