[Insight-users] DICOM Registration

[uclao]

Hi,

I'm trying to register two sets of DICOM CT images locates in two separated
directories. What I did is taking the sample code ImageRegistration8.cxx and
replacing the image readers with image series readers, copied from
http://www.itk.org/Wiki/ITK/Examples/DICOM/ResampleDICOM. Everything seems
to be working until it gets to the registration part. Once it get to that
line the program simply crushes with exit message "R6010 -abort() has been
called". When I do fixedImageReader -> Print(std::cout) , I have the
following printed:

ImageSeriesReader (0431BC50)
  RTTI typeinfo:   class itk::ImageSeriesReader<class itk::Image<float,3> >
  Reference Count: 1
  Modified Time: 116
  Debug: Off
  Observers: 
    none
  No Inputs
  Outputs: 
    Primary: (04125EA8)
  NumberOfIndexedInputs: 0
  NumberOfRequiredInputs: 0
  NumberOfIndexedOutputs: 1
  NumberOfRequiredOutputs: 1
  Number Of Threads: 4
  ReleaseDataFlag: Off
  ReleaseDataBeforeUpdateFlag: Off
  AbortGenerateData: Off
  Progress: 1
  Multithreader: 
    RTTI typeinfo:   class itk::MultiThreader
    Reference Count: 1
    Modified Time: 98
    Debug: Off
    Observers: 
      none
    Thread Count: 4
    Global Maximum Number Of Threads: 128
    Global Default Number Of Threads: 4
  ReverseOrder: 0
  UseStreaming: 1
  ImageIO: 
    GDCMImageIO (0411D5D8)
      RTTI typeinfo:   class itk::GDCMImageIO
      Reference Count: 3
      Modified Time: 2977
      Debug: Off
      Observers: 
        none
      AbortGenerateData: Off
      Progress: 0
      FileName: SS0001/S02A01_16/1.2.826.0.1.3680043.2.656.4.1.1.19.1.dcm
      FileType: Binary
      ByteOrder: LittleEndian
      IORegion: 
        ImageIORegion (0411D650)
          Dimension: 3
          Index: 0 0 0 
          Size: 512 512 1 
      Number of Components/Pixel: 1
      Pixel Type: scalar
      Component Type: int
      Dimensions: ( 512 512 1 )
      Origin: ( -200.6 -160 -29.5 )
      UseCompression: Off
      UseStreamedReading: On
      UseStreamedWriting: Off
      Internal Component Type: unknown
      RescaleSlope: 1
      RescaleIntercept: -1024
      KeepOriginalUID:Off
      UIDPrefix: 1.2.826.0.1.3680043.2.1125.1
      StudyInstanceUID: 
      SeriesInstanceUID: 
      FrameOfReferenceInstanceUID: 
      CompressionType:1
      Patient Name:BX0002
      Patient ID:A16 
      ...
      Study Date:20000101
      Modality:CT
      Manufacturer:GE MEDICAL SYSTEMS
      Institution Name:
      Model:LightSpeed Ultra
      Scan Options:HELICAL MODE
  MetaDataDictionaryArrayMTime: 2989
  MetaDataDictionaryArrayUpdate: 1

Similar for the fixedImageReader print. The ITK version I am using is 4.1
and compiled with MS Visual C++ 2010.

I went through the user guide Ch. 8 but couldn't find any example related to
the problem directly. It would be nice if anyone can also give one or two
sample code on DICOM registration.

Appreciate your help in this matter. 

Regards, 
Jeff

--
View this message in context: http://itk-insight-users.2283740.n2.nabble.com/DICOM-Registration-tp7519686.html
Sent from the ITK Insight Users mailing list archive at Nabble.com.