X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=vtk%2FvtkGdcmWriter.cxx;h=333be8f232ba5e4be1a472528bd1319ae05c8cf3;hb=4297b237d3010e089a89e0242d4151168553778f;hp=e4ba6795cbf9b980e326e1844f9b664a04febb30;hpb=446ab2af0e3045cd559b086ced08a8cc5a61d2e2;p=gdcm.git diff --git a/vtk/vtkGdcmWriter.cxx b/vtk/vtkGdcmWriter.cxx index e4ba6795..333be8f2 100644 --- a/vtk/vtkGdcmWriter.cxx +++ b/vtk/vtkGdcmWriter.cxx @@ -3,8 +3,8 @@ Program: gdcm Module: $RCSfile: vtkGdcmWriter.cxx,v $ Language: C++ - Date: $Date: 2004/12/09 11:31:52 $ - Version: $Revision: 1.5 $ + Date: $Date: 2005/10/25 14:52:37 $ + Version: $Revision: 1.26 $ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de l'Image). All rights reserved. See Doc/License.txt or @@ -16,9 +16,10 @@ =========================================================================*/ -#include "gdcmHeader.h" #include "gdcmFile.h" +#include "gdcmFileHelper.h" #include "gdcmDebug.h" +#include "gdcmUtil.h" #include "vtkGdcmWriter.h" #include @@ -26,14 +27,20 @@ #include #include -vtkCxxRevisionMacro(vtkGdcmWriter, "$Revision: 1.5 $"); -vtkStandardNewMacro(vtkGdcmWriter); +#ifndef vtkFloatingPointType +#define vtkFloatingPointType float +#endif + +vtkCxxRevisionMacro(vtkGdcmWriter, "$Revision: 1.26 $") +vtkStandardNewMacro(vtkGdcmWriter) //----------------------------------------------------------------------------- // Constructor / Destructor vtkGdcmWriter::vtkGdcmWriter() { this->LookupTable = NULL; + this->FileDimensionality = 3; + this->WriteType = VTK_GDCM_WRITE_TYPE_EXPLICIT_VR; } vtkGdcmWriter::~vtkGdcmWriter() @@ -42,90 +49,125 @@ 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(); } //----------------------------------------------------------------------------- // 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 /** * Copy the image and reverse the Y axis */ -// The output datas must be deleted by the user of the method !!! -size_t ReverseData(vtkImageData *img,unsigned char **data) +// The output data must be deleted by the user of the method !!! +size_t ReverseData(vtkImageData *image,unsigned char **data) { - int *dim = img->GetDimensions(); - size_t lineSize = dim[0] * img->GetScalarSize() - * img->GetNumberOfScalarComponents(); + int inc[3]; + int *extent = image->GetUpdateExtent(); + int dim[3] = {extent[1]-extent[0]+1, + extent[3]-extent[2]+1, + extent[5]-extent[4]+1}; + + size_t lineSize = dim[0] * image->GetScalarSize() + * image->GetNumberOfScalarComponents(); size_t planeSize = dim[1] * lineSize; size_t size = dim[2] * planeSize; - *data = new unsigned char[size]; - - unsigned char *src = (unsigned char *)img->GetScalarPointer(); - unsigned char *dst = *data + planeSize - lineSize; - for (int plane = 0; plane < dim[2]; plane++) + if( size>0 ) { - for (int line = 0; line < dim[1]; line++) - { - // Copy one line at proper destination: - memcpy((void*)dst, (void*)src, lineSize); + *data = new unsigned char[size]; - src += lineSize; - dst -= lineSize; + 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++) + { + 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; + } + dst += 2 * planeSize; } - dst += 2 * planeSize; + } + else + { + *data = NULL; } return size; } /** - * Set the datas informations in the file + * Set the data informations in the file */ -void SetImageInformation(gdcm::File *file,vtkImageData *image) +void SetImageInformation(gdcm::FileHelper *file, vtkImageData *image) { std::ostringstream str; // Image size - int *dim = image->GetDimensions(); + int *extent = image->GetUpdateExtent(); + int dim[3] = {extent[1]-extent[0]+1, + extent[3]-extent[2]+1, + extent[5]-extent[4]+1}; str.str(""); str << dim[0]; - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0011); + file->InsertEntryString(str.str(),0x0028,0x0011); // Columns str.str(""); str << dim[1]; - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0010); + file->InsertEntryString(str.str(),0x0028,0x0010); // Rows if(dim[2]>1) { str.str(""); str << dim[2]; - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0012); + //file->Insert(str.str(),0x0028,0x0012); // Planes + file->InsertEntryString(str.str(),0x0028,0x0008); // Number of Frames } // Pixel type str.str(""); str << image->GetScalarSize()*8; - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0100); - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0101); + file->InsertEntryString(str.str(),0x0028,0x0100); // Bits Allocated + file->InsertEntryString(str.str(),0x0028,0x0101); // Bits Stored str.str(""); str << image->GetScalarSize()*8-1; - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0102); + file->InsertEntryString(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.str(""); - if( image->GetScalarType() == VTK_UNSIGNED_CHAR || + if( image->GetScalarType() == VTK_UNSIGNED_CHAR || image->GetScalarType() == VTK_UNSIGNED_SHORT || - image->GetScalarType() == VTK_UNSIGNED_INT || + image->GetScalarType() == VTK_UNSIGNED_INT || image->GetScalarType() == VTK_UNSIGNED_LONG ) { str << "0"; // Unsigned @@ -134,97 +176,198 @@ void SetImageInformation(gdcm::File *file,vtkImageData *image) { str << "1"; // Signed } - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0103); + file->InsertEntryString(str.str(),0x0028,0x0103); // Pixel Representation // Samples per pixel str.str(""); str << image->GetNumberOfScalarComponents(); - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0002); + file->InsertEntryString(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.str(""); - str << sp[0] << "\\" << sp[1]; - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x0030); + // 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[1] << "\\" << sp[0]; + file->InsertEntryString(str.str(),0x0028,0x0030); // Pixel Spacing str.str(""); str << sp[2]; - file->ReplaceOrCreateByNumber(str.str(),0x0018,0x0088); + file->InsertEntryString(str.str(),0x0018,0x0088); // Spacing Between Slices // Origin - double *org = image->GetOrigin(); + vtkFloatingPointType *org = image->GetOrigin(); + + /// \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->ReplaceOrCreateByNumber(str.str(),0x0020,0x0032); + file->InsertEntryString(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.str(""); str << rng[1]-rng[0]; - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x1051); + file->InsertEntryString(str.str(),0x0028,0x1051); // Window Width str.str(""); str << (rng[1]+rng[0])/2.0; - file->ReplaceOrCreateByNumber(str.str(),0x0028,0x1050); + file->InsertEntryString(str.str(),0x0028,0x1050); // Window Center // Pixels unsigned char *data; size_t size = ReverseData(image,&data); - file->SetImageData(data,size); + file->SetUserData(data,size); } /** * 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) { - (void)axis; // To avoid warning 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; } - WriteDcmFile(this->FileName,image); + RecursiveWrite(axis,image, image, file); + //WriteDcmFile(this->FileName,image); } -void vtkGdcmWriter::RecursiveWrite(int axis, vtkImageData *image, - vtkImageData *cache, ofstream *file) +void vtkGdcmWriter::RecursiveWrite(int axis, vtkImageData *cache, + vtkImageData *image, ofstream *file) { - (void)axis; // To avoid warning - (void)image; // To avoid warning - (void)cache; // To avoid warning - (void)file; // To avoid warning + int idx, min, max; + + // if the file is already open then just write to it + if( file ) + { + vtkErrorMacro( << "File musn't be open"); + return; + } + + // if we need to open another slice, do it + if( (axis + 1) == this->FileDimensionality ) + { + // determine the name + if (this->FileName) + { + sprintf(this->InternalFileName, "%s", this->FileName); + } + else + { + if (this->FilePrefix) + { + sprintf(this->InternalFileName, this->FilePattern, + this->FilePrefix, this->FileNumber); + } + 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; + } + else if (this->FileNumber > this->MaximumFileNumber) + { + this->MaximumFileNumber = this->FileNumber; + } +#endif + } + + // Write the file + WriteDcmFile(this->InternalFileName,image); + ++this->FileNumber; + return; + } + + // if the current region is too high a dimension for the file + // the we will split the current axis + cache->GetAxisUpdateExtent(axis, min, max); + + // if it is the y axis then flip by default + if (axis == 1 && !this->FileLowerLeft) + { + for(idx = max; idx >= min; idx--) + { + cache->SetAxisUpdateExtent(axis, idx, idx); + this->RecursiveWrite(axis - 1, cache, image, file); + } + } + else + { + for(idx = min; idx <= max; idx++) + { + cache->SetAxisUpdateExtent(axis, idx, idx); + this->RecursiveWrite(axis - 1, cache, image, file); + } + } + + // restore original extent + 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::File *dcmFile = new gdcm::File(); + gdcm::FileHelper *dcmFile = gdcm::FileHelper::New(); // Set the image informations - SetImageInformation(dcmFile,image); + 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 " + vtkErrorMacro( << "File " << this->FileName << "cannot be written by " << " the gdcm library"); - std::cerr << "not written \n"; } - delete dcmFile; + // Clean up + if( dcmFile->GetUserData() && dcmFile->GetUserDataSize()>0 ) + { + delete[] dcmFile->GetUserData(); + } + dcmFile->Delete(); } //-----------------------------------------------------------------------------