[Insight-users] Reading DICOM time series study

[Mathieu Malaterre]

Martin,

	I'll repeat what Mark just told you: I never saw those tags before. 
I'll double check what 3.3 says about wether or not those tags are 
compulsory, but I would rather order the files using the usual Date/Time 
tags. Those I am pretty sure to find everywhere and I am also pretty 
sure that they are correct :)

	BTW is this a dual echo sequences ? I'd would like to know how they 
increment the number in this case ?

Thanks,
Mathieu

Martin Kavec wrote:
> Thanks for the reply, Karthik.
> 
> Your code may, for sure, be a valid solution. However, after I read your 
> message, I told myself there must be an information related to the time 
> series in the DICOM header of an image. So I randomly grabbed an image from 
> an fMRI time series and found following two (correct) tags
> 
> 0020,0100  Temporal Position Identifier: 52
> 0020,0105  Number of Temporal Positions: 120 
> 
> Everything what's needed to construct the 4D image is there, so why not to use 
> it? Quoting myself: Is this a limitation of ITK's DICOM image reader 
> implementation? I would still be tempted to say: yes.
> 
> What do you say?
> 
> Regards,
> Martin
> 
> On Monday 11 July 2005 14:45, Karthik Krishnan wrote:
> 
>>I don't think it is not a limitation of the dicom image reader. You need
>>to do that with your code.
>>
>>You could for instance read all the dicom series in the directory and
>>use the JoinSeriesImageFilter to create a 4D image from the 3D images
>>you just read. It is not a limitation of the dicom reader because a
>>dicom file is  a 2D file. So the ImageSeriesReader which passes the set
>>of filenames to the GDCMImageIO has to be instantiated over an image of
>>dimension 2+1.
>>
>>You could do something like
>>    const SeriesIdContainer & seriesUID = nameGenerator->GetSeriesUIDs();
>>    SeriesIdContainer::const_iterator seriesItr = seriesUID.begin();
>>    SeriesIdContainer::const_iterator seriesEnd = seriesUID.end();
>>    while( seriesItr != seriesEnd )
>>      {
>>      seriesReader->SetFileNames( nameGenerator->GetFileNames(
>>seriesItr->c_str() ) );
>>      seriesReader->Update();
>>      joinSeriesFilter->SetInput( i++, seriesReader->GetOutput() );
>>      seriesItr++;
>>      }
>>
>>my4DImage = joinSeriesFilter->GetOutput();
>>
>>Thanks
>>karthik
>>
>>Martin Kavec wrote:
>>
>>>Hi,
>>>
>>>I have a time series (fMRI, bolus tracking) studies in DICOM format, which
>>>I am trying to read in. I gave a brief trial test to
>>>DicomSeriesReadImageWrite2.cxx from Examples/IO. Although I changed the
>>>Dimension = 4, the output Analyze format image is 3D only. I assume that
>>>internal itkImage dimension is also 3D.
>>>
>>>Is this a limitation of ITK's DICOM image reader implementation?
>>>
>>>Regards.
>>>
>>>Martin
>>>_______________________________________________
>>>Insight-users mailing list
>>>Insight-users at itk.org
>>>http://www.itk.org/mailman/listinfo/insight-users
> 
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users
>