BUG: fixed ALive-research#248, adjusted region growing and surface re…

…solution settings
RuoyanMeng · Sep 25, 2023 · 2ccac63 · 2ccac63
1 parent c01526c
commit 2ccac63
Show file tree

Hide file tree

Showing 3 changed files with 13 additions and 4 deletions.
diff --git a/Liver/Resources/UI/ResectionsWidget.ui b/Liver/Resources/UI/ResectionsWidget.ui
@@ -543,7 +543,7 @@
             <double>1.000000000000000</double>
            </property>
            <property name="minimum">
-            <double>5.000000000000000</double>
+            <double>1.000000000000000</double>
            </property>
            <property name="maximum">
             <double>10.000000000000000</double>

diff --git a/LiverVolumetry/LiverVolumetry.py b/LiverVolumetry/LiverVolumetry.py
@@ -387,8 +387,14 @@ def computeVolume(self, segmentsVolumeNode, targetSegmentVolumeNode, segmentatio
         segStatLogic.computeStatistics()
         stats = segStatLogic.getStatistics()
         for segmentId in stats["SegmentIDs"]:
-          voxel_count = stats[segmentId,"LabelmapSegmentStatisticsPlugin.voxel_count"]
-          volume_cm3 = stats[segmentId,"LabelmapSegmentStatisticsPlugin.volume_cm3"]
+          voxel_count = 0
+          volume_cm3 = 0
+          if stats[segmentId,"LabelmapSegmentStatisticsPlugin.voxel_count"]:
+            voxel_count = stats[segmentId,"LabelmapSegmentStatisticsPlugin.voxel_count"]
+            volume_cm3 = stats[segmentId,"LabelmapSegmentStatisticsPlugin.volume_cm3"]
+          elif stats[segmentId,"ScalarVolumeSegmentStatisticsPlugin.voxel_count"]:
+            voxel_count = stats[segmentId,"ScalarVolumeSegmentStatisticsPlugin.voxel_count"]
+            volume_cm3 = stats[segmentId,"ScalarVolumeSegmentStatisticsPlugin.volume_cm3"]
           segmentName = segmentationNode.GetSegmentation().GetSegment(segmentId).GetName()
           statistics[segmentId] = [segmentName, voxel_count, volume_cm3]
           self.scl.VolumetryTable(segmentName, targetSegmentVolume, voxel_count, volume_cm3,outputTable)

diff --git a/LiverVolumetry/Logic/vtkLiverVolumetryLogic.cxx b/LiverVolumetry/Logic/vtkLiverVolumetryLogic.cxx
@@ -84,7 +84,7 @@ void vtkLiverVolumetryLogic::PrintSelf(ostream &os, vtkIndent indent)
 void vtkLiverVolumetryLogic::ComputeAdvancedPlanningVolumetry(vtkMRMLLabelMapVolumeNode* SelectedSegmentsLabelMap, vtkMRMLTableNode* OutputTableNode, vtkMRMLMarkupsFiducialNode* ROIMarkersList, vtkCollection* ResectionNodes, double TargetSegmentationVolume){
   vtkLabelMapHelper::LabelMapType::RegionType TargetSegmentsBoundingBox;
   vtkLabelMapHelper::LabelMapType::Pointer TargetSegmentsITKImage;
-  int baseValue = 7;
+  int baseValue = 100;
   double spacing[3];
   SelectedSegmentsLabelMap->GetSpacing(spacing);
 
@@ -116,6 +116,9 @@ void vtkLiverVolumetryLogic::ComputeAdvancedPlanningVolumetry(vtkMRMLLabelMapVol
         {
         auto bezierSurfaceNode = vtkMRMLMarkupsBezierSurfaceNode::SafeDownCast(this->resectionNodes->GetItemAsObject(i));
         auto Res =  GetRes(bezierSurfaceNode, spacing, 300);
+        if(Res < 500){
+          Res  =  500;
+        }
         BezierHR = GenerateBezierSurface(Res, bezierSurfaceNode);
         if(i == 0){
           this->ProjectedTargetSegmentImage = vtkLabelMapHelper::VolumeNodeToItkImage(TargetSegmentLabelMapCopy, true, false);