Station8 almost ok.

[clitk.git] / segmentation / clitkExtractLymphStation_3P.txx

diff --git a/segmentation/clitkExtractLymphStation_3P.txx b/segmentation/clitkExtractLymphStation_3P.txx
new file mode 100644 (file)
index 0000000..effcd4f
--- /dev/null
+++ b/segmentation/clitkExtractLymphStation_3P.txx
@@ -0,0 +1,55 @@
+
+#include <itkBinaryDilateImageFilter.h>
+#include <itkMirrorPadImageFilter.h>
+
+//--------------------------------------------------------------------
+template <class ImageType>
+void 
+clitk::ExtractLymphStationsFilter<ImageType>::
+ExtractStation_3P_SetDefaultValues()
+{
+}
+//--------------------------------------------------------------------
+
+//--------------------------------------------------------------------
+template <class ImageType>
+void 
+clitk::ExtractLymphStationsFilter<ImageType>::
+ExtractStation_3P_SI_Limits() 
+{
+  /*
+    Apex of the chest & Carina.
+   */
+  StartNewStep("[Station 3P] Inf/Sup limits with apex of the chest and carina");
+
+  writeImage<MaskImageType>(m_Working_Support, "support.mhd");
+
+  // Get Carina position (has been determined in Station8)
+  m_CarinaZ = GetAFDB()->GetDouble("CarinaZ");
+  DD(m_CarinaZ);
+  
+  // Get Apex of the Chest
+  MaskImagePointer Lungs = GetAFDB()->template GetImage<MaskImageType>("Lungs");
+  DD("lung ok");
+  MaskImagePointType p;
+  bool found = clitk::FindExtremaPointInAGivenDirection<MaskImageType>(Lungs, 
+                                                                       GetBackgroundValue(), 
+                                                                       1, true, p);
+  DD(found);
+  DD(p);
+  double m_ApexOfTheChest = p[2];
+  DD(m_ApexOfTheChest);
+
+  /* Crop support :
+     Superior limit = carina
+     Inferior limit = Apex of the chest */
+  m_Working_Support = 
+    clitk::CropImageAlongOneAxis<MaskImageType>(m_Working_Support, 2, 
+                                                m_ApexOfTheChest,
+                                                m_CarinaZ, true,
+                                                GetBackgroundValue());
+
+  StopCurrentStep<MaskImageType>(m_Working_Support);
+  m_ListOfStations["3P"] = m_Working_Support;
+}
+//--------------------------------------------------------------------