ing">'_airway.nii.gz'</span>];
skel_name = [dataname, <span class="string">'_airway_PTKskel.nii.gz'</span>];

<span class="comment">% Load data into relevant formations from nifti images.</span>
<span class="comment">% Load the source CT metadata information</span>
meta = niftiinfo(CT_name);
<span class="comment">% Load CT data as double</span>
source = double(niftiread(meta));

<span class="comment">% Load Airway segmentation and its skeleton as logicals</span>
seg = logical(niftiread(seg_name));
skel = lo