GDCM 2.2.0

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

iU22tomultisc.cxx

    /*=========================================================================

      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.

    =========================================================================*/
    /*
     * iU22 Raw Data extractor
     */
    #include "gdcmReader.h"
    #include "gdcmImageWriter.h"
    #include "gdcmAttribute.h"
    #include "gdcmPrivateTag.h"

    #include <math.h>

    int main(int argc, char *argv [])
    {
      if( argc < 2 ) return 1;
      // IM_001
      const char *filename = argv[1];

      gdcm::Reader reader; // Do not use ImageReader
      reader.SetFileName( filename );
      if( !reader.Read() )
        {
        std::cerr << "Failed to read: " << filename << std::endl;
        return 1;
        }

     // * The data is simply 8-bit unsigned in the obvious x/y/z order
     // * 200D,300B contains the data
     // * 200D,3001 contains the no. of voxels (416,412,256 in this case)
     // * 200D,3003 contains the voxel sizes (0.156184527398215 /
     // 0.1223749613981957 / 0.328479990704639 in this case)

      const gdcm::File &file = reader.GetFile();
      const gdcm::DataSet &ds = file.GetDataSet();
      const gdcm::PrivateTag trawdataus( 0x200d, 0x0b, "Philips US Imaging DD 033" );
      const gdcm::DataElement &rawdataus = ds.GetDataElement( trawdataus );

      const gdcm::PrivateTag tcolsrowsframes( 0x200d, 0x01, "Philips US Imaging DD 036" );
      const gdcm::DataElement &colsrowsframes = ds.GetDataElement( tcolsrowsframes );
      // const gdcm::PrivateTag tcolsrowsframes( 0x200d, 0x02, "Philips US Imaging DD 036" );
      // this is just a duplicate previous tag.
      const gdcm::PrivateTag tvoxelspacing( 0x200d, 0x03, "Philips US Imaging DD 036" );
      const gdcm::DataElement &voxelspacing = ds.GetDataElement( tvoxelspacing );

      gdcm::Element<gdcm::VR::DS,gdcm::VM::VM3> dims; // Use DS to interpret value stored in LO
      dims.SetFromDataElement( colsrowsframes );

      gdcm::Element<gdcm::VR::DS,gdcm::VM::VM3> spacing;
      spacing.SetFromDataElement( voxelspacing );

      gdcm::ImageWriter writer;

      gdcm::Image &image = writer.GetImage();
      image.SetNumberOfDimensions( 3 ); // good default
      image.SetDimension(0, dims[0] );
      image.SetDimension(1, dims[1] );
      image.SetDimension(2, dims[2] );
      image.SetSpacing(0, spacing[0] );
      image.SetSpacing(1, spacing[1] );
      image.SetSpacing(2, spacing[2] );
      gdcm::PixelFormat pixeltype = gdcm::PixelFormat::UINT8;

      gdcm::PhotometricInterpretation pi;
      pi = gdcm::PhotometricInterpretation::MONOCHROME2;
      image.SetPhotometricInterpretation( pi );
      image.SetPixelFormat( pixeltype );

      image.SetDataElement( rawdataus );

      std::string outfilename = "outiu22.dcm";

      gdcm::DataElement de( gdcm::Tag(0x8,0x16) ); // SOP Class UID
      de.SetVR( gdcm::VR::UI );
      gdcm::MediaStorage ms(
        gdcm::MediaStorage::UltrasoundMultiFrameImageStorage );
    //    gdcm::MediaStorage::MultiframeGrayscaleByteSecondaryCaptureImageStorage );
      de.SetByteValue( ms.GetString(), strlen(ms.GetString()));
      writer.GetFile().GetDataSet().Replace( de );

      writer.SetFileName( outfilename.c_str() );
      if( !writer.Write() )
        {
        std::cerr << "could not write: " << outfilename << std::endl;
        return 1;
        }

      return 0;
    }