Add post limits with VertebralBody and ant limits with Sternum.

[clitk.git] / segmentation / clitkExtractLymphStation_3A.txx

diff --git a/segmentation/clitkExtractLymphStation_3A.txx b/segmentation/clitkExtractLymphStation_3A.txx
index 3e7005906131aafe6eb42fa1ca6aa2436a5fc20e..a55071c19e5874d3f9804e03befa2a1509bff152 100644 (file)
--- a/segmentation/clitkExtractLymphStation_3A.txx
+++ b/segmentation/clitkExtractLymphStation_3A.txx
@@ -1,13 +1,12 @@
 
-#include <itkBinaryDilateImageFilter.h>
-#include <itkMirrorPadImageFilter.h>
-
 //--------------------------------------------------------------------
 template <class ImageType>
 void 
 clitk::ExtractLymphStationsFilter<ImageType>::
 ExtractStation_3A_SetDefaultValues()
 {
+  SetFuzzyThreshold("3A", "Sternum", 0.5);
+  SetFuzzyThreshold("3A", "SubclavianArtery", 0.5);
 }
 //--------------------------------------------------------------------
 
@@ -56,7 +55,27 @@ ExtractStation_3A_Ant_Limits()
   MaskImagePointer Sternum = GetAFDB()->template GetImage<MaskImageType>("Sternum");
   m_Working_Support = 
     clitk::SliceBySliceRelativePosition<MaskImageType>(m_Working_Support, Sternum, 2, 
-                                                       0.5, "PostTo", 
+                                                       GetFuzzyThreshold("3A", "Sternum"), "PostTo", 
+                                                       false, 3, true, false);
+  StopCurrentStep<MaskImageType>(m_Working_Support);
+  m_ListOfStations["3A"] = m_Working_Support;
+}
+//--------------------------------------------------------------------
+
+
+//--------------------------------------------------------------------
+template <class ImageType>
+void 
+clitk::ExtractLymphStationsFilter<ImageType>::
+ExtractStation_3A_Post_Limits() 
+{
+  StartNewStep("[Station 3A] Post limits with SubclavianArtery");
+
+  // Get Sternum, keep posterior part.
+  MaskImagePointer SubclavianArtery = GetAFDB()->template GetImage<MaskImageType>("SubclavianArtery");
+  m_Working_Support = 
+    clitk::SliceBySliceRelativePosition<MaskImageType>(m_Working_Support, SubclavianArtery, 2, 
+                                                       GetFuzzyThreshold("3A", "SubclavianArtery"), "AntTo", 
                                                        false, 3, true, false);
   StopCurrentStep<MaskImageType>(m_Working_Support);
   m_ListOfStations["3A"] = m_Working_Support;