separate.m

function [VisRef,OccRef] = separate(Ref,Mask)
%SEPARATE Summary of this function goes here
%   Detailed explanation goes here

[Offset,ColVis,LabVis] = SR.SetToFields(Ref);
numLab = numel(LabVis);

ColOcc = ColVis;
LabOcc = LabVis;

for lab = 1:numLab
    LabOcc{lab} = LabVis{lab}==lab & ~Mask{lab};
    LabVis{lab} = LabVis{lab}==lab &  Mask{lab};
    
    ColVis  {lab}( ~LabVis{lab}) = nan;
    ColOcc  {lab}(~~LabVis{lab}) = nan;
    
    LabVis{lab} = lab*double(LabVis{lab});
    LabOcc{lab} = lab*double(LabOcc{lab});
end

VisRef = SR.FieldsToSet(Offset,ColVis,LabVis);
OccRef = SR.FieldsToSet(Offset,ColOcc,LabOcc);

end