function Xc=FecgImpArtCanc(X,fs,cName,graph,dbFlag)
% -------------------------------------------------------------------------------------
% Impulsive artifact removal from each ECG channel,
% applying the function "ImpArtElimS" to each channel.
%
% Xc=FecgImpArtCanc(X,fs,cName,graph,dbFlag)
%
% X       : input signal matrix (one signal per column)
% fs      : sampling frequency
% cName   : record name
% graph   : flag enabling figure drawing
% dbFlag  : flag enabling figure drawing for debugging
%
% Xc      : output cleaned signal matrix (one signal per column)
%
%   Maurizio Varanini, Clinical Physiology Institute, CNR, Pisa, Italy
%   For any comment or bug report, please send e-mail to: maurizio.varanini@ifc.cnr.it
% -------------------------------------------------------------------------------------

if(nargin<3), cName=''; end
if(nargin<4), graph=0; end
if(nargin<5), dbFlag=0; end

[ndt, ns]=size(X);
%-------------------------------------------------------------------------------
% Set the first ten samples to the median of the following three
npti=10; nptim=3;
for is=1:ns
    X(1:npti,is)=median(X(npti+1:npti+nptim,is));
end
%-------------------------------------------------------------------------------
thE=4;               %  threshold on absolute error
wm=fix(0.06*fs);   %  window length for "deviation" estimation (60ms)
pvsc=2;
Xc=zeros(size(X));
for is=1:ns
    Xc(:,is)=ImpArtElimS(X(:,is),thE,wm,pvsc,dbFlag);
end

if(graph)
    vtime= [1:ndt]/fs;
    figure, set(gcf,'Color','white');
    for is=1:ns,
        subplot(ns,1,is), hold on, plot(vtime,X(:,is),'r'), plot(vtime,Xc(:,is),'b');
        wgmimaV=mimaxscG(Xc(:,is),0,0,.1);
        ylim(wgmimaV);
        if(is~=ns), set(gca,'XTickLabel',''); end
        if(is==1), title([cName,': original (r) & cleaned (b) ECG']); end
        shg
    end
end
end %== function ================================================================
%