matRad_getPhotonLQMParameters

 matRad function to receive the photon LQM reference parameter 
 
 call
   [ax,bx] = matRad_getPhotonLQMParameters(cst,numVoxel,ctScen,VdoseGrid)

 input
   cst:        matRad cst struct
   numVoxel:   number of voxels of the dose cube
   ctScen:     CT scenarios for alpha_x and beta_x should be calculated
   VdoseGrid:  optional linear index vector that allows to specify subindices
               for which ax and bx will be computed

 output
   ax:         vector containing for each linear voxel index alpha_x
   bx:         vector containing for each linear voxel index beta_x

 References
   -

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 Copyright 2018 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.

 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

0001 function [ax,bx] = matRad_getPhotonLQMParameters(cst,numVoxel,ctScen,VdoseGrid)
0002 % matRad function to receive the photon LQM reference parameter
0003 %
0004 % call
0005 %   [ax,bx] = matRad_getPhotonLQMParameters(cst,numVoxel,ctScen,VdoseGrid)
0006 %
0007 % input
0008 %   cst:        matRad cst struct
0009 %   numVoxel:   number of voxels of the dose cube
0010 %   ctScen:     CT scenarios for alpha_x and beta_x should be calculated
0011 %   VdoseGrid:  optional linear index vector that allows to specify subindices
0012 %               for which ax and bx will be computed
0013 %
0014 % output
0015 %   ax:         vector containing for each linear voxel index alpha_x
0016 %   bx:         vector containing for each linear voxel index beta_x
0017 %
0018 % References
0019 %   -
0020 %
0021 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0022 %
0023 % Copyright 2018 the matRad development team.
0024 %
0025 % This file is part of the matRad project. It is subject to the license
0026 % terms in the LICENSE file found in the top-level directory of this
0027 % distribution and at https://github.com/e0404/matRad/LICENSES.txt. No part
0028 % of the matRad project, including this file, may be copied, modified,
0029 % propagated, or distributed except according to the terms contained in the
0030 % LICENSE file.
0031 %
0032 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0033 
0034 ax = zeros(numVoxel,ctScen);
0035 bx = zeros(numVoxel,ctScen);
0036 
0037 for i = 1:size(cst,1)
0038    if isequal(cst{i,3},'OAR') || isequal(cst{i,3},'TARGET')
0039       for ctScen = 1:numel(cst{i,4})
0040          
0041          if exist('VdoseGrid','var')
0042             isInVdoseGrid = ismember(VdoseGrid,cst{i,4}{ctScen});
0043             ax(VdoseGrid(isInVdoseGrid)) = cst{i,5}.alphaX;
0044             bx(VdoseGrid(isInVdoseGrid)) = cst{i,5}.betaX;
0045          else
0046             ax(cst{i,4}{ctScen},ctScen) = cst{i,5}.alphaX;
0047             bx(cst{i,4}{ctScen},ctScen) = cst{i,5}.betaX;
0048          end
0049          
0050       end
0051    end
0052 end
0053 
0054 
0055 
0056