/*========================================================================= Program: wxMaracas Module: $RCSfile: marSimpleDicom.cpp,v $ Language: C++ Date: $Date: 2008/10/31 16:32:56 $ Version: $Revision: 1.1 $ Copyright: (c) 2002, 2003 License: 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. =========================================================================*/ #ifdef _MSC_VER #pragma warning ( disable : 4786 ) #pragma warning ( disable : 4251 ) #endif //_MSC_VER // PS -> #include "gdcm.h" #include "marSimpleDicom.h" #include #include #include #include #include // PS -> #include "vtkGdcmReader.h" #include #include #include #include // PS -> #include #include #include #include #include #include //---------------------------------------------------------------------------- // Callback for the interaction class marProgressObserver : public vtkCommand { public: static marProgressObserver *New() { return new marProgressObserver; } marProgressObserver() { this->FrameGauge = NULL; } virtual void Execute(vtkObject *wdg, unsigned long event, void* calldata) { if ( this->FrameGauge ) { int val = this->FrameGauge->GetValue(); this->FrameGauge->SetValue( val + 1); //wxYield(); //wxSafeYield(); //wxYieldIfNeeded(); //wxWakeUpIdle(); //wxSafeYield( this->FrameGauge, FALSE); } } wxGauge *FrameGauge; }; //---------------------------------------------------------------------------- static char* DicomTagsNames[] = { "ID_File_Name", "ID_SOP_Class_UID", //!0008 0016 UI ID SOP Class UID "ID_SOP_Instance_UID", //!0008 0018 UI ID SOP Instance UID "ID_Study_Date", //!0008 0020 DA ID Study Date "ID_Series_Date", //!0008 0021 DA ID Series Date "ID_Acquisition_Date", //!0008 0022 DA ID Acquisition Date "ID_Image_Date", //!0008 0023 DA ID Image Date "ID_Study_Time", //!0008 0030 TM ID Study Time "ID_Series_Time", //!0008 0031 TM ID Series Time "ID_Acquisition_Time", //!0008 0032 TM ID Acquisition Time "ID_Image_Time", //!0008 0033 TM ID Image Time "ID_Modality", //!0008 0060 CS ID Modality "ID_Manufacturer", //!0008 0070 LO ID Manufacturer "ID_Institution_Name", //!0008 0080 LO ID Institution Name "ID_Study_Description", //!0008 1030 LO ID Study Description "ID_Series_Description", //!0008 103e LO ID Series Description "ID_Admitting_Diagnoses_Description", //!0008 1080 LO ID Admitting Diagnoses Description "ID_Patient_Name", //!0010 0010 PN PAT Patient Name "ID_Patient_ID", //!0010 0020 LO PAT Patient ID "ID_Body_Part_Examined", //!0018 0015 CS ACQ Body Part Examined "ID_Scanning_Sequence", //!0018 0020 CS ACQ Scanning Sequence "ID_Sequence_Variant", //!0018 0021 CS ACQ Sequence Variant "ID_Scan_Options", //!0018 0022 CS ACQ Scan Options "ID_MR_Acquisition_Type", //!0018 0023 CS ACQ MR Acquisition Type "ID_Sequence_Name", //!0018 0024 SH ACQ Sequence Name "ID_Slice_Thickness", //!0018 0050 DS ACQ Slice Thickness "ID_Repetition_Time", //!0018 0080 DS ACQ Repetition Time "ID_Echo_Time", //!0018 0081 DS ACQ Echo Time "ID_Inversion_Time", //!0018 0082 DS ACQ Inversion Time "ID_Number_of_Averages", //!0018 0083 DS ACQ Number of Averages "ID_Imaging_Frequency", //!0018 0084 DS ACQ Imaging Frequency "ID_Imaged_Nucleus", //!0018 0085 SH ACQ Imaged Nucleus "ID_Echo_Number", //!0018 0086 IS ACQ Echo Number "ID_Magnetic_Field_Strength", //!0018 0087 DS ACQ Magnetic Field Strength "ID_Spacing_Between_Slices", //!0018 0088 DS ACQ Spacing Between Slices "ID_Echo_Train_Length", //!0018 0091 IS ACQ Echo Train Length "ID_Percent_Sampling", //!0018 0093 DS ACQ Percent Sampling "ID_Percent_Phase_Field_of_View", //!0018 0094 DS ACQ Percent Phase Field of View "ID_Receiving_Coil", //!0018 1250 SH ACQ Receiving Coil "ID_Patient_Position", //!0018 5100 CS ACQ Patient Position "ID_Study_Instance_UID", //!0020 000d UI REL Study Instance UID "ID_Series_Instance_UID", //!0020 000e UI REL Series Instance UID "ID_Study_ID", //!0020 0010 SH REL Study ID "ID_Series_Number", //!0020 0011 IS REL Series Number "ID_Acquisition_Number", //!0020 0012 IS REL Acquisition Number "ID_Image_Number", //!0020 0013 IS REL Image Number "ID_Patient_Orientation", //!0020 0020 CS REL Patient Orientation "ID_Image_Position", //!0020 0030 RET REL Image Position "ID_Image_Position_Patient", //!0020 0032 DS REL Image Position Patient "ID_Image_Orientation", //!0020 0035 RET REL Image Orientation "ID_Image_Orientation_Patient", //!0020 0037 DS REL Image Orientation (Patient) "ID_Location", //!0020 0050 RET REL Location "ID_Frame_of_Reference_UID", //!0020 0052 UI REL Frame of Reference UID "ID_Slice_Location", //!0020 1041 DS REL Slice Location "ID_Image_Comments", //!0020 4000 LT REL Image Comments "ID_Pixel_Spacing", //!0028 0030 DS IMG Pixel Spacing "ID_Window_Center", //!0028 1050 DS IMG Window Center "ID_Window_Width", //!0028 1051 DS IMG Window Width }; // ------------------------------------------------------------------------- marSimpleDicom::marSimpleDicom( marParameters* p ) : marFilesBase( p ), m_pSlicesPixels(NULL), m_ActualStudy(0), m_ActualSerie(0) // PS -> ,helper(NULL) { // FillDicomInfo( ); } // ------------------------------------------------------------------------- marSimpleDicom::~marSimpleDicom(){ reset( ); } /** Conversion de wxString en string. */ //#define ws2s(as) (string(as.GetData())) /** Conversion de string en wxString. */ #define s2ws(s) (wxString((s).c_str())) /** Conversion de C-string en wxString. */ //#define cs2ws(s) (wxString(s)) // ------------------------------------------------------------------------- // PS -> wxArrayString marGdcmDicom::gdcmGetExamInfo( wxString directoryname ) // PS -> { // PS -> wxArrayString temparray; // PS -> temparray.Alloc( 58 ); // PS -> // PS -> wxDir dir(directoryname); // PS -> wxString filename; // PS -> bool cont = dir.GetFirst( &filename, "*", wxDIR_FILES ); // PS -> filename = directoryname + "/" + filename; // PS -> // PS -> std::string val; // PS -> // PS -> gdcmHeader GdcmHeader( filename.c_str() ); // PS -> temparray.Add( filename ); // "ID_File_Name", // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0016).c_str() );// "ID_SOP_Class_UID", //0008 0016 UI ID SOP Class UID // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0018).c_str() );// "ID_SOP_Instance_UID", //0008 0018 UI ID SOP Instance UID // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0020).c_str() );// "ID_Study_Date", //0008 0020 DA ID Study Date // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0021).c_str() );// "ID_Series_Date", //0008 0021 DA ID Series Date // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0022).c_str() );// "ID_Acquisition_Date", //0008 0022 DA ID Acquisition Date // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0023).c_str() );// "ID_Image_Date", //0008 0023 DA ID Image Date // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0030).c_str() );// "ID_Study_Time", //0008 0030 TM ID Study Time // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0031).c_str() );// "ID_Series_Time", //0008 0031 TM ID Series Time // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0032).c_str() );// "ID_Acquisition_Time", //0008 0032 TM ID Acquisition Time // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0033).c_str() );// "ID_Image_Time", //0008 0033 TM ID Image Time // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0060).c_str() );// "ID_Modality", //0008 0060 CS ID Modality // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0070).c_str() );// "ID_Manufacturer", //0008 0070 LO ID Manufacturer // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x0080).c_str() );// "ID_Institution_Name", //0008 0080 LO ID Institution Name // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x1030).c_str() );// "ID_Study_Description", //0008 1030 LO ID Study Description // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x103e).c_str() );// "ID_Series_Description", //0008 103e LO ID Series Description // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0008,0x1080).c_str() );// "ID_Admitting_Diagnoses_Description", //0008 1080 LO ID Admitting Diagnoses Description // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0010,0x0010).c_str() );// "ID_Patient_Name", //0010 0010 PN PAT Patient Name // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0010,0x0020).c_str() );// "ID_Patient_ID", //0010 0020 LO PAT Patient ID // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0015).c_str() );// "ID_Body_Part_Examined", //0018 0015 CS ACQ Body Part Examined // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0020).c_str() );// "ID_Scanning_Sequence", //0018 0020 CS ACQ Scanning Sequence // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0021).c_str() );// "ID_Sequence_Variant", //0018 0021 CS ACQ Sequence Variant // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0022).c_str() );// "ID_Scan_Options", //0018 0022 CS ACQ Scan Options // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0023).c_str() );// "ID_MR_Acquisition_Type", //0018 0023 CS ACQ MR Acquisition Type // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0024).c_str() );// "ID_Sequence_Name", //0018 0024 SH ACQ Sequence Name // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0050).c_str() );// "ID_Slice_Thickness", //0018 0050 DS ACQ Slice Thickness // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0080).c_str() );// "ID_Repetition_Time", //0018 0080 DS ACQ Repetition Time // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0081).c_str() );// "ID_Echo_Time", //0018 0081 DS ACQ Echo Time // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0082).c_str() );// "ID_Inversion_Time", //0018 0082 DS ACQ Inversion Time // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0083).c_str() );// "ID_Number_of_Averages", //0018 0083 DS ACQ Number of Averages // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0084).c_str() );// "ID_Imaging_Frequency", //0018 0084 DS ACQ Imaging Frequency // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0085).c_str() );// "ID_Imaged_Nucleus", //0018 0085 SH ACQ Imaged Nucleus // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0086).c_str() );// "ID_Echo_Number", //0018 0086 IS ACQ Echo Number // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0087).c_str() );// "ID_Magnetic_Field_Strength", //0018 0087 DS ACQ Magnetic Field Strength // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0088).c_str() );// "ID_Spacing_Between_Slices", //0018 0088 DS ACQ Spacing Between Slices // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0091).c_str() );// "ID_Echo_Train_Length", //0018 0091 IS ACQ Echo Train Length // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0093).c_str() );// "ID_Percent_Sampling", //0018 0093 DS ACQ Percent Sampling // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x0094).c_str() );// "ID_Percent_Phase_Field_of_View", //0018 0094 DS ACQ Percent Phase Field of View // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x1250).c_str() );// "ID_Receiving_Coil", //0018 1250 SH ACQ Receiving Coil // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0018,0x5100).c_str() );// "ID_Patient_Position", //0018 5100 CS ACQ Patient Position // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x000d).c_str() );// "ID_Study_Instance_UID", //0020 000d UI REL Study Instance UID // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x000e).c_str() );// "ID_Series_Instance_UID", //0020 000e UI REL Series Instance UID // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0010).c_str() );// "ID_Study_ID", //0020 0010 SH REL Study ID // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0011).c_str() );// "ID_Series_Number", //0020 0011 IS REL Series Number // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0012).c_str() );// "ID_Acquisition_Number", //0020 0012 IS REL Acquisition Number // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0013).c_str() );// "ID_Image_Number", //0020 0013 IS REL Image Number // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0020).c_str() );// "ID_Patient_Orientation", //0020 0020 CS REL Patient Orientation // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0030).c_str() );// "ID_Image_Position", //0020 0030 RET REL Image Position // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0032).c_str() );// "ID_Image_Position_Patient", //0020 0032 DS REL Image Position Patient // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0035).c_str() );// "ID_Image_Orientation", //0020 0035 RET REL Image Orientation // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0037).c_str() );// "ID_Image_Orientation_Patient", //0020 0037 DS REL Image Orientation (Patient) // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0050).c_str() );// "ID_Location", //0020 0050 RET REL Location // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x0052).c_str() );// "ID_Frame_of_Reference_UID", //0020 0052 UI REL Frame of Reference UID // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x1041).c_str() );// "ID_Slice_Location", //0020 1041 DS REL Slice Location // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0020,0x4000).c_str() );// "ID_Image_Comments", //0020 4000 LT REL Image Comments // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0028,0x0030).c_str() );// "ID_Pixel_Spacing", //0028 0030 DS IMG Pixel Spacing // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0028,0x1050).c_str() );// "ID_Window_Center", //0028 1050 DS IMG Window Center // PS -> temparray.Add( GdcmHeader.GetPubElValByNumber(0x0028,0x1051).c_str() );// "ID_Window_Width", //0028 1051 DS IMG Window Width // PS -> // PS -> return temparray; // PS -> // PS -> } #undef s2ws // ---------------------------------------------------------------------------- // PS -> void marGdcmDicom::FillDicomInfo() // PS -> { // PS -> reset(); // PS -> // PS -> //!\todo FIXME : m_DicomInfo.size() == 0 // PS -> //if( m_DicomInfo.size() == 0 ) // PS -> { // PS -> wxString file,subsubdir, subdir, rootdir = getParameters( )->getStringParam( marParameters::e_dicom_images_directory ); // PS -> wxDir dir(rootdir); // PS -> bool cont = dir.GetFirst( &subdir, "*", wxDIR_DIRS ); // PS -> while ( cont ) // PS -> { // PS -> subdir = rootdir + "/" + subdir; // PS -> // PS -> SerieInfo tempserie; // PS -> wxDir dir2( subdir ); // PS -> bool cont2 = dir2.GetFirst( &subsubdir, "*", wxDIR_DIRS ); // PS -> while ( cont2 ) // PS -> { // PS -> subsubdir = subdir + "/" + subsubdir; // PS -> wxArrayString temptemp(gdcmGetExamInfo(subsubdir)); // PS -> tempserie.push_back(temptemp); // PS -> cont2 = dir2.GetNext(&subsubdir); // PS -> } // PS -> // PS -> m_DicomInfo.push_back( tempserie ); // PS -> cont = dir.GetNext(&subdir); // PS -> } // PS -> } // PS -> } //------------------------------------------------------------------------- void marSimpleDicom::loadActualSerie( wxGauge* gauge) { // PS -> #ifdef DXMM loadVolumeDXMM( true , gauge); // PS -> #else // PS -> loadVolume( true , gauge); // PS -> #endif } // ------------------------------------------------------------------------- // PS -> void marGdcmDicom::loadVolume( bool force , wxGauge *gauge) // PS -> { // PS -> float spacing[3]; // PS -> vtkGdcmReader *gdcmReader = NULL; // PS -> // PS -> if( force || !_volume ) // PS -> { // PS -> freeVolume( ); // PS -> // PS -> wxFileName temp(m_DicomInfo[m_ActualStudy][m_ActualSerie][ID_File_Name]); // PS -> // PS -> helper = new gdcmSerieHeaderHelper(); // PS -> helper->SetDirectory( temp.GetPath(wxPATH_GET_VOLUME | wxPATH_GET_SEPARATOR ).c_str() ); // PS -> helper->OrderGdcmFileList(); // PS -> // PS -> std::list flist = helper->GetGdcmFileList(); // PS -> if( flist.size() > 0) // PS -> { // PS -> gdcmReader = vtkGdcmReader::New(); // PS -> for(std::list::iterator it = flist.begin(); it != flist.end(); it++ ) // PS -> { // PS -> gdcmReader->AddFileName( (*it)->GetFileName().c_str() ); // PS -> } // PS -> } // fi // PS -> } // fi // PS -> // PS -> marProgressObserver *progressEvent = marProgressObserver::New(); // PS -> progressEvent->FrameGauge = gauge; // PS -> //gdcmReader->AddObserver( vtkCommand::ProgressEvent, progressEvent); // PS -> progressEvent->Delete(); // PS -> // PS -> //To simplyfy calculation later I turn Origin to (0, 0, 0) // PS -> //this shouldn't be 'trop grave' // PS -> vtkImageChangeInformation *change = vtkImageChangeInformation::New(); // PS -> change->SetInput( gdcmReader->GetOutput() ); // PS -> change->SetOutputOrigin( 0, 0, 0); // PS -> change->Update(); // PS -> // PS -> float *range = change->GetOutput()->GetScalarRange(); // PS -> // PS -> //Our algorithm are made based on unsigned short pixel/voxel: // PS -> //We need to cast our imagedata // PS -> vtkImageCast *cast = vtkImageCast::New(); // PS -> cast->SetInput( change->GetOutput() ); // PS -> cast->SetOutputScalarTypeToUnsignedShort(); // PS -> cast->Update(); //important // PS -> // PS -> //As I need to shift & scale I don't need anymore vtkImageCast: // PS -> /* vtkImageShiftScale *shift = vtkImageShiftScale::New(); // PS -> shift->SetInput( change->GetOutput() ); // PS -> shift->SetShift( -range[0] ); // PS -> shift->SetScale( VTK_UNSIGNED_SHORT_MAX / (range[1] - range[0]) ); // PS -> shift->SetOutputScalarTypeToUnsignedShort(); // PS -> shift->Update();*/ // PS -> // PS -> cast->GetOutput()->GetSpacing( spacing ); // PS -> getParameters( )->setDoubleParam( marParameters::e_voxel_x_dimension, spacing[0] ); // PS -> getParameters( )->setDoubleParam( marParameters::e_voxel_y_dimension, spacing[1] ); // PS -> getParameters( )->setDoubleParam( marParameters::e_voxel_z_dimension, spacing[2] ); // PS -> // PS -> //FIXME: kVolume doesn't need to be -again- matrix rotated (done in vtkGdcmReader)!! // PS -> _volume = new kVolume( cast->GetOutput() ); // PS -> // PS -> //!\todo thanks // PS -> if(gdcmReader) gdcmReader->Delete(); // PS -> cast->Delete(); // PS -> change->Delete(); // PS -> } void marSimpleDicom::loadVolumeDXMM( bool force , wxGauge *gauge) { // PS -> // PS -> FILE *stream; // PS -> char string[20]; // PS -> wxString prefix, fileParam, fileVol, rootdir; // PS -> rootdir=wxString("c:/temp"); // PS -> wxDir dir(rootdir); // PS -> int dimX, dimY, dimZ; // PS -> float spacing[3]; // PS -> vtkImageReader *volVTK = NULL; // PS -> // PS -> fileParam = rootdir + "/" + "paramDxMM.txt"; // PS -> stream = fopen(fileParam.c_str(),"r"); // PS -> // PS -> wxString errorMsg; // PS -> errorMsg= "Cannot open file "+fileParam; // PS -> wxASSERT_MSG(stream!=NULL,errorMsg); // PS -> // PS -> fscanf(stream,"%s", string); // PS -> dimX = atoi(string); // PS -> fscanf(stream,"%s", string); // PS -> dimY = atoi(string); // PS -> fscanf(stream,"%s", string); // PS -> dimZ = atoi(string); // PS -> fscanf(stream,"%s", string); // PS -> spacing[0] = atof(string); // PS -> fscanf(stream,"%s", string); // PS -> spacing[1] = atof(string); // PS -> fscanf(stream,"%s", string); // PS -> spacing[2] = atof(string); // PS -> fclose(stream); // PS -> // PS -> getParameters( )->setDoubleParam( marParameters::e_voxel_x_dimension, spacing[0] ); // PS -> getParameters( )->setDoubleParam( marParameters::e_voxel_y_dimension, spacing[1] ); // PS -> getParameters( )->setDoubleParam( marParameters::e_voxel_z_dimension, spacing[2] ); // PS -> // PS -> if( force || !_volume ) // PS -> { // PS -> freeVolume( ); // PS -> } // fi // PS -> // PS -> marProgressObserver *progressEvent = marProgressObserver::New(); // PS -> progressEvent->FrameGauge = gauge; // PS -> progressEvent->Delete(); // PS -> // PS -> prefix = rootdir + "/" + "imaDxMM"; // PS -> volVTK = vtkImageReader::New(); // PS -> volVTK->SetDataByteOrderToLittleEndian(); // PS -> volVTK->SetDataScalarTypeToUnsignedShort(); // PS -> volVTK->SetDataExtent(0, dimX-1, 0, dimY-1, 0, dimZ-1); // PS -> volVTK->SetFilePrefix( prefix.c_str() ); // PS -> volVTK->SetFilePattern( "%s%d.raw"); // PS -> volVTK->SetDataSpacing(spacing[0],spacing[1],spacing[2]); // PS -> volVTK->Update(); // PS -> vtkImageResample* ir = vtkImageResample::New( ); // PS -> ir->SetInput(volVTK->GetOutput()); // PS -> ir->SetDimensionality( 3 ); // PS -> ir->SetAxisOutputSpacing( 0, getParameters( )->getVoxelSize( ) ); // PS -> ir->SetAxisOutputSpacing( 1, getParameters( )->getVoxelSize( ) ); // PS -> ir->SetAxisOutputSpacing( 2, getParameters( )->getVoxelSize( ) ); // PS -> ir->InterpolateOn( ); // PS -> ir->Update( ); //important // PS -> // PS -> _volume = new kVolume( ir->GetOutput( ) ); // PS -> // PS -> //!\todo thanks // PS -> if (ir) ir->Delete( ); // PS -> if(volVTK) volVTK->Delete(); if (m_pSlicesPixels!=NULL) { vtkUnsignedShortArray * dataArray=vtkUnsignedShortArray::New(); int iRow,iColumn,iSlice,iPixel; int iiSlice; for (iiSlice=0;iiSliceInsertNextValue(m_pSlicesPixels[iPixel]); } } } vtkImageData*imgData=vtkImageData::New(); imgData->SetDimensions(m_nX,m_nY,m_nSlices); imgData->SetSpacing(m_PixelSpacingRow,m_PixelSpacingColumn,m_PixelSpacingSlice); imgData->SetOrigin(0,0,0); imgData->SetScalarTypeToUnsignedShort(); imgData->GetPointData()->SetScalars(dataArray); // TO DO : remplir le vtkImageResample vtkImageResample* ir = vtkImageResample::New(); ir->SetInput(imgData); ir->SetDimensionality( 3 ); ir->SetAxisOutputSpacing( 0, getParameters()->getVoxelSize() ); ir->SetAxisOutputSpacing( 1, getParameters()->getVoxelSize() ); ir->SetAxisOutputSpacing( 2, getParameters()->getVoxelSize() ); ir->InterpolateOn( ); ir->Update( ); //important // EED // _volume = new kVolume( ir->GetOutput( ) ); SetVolume( new kVolume( ir->GetOutput() ) ); //!\todo thanks if (ir) ir->Delete(); if(imgData) imgData->Delete(); if(dataArray) dataArray->Delete(); getParameters( )->setDoubleParam( marParameters::e_RescaleSlope , m_RescaleSlope ); getParameters( )->setDoubleParam( marParameters::e_RescaleIntercept, m_RescaleIntercept ); } } // ------------------------------------------------------------------------- void marSimpleDicom::copyFrom( const marObject& from ) { //!\todo TODO } // ------------------------------------------------------------------------- bool marSimpleDicom::save( std::ofstream& os ) { /* int s = _actualStudy.length( ); os.write( ( const char* )&s, sizeof( int ) ); os.write( ( char* )_actualStudy.c_str( ), s * sizeof( char ) ); s = _actualSerie.length( ); os.write( ( const char* )&s, sizeof( int ) ); os.write( ( char* )_actualSerie.c_str( ), s * sizeof( char ) );*/ for(std::vector::iterator it1 = m_DicomInfo.begin(); it1 != m_DicomInfo.end(); it1++ ) { for(std::vector::iterator it2 = it1->begin(); it2 != it1->end(); it2++ ) { for(int i=0; iGetCount(); i++) { //wxString temp( it2->Item(i) ); os << it2->Item(i); //temp; } } } return( true ); } // ------------------------------------------------------------------------- bool marSimpleDicom::load( std::ifstream& is ) { /* int s; reset( ); is.read( ( char* )&s, sizeof( int ) ); _actualStudy.resize( s ); is.read( ( char* )_actualStudy.c_str( ), s * sizeof( char ) ); is.read( ( char* )&s, sizeof( int ) ); _actualSerie.resize( s ); is.read( ( char* )_actualSerie.c_str( ), s * sizeof( char ) );*/ for(std::vector::iterator it1 = m_DicomInfo.begin(); it1 != m_DicomInfo.end(); it1++ ) { for(std::vector::iterator it2 = it1->begin(); it2 != it1->end(); it2++ ) { for(int i=0; iGetCount(); i++) { std::string temp; is >> temp ; it2->Add( temp.c_str() ); } } } return( true ); } // ------------------------------------------------------------------------- void marSimpleDicom::SetVolumeData( int dimX, int dimY, int dimZ, float spacingX, float spacingY, float spacingZ, float rescaleSlope, float rescaleIntercept, unsigned short *pixels) { m_pSlicesPixels = pixels; m_nX = dimX; m_nY = dimY; m_nSlices = dimZ; m_PixelSpacingColumn= spacingX; m_PixelSpacingRow = spacingY; m_PixelSpacingSlice = spacingZ; m_RescaleSlope = rescaleSlope; m_RescaleIntercept = rescaleIntercept; } // ------------------------------------------------------------------------- void marSimpleDicom::reset( ){ // freeVolume( ); //_imageFileNames.Clear( ); int n = m_DicomInfo.size(); ///\todo does this clear evrythg ? for(int i=0; i if(helper) // PS -> delete helper; // PS -> helper = NULL; m_pSlicesPixels = NULL; m_nX = 0; m_nY = 0; m_nSlices = 0; m_PixelSpacingColumn= 0; m_PixelSpacingRow = 0; m_PixelSpacingSlice = 0; m_RescaleSlope = 1; m_RescaleIntercept = 0; } // ------------------------------------------------------------------------- void marSimpleDicom::SetInvestSliceOrder(bool investSliceOrder){ _investSliceOrder=investSliceOrder; } // -------------------------------------------------------------------------