function ecgpuwave(varargin)
ecgpuwave(recordName,annFileName,startTime,stopTime,qrsAnn,pflag,signalList)
Wrapper to the ECGPUWAVE binary written by Pablo Laguna (laguna@posta.unizar.es), Raimon Jan[�e], Eudald Bogatell, and David Vigo Anglada:
http://www.physionet.org/physiotools/wag/ecgpuw-1.htm
ECGPUWAVE analyses an ECG signal from the specified record 'recordName', detecting the QRS complexes
and locating the beginning, peak, and end of the P, QRS, and ST-T waveforms. The output of ecgpuwave
is written as a standard WFDB-format annotation file associated
with the specified annotator 'outputAnnFileName'.
The QRS detector is based on the algorithm of Pan and Tompkins (reference 1) with some improvements that
make use of slope information (reference 2). Optionally, QRS annotations can be provided as input (see option 'qrsAnn'),
permitting the use of external QRS detectors such as sqrs(1) or manually-edited annotations
(which can be created using wave(1) ). The waveform limit locator is based on the algorithm
described in reference 3 and evaluated in references 3 and 4.
The output annotation file contains PWAVE ("p") and TWAVE ("t") annotations
that indicate the P- and T-wave peaks, as well as QRS annotations (NORMAL ("N") if generated
by the built-in QRS detector, or copies of the input QRS annotations if these were supplied).
ECGPUWAVE classifies each T wave as type 0 (normal), 1 (inverted), 2 (positive monophasic),
3 (negative monophasic), 4 (biphasic negative-positive), or 5 (biphasic positive-negative);
this numeric classification is written into the num field of each TWAVE annotation.
The P, QRS, and T waveform onsets and ends are marked in the output annotation file using
WFON ("(") and WFOFF (")") annotations. The num field of each WFON and WFOFF annotation designates
the type of waveform with which it is associated: 0 for a P wave, 1 for a QRS complex, or 2 for a T wave.
WARNING:
If ECGPUWAVE is used without providing an annotator file for the R peaks it relies on a Pan-Tompkins detector.
This detector may not be optimized for your particular signal, and thus may need to be tuned, yielding anempty annotation
file. In this case, you may want to first run a QRS detector to generate an appropriate QRS annotation file,
which you can then check manually prior to calling ECGPUWAVE with the
'qrsAnn' option described below.
IMPORTANT NOTE:
A patch to the original Fortran source files was provided by Roberto Sassi in
order to compile ECGPUWAVE in 64-bit arch. This patch requires
compiling the source code with gfortran (g77 is not supported in gcc 4.x), which can
yield slighly different results. For more information please see:
http://www.dti.unimi.it/~sassi/software/wfdb64HowTo.htm
References:
1. Pan J and Tompkins WJ. A Real-Time QRS Detection Algorithm. IEEE Transactions on Biomedical Engineering 32(3):230-236, 1985.
2. Laguna P. New Electrocardiographic Signal Processing Techniques: Application to Long-term Records. Ph. D. dissertation, Science Faculty, University of Zaragoza, 1990.
3. Laguna P, Jan[�e] R, Caminal P. Automatic Detection of Wave Boundaries in Multilead ECG Signals: Validation with the CSE Database. Computers and Biomedical Research 27(1):45-60, 1994.
4. Jan[�e] R, Blasi A, Garc[�i]a J, and Laguna P. Evaluation of an automatic threshold based detector of waveform limits in Holter ECG with the QT database. Computers in Cardiology 24:295-298 (1997; available at http://www.physionet.org/physiobank/database/qtdb/eval/ )
Required Parameters:
recordName
String specifying the name of the record to annotate with ECGPUWAVE.
annFileName
String specifying the output file name that ECPUWAVE will store its
output.
ecgpuwave(recordName,outputAnnFileName,startTime,stopTime,qrsAnn,pflag,signalList)
Optional Parameters are:
startTime
Either WFDB Time String or 1x1 integer specifying the sample number
specifying the starting location.
stopTime
Either WFDB Time String or 1x1 integer specifying the sample number
specifying the ending location.
qrsAnn
A String represengint the QRS annotation file that ECGPUWAVE should
use (must be in the WFDP PATH or within the current directory).
pflag
A boolean (default = 0), if true all beats are processed, if false,
only normal beats are used.
signalList
A Nx1 vector of integers specifying which signals in 'recordName' are
to be analyzed (default = [0].
Wrapper written by Ikaro Silva, 2013
Last Modified: May 29 , 2014
Version 0.0.1
%Example - Will go into the WFDB directory of examples to load data
curdir=pwd;
[~,config]=wfdbloadlib;eval(['cd ' config.WFDB_JAVA_HOME filesep 'example'])
ecgpuwave('100s','test');
[tm,signal]=rdsamp('100s');
pwaves=rdann('100s','test',[],[],[],'p');
plot(tm,signal(:,1));hold on;grid on
plot(tm(pwaves),signal(pwaves),'or')
cd(curdir) %CD back to current directory
Since 0.9.5
See also WFDBDESC, PHYSIONETDB, RDANN, WRANN, SQRS, WQRS, WFDBEXEC