CN114287919A - J wave positioning method, device, equipment and medium based on cardiac shock signal - Google Patents
J wave positioning method, device, equipment and medium based on cardiac shock signal Download PDFInfo
- Publication number
- CN114287919A CN114287919A CN202111529493.4A CN202111529493A CN114287919A CN 114287919 A CN114287919 A CN 114287919A CN 202111529493 A CN202111529493 A CN 202111529493A CN 114287919 A CN114287919 A CN 114287919A
- Authority
- CN
- China
- Prior art keywords
- wave
- peak
- signal
- sampling
- preset
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000000747 cardiac effect Effects 0.000 title claims abstract description 47
- 230000035939 shock Effects 0.000 title claims abstract description 46
- 238000005070 sampling Methods 0.000 claims abstract description 176
- 238000012216 screening Methods 0.000 claims abstract description 20
- 238000007781 pre-processing Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims description 24
- 230000004807 localization Effects 0.000 claims description 20
- 230000010354 integration Effects 0.000 claims description 9
- 208000010125 myocardial infarction Diseases 0.000 claims description 6
- 210000001112 cardioblast Anatomy 0.000 description 14
- 238000001914 filtration Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000009499 grossing Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 238000009610 ballistocardiography Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The embodiment of the application relates to the field of smart home, and discloses a J-wave positioning method, device, equipment and medium based on a cardiac shock signal. Wherein the method comprises the following steps: acquiring a cardiac shock sampling signal with a preset sampling frequency; preprocessing the cardioimpact sampling signal to obtain preprocessed signal data, wherein the amplitude difference between any peak of the preprocessed signal data and a first reference point is larger than the amplitude difference between a corresponding peak of the cardioimpact sampling signal and a second reference point; acquiring a plurality of wave crests and a plurality of wave troughs in the preprocessed signal data; screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs; and positioning the J wave in the ballistocardiographic sampling signal according to the effective wave crest. The positioning accuracy of the J wave of BCG that places thicker intelligent house mattress in can improving in this application and gather.
Description
Technical Field
The embodiment of the application relates to the field of smart home, in particular to a J-wave positioning method, device, equipment and medium based on a cardiac shock signal.
Background
Ballistocardiogram (BCG) can be measured by a non-contact piezoelectric sensor in a portable manner without wearing cumbersome electrodes. BCG is the instantaneous impact force generated by the heart beat on the contact surface collected by a piezoelectric sensor, and it contains abundant physiological information related to sleep, such as respiration, heart rate, etc. In the waveform of the heart attack signal, the main wave peaks of the BCG have H, I, J, K wave peaks and the formed HIJK wave group has close relation with the heart beat.
In the process of implementing the embodiment of the present application, the inventors of the present application find that: generally, BCG can calculate heart rate and heart rate variability, and therefore, J-waves of BCG need to be located to find out the position of J-waves generated by the heart beat, and further calculate the detection index related to the heart rate. However, if the piezoelectric sensor for acquiring the BCG is embedded in a thick smart home mattress, the position of the acquired J wave of the BCG is not obvious.
Disclosure of Invention
The embodiment of the application aims to provide a method, a device, equipment and a medium for positioning J waves based on a cardiac shock signal, and the positioning accuracy of the J waves of BCG acquired by a piezoelectric sensor arranged in a thick smart home mattress can be improved.
In order to solve the technical problem, the embodiment of the application adopts the following technical scheme:
in a first aspect, an embodiment of the present application provides a method for J-wave localization based on a ballistocardiogram signal, including:
acquiring a cardiac shock sampling signal with a preset sampling frequency;
preprocessing the cardioimpact sampling signal to obtain preprocessed signal data, wherein the amplitude difference between any peak of the preprocessed signal data and a first reference point is larger than the amplitude difference between a corresponding peak of the cardioimpact sampling signal and a second reference point, the first reference point is a point on the preprocessed signal data except for the peak, the second reference point is a point on the cardioimpact sampling signal except for the peak, and the second reference point corresponds to the first reference point;
acquiring a plurality of wave crests and a plurality of wave troughs in the preprocessed signal data;
screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs;
and positioning the J wave in the ballistocardiographic sampling signal according to the effective wave crest.
In some embodiments, the preprocessing the ballistocardiograph sampling signal to obtain preprocessed signal data includes:
performing first-order difference processing on the cardiac shock sampling signal to obtain a first-order difference signal;
squaring the first-order difference signal to obtain a first signal;
and performing integration processing on the first signal to obtain preprocessed signal data.
In some embodiments, the integrating the first signal to obtain the pre-processed signal data includes:
performing integration processing on the first signal to obtain a second signal;
and denoising the second signal to obtain the preprocessed signal data.
In some embodiments, the screening out valid peaks from the preprocessed signal data according to the preset sampling frequency, the peaks and the troughs includes:
obtaining a plurality of first wave crests from the plurality of wave crests according to the plurality of wave crests and the plurality of wave troughs so as to obtain a first wave crest set;
traversing a plurality of first wave crests in the first wave crest set, and calculating the position distance between the first wave crests and second wave crests, wherein the first wave crests are currently traversed wave crests, and the second wave crests are adjacent to the first wave crests;
judging whether the position interval is larger than a preset interval or not, wherein the preset interval is obtained based on the preset sampling frequency;
if the position interval is larger than a preset interval, recording the first peak traversed currently, and taking the recorded first peak as an effective peak;
and if the position interval is not larger than the preset interval, rejecting the first peak traversed currently.
In some embodiments, said obtaining, from said number of peaks and said number of troughs, a number of first peaks among said number of peaks to obtain a first set of peaks comprises:
calculating a first difference between the amplitude of the first peak and the amplitude of the first trough and a second difference between the amplitude of the first peak and the amplitude of the second trough, wherein the first trough and the second trough are respectively points adjacent to the first peak;
judging whether the first difference value and the second difference value are both larger than a preset amplitude value;
and if the first difference value and the second difference value are both larger than a preset amplitude value, recording the first peak to obtain a first peak set.
In some embodiments, the method further comprises:
calculating the average value of the amplitudes of the plurality of wave peaks according to the plurality of wave peaks;
and obtaining the preset amplitude according to the average value.
In some embodiments, after the determining whether the first difference and the second difference are both greater than a preset magnitude, the method further includes:
and if the first difference is not larger than the preset amplitude and/or the second difference is not larger than the preset amplitude, rejecting the first peak.
In some embodiments, said locating a J-wave in said ballistocardiographic sampling signal from said effective peak comprises:
acquiring a reference peak corresponding to the effective peak from the heart attack sampling signal;
and taking the reference peak as a center, and acquiring the reference peak with the maximum amplitude as a J wave within a preset range of the reference peak.
In a second aspect, an embodiment of the present application further provides a J-wave localization apparatus based on a ballistocardiographic signal, where the apparatus includes:
the sampling signal acquisition module is used for acquiring a cardiac shock sampling signal with a preset sampling frequency;
the preprocessing module is used for preprocessing the cardioimpact sampling signal to obtain preprocessed signal data, wherein the amplitude difference between any peak of the preprocessed signal data and a first reference point is larger than the amplitude difference between a corresponding peak of the cardioimpact sampling signal and a second reference point, the first reference point is a point on the preprocessed signal data except for the peak, the second reference point is a point on the cardioimpact sampling signal except for the peak, and the second reference point corresponds to the position of the first reference point;
a peak and trough obtaining module, configured to obtain a plurality of peaks and a plurality of troughs in the preprocessed signal data;
the screening module is used for screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs;
and the positioning module is used for positioning the J wave in the cardiac shock sampling signal according to the effective wave crest.
In a third aspect, the present application further provides a J-wave localization apparatus based on ballistocardiogram signals, including:
at least one processor, and
a memory communicatively connected to the processor, the memory storing instructions executable by the at least one processor to enable the at least one processor to perform the method according to the first aspect.
In a fourth aspect, the present application further provides a non-transitory computer-readable storage medium, wherein the computer-readable storage medium stores computer-executable instructions, which, when executed by a ballistocardiograph-based J-wave localization apparatus, cause the ballistocardiograph-based J-wave localization apparatus to perform the method according to any one of the first aspect.
The beneficial effects of the embodiment of the application are as follows: different from the situation of the prior art, the J-wave positioning method, device, equipment and medium based on the cardiac shock signal provided in the embodiment of the present application obtain the cardiac shock sampling signal with the preset sampling frequency, then pre-process the cardiac shock sampling signal to obtain the pre-processed signal data, where an amplitude difference between any peak of the pre-processed signal data and a first reference point is greater than an amplitude difference between a corresponding peak of the cardiac shock sampling signal and a second reference point, where the first reference point is a point on the pre-processed signal data except for the peak, the second reference point is a point on the cardiac shock sampling signal except for the peak, and the second reference point corresponds to the position of the first reference point, so as to achieve the salient processing of the peak position in the cardiac shock sampling signal, so that the peak position is more obvious, the sensitivity of the wave crest in the amplitude of the whole cardiac shock sampling signal is improved, and J wave confirmation is facilitated; acquiring a plurality of wave crests and a plurality of wave troughs in the preprocessed signal data; screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs to obtain wave crests matched with a J wave band; and positioning the J wave in the cardiac shock sampling signal according to the effective wave peak, namely, taking the effective wave peak as a reference, accurately positioning the J wave in the cardiac shock sampling signal and improving the accuracy of J wave positioning.
Drawings
One or more embodiments are illustrated by way of example in the accompanying drawings, which correspond to the figures in which like reference numerals refer to similar elements and which are not to scale unless otherwise specified.
FIG. 1 is a schematic flow chart diagram illustrating an embodiment of a method for J-wave localization based on ballistocardiographic signals according to the present application;
FIG. 2 is a waveform diagram of a ballistocardiograph signal according to an embodiment of the method for J-wave localization based on ballistocardiograph signals;
FIG. 3 is a waveform diagram of the ballistocardiogram signal after the filtering process of FIG. 2;
FIG. 4 is a waveform diagram for the pre-processed signal data of FIG. 2;
FIG. 5 is a schematic flow chart diagram illustrating a further embodiment of the method for J-wave localization based on ballistocardiographic signals according to the present application;
FIG. 6 is a schematic structural diagram of an embodiment of the J-wave localization apparatus based on ballistocardiographic signals according to the present application;
fig. 7 is a hardware structure diagram of a controller in an embodiment of the J-wave localization apparatus based on ballistocardiogram signals.
Detailed Description
The present application will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the present application, but are not intended to limit the present application in any way. It should be noted that various changes and modifications can be made by one skilled in the art without departing from the spirit of the application. All falling within the scope of protection of the present application.
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the present application and are not intended to limit the present application.
It should be noted that, if not conflicted, the various features of the embodiments of the present application may be combined with each other within the scope of protection of the present application. Additionally, while functional block divisions are performed in apparatus schematics, with logical sequences shown in flowcharts, in some cases, steps shown or described may be performed in sequences other than block divisions in apparatus or flowcharts. Further, the terms "first," "second," "third," and the like, as used herein, do not limit the data and the execution order, but merely distinguish the same items or similar items having substantially the same functions and actions.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used in the description of the present application is for the purpose of describing particular embodiments only and is not intended to be limiting of the present application. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
In addition, the technical features mentioned in the embodiments of the present application described below may be combined with each other as long as they do not conflict with each other.
The method and the device for J-wave positioning based on the cardioblast signal can be applied to J-wave positioning equipment based on the cardioblast signal, and the J-wave positioning equipment based on the cardioblast signal comprises a controller, a piezoelectric sensor, a low-pass filter and a band-pass filter. The controller is respectively and electrically connected with the piezoelectric sensor, the low-pass filter and the band-pass filter, and the controller serves as a master control center. The piezoelectric sensor is used for collecting the instant impact force generated by the heart beat of the user on the contact surface of the piezoelectric sensor to obtain the BCG. The method comprises the steps that after a BCG is input into a low-pass filter to be filtered, the BCG is input into a band-pass filter to be filtered, a cardioshock sampling signal with a preset sampling frequency is obtained, the cardioshock sampling signal is preprocessed, preprocessed signal data is obtained, the amplitude difference between any peak of the preprocessed signal data and a first reference point is larger than the amplitude difference between a corresponding peak of the cardioshock sampling signal and a second reference point, the first reference point is a point on the preprocessed signal data except for the peak, the second reference point is a point on the cardioshock sampling signal except for the peak, the position of the second reference point corresponds to the position of the first reference point, and a plurality of peaks and a plurality of troughs are obtained from the preprocessed signal data; screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs; and positioning the J wave in the ballistocardiographic sampling signal according to the effective wave crest.
J ripples locating device based on heart impact signal can be great intelligent house mattress of thickness, and intelligent house mattress embeds piezoelectric sensor, gathers human heart impact signal (BCG) on the mattress.
The following J-wave positioning devices based on the cardiac shock signals are all exemplified by smart home mattresses.
Please refer to fig. 1, which is a flowchart illustrating an embodiment of a method for determining J-wave based on ballistocardiogram according to the present application, where the method may be executed by a controller in a J-wave determining apparatus based on ballistocardiogram, and the method includes steps S101 to S105.
S101: and acquiring a cardioblast sampling signal with a preset sampling frequency.
The user is located intelligent house mattress, and the built-in piezoelectric sensor of intelligence house mattress can gather the impact force in the twinkling of an eye that user's heart beats the production to the piezoelectric sensor contact surface in real time to obtain real-time BCG. The piezoelectric sensor can transmit the BCG to other devices for processing, and the other devices transmit the processed BCG to the controller, so that the controller acquires a heart attack sampling signal with preset sampling frequency. The controller may set the sampling frequency of the signals acquired from the piezoelectric sensor or other device, for example, the controller may set the sampling frequency of BCG acquired from the piezoelectric sensor to 120 Hz. As shown in fig. 2, fig. 2 is a waveform diagram of the acquired ballistocardiogram signal, and the ballistocardiogram signal is a BCG of 100s-160s of the user.
In some embodiments, in order to denoise the ballistocardiograph sampling signal and save the memory and the operation time of the J-wave positioning device based on the ballistocardiograph signal, the obtaining of the ballistocardiograph sampling signal with the preset sampling frequency specifically includes:
sampling the impact cardiac signal at a first frequency to obtain a first sampling signal;
inputting the first sampling signal into a low-pass filter for filtering to obtain a first filtering signal;
sampling the first filtering signal at a second frequency to obtain a second sampling signal;
and inputting the second sampling signal into a band-pass filter for filtering to obtain a cardioimpact sampling signal, and taking the second frequency as a preset sampling frequency, wherein the first frequency is greater than the second frequency.
Firstly, the controller samples a period of at least 4 seconds of the ballistocardiograph signal at a first frequency to obtain a first sampled signal, for example, the first frequency may be 120Hz, that is, the ballistocardiograph signal is sampled at a sampling frequency of 120Hz to obtain a first sampled signal, and then the first sampled signal is input to a low-pass filter for filtering to obtain a first filtered signal, the cut-off frequency of the low-pass filter may be 12Hz, and the low-pass filter may be a low-pass filter of 4-order butterworth, and by filtering the first sampled signal, the ballistocardiograph signal can be effectively denoised.
Then, sampling the first filtering signal at a second frequency to obtain a second sampling signal; and taking the second frequency as a preset sampling frequency, wherein the first frequency is greater than the second frequency. The second frequency may be 40Hz, after the second sampling signal is obtained, the second sampling signal is input to a band-pass filter for filtering, so as to obtain a cardioblast sampling signal, because the components of the heart rate signal are mainly concentrated in 2Hz to 12Hz, the cut-off frequencies of the band-pass filter are respectively 2Hz and 12Hz, after the second sampling signal is filtered by using the band-pass filters with the cut-off frequencies respectively 2Hz and 12Hz, the signals with the frequencies less than 2Hz and greater than 12Hz in the second sampling signal are filtered, and the signals with the frequencies between 2Hz and 12Hz are retained, so as to obtain the cardioblast sampling signal. As shown in fig. 3, after the ballistocardiogram signal of fig. 2 is sampled and filtered, the ballistocardiogram sampled signal shown in fig. 3 is obtained, it can be understood that the J wave in fig. 3 is only obtained after the subsequent positioning, and therefore, the ballistocardiogram sampled signal obtained after the ballistocardiogram signal of fig. 2 is sampled and filtered is substantially filtered data of the non-labeled J wave.
S102, preprocessing the cardioimpact sampling signal to obtain preprocessed signal data, wherein the amplitude difference between any peak of the preprocessed signal data and a first reference point is larger than the amplitude difference between a corresponding peak of the cardioimpact sampling signal and a second reference point, the first reference point is a point on the preprocessed signal data except for the peak, the second reference point is a point on the cardioimpact sampling signal except for the peak, and the second reference point corresponds to the position of the first reference point.
In order to highlight the peak, to highlight the position of the J wave, the cardiac shock sampling signal needs to be preprocessed first, and the preprocessing is performed on the cardiac shock sampling signal to obtain preprocessed signal data, which may include:
performing first-order difference processing on the cardiac shock sampling signal to obtain a first-order difference signal;
squaring the first-order difference signal to obtain a first signal;
and performing integration processing on the first signal to obtain preprocessed signal data.
Specifically, the first-order difference processing is performed on the cardiac shock sampling signal x [ n ], and for performing the first-order difference processing on the amplitude of each point in the cardiac shock sampling signal x [ n ], a first-order difference signal y [ n ] ═ x [ n ] -x [ n-1 is obtained.
Then, for the first order difference signal y [ n ]]The squaring process is performed, again on the first order difference signal y [ n ]]The amplitude of each point is squared to obtain a first signal z [ n ]]Namely: z [ n ]]＝(y[n])2The peak position is further highlighted by a squaring process.
Performing integration processing on the first signal to obtain a second signal s [ n ], which is obtained by calculating according to the following formula 1:
for example, when n is 3, the peak positions can be emphasized again by the integration processing for the points 1 to 5 in the first signal, so that the peak positions are more conspicuous.
As shown in fig. 4, fig. 4 is a waveform diagram of the preprocessed signal data, and the preprocessed signal data obtained by preprocessing the ballistocardiograph sampling signal of fig. 2 is shown in fig. 4. Obviously, the amplitude difference between any peak of the preprocessed signal data and a first reference point is greater than the amplitude difference between a corresponding peak of the cardioverter sampled signal and a second reference point, where the first reference point is a point on the preprocessed signal data other than the peak, the second reference point is a point on the cardioverter sampled signal other than the peak, and the second reference point corresponds to the position of the first reference point.
That is, in the preprocessed signal data, each point on the waveform corresponds to each point of the waveform of the non-preprocessed ballistocardiographic sampling signal, and includes that a certain peak of the preprocessed signal data corresponds to a certain peak of the ballistocardiographic sampling signal, and the position of the first reference point of the preprocessed signal data corresponds to the position of the second reference point of the ballistocardiographic sampling signal.
After the preprocessing, the peak position of the preprocessed signal data is highlighted relative to the amplitudes of other points (first reference points) in the preprocessed signal data, and the amplitude difference between any peak in the preprocessed signal data and the first reference point is larger than the amplitude difference between the corresponding peak in the non-preprocessed ballistocardiographic sampling signal and the second reference point; therefore, the peak of the preprocessed signal data is more pronounced than for the ballistocardiographic sampling signal.
After the second signal s [ n ] is obtained, in order to remove the influence of Gaussian noise on the waveform, denoising processing is carried out on the second signal, and the preprocessed signal data h [ n ] is obtained. The method specifically comprises the following steps: and performing average smoothing processing on the second signal, and calculating by adopting the following formula 2:
according to the formula 2, the average smoothing processing on the second signal is denoising processing by adopting a 7-point average smoothing processing mode, and gaussian noise in the second signal is eliminated through the denoising processing to obtain preprocessed signal data h [ n ].
S103: and acquiring a plurality of wave crests and a plurality of wave troughs in the preprocessed signal data.
In the preprocessed signal data, all maximum value points and corresponding positions thereof are obtained and are respectively recorded as a maximum value point amplitude value set Y ═ Y1,y2,y3……yn]And a maximum point position set P [ [ P ]1,p2,p3……pn]. Wherein, the maximum value point is a peak.
Correspondingly, all minimum value points and corresponding positions thereof are obtained and are respectively recorded as a minimum value point amplitude value set W ═ W1,w2,w3……wn]And a minimum point position set J ═ J1,j2,j3……jn]. Wherein, the minimum value point is a wave trough.
S104: and screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs.
In some embodiments, as shown in fig. 5, fig. 5 is a schematic flowchart of the process after being refined in step S104, and in step S104, the screening out effective peaks from the preprocessed signal data according to the preset sampling frequency, the peaks and the troughs may include:
s41: obtaining a plurality of first wave crests from the plurality of wave crests according to the plurality of wave crests and the plurality of wave troughs so as to obtain a first wave crest set;
s42: traversing a plurality of first wave crests in the first wave crest set, and calculating the position distance between the first wave crests and second wave crests, wherein the first wave crests are currently traversed wave crests, and the second wave crests are adjacent to the first wave crests;
s43: judging whether the position interval is larger than a preset interval or not, wherein the preset interval is obtained based on the preset sampling frequency;
s44: if the position interval is larger than a preset interval, recording the first peak traversed currently, and taking the recorded first peak as an effective peak;
s45: and if the position interval is not larger than the preset interval, rejecting the first peak traversed currently.
Specifically, after obtaining the peaks and the troughs, in order to remove the pseudo peaks from the peaks, obtaining a plurality of first peaks from the plurality of peaks according to the plurality of peaks and the plurality of troughs to obtain a first peak set, which may include:
calculating a first difference between the amplitude of the first peak and the amplitude of the first trough and a second difference between the amplitude of the first peak and the amplitude of the second trough, wherein the first trough and the second trough are respectively points adjacent to the first peak;
judging whether the first difference value and the second difference value are both larger than a preset amplitude value;
and if the first difference value and the second difference value are both larger than a preset amplitude value, recording the first peak to obtain a first peak set.
Specifically, traversing all maximum value points, calculating a first difference between the amplitude of the first peak and the amplitude of the first valley, and a second difference between the amplitude of the first peak and the amplitude of the second valley, where the first valley and the second valley are respectively adjacent points to the first peak, that is, respectively subtracting the amplitude of the currently traversed maximum value point and the amplitudes of the two adjacent valleys to obtain a first difference D [ n [ ]]And a second difference D [ n + 1]]I.e. D [ n ]]＝yn-wn，D[n+1]＝yn-wn+1。
Then, judging whether the first difference and the second difference are both larger than a preset amplitude, wherein the obtaining of the preset amplitude comprises:
calculating the average value of the amplitudes of the plurality of wave peaks according to the plurality of wave peaks;
and obtaining the preset amplitude according to the average value.
Obtaining maximum point amplitude value set Y ═ Y1,y2,y3……yn]Then, the maximum point amplitude value set Y is calculated as [ Y ═ Y1,y2,y3……yn]And obtaining an average value of the amplitudes of all the maximum value points (wave crests), recording the average value as M, and obtaining a preset amplitude value according to the average value M, wherein the preset amplitude value THR is 0.2M.
And judging whether the first difference and the second difference are both greater than a preset amplitude THR, if so, indicating that the first peak is not a pseudo peak, recording the first peak, and after traversing the peaks, recording all the first peaks which are not pseudo peaks to obtain a first peak set. It will be appreciated that recording the first peak includes recording the amplitude and position of the peak.
Correspondingly, if the first difference is not greater than the preset amplitude and/or the second difference is not greater than the preset amplitude, it is determined that the first peak is a pseudo peak, and the first peak is rejected as a pseudo peak.
After obtaining a first peak set P, obtaining all real peaks as first peaks, then traversing all the first peaks in the first peak set, and calculating the position distance between the first peaks and a second peak, wherein the first peaks are currently traversed peaks, the second peaks are peaks adjacent to the first peaks, and the position distance is recorded as beta [ n ]]＝pn-pn-1。
The preset interval is obtained based on said preset sampling frequency, in particular, since the maximum tolerated value of a typical adult heart beat is 300bpm, i.e. the heart rate value is maximally 300 beats per minute, the interval between adjacent J-waves is 40/(300/60) ═ 8 within the preset sampling rate with the second frequency of 40Hz, and thus the preset interval is set to 8.
For example, if a point adjacent to a currently traversed first peak is a second peak, and if a position interval between the first peak and the second peak is greater than a preset interval 8, it indicates that the first peak may be a J wave, and needs to be preserved, the currently traversed first peak is recorded, and the recorded first peak is taken as a valid peak.
Correspondingly, if the point adjacent to the currently traversed first peak is the second peak, if the position interval between the first peak and the second peak is not greater than the preset interval 8, then the first peak may not be a J wave and does not need to be kept, and therefore, the currently traversed first peak is eliminated.
And all the wave crests are effective wave crests after the situation that the position intervals between the currently traversed first wave crest and the adjacent second wave crest of the first wave crest are larger than the preset intervals is determined.
And S105, positioning the J wave in the heart shock sampling signal according to the effective wave crest.
Because the effective wave crest is obtained from the preprocessed cardiac shock sampling signal, the effective wave crest and the real J wave in the cardiac shock sampling signal have deviation, and the deviation distance is not large usually, the effective wave crest can be used as a reference to perform J point correction on the cardiac shock sampling signal.
In some embodiments, locating a J-wave in the ballistocardiographic sampling signal according to the effective peak may include:
acquiring a reference peak corresponding to the effective peak from the heart attack sampling signal;
and taking the reference peak as a center, and acquiring the reference peak with the maximum amplitude as a J wave within a preset range of the reference peak.
Specifically, in the ballistocardiograph sampling signal, a reference peak corresponding to the effective peak is obtained, that is, a reference peak position corresponding to the effective peak in the ballistocardiograph sampling signal x [ n ], for example, a 121 th sampling point of the effective peak in fig. 4, corresponds to a 120 th sampling point in the ballistocardiograph sampling signal x [ n ] in fig. 3, and therefore, when the effective peak is 121 sampling points, the reference peak corresponding to the effective peak being 121 sampling points in the ballistocardiograph sampling signal is obtained as the 120 th sampling point in fig. 3. Then, in the ballistocardiogram sampling signal, taking the reference peak as a center (taking the 120 th sampling point as an example), in a preset range (n-4, n +4) of the reference peak, finding out the reference peak with the largest amplitude as a J wave, thereby realizing real J wave position positioning, and as shown in fig. 3, finding out the reference peak with the largest amplitude among 7 points of the 117 th sampling point, the 118 th sampling point, the 119 th sampling point, the 120 th sampling point, the 121 th sampling point, the 122 th sampling point, and the 123 th sampling point as the J wave.
It is understood that the predetermined range may be selected as a predetermined range with the reference peak x [ n ] as the center and the neighborhood region as (n-4, n +4), the neighborhood region being a specific region, and any open region with the point n as the center, called the neighborhood region of the point a, denoted as u (n). δ neighborhood of point n: assuming that δ is a positive number, the open interval (n- δ, n + δ) is referred to as a δ neighborhood of the point n, the point n is referred to as the center of the neighborhood, and δ is referred to as the radius of the neighborhood, in this embodiment, δ is 4, so that there are 7 points in the 4 neighborhoods (n-4, n +4) of the reference peak n, and 7 points in the neighborhood are converted into time points, thereby obtaining a duration corresponding to a preset range, and the duration corresponding to the preset range corresponds to the duration of the HIJK wave group, where the HIJK wave group is a heartbeat wave group, and the duration of the HIJK wave group is the maximum duration of the heartbeat wave group.
As shown in FIG. 3, after the J wave is positioned, the circled points on the wave crest are all J points, so that the J wave can be accurately positioned in the ballistocardiographic sampling signal.
According to the embodiment of the application, a cardiac shock sampling signal with a preset sampling frequency is obtained, then the cardiac shock sampling signal is preprocessed, preprocessed signal data are obtained, the amplitude difference between any peak of the preprocessed signal data and a first reference point is larger than the amplitude difference between a corresponding peak of the cardiac shock sampling signal and a second reference point, wherein the first reference point is a point on the preprocessed signal data except for the peak, the second reference point is a point on the cardiac shock sampling signal except for the peak, and the second reference point corresponds to the position of the first reference point, so that the peak position in the cardiac shock sampling signal is highlighted, the position of the peak is more obvious, the sensitivity of the peak detected on the amplitude of the whole cardiac shock sampling signal is improved, and J waves are confirmed; acquiring a plurality of wave crests and a plurality of wave troughs in the preprocessed signal data; screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs to obtain wave crests matched with a J wave band; and positioning the J wave in the cardiac shock sampling signal according to the effective wave peak, namely, taking the effective wave peak as a reference, accurately positioning the J wave in the cardiac shock sampling signal and improving the accuracy of J wave positioning.
The embodiment of the present application further provides a J-wave positioning device based on a cardioblast signal, please refer to fig. 6, which shows a structure of the J-wave positioning device based on a cardioblast signal provided in the embodiment of the present application, and the J-wave positioning device 600 based on a cardioblast signal includes:
a sampling signal obtaining module 601, configured to obtain a cardiac shock sampling signal with a preset sampling frequency;
a preprocessing module 602, configured to preprocess the cardioverter sampling signal to obtain preprocessed signal data, where an amplitude difference between any peak of the preprocessed signal data and a first reference point is greater than an amplitude difference between a corresponding peak of the cardioverter sampling signal and a second reference point, where the first reference point is a point on the preprocessed signal data except for the peak, the second reference point is a point on the cardioverter sampling signal except for the peak, and the second reference point corresponds to a position of the first reference point;
a peak and trough obtaining module 603, configured to obtain a plurality of peaks and a plurality of troughs in the preprocessed signal data;
a screening module 604, configured to screen an effective peak from the preprocessed signal data according to the preset sampling frequency, the multiple peaks, and the multiple troughs;
and a positioning module 605, configured to position a J wave in the cardiac shock sampling signal according to the effective peak.
According to the embodiment of the application, a cardioblast sampling signal with a preset sampling frequency is obtained, then the cardioblast sampling signal is preprocessed, preprocessed signal data are obtained, the amplitude difference between any peak of the preprocessed signal data and a first reference point is larger than the amplitude difference between a corresponding peak of the cardioblast sampling signal and a second reference point, wherein the first reference point is a point on the preprocessed signal data except for the peak, the second reference point is a point on the cardioblast sampling signal except for the peak, and the second reference point corresponds to the first reference point in position, so that the peak is more obvious in position, the sensitivity of the peak detected on the amplitude of the whole cardioblast sampling signal is improved, and J waves are favorably confirmed; acquiring a plurality of wave crests and a plurality of wave troughs in the preprocessed signal data; screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs to obtain wave crests matched with a J wave band; and positioning the J wave in the cardiac shock sampling signal according to the effective wave peak, namely, taking the effective wave peak as a reference, accurately positioning the J wave in the cardiac shock sampling signal and improving the accuracy of J wave positioning.
In some embodiments, the pre-processing module 602 is further configured to:
performing first-order difference processing on the cardiac shock sampling signal to obtain a first-order difference signal;
squaring the first-order difference signal to obtain a first signal;
and performing integration processing on the first signal to obtain preprocessed signal data.
In some embodiments, the pre-processing module 602 is further configured to:
performing integration processing on the first signal to obtain a second signal;
and denoising the second signal to obtain the preprocessed signal data.
In some embodiments, the screening module 604 is further configured to:
obtaining a plurality of first wave crests from the plurality of wave crests according to the plurality of wave crests and the plurality of wave troughs so as to obtain a first wave crest set;
traversing a plurality of first wave crests in the first wave crest set, and calculating the position distance between the first wave crests and second wave crests, wherein the first wave crests are currently traversed wave crests, and the second wave crests are adjacent to the first wave crests;
judging whether the position interval is larger than a preset interval or not, wherein the preset interval is obtained based on the preset sampling frequency;
if the position interval is larger than a preset interval, recording the first peak traversed currently, and taking the recorded first peak as an effective peak;
and if the position interval is not larger than the preset interval, rejecting the first peak traversed currently.
In some embodiments, the screening module 604 is further configured to:
calculating a first difference between the amplitude of the first peak and the amplitude of the first trough and a second difference between the amplitude of the first peak and the amplitude of the second trough, wherein the first trough and the second trough are respectively points adjacent to the first peak;
judging whether the first difference value and the second difference value are both larger than a preset amplitude value;
and if the first difference value and the second difference value are both larger than a preset amplitude value, recording the first peak to obtain a first peak set.
In some embodiments, the screening module 604 is further configured to:
calculating the average value of the amplitudes of the plurality of wave peaks according to the plurality of wave peaks;
and obtaining the preset amplitude according to the average value.
In some embodiments, the screening module 604 is further configured to:
and if the first difference is not larger than the preset amplitude and/or the second difference is not larger than the preset amplitude, rejecting the first peak.
In some embodiments, the positioning module 605 is further configured to:
acquiring a reference peak corresponding to the effective peak from the heart attack sampling signal;
and taking the reference peak as a center, and acquiring the reference peak with the maximum amplitude as a J wave within a preset range of the reference peak.
In some embodiments, the sampling signal acquisition module 601 is configured to:
sampling the impact cardiac signal at a first frequency to obtain a first sampling signal;
inputting the first sampling signal into a low-pass filter for filtering to obtain a first filtering signal;
sampling the first filtering signal at a second frequency to obtain a second sampling signal;
and inputting the second sampling signal into a band-pass filter for filtering to obtain a cardioimpact sampling signal, and taking the second frequency as a preset sampling frequency, wherein the first frequency is greater than the second frequency.
It should be noted that the above-mentioned apparatus can execute the method provided by the embodiments of the present application, and has corresponding functional modules and beneficial effects for executing the method. For technical details which are not described in detail in the device embodiments, reference is made to the methods provided in the embodiments of the present application.
Fig. 7 is a schematic diagram of a hardware structure of a controller in an embodiment of a J-wave localization apparatus based on a ballistocardiogram signal, and as shown in fig. 7, the controller 11 includes:
one or more processors 111, memory 112. Fig. 7 illustrates an example of one processor 111 and one memory 112.
The processor 111 and the memory 112 may be connected by a bus or other means, and fig. 7 illustrates the connection by a bus as an example.
The memory 112, as a non-volatile computer-readable storage medium, may be used to store non-volatile software programs, non-volatile computer-executable programs, and modules, such as program instructions/modules corresponding to the J-wave localization method based on ballistocardiography signals in the embodiments of the present application (for example, the sampling signal obtaining module 601, the preprocessing module 602, the peak-valley obtaining module 603, the filtering module 604, and the localization module 605 shown in fig. 6). The processor 111 executes various functional applications and data processing of the controller by running the nonvolatile software programs, instructions and modules stored in the memory 112, that is, the method for J-wave localization based on ballistocardiography according to the above-described method embodiment is realized.
The memory 112 may include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function; the storage data area may store data created according to use of the person entry and exit detection apparatus, and the like. Further, the memory 112 may include high speed random access memory, and may also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other non-volatile solid state storage device. In some embodiments, the memory 112 optionally includes memory remotely located from the processor 111, which may be connected to a ballistocardiographic signal-based J-wave localization apparatus via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
The one or more modules are stored in the memory 112 and, when executed by the one or more processors 111, perform the method for shockwave based J-wave localization in any of the above-described method embodiments, e.g., performing the above-described method steps S101-S105 in fig. 1; the functions of the modules 601 and 605 in fig. 6 are realized.
The product can execute the method provided by the embodiment of the application, and has the corresponding functional modules and beneficial effects of the execution method. For technical details that are not described in detail in this embodiment, reference may be made to the methods provided in the embodiments of the present application.
The present application provides a non-transitory computer-readable storage medium, which stores computer-executable instructions, which are executed by one or more processors, such as one processor 111 in fig. 7, to enable the one or more processors to perform the method for locating J-waves based on ballistocardiograms in any of the above method embodiments, such as performing the above-described method steps S101 to S105 in fig. 1; the functions of the modules 601 and 605 in fig. 6 are realized.
The above-described embodiments of the apparatus are merely illustrative, and the units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment.
Through the above description of the embodiments, those skilled in the art will clearly understand that the embodiments may be implemented by software plus a general hardware platform, and may also be implemented by hardware. It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware related to instructions of a computer program, which can be stored in a computer readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), or the like.
Finally, it should be noted that: the above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; within the idea of the invention, also technical features in the above embodiments or in different embodiments may be combined, steps may be implemented in any order, and there are many other variations of the different aspects of the invention as described above, which are not provided in detail for the sake of brevity; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (11)
1. A J-wave positioning method based on a ballistocardiographic signal, the method comprising:
acquiring a cardiac shock sampling signal with a preset sampling frequency;
preprocessing the cardioimpact sampling signal to obtain preprocessed signal data, wherein the amplitude difference between any peak of the preprocessed signal data and a first reference point is larger than the amplitude difference between a corresponding peak of the cardioimpact sampling signal and a second reference point, the first reference point is a point on the preprocessed signal data except for the peak, the second reference point is a point on the cardioimpact sampling signal except for the peak, and the second reference point corresponds to the first reference point;
acquiring a plurality of wave crests and a plurality of wave troughs in the preprocessed signal data;
screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs;
and positioning the J wave in the ballistocardiographic sampling signal according to the effective wave crest.
2. The method of claim 1, wherein the pre-processing the ballistocardiographic sampling signal to obtain pre-processed signal data comprises:
performing first-order difference processing on the cardiac shock sampling signal to obtain a first-order difference signal;
squaring the first-order difference signal to obtain a first signal;
and performing integration processing on the first signal to obtain preprocessed signal data.
3. The method of claim 2, wherein the integrating the first signal to obtain the pre-processed signal data comprises:
performing integration processing on the first signal to obtain a second signal;
and denoising the second signal to obtain the preprocessed signal data.
4. The method of claim 1, wherein the screening the pre-processed signal data for valid peaks according to the predetermined sampling frequency, the peaks and the troughs comprises:
obtaining a plurality of first wave crests from the plurality of wave crests according to the plurality of wave crests and the plurality of wave troughs so as to obtain a first wave crest set;
traversing a plurality of first wave crests in the first wave crest set, and calculating the position distance between the first wave crests and second wave crests, wherein the first wave crests are currently traversed wave crests, and the second wave crests are adjacent to the first wave crests;
judging whether the position interval is larger than a preset interval or not, wherein the preset interval is obtained based on the preset sampling frequency;
if the position interval is larger than a preset interval, recording the first peak traversed currently, and taking the recorded first peak as an effective peak;
and if the position interval is not larger than the preset interval, rejecting the first peak traversed currently.
5. The method of claim 4, wherein obtaining a first number of peaks from the number of peaks and a number of valleys to obtain a first set of peaks comprises:
calculating a first difference between the amplitude of the first peak and the amplitude of the first trough and a second difference between the amplitude of the first peak and the amplitude of the second trough, wherein the first trough and the second trough are troughs adjacent to the first peak respectively;
judging whether the first difference value and the second difference value are both larger than a preset amplitude value;
and if the first difference value and the second difference value are both larger than a preset amplitude value, recording the first peak to obtain a first peak set.
6. The method of claim 5, further comprising:
calculating the average value of the amplitudes of the plurality of wave peaks according to the plurality of wave peaks;
and obtaining the preset amplitude according to the average value.
7. The method of claim 5, wherein after determining whether the first difference and the second difference are both greater than a predetermined magnitude, the method further comprises:
and if the first difference is not larger than the preset amplitude and/or the second difference is not larger than the preset amplitude, rejecting the first peak.
8. The method according to any one of claims 1 to 7, wherein said locating J-waves in said ballistocardiographic sampling signal according to said effective peak comprises:
acquiring a reference peak corresponding to the effective peak from the heart attack sampling signal;
and taking the reference peak as a center, and acquiring the reference peak with the maximum amplitude as a J wave within a preset range of the reference peak.
9. A J-wave localization apparatus based on ballistocardiographic signals, the apparatus comprising:
the sampling signal acquisition module is used for acquiring a cardiac shock sampling signal with a preset sampling frequency;
the preprocessing module is used for preprocessing the cardioimpact sampling signal to obtain preprocessed signal data, wherein the amplitude difference between any peak of the preprocessed signal data and a first reference point is larger than the amplitude difference between a corresponding peak of the cardioimpact sampling signal and a second reference point, the first reference point is a point on the preprocessed signal data except for the peak, the second reference point is a point on the cardioimpact sampling signal except for the peak, and the second reference point corresponds to the position of the first reference point;
a peak and trough obtaining module, configured to obtain a plurality of peaks and a plurality of troughs in the preprocessed signal data;
the screening module is used for screening effective wave crests from the preprocessed signal data according to the preset sampling frequency, the wave crests and the wave troughs;
and the positioning module is used for positioning the J wave in the cardiac shock sampling signal according to the effective wave crest.
10. A ballistocardiograph signal-based J-wave localization apparatus, comprising:
at least one processor, and
a memory communicatively coupled to the processor, the memory storing instructions executable by the at least one processor to enable the at least one processor to perform the method of any of claims 1-8.
11. A non-transitory computer-readable storage medium having stored thereon computer-executable instructions that, when executed by a ballistocardiograph-based J-wave localization device, cause the ballistocardiograph-based J-wave localization device to perform the method of any one of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111529493.4A CN114287919A (en) | 2021-12-14 | 2021-12-14 | J wave positioning method, device, equipment and medium based on cardiac shock signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111529493.4A CN114287919A (en) | 2021-12-14 | 2021-12-14 | J wave positioning method, device, equipment and medium based on cardiac shock signal |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114287919A true CN114287919A (en) | 2022-04-08 |
Family
ID=80968305
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111529493.4A Pending CN114287919A (en) | 2021-12-14 | 2021-12-14 | J wave positioning method, device, equipment and medium based on cardiac shock signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114287919A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116172527A (en) * | 2023-03-15 | 2023-05-30 | 康原(江苏)科技有限公司 | Intelligent health monitoring method, system and device |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105662375A (en) * | 2016-03-17 | 2016-06-15 | 广州中科新知科技有限公司 | Method and device for non-contact detecting vital sign signals |
| CN106691424A (en) * | 2015-12-18 | 2017-05-24 | 深圳市汇顶科技股份有限公司 | Method and device for detecting heart rate |
| CN107137071A (en) * | 2017-04-26 | 2017-09-08 | 可瑞尔科技(扬州)有限公司 | It is a kind of to analyze the method that heart impact signal is used for calculating short-term heart beat value |
| CN110115574A (en) * | 2018-02-07 | 2019-08-13 | 普天信息技术有限公司 | The method and apparatus of rhythm of the heart |
| CN110916636A (en) * | 2019-11-22 | 2020-03-27 | 新绎健康科技有限公司 | BCG signal heart rate calculation method and system based on dynamic second-order differential threshold |
| CN111084614A (en) * | 2019-04-16 | 2020-05-01 | 深圳市格兰莫尔科技有限公司 | Multi-sensor blood pressure detection device based on BCG (BCG-shaped BCG) signals |
| CN111938621A (en) * | 2020-06-29 | 2020-11-17 | 杨铭轲 | Heart rate detection method and device, computer equipment and readable storage medium |
| CN113499059A (en) * | 2021-06-01 | 2021-10-15 | 武汉理工大学 | BCG signal processing system and method based on optical fiber sensing non-contact |
| CN113647925A (en) * | 2021-07-05 | 2021-11-16 | 新绎健康科技有限公司 | Heart rate determination method and device based on heart attack signal |
-
2021
- 2021-12-14 CN CN202111529493.4A patent/CN114287919A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106691424A (en) * | 2015-12-18 | 2017-05-24 | 深圳市汇顶科技股份有限公司 | Method and device for detecting heart rate |
| CN105662375A (en) * | 2016-03-17 | 2016-06-15 | 广州中科新知科技有限公司 | Method and device for non-contact detecting vital sign signals |
| CN107137071A (en) * | 2017-04-26 | 2017-09-08 | 可瑞尔科技(扬州)有限公司 | It is a kind of to analyze the method that heart impact signal is used for calculating short-term heart beat value |
| CN110115574A (en) * | 2018-02-07 | 2019-08-13 | 普天信息技术有限公司 | The method and apparatus of rhythm of the heart |
| CN111084614A (en) * | 2019-04-16 | 2020-05-01 | 深圳市格兰莫尔科技有限公司 | Multi-sensor blood pressure detection device based on BCG (BCG-shaped BCG) signals |
| CN110916636A (en) * | 2019-11-22 | 2020-03-27 | 新绎健康科技有限公司 | BCG signal heart rate calculation method and system based on dynamic second-order differential threshold |
| CN111938621A (en) * | 2020-06-29 | 2020-11-17 | 杨铭轲 | Heart rate detection method and device, computer equipment and readable storage medium |
| CN113499059A (en) * | 2021-06-01 | 2021-10-15 | 武汉理工大学 | BCG signal processing system and method based on optical fiber sensing non-contact |
| CN113647925A (en) * | 2021-07-05 | 2021-11-16 | 新绎健康科技有限公司 | Heart rate determination method and device based on heart attack signal |
Non-Patent Citations (2)
| Title |
|---|
| 吴玲燕;杨涛;杨卫;吴小玲;: "基于差分方法的QRS波精确定位", 南京医科大学学报(自然科学版), no. 06, 15 June 2010 (2010-06-15), pages 56 - 61 * |
| 林海平;向馗;陈静;: "心电信号QRS波精确定位算法研究", 航天医学与医学工程, no. 04, 15 August 2009 (2009-08-15), pages 56 - 60 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116172527A (en) * | 2023-03-15 | 2023-05-30 | 康原(江苏)科技有限公司 | Intelligent health monitoring method, system and device |
| CN116172527B (en) * | 2023-03-15 | 2023-09-22 | 康原(江苏)科技有限公司 | Intelligent health monitoring method, system and device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rodríguez et al. | Feature extraction of electrocardiogram signals by applying adaptive threshold and principal component analysis | |
| Sadhukhan et al. | R-peak detection algorithm for ECG using double difference and RR interval processing | |
| Strasser et al. | Motion artifact removal in ECG signals using multi-resolution thresholding | |
| CN114159050A (en) | Respiration rate monitoring method, device, equipment and medium | |
| WO2017058595A1 (en) | Method, apparatus and system for biometric identification | |
| CN104616001A (en) | Fingerprint recognition system and fingerprint recognition method | |
| EP2459065A1 (en) | Method and apparatus for the analysis of a ballistocardiogram signal | |
| US11337638B2 (en) | ECG signal processing method and apparatus | |
| CN106618560B (en) | Method and device for processing brain wave signals | |
| CN103126669B (en) | Method of improving wave form drawing accuracy degree and drawing method and device applied method | |
| Matos et al. | Embedded system for individual recognition based on ECG Biometrics | |
| CN109009073A (en) | Atrial fibrillation detection device and storage medium | |
| KR20180072290A (en) | Method for generating electrocardiogram for identifying a person and method for identifying a person using the electrocardiogram | |
| CN107997754A (en) | A kind of intelligence mattress system and human body physiological characteristics data extraction method | |
| CN114287919A (en) | J wave positioning method, device, equipment and medium based on cardiac shock signal | |
| Sadhukhan et al. | Detection of ECG characteristic features using slope thresholding and relative magnitude comparison | |
| CN112190268A (en) | Physiological signal processing device | |
| Illanes-Manriquez et al. | An algorithm for robust detection of QRS onset and offset in ECG signals | |
| CN111582138A (en) | Electrocardio identity recognition method and system based on frequency domain cepstrum coefficient characteristics | |
| CN109077721A (en) | atrial fibrillation detection device and storage medium | |
| CN103040524A (en) | Device and method for reducing interference of physiological activities to medical imaging or measuring results | |
| CN114916940A (en) | Arrhythmia signal detection method, system and terminal based on expert knowledge | |
| Bashar et al. | Smartwatch based atrial fibrillation detection from photoplethysmography signals | |
| CN114557685A (en) | Non-contact motion robust heart rate measuring method and measuring device | |
| Rodríguez et al. | Adaptive threshold and principal component analysis for features extraction of electrocardiogram signals |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |