CN105411575B - A kind of electrocardioelectrode position monitoring method - Google Patents
A kind of electrocardioelectrode position monitoring method Download PDFInfo
- Publication number
- CN105411575B CN105411575B CN201510865936.5A CN201510865936A CN105411575B CN 105411575 B CN105411575 B CN 105411575B CN 201510865936 A CN201510865936 A CN 201510865936A CN 105411575 B CN105411575 B CN 105411575B
- Authority
- CN
- China
- Prior art keywords
- ecg
- formula
- data
- numerical value
- electrocardiogram
- 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.)
- Expired - Fee Related
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The invention discloses a kind of electrocardioelectrode position monitoring method, for judging three electrode positions led during ECG detecting whether wrong, described method includes that simplifying electrode position judges step, raw data acquisition step, data reconstruction step, decision function construction step and last formal detection and the step judging electrode position whether wrong.The electrocardioelectrode position monitoring method of the present invention, at most have only to change six location status and be assured that correct electrode position, it is hereby achieved that be available for the electrocardiogram (ECG) data of practicality accurately, layman can be avoided cannot to judge the inaccurate defect of electrocardiogram (ECG) data obtained, the diagnosis and treatment avoiding the heart disease to detected personnel are forbidden, or even the serious consequence that delays treatment is brought.
Description
Technical field
The present invention relates to a kind of monitor the method that medical treatment detection device is the most proper use of, especially monitor the heart
The method whether position of electrical detection device in use each electrode has been placed into correct position,
Particularly a kind of electrocardioelectrode position monitoring method.
Background technology
Along with urbanization and the development of aging, existing medical system can't bear the heavy load, novel community
Medical treatment, community endowment, remotely see a doctor the most popular, but thing followed problem is health care professional
Lacking, considerable medical diagnosis, detection not only need specialized equipment, but also need specialization
Personnel correctly operate these equipment.
For as a example by the heart disease being common in middle-aged and elderly people, in order to prevent to diagnose early, typically in advance
Be required for using the electrocardiogram acquisition equipment of specialty to detect electrocardiogram (ECG) data, namely common people visual understanding
So-called thought-read electrograph, its most basic operation is accurately to install electrocardioelectrode with detected object.
Fig. 1 is shown that the three electrode position schematic diagrams leading in ECG detecting that prior art is conventional, i.e.
Such as figure, three ECG detecting that lead include seven electrodes, and wherein, the first positive pole led is expressed as CH1+
(electrode sequence number in standard 12 lead specification is expressed as V1), negative pole is expressed as CH1-(standard 12
The electrode sequence number led in specification is expressed as V2)；Second positive pole led is expressed as CH2+ (standard 12
The electrode sequence number led in specification is expressed as V3), negative pole is expressed as CH2-(in standard 12 lead specification
Electrode sequence number be expressed as V4)；3rd positive pole led is expressed as CH3+ (in standard 12 lead specification
Electrode sequence number be expressed as V5), negative pole is expressed as the CH3-(electrode sequence number in standard 12 lead specification
It is expressed as V6)；7th electrode RL is the ground electrode (electrode sequence number in standard 12 lead specification
It is expressed as V7).
The accurate location of V1-V7 electrode is respectively as follows: V1 electrode in intercostal space, left anterior axillary line the 5th；V2 electricity
Pole is at right clavicle and breastbone intersection；V3 electrode is in right border of sternum the 4th intercostal space；V4 electrode is at left lock
Bone and breastbone intersection；V5 electrode is at left side the 5th rib midline position；V6 electrode on presternum,
It is positioned under CH1-electrode and CH2-electrode；V7 electrode is in arcus costarum lower edge position, right side.
In existing universal standard specification, the color for the electrode wires of each electrode also has clear and definite regulation.
According to the standard of AHA (American Heart Association), the electrode wires color of V1-V7 electrode is respectively as follows: redness,
White, brown, black, orange, blue, green.Mark according to IEC (International Electrotechnical Commission)
Standard, the electrode wires color of V1-V7 electrode is respectively as follows: green, red, white, and yellow is orange, blue
Color, black.
The electrode position that shows from Fig. 1 it can be seen that the color of each electrode, position are different,
Needing suitable Professional knowledge could correctly operate the positioning of electrode of complexity, owing to circuit is more, location is multiple
Miscellaneous, non-professional detection doctor cannot be competent at, and therefore, ordinary individual has been difficult to the ECG detecting of specialty.
Although occurring in that some aim at the ECG detecting equipment of individual's design in the market, but structure being complicated, behaviour
Make bother the most very much, it is often more important that once electrode position place mistake, it is thus achieved that electrocardiogram (ECG) data be exactly not
Accurately, unpredictable serious consequence will be brought in the diagnosis and treatment basis in this, as heart disease.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of electrocardioelectrode position monitoring method, reducing or
Avoid problem noted earlier.
For solving above-mentioned technical problem, the present invention proposes a kind of electrocardioelectrode position monitoring method, is used for
Judging three electrode positions led during ECG detecting whether wrong, described three ECG detecting that lead include
CH1+, CH1-, CH2+, CH2-, CH3+, CH3-, RL totally seven electrodes, wherein said electrode
CH1+ and CH1-constitutes the first both positive and negative polarity led；Described electrode CH2+ and CH2-constitutes second and leads
Both positive and negative polarity；Described electrode CH3+ and CH3-constitutes the 3rd both positive and negative polarity led；Electrode RL is ground connection
Electrode；Described method comprises the steps:
Step one: simplify electrode position and judge, get rid of the described ground electrode RL away from remaining six electrode
Position wrong, the position of described negative electrode CH1-, CH2-and CH3-getting rid of three close positions connects
Wrong；
Step 2: raw data acquisition, according to correct electrode position connected mode, for different surveys
Examination object, uses same ECG detecting equipment to gather the correct electrocardiogram (ECG) data of multiple seasonal effect in time series and is stored in
In one raw data base, each time series comprises a plurality of described of one group of same time collection in interval
Correct electrocardiogram (ECG) data, each described correct electrocardiogram (ECG) data include described first leading, second lead and
3rd electrocardio test voltage CV1 led, CV2 and CV3；
Step 3: assuming that electrocardio test voltage CV1 in each described correct electrocardiogram (ECG) data, CV2 and
CV3, can carry out rebuilding acquisition by following formula 1, and described formula 1 is:
CV1=b11*1+b12*CV2+b13*CV3
CV2=b21*1+b22*CV1+b23*CV3
CV3=b31*1+b32*CV1+b33*CV2
Described correct electrocardiogram (ECG) data in described raw data base is substituted into described formula 1, calculates and obtain institute
State the matrix numerical value of coefficient bk in formula 1:
Step 4: the matrix numerical value calculating the described coefficient bk of acquisition in step 3 is substituted into described formula 1,
By described electrocardio test voltage CV1 of described for each in described raw data base correct electrocardiogram (ECG) data,
CV2 and CV3 substitutes into described formula 1 equally, correspondence can obtain one group of virtual electrocardio voltage DV1, DV2
And DV3；
Calculate electrocardio test voltage CV1 of correct electrocardiogram (ECG) data, CV2 described in each seasonal effect in time series one group
Described virtual correlation coefficient f1 electrocardio voltage DV1, DV2 and DV3 between corresponding with CV3,
F2 and f3；
Define a linear function formula 2:
Z=T0+T1*f1+T2*f2+T3*f3
Each group is calculated described correlation coefficient f1, f2 and the f3 substitution formula 2 obtained and all can obtain one
The function Z of individual correspondence, by each described function Z one decision-making formula 3 of substitution:
Determined equal to 1 by the described decision function g (Z) corresponding to correct electrocardiogram (ECG) data, solution formula
3, substitute into formula 2 by formula 3 by calculating each described function Z obtained, calculate and obtain in formula 2
T coefficient matrix numerical value:
T=[T0 T1 T2 T3]
Step 5: use ECG detecting equipment official testing electrocardiogram (ECG) data same in step 2, adopt equally
Step is judged, it is thus achieved that a plurality of official testing of multiple seasonal effect in time series with the described simplification electrode position of step one
Described the first of electrocardiogram (ECG) data lead, second lead and the 3rd formal electrocardio test voltage CV1 led,
CV2 and CV3；
The described formal heart of a plurality of official testing electrocardiogram (ECG) data described in multiple seasonal effect in time series that test is obtained
Electrical testing voltage CV1, CV2 and CV3 and step 3 calculate the matrix of the described coefficient bk of acquisition
Numerical value substitutes into described formula 1, and each time series correspondence obtains one group of virtual formal voltage DV1, DV2
And DV3；
Calculate formal electrocardio test voltage CV1 described in each seasonal effect in time series one group, CV2 and CV3 and its
Corresponding virtual formal correlation coefficient f1, f2 and f3 between voltage DV1, DV2 and DV3；Will
Calculate described correlation coefficient f1, f2 and the f3 obtained and step 4 calculates the described T coefficient matrix obtained
Numerical value substitutes into formula 2 and obtains the numerical value of linear function Z, and the numerical value of described function Z is substituted into decision-making formula
3, calculate the numerical value obtaining described decision function g (Z)；
If the demarcation numerical value that the numerical value calculating the described decision function g (Z) obtained sets more than or equal to one
T, then judge that electrode position does not has wrong, then by described formal electrocardio test voltage CV1, CV2 and
CV3 storage is standby；
If the numerical value calculating the described decision function g (Z) obtained is less than described demarcation numerical value t, then judge
Electrode position wrong also sends the prompting that reports an error, and repeats step 5, until obtaining correct described formal electrocardio
Test voltage CV1, CV2 and CV3 storage is standby.
Preferably, described step 2 farther includes data base and increases step: by described raw data base
In each described correct electrocardiogram (ECG) data in described electrocardio test voltage CV1, CV2 and CV3
Carry out permutation and combination, form five new wrong electrocardiogram (ECG) datas and be stored in described raw data base.
Preferably, described for each in described raw data base mistake electrocardiogram (ECG) data is substituted in step 4,
Determined equal to 0 by the described decision function g (Z) corresponding to mistake electrocardiogram (ECG) data, solution formula 3,
Substitute into formula 2 by formula 3 by calculating each described function Z obtained, calculate the T obtained in formula 2
Coefficient matrix numerical value.
Preferably, the described correct electrocardiogram (ECG) data in described raw data base and described mistake electrocardio number are utilized
Described T coefficient matrix numerical value is obtained according to using gradient descent method to calculate.
Preferably, the step gathering correct electrocardiogram (ECG) data described in described step 2 is: each test object
Follow-on test 24 hours, in these 24 hours, described correctly every test in 1 hour acquisition in 10 seconds
Electrocardiogram (ECG) data.
Preferably, in described step 3, the matrix numerical value of the coefficient bk in described formula 1 uses a young waiter in a wineshop or an inn
Multiplication calculates and obtains.
Preferably, in described step 3, the T coefficient matrix numerical value in described formula 2 uses gradient to decline
Method calculates and obtains.
Preferably, in described step 4, the demarcation numerical value t of described setting is more than or equal to 0.5.
Preferably, in described step 3, the matrix numerical value of the coefficient bk in described formula 1 is:
Preferably, in described step 4, the T coefficient matrix numerical value in described formula 2 is:
T=[-3.4115 7.7507 4.1454-4.5733].
The electrocardioelectrode position monitoring method of the present invention, at most has only to change six location status the most permissible
Determine correct electrode position, it is hereby achieved that be available for the electrocardiogram (ECG) data of practicality accurately, can avoid
Layman cannot judge the inaccurate defect of electrocardiogram (ECG) data obtained, it is to avoid the heart to detected personnel
The diagnosis and treatment of dirty disease are forbidden, or even the serious consequence that delays treatment is brought.
Accompanying drawing explanation
The following drawings is only intended to, in schematically illustrating the present invention and explaining, not limit the model of the present invention
Enclose.Wherein,
Fig. 1 is shown that the three electrode position schematic diagrams leading in ECG detecting that prior art is conventional；
Fig. 2 is shown that the curve synoptic diagram of the decision function that the decision model according to the present invention is drawn.
Detailed description of the invention
In order to the technical characteristic of the present invention, purpose and effect are more clearly understood from, now compare accompanying drawing
The detailed description of the invention of the present invention is described.Wherein, identical parts use identical label.
The most as described in the background section, lead for ECG detecting process even for three, be required for relating to
And the correct connectivity problem of seven electrodes, ordinary people is difficult to accurately complete the correct detection of electrocardiogram (ECG) data, because of
This is badly in need of a kind of method, knows electrode position whether wrong in order to monitor.
Based on the problems referred to above, the invention provides a kind of electrocardioelectrode position monitoring method, described method can
Lead the whether wrong of the electrode position during ECG detecting for monitoring three.Certainly, this is being described in detail
Before the method for invention, it is still desirable to be similar to background section and lead to three like that ECG detecting relates to
And each electrode re-define, owing to background parts has been carried out the most detailed with reference to Fig. 1
Illustrate, carry out the most simply defining only for electrode closely-related with the method for the present invention the most herein:
The three of the present invention ECG detecting that leads includes CH1+, CH1-, CH2+, CH2-, CH3+, CH3-,
RL totally seven electrodes, wherein electrode CH1+ and CH1-constitutes the first both positive and negative polarity led；Electrode CH2+
The second both positive and negative polarity led is constituted with CH2-；Electrode CH3+ and CH3-constitutes the 3rd both positive and negative polarity led；
Electrode RL is ground electrode.
Step one: first, before judging three electrode positions led during ECG detecting whether wrong,
Complex array combined result in view of seven malposition of electrode is the biggest, it is therefore desirable to carry out simplifying electrode position
Judge.I.e., first, get rid of the position wrong of the ground electrode RL away from remaining six electrode, with reference to Fig. 1
Visible, the 7th electrode RL, away from remaining six electrode, is very easy to connect, and the probability of wrong is very
Low, therefore the probability of this electrode RL position wrong is got rid of, subsequent step is made at eliminating equally
Reason can obtain same effect.
Then, the position wrong of negative electrode CH1-, CH2-and the CH3-of three close positions is got rid of.Row
After the situation of the 7th electrode position wrong, it is still desirable to consider three positions to totally six electrodes of leading
Put the situation of dislocation, it is assumed that six electrode positions arbitrarily exchange, and may result in 6！=720 kinds of possible schemes,
Judgement is got up highly difficult.But, fortunately, as it is shown in figure 1, three negative electrode CH1-, CH2-
Close with the position of CH3-, electric potential difference is the least, it can be assumed that they are essentially equipotential, and
If there is exchange three negative electrode positions, because their position is close, can ignore therefore make completely
The small electric potential difference become.Therefore simplify the negative pole position of three close positions of eliminating during electrode position judges to connect
Wrong situation is also rational, same, makees eliminating process equally and also can obtain same in subsequent step
Effect.
After carrying out simplifying electrode position judgement process, it is referred to Fig. 1 and finds out, due to three anelectrodes
CH1+, CH2+ and CH3+ are arranged in far-end, and apart from remote cause, electric potential difference is big each other, therefore, and this
Invent the determination methods about electrode position wrong, it is thus only necessary to consider when three positive electrode position are exchanged
Situation, the location status that i.e. three anelectrodes CH1+, CH2+ and CH3+ may possess is listed as follows:
Wherein, P1, P2, P3 represent three anelectrode CH1+, and the position at CH2+, CH3+ place is suitable
In sequence, such as location status 1, P1 position is the position at CH1+ place in Fig. 1, and P2 position is Fig. 1
The position at middle CH2+ place, P3 position is the position at CH3+ place in Fig. 1.Wherein, a total of six
Plant sequence of positions, it is assumed that it is correct for wherein only having location status 1, i.e. electrode CH1+, CH2+, CH3+
Respectively in P1 position, P2 position, P3 position.
Step 2: below it is envisaged that judge that, in six kinds of possible location status, any is correct
, remaining is all the position of mistake.
Therefore, step 2 first has to store the data of substantial amounts of tram state, i.e. passes through initial data
Gather, in order to as the basis judging wrong.I.e., in a specific embodiment, according to correct electricity
Position, pole connected mode, for different test objects, uses same ECG detecting equipment to gather multiple
The correct electrocardiogram (ECG) data of seasonal effect in time series is stored in a raw data base, and each time series comprises one group
Being spaced a plurality of correct electrocardiogram (ECG) data of same time collection, each correct electrocardiogram (ECG) data includes that first leads
Connection, second lead and the 3rd electrocardio test voltage CV1 led, CV2 and CV3.
Concrete, the step gathering correct electrocardiogram (ECG) data is: each test object follow-on test 24 hours,
In these 24 hours, within 10 seconds, obtain a correct electrocardiogram (ECG) data, the most each survey every test in 1 hour
Examination object follow-on test can obtain 24 correct electrocardiogram (ECG) datas for one day.Certainly, those skilled in the art
Even if it can be seen that such follow-on test a whole day, a test object also can only obtain little number
According to, in view of the surging situation of cost of labor, wanting to obtain abundant correct electrocardiogram (ECG) data is also one ten
Divide expensive work.
Therefore, invention further provides the step that a kind of data base increases, in order to reduce initial data
Collection capacity, reduce cost of labor.Particularly as follows:
By electrocardio test voltage CV1 in each correct electrocardiogram (ECG) data in raw data base, CV2 with
And CV3 carries out permutation and combination, form five new data and be stored in raw data base.
Such as, the original alignment order of a correct electrocardiogram (ECG) data in raw data base is, CV1, CV2
And CV3, therefore carry out permutation and combination according to according to optional position, six data altogether can be formed,
Wherein Article 1 data are original correct electrocardiogram (ECG) datas, and remaining five data is the positional fault newly increased
Wrong electrocardiogram (ECG) data, such as shown in following table:
CV1, CV2, CV3
CV2, CV1, CV3
CV3, CV2, CV1
CV1, CV3, CV2
CV3, CV1, CV2
CV2, CV3, CV1
The new data that combinations thereof is formed are also stored in raw data base, therefore pass through data base
Increase step and the data volume of raw data base can be expanded as original six times, greatly reduce collection number
According to workload, reduce cost.Certainly, raw data base needs which marks is correct electrocardio number
According to, which is mistake electrocardiogram (ECG) data, in order to the correct electrocardiogram (ECG) data of later use carries out gradient descent method instruction
Practice when solve coefficient matrix numerical value, after will to this further describe.
Step 3: data reconstruction.Assuming that the electrocardio in each correct electrocardiogram (ECG) data in raw data base
Test voltage CV1, CV2 and CV3, can carry out rebuilding acquisition, institute by following formula 1
Stating formula 1 is:
CV1=b11*1+b12*CV2+b13*CV3
CV2=b21*1+b22*CV1+b23*CV3
CV3=b31*1+b32*CV1+b33*CV2
The above-mentioned hypothesis rebuild is a kind of protocol step proposed based on the invention thinking, logical
Cross follow-up and true inspection, shown that this supposition possesses the value of practicality, invent as its principle
People is also difficult to complete analysis explanation, and those skilled in the art's hypothesis based on the present invention carries out operation and can obtain
Obtain result correctly.
Afterwards, the correct electrocardiogram (ECG) data in raw data base is substituted into described formula 1, calculate and obtain formula 1
In the matrix numerical value of coefficient bk:
From formula 1 it can be seen that total of three equation comprises nine coefficients, these nine coefficients be solved
Calculating, at least need three official testing electrocardiogram (ECG) datas, each data has three formal electrocardios to survey
Examination voltage, substitutes into formula 1 and can form nine equations, and simple solution can obtain the matrix function of coefficient bk
Value.Certainly, in order to reduce the error that individual data brings, it is possible to use such as 100 in raw data base
Bar, so 1000 data (include data base increase after data), use method of least square mode
The Fitting Calculation obtains the matrix numerical value of coefficient bk respectively, and the matrix numerical value of the coefficient bk being achieved in that will
Relatively accurately a lot.
It is the public affairs of art of mathematics about the method being obtained more exact computation results by least square fitting
Know general knowledge, by textbook or internet checking, those skilled in the art can know that matching of the present invention obtains
Ultimate principle to the matrix numerical value of coefficient bk: find the optimal of data by the quadratic sum minimizing error
Function mates so that between data and the real data tried to achieve, the quadratic sum of error is minimum.Due to minimum
Square law is a kind of method of known fitted data, and its point principle substantially is not the protection model of the present invention
Enclosing, therefore this is no longer going to repeat them.
In a specific embodiment, using same ECG detecting equipment collection, each test object is even
Continuous test 24 hours, in these 24 hours, obtains a correct electrocardio number in 10 seconds every test in 1 hour
According to, such as after test object is more than 100 people, calculated in the formula 1 obtained by above-mentioned steps
The matrix numerical value of coefficient bk is:
Step 4: decision function builds.The matrix numerical value calculating the coefficient bk of acquisition in step 3 is substituted into
Described formula 1, simultaneously by the electrocardio test voltage of each correct electrocardiogram (ECG) data in raw data base
CV1, CV2 and CV3 substitute into formula 1 equally, correspondence can obtain one group of virtual electrocardio voltage DV1,
DV2 and DV3.
Certainly, increase in step aforesaid data base, former by rearranging being stored in of combination acquisition
Each wrong electrocardiogram (ECG) data in beginning data base also is able to substitute into formula 1 equally, it is also possible to corresponding acquisition
One group of virtual electrocardio voltage, to be ready for use on the use of subsequent flows journey.
It should be noted that, the matrix numerical value of the coefficient bk owing to calculating acquisition in step 3 is many data
The result obtained after utilizing the mathematical method Optimal Fitting of method of least square etc is therefore correct by each
Electrocardio test voltage CV1 of electrocardiogram (ECG) data, the result that CV2 and CV3 obtains after substituting in formula 1 with
Small deviation can be there is in initial data, therefore through the calculated one group of numerical value of formula 1 be not true
Real numerical value, it may be that close to the virtual numerical value of actual value, be defined as the virtual heart the most in the present invention
Piezoelectric voltage DV1, DV2 and DV3, in order to make a distinction with actual value, certain those skilled in the art
Can also be defined by other title, the most only a title code name, people in the art
A virtual word therein should be interpreted as other implication any by member.
Subsequently, calculate electrocardio test voltage CV1 of each one group of correct electrocardiogram (ECG) data of seasonal effect in time series,
The virtual correlation coefficient f1 between electrocardio voltage DV1, DV2 and DV3 that CV2 and CV3 is corresponding,
F2 and f3, in order to the difference between these two groups of data of comparison.That is, it is possible to use known to art of mathematics
The formula calculating correlation coefficient calculates the phase between one group of DV1 and CV1 under each time series respectively
Close coefficient f1, the correlation coefficient f2 between one group of DV2 and CV2 under each time series and each
The correlation coefficient f3 between one group of DV3 and CV3 under time series.About the meter solving correlation coefficient
Calculation method is the common knowledge of art of mathematics, and its ultimate principle is not the scope of protection of present invention, this
Skilled person searches mathematical textbooks and can be readily available, and this is no longer going to repeat them.
Afterwards, one linear function formula 2 of definition:
Z=T0+T1*f1+T2*f2+T3*f3
By each group calculate obtain correlation coefficient f1, f2 and f3 substitute into formula 2 all can obtain one right
The function Z answered, certainly, due to the coefficient T 0, T1 in now formula 2, T2 and T3 is or unknown
, now function Z is not the numerical value determined, needs to solve acquisition coefficient T 0, T1 by following step,
After the concrete numerical value of T2 and T3, formula 2 could use in follow-up official testing.
Above-mentioned formula 2 is a kind of Scenario proposed based on the invention thinking, i.e. assumes each
Group meets linear relationship between correlation coefficient f1, f2 and f3, by follow-up and true inspection, shows
This supposition possesses the value of practicality, is also difficult to complete analysis explanation, this area as its principle inventor
Technical staff's hypothesis based on the present invention carries out operation can obtain correct result.
Then, each described function Z is substituted into a decision-making formula 3:
Now, if it is known that the numerical value of decision function g (Z), then the equation left side of formula 3 can be substituted into,
The concrete numerical value of the function Z in obtaining on the right of equation the most just can be reversely solved by formula 3, will meter
Calculate each described function Z obtained and substitute into formula 2, the T coefficient matrix numerical value in calculating acquisition formula 2:
T=[T0 T1 T2 T3]
About decision-making formula 3, it is that the one that inventor chooses from different decision schemes is applicable to this
Bright decision model, as in figure 2 it is shown, it is shown that the decision-making that the decision model according to the present invention is drawn
The curve synoptic diagram of function, wherein, using linear function Z as decision boundary function in Fig. 2, with transverse axis
Coordinate representation, ordinate of orthogonal axes is expressed as decision function g (Z).It can be seen in fig. 2 that when making with formula 3
Judge the decision function of electrode position whether wrong for the present invention in the case of, when decision function g (Z) is 1
Time, electrode position should be right-on, and when decision function g (Z) is 0, then electrode position should
It is full of prunes.
Therefore, principle based on decision-making formula 3 understands, by each correct electrocardio in raw data base
Data are after abovementioned steps converts, and the corresponding determining function g (Z) obtained just should be equal to 1；With
Sample, by each wrong electrocardiogram (ECG) data in raw data base after abovementioned steps converts, corresponding
The determining function g (Z) obtained just should be equal to 0, and (mistake electrocardiogram (ECG) data is increased step by aforesaid data base
Suddenly obtained).
Therefore, the functional value feature equal to 1 when utilizing decision-making formula 3 position correct, can be by original number
According to the correct electrocardiogram (ECG) data solution formula 3 in storehouse, then the concrete numerical value generation of the function Z obtained will be solved
Enter formula 2, it is possible to obtain T coefficient matrix numerical value.Same, it is also possible to utilize decision-making formula 3
The functional value feature equal to 0 during positional fault, by the wrong electrocardiogram (ECG) data solution formula in raw data base
3, then the concrete numerical value solving the function Z obtained is substituted into formula 2, equally obtain T coefficient square
Battle array numerical value.
Owing to, in formula 2, correlation coefficient f1, f2 and f3 are known, coefficient T 0, T1, T2 and
T3 is unknown, and now formula 2 is equal to a quaternary linear equation, minimum has only to four groups of phase relations
Number substitutes into formula 2 and is obtained with coefficient T 0, T1, the T coefficient matrix numerical value of T2 and T3.
Certainly, owing to the correct electrocardiogram (ECG) data in the raw data base of the present invention will not be only limitted to four certainly,
Meanwhile, the wrong data of five haplotype data amounts can be obtained by data base's expansion step, therefore can pass through
The T coefficient matrix number of the method acquisition present invention that the matrix numerical value of the coefficient bk in solution formula 1 is similar
Value.In a specific embodiment, present invention preferably employs gradient descent method to calculate in acquisition formula 2
T coefficient matrix numerical value.
Being similar to aforesaid method of least square, gradient descent method is also mathematic calculation known to one,
By textbook or internet checking, those skilled in the art can know that the present invention passes through gradient descent method
Calculate the ultimate principle obtaining T coefficient matrix numerical value: utilize negative gradient direction to determine the new of each iteration
The direction of search so that every time iteration can make object function to be optimized progressively reduce, and gradient descent method is
By continuous iteration to obtain the convergence of function parameter, within the specific limits, the scale of input data is big
The degree of convergence of little its result of calculation of decision, that is the data volume in raw data base is the biggest, the T obtained
Coefficient matrix numerical value closer to theoretical value (when accuracy exceedes certain threshold value, such as 99%, time then need not
Carry out more data again, because too much data can increase human cost, be likely to be due to input number simultaneously
Convergence concussion is caused) according to itself there is error.Owing to gradient descent method is the side of a kind of known fitted data
Method, its ultimate principle is not protection scope of the present invention, and therefore this is no longer going to repeat them.
In a specific embodiment, using same ECG detecting equipment collection, each test object is even
Continuous test 24 hours, in these 24 hours, obtains a correct electrocardio number in 10 seconds every test in 1 hour
According to, such as after test object is more than 100 people, calculated in the formula 2 obtained by above-mentioned steps
T coefficient matrix numerical value is:
T=[-3.4115 7.7507 4.1454-4.5733]
Step 5: formal detection with judge electrode position whether wrong.When operator formally start to detect,
It can use ECG detecting equipment official testing electrocardiogram (ECG) data same in step 2, uses step equally
The simplification electrode position of one judges step, it is thus achieved that multiple seasonal effect in time series a plurality of official testing electrocardiogram (ECG) data
First lead, second lead and the 3rd formal electrocardio test voltage CV1 led, CV2 and CV3.
The when of the most formally starting to detect, need to keep initial condition constant, i.e. assume the 7th electrode,
The position of whole three negative poles does not has wrong.
Then, the formal electrocardio of multiple seasonal effect in time series a plurality of official testing electrocardiogram (ECG) data test obtained
Test voltage CV1, CV2 and CV3 and step 3 calculate the matrix numerical value generation of the coefficient bk of acquisition
Enter described formula 1, each time series correspondence obtain one group of virtual formal voltage DV1, DV2 and
DV3。
Afterwards, formal electrocardio test voltage CV1 of each seasonal effect in time series one group, CV2 and CV3 are calculated
Corresponding virtual formal correlation coefficient f1, f2 and f3 between voltage DV1, DV2 and DV3；
Correlation coefficient f1, f2 and the f3 obtained and the T coefficient matrix numerical value generation of step 4 calculating acquisition will be calculated
Entering formula 2 and obtain the numerical value of linear function Z, the numerical value of function Z is substituted into decision-making formula 3, calculating obtains
Obtain the numerical value of decision function g (Z).
If the demarcation numerical value t that the numerical value calculating the decision function g (Z) obtained sets more than or equal to,
Then judging that electrode position does not has wrong, then by formal electrocardio test voltage CV1, CV2 and CV3 stores
Standby, the demarcation numerical value t of the most described setting can be more than or equal to 0.5.
If the numerical value calculating the decision function g (Z) obtained is less than described demarcation numerical value t, then judge electrode
Position wrong also sends the prompting that reports an error, and repeats step 5, until obtaining correct formal electrocardio test voltage
CV1, CV2 and CV3 storage is standby.
Such as, in the case of setting demarcation numerical value t is more than or equal to 0.5, if calculating the decision-making letter obtained
The numerical value of number g (Z) more than or equal to 0.5, then judges that electrode position does not has wrong.This is because it is public
T coefficient matrix numerical value in formula 2 is data or the number of whole errors present utilizing whole tram
According to calculating acquisition, thus the data of formal detection substitute into formula 2, then substitute into the g (Z) that formula 3 obtains
Only it is likely to occur two kinds of probabilities, close to 0 or close to 1, if greater than equal to 0.5, then according to certainly
It is correct that the curve law of plan function may determine that electrode position is particularly likely that.Actually if position
The numerical result of the decision function g (Z) that mistake obtains, it may be that be nearly equal to 0, sets and demarcates number
The threshold value that value t the most greatly relaxes judgement equal to 0.5 so that judge that the correct accuracy rate in position is higher,
Demarcating numerical value t equal to 0.3,0.4 being also feasible if namely set, being also possess practicality accurate
Rate, the accuracy rate only with respect to 0.5 is lower, same, and demarcating numerical value t can also set
Equal to 0.6,0.7, accuracy rate can be higher.Simply along with the raising of accuracy rate, correct for details of operation
Requirement higher, it is possible to need repeated multiple times measurement with comparison closest to that position accurately.
Whereas if judge electrode position wrong, then can be at the equipment by computer, mobile phone etc
The numerical value of comparison g (Z), less than in the case of t, controls speaker or light or display screen sends report
Wrong prompting, after tester experiences the prompting that reports an error, can change electrode position, retest and compare
To, judge electrode position whether wrong.
In a specific embodiment, testing staff can adjust one by one according to an electrode position status list
Whole electrode position, such as, the location status that three anelectrodes CH1+, CH2+ and CH3+ may possess
It is listed as follows:
If judging that location status 1 exists mistake after formal detection, then can adjust to location status 2 weight
Repetition measurement examination comparison, if it is determined that position is correct, then storage test data are standby.If position wrong, then
Changing to location status 3 again, the rest may be inferred, at most has only to change six location status and is assured that
Correct electrode position, it is hereby achieved that be available for the electrocardiogram (ECG) data of practicality accurately, can avoid background
The layman that technology segment is mentioned cannot judge the inaccurate defect of electrocardiogram (ECG) data obtained, it is to avoid right
The diagnosis and treatment of the heart disease of detected personnel are forbidden, or even the serious consequence that delays treatment is brought.
Although it will be appreciated by those skilled in the art that the present invention is to retouch according to the mode of multiple embodiments
State, but the most each embodiment only comprises an independent technical scheme.Description so describes
For the sake of being only used to understand, description should be understood by those skilled in the art as an entirety,
And technical scheme involved in each embodiment regarded as can be in the way of being mutually combined into different embodiment
Understand protection scope of the present invention.
The foregoing is only the schematic detailed description of the invention of the present invention, be not limited to the model of the present invention
Enclose.Any those skilled in the art, is made on the premise of without departing from the design of the present invention and principle
Equivalent variations, revise and combine, the scope of protection of the invention all should be belonged to.
Claims (9)
1. an electrocardioelectrode position monitoring method, the electrode during ECG detecting is led in monitoring three
Position whether wrong, described three ECG detecting that lead include CH1+, CH1-, CH2+, CH2-, CH3+,
CH3-, RL totally seven electrodes, wherein said electrode CH1+ and CH1-constitutes the first both positive and negative polarity led；
Described electrode CH2+ and CH2-constitutes the second both positive and negative polarity led；Described electrode CH3+ and CH3-is constituted
3rd both positive and negative polarity led；Electrode RL is ground electrode；It is characterized in that, described method includes as follows
Step:
Step one: simplify electrode position and judge, get rid of the described ground electrode RL away from remaining six electrode
Position wrong, the position of described negative electrode CH1-, CH2-and CH3-getting rid of three close positions connects
Wrong；
Step 2: raw data acquisition, according to correct electrode position connected mode, for different surveys
Examination object, uses same ECG detecting equipment to gather the correct electrocardiogram (ECG) data of multiple seasonal effect in time series and is stored in
In one raw data base, each time series comprises a plurality of described of one group of same time collection in interval
Correct electrocardiogram (ECG) data, each described correct electrocardiogram (ECG) data include described first leading, second lead and
3rd electrocardio test voltage CV1 led, CV2 and CV3；
Step 3: assuming that electrocardio test voltage CV1 in each described correct electrocardiogram (ECG) data, CV2 and
CV3, can carry out rebuilding acquisition by following formula 1, and described formula 1 is:
CV1=b11*1+b12*CV2+b13*CV3
CV2=b21*1+b22*CV1+b23*CV3
CV3=b31*1+b32*CV1+b33*CV2
Described correct electrocardiogram (ECG) data in described raw data base is substituted into described formula 1, calculates and obtain institute
State the matrix numerical value of coefficient bk in formula 1:
Step 4: the matrix numerical value calculating the described coefficient bk of acquisition in step 3 is substituted into described formula 1,
By described electrocardio test voltage CV1 of described for each in described raw data base correct electrocardiogram (ECG) data,
CV2 and CV3 substitutes into described formula 1 equally, correspondence can obtain one group of virtual electrocardio voltage DV1, DV2
And DV3；
Calculate electrocardio test voltage CV1 of correct electrocardiogram (ECG) data, CV2 described in each seasonal effect in time series one group
Described virtual correlation coefficient f1 electrocardio voltage DV1, DV2 and DV3 between corresponding with CV3,
F2 and f3；
Define a linear function formula 2:
Z=T0+T1*f1+T2*f2+T3*f3
Each group is calculated described correlation coefficient f1, f2 and the f3 substitution formula 2 obtained and all can obtain one
The function Z of individual correspondence, by each described function Z one decision-making formula 3 of substitution:
Determined equal to 1 by the described decision function g (Z) corresponding to correct electrocardiogram (ECG) data, solution formula
3, substitute into formula 2 by formula 3 by calculating each described function Z obtained, calculate and obtain in formula 2
T coefficient matrix numerical value:
T=[T0 T1 T2 T3]
Step 5: use ECG detecting equipment official testing electrocardiogram (ECG) data same in step 2, adopt equally
Step is judged, it is thus achieved that a plurality of official testing of multiple seasonal effect in time series with the described simplification electrode position of step one
Described the first of electrocardiogram (ECG) data lead, second lead and the 3rd formal electrocardio test voltage CV1 led,
CV2 and CV3；
The described formal heart of a plurality of official testing electrocardiogram (ECG) data described in multiple seasonal effect in time series that test is obtained
Electrical testing voltage CV1, CV2 and CV3 and step 3 calculate the matrix of the described coefficient bk of acquisition
Numerical value substitutes into described formula 1, and each time series correspondence obtains one group of virtual formal voltage DV1, DV2
And DV3；
Calculate formal electrocardio test voltage CV1 described in each seasonal effect in time series one group, CV2 and CV3 and its
Corresponding virtual formal correlation coefficient f1, f2 and f3 between voltage DV1, DV2 and DV3；Will
Calculate described correlation coefficient f1, f2 and the f3 obtained and step 4 calculates the described T coefficient matrix obtained
Numerical value substitutes into formula 2 and obtains the numerical value of linear function Z, and the numerical value of described function Z is substituted into decision-making formula
3, calculate the numerical value obtaining described decision function g (Z)；
If the demarcation numerical value that the numerical value calculating the described decision function g (Z) obtained sets more than or equal to one
T, then judge that electrode position does not has wrong, then by described formal electrocardio test voltage CV1, CV2 and
CV3 storage is standby；
If the numerical value calculating the described decision function g (Z) obtained is less than described demarcation numerical value t, then judge
Electrode position wrong also sends the prompting that reports an error, and repeats step 5, until obtaining correct described formal electrocardio
Test voltage CV1, CV2 and CV3 storage is standby.
2. the method for claim 1, it is characterised in that described step 2 farther includes number
Step is increased: by the described heart in described for each in described raw data base correct electrocardiogram (ECG) data according to storehouse
Electrical testing voltage CV1, CV2 and CV3 carry out permutation and combination, form five new wrong electrocardiogram (ECG) datas
And be stored in described raw data base.
3. method as claimed in claim 2, it is characterised in that each by described raw data base
Described in bar, mistake electrocardiogram (ECG) data substitutes in step 4, by the described decision function corresponding to mistake electrocardiogram (ECG) data
G (Z) is determined equal to 0, solution formula 3, will calculate each described function Z obtained by formula 3
Substitute into formula 2, calculate the T coefficient matrix numerical value obtained in formula 2.
4. method as claimed in claim 3, it is characterised in that utilize the institute in described raw data base
State correct electrocardiogram (ECG) data and described mistake electrocardiogram (ECG) data uses gradient descent method to calculate the described T coefficient square of acquisition
Battle array numerical value.
5. the method for claim 1, it is characterised in that gather correct described in described step 2
The step of electrocardiogram (ECG) data is: each test object follow-on test 24 hours, in these 24 hours, every
Test in 1 hour obtains a described correct electrocardiogram (ECG) data in 10 seconds.
6. the method for claim 1, it is characterised in that in described step 3, described formula 1
In coefficient bk matrix numerical value use method of least square calculate obtain.
7. the method for claim 1, it is characterised in that in described step 3, described formula 2
In T coefficient matrix numerical value use gradient descent method calculate obtain.
8. the method as described in one of claim 3-6, it is characterised in that in described step 3, described
The matrix numerical value of the coefficient bk in formula 1 is:
9. the method as described in one of claim 3-6, it is characterised in that in described step 4, described
T coefficient matrix numerical value in formula 2 is:
T=[-3.4115 7.7507 4.1454-4.5733].
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510865936.5A CN105411575B (en) | 2015-12-01 | 2015-12-01 | A kind of electrocardioelectrode position monitoring method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510865936.5A CN105411575B (en) | 2015-12-01 | 2015-12-01 | A kind of electrocardioelectrode position monitoring method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105411575A CN105411575A (en) | 2016-03-23 |
| CN105411575B true CN105411575B (en) | 2016-08-17 |
Family
ID=55490465
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510865936.5A Expired - Fee Related CN105411575B (en) | 2015-12-01 | 2015-12-01 | A kind of electrocardioelectrode position monitoring method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105411575B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105943020B (en) * | 2016-05-03 | 2018-10-23 | 东北大学 | A kind of 6 lead misconnection detecting system of cardiac electric axis and method |
| CN107374617A (en) * | 2017-08-24 | 2017-11-24 | 复旦大学附属中山医院 | A kind of automatic error-correcting electrocardiogram system and method |
| CN111491557B (en) * | 2017-12-15 | 2024-04-05 | 深圳迈瑞生物医疗电子股份有限公司 | Lead electrode identification device, method, storage medium and medical monitoring equipment |
| CN111345897B (en) * | 2018-12-20 | 2021-02-09 | 四川锦江电子科技有限公司 | Method and device for determining entrance and exit of catheter electrode into and out of sheath |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001204701A (en) * | 1999-12-31 | 2001-07-31 | Ge Marquette Medical Systems Inc | Method for identifying electrode arrangement |
| CN102793539A (en) * | 2012-08-31 | 2012-11-28 | 深圳市理邦精密仪器股份有限公司 | Vector cardiogram detection and analysis method and system |
| CN103431856A (en) * | 2013-08-30 | 2013-12-11 | 深圳市理邦精密仪器股份有限公司 | Method and device for selecting electrocardiogram lead in multiple lead synchronous electrocardiographic signals |
| CN103800004A (en) * | 2013-11-19 | 2014-05-21 | 邱磊 | Automatic detection method used for wrong placement of electrocardio-electrodes and based on kernel function classification algorithm |
| CN103908244A (en) * | 2014-04-03 | 2014-07-09 | 深圳市理邦精密仪器股份有限公司 | ECG lead wrong connection judging method and device |
| CN104856669A (en) * | 2015-04-30 | 2015-08-26 | 厦门纳龙科技有限公司 | Method for correcting electrocardiographic lead misconnection |
-
2015
- 2015-12-01 CN CN201510865936.5A patent/CN105411575B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001204701A (en) * | 1999-12-31 | 2001-07-31 | Ge Marquette Medical Systems Inc | Method for identifying electrode arrangement |
| CN102793539A (en) * | 2012-08-31 | 2012-11-28 | 深圳市理邦精密仪器股份有限公司 | Vector cardiogram detection and analysis method and system |
| CN103431856A (en) * | 2013-08-30 | 2013-12-11 | 深圳市理邦精密仪器股份有限公司 | Method and device for selecting electrocardiogram lead in multiple lead synchronous electrocardiographic signals |
| CN103800004A (en) * | 2013-11-19 | 2014-05-21 | 邱磊 | Automatic detection method used for wrong placement of electrocardio-electrodes and based on kernel function classification algorithm |
| CN103908244A (en) * | 2014-04-03 | 2014-07-09 | 深圳市理邦精密仪器股份有限公司 | ECG lead wrong connection judging method and device |
| CN104856669A (en) * | 2015-04-30 | 2015-08-26 | 厦门纳龙科技有限公司 | Method for correcting electrocardiographic lead misconnection |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105411575A (en) | 2016-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105411575B (en) | A kind of electrocardioelectrode position monitoring method | |
| CN105997072B (en) | A kind of electrode contact condition detection method of the electrical impedance detection based on measurement accuracy | |
| CN105496392B (en) | A kind of have electrode and connect the three of discrimination function and lead cardioelectric monitor system | |
| CN105380643B (en) | A kind of three connection method of discrimination leading electrocardioelectrode | |
| CN105286854B (en) | A kind of remote ecg monitoring method | |
| CN103886328A (en) | Functional magnetic resonance image data classification method based on brain network modular structure characteristics | |
| CN102670190A (en) | Heart rate variability nonlinear characteristic-based automatic diagnosis method for congestive heart failure | |
| CN106443326A (en) | Test system and method for fault locator of power distribution network | |
| CN105184107A (en) | Health risk pre-warning system for human body circulatory system | |
| CN106257511A (en) | A kind of grid faults characteristics quality testing method | |
| CN105310687B (en) | Dynamic electrocardiogram method for real-time monitoring based on mobile Internet | |
| CN117238226B (en) | LED display screen fault self-checking method and system | |
| Effendi et al. | The impact of EQ, SQ and IQ towards AQ using path analysis and Rasch logit for modelling purpose among mechanical engineering students | |
| CN105528857A (en) | Intelligent remote body sign data acquisition device | |
| CN106971361A (en) | A kind of professional ability remote test system and professional ability remote test method | |
| CN109886428A (en) | A kind of power equipment safety configuration check method | |
| CN107273984A (en) | A kind of oil for electric power chemical supervision management method and system | |
| CN110956330A (en) | Method and system for predicting line loss of power transmission line based on multidimensional influence quantity | |
| CN106503417A (en) | A kind of electrical equipment fault rate computational methods based on state evaluation | |
| CN106618551B (en) | A kind of intelligent terminal for being used for the connection of three lead electrocardioelectrodes and differentiating | |
| CN106725447B (en) | A kind of three lead electrocardioelectrodes based on feedforward neural network fitting connect method of discrimination | |
| CN105662396A (en) | Method for acquiring human body life data | |
| CN105092980B (en) | A kind of input and output impedance intelligent testing method | |
| CN211633227U (en) | Novel electronic visual chart | |
| Sidik et al. | The identification of outliers in wrapped normal data by using Ga statistics |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180514 Address after: 312000 ninth floor of Pingjiang Road 328, Yuecheng District, Shaoxing, Zhejiang. Patentee after: ZHEJIANG MEDZONE MEDICAL EQUIPMENT CO.,LTD. Address before: 312099 Pingjiang Road, Shaoxing, Shaoxing, Zhejiang, No. ninth building Co-patentee before: ZHEJIANG MEDZONE BIOMEDICAL MATERIAL AND DEVICE Research Institute Patentee before: ZHEJIANG MINGZHONG TECHNOLOGY Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200918 Address after: Room 703, 7 / F, No. 3766, Nanhuan Road, Binjiang District, Hangzhou City, Zhejiang Province Patentee after: ZHEJIANG MINGZHONG TECHNOLOGY Co.,Ltd. Address before: 312000, Zhejiang, Shaoxing City, Pingjiang Road, No. 328, ninth, one floor Patentee before: ZHEJIANG MEDZONE MEDICAL EQUIPMENT Co.,Ltd. |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 Termination date: 20211201 |