Pareto navigation

MatRad_Config.m

@@ -164,13 +164,14 @@ function setDefaultProperties(obj)
             %setDefaultProperties set matRad's default computation
             %   properties
             %  input
-
+
+            %Geometric stf options
             obj.propStf.defaultLongitudinalSpotSpacing = 2;
             obj.propStf.defaultAddMargin = true; %expand target for beamlet finding
-
             obj.propStf.defaultBixelWidth = 5;
 
             %Dose Calculation Options
+            obj.propDoseCalc.defaultDoseEngines = {'SVDPB','HongPB'}; %Names for default engines used when no other is given
             obj.propDoseCalc.defaultResolution = struct('x',3,'y',3,'z',3); %[mm]
             obj.propDoseCalc.defaultDosimetricLateralCutOff = 0.995; %[rel.]
             obj.propDoseCalc.defaultGeometricLateralCutOff = 50; %[mm]
@@ -183,8 +184,6 @@ function setDefaultProperties(obj)
             obj.propDoseCalc.defaultSelectVoxelsInScenarios = 'all';
             obj.propDoseCalc.defaultAirOffsetCorrection = true;
 
-            obj.propDoseCalc.defaultDoseEngines = {'SVDPB','HongPB'}; %Names for default engines used when no other is given
-
             % default properties for fine sampling calculation
             obj.propDoseCalc.defaultFineSamplingProperties.sigmaSub = 1;
             obj.propDoseCalc.defaultFineSamplingProperties.N = 2;
@@ -194,12 +193,7 @@ function setDefaultProperties(obj)
             obj.propDoseCalc.defaultNumHistoriesPerBeamlet = 2e4;
             obj.propDoseCalc.defaultNumHistoriesDirect = 1e6;
             obj.propDoseCalc.defaultOutputMCvariance = true;
-
-            obj.propOpt.defaultMaxIter = 500;
-            obj.propOpt.defaultRunDAO = 0;
-            obj.propOpt.defaultRunSequencing = 0;
-            obj.propOpt.defaultClearUnusedVoxels = false;
-
+
             %deprecated monte carlo options
             obj.propMC.ompMC_defaultHistories = 1e6;
             obj.propMC.ompMC_defaultOutputVariance = false;
@@ -218,6 +212,15 @@ function setDefaultProperties(obj)
             obj.propMC.defaultExternalCalculation = false;
             obj.propMC.defaultCalcDij = false;
 
+
+            %Default Optimization Options
+            obj.propOpt.defaultMaxIter = 10000;
+            obj.propOpt.defaultRunDAO = 0;
+            obj.propOpt.defaultRunSequencing = 0;
+            obj.propOpt.defaultClearUnusedVoxels = false;
+
+
+
             obj.disableGUI = false;
 
             obj.defaults.samplingScenarios = 25;
@@ -234,40 +237,23 @@ function setDefaultPropertiesForTesting(obj)
 
             obj.setDefaultProperties();
 
-            obj.logLevel   = 3; %Omit output except errors
-
+            obj.logLevel = 1; %Omit output except errors
+
+            %Geomteric stf options
             obj.propStf.defaultLongitudinalSpotSpacing = 20;
-            obj.propStf.defaultAddMargin = true; %expand target for beamlet finding
-
             obj.propStf.defaultBixelWidth = 20;
 
+            %Dose Calculation Options
             obj.propDoseCalc.defaultResolution = struct('x',5,'y',6,'z',7); %[mm]
             obj.propDoseCalc.defaultGeometricLateralCutOff = 20;
             obj.propDoseCalc.defaultDosimetricLateralCutOff = 0.8;
             obj.propDoseCalc.defaultKernelCutOff = 20; %[mm]
             obj.propDoseCalc.defaultSsdDensityThreshold = 0.05;
-            obj.propDoseCalc.defaultUseGivenEqDensityCube = false; %Use the given density cube ct.cube and omit conversion from cubeHU.
-            obj.propDoseCalc.defaultIgnoreOutsideDensities = true;
-            obj.propDoseCalc.defaultUseCustomPrimaryPhotonFluence = false; %Use a custom primary photon fluence
-            obj.propDoseCalc.defaultCalcLET = true; %calculate LET for particles
-            obj.propDoseCalc.defaultSelectVoxelsInScenarios = 'all';
-
-            % default properties for fine sampling calculation
-            obj.propDoseCalc.defaultFineSamplingProperties.sigmaSub = 1;
-            obj.propDoseCalc.defaultFineSamplingProperties.N = 2;
-            obj.propDoseCalc.defaultFineSamplingProperties.method = 'fitCircle';
-
-            obj.propOpt.defaultMaxIter = 10;
-            obj.propOpt.defaultClearUnusedVoxels = false;
-
-            obj.propDoseCalc.defaultDoseEngines = {'SVDPB','HongPB'}; %Names for default engines used when no other is given
-
+
             %Monte Carlo options
             obj.propDoseCalc.defaultNumHistoriesPerBeamlet = 100;
             obj.propDoseCalc.defaultNumHistoriesDirect = 100;
             obj.propDoseCalc.defaultOutputMCvariance = true;
-
-            obj.propOpt.defaultMaxIter = 10;
 
             %Deprecated Monte Carlo options
             obj.propMC.ompMC_defaultHistories = 100;
@@ -276,10 +262,8 @@ function setDefaultPropertiesForTesting(obj)
             % Set default histories for MonteCarlo
             obj.propMC.defaultNumHistories = 100;
 
-            obj.propMC.default_photon_engine = 'matRad_OmpConfig';
-            %             obj.propMC.default_photon_engine = 'matRad_TopasConfig';
-            obj.propMC.default_proton_engine = 'matRad_MCsquareConfig';
-            obj.propMC.default_carbon_engine = 'matRad_TopasConfig';
+            % default optimization options
+            obj.propOpt.defaultMaxIter = 10;
 
             obj.defaults.samplingScenarios = 2;
 

Pareto/NavigationUI/matRad_UIData.m

@@ -0,0 +1,66 @@
+classdef matRad_UIData < handle
+    % matRad_UIData implements a class that allows easy storing of
+    % variables related to the pareto Navigation
+    %
+    % References
+    %   -
+    %
+    % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+    %
+    % Copyright 2020 the matRad development team. 
+    % 
+    % This file is part of the matRad project. It is subject to the license 
+    % terms in the LICENSE file found in the top-level directory of this 
+    % distribution and at https://github.com/e0404/matRad/LICENSES.txt. No part 
+    % of the matRad project, including this file, may be copied, modified, 
+    % propagated, or distributed except according to the terms contained in the 
+    % LICENSE file.
+    %
+    % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+    properties
+        wRef %weight vector of last calculated plan
+        fRef %objective function values of last calculated plan
+        fIndsAll % stores all objective function values
+        fIndsRed % stores only the "reduced points"
+        upboundInit
+        upboundRed
+        upboundSlider
+        lowboundSliderInit
+        lowboundSlider
+        linestyle
+    end
+
+    methods
+        function obj = matRad_UIData(wRef,fRef,fInds)
+            obj.wRef = wRef; 
+            obj.fRef = fRef; 
+            obj.fIndsAll = fInds;
+            obj.fIndsRed = fInds;
+            obj.upboundInit = max(fInds,[],1);
+            obj.upboundRed = max(fInds,[],1);
+            obj.upboundSlider= max(fInds,[],1);
+            obj.lowboundSliderInit = min(fInds,[],1);
+            obj.lowboundSlider = min(fInds,[],1);
+            obj.linestyle = 2;
+        end       
+
+        function [sliderLowBound,sliderUpBound] = restrictObjective(obj,i,bound)
+            obj.upboundRed(i) = bound;
+            obj.fIndsRed = obj.fIndsRed(obj.fIndsRed(:,i) <= bound,:);
+            obj.upboundSlider= max([obj.fIndsRed;obj.fRef],[],1);
+            obj.upboundSlider(i) = bound;
+            obj.lowboundSlider = min([obj.fIndsRed;obj.fRef],[],1);
+            sliderLowBound = obj.lowboundSlider;
+            sliderUpBound = obj.upboundSlider;
+        end
+
+        function releaseObjectiveBounds(obj)
+            obj.upboundRed = obj.upboundInit;
+            obj.fIndsRed = obj.fIndsAll;    
+            obj.upboundSlider = obj.upboundInit;
+            obj.lowboundSlider = obj.lowboundSliderInit;
+        end
+
+    end
+end