2 #include <itkBinaryDilateImageFilter.h>
3 #include <itkMirrorPadImageFilter.h>
5 //--------------------------------------------------------------------
6 template <class ImageType>
8 clitk::ExtractLymphStationsFilter<ImageType>::
9 ExtractStation_3A_SetDefaultValues()
12 //--------------------------------------------------------------------
15 //--------------------------------------------------------------------
16 template <class ImageType>
18 clitk::ExtractLymphStationsFilter<ImageType>::
19 ExtractStation_3A_SI_Limits()
21 // Apex of the chest or Sternum & Carina.
22 StartNewStep("[Station 3A] Inf/Sup limits with Sternum and Carina");
24 // Get Carina position (has been determined in Station8)
25 m_CarinaZ = GetAFDB()->GetDouble("CarinaZ");
27 // Get Sternum and search for the upper position
28 MaskImagePointer Sternum = GetAFDB()->template GetImage<MaskImageType>("Sternum");
30 // Search most sup point
31 MaskImagePointType ps = Sternum->GetOrigin(); // initialise to avoid warning
32 clitk::FindExtremaPointInAGivenDirection<MaskImageType>(Sternum, GetBackgroundValue(), 2, false, ps);
33 double m_SternumZ = ps[2]+Sternum->GetSpacing()[2]; // One more slice, because it is below this point
37 clitk::CropImageAlongOneAxis<MaskImageType>(m_Working_Support, 2,
38 m_CarinaZ, m_SternumZ, true,
39 GetBackgroundValue());
41 StopCurrentStep<MaskImageType>(m_Working_Support);
42 m_ListOfStations["3A"] = m_Working_Support;
44 //--------------------------------------------------------------------
47 //--------------------------------------------------------------------
48 template <class ImageType>
50 clitk::ExtractLymphStationsFilter<ImageType>::
51 ExtractStation_3A_Ant_Limits()
53 StartNewStep("[Station 3A] Ant limits with Sternum");
55 // Get Sternum, keep posterior part.
56 MaskImagePointer Sternum = GetAFDB()->template GetImage<MaskImageType>("Sternum");
58 clitk::SliceBySliceRelativePosition<MaskImageType>(m_Working_Support, Sternum, 2,
60 false, 3, true, false);
61 StopCurrentStep<MaskImageType>(m_Working_Support);
62 m_ListOfStations["3A"] = m_Working_Support;
64 //--------------------------------------------------------------------