matRad inverse planning wrapper function call [vAlpha, vBeta] = matRad_calcLQParameter(vRadDepths,mTissueClass,baseData) input vRadDepths: radiological depths of voxels mTissueClass: tissue classes of voxels baseData: biological base data output vAlpha: alpha values for voxels interpolated from base data vBeta: beta values for voxels interpolated from base data References - %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Copyright 2015 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 [vAlpha, vBeta] = matRad_calcLQParameter(vRadDepths,mTissueClass,baseData) 0002 % matRad inverse planning wrapper function 0003 % 0004 % call 0005 % [vAlpha, vBeta] = matRad_calcLQParameter(vRadDepths,mTissueClass,baseData) 0006 % 0007 % input 0008 % vRadDepths: radiological depths of voxels 0009 % mTissueClass: tissue classes of voxels 0010 % baseData: biological base data 0011 % 0012 % output 0013 % vAlpha: alpha values for voxels interpolated from base data 0014 % vBeta: beta values for voxels interpolated from base data 0015 % 0016 % References 0017 % - 0018 % 0019 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0020 % 0021 % Copyright 2015 the matRad development team. 0022 % 0023 % This file is part of the matRad project. It is subject to the license 0024 % terms in the LICENSE file found in the top-level directory of this 0025 % distribution and at https://github.com/e0404/matRad/LICENSES.txt. No part 0026 % of the matRad project, including this file, may be copied, modified, 0027 % propagated, or distributed except according to the terms contained in the 0028 % LICENSE file. 0029 % 0030 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0031 0032 vAlpha = NaN*ones(numel(vRadDepths),1); 0033 vBeta = NaN*ones(numel(vRadDepths),1); 0034 0035 numOfTissueClass = size(baseData(1).alpha,2); 0036 0037 % range shift 0038 depths = baseData.depths + baseData.offset; 0039 0040 for i = 1:numOfTissueClass 0041 mask = mTissueClass == i; 0042 if any(mask) 0043 vAlpha(mask) = matRad_interp1(depths,baseData.alpha(:,i),vRadDepths(mask)); 0044 vBeta(mask) = matRad_interp1(depths,baseData.beta(:,i), vRadDepths(mask)); 0045 end 0046 end 0047