and header variables.

% get CT, segmentation and skeleton images of the case
CT_name = [dataname, '_source.nii.gz'];
seg_name = [dataname, '_airway.nii.gz'];
skel_name = [dataname, '_airway_PTKskel.nii.gz'];

% Load data into relevant formations from nifti images.
% Load the source CT metadata information
meta = niftiinfo(CT_name);
% Load CT data as double
source = double(niftiread(meta));

% Load Airway segmentation and its skeleton as logicals
seg = logical(niftiread(seg_name));
skel = logical(niftiread(