function [stat, res] = apply_ROIrules(ops, stat0, res0, clustrules)

stat0    = get_regions(stat0, res0);

%%
if ops.splitROIs
    [stat, res] = get_validregions(stat0,res0, clustrules);
else
    stat = stat0;
    res = res0;
    
    for j = 1:length(stat)
        stat(j).igood = 1;
    end
    
    % trim rogue pixels not in any compact subregion. Assign to Nk+1th cluster
    [stat, res] = trimPixels(stat, res);           
    
    
end


if isfield(ops, 'ResultsSavePath')
    if ~exist(ops.ResultsSavePath, 'dir')
        mkdir(ops.ResultsSavePath)
    end
    
    iplane = ops.iplane;
    Nk = ops.Nk;
    
    
    save(sprintf('%s/F_%s_%s_plane%d_Nk%d.mat', ops.ResultsSavePath, ...
        ops.mouse_name, ops.date, iplane, Nk),  'ops', 'res', 'stat', 'stat0', 'res0', 'clustrules')
end
end