[Insight-users] Reading 4D dicom CT with ITK

[Thibaut Merlin]

Hi,

I have the same kind of problem with 4D dicom PET data. I used the
example code from "DicomSeriesReadImageWrite2". I only add a series
restriction in the code using unique fields for each 3D subset, since
the tags used in setUseSeriesDetails(true) were not sufficient to
separate each subset :

(...)
  typedef itk::GDCMSeriesFileNames NamesGeneratorType;
  NamesGeneratorType::Pointer nameGenerator = NamesGeneratorType::New();

  nameGenerator->SetUseSeriesDetails( true );

  // Add restriction using tag specific to each 3D subset.
  nameGenerator->AddSeriesRestriction("0008|0033" );

  nameGenerator->SetDirectory( argv[1] );
(...)

Unfortunately I still get the whole data in one serie UID instead of
30 separate ones. Is there something to add in order to sub-refine the
4D dataset ?

Thanks in advance for any hints.

Best regards

Thibaut




2011/11/30, Zein Salah <zeinsalah at gmail.com>:
> Hi,
>
> I am trying to read -with ITK- a dicom study in which a set of dicom
> series, acquired at successive time steps, are stored. In other words,
> The whole dataset is a 4D CT (3D+t).
>
> The SeriesNamesGenerator already generates the list of series, BUT
> the the names represent series UID', which do NOT say anything about
> timestamps of the series (acquisition time).
>
> Has someone done such thing before and can tell me how to go
> around the problem?
>
>
> Much thanks in advance,
> Zein
> _____________________________________