% Author: Dr. Joshua Morgan and Isabella Bagdasarian, TIME lab, UCR
% Filter and enhance the phalloidin image. 
%PARAMETERS
%           P1 - NxN image (normally already filtered)

function [Po, Pth] = CytCleanUp(P1)
%Clean up cytosolic channel
%   Detailed explanation goes here
    Poc = imocbr(P1,strel('disk',2));%soft denoising
    Poc2 = imocbr(Poc,strel('disk',10)); %smoothing
    Po = imopen(Poc2,strel('disk',2)); %edge smoothing
    Pth = imtophat(Po, strel('disk',10)); %even out background
end