<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.2 20190208//EN" "https://jats.nlm.nih.gov/publishing/1.2/JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:mml="http://www.w3.org/1998/Math/MathML" article-type="brief-report" dtd-version="1.2" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">IJSP</journal-id>
<journal-title-group><journal-title>International Journal of Signal Processing</journal-title></journal-title-group>
<publisher><publisher-name>Thomson &amp; Ryberg Publications</publisher-name></publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.62715/ijsp.2026.v1.i1.pp21-25</article-id>
<article-id pub-id-type="publisher-id">e2e33ceed3</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Short Communication</subject></subj-group></article-categories>
<title-group><article-title>An Efficient Wavelet Thresholding Rule for Denoising Wearable Electrocardiogram Recordings</article-title></title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes"><name><surname>Mansour</surname><given-names>Rania</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1">*</xref></contrib>
<contrib contrib-type="author"><name><surname>Yamamoto</surname><given-names>Aditya</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib>
<contrib contrib-type="author"><name><surname>Roy</surname><given-names>Hirut</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib>
<aff id="aff1"><label>1</label>Department of Electrical Electronics, Pannonian Academy of Electrical Electronics, Czech Republic</aff>
<aff id="aff2"><label>2</label>Department of Electrical Electronics, Rhone University Centre for Electrical Electronics, Indonesia</aff>
<aff id="aff3"><label>3</label>Department of Electrical Electronics, Han River Academy of Electrical Electronics, Spain</aff>
</contrib-group>
<author-notes><corresp id="cor1">* Correspondence: <email>rania.mansour@pannonian.edu</email></corresp></author-notes>
<pub-date pub-type="epub" iso-8601-date="2026-07-23"><day>23</day><month>07</month><year>2026</year></pub-date>
<volume>1</volume>
<issue>1</issue>
<fpage>21</fpage>
<lpage>25</lpage>
<history>
<date date-type="received" iso-8601-date="2026-03-20"><day>20</day><month>03</month><year>2026</year></date>
<date date-type="accepted" iso-8601-date="2026-05-19"><day>19</day><month>05</month><year>2026</year></date>
</history>
<permissions><copyright-statement>© 2026 The Author(s)</copyright-statement><copyright-year>2026</copyright-year><license license-type="open-access" xlink:href="https://creativecommons.org/licenses/by/4.0/"><license-p>This is an open access article distributed under the Creative Commons Attribution 4.0 International License (CC BY 4.0).</license-p></license></permissions>
<self-uri xlink:href="https://trpub.org/article.php?id=e2e33ceed3"/>
<abstract><p>Electrocardiograms recorded by patches, chest straps and wrist devices are corrupted by muscle activity, electrode motion and baseline drift to a degree seldom met at the bedside, and they must be cleaned on processors that run for days on a small battery. Wavelet shrinkage is a natural candidate for the task, yet the textbook rule, one universal threshold applied to every detail coefficient, was derived for white Gaussian noise and tends to flatten the QRS complex when applied to ambulatory recordings. This technical note describes a thresholding rule designed for the wearable setting. It combines four ingredients that are each well established: an integer lifting implementation of the discrete wavelet transform, a noise scale estimated separately at every decomposition level from the median absolute deviation, a level-dependent threshold that is relaxed in a short window around detected R peaks, and the non-negative garrote shrinkage function in place of hard or soft thresholding. The note explains the reasoning behind each choice and the arithmetic cost of the whole. The rule is offered as a design for others to test. No recordings were processed for this note, and the closing section sets out how a fair assessment could be made with public databases and standard noise records.</p></abstract>
<kwd-group kwd-group-type="author">
<kwd>Electrocardiogram</kwd>
<kwd>Wavelet shrinkage</kwd>
<kwd>Level-dependent threshold</kwd>
<kwd>Wearable sensors</kwd>
<kwd>Lifting scheme</kwd>
<kwd>Motion artefact</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="sec1"><title>1. Introduction</title>
<p>The diagnostic content of an electrocardiogram lies in a few brief waveforms per heartbeat: the P wave, the QRS complex and the T wave. In a recording from a wearable device these sit among disturbances that overlap them in both time and frequency. Muscle noise is broadband. Electrode motion produces slow excursions and abrupt steps, and respiration and posture move the baseline. Linear filters with fixed passbands remove what lies outside the cardiac band and can do nothing about what lies inside it.</p>
<p>Wavelet shrinkage acts differently. The signal is expanded in a wavelet basis, small coefficients are suppressed and the remainder is transformed back. Because a QRS complex is represented by a few large coefficients while noise spreads thinly over all of them, thresholding can separate the two even where their spectra overlap. Donoho and Johnstone (1994) proved that a single threshold proportional to the noise standard deviation and to the square root of twice the logarithm of the sample count comes within a logarithmic factor of ideal performance, and Donoho (1995) analysed the soft-thresholding form. The result assumes white Gaussian noise of known level. Wearable recordings violate that assumption, and this note describes a rule adapted to them.</p>
</sec>
<sec id="sec2"><title>2. Why the Universal Threshold Falls Short</title>
<p>Three properties of ambulatory noise defeat the textbook rule. The noise is coloured: motion artefact concentrates at low frequencies and muscle noise at high ones, so its variance differs from one decomposition level to the next and a threshold fixed from the finest level is too small for the coarse levels. It is also non-stationary, rising when the wearer moves and falling at rest. And the universal threshold is conservative by design, which with soft thresholding shrinks every retained coefficient by the full threshold and so lowers the amplitude of the R wave. Since clinical measurements and beat detectors both depend on QRS amplitude and slope, that bias matters more here than in generic denoising.</p>
<p>The statistical literature supplies the remedies. Donoho and Johnstone (1995) showed that thresholds chosen level by level adapt to unknown smoothness, and the practice of estimating the noise scale from the median absolute deviation of the coefficients, divided by 0.6745, makes the estimate insensitive to the few large coefficients that carry signal (Donoho and Johnstone, 1994). Gao (1998) proposed the non-negative garrote as a shrinkage function that is continuous like soft thresholding yet approaches the identity for large coefficients like hard thresholding. The rule below assembles these with one element specific to the electrocardiogram.</p>
</sec>
<sec id="sec3"><title>3. The Proposed Rule</title>
<p>The rule has four parts and a closing step for the baseline, summarised in Table 1. Figure 1 shows how the treatment differs across the sub-bands of the decomposition.</p>
<table-wrap id="tbl1"><label>Table 1</label><caption><p>The parts of the proposed thresholding rule, with the purpose and the per-sample cost of each</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Part</th><th>What is done</th><th>Purpose</th><th>Cost</th></tr></thead><tbody><tr><td>Transform</td><td>Decimated wavelet transform by integer lifting steps</td><td>Multiresolution analysis without floating point</td><td>A few additions and shifts</td></tr><tr><td>Noise scale</td><td>Median absolute deviation per level, per block</td><td>Follow coloured, time-varying noise</td><td>One small sort per level</td></tr><tr><td>Threshold</td><td>Level-wise universal threshold, relaxed around R peaks</td><td>Protect QRS amplitude and slope</td><td>One comparison</td></tr><tr><td>Shrinkage</td><td>Non-negative garrote function</td><td>Continuity with little bias on large coefficients</td><td>One division</td></tr><tr><td>Approximation</td><td>Attenuate the final approximation band</td><td>Remove baseline drift</td><td>One multiplication</td></tr></tbody></table></table-wrap>
<fig id="fig1"><label>Figure 1</label><caption><p>Treatment applied to each group of sub-bands under the proposed rule</p></caption><graphic xlink:href="IJSP_V1I1_2026_04_Figure1.png" mimetype="image" mime-subtype="png"/></fig>
<sec id="sec3-2"><title>3.1. Transform</title>
<p>The discrete wavelet transform is computed by the pyramid algorithm of Mallat (1989), whose cost is proportional to the number of samples. On a microcontroller it is best realised by lifting. Daubechies and Sweldens (1998) showed that any finite wavelet filter pair factors into a sequence of simple prediction and update steps performed in place, and Calderbank et al. (1998) showed that rounding inside those steps yields transforms that map integers to integers and remain exactly invertible. A short biorthogonal wavelet with a smooth, symmetric analysis function is recommended, since symmetry avoids phase distortion of the waveform. The depth is chosen from the sampling rate so that the final approximation band lies below the lowest frequencies of clinical interest.</p>
</sec>
<sec id="sec3-3"><title>3.2. Noise scale and threshold</title>
<p>The recording is processed in blocks of a few seconds. Within each block and at each level j, the noise scale σj is the median absolute deviation of the detail coefficients divided by 0.6745, and the base threshold is σj multiplied by the square root of 2 ln nj, where nj is the number of coefficients at that level in the block. Estimating the scale afresh for each block lets the threshold rise during movement and fall at rest.</p>
<p>The element specific to the electrocardiogram is a relaxation of the threshold where the QRS complex is expected. Most QRS energy lies between roughly 5 and 15 Hz, the band on which the detector of Pan and Tompkins (1985) is built. A simple detector of that kind, run on the detail levels covering this band, marks the R peaks. In a short window around each mark, and only at the QRS-band levels, the threshold is multiplied by a factor below one. Elsewhere the full threshold applies. The coarser detail levels, where P and T waves reside, receive a lighter threshold throughout because those waves are low in amplitude and easily erased. The relaxation factor and window length are the two free parameters of the rule.</p>
</sec>
<sec id="sec3-4"><title>3.3. Shrinkage and baseline</title>
<p>Each detail coefficient d with threshold t is replaced by zero if its magnitude is below t and by d - t2/d otherwise, which is the garrote function of Gao (1998). The subtracted term decays as the coefficient grows, so the tall R-wave coefficients pass almost unaltered while coefficients near the threshold are shrunk smoothly, without the discontinuity that makes hard thresholding produce spurious oscillations. Baseline drift is handled at no extra cost by attenuating the final approximation coefficients before reconstruction, the wavelet equivalent of a high-pass filter.</p>
</sec>
</sec>
<sec id="sec4"><title>4. Why the Rule Is Cheap</title>
<p>Efficiency follows from the structure and can be stated without measurement. The lifting transform and its inverse need a fixed, small number of integer additions and shifts per sample. Thresholding adds one comparison and, for retained coefficients only, one division. The median is the single operation that is not linear in block length when done by sorting, and it is needed once per level per block on arrays that halve in size at each level; a selection algorithm or a running histogram of the integer coefficients brings it to linear time as well. Memory is one block of samples, since lifting works in place. The decimated transform is not translation invariant, and the cycle-spinning remedy of Coifman and Donoho (1995) would multiply the cost by the number of shifts averaged. The rule leaves it out and relies on the garrote function to limit the artefacts that cycle spinning was designed to remove.</p>
</sec>
<sec id="sec5"><title>5. How the Rule Could Be Assessed</title>
<p>No recordings were processed for this note, and the rule should be regarded as untested. A fair assessment is nonetheless easy to specify, because the field has shared materials. The MIT-BIH Arrhythmia Database provides annotated two-channel recordings (Moody and Mark, 2001), and the noise stress test of Moody et al. (1984) supplies separate records of baseline wander, muscle artefact and electrode motion artefact that can be added to clean segments at chosen signal-to-noise ratios. With a known clean reference, the improvement in signal-to-noise ratio can be computed directly.</p>
<p>Waveform fidelity is not the only criterion. What the user of a wearable needs is that beats are still detected and that intervals and amplitudes are still measured correctly. An assessment should therefore also report the sensitivity and positive predictivity of a standard QRS detector before and after denoising, and the shift in wave boundaries found by a delineator such as that of Martínez et al. (2004). The rule should be compared with the universal soft threshold, with a level-dependent threshold lacking the R-peak relaxation, and with a conventional band-pass filter, so that the contribution of each part can be seen. Records added from noise databases remain an approximation to the artefacts of a specific device, and final judgement needs recordings from that device.</p>
</sec>
<sec id="sec6"><title>6. Conclusion</title>
<p>A thresholding rule suited to wearable electrocardiography can be put together from components that are individually well understood: integer lifting, level-wise robust noise estimation, a threshold relaxed around the R peak and garrote shrinkage. The combination respects the coloured, changing noise of ambulatory recording, guards the QRS complex against the amplitude loss of the universal rule and stays within the arithmetic budget of a small processor. Its two free parameters and its dependence on a preliminary R-peak detector are the points at which it may fail, and they are the first things an experimental study should examine.</p>
</sec>
<sec id="sec7"><title>Conflict of interest</title>
<p>To be declared by the author(s) before publication.</p>
</sec>
<sec id="sec8"><title>Funding</title>
<p>To be declared by the author(s) before publication.</p>
</sec>
</body>
<back>
<ref-list><title>References</title>
<ref id="ref1"><label>1</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Calderbank, A.R., Daubechies, I., Sweldens, W. and Yeo, B.-L.</string-name></person-group> (<year>1998</year>). <article-title>Wavelet transforms that map integers to integers</article-title>. <source>Applied and Computational Harmonic Analysis</source>, <volume>5(3)</volume>, <fpage>332-369</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1006/acha.1997.0238">https://doi.org/10.1006/acha.1997.0238</ext-link></mixed-citation></ref>
<ref id="ref2"><label>2</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Coifman, R.R. and Donoho, D.L.</string-name></person-group> (<year>1995</year>). <article-title>Translation-invariant de-noising</article-title>. <source>Wavelets and Statistics (Lecture Notes in Statistics, vol. 103), Springer, New York</source>, <fpage>125-150</fpage></mixed-citation></ref>
<ref id="ref3"><label>3</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Daubechies, I. and Sweldens, W.</string-name></person-group> (<year>1998</year>). <article-title>Factoring wavelet transforms into lifting steps</article-title>. <source>Journal of Fourier Analysis and Applications</source>, <volume>4(3)</volume>, <fpage>247-269</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1007/BF02476026">https://doi.org/10.1007/BF02476026</ext-link></mixed-citation></ref>
<ref id="ref4"><label>4</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Donoho, D.L.</string-name></person-group> (<year>1995</year>). <article-title>De-noising by soft-thresholding</article-title>. <source>IEEE Transactions on Information Theory</source>, <volume>41(3)</volume>, <fpage>613-627</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/18.382009">https://doi.org/10.1109/18.382009</ext-link></mixed-citation></ref>
<ref id="ref5"><label>5</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Donoho, D.L. and Johnstone, I.M.</string-name></person-group> (<year>1994</year>). <article-title>Ideal spatial adaptation by wavelet shrinkage</article-title>. <source>Biometrika</source>, <volume>81(3)</volume>, <fpage>425-455</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1093/biomet/81.3.425">https://doi.org/10.1093/biomet/81.3.425</ext-link></mixed-citation></ref>
<ref id="ref6"><label>6</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Donoho, D.L. and Johnstone, I.M.</string-name></person-group> (<year>1995</year>). <article-title>Adapting to unknown smoothness via wavelet shrinkage</article-title>. <source>Journal of the American Statistical Association</source>, <volume>90(432)</volume>, <fpage>1200-1224</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1080/01621459.1995.10476626">https://doi.org/10.1080/01621459.1995.10476626</ext-link></mixed-citation></ref>
<ref id="ref7"><label>7</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Gao, H.-Y.</string-name></person-group> (<year>1998</year>). <article-title>Wavelet shrinkage denoising using the non-negative garrote</article-title>. <source>Journal of Computational and Graphical Statistics</source>, <volume>7(4)</volume>, <fpage>469-488</fpage></mixed-citation></ref>
<ref id="ref8"><label>8</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Mallat, S.G.</string-name></person-group> (<year>1989</year>). <article-title>A theory for multiresolution signal decomposition: The wavelet representation</article-title>. <source>IEEE Transactions on Pattern Analysis and Machine Intelligence</source>, <volume>11(7)</volume>, <fpage>674-693</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/34.192463">https://doi.org/10.1109/34.192463</ext-link></mixed-citation></ref>
<ref id="ref9"><label>9</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Martínez, J.P., Almeida, R., Olmos, S., Rocha, A.P. and Laguna, P.</string-name></person-group> (<year>2004</year>). <article-title>A wavelet-based ECG delineator: Evaluation on standard databases</article-title>. <source>IEEE Transactions on Biomedical Engineering</source>, <volume>51(4)</volume>, <fpage>570-581</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/TBME.2003.821031">https://doi.org/10.1109/TBME.2003.821031</ext-link></mixed-citation></ref>
<ref id="ref10"><label>10</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Moody, G.B. and Mark, R.G.</string-name></person-group> (<year>2001</year>). <article-title>The impact of the MIT-BIH Arrhythmia Database</article-title>. <source>IEEE Engineering in Medicine and Biology Magazine</source>, <volume>20(3)</volume>, <fpage>45-50</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/51.932724">https://doi.org/10.1109/51.932724</ext-link></mixed-citation></ref>
<ref id="ref11"><label>11</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Moody, G.B., Muldrow, W.K. and Mark, R.G.</string-name></person-group> (<year>1984</year>). <article-title>A noise stress test for arrhythmia detectors</article-title>. <source>Computers in Cardiology</source>, <volume>11</volume>, <fpage>381-384</fpage></mixed-citation></ref>
<ref id="ref12"><label>12</label><mixed-citation publication-type="other"><person-group person-group-type="author"><string-name>Pan, J. and Tompkins, W.J.</string-name></person-group> (<year>1985</year>). <article-title>A real-time QRS detection algorithm</article-title>. <source>IEEE Transactions on Biomedical Engineering</source>, <volume>BME-32(3)</volume>, <fpage>230-236</fpage> <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1109/TBME.1985.325532">https://doi.org/10.1109/TBME.1985.325532</ext-link></mixed-citation></ref>
</ref-list>
</back>
</article>