c++ - GDCM を使用した単純な CT DICOM 画像の生成

Question

GDCMライブラリを使用して一連の CT DICOM 画像 (512x512、MONOCHROME) を生成する必要があります。私が取り組んでいるアプリケーションの構造上、イメージごとにgdcm::Imageとインスタンスを個別に準備し、最後にとメソッドを使用してそれらを渡すことが理想的です。gdcm::FilegdcmImageWriterSetFile()SetImage()

は、gdcmImageWriter::write()DICOM イメージの書き込みに失敗します。ヒントとして、 VTKGDCMImageWriterクラスのソース コードとこの GDCM の例を既に確認しました。この問題は、gdcm::Imageとgdcm::Fileインスタンスの間の非互換性 (たとえば、両方のクラスにDICOM 画像の測光解釈gdcmImageWriterに関する情報が含まれている) に起因するようですが、それらが何であるかは正確にはわかりません。

次の簡略化された C++ コードは、私が得たエラーを再現しています。

#include <gdcmFile.h>
#include <gdcmImage.h>
#include <gdcmTag.h>
#include "gdcmDataElement.h"
#include "gdcmByteValue.h"
#include "gdcmPixelFormat.h"
#include "gdcmImageWriter.h"
#include "gdcmAnonymizer.h"

#include <iostream>
#include <string>

InitializeGDCMImage(gdcm::Image*  imagePtr);
InitializeGDCMFile( gdcm::File*   filePtr );
WriteFile( gdcm::Image* imagePtr , gdcm::File* filePtr , std::string  fileName );

int main()
{
    std::string  fileName = "./TEST.dcm";

    // Will be deleted by gdcm::writer.
    gdcm::Image* imagePtr = new gdcm::Image;
    gdcm::File*  filePtr  = new gdcm::File;   

    InitializeGDCMImage(imagePtr);

    InitializeGDCMFile( filePtr );

    WriteFile( imagePtr , filePtr , fileName );

    return 0;
}

どこ、

void WriteFile( gdcm::Image* imagePtr , gdcm::File* filePtr , std::string fileName )
{
    gdcm::ImageWriter* writer = new gdcm::ImageWriter; 
    writer->SetFileName(  fileName.c_str() );
    writer->SetImage(    *imagePtr         );
    writer->SetFile(     *filePtr          );
    if( !writer->Write() ){
      std::cerr << "ERROR: Could not write to \"" << fileName <<  "\"" << std::endl;
    }
    delete writer;
}

WriteFile()gdcm::Fileおよびを設定し、gdcm::Imageを呼び出しますwrite()。 InitializeGDCMFile()を使用していくつかの DICOM タグを変更し、実際のピクセル値gdcm::AnonymizerをInitializeGDCMImage()設定します (簡潔にするためにここではスキップします)。

void InitializeGDCMFile( gdcm::File* origianlFilePtr )
{
  gdcm::File* filePtr = new gdcm::File;
  gdcm::Anonymizer anon;
  anon.SetFile( *filePtr );  

  anon.Replace( gdcm::Tag(0x0028,0x0002) , "1" );                       //SamplesperPixel
  anon.Replace( gdcm::Tag(0x0028,0x0004) , "MONOCHROME2" );             //PhotometricInterpretation
  anon.Replace( gdcm::Tag(0x0028,0x0010) , "512" );                     //Rows
  anon.Replace( gdcm::Tag(0x0028,0x0011) , "512" );                     //Columns
  anon.Replace( gdcm::Tag(0x0028,0x0030) , "0.781\\0.781" );            //PixelSpacing

  anon.Replace( gdcm::Tag(0x0028,0x1050) , "0" );                        //WindowCenter
  anon.Replace( gdcm::Tag(0x0028,0x1051) , "100" );                      //WindowWidth
  anon.Replace( gdcm::Tag(0x0028,0x1052) , "0" );                        //RescaleIntercept
  anon.Replace( gdcm::Tag(0x0028,0x1053) , "1" );                        //RescaleSlope

  *origianlFilePtr = *filePtr;
}


void InitializeGDCMImage(gdcm::Image* imagePtr)
{
  imagePtr->SetPhotometricInterpretation( gdcm::PhotometricInterpretation::MONOCHROME2 );

  imagePtr->SetNumberOfDimensions(2);
  unsigned int dims[2]={512,512};
  imagePtr->SetDimensions( dims );

  imagePtr->SetSpacing( 0 , 1. );
  imagePtr->SetSpacing( 1 , 1. );

  imagePtr->SetIntercept(0.);
  imagePtr->SetSlope(1.);

  double dirCos[6]={1.,0.,0.,0.,1.,0.};
  imagePtr->SetDirectionCosines( dirCos );

  char *buffer = new char[512*512*sizeof(int16_t)];
  imagePtr->SetPixelFormat( gdcm::PixelFormat::INT16 );
  imagePtr->GetDataElement().SetByteValue( buffer , 512*512*sizeof(int16_t) );
  imagePtr->GetPixelFormat().SetSamplesPerPixel(1);
  delete[] buffer;
}