[vtkusers] Why vtkDICOMFileSorter only find one series

[David Gobbi]

Yes, there is a roundabout way that you can do this, because the reader can
sort the files into stacks.

  vtkSmartPointer<vtkDICOMReader> reader = vtkSmartPointer<
vtkDICOMReader>::New();
  reader->SetFileNames(fileArray);
  reader->UpdateInformation();

  vtkStringArray *stacks = reader->GetStackIDs();
  int numberOfStacks = stacks->GetNumberOfValues();
  int stackNumber = 1;  // for example

  reader->SetDesiredStackID(stacks->GetValue(stackNumber));
  reader->UpdateInformation();

  vtkIntArray *filesInStack = reader->GetFileIndexArray();
  for (int fileNumber = 0; fileNumber < filesInStack->GetNumberOfValues();
fileNumber++)
  {
    std::cout << fileArray->GetValue(filesInStack->GetValue(jj)).c_str() <<
std::endl;
  }

So this code will use the FileIndexArray() to get a sorted list of the
files in the stack.  Once you have the files you want, you can throw away
the reader.  Note that I only call UpdateInformation() on the reader, so it
is only reading the file headers.

 - David


On Wed, Mar 1, 2017 at 7:54 AM, chenjianyyzz <chenjianyyzz at 163.com> wrote:

> Hello David,
>
> greak, thank you very much. further question: in the attached files there
> are 2 stacks, and now I want to find which DICOM files belongs to stack0,
> and which belongs stack1, how can I find it?
> is there any function that I can get the files list for stack0, just like
> how I get the file list for somes series?
>      vtkStringArray *seriesFilenames = sorter->GetFileNamesForSeries(i);
>
>
> Best Regards
> James
>
>
>
>
>
> At 2017-03-01 21:45:47, "David Gobbi" <david.gobbi at gmail.com> wrote:
>
> Hi James,
>
> All of the images in that zip file have the same SeriesInstanceUID.  So
> this is, in fact, a single series that contains two stacks. This is fairly
> common.
>
> The vtkDICOMReader can be used to find the number of stacks:
>
>   vtkSmartPointer<vtkDICOMReader> reader = vtkSmartPointer<vtkDICOMReader
> >::New();
>   reader->SetFileNames(sorter->GetFileNamesForSeries(0));
>   reader->UpdateInformation();
>   vtkStringArray *stacks = reader->GetStackIDs();
>   int numberOfStacks = stacks->GetNumberOfValues();
>
> The reader can also be told which stack to read:
>
>   reader->SetDesiredStackID(stacks->GetValue(0));
>   reader->Update();
>
> With the "Enhanced MR Image IOD" and the "Enhanced CT Image IOD", it is
> even possible for a single _file_ to contain multiple stacks.
>
>  - David
>
>
>
>
>
>
>
> On Wed, Mar 1, 2017 at 3:50 AM, chenjianyyzz <chenjianyyzz at 163.com> wrote:
>
>> Hello David,
>>
>> I have a set  of images (it includes 2 series), but the
>> vtkDICOMFileSorter can only 1 series?
>> I don't know what's the problem, will any body help me out?
>>
>> attached please find the data, and below is the code:
>> ========================================================
>> vtkSmartPointer<vtkStringArray> filenames =
>> vtkSmartPointer<vtkStringArray>::New();
>> filenames->InsertNextValue("A_10_0001_0325.dcm");
>> filenames->InsertNextValue("A_10_0001_0326.dcm");
>> filenames->InsertNextValue("A_10_0001_0327.dcm");
>> filenames->InsertNextValue("A_10_0001_0328.dcm");
>> filenames->InsertNextValue("A_10_0001_0329.dcm");
>> filenames->InsertNextValue("A_10_0001_0330.dcm");
>> filenames->InsertNextValue("A_10_0001_0331.dcm");
>> filenames->InsertNextValue("A_10_0001_0332.dcm");
>> filenames->InsertNextValue("A_10_0001_0333.dcm");
>> filenames->InsertNextValue("A_10_0001_0334.dcm");
>> filenames->InsertNextValue("A_10_0001_0335.dcm");
>> vtkSmartPointer <vtkDICOMFileSorter> sorter = vtkSmartPointer
>> <vtkDICOMFileSorter>::New();
>> sorter->SetInputFileNames(filenames);
>> sorter->Update();
>>
>> // Get number of series
>> int numSeries = sorter->GetNumberOfSeries(); //==> here I expect to get
>> 2, but only 1 got
>> ========================================================
>>
>> actually "A_10_0001_0325.dcm" should be seperate into another series,
>> will any body help me out?
>> thanks in advance
>>
>>
>>
>>
>
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/vtkusers/attachments/20170301/fc34d82c/attachment.html>