Program: gdcm
Module: $RCSfile: vtkGdcmWriter.cxx,v $
Language: C++
- Date: $Date: 2005/03/01 16:12:20 $
- Version: $Revision: 1.18 $
+ Date: $Date: 2005/05/11 14:40:58 $
+ Version: $Revision: 1.22 $
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.18 $");
+#ifndef vtkFloatingPointType
+#define vtkFloatingPointType float
+#endif
+
+vtkCxxRevisionMacro(vtkGdcmWriter, "$Revision: 1.22 $");
vtkStandardNewMacro(vtkGdcmWriter);
//-----------------------------------------------------------------------------
this->LookupTable = NULL;
this->FileDimensionality = 3;
this->WriteType = VTK_GDCM_WRITE_TYPE_EXPLICIT_VR;
-
- UIDPrefix = "";
- StudyInstanceUID = "";
- SeriesInstanceUID = "";
- FrameOfReferenceInstanceUID = "";
}
vtkGdcmWriter::~vtkGdcmWriter()
//-----------------------------------------------------------------------------
// Print
-void vtkGdcmWriter::PrintSelf(ostream& os, vtkIndent indent)
+void vtkGdcmWriter::PrintSelf(ostream &os, vtkIndent indent)
{
- this->Superclass::PrintSelf(os,indent);
+ this->Superclass::PrintSelf(os, indent);
os << indent << "Write type : " << this->GetWriteTypeAsString();
}
}
}
-void vtkGdcmWriter::SetUIDPrefix(const char *prefix)
-{
- UIDPrefix = prefix;
-}
-
-const char *vtkGdcmWriter::GetUIDPrefix()
-{
- return UIDPrefix.c_str();
-}
-
-void vtkGdcmWriter::NewStudyInstanceUID()
-{
- StudyInstanceUID = "";
-}
-
-void vtkGdcmWriter::NewSeriesInstanceUID()
-{
- SeriesInstanceUID = "";
-}
-
-void vtkGdcmWriter::NewFrameOfReferenceInstanceUID()
-{
- FrameOfReferenceInstanceUID = "";
-}
-
//-----------------------------------------------------------------------------
// Protected
/**
* Copy the image and reverse the Y axis
*/
-// The output datas must be deleted by the user of the method !!!
+// The output data must be deleted by the user of the method !!!
size_t ReverseData(vtkImageData *image,unsigned char **data)
{
int inc[3];
}
/**
- * Set the datas informations in the file
+ * Set the data informations in the file
*/
-void SetImageInformation(gdcm::FileHelper *file,vtkImageData *image)
+void SetImageInformation(gdcm::FileHelper *file, vtkImageData *image)
{
std::ostringstream str;
extent[3]-extent[2]+1,
extent[5]-extent[4]+1};
- str.seekp(0);
+ str.str("");
str << dim[0];
file->InsertValEntry(str.str(),0x0028,0x0011); // Columns
- str.seekp(0);
+ str.str("");
str << dim[1];
file->InsertValEntry(str.str(),0x0028,0x0010); // Rows
if(dim[2]>1)
{
- str.seekp(0);
+ str.str("");
str << dim[2];
//file->Insert(str.str(),0x0028,0x0012); // Planes
file->InsertValEntry(str.str(),0x0028,0x0008); // Number of Frames
}
// Pixel type
- str.seekp(0);
+ str.str("");
str << image->GetScalarSize()*8;
file->InsertValEntry(str.str(),0x0028,0x0100); // Bits Allocated
file->InsertValEntry(str.str(),0x0028,0x0101); // Bits Stored
- str.seekp(0);
+ str.str("");
str << image->GetScalarSize()*8-1;
file->InsertValEntry(str.str(),0x0028,0x0102); // High Bit
// Pixel Repr
// FIXME : what do we do when the ScalarType is
// VTK_UNSIGNED_INT or VTK_UNSIGNED_LONG
- str.seekp(0);
+ str.str("");
if( image->GetScalarType() == VTK_UNSIGNED_CHAR ||
image->GetScalarType() == VTK_UNSIGNED_SHORT ||
image->GetScalarType() == VTK_UNSIGNED_INT ||
file->InsertValEntry(str.str(),0x0028,0x0103); // Pixel Representation
// Samples per pixel
- str.seekp(0);
+ str.str("");
str << image->GetNumberOfScalarComponents();
file->InsertValEntry(str.str(),0x0028,0x0002); // Samples per Pixel
+ /// \todo : Spacing Between Slices is meaningfull ONLY for CT an MR modality
+ /// We should perform some checkings before forcing the Entry creation
+
// Spacing
- double *sp = image->GetSpacing();
+ vtkFloatingPointType *sp = image->GetSpacing();
- str.seekp(0);
+ str.str("");
// We are about to enter floating point value. By default ostringstream are smart and don't do fixed point
// thus forcing to fixed point value
str.setf( std::ios::fixed );
- str << sp[0] << "\\" << sp[1];
+ str << sp[1] << "\\" << sp[0];
file->InsertValEntry(str.str(),0x0028,0x0030); // Pixel Spacing
- str.seekp(0);
+ str.str("");
str << sp[2];
file->InsertValEntry(str.str(),0x0018,0x0088); // Spacing Between Slices
// Origin
- double *org = image->GetOrigin();
+ vtkFloatingPointType *org = image->GetOrigin();
- str.seekp(0);
+ /// \todo : Image Position Patient is meaningfull ONLY for CT an MR modality
+ /// We should perform some checkings before forcing the Entry creation
+
+ str.str("");
str << org[0] << "\\" << org[1] << "\\" << org[2];
file->InsertValEntry(str.str(),0x0020,0x0032); // Image Position Patient
str.unsetf( std::ios::fixed ); //done with floating point values
// Window / Level
- double *rng=image->GetScalarRange();
+ vtkFloatingPointType *rng = image->GetScalarRange();
- str.seekp(0);
+ str.str("");
str << rng[1]-rng[0];
file->InsertValEntry(str.str(),0x0028,0x1051); // Window Width
- str.seekp(0);
+ str.str("");
str << (rng[1]+rng[0])/2.0;
file->InsertValEntry(str.str(),0x0028,0x1050); // Window Center
* Write of the files
* The call to this method is recursive if there is some files to write
*/
-void vtkGdcmWriter::RecursiveWrite(int axis, vtkImageData *image, ofstream *file)
+void vtkGdcmWriter::RecursiveWrite(int axis, vtkImageData *image,
+ ofstream *file)
{
if(file)
{
- vtkErrorMacro( << "File musn't be opened");
+ vtkErrorMacro( << "File must not be open");
return;
}
if( image->GetScalarType() == VTK_FLOAT ||
image->GetScalarType() == VTK_DOUBLE )
{
- vtkErrorMacro(<< "Bad input type. Scalar type musn't be of type "
+ vtkErrorMacro(<< "Bad input type. Scalar type must not be of type "
<< "VTK_FLOAT or VTKDOUBLE (found:"
<< image->GetScalarTypeAsString());
return;
}
- RecursiveWrite(axis,image,image,file);
+ RecursiveWrite(axis,image, image, file);
//WriteDcmFile(this->FileName,image);
}
// if the file is already open then just write to it
if( file )
{
- vtkErrorMacro( << "File musn't be opened");
+ vtkErrorMacro( << "File musn't be open");
return;
}
// determine the name
if (this->FileName)
{
- sprintf(this->InternalFileName,"%s",this->FileName);
+ sprintf(this->InternalFileName, "%s", this->FileName);
}
else
{
{
sprintf(this->InternalFileName, this->FilePattern,this->FileNumber);
}
+// Remove this code in case user is using VTK 4.2...
+#if !(VTK_MAJOR_VERSION == 4 && VTK_MINOR_VERSION == 2)
if (this->FileNumber < this->MinimumFileNumber)
{
this->MinimumFileNumber = this->FileNumber;
{
this->MaximumFileNumber = this->FileNumber;
}
+#endif
}
// Write the file
return;
}
- // if the current region is too high a dimension forthe file
+ // if the current region is too high a dimension for the file
// the we will split the current axis
cache->GetAxisUpdateExtent(axis, min, max);
cache->SetAxisUpdateExtent(axis, min, max);
}
-void vtkGdcmWriter::WriteDcmFile(char *fileName,vtkImageData *image)
+void vtkGdcmWriter::WriteDcmFile(char *fileName, vtkImageData *image)
{
// From here, the write of the file begins
gdcm::FileHelper *dcmFile = new gdcm::FileHelper();
- // Set the image UID
- if( StudyInstanceUID.empty() )
- StudyInstanceUID = gdcm::Util::CreateUniqueUID( UIDPrefix );
- if( SeriesInstanceUID.empty() )
- SeriesInstanceUID = gdcm::Util::CreateUniqueUID( UIDPrefix );
- if( FrameOfReferenceInstanceUID.empty() )
- FrameOfReferenceInstanceUID = gdcm::Util::CreateUniqueUID( UIDPrefix );
- std::string uid = gdcm::Util::CreateUniqueUID( UIDPrefix );
-
- dcmFile->InsertValEntry(uid,0x0008,0x0018); //[SOP Instance UID]
- dcmFile->InsertValEntry(uid,0x0002,0x0003); //[Media Stored SOP Instance UID]
- dcmFile->InsertValEntry(StudyInstanceUID,0x0020,0x000d); //[Study Instance UID]
- dcmFile->InsertValEntry(SeriesInstanceUID,0x0020,0x000e); //[Series Instance UID]
- dcmFile->InsertValEntry(FrameOfReferenceInstanceUID,0x0020, 0x0052); //[Frame of Reference UID]
-
// Set the image informations
- SetImageInformation(dcmFile,image);
+ SetImageInformation(dcmFile, image);
// Write the image
switch(this->WriteType)
if(!dcmFile->Write(fileName))
{
- vtkErrorMacro( << "File " << this->FileName << "couldn't be written by "
+ vtkErrorMacro( << "File " << this->FileName << "cannot be written by "
<< " the gdcm library");
}