CN107374610A - Magnetocardiogram generation method and generation system - Google Patents
Magnetocardiogram generation method and generation system Download PDFInfo
- Publication number
- CN107374610A CN107374610A CN201710564362.7A CN201710564362A CN107374610A CN 107374610 A CN107374610 A CN 107374610A CN 201710564362 A CN201710564362 A CN 201710564362A CN 107374610 A CN107374610 A CN 107374610A
- Authority
- CN
- China
- Prior art keywords
- passage
- coordinate
- heart magnetic
- data
- magnetic data
- 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.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/242—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents
- A61B5/243—Detecting biomagnetic fields, e.g. magnetic fields produced by bioelectric currents specially adapted for magnetocardiographic [MCG] signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The present invention provides a kind of magnetocardiogram generation method and system applied to magnetocardiograph, including：1）The data of each passage a cycle are intercepted first, and a certain moment is chosen from the data of interception, take out the data at each passage moment, then calculate each passage corresponding position coordinates above chest；2）Allow the data of each passage sometime should with coordinate pair；3）The two-dimentional heart magnetic data arranged to opsition dependent carries out two-dimensional interpolation；4）Color table is inquired about to the 2-D data after interpolation；5）Draw magnetocardiogram.The magnetocardiogram that the present invention generates, image clearly, distribution of color is uniform, and+two dipoles can be clearly visible that in T wave bands.In addition, if the magnetocardiogram from the S sections of QRS wave to T sections by setting step-length picture at different moments, then it can be seen that the temporal-spatial evolution process of cardiomotility electric current, doctor can diagnose whether the patient suffers from heart disease according to the dipole number of magnetocardiogram and the evolution direction of cardiomotility electric current.
Description
Technical field
The invention belongs to biomedicine signals analysis field, more particularly to a kind of magnetocardiogram life being applied in magnetocardiograph
Into method and system.
Background technology
Cardiac muscle cell's periodically polarization forms electrocardiogram (Electrocardiograph, ECG) with process of depolarization,
But the loss of signal is larger after tissue is spread out of due to electrocardiosignal, and the conduction of magnetic signal is not by the shadow of the tissues such as muscle
Ring, therefore magnetocardiogram (magnetocardiogram, MCG) turns into the focus studied recently.
Magnetocardiogram is the magnetic field intensity that different zones above human heart are measured using multichannel magnetocardiograph, and these
The magnetic field value opsition dependent arrangement measured is shown with different colors.Research shows that magnetocardiogram is in coronary heart disease, arrhythmia cordis etc.
Obvious advantage be present in the diagnosis of heart disease, risk stratification etc..Brisinda et al. has found that magnetocardiogram can be used for differentiating
Patients with left bundle branch block (their the R peak positions magnetic field evolution directions of magnetocardiogram are opposite with normal person)[1].Park's et al.
Research shows that, for prediction acute chest pain and the patients with coronary heart disease raised without ST sections, MCG is better than ECG, echocardiogram and flesh calcium
Albumen test；In addition, MCG can accurately detect the bundle-branch block patient with acute chest pain, and ECG can not then be examined
It is disconnected[2].Fainziberg et al. have studied 123 patients with coronary heart disease and 124 normal persons using 4 passage magnetocardiographs, and to by
The magnetocardiogram of examination person is classified (0 normal -4 is completely abnormal), if finding, with classification value 1.75 be threshold value, and MCG classifies quick
Perception and uniformity have respectively reached 76% and 81%[3]。
The studies above shows that magnetocardiogram has potential advantages in diagnosis, prediction, risk assessment of heart disease etc., because
How this, turn into a focus of clinical research by the magnetocardiogram of magnetocardiogram algorithm generation high quality.Chinese patent literature number
CN1552285A discloses (bulletin) day 2004.12.08, discloses a kind of heart magnetic wave group figure and its analysis method, but this method is only
It is that the mcg-signalses that diverse location synchronization measures are drawn in same reference axis, diverse location above heart can not be reflected
Distribution of Magnetic Field situation.Chinese patent literature CN102682425B authorizes (bulletin) day 2014.10.15, discloses a kind of heart
Magnetic chart system and the method for creating magnetocardiogram image, this method use a kind of algorithm based on model learning, i.e., using based on complete
A large amount of high-resolution MCG images that Ao-Sa cuts down that law and randomly generated carry out tectonic model, by being equipped with sparse survey for model
Amount, create high-resolution MCG images.But this method will could produce the heart by creating model, and being equipped with sparseness measuring for model
Magnetic chart image, algorithm is excessively complicated, processing time length, can not meet the requirement for clinically analyzing magnetocardiogram collection of illustrative plates in real time.Therefore,
Field strength values such as how finite point position produce high-resolution, real-time good magnetocardiogram is we need to solve one
Problem.
Bibliography：
[1]D.Brisinda and R.Fenici,"Magnetocardiographic study of patients
with right and left bundle branch blocks,"in International Congress Series,
2007,pp.451-454.
[2]J.W.Park,P.M.Hill,N.Chung,P.G.Hugenholtz,and F.Jung,"
Magnetocardiography predicts coronary artery disease in patients with acute
chest pain,"Annals of noninvasive electrocardiology,vol.10,pp.312-323,2005.
[3]L.Fainzilberg,I.Chaikovsky,S.Auth-Eisernitz,B.Awolin,D.Ivaschenko,
and B.Hailer,"Sensitivity and specificity of magnetocardiography,using
computerized classification of current density vectors maps,in ischemic
patients with normal ECG and echocardiogram,"in International Congress
Series,2007,pp.468-471.
The content of the invention
In view of the above the shortcomings that prior art, it is an object of the invention to provide a kind of heart applied to magnetocardiograph
Magnetic chart generation method and system, for solving above mentioned problem of the prior art, so as to provide ginseng for clinical diagnosis heart disease
Examine value.
In order to achieve the above objects and other related objects, the present invention provides a kind of generation method of magnetocardiogram, including：From the heart
The heart magnetic data of a cycle is intercepted in each passage of magnetic chart instrument respectively, a certain moment is chosen from the heart magnetic data of interception, takes
Go out heart magnetic data of each passage at the moment；Each passage is calculated in position coordinates corresponding to the surface of thoracic cavity；Each passage is made at this
The heart magnetic data at quarter is corresponding with the position coordinates being calculated, and each coordinate opsition dependent after correspondence is arranged to form Two-Dimensional Moment
Battle array；Two-dimensional interpolation processing is carried out to the heart magnetic data of the two-dimensional matrix；Found in color table and the heart magnetic number after each interpolation
Shown according to the corresponding color of value, and in the correspondence position of magnetocardiogram to draw out the magnetocardiogram of colour.
In one embodiment of the invention, the heart magnetic data of each passage a cycle of interception, including：Calculate all passages
The average period of heart magnetic data Wherein, RtiAt the time of for where i-th of R peak, N is
The number of R ripples in heart magnetic data, K are channel number；Using the moment where the R peaks in each passage heart magnetic data as reference point, respectively
The data of fraction times average period are taken forward, the data of 1- fractions times individual average period are taken backward, so as to form a cycle
Heart magnetic data.
It is described to calculate each passage in position coordinates corresponding to the surface of thoracic cavity in one embodiment of the invention, including：Order：Institute
There is passage on an X/Y plane, the coordinate of first passage is (0,0), and the distance of each interchannel is equal and is S；Then：From original
Point (0,0) starts, and the coordinate of first row passage is (0,0+iS), and the coordinate of secondary series passage is (S, 0+iS), the 3rd row passage
Coordinate be (2S, 0+iS) ..., the coordinate of the n-th row passage is ((n-1) S, 0+iS)；Wherein, i is the step number moved along Y-axis.
In one embodiment of the invention, each coordinate opsition dependent by after correspondence is arranged to form two-dimensional matrix, including：
The distribution of its coordinate is expressed as C { (x by the heart magnetic data after the completion of order is corresponding along Y-axis negative sense and the positive arrangement of X-axis1,y1)、(x2,
y2)、(x3,y3)...(xK,yK)}；Search and coordinate (X in C are distributed in coordinatei, Yi) passage corresponding to equal coordinate, its
In, (Xi, Yi) it is arrangement coordinates of the passage i on X/Y plane；Search circulation is jumped out when finding corresponding coordinate, and by passage i
Heart magnetic data MCGiK-th of position of one-dimension array is put into, wherein, k is the execution number of search circulation；Repeat, directly
It is put into by all heart magnetic datas corresponding to one-dimension array on position；N data are put before being taken out from the one-dimension array of acquisition
In the first row, further take out (n+1)th to the 2n data and be placed on the second row, by that analogy, until by one-dimension array by Y-axis negative sense and
X-axis forward direction is arranged in two-dimensional matrix.
In one embodiment of the invention, the heart magnetic data to two-dimensional matrix carries out two-dimensional interpolation processing, including：Order：
A series of coordinate points (x on X/Y planei, yj) (i=0,1,2 ..., n；J=0,1,2 ..., m) on heart magnetic data value be
aij=f (xi,yj)；Binary function z=z (x, y) is constructed, makes z (xi,yj) value approach aij, now, z (xi,yj) it is required
Coordinate points (xi, yj) place interpolation.
In one embodiment of the invention, the binary function z of construction is：Wherein,For X-axis coordinate section { x0＜ x1＜ x2＜ ... ＜ xnOn non-zero Gene radiotherapy three times, ψj(y) it is Y-axis coordinate area
Between { y0＜ y1＜ y2＜ ... ＜ ymOn non-zero Gene radiotherapy three times, uijFor undetermined coefficient, by being obtained to z (x, y) derivation
Arrive.
In one embodiment of the invention, the magnetocardiogram for drawing colour, in addition to：Find out maximum in image and most
Small value, drawn respectively at maximum+number, drawn at minimum value-, and maximum has been connected with minimum value with straight line with the arrow
Come, wherein, arrow points to minimum value.
In order to achieve the above objects and other related objects, the present invention provides a kind of generation system of magnetocardiogram, including：Data
Acquisition module, for intercepting the heart magnetic data of a cycle respectively from each passage of magnetocardiograph, from the heart magnetic data of interception
It is middle to choose a certain moment, take out heart magnetic data of each passage at the moment；Coordinate calculation module, for calculating each passage in thoracic cavity
Position coordinates corresponding to surface；Alignment processing module, for making heart magnetic data and the position that is calculated of each passage at the moment
Put coordinate pair to answer, and each coordinate opsition dependent after correspondence is arranged to form two-dimensional matrix；Interpolation processing module, for described
The heart magnetic data of two-dimensional matrix carries out two-dimensional interpolation processing；Image generation module, for found in color table with after each interpolation
Heart magnetic data be worth corresponding color, and shown in the correspondence position of magnetocardiogram to draw out the magnetocardiogram of colour.
In one embodiment of the invention, the heart magnetic data of each passage a cycle of interception, including：Calculate all passages
The average period of heart magnetic data Wherein, RtiAt the time of for where i-th of R peak, N is
The number of R ripples in heart magnetic data, K are channel number；Using the moment where the R peaks in each passage heart magnetic data as reference point, respectively
The data of fraction times average period are taken forward, the data of 1- fractions times individual average period are taken backward, so as to form a cycle
Heart magnetic data.
It is described to calculate each passage in position coordinates corresponding to the surface of thoracic cavity in one embodiment of the invention, including：Order：Institute
There is passage on an X/Y plane, the coordinate of first passage is (0,0), and the distance of each interchannel is equal and is S；Then：From original
Point (0,0) starts, and the coordinate of first row passage is (0,0+iS), and the coordinate of secondary series passage is (S, 0+iS), the 3rd row passage
Coordinate be (2S, 0+iS) ..., the coordinate of the n-th row passage is ((n-1) S, 0+iS)；Wherein, i is the step number moved along Y-axis.
In one embodiment of the invention, each coordinate opsition dependent by after correspondence is arranged to form two-dimensional matrix, including：
The distribution of its coordinate is expressed as C { (x by the heart magnetic data after the completion of order is corresponding along Y-axis negative sense and the positive arrangement of X-axis1,y1)、(x2,
y2)、(x3,y3)...(xK,yK)}；Search and coordinate (X in C are distributed in coordinatei, Yi) passage corresponding to equal coordinate, its
In, (Xi, Yi) it is arrangement coordinates of the passage i on X/Y plane；Search circulation is jumped out when finding corresponding coordinate, and by passage i
Heart magnetic data MCGiK-th of position of one-dimension array is put into, wherein, k is the execution number of search circulation；Repeat, directly
It is put into by all heart magnetic datas corresponding to one-dimension array on position；N data are put before being taken out from the one-dimension array of acquisition
In the first row, further take out (n+1)th to the 2n data and be placed on the second row, by that analogy, until by one-dimension array by Y-axis negative sense and
X-axis forward direction is arranged in two-dimensional matrix.
In one embodiment of the invention, the heart magnetic data to two-dimensional matrix carries out two-dimensional interpolation processing, including：Order：
A series of coordinate points (x on X/Y planei, yj) (i=0,1,2 ..., n；J=0,1,2 ..., m) on heart magnetic data value be
aij=f (xi,yj)；Binary function z=z (x, y) is constructed, makes z (xi,yj) value approach aij, now, z (xi,yj) it is required
Coordinate points (xi, yj) place interpolation.
In one embodiment of the invention, the binary function z (x, y) of construction is：Its
In,For X-axis coordinate section { x0＜ x1＜ x2＜ ... ＜ xnOn non-zero Gene radiotherapy three times, ψj(y) it is Y-axis coordinate
Section { y0＜ y1＜ y2＜ ... ＜ ymOn non-zero Gene radiotherapy three times, uijFor undetermined coefficient, by z (x, y) derivation
Obtain.
In one embodiment of the invention, described image generation module is additionally operable to：The maximum and minimum value in image are found out,
Drawn respectively at maximum+number, drawn at minimum value-, and maximum and minimum value are connected with straight line with the arrow, its
In, arrow points to minimum value.
As described above, the generation method and generation system of the magnetocardiogram of the present invention, the magnetocardiogram image clearly of generation, color
It is evenly distributed ,+- two dipoles can be clearly visible that in T wave bands.If draw difference from the S sections of QRS wave to T sections by setting step-length
The magnetocardiogram at moment, then it can be seen that the temporal-spatial evolution of cardiomotility electric current, doctor can according to the dipole number of magnetocardiogram and
The evolution direction of cardiomotility electric current diagnoses whether the patient suffers from heart disease.
Brief description of the drawings
Fig. 1 is shown as the flow chart of the generation method of the magnetocardiogram in one embodiment of the invention.
Fig. 2 is shown as during step-length i=5 in one embodiment of the invention K passage in the coordinate distribution situation of X/Y plane.
Fig. 3 is shown as the schematic diagram of the interpolation algorithm of the heart magnetic data in one embodiment of the invention.
Fig. 4 is shown as the magnetocardiogram of the normal person generated with the magnetocardiogram algorithm of the present invention.
Fig. 5 is shown as the magnetocardiogram of the cardiac generated with the magnetocardiogram algorithm of the present invention.
Fig. 6 is shown as the module map of the generation system of the magnetocardiogram in one embodiment of the invention.
Embodiment
Illustrate embodiments of the present invention below by way of specific instantiation, those skilled in the art can be by this specification
Disclosed content understands other advantages and effect of the present invention easily.The present invention can also pass through specific realities different in addition
The mode of applying is embodied or practiced, the various details in this specification can also be based on different viewpoints with application, without departing from
Various modifications or alterations are carried out under the spirit of the present invention.It should be noted that in the case where not conflicting, following examples and implementation
Feature in example can be mutually combined.
It should be noted that the diagram provided in following examples only illustrates the basic structure of the present invention in a schematic way
Think, only show the component relevant with the present invention in schema then rather than according to component count, shape and the size during actual implement
Draw, kenel, quantity and the ratio of each component can be a kind of random change during its actual implementation, and its assembly layout kenel
It is likely more complexity.
As shown in figure 1, the magnetocardiogram generation method of the present invention, specifically includes 5 steps：
Step S1：Intercept the heart magnetic data of each passage sometime.Concrete methods of realizing is：First, intercept each logical
(number of the heart magnetic data for a cycle taken out from each passage is heart magnetic data MCG1~MCGM of road a cycle
M, for 36 passage heart magnetic acquisition systems, total number of channels 36, the data of each passage a cycle probably have more than 800),
And a certain moment is chosen from the data of interception, take out the data at each passage moment；Then, each passage is calculated in chest table
The position coordinates of measurement point corresponding to face.
The data of the described each passage a cycle of interception, including：First, being averaged for all passage heart magnetic datas is calculated
CycleThen, using the moment where R peaks in each passage heart magnetic data as reference point, 2/5 average period is respectively forwardly taken
Data, the data of 3/5 average period are taken backward, form a cycle.
Described calculation formula average period is,Wherein, RtiFor i-th of R peaks institute
At the time of, N be heart magnetic data in R ripples number, K is total channel number.
Each passage of described calculating corresponding position coordinates above chest：Assuming that all passages (K) are put down in an XY
On face, the coordinate of first passage is (0,0), and the distance of each interchannel is equal and is S, then along X-direction since origin
The coordinate of first row passage is (0,0+iS), and i is the step number moved up along Y direction, and i value is whole in the range of 1~K
Number；The coordinate of secondary series passage is (S, 0+iS), by that analogy, can calculate X/Y plane coordinate corresponding to other all passages.Fig. 2
When showing step-length i=5, the coordinate distribution situation of K passage.
Step S2：The heart magnetic data of each passage sometime is allowed to be answered with coordinate pair, i.e. the heart magnetic each passage sometime
Data are put on corresponding position coordinates, form two-dimensional matrix.Concrete methods of realizing is：Data after the completion of false coordinate is corresponding
Arranged along Y-axis is positive, X-axis is positive, its coordinate is distributed as C { (x1,y1)、(x2,y2)、(x3,y3)...(xK,yK)}；Assuming that passage i
Arrangement coordinate on XY two dimensional surfaces is (Xi, Yi).First, it is distributed in C and finds and coordinate (X in coordinatei, Yi) equal logical
Road, when finding corresponding coordinate, jump out search circulation, it is assumed that now cycle-index is k, then passage i heart magnetic data
MCGi is put into k-th of position of one-dimension array；Then, all heart magnetic datas are put into corresponding to one-dimension array as stated above
On position；Such as：36 passage heart magnetic systems, each passage of synchronization take with regard to each data for taking out this moment, i.e., each passage
Go out a data, take out 36 data altogether, found and each data coordinates (X on Ci,Yi) duplicate coordinate, search out
Circulation is exited, cycle-index is corresponding location number 1...36 in 36 positions；Finally, above-mentioned one-dimension array is pressed Y-axis
Negative sense, X-axis forward direction are arranged in two-dimensional matrix, that is, assume that two-dimensional matrix columns is n, then take out the preceding n number of one-dimension array first
According to the first row is placed on, (n+1)th to the 2n data for then taking out one-dimension array are placed on the second row, by that analogy, all numbers
According to the relevant position for being placed on two-dimensional matrix.
Step S3：The two-dimentional heart magnetic data arranged to opsition dependent carries out two-dimensional interpolation.Because the signal of magnetocardiograph leads to
Road number is limited, such as：36 passage magnetocardiograph a certain moment only had 36 heart magnetic datas, and heart magnetic is carried out with this 36 heart magnetic datas
Figure imaging, image resolution ratio is too low, therefore, it is necessary to the two-dimentional heart magnetic data arranged to opsition dependent carries out two-dimensional interpolation.
Referring to Fig. 3, the present invention enters row interpolation using bicubic spline interpolation algorithm to two-dimentional heart magnetic data, it is embodied
Method is：A series of assuming that coordinate points (x on X/Y planei, yj) (i=0,1,2 ..., n；J=0,1,2 ..., m；N's and m
Value is the integer in the range of 1~K) on heart magnetic data value be aij=f (xi,yj), construction binary function z=z (x, y), make z
(xi,yj) value approach aij, i.e. z (xi,yj) it is required coordinate points (xi, yj) place interpolation.Fig. 3 shows adopting before interpolation
Collect data zi,j(xi, yj) and interpolated data zim,n(xim,n, yim,n) distribution situation.
Constructed in the present inventionWherein,For X-axis coordinate section { x0＜ x1＜ x2
＜ ... ＜ xnOn non-zero Gene radiotherapy three times, ψj(y) it is Y-axis coordinate section { y0＜ y1＜ y2＜ ... ＜ ymOn it is non-
Zero Gene radiotherapy three times, uij, can be by being obtained to z (x, y) derived function for undetermined coefficient.
Step S4：Color table is inquired about to the 2-D data after interpolation.According to the heart magnetic data value after interpolation in color table
Corresponding color is found to show in magnetocardiogram correspondence position.To improve the image resolution ratio of magnetocardiogram, color table uses for
Circulation 256 colo(u)r streak uniform color tables of generation, the difference of each color can sets itself, to reach preferable uniformity and continuous
Property.
Step S5：Draw magnetocardiogram.Concrete methods of realizing is：First, the size of image is set；Then, in the correspondence of image
The value that heart magnetic data is pressed on position searches color table, and corresponding color value is drawn on image；Finally, the maximum in image is found out
Value and minimum value, drawn respectively at maximum+number, drawn at minimum value-, and maximum and minimum value are connected with straight line with the arrow
Pick up and, arrow points to minimum value.
To verify the validity of magnetocardiogram generating algorithm, we are using 9 passage magnetocardiographs in Xuhui District central hospital point
The heart magnetic data for not acquiring normal person and cardiac carries out proof of algorithm.First, using 9 passage magnetocardiographs in thoracic cavity
Four fixing point station acquisitions, four heart magnetic datas (data of 9 passages of collection every time) of top setting, it is logical to form 36 altogether
Road heart magnetic data；Then, according to the magnetocardiogram generating algorithm of the present invention, 1. the heart magnetic data of this 36 passage is cut successively
The data of each passage a cycle are taken, and a certain moment is chosen from the data of interception, take out the number at each passage moment
According to；Then each passage corresponding position coordinates above chest is calculated.2. allow the data of each passage sometime should with coordinate pair.
3. the two-dimentional heart magnetic data arranged to opsition dependent carries out two-dimensional interpolation.4. color table is inquired about to the 2-D data after interpolation.⑤
Draw 5 steps such as magnetocardiogram, normal person's magnetocardiogram of generation as shown in Figure 4, the magnetocardiogram such as institute of accompanying drawing 5 of cardiac
Show.
As shown in fig. 6, with above method embodiment principle similarly, the present invention also provides a kind of generation system of magnetocardiogram
System, including：Data acquisition module 601, coordinate calculation module 602, alignment processing module 603, interpolation processing module 604, image
Generation module 605.Because the technical characteristic in preceding method embodiment can apply to the system embodiment, thus it is not repeated
Repeat.
Data acquisition module 604 intercepts the heart magnetic data of a cycle respectively from each passage of magnetocardiograph, from interception
Heart magnetic data in choose a certain moment, take out heart magnetic data of each passage at the moment.Coordinate calculation module 602 calculates each logical
Road is in position coordinates corresponding to the surface of thoracic cavity.Alignment processing module 603 make each passage the moment heart magnetic data and calculate
The position coordinates arrived is corresponding, and each coordinate opsition dependent after correspondence is arranged to form two-dimensional matrix.Interpolation processing module 604 is right
The heart magnetic data of the two-dimensional matrix carries out two-dimensional interpolation processing.Image generation module 605 is found and each interpolation in color table
Heart magnetic data afterwards is worth corresponding color, and is shown in the correspondence position of magnetocardiogram to draw out the magnetocardiogram of colour.
In summary, magnetocardiogram generation method of the invention and generation system, effectively overcome of the prior art a variety of
Shortcoming and have high industrial utilization.
The above-described embodiments merely illustrate the principles and effects of the present invention, not for the limitation present invention.It is any ripe
Know the personage of this technology all can carry out modifications and changes under the spirit and scope without prejudice to the present invention to above-described embodiment.Cause
This, those of ordinary skill in the art is complete without departing from disclosed spirit and institute under technological thought such as
Into all equivalent modifications or change, should by the present invention claim be covered.
Claims (14)
- A kind of 1. generation method of magnetocardiogram, it is characterised in that including：Intercept the heart magnetic data of a cycle respectively from each passage of magnetocardiograph, chosen from the heart magnetic data of interception a certain At the moment, take out heart magnetic data of each passage at the moment；Each passage is calculated in position coordinates corresponding to the surface of thoracic cavity；Make each passage corresponding with the position coordinates being calculated in the heart magnetic data at the moment, and by each coordinate step-by-step after correspondence Arrangement is put to form two-dimensional matrix；Two-dimensional interpolation processing is carried out to the heart magnetic data of the two-dimensional matrix；The color corresponding with the heart magnetic data value after each interpolation is found in color table, and is shown in the correspondence position of magnetocardiogram Out with draw out colour magnetocardiogram.
- 2. according to the method for claim 1, it is characterised in that the heart magnetic data of each passage a cycle of interception, bag Include：Calculate the average period of all passage heart magnetic datas Wherein, RtiFor i-th of R At the time of where peak, N is the number of R ripples in heart magnetic data, and K is channel number；Using the moment where the R peaks in each passage heart magnetic data as reference point, the data of fraction times average period are respectively forwardly taken, The data of 1- fractions times individual average period are taken backward, so as to form the heart magnetic data of a cycle.
- 3. according to the method for claim 1, it is characterised in that each passage of calculating is sat in position corresponding to the surface of thoracic cavity Mark, including：Order：All passages are on an X/Y plane, and the coordinate of first passage is (0,0), and the distance of each interchannel is equal and is S；Then：Starting from origin (0,0), the coordinate of first row passage is (0,0+iS), and the coordinate of secondary series passage is (S, 0+iS), The coordinate of 3rd row passage is (2S, 0+iS) ..., and the coordinate of the n-th row passage is ((n-1) S, 0+iS)；Wherein, i is to be moved along Y-axis Dynamic step number.
- 4. according to the method for claim 3, it is characterised in that each coordinate opsition dependent by after correspondence is arranged to be formed Two-dimensional matrix, including：The distribution of its coordinate is expressed as C { (x by the heart magnetic data after the completion of order is corresponding along Y-axis forward direction and the positive arrangement of X-axis1,y1)、 (x2,y2)、(x3,y3)...(xK,yK)}；Search and coordinate (X in C are distributed in coordinatei, Yi) passage corresponding to equal coordinate, wherein, (Xi, Yi) exist for passage i Arrangement coordinate on X/Y plane；Search circulation is jumped out when finding corresponding coordinate, and by passage i heart magnetic data MCGiIt is put into K-th of position of one-dimension array, wherein, k is the execution number of search circulation；Repeat, until by all heart magnetic datas It is put into corresponding to one-dimension array on position；N data are placed on the first row before being taken out from the one-dimension array of acquisition, further take out (n+1)th to the 2n data and are placed on second OK, by that analogy, until one-dimension array is arranged in into two-dimensional matrix by Y-axis negative sense and X-axis forward direction.
- 5. according to the method for claim 1, it is characterised in that the heart magnetic data to two-dimensional matrix carries out two-dimensional interpolation Processing, including：Order：A series of coordinate points (x on X/Y planei, yj) (i=0,1,2 ..., n；J=0,1,2 ..., m) on heart magnetic number It is a according to valueij=f (xi,yj)；Binary function z=z (x, y) is constructed, makes z (xi,yj) value approach aij, now, z (xi,yj) it is required coordinate points (xi, yj) place interpolation.
- 6. according to the method for claim 5, it is characterised in that the binary function z of construction is：Wherein,For X-axis coordinate section { x0＜ x1＜ x2＜ ... ＜ xnOn non-zero Gene radiotherapy three times, ψj(y) it is Y Axial coordinate section { y0＜ y1＜ y2＜ ... ＜ ymOn non-zero Gene radiotherapy three times, uijFor undetermined coefficient, by z (x, Y) derivation obtains.
- 7. according to the method for claim 1, it is characterised in that the magnetocardiogram for drawing colour, in addition to：Find out image In maximum and minimum value, drawn respectively at maximum+number, drawn at minimum value-, and maximum and minimum value with the arrow Straight line connect, wherein, arrow point to minimum value.
- A kind of 8. generation system of magnetocardiogram, it is characterised in that including：Data acquisition module, for intercepting the heart magnetic data of a cycle respectively from each passage of magnetocardiograph, from interception The a certain moment is chosen in heart magnetic data, takes out heart magnetic data of each passage at the moment；Coordinate calculation module, for calculating each passage in position coordinates corresponding to the surface of thoracic cavity；Alignment processing module, for making each passage corresponding with the position coordinates being calculated in the heart magnetic data at the moment, and will Each coordinate opsition dependent after corresponding is arranged to form two-dimensional matrix；Interpolation processing module, for carrying out two-dimensional interpolation processing to the heart magnetic data of the two-dimensional matrix；Image generation module, for finding the color corresponding with the heart magnetic data value after each interpolation in color table, and in the heart The correspondence position of magnetic chart is shown to draw out the magnetocardiogram of colour.
- 9. system according to claim 8, it is characterised in that the heart magnetic data of each passage a cycle of interception, bag Include：Calculate the average period of all passage heart magnetic datas Wherein, RtiFor i-th of R At the time of where peak, N is the number of R ripples in heart magnetic data, and K is channel number；Using the moment where the R peaks in each passage heart magnetic data as reference point, the data of fraction times average period are respectively forwardly taken, The data of 1- fractions times individual average period are taken backward, so as to form the heart magnetic data of a cycle.
- 10. system according to claim 8, it is characterised in that described to calculate each passage in position corresponding to the surface of thoracic cavity Coordinate, including：Order：All passages are on an X/Y plane, and the coordinate of first passage is (0,0), and the distance of each interchannel is equal and is S；Then：Starting from origin (0,0), the coordinate of first row passage is (0,0+iS), and the coordinate of secondary series passage is (S, 0+iS), The coordinate of 3rd row passage is (2S, 0+iS) ..., and the coordinate of the n-th row passage is ((n-1) S, 0+iS)；Wherein, i is to be moved along Y-axis Dynamic step number.
- 11. system according to claim 10, it is characterised in that each coordinate opsition dependent by after correspondence is arranged with shape Into two-dimensional matrix, including：The distribution of its coordinate is expressed as C { (x by the heart magnetic data after the completion of order is corresponding along Y-axis forward direction and the positive arrangement of X-axis1,y1)、 (x2,y2)、(x3,y3)...(xK,yK)}；Search and coordinate (X in C are distributed in coordinatei, Yi) passage corresponding to equal coordinate, wherein, (Xi, Yi) exist for passage i Arrangement coordinate on X/Y plane；Search circulation is jumped out when finding corresponding coordinate, and by passage i heart magnetic data MCGiIt is put into K-th of position of one-dimension array, wherein, k is the execution number of search circulation；Repeat, until by all heart magnetic datas It is put into corresponding to one-dimension array on position；N data are placed on the first row before being taken out from the one-dimension array of acquisition, further take out (n+1)th to the 2n data and are placed on second OK, by that analogy, until one-dimension array is arranged in into two-dimensional matrix by Y-axis negative sense and X-axis forward direction.
- 12. system according to claim 8, it is characterised in that the heart magnetic data to two-dimensional matrix carries out two dimension and inserted Value processing, including：Order：A series of coordinate points (x on X/Y planei, yj) (i=0,1,2 ..., n；J=0,1,2 ..., m) on heart magnetic number It is a according to valueij=f (xi,yj)；Binary function z=z (x, y) is constructed, makes z (xi,yj) value approach aij, now, z (xi,yj) it is required coordinate points (xi, yj) place interpolation.
- 13. system according to claim 12, it is characterised in that the binary function z (x, y) of construction is：Wherein,For X-axis coordinate section { x0＜ x1＜ x2＜ ... ＜ xnOn non-zero Gene radiotherapy three times, ψj(y) it is Y Axial coordinate section { y0＜ y1＜ y2＜ ... ＜ ymOn non-zero Gene radiotherapy three times, uijFor undetermined coefficient, by z (x, Y) derivation obtains.
- 14. system according to claim 8, it is characterised in that described image generation module is additionally operable to：Find out in image Maximum and minimum value, drawn respectively at maximum+number, drawn at minimum value-, and maximum and minimum value with the arrow straight Line connects, wherein, arrow points to minimum value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710564362.7A CN107374610B (en) | 2017-07-12 | 2017-07-12 | Magnetocardiogram generation method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710564362.7A CN107374610B (en) | 2017-07-12 | 2017-07-12 | Magnetocardiogram generation method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107374610A true CN107374610A (en) | 2017-11-24 |
| CN107374610B CN107374610B (en) | 2020-05-08 |
Family
ID=60339073
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710564362.7A Active CN107374610B (en) | 2017-07-12 | 2017-07-12 | Magnetocardiogram generation method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107374610B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109998519A (en) * | 2019-04-19 | 2019-07-12 | 北京航空航天大学 | A kind of magnetocardiogram measurement and magnetocardiogram generation system based on SERF atom magnetometer |
| CN116965822A (en) * | 2023-09-25 | 2023-10-31 | 合肥工业大学 | Magnetocardiogram color space circle diagram generation and wave band time identification method and storage medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5211178A (en) * | 1990-09-30 | 1993-05-18 | Agency Of Industrial Science & Technology | Method and apparatus for synchronously detecting a magnetocardiogram and method and apparatus for synchronously adding magnetocardiograms |
| DE602005002326T2 (en) * | 2004-09-29 | 2008-05-29 | Hitachi High-Technologies Corp. | Biomagnetic measuring device |
| CN101836860A (en) * | 2010-06-10 | 2010-09-22 | 北京大学 | Method for determining point positions of maximum value and minimum value of cardiac magnetic field |
| CN104825133A (en) * | 2015-05-04 | 2015-08-12 | 河南理工大学 | Colored Doppler 3D (three-dimensional) imaging based quasistatic ventricle-heart magnetic field model |
-
2017
- 2017-07-12 CN CN201710564362.7A patent/CN107374610B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5211178A (en) * | 1990-09-30 | 1993-05-18 | Agency Of Industrial Science & Technology | Method and apparatus for synchronously detecting a magnetocardiogram and method and apparatus for synchronously adding magnetocardiograms |
| DE602005002326T2 (en) * | 2004-09-29 | 2008-05-29 | Hitachi High-Technologies Corp. | Biomagnetic measuring device |
| CN101836860A (en) * | 2010-06-10 | 2010-09-22 | 北京大学 | Method for determining point positions of maximum value and minimum value of cardiac magnetic field |
| CN104825133A (en) * | 2015-05-04 | 2015-08-12 | 河南理工大学 | Colored Doppler 3D (three-dimensional) imaging based quasistatic ventricle-heart magnetic field model |
Non-Patent Citations (2)
| Title |
|---|
| 宋爱斌等: "双立方插值和二元三次样条插值在心磁图像处理中的应用", 《数学的实践与认识》 * |
| 李明等: "９通道心磁图仪性能测试与临床应用研究", 《低温物理学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109998519A (en) * | 2019-04-19 | 2019-07-12 | 北京航空航天大学 | A kind of magnetocardiogram measurement and magnetocardiogram generation system based on SERF atom magnetometer |
| CN116965822A (en) * | 2023-09-25 | 2023-10-31 | 合肥工业大学 | Magnetocardiogram color space circle diagram generation and wave band time identification method and storage medium |
| CN116965822B (en) * | 2023-09-25 | 2023-12-29 | 合肥工业大学 | Magnetocardiogram color space circle diagram generation and wave band time identification method and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107374610B (en) | 2020-05-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9125581B2 (en) | Continuous modeling for dipole localization from 2D MCG images with unknown depth | |
| EP2020914B1 (en) | Methods and apparatus of three dimensional cardiac electrophysiological imaging | |
| Disbrow et al. | Evidence for interhemispheric processing of inputs from the hands in human S2 and PV | |
| Li et al. | Localization of the site of origin of cardiac activation by means of a heart-model-based electrocardiographic imaging approach | |
| He et al. | Body surface Laplacian ECG mapping | |
| Brookes et al. | Measuring functional connectivity using MEG: methodology and comparison with fcMRI | |
| Habermehl et al. | Somatosensory activation of two fingers can be discriminated with ultrahigh-density diffuse optical tomography | |
| US8553956B2 (en) | 3D current reconstruction from 2D dense MCG images | |
| US5782773A (en) | Three-dimensional electrocardiogram display method | |
| Diwakar et al. | Dual-core beamformer for obtaining highly correlated neuronal networks in MEG | |
| Numminen et al. | Transformation of multichannel magnetocardiographic signals to standard grid form | |
| JP2012157696A (en) | Magnetocardiogram system, and method for creating magnetocardiogram image | |
| CN106073753A (en) | 18 lead holographic dynamic and static ecg analysis method and system | |
| CN105212895A (en) | Dynamic brain source localization method | |
| JP2013066708A (en) | Noise elimination of mcg measurement | |
| CN107374610A (en) | Magnetocardiogram generation method and generation system | |
| CN106901717A (en) | Three-dimensional heart source localization method, system and server | |
| Loeffler et al. | Reaction-diffusion informed approach to determine myocardial ischemia using stochastic in-silico ECGs and CNNs | |
| He et al. | A practical method for quickly determining electrode positions in high-density EEG studies | |
| Dossel et al. | Optimization of electrode positions for multichannel electrocardiography with respect to electrical imaging of the heart | |
| Popescu et al. | Reconstruction of fetal cardiac vectors from multichannel fMCG data using recursively applied and projected multiple signal classification | |
| Marhl et al. | Application of source localization algorithms in magnetoencephalography: test on a new generation of magnetometers | |
| Vilchez-Monge et al. | Image reconstruction of the human foearm by electrical impedance tomography | |
| Tsutada et al. | Detecting functional asymmetries through the dipole moment of magnetoencephalography | |
| Bowman et al. | Modelling of cardiac imaging data with spatial correlation |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180911 Address after: 201108 221, room 1, 3688 Kim do road, Minhang District, Shanghai. Applicant after: Man Di medical instruments (Shanghai) Co., Ltd. Address before: 200050 865 Changning Road, Changning District, Changning District, Shanghai. Applicant before: Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |