% Author: Dr. Joshua Morgan and Isabella Bagdasarian, TIME lab, UCR
% Filter 2D nuclei stain

%PARAMETERS
%           D - NxN image

function [Do] = NucCleanUp(D)
%clean up nuclei channel

    Doc = imocbr(D,strel('disk',2)); %soft denoising
    Dth = imthbr(Doc,strel('disk',50)); %tophat to remove uneven illumination
    Doc2 = imocbr(Dth,strel('disk',2)); %smoothing
    Do = imopen(Doc2,strel('disk',2)); %edge smoothing
    
end