Program: gdcm
Module: $RCSfile: vtkGdcmWriter.cxx,v $
Language: C++
- Date: $Date: 2004/12/10 08:34:08 $
- Version: $Revision: 1.7 $
+ Date: $Date: 2005/01/18 18:24:07 $
+ Version: $Revision: 1.10 $
Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
l'Image). All rights reserved. See Doc/License.txt or
#include <vtkPointData.h>
#include <vtkLookupTable.h>
-vtkCxxRevisionMacro(vtkGdcmWriter, "$Revision: 1.7 $");
+vtkCxxRevisionMacro(vtkGdcmWriter, "$Revision: 1.10 $");
vtkStandardNewMacro(vtkGdcmWriter);
//-----------------------------------------------------------------------------
{
this->LookupTable = NULL;
this->FileDimensionality = 3;
+ this->WriteType = VTK_GDCM_WRITE_TYPE_EXPLICIT_VR;
}
vtkGdcmWriter::~vtkGdcmWriter()
void vtkGdcmWriter::PrintSelf(ostream& os, vtkIndent indent)
{
this->Superclass::PrintSelf(os,indent);
+
+ os << indent << "Write type : " << this->GetWriteTypeAsString();
}
//-----------------------------------------------------------------------------
// Public
+const char *vtkGdcmWriter::GetWriteTypeAsString()
+{
+ switch(WriteType)
+ {
+ case VTK_GDCM_WRITE_TYPE_EXPLICIT_VR :
+ return "Explicit VR";
+ case VTK_GDCM_WRITE_TYPE_IMPLICIT_VR :
+ return "Implicit VR";
+ case VTK_GDCM_WRITE_TYPE_ACR :
+ return "ACR";
+ case VTK_GDCM_WRITE_TYPE_ACR_LIBIDO :
+ return "ACR Libido";
+ default :
+ return "Unknow type";
+ }
+}
//-----------------------------------------------------------------------------
// Protected
size_t planeSize = dim[1] * lineSize;
size_t size = dim[2] * planeSize;
- *data = new unsigned char[size];
-
- image->GetIncrements(inc);
- unsigned char *src = (unsigned char *)image->GetScalarPointerForExtent(extent);
- unsigned char *dst = *data + planeSize - lineSize;
- for (int plane = extent[4]; plane <= extent[5]; plane++)
+ if( size>0 )
{
- for (int line = extent[2]; line <= extent[3]; line++)
+ *data = new unsigned char[size];
+
+ image->GetIncrements(inc);
+ unsigned char *src = (unsigned char *)image->GetScalarPointerForExtent(extent);
+ unsigned char *dst = *data + planeSize - lineSize;
+ for (int plane = extent[4]; plane <= extent[5]; plane++)
{
- // Copy one line at proper destination:
- memcpy((void*)dst, (void*)src, lineSize);
+ for (int line = extent[2]; line <= extent[3]; line++)
+ {
+ // Copy one line at proper destination:
+ memcpy((void*)dst, (void*)src, lineSize);
- src += inc[1]*image->GetScalarSize();
- dst -= lineSize;
+ src += inc[1] * image->GetScalarSize();
+ dst -= lineSize;
+ }
+ dst += 2 * planeSize;
}
- dst += 2 * planeSize;
+ }
+ else
+ {
+ *data = NULL;
}
return size;
str.str("");
str << dim[0];
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0011); // Columns
+ file->ReplaceOrCreate(str.str(),0x0028,0x0011); // Columns
str.str("");
str << dim[1];
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0010); // Rows
+ file->ReplaceOrCreate(str.str(),0x0028,0x0010); // Rows
if(dim[2]>1)
{
str.str("");
str << dim[2];
- //file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0012); // Planes
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0008); // Number of Frames
+ //file->ReplaceOrCreate(str.str(),0x0028,0x0012); // Planes
+ file->ReplaceOrCreate(str.str(),0x0028,0x0008); // Number of Frames
}
// Pixel type
str.str("");
str << image->GetScalarSize()*8;
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0100); // Bits Allocated
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0101); // Bits Stored
+ file->ReplaceOrCreate(str.str(),0x0028,0x0100); // Bits Allocated
+ file->ReplaceOrCreate(str.str(),0x0028,0x0101); // Bits Stored
str.str("");
str << image->GetScalarSize()*8-1;
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0102); // High Bit
+ file->ReplaceOrCreate(str.str(),0x0028,0x0102); // High Bit
// Pixel Repr
// FIXME : what do we do when the ScalarType is
{
str << "1"; // Signed
}
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0103); // Pixel Representation
+ file->ReplaceOrCreate(str.str(),0x0028,0x0103); // Pixel Representation
// Samples per pixel
str.str("");
str << image->GetNumberOfScalarComponents();
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0002); // Samples per Pixel
+ file->ReplaceOrCreate(str.str(),0x0028,0x0002); // Samples per Pixel
// Spacing
double *sp = image->GetSpacing();
str.str("");
str << sp[0] << "\\" << sp[1];
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0030); // Pixel Spacing
+ file->ReplaceOrCreate(str.str(),0x0028,0x0030); // Pixel Spacing
str.str("");
str << sp[2];
- file->ReplaceOrCreateByNumber(str.str(),0x0018,0x0088); // Spacing Between Slices
+ file->ReplaceOrCreate(str.str(),0x0018,0x0088); // Spacing Between Slices
// Origin
double *org = image->GetOrigin();
str.str("");
str << org[0] << "\\" << org[1] << "\\" << org[2];
- file->ReplaceOrCreateByNumber(str.str(),0x0020,0x0032); // Image Position Patient
+ file->ReplaceOrCreate(str.str(),0x0020,0x0032); // Image Position Patient
// Window / Level
double *rng=image->GetScalarRange();
str.str("");
str << rng[1]-rng[0];
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x1051); // Window Width
+ file->ReplaceOrCreate(str.str(),0x0028,0x1051); // Window Width
str.str("");
str << (rng[1]+rng[0])/2.0;
- file->ReplaceOrCreateByNumber(str.str(),0x0028,0x1050); // Window Center
+ file->ReplaceOrCreate(str.str(),0x0028,0x1050); // Window Center
// Pixels
unsigned char *data;
size_t size = ReverseData(image,&data);
- file->SetImageData(data,size);
+ file->SetUserData(data,size);
}
/**
SetImageInformation(dcmFile,image);
// Write the image
+ switch(this->WriteType)
+ {
+ case VTK_GDCM_WRITE_TYPE_EXPLICIT_VR :
+ dcmFile->SetWriteTypeToDcmExplVR();
+ break;
+ case VTK_GDCM_WRITE_TYPE_IMPLICIT_VR :
+ dcmFile->SetWriteTypeToDcmImplVR();
+ break;
+ case VTK_GDCM_WRITE_TYPE_ACR :
+ dcmFile->SetWriteTypeToAcr();
+ break;
+ case VTK_GDCM_WRITE_TYPE_ACR_LIBIDO :
+ dcmFile->SetWriteTypeToAcrLibido();
+ break;
+ default :
+ dcmFile->SetWriteTypeToDcmExplVR();
+ }
+
if(!dcmFile->Write(fileName))
{
vtkErrorMacro( << "File " << this->FileName << "couldn't be written by "
<< " the gdcm library");
- std::cerr << "not written \n";
}
+ // Clean up
+ if( dcmFile->GetUserData() && dcmFile->GetUserDataSize()>0 )
+ {
+ delete[] dcmFile->GetUserData();
+ }
delete dcmFile;
}