Indonesian
J
our
nal
of
Electrical
Engineering
and
Computer
Science
V
ol.
36,
No.
2,
No
v
ember
2024,
pp.
837
∼
845
ISSN:
2502-4752,
DOI:
10.11591/ijeecs.v36.i2.pp837-845
❒
837
Electr
ocardiogram
r
econstruction
based
on
Hermite
inter
polating
polynomial
with
Chebyshe
v
nodes
Shashwati
Ray
1
,
V
andana
Chouhan
2
1
Department
of
Electrical
Engineering,
Bhilai
Institute
of
T
echnology
,
Dur
g
(C.G.),
India
2
Department
of
Electronics
and
T
elecommunication
Engineering,
Bhilai
Institute
of
T
echnology
,
Dur
g
(C.G.),
India
Article
Inf
o
Article
history:
Recei
v
ed
No
v
6,
2023
Re
vised
Jul
5,
2024
Accepted
Jul
14,
2024
K
eyw
ords:
Chebyshe
v
nodes
Electrocardiogram
Equally
spaced
nodes
Hermite
Noise
Polynomial
ABSTRA
CT
Electrocardiogram
(ECG)
signals
generate
massi
v
e
v
olume
of
digital
data,
so
the
y
need
to
be
suitably
compressed
for
ef
cient
transmission
and
storage.
Poly-
nomial
approximations
and
polynomial
interpolation
ha
v
e
been
used
for
ECG
data
compression
where
the
data
signal
is
described
by
polynomial
coef
cients
only
.
Here,
we
propose
approximation
using
hermite
polynomial
interpolation
with
chebyshe
v
nodes
for
compressing
ECG
signals
that
consequently
denoises
them
too.
Recommended
algorithm
is
applied
on
v
arious
ECG
signals
tak
en
from
MIT
-BIH
arrh
ythmia
database
without
an
y
additional
noise
as
the
signals
are
already
contaminated
with
noise.
Performance
of
the
proposed
algorithm
is
e
v
aluated
using
v
arious
performance
metrics
and
compared
with
some
recent
compression
techniques.
Experimental
results
pro
v
e
that
the
proposed
method
ef
ciently
compresses
the
ECG
signals
while
preserving
the
minute
details
of
important
morphological
features
of
ECG
signal
required
for
clinical
diagnosis.
This
is
an
open
access
article
under
the
CC
BY
-SA
license
.
Corresponding
A
uthor:
Shashw
ati
Ray
Department
of
Electrical
Engineering,
Bhilai
Institute
of
T
echnology
Dur
g
(C.G.),
India
Email:
shashw
atiray@yahoo.com
1.
INTR
ODUCTION
Electrocardiogram
(ECG
or
EKG)
is
a
recording
of
1-D
time
series
data
sequence
generated
by
cardiac
muscles.
The
y
help
to
track
and
detect
abnormalities
in
the
heart
rh
ythm
based
on
the
morphology
and
frequenc
y
of
heartbeat
[1].
24
hours
ECG
record
with
the
sampling
rate
of
360
Hz
and
11
bit/sample
data
resolution
requires
about
43
MBytes
per
channel
[2].
Therefore,
an
ef
fecti
v
e
data
compression
scheme
is
often
required
for
ef
cient
ECG
data
storage
and
transmission
o
v
er
telephone
line
or
digital
telecommunication
netw
ork
[3].
ECG
compression
techniques
can
be
broadly
classied
i
nto
three
major
cate
gories:
direct
time
domain,
parameter
e
xtraction
and
transform
domain
method
[4].
In
direct
ti
me
domain
method
compression
is
achie
v
ed
by
nding
correlation
between
the
adjace
nt
samples,
i.e.,
intra-beat
redundanc
y
in
a
group
and
encode
them
into
a
sma
ller
sub-group
[5].
Some
popu-
lar
algorithms
of
direct
method
are:
amplitude
zone
time
epoch
coding
(AZTECH)
[6],
coordinate
reduction
time
encoding
system
(COR
TES)
[7],
entrop
y
coding
[8],
scan-along
polygonal
approximation
(SAP
A)
[9]
and
long–ter
m
prediction
(L
TP)
[10].
P
arameter
e
xtraction
methods
e
xtract
important
morphological
features
from
the
signal
and
encode
these
features
to
achie
v
e
desired
compression
[11],
e.g.,
linear
prediction
[12]
and
residual
encoding
methods
[13].
In
transform
domain
method
signals
are
transformed
to
another
space
using
v
arious
transform
schemes
[14],
in
the
subsequent
steps
small
transform
coef
cients
are
discarded
J
ournal
homepage:
http://ijeecs.iaescor
e
.com
Evaluation Warning : The document was created with Spire.PDF for Python.
838
❒
ISSN:
2502-4752
and
only
critical
information
is
encoded
to
achie
v
e
the
desired
compression.
T
ransform
techniques
include
fourier
transforms
[15],
discrete
w
a
v
elet
transform
(D
WT)
[16],
discrete
cosine
transform
(DCT)
[17],
discrete
Le
gendre
transform
[18]
and
Karhunen
Louv
e
transform
(KL
T)
[19].
Polynomials
nd
application
in
signal
processing,
viz.,
for
ltering
noisy
signals,
interpolating
data
se-
quence,
and
for
data
compression
[20].
High
compression
ratio
can
be
achie
v
ed
by
polynomial
approximation,
since,
polynomial
coef
cients
describing
the
signal
are
only
required
for
compression
[21].
Le
gendre
poly-
nomials
is
used
by
[22]
for
image
compression
and
reconstruction.
F
or
higher
order
of
polynomials
the
V
an-
dermonde
matrix
gets
ill
conditioned
which
is
a
v
oided
by
di
viding
image
matrix
into
sub-matrices
to
achie
v
e
desired
compression
ratio.
ECG
data
compression
utilizing
Jacobi
polynomials
is
proposed
by
[23].
ECG
signals
are
rst
se
gmented
into
blocks
that
match
with
cardiac
c
ycles
before
being
decomposed
in
Jacobi
poly-
nomials
bases.
Gauss
quadratures
mechanism
for
numerical
inte
gration
is
used
to
compute
Jacobi
transforms
coef
cients.
T
o
achie
v
e
desired
compression,
coef
cients
of
small
v
alues
are
discarded
in
the
reconstruction
stage.
As
the
deri
v
ed
polynomials
use
recurrence
formula,
rounding
errors
pile
up
during
the
computation
and
polynomials
gradually
denature
with
increase
in
order
of
the
polynomial.
ECG
data
compression
based
on
B-spline
basis
functions
is
proposed
by
[24].
The
position
of
the
knots
are
computed
using
run-length
coding.
Ho
we
v
er
,
for
f
alse
R-R
comple
x
ne
w
sequence
of
knots
has
to
be
sent,
thereby
increasing
the
o
v
erhead
data
resulting
in
increasing
computational
comple
xity
.
Nyg
aard
and
Haugland
[25]
proposed
piece
wise
polynomial
approximation
for
reconstructing
the
ECG
signal
by
second
order
polynomials.
Khetk
eeree
and
Chansamorn
[26]
introduced
signal
reconstruction
based
on
the
second
order
tetration
polynomial.
Fundamental
signals
such
as
square,
sa
w-tooth
and
sine
w
a
v
e
with
v
arious
sampling
resolution
were
applied
to
test
the
interpolation
performance.
High
v
alues
of
peak
signal-to-noise
ratio
(PSNR)
were
obtai
ned
for
square
w
a
v
e
and
sa
w-tooth
w
a
v
e.
No
w
a
days,
long
term
ECG
monitoring
is
applied
for
management
of
cardio
v
ascular
diseases
where
wireless
technology
is
used
to
transmit
ECG
data
through
communication
channels.
Channel
bandwidth
can
be
optimized
by
performing
ECG
compression.
ECG
signals
can
be
compressed
using
polynomial
interpolation.
Here
high
compression
ratio
can
be
achie
v
ed,
since
ECG
signals
can
be
reconstructed
using
fe
w
sampling
points.
W
e
propose
here
an
algorithm
to
obtain
an
ECG
approximating
model
based
on
lagrange
form
of
hermite
interpolating
polynomial
with
chebyshe
v
nodes.
Hermite
interpolating
polynomials
are
more
rob
ust,
since
we
ha
v
e
higher
de
gree
of
freedom
with
deri
v
ati
v
e
v
alues
as
additional
information
which
is
equi
v
alent
to
almost
twice
the
order
of
the
interpolating
polynomial.
This
polynomial
smoothly
interpolates
between
the
k
e
y-points
and
compresses
the
ECG
data,
thus
f
acilitating
less
data
for
storage
and
transmission.
In
the
process
it
also
denoises
the
ECG
signal
in
an
ef
cient
w
ay
while
preserving
the
morphologi
cal
features
as
required
by
the
cardiologist.
The
or
g
anization
of
the
paper
is
as
follo
ws:
in
section
2
we
e
xplain
the
research
methodology
of
this
w
ork,
in
section
3
we
pro
vide
e
xperimental
details
and
in
section
4
we
present
the
conclusions
dra
wn
out
of
our
w
ork.
2.
RESEARCH
METHOD
In
real
life
applications
e
xperimental
data
are
in
the
form
of
set
of
discrete
data
points
and
functional
relation
between
input
and
output
is
nondeterministic.
In
such
situations
polynomial
interpolation
plays
an
important
role
in
determining
a
polynomial
matching
the
points.
In
our
research
w
ork,
we
approximate
an
ECG
signal
f
consisting
of
N
ECG
samples
with
lagrange
form
of
hermite
interpolating
polynomial
H
p
n
(
x
)
using
n
chebyshe
v
nodes.
Our
research
methodology
comprises
mainly
of
v
e
stages:
-
Obtain
ra
w
discrete
ECG
data
as
.mat
le
from
the
MIT
-BIH
arrh
ythmia
database
no
w
freely
a
v
ailable
on
Ph
ysioNet.
-
At
this
sta
g
e
preprocess
and
normalize
the
signals.
The
ECG
data
is
con
v
erted
into
ph
ysical
units
(mV),
then
normalised
by
reducing
the
g
ain
and
re-scaled
to
limit
the
range
within
[
−
1
,
1]
.
-
Map
the
‘
n
’
chebyshe
v
nodes,
x
k
,
k
=
1
,
...,
n
on
the
abscissa
of
time
(seconds)
with
the
equi
v
alent
ECG
data
to
obtain
them
as
function
v
alues
f
x
k
,
k
=
1
,
...,
n
.
-
Obtain
the
deri
v
ati
v
es
f
′
x
k
at
all
k
=
1
,
...,
n
points
using
numerical
dif
ferentiation
method.
-
Construct
the
lagrange
form
of
hermite
polynomial
with
f
x
k
and
f
′
x
k
,
k
=
1
,
...,
n
.
Indonesian
J
Elec
Eng
&
Comp
Sci,
V
ol.
36,
No.
2,
No
v
ember
2024:
837–845
Evaluation Warning : The document was created with Spire.PDF for Python.
Indonesian
J
Elec
Eng
&
Comp
Sci
ISSN:
2502-4752
❒
839
2.1.
Function
v
alues
at
inter
polating
nodes
Since
the
ECG
samples
are
obtained
for
an
arbitrary
length,
we
require
to
transform
them
to
the
designated
interv
al
and
compute
the
function
v
alues
(we
consider
the
ECG
signal
as
a
function)
at
all
the
interpolating
nodes.
Let
the
ECG
signal
be
sampled
at
a
frequenc
y
F
s
and
the
sampled
v
alues
be
dened
as
function
v
alues
f
,
thus
comprising
of
total
N
samples.
W
ith
spacing
h
=
1
/F
s
,
compute
the
end
points
of
[
a,
b
]
on
the
abscissa
of
time
as,
a
=
1
/F
s
,
b
=
N
/F
s
compute
the
n
chebyshe
v
nodes
x
k
,
k
=
1
,
...,
n
on
[
a,
b
]
as,
x
k
=
a
+
b
2
−
b
−
a
2
cos
2
k
−
1
2
n
π
,
k
=
1
,
·
·
·
,
n
nd
all
the
equi
v
alent
function
v
alues
as,
f
x
k
=
f
(
x
k
[1
:
n
])
an
y
missing
function
v
alue
is
e
v
aluated
with
linear
interpolation
using
the
adjacent
sampled
v
alues.
No
w
,
we
ha
v
e
the
data
points
in
the
form
(
x
k
,
f
x
k
)
,
k
=
1
,
...,
n
.
2.2.
Lagrange
f
orm
of
hermite
inter
polation
Hermite
interpolating
polynomials
require
the
kno
wledge
of
the
deri
v
ati
v
es
at
the
interpolating
n
odes
.
Since
the
function
v
alues
are
discrete,
the
deri
v
ati
v
es
are
computed
applying
numerical
dif
ferentiation
meth-
ods.
The
numerical
deri
v
ati
v
es
are
computed
using
forw
ard,
central
and
bac
k
w
a
rd
dif
ferences.
Use
forw
ard
dif
ference
to
compute
f
′
x
k
at
lo
wer
points
of
[
a,
b
]
,
f
′
x
k
=
1
2
h
[
−
3
f
x
k
+
4
f
x
k
+1
−
f
x
k
+2
]
use
backw
ard
dif
ference
to
compute
f
′
x
k
at
upper
points
of
[
a,
b
]
,
f
′
x
k
=
1
2
h
[
−
3
f
x
k
−
4
f
x
k
−
1
+
f
x
k
−
2
]
use
central
dif
ference
to
compute
f
′
x
k
at
intermediate
points
of
[
a,
b
]
.
f
′
x
k
=
1
2
h
[
f
x
k
+1
−
f
x
k
−
1
]
F
or
the
data
of
the
form
(
x
k
,
f
x
k
)
,
(
x
k
,
f
′
x
k
)
,
k
=
1
,
...,
n
,
the
unique
lagrange
form
of
hermite
polynomial
H
p
n
(
x
)
of
de
gree
2
n
+
1
that
agrees
with
f
x
k
and
f
′
x
k
is
gi
v
en
by:
H
p
n
(
x
)
=
n
X
k
=1
f
x
k
A
n,k
(
x
)
+
n
X
k
=1
f
′
x
k
B
n,k
(
x
)
where,
A
n,k
(
x
)
=
[1
−
2(
x
−
x
k
)
L
′
n,k
(
x
k
)]
L
2
n,k
(
x
)
and
B
n,k
(
x
)
=
(
x
−
x
k
)
L
2
n,k
(
x
)
where,
L
n,k
(
x
)
denotes
lagrange
basis
function
of
order
n
dened
by
,
L
n,k
(
x
)
=
n
Y
i
=1
,i
̸
=
k
(
x
−
x
i
)
(
x
k
−
x
i
)
the
error
using
lagrange
form
of
hermite
interpolation
with
chebyshe
v
nodes
is
gi
v
en
by
,
E
(
x
)
=
|
f
(
x
)
−
H
p
n
(
x
)
|
≤
1
2
n
(
n
+
1)!
(
b
−
a
)
2
(
n
+1)
max
a
≤
ξ
≤
b
|
f
(
n
+1)
(
ξ
)
|
if
E
≥
ε
,
where
the
tolerance
ε
=
10
−
2
,
then
n
is
increased
by
10
and
the
entire
procedure
is
repeated.
Electr
ocar
dio
gr
am
r
econstruction
based
on
Hermite
interpolating
...
(Shashwati
Ray)
Evaluation Warning : The document was created with Spire.PDF for Python.
840
❒
ISSN:
2502-4752
3.
RESUL
TS
AND
AN
AL
YSIS
The
proposed
algorithms
are
implemented
in
MA
TLAB
(R2013b)
v
ersion.
Each
part
of
the
proposed
ECG
approximation
algorithm
is
written
in
the
.m
le
as
a
subroutine
module.
All
the
computations
are
carried
out
on
ECG
signals
tak
en
from
MIT
-BIH
arrh
ythmia
database
a
v
ailable
on
Ph
ysioNet
[27].
W
e
consider
here
v
arious
signals
of
channel
1,
sampled
at
360
Hz
with
a
resolution
of
11
bits
per
sample
with
duration
of
5
seconds
resulting
in
1,800
s
amples.
These
sample
points
are
the
ECG
signal
magnitudes
obtained
at
equal
interv
als
of
‘
1
/
360
’
second.
W
e
perform
the
delity
assessment
of
the
propos
ed
approximation
method
using
the
performance
or
error
measures
as
-
root
m
ean
square
error
(RMS),
percentage
root
mean
dif
ference
(PRD),
signal
to
noise
ratio
(SNR),
and
compression
ratio
(CR).
Here,
we
consider
CR
as
the
ratio
of
the
number
of
bytes
in
the
uncompressed
representation
to
the
number
of
bytes
in
the
compressed
representation.
T
o
test
the
ef
cac
y
of
our
proposed
method
we
apply
the
de
v
eloped
algorithms
on
12
records
sho
wn
in
Figures
1
and
2
with
duration
of
5
seconds
and
approximate
them
in
the
form
of
respecti
v
e
polynomials
with
300
chebyshe
v
nodes.
All
the
12
obtained
results
are
illustrated
in
Figures
1
and
2,
Figure
1(a)
#100,
Figure
1(b)
#104,
Figure
1(c)
#108,
Figure
1(d)
#112,
Figure
1(
e)
#115,
Figure
1(f)
#117,
and
Figure
2(a)
#122,
Figure
2(b)
#201,
Figure
2(c)
#205,
Figure
2(d)
#207,
Figure
2(e)
#214,
Figure
2(f)
#220
depicting
the
original
signals
and
the
reconstructed
polynomials
in
red
and
blue
colours
respecti
v
ely
.
(a)
(b)
(c)
(d)
(e)
(f)
Figure
1.
Original
noisy
ECG
signals
(pink)
and
reconstructed
samples
by
proposed
method
(blue):
(a)
ECG
record
#
100
,
(b)
ECG
record
#
104
,
(c)
ECG
record
#
108
,
(d)
ECG
record
#
112
,
(e)
ECG
record
#
115
,
and
(f)
ECG
record
#
117
Indonesian
J
Elec
Eng
&
Comp
Sci,
V
ol.
36,
No.
2,
No
v
ember
2024:
837–845
Evaluation Warning : The document was created with Spire.PDF for Python.
Indonesian
J
Elec
Eng
&
Comp
Sci
ISSN:
2502-4752
❒
841
(a)
(b)
(c)
(d)
(e)
(f)
Figure
2.
Original
noisy
ECG
signals
(pink)
and
reconstructed
samples
by
proposed
method
(blue):
(a)
ECG
record
#
122
,
(b)
ECG
record
#
201
,
(c)
ECG
record
#
205
,
(d)
ECG
record
#
207
,
(e)
ECG
record
#
214
,
and
(f)
ECG
record
#
220
Since
the
proposed
method
is
an
e
xtension
of
the
method
propose
d
by
Y
ada
v
and
Ray
[21],
we
deem
it
t
to
compare
the
performance
statistics
of
the
proposed
method
with
the
latter
.
Y
ada
v
and
Ray
[21],
ha
v
e
approximated
all
these
records
of
same
duration
with
assorted
noise
le
v
els
using
lagrange
interpolating
poly-
nomial
with
400
chebyshe
v
nodes
respecti
v
ely
.
The
performance
statistics
of
the
method
by
Y
ada
v
and
Ray
,
and
the
proposed
method
are
reported
in
T
able
1.
It
is
w
orth
mentioning
that
in
the
process
of
approximating,
both
the
methods
e
v
entually
denoise
the
ECG
signals.
Comparing
the
respecti
v
e
entries
of
T
able
1
for
each
ECG
record,
we
observ
e
that
lo
wer
order
hermit
e
form
of
interpolati
n
g
polynomial
outperforms
the
lagrange
form
in
all
performance
metrics.
Lo
wer
v
alues
of
RMS
are
indicati
v
e
of
least
distortion
and
better
approximation.
The
v
alues
of
CR
in
both
the
methods
remain
constant
for
all
the
signals,
because
the
number
of
samples
and
the
respecti
v
e
number
of
interpolating
nodes
are
persistent
in
all
the
signals.
Electr
ocar
dio
gr
am
r
econstruction
based
on
Hermite
interpolating
...
(Shashwati
Ray)
Evaluation Warning : The document was created with Spire.PDF for Python.
842
❒
ISSN:
2502-4752
The
tw
o
important
features
of
a
compression
algorithm
are
the
compression
measure
and
the
recon-
struction
error
.
Not
man
y
approximating
methods
are
a
v
ailable
in
the
e
xisting
literature
for
ECG
signals.
T
o
ha
v
e
a
comprehensi
v
e
re
vie
w
of
the
proposed
method,
we
abstractly
compare
the
performance
of
the
proposed
method
with
tw
o
e
xisting
recent
w
orks
on
ECG
compression.
The
identied
methods
are:
Y
ang
et
al.
[14]
us-
ing
empirical
mode
decomposition
(EMD)
and
Hamza
et
al.
[28]
based
on
discrete
w
a
v
elet
transform
(D
WT)
and
dual
encoding
technique.
T
able
1.
Comparison
of
the
proposed
method
with
Y
ada
v
and
Ray
[21]
for
signal
length
of
5
sec
Record
Y
ada
v
and
Ray
[21]
with
n
=
400
Proposed
method
with
n
=
300
RMS
PRD
SNR
CR
RMS
PRD
SNR
CR
100
0.11
15.75
8.45
4.49
0.04
12.27
11.44
6.00
104
0.06
18.92
12.34
4.49
0.07
19.85
11.86
6.00
108
0.02
6.39
16.76
4.49
0.03
7.81
15.33
6.00
112
0.02
2.55
17.09
4.49
0.02
1.70
21.35
6.00
115
0.09
14.72
10.57
4.49
0.04
6.21
18.07
6.00
117
0.03
3.71
16.81
4.49
0.02
2.73
19.71
6.00
122
0.04
4.99
17.85
4.49
0.03
3.53
21.21
6.00
201
0.02
9.22
17.29
4.49
0.02
8.14
19.68
6.00
205
0.05
11.43
10.51
4.49
0.02
4.88
18.33
6.00
207
0.02
5.10
23.56
4.49
0.02
6.97
21.93
6.00
214
0.04
8.79
19.57
4.49
0.03
7.67
21.96
6.00
220
0.11
15.75
8.70
4.49
0.06
8.60
14.07
6.00
F
or
comparison
with
Y
ang
et
al.
[14]
method,
we
choose
8
MIT
-BIH
arrh
ythmia
data
sets
[27]
as
test
signals
with
time
period
as
4.2
seconds
and
sampling
rate
as
360
Hz
for
all
the
signals.
The
8
records
are
referred
as
#
100
,
#103,
#
107
,
#
109
,
#
116
,
#
117
,
#
119
,
and
#
200
.
T
o
e
v
aluate
the
quality
of
the
pro-
posed
algorithm
we
use
RMS
and
CR
as
the
perf
ormance
measures
and
illustrate
the
results
as
bar
graphs
in
Figure
3
and
Figure
4.
From
the
comparisons
we
can
easily
infer
that
the
proposed
method
f
airs
v
ery
well
in
RMS
metric
and
compares
well
in
CR
metric.
F
or
comparison
with
Hamza
et
al.
[28]
method,
we
choos
e
5
MIT
-BIH
arrh
ythmia
data
sets
as
test
signals
with
time
period
as
10
seconds
and
sampling
rate
as
360
Hz
for
all
the
signals
resulting
in
3,600
samples.
The
5
records
are
referred
as
#
100
,
#
109
,
#
115
,
#
119
,
and
#
200
.
T
o
e
v
aluate
the
quality
of
the
proposed
algorithm
we
use
RMS
as
the
performance
measure
and
depict
the
results
in
Figure
5
fr
om
where
we
observ
e
that
the
proposed
method
performs
v
ery
well.
Here,
we
ha
v
en’
t
considered
the
comparison
of
CR
since
we
dif
fer
in
our
denitions.
Figure
3.
Bar
graph
of
RMS
v
alues
of
v
arious
records
with
signal
length
of
4
.
2
sec
obtained
by
Y
ang
et
al.
[14]
and
the
proposed
method
Indonesian
J
Elec
Eng
&
Comp
Sci,
V
ol.
36,
No.
2,
No
v
ember
2024:
837–845
Evaluation Warning : The document was created with Spire.PDF for Python.
Indonesian
J
Elec
Eng
&
Comp
Sci
ISSN:
2502-4752
❒
843
Figure
4.
Bar
graph
of
CR
v
alues
of
v
arious
records
with
signal
length
of
4
.
2
sec
obtained
by
Y
ang
et
al.
[14]
and
the
proposed
method
Figure
5.
Bar
graph
of
RMS
v
alues
of
v
arious
records
with
signal
length
of
10
sec
obtained
by
Hamza
et
al.
[28]
and
the
proposed
method
4.
CONCLUSION
In
this
w
ork,
the
superiority
of
the
proposed
approximation
model
of
lagrange
form
of
hermite
poly-
nomial
interpolation
with
chebyshe
v
nodes
is
established
by
applying
on
v
arious
ECG
signals
tak
en
from
MIT/BIH
arrh
ythmia
database
and
comparing
with
fe
w
e
xisti
ng
methods
tak
en
from
recent
literature.
From
all
the
analysis
we
infer
that
the
proposed
method
of
approximation
outperforms
all
the
methods
in
most
of
the
metrics.
The
proposed
method
compresses
ECG
signal
thus
reducing
the
memory
requir
ement.
Apart
from
this,
the
proposed
scheme
not
only
eliminates
noise,
b
ut
also
preserv
es
important
morphological
features
re-
quired
for
analysis
of
v
arious
conditions
lik
e
arrh
ythmias,
inadequate
coronary
artery
blood
o
w
,
electrolyte
disturbances,
and
cardiomyopath
y
.
Most
signicant
is
that
the
proposed
method
is
able
to
con
v
ert
the
ECG
signal
into
a
polynomial;
and
all
polynomial
operations
emphasize
can
be
applied
to
e
xtract
v
arious
morpho-
logical
features
for
the
diagnosis
of
v
arious
diseases
that
are
re
ected
in
the
ECG.
The
proposed
model
can
also
be
e
xtended
in
approximating
other
time
series
data
such
as
economic
and
sales
forecasting,
b
udgetary
and
stock
mark
et
analysis,
yield
projections,
process
and
quality
control,
to
predict
the
future
price
of
the
stock
mark
et,
and
e
xchange
rate
forecast.
Furthermore,
the
proposed
method
is
riddled
with
certain
challenges.
In
case
of
critical
base
line
w
ander
additional
preprocessing
step
has
to
be
applied.
Moreo
v
er
,
detrending
of
time
series
data
is
necessary
whene
v
er
there
is
a
base
line
drift
in
the
signal.
Electr
ocar
dio
gr
am
r
econstruction
based
on
Hermite
interpolating
...
(Shashwati
Ray)
Evaluation Warning : The document was created with Spire.PDF for Python.
844
❒
ISSN:
2502-4752
REFERENCES
[1]
T
.
A.
V
u,
H.
Q.
Huy
,
P
.
D.
Khanh,
N.
T
.
M.
Huyen,
T
.
T
.
T
.
Uyen,
and
P
.
T
.
V
.
Huong,
“Classify
arrh
ythmia
by
using
2D
spectral
images
and
deep
neural
netw
ork,
”
Indonesian
Journal
of
Electrical
Engineering
and
Computer
Science
(IJEECS)
,
v
ol.
25,
no.
2,
pp.
931–940,
2022,
doi:
10.11591/ijeecs.v25.i2.pp931-940.
[2]
M
.
M.
Goudarzi,
A.
T
aheri,
and
M.
Poo
yan,
“Ef
cient
method
for
ECG
compression
using
tw
o
dimensional
multiw
a
v
elet
transform,
”
International
Journal
of
Information
T
echnology
,
v
ol.
1,
no.
3,
pp.
226–232,
2004.
[3]
Z.
Lu,
D.
Y
.
Kim,
and
W
.
A.
Pearlman,
“W
a
v
elet
compression
of
ECG
signals
by
the
set
partitioning
in
hierarchical
trees
algorithm,
”
IEEE
T
ransactions
on
Biomedical
Engineering
,
v
ol.
47,
no.
7,
pp.
849–856,
2000,
doi:
10.1109/10.846678.
[4]
C
.
K.
Jha
and
M.
H.
K
olekar
,
“Electrocardiogram
data
compression
using
DCT
based
discrete
orthogonal
Stockwell
transform,
”
Biomedical
Signal
Processing
and
Control
,
v
ol.
46,
pp.
174–181,
2018,
doi:
10.1016/j.bspc.2018.06.009.
[5]
R
.
Gupta,
“Quality
a
w
are
compress
ion
of
electrocardiogram
using
principal
component
analysis,
”
Journal
of
Medical
Systems
,
v
ol.
40,
no.
5,
p.112,
2016,
doi:
10.1007/s10916-016-0468-7.
[6]
V
.
K
umar
,
S.
C.
Sax
ena,
V
.
K.
Giri,
and
D.
Singh,
“Impro
v
ed
modied
AZTEC
technique
for
ECG
data
compression:
Ef
fect
of
length
of
parabolic
lter
on
reconstructed
signal,
”
Computers
and
Electrical
Engineering
,
v
ol.
31,
no.
4–5,
pp.
334–344,
2005,
doi:
10.1016/j.compeleceng.2005.02.002.
[7]
J.
P
.
Abenstein
and
W
.
J.
T
ompkins,
“
A
ne
w
data-reduction
algorithm
for
real-time
ECG
analysis,
”
IEEE
T
ransactions
on
Biomedical
Engineering
,
v
ol.
29,
no.
1,
pp.
43–48,
1982,
doi:
10.1109/TBME.1982.324962.
[8]
U.
E.
Ruttimann
and
H.
V
.
Pipber
ger
,
“Compression
of
the
ECG
by
prediction
or
interpolation
and
entrop
y
encoding,
”
IEEE
T
rans-
actions
on
Biomedical
Engineering
,
v
ol.
26,
no.
11,
pp.
613–623,
1979,
doi:
10.1109/TBME.1979.326543.
[9]
R.
C.
Barr
,
S.
M.
Blanchard,
and
D.
A.
Dipersi
o,
“SAP
A-2
is
the
F
an,
”
IEEE
T
ransactions
on
Biomedical
Engineering
,
v
ol.
32,
no.
5,
p.
337,
1985,
doi:
10.1109/TBME.1985.325548.
[10]
G.
Na
v
e
and
A.
Cohen,
“EC
G
compression
using
long-term
prediction,
”
IEEE
T
ransactions
on
Biomedical
Engineering
,
v
ol.
40,
no.
9,
pp.
877–885,
1993,
doi:
10.1109/10.245608.
[11]
F
.
W
ang
et
al.
,
“
A
no
v
el
ECG
signal
compression
m
ethod
using
spindle
con
v
olutional
auto-encoder
,
”
Computer
Methods
and
Pro-
grams
in
Biomedicine
,
v
ol.
175,
pp.
139–150,
2019,
doi:
10.1016/j.cmpb
.2019.03.019.
[12]
C.
J.
Deepu
and
Y
.
Lian,
“
A
Joint
QRS
detection
and
data
compression
scheme
for
wearable
sensors,
”
IEEE
T
ransactions
on
Biomedical
Engineering
,
v
ol.
62,
no.
1,
pp.
165–175,
2015,
doi:
10.1109/TBME.2014.2342879.
[13]
P
.
S.
Hamilton
and
W
.
J.
T
ompkins,
“Compression
of
the
amb
ulatory
ECG
by
a
v
erage
beat
subtraction
and
residual
dif
ferencing,
”
IEEE
T
ransactions
on
Biomedical
Engineering
,
v
ol.
38,
no.
3,
pp.
253–259,
1991,
doi:
10.1109/10.133206.
[14]
D.
Y
ang,
M.
Z.
Qin,
and
B.
Xu,
“ECG
compression
algorithm
based
on
empirical
mode
decomposition,
”
International
Journal
of
Signal
Processing,
Image
Processing,
and
P
attern
Recognition
,
v
ol.
8,
no.
2,
pp.
165–174,
2015.
[15]
B.
R.
S.
Reddy
and
I.
S.
N.
Murth
y
,
“ECG
data
compression
using
fourier
descriptors,
”
IEEE
T
ransactions
on
Biomedical
Engineer
-
ing
,
v
ol.
33,
no.
4,
pp.
428–434,
1986,
doi:
10.1109/TBME.1986.325799.
[16]
O.
El
B’charri,
R.
Latif,
W
.
Jenkal,
and
A.
Abenaou,
“The
ECG
signal
compression
using
an
ef
cient
algorithm
based
on
the
D
WT
,
”
International
Journal
of
Adv
anced
Computer
Science
and
Applications
,
v
ol.
7,
no.
3,
pp.
181-187,
2016,
doi:
10.14569/ijacsa.2016.070325.
[17]
L.
V
.
Batista,
E.
U.
K.
Melcher
,
and
L.
C.
Carv
alho,
“Compression
of
ECG
signals
by
optimized
quantization
of
discrete
cosine
transform
coef
cients,
”
Medical
Engineering
and
Ph
ysics
,
v
ol.
23,
no.
2,
pp.
127–134,
2001,
doi:
10.1016/S1350-4533(01)00030-3.
[18]
A.
Albiol
Colomer
and
A.
Albiol
Colomer
,
“
Adapti
v
e
ECG
data
compression
using
discrete
le
gendre
transform,
”
Digital
Signal
Processing
,
v
ol.
7,
no.
4,
pp.
222–228,
1997,
doi:
10.1006/dspr
.1997.0295.
[19]
T
.
Blanchett,
G.
C.
K
ember
,
and
G.
A.
Fenton,
“KL
T
-Based
quality
controlled
compression
of
single-lead
ECG,
”
IEEE
T
ransactions
on
Biomedical
Engineering
,
v
ol.
45,
no.
7,
pp.
942–945,
1998,
doi:
10.1109/10.686803.
[20]
W
.
Philips
and
G.
De
Jonghe,
“Data
compression
of
ECG’
s
by
high-de
gree
polynomial
approximation,
”
IEEE
T
ransactions
on
Biomedical
Engineering
,
v
ol.
39,
no.
4,
pp.
330–337,
1992,
doi:
10.1109/10.126605.
[21]
O.
P
.
Y
ada
v
and
S.
Ray
,
“
A
no
v
el
method
of
preprocessing
and
modeling
ECG
signals
with
Lagrange–Chebyshe
v
interpolat-
ing
polynomials,
”
International
Journal
of
System
Assurance
Engineering
and
Management
,
v
ol.
12,
no.
3,
pp.
377–390,
2021,
doi:
10.1007/s13198-021-01077-z.
[22]
G.
Li
and
C.
W
en,
“Le
gendre
polynomials
in
signal
reconstruction
and
compression,
”
in
Proceedings
of
the
2010
5th
IEEE
Confer
-
ence
on
Industrial
Electronics
and
Applications,
ICIEA
2010
,
Jun.
2010,
pp.
1636–1640,
doi:
10.1109/ICIEA.2010.5514776.
[23]
D.
Tchiotsop,
D.
W
olf,
V
.
Louis-Dorr
,
and
R.
Husson,
“ECG
data
compression
using
Jacobi
polynomials,
”
in
An-
nual
International
Conference
of
the
IEEE
Engineering
in
Medicine
and
Biology
-
Proceedings
,
2007,
pp.
1863–1867,
doi:
10.1109/IEMBS.2007.4352678.
[24]
M.
Karcze
wicz
and
M.
Gabbouj,
“ECG
data
compression
by
spline
approximation,
”
Signal
Processing
,
v
ol.
59,
no.
1,
pp.
43–59,
May
1997,
doi:
10.1016/s0165-1684(97)00037-6.
[25]
R.
Nyg
aard
and
D.
Haugland,
“Compressing
ECG
signals
by
piece
wise
polynomial
approximation,
”
in
ICASSP
,
IEEE
International
Conference
on
Acoustics,
Speech
and
Signal
Processing
-
Proceedings
,
1998,
v
ol.
3,
pp.
1809–1812,
doi:
10.1109/ICASSP
.1998.681812.
[26]
S.
Khetk
eeree
and
C.
Chansamorn,
“Signal
reconstruction
using
second
order
tetration
polynomial,
”
in
34th
International
T
ech-
nical
Conference
on
Circuits/Systems,
Computers
and
Communications,
ITC-CSCC
2019
,
2019,
pp.
1–4,
doi:
10.1109/ITC-
CSCC.2019.8793435.
[27]
A.
L.
Goldber
ger
et
al.
,
“Ph
ysioBank,
Ph
ysioT
oolkit,
and
Ph
ysioNet:
components
of
a
ne
w
research
resource
for
comple
x
ph
ysio-
logic
signals.,
”
Circulation
,
v
ol.
101,
no.
23,
2000,
doi:
10.1161/01.cir
.101.23.e215.
[28]
R.
Q.
Hamza,
K.
S.
Rijab,
and
M.
A.
R.
Hussien,
“Ef
cient
electrocardiogram
signal
compression
algorithm
using
dual
encoding
technique,
”
Indonesian
Journal
of
Electrical
Engineering
and
Computer
Science
(IJEECS)
,
v
ol.
25,
no.
3,
pp.
1529-1538,
2022,
doi:
10.11591/ijeecs.v25.i3.pp1529-1538.
Indonesian
J
Elec
Eng
&
Comp
Sci,
V
ol.
36,
No.
2,
No
v
ember
2024:
837–845
Evaluation Warning : The document was created with Spire.PDF for Python.
Indonesian
J
Elec
Eng
&
Comp
Sci
ISSN:
2502-4752
❒
845
BIOGRAPHIES
OF
A
UTHORS
Dr
.
Shashwati
Ray
recei
v
ed
her
B.Sc.
(Engg.)
de
gree
in
Electrical
Engineering
and
M.T
ech.
de
gree
in
Control
Systems
from
National
Institute
of
T
echnol
ogy
,
K
urukshetra,
India,
and
her
Ph.D.
de
gree
from
Indian
Institute
of
T
echnology
,
Bombay
,
India
in
2007.
She
has
been
a
professor
in
the
Department
of
Electrical
Engineering,
Bhilai
Institute
of
T
echnology
,
Dur
g,
India.
Her
research
interests
include
interv
al
analysis
techniques,
numerical
analysis,
optimization,
po
wer
system
control,
rene
w
able
ener
gy
sources,
rob
ust
control,
and
signal
processing.
She
can
be
contacted
at
email:
shashw
atiray@yahoo.com.
V
andana
Chouhan
recei
v
ed
her
B.E.
de
gree
in
electronics
engineering
from
CEC,
Chandrapur
,
India.
M.T
ech.
de
gree
in
Instrumentation
and
M.E
in
En
vironmental
from
CSVTU,
Bhilai.
She
is
pursuing
her
Ph.D.
de
gree
from
CSVTU,
Bhilai.
Her
research
interests
include
biomed-
ical
and
digital
signal
processing.
She
ca
n
be
contacted
at
email:
v
andanachouhan2212@yahoo.com.
Electr
ocar
dio
gr
am
r
econstruction
based
on
Hermite
interpolating
...
(Shashwati
Ray)
Evaluation Warning : The document was created with Spire.PDF for Python.