发明名称 Technique for real-time removal of power line interference in ECG
摘要 A technique for real-time removal of power line interference (PLI) in electrocardiograms (ECG) with sampling rate a relevant integer multiple of the nominal power line frequency (NPLF), comprising band-pass PLI filtering of the ECG, determining expedient weights and weighted Least-Squares parabolic approximating downsampled to the NPLF series of band-pass PLI filtered values. The disclosed technique utilizes only integer arithmetic, includes pre-adjusting for the routine ECG sampling rates, causes negligible ringing and is still accurate in the common occurrences of PLI with deviating from the nominal 60 or 50 Hz frequency and varying amplitude.
申请公布号 US9585621(B2) 申请公布日期 2017.03.07
申请号 US201313849666 申请日期 2013.03.25
申请人 Mitov Iliya 发明人 Mitov Iliya
分类号 A61B5/0402;A61B5/00 主分类号 A61B5/0402
代理机构 代理人 Coblenz Michael
主权项 1. A technique for real-time removal of power line interference (PLI) in electrocardiographic (ECG) devices with sampling rate sf set equal to an integer multiple of the nominal power line frequency (NPLF, given by pf=60 or 50 Hz) in a range from 240 Hz to 6*pf Hz, the technique comprising the following steps performed in each of said ECG devices after every sampling: (a) storing into a part of the ECG device memory a current sample obtained by analog-to-digital conversion of the treated ECG signal as x(1) in a First-In-First-Out (FIFO) buffer of 16-bit integers x(j), where j=1, . . . ,nx, nx=2*k−1+sf/5 and k=sf/pf is 4, 5 or 6, computing a 32-bit difference d(.)=x(1)−x(1+k) of phased with the NPLF ECG samples x(1) and x(1+k) for said k=4, 5 or 6 in said FIFO buffer x(j) and storing said d(.) as d(1) in a FIFO buffer of 32-bit integers d(j) with j=1 , . . . ,4; (b) filtering ECG samples x(1), . . . ,x(4*k−3) with said k=4, 5 or 6 in said FIFO buffer x(j) by a series of two identical band-pass PLI filters that are based on double moving average and have scaled-up to obviate fractional numbers impulse responses [−1, −2, . . . , −(k−1), k2−k, −(k−1), . . . , −2, −1] for said k=4, 5 or 6; (c) determining a 64-bit integer error e(.)=(3*d(1)−4*d(2)−d(3)+2*d(4))2+(4*d(1)−7*d(2) +2*d(3)+d(4))2+(d(1)+2*d(2)−7*d(3)+4*d(4))2+(2*d(1)−d(2)−4*d(3)+3*d(4))2 in Least-Squares (LS) linear approximation of said four 32-bit differences d(j) in said FIFO buffer, storing said e(.) as e(1) in a FIFO buffer of two 64-bit integers [e(1), e(2)], and storing f(.)=e(1)+e(2) as a new smoothed error f(1) in a FIFO buffer f(j) with j=1, . . . ,nf, where nf and a first index for use ff are pre-adjusted to dependent on said pf and sf values specified below; (d) determining a maximum mf of smoothed errors f(ff), . . . ,f(nj) in said FIFO buffer f(j), computing a ranged from 1 to 1000 integer weight w(.)=(mf+1000*da)/(mf+da), where da is parameter that for said ECG devices with analog-to-digital converters of 11- or 12-bit resolution over a 10 mV range is experimentally adjusted to da=3000, and storing said w(.) and its 64-bit integer product v(.) with a current output value of said band-pass PLI filters as w(1) and respectively v(1) in two FIFO buffers w(j) and v(j) with j=1, . . . ,nw, nw=2*sf/5+1 and lengths corresponding to 0.4 s; (e) estimating PLI by weighted LS parabolic approximating that uses a pair of current sub-series [w(1), w(k+1), . . . , w(nw)] and [v(1), v(k+1), . . . , v(nw)] for said k=4, 5 or 6 that at said pf=60 or 50 Hz have LS indices [−m, . . . ,m] with m=pf/5 and subsume each 2*m+1 of downsampled to the NPLF weights w(j) and products v(j) in said FIFO buffers, and subtracting thus obtained PLI estimate from ECG sample x(nx) that is last in said FIFO buffer x(j), thereby producing a cleaned of power line interference electrocardiogram for usable printing or displaying by said ECG devices.
地址 Vidin BG