Noninvasive Fetal ECG: The PhysioNet/Computing in Cardiology Challenge 2013 1.0.0

function [fetal_QRSAnn_est,QT_Interval] = physionet2013(tm,ECG,cName,qrsAf,graphDflags)
% Algorithm for Physionet/CinC competition 2013.
%
% [fetal_QRSAnn_est,QT_Interval] = physionet2013(tm,ECG,cName,qrsAf)
%
% inputs:
%   ECG:    60000x4 (4 channels and 1min of signal at 1000Hz) matrix of
%           abdominal ECG channels.
%   tm :    Nx1 vector of time in milliseconds
%   cName:  record name (optional)
%   qrsAf:  QRS markers in seconds (optional, learning set only)
% output:
%   FQRS:  FQRS markers in seconds. Each marker indicates the position of one
%          of the FQRS detected by the algorithm.
%   QT_Interval:   1x1 estimated fetal QT duration (enter NaN or 0 if you do wish to calculate)
%

try
    if(nargin<3), cName=''; end
    if(nargin<4), qrsAf=[]; end
    if(nargin<5)
        dbFlag=0;                   % debug flag
        graph=0;                    % enable/disable graphical representation
        saveFig=0;                  % =1 => save figures of the processing phases
        saveFigRRf=0;               % =1 => save estimated fetal RR figures
    else
        dbFlag=graphDflags.dbFlag;
        graph=graphDflags.graph;
        saveFig=graphDflags.saveFig;
        saveFigRRf=graphDflags.saveFigRRf;
    end
           
    % ---- check size of ECG ----
    if size(ECG,2)>size(ECG,1)
        ECG = ECG';
    end
        
    fs = 1000;             % sampling frequency
    
    % ---- Artifact canceling ----
    % X=FecgFecgImpArtCanc(ECG,fs,cName,graph,dbFlag);
    X=FecgImpArtCanc(ECG,fs,cName,0,0);
    
    % ---- detrending  ----
    % Xd=FecgDetrFilt(X,fs,cName,graph,dbFlag);
    Xd=FecgDetrFilt(X,fs,cName,0,0);
    
    % ---- Power line interference removal by notch filtering ----
    % Xf=FecgNotchFilt(Xd,fs,cName,graph,dbFlag);
    Xf=FecgNotchFilt(Xd,fs,cName,0,0);
    
    % ---- Independent Component Analysis ----
    % Xm=FecgICAm(Xf,fs,cName,graph,dbFlag,saveFig);
    Se=FecgICAm(Xf,fs,cName,graph,dbFlag,saveFig);
    
    % ---- Signal interpolation
    % Xi=FecgInterp(X,fs,interpFact,cName,graph);
    [Se,fs]=FecgInterp(Se,fs,4,cName,0);
    
    % ---- Channel selection and Mother QRS detection
    qrsM=FecgQRSmDet(Se,fs,cName,graph,dbFlag,saveFig,qrsAf);
    
    % ---- Mother QRS cancelling
    Xr=FecgQRSmCanc(Se,qrsM,fs,cName,graph,dbFlag,saveFig,qrsAf);
    
    % ---- Source separation by ICA on residual signals
    Ser=FecgICAf(Xr,fs,cName,graph,dbFlag,saveFig);
    
    % ---- Channel selection and Fetal QRS detection
    % qrsF=FecgQRSfDniAdf(Ser,fs,cName,qrsM,graph,dbFlag,saveFig,saveFigRRf,qrsAf);
    qrsF=FecgQRSfDet(Ser,fs,cName,qrsM,graph,dbFlag,saveFig,saveFigRRf,qrsAf);
    
    fetal_QRSAnn_est=qrsF;
    QT_Interval         = [];
    
catch
    fetal_QRSAnn_est    = [];
    QT_Interval         = [];
end

end %== function ================================================================
%