-//--------------------------------------------------------------------
-template <class ImageType>
-void
-clitk::ExtractLymphStationsFilter<ImageType>::
-ExtractStation_3P_LR_inf_Limits()
-{
- /*
- "On the right side, the limit is defined by the air–soft-tissue
- interface. On the left side, it is defined by the air–tissue
- interface superiorly (Fig. 2A–C) and the aorta inferiorly
- (Figs. 2D–I and 3A–C)."
-
- inf : not Right to Azygousvein
- */
- StartNewStep("[Station 3P] Left/Right limits (inferior part) with Azygousvein and Aorta");
-
- // Load structures
- MaskImagePointer AzygousVein = GetAFDB()->template GetImage<MaskImageType>("AzygousVein");
- MaskImagePointer Aorta = GetAFDB()->template GetImage<MaskImageType>("Aorta");
-
- DD("ici");
- writeImage<MaskImageType>(m_Working_Support, "ws.mhd");
-
- typedef clitk::AddRelativePositionConstraintToLabelImageFilter<MaskImageType> RelPosFilterType;
- typename RelPosFilterType::Pointer relPosFilter = RelPosFilterType::New();
- relPosFilter->VerboseStepFlagOff();
- relPosFilter->WriteStepFlagOff();
- relPosFilter->SetBackgroundValue(GetBackgroundValue());
- relPosFilter->SetInput(m_Working_Support);
- relPosFilter->SetInputObject(AzygousVein);
- relPosFilter->AddOrientationTypeString("R");
- relPosFilter->SetInverseOrientationFlag(true);
- relPosFilter->SetIntermediateSpacing(3);
- relPosFilter->SetIntermediateSpacingFlag(false);
- relPosFilter->SetFuzzyThreshold(0.7);
- relPosFilter->AutoCropFlagOn();
- relPosFilter->Update();
- m_Working_Support = relPosFilter->GetOutput();
-
- DD("la");
- writeImage<MaskImageType>(m_Working_Support, "before-L-aorta.mhd");
- relPosFilter->SetInput(m_Working_Support);
- relPosFilter->SetInputObject(Aorta);
- relPosFilter->AddOrientationTypeString("L");
- relPosFilter->SetInverseOrientationFlag(true);
- relPosFilter->SetIntermediateSpacing(3);
- relPosFilter->SetIntermediateSpacingFlag(false);
- relPosFilter->SetFuzzyThreshold(0.7);
- relPosFilter->AutoCropFlagOn();
- relPosFilter->Update();
- m_Working_Support = relPosFilter->GetOutput();
- writeImage<MaskImageType>(m_Working_Support, "after-L-aorta.mhd");
-
- StopCurrentStep<MaskImageType>(m_Working_Support);
- m_ListOfStations["3P"] = m_Working_Support;
-}
-//--------------------------------------------------------------------