GDCM 2.2.1

- Main Page
- Related Pages
- Namespaces
- Classes
- Files
- Examples

PatchFile.cxx

This is a C++ example on how to use gdcm::Attribute

/*=========================================================================
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;
}