CN116965822B - Magnetocardiogram color space circle diagram generation and wave band time identification method and storage medium - Google Patents
Magnetocardiogram color space circle diagram generation and wave band time identification method and storage medium Download PDFInfo
- Publication number
- CN116965822B CN116965822B CN202311234687.0A CN202311234687A CN116965822B CN 116965822 B CN116965822 B CN 116965822B CN 202311234687 A CN202311234687 A CN 202311234687A CN 116965822 B CN116965822 B CN 116965822B
- Authority
- CN
- China
- Prior art keywords
- magnetocardiogram
- color space
- data
- circle
- color
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000003860 storage Methods 0.000 title claims abstract description 9
- 238000010586 diagram Methods 0.000 title claims description 60
- 239000000523 sample Substances 0.000 claims abstract description 31
- 238000003708 edge detection Methods 0.000 claims abstract description 14
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 13
- 241000238366 Cephalopoda Species 0.000 claims abstract description 12
- 238000001514 detection method Methods 0.000 claims abstract description 12
- 238000005520 cutting process Methods 0.000 claims abstract description 10
- 238000010606 normalization Methods 0.000 claims abstract description 10
- 230000009467 reduction Effects 0.000 claims abstract description 9
- 230000000007 visual effect Effects 0.000 claims abstract description 7
- 108010076504 Protein Sorting Signals Proteins 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 30
- 238000012545 processing Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 230000000747 cardiac effect Effects 0.000 claims description 8
- 238000003672 processing method Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 5
- 101100391182 Dictyostelium discoideum forI gene Proteins 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 239000000725 suspension Substances 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 230000003628 erosive effect Effects 0.000 claims 1
- 238000006243 chemical reaction Methods 0.000 abstract description 5
- 230000008859 change Effects 0.000 description 12
- 230000006870 function Effects 0.000 description 12
- 238000004891 communication Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000000694 effects Effects 0.000 description 7
- 230000009471 action Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000036541 health Effects 0.000 description 3
- 230000007797 corrosion Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 210000000170 cell membrane Anatomy 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 208000019622 heart disease Diseases 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 150000002500 ions Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000008807 pathological lesion Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000008733 trauma Effects 0.000 description 1
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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
- A61B5/7285—Specific aspects of physiological measurement analysis for synchronising or triggering a physiological measurement or image acquisition with a physiological event or waveform, e.g. an ECG signal
- A61B5/7289—Retrospective gating, i.e. associating measured signals or images with a physiological event after the actual measurement or image acquisition, e.g. by simultaneously recording an additional physiological signal during the measurement or image acquisition
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Signal Processing (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The invention relates to a magnetocardiogram color space circle graph generation and wave band time identification method and a storage medium, wherein the method comprises the steps of obtaining 36-channel magnetocardiogram signals on the surface of a human body based on SQUID probe measurement, designing a window function, intercepting to obtain a complete magnetocardiogram signal period sample, and then carrying out R peak alignment operation and integral normalization operation to realize pretreatment of data; drawing concentric rings to construct a circle graph of the preprocessed signals, calculating an angle value according to the length of a magnetocardiogram signal sequence, carrying out sector cutting on the concentric rings in the circle graph, and carrying out color filling on sectors at corresponding positions based on magnetocardiogram signal data values to realize 36-channel magnetocardiogram signal magnetocardiogram color space circle graph conversion; and performing binarization and noise reduction treatment on the magnetocardiogram color space circular graph, performing edge detection on the treated image, performing Hough transformation linear detection, performing visual display on the linear, and calculating the P-band and QRS-band time of the magnetocardiogram signal according to the linear angle to realize the identification of the band time.
Description
Technical Field
The invention relates to the technical field of magnetocardiogram signal image generation and magnetocardiogram signal wave band time identification, in particular to a magnetocardiogram color space circle graph generation and wave band time identification method and a storage medium.
Background
The heart of a human body can electrically excite myocardial cells during the activity process, so that weak current is generated to form a potential difference on the cell membrane of the heart of the human body. The change of the electric signal can generate corresponding magnetic field change, namely a magnetocardiogram signal, in the body surface space, and the magnetocardiogram signal can show almost periodic transformation due to systole and comfort. The magnetocardiogram signal is mainly formed by collecting magnetic field data of the heart through a superconducting quantum interferometer, and the non-contact measurement mode without external excitation can effectively reduce the trauma and risk caused by detection. Meanwhile, most pathological lesions of the heart are related to the change of the transmembrane ions of the heart, and compared with an Electrocardiosignal (ECG) measured by adopting a potential method, the electrocardiosignal is sensitive to the change, so that abnormal activity of the heart can be found earlier. It contains a great deal of regular, physiological and pathological information of heart activity, and has important significance for clinical diagnosis of heart diseases.
Magnetocardiography (MCG) is an effective method for recording and displaying human magnetocardiography signals, which displays the collected magnetocardiography signals of a superconducting quantum interferometer in different colors according to spatial position arrangement. A doctor can acquire the distribution condition of the magnetic field of the human heart at a certain moment through the color change of the magnetocardiogram, and if a plurality of magnetocardiograms are played in time sequence, the magnetic field change trend of the human heart can be seen. However, one magnetocardiogram can only reflect the respective regularity of the human heart magnetic field at a certain moment, that is, only the spatial distribution of the heart magnetic field at a certain moment can be observed on one magnetocardiogram, and the time change trend at the corresponding position cannot be observed on one magnetocardiogram. How to effectively display the spatial distribution and time variation trend of the cardiac magnetic field on one picture at the same time is a problem to be solved. In addition, after the trend of time transformation is displayed, how to identify the Q band and QRS band in the time transformation process is also a problem to be solved together. The magnetocardiogram color space circle diagram of the invention corresponds to the surface of a human body, 36-channel magnetocardiogram signals are obtained based on SQUID probe measurement, the time variation trend of single-channel magnetocardiogram signals is displayed through the sector on a single circle on the circle diagram, and the spatial position up-conversion trend of 36-channel magnetocardiogram signals is displayed through the sector between the circles on the circle diagram. A doctor can observe the space information and the time information of the 36-channel magnetocardiogram signals through the magnetocardiogram color space circular graph, and simultaneously obtain the activity condition of the heart of the patient in one period according to the calculated P wave band and QRS wave band time, so as to help the doctor to judge the health state of the heart of the patient.
Disclosure of Invention
The magnetocardiogram color space circle graph generation and wave band time identification method provided by the invention can at least solve one of the technical problems in the background technology.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
a characteristic image generating and recognizing method of electrocardiogram data includes the following steps,
step A: human body surface 36 channel magnetocardiogram signal set obtained based on SQUID probe measurementThe unit is millivolts; />The time is in seconds, a window function is designed to intercept to obtain a complete magnetocardiogram signal period sample, and R peak alignment operation and integral normalization operation are carried out on the intercepted magnetocardiogram signal to obtain a normalized magnetocardiogram signal set->Preprocessing of data is achieved;
and (B) step (B): for normalized magnetocardiogram signal setDrawing concentric rings to construct a circle graph according to normalized magnetocardiogram signal set +.>Cutting concentric circles in the circle graph into sectors to give sectors and +.>The position relation of the data is based on +.>Performing color filling operation on the data to obtain a magnetocardiogram color space circle diagram>；
Step C: centering magnetic color space circle diagramDrawing an edge magnetocardiogram space circle graph for binarization and noise reduction processing and edge detection of the processed image>Then->And carrying out Hough transformation linear detection, and calculating the time of the P wave band and the QRS wave band of the magnetocardiogram signal according to the linear angle value.
Further said step a: for 36-channel magnetocardiogram signal set obtained based on SQUID probe measurementThe unit is millivolts; />The time is in seconds, a window function is designed to intercept to obtain a complete magnetocardiogram signal period sample, and R peak alignment operation and integral normalization operation are carried out on the intercepted magnetocardiogram signal to obtain a normalized magnetocardiogram signal set->The preprocessing of the data is realized, which comprises the following specific steps:
s101: 36-channel magnetocardiogram signal set based on SQUID probe measurementExpressed in the form of a sequence, i.e. 36 channel magnetocardiogram signal set +.>Designing a rectangular window function to intercept to obtain a complete magnetocardiogram signal period sample +.>；
Wherein,for the sampling amount +.>Is 36-channel magnetocardiogram signal set->Is a sampling frequency of (a);for the window function +.>、/>Respectively a window function start point and an end point, which correspond to the pair +.>When intercepting each line of (a)>、/>Are all different in value; />Is 36-channel magnetocardiogram signal set sequence, the unit is millivolt, is 36-channel magnetocardiogram signal set +.>In the form of a sequence representation of (a),it is +.>A two-dimensional matrix of size, each row representing a magnetocardiogram signal sequence of a channel; />For a complete magnetocardiogram signal period sample, the unit is millivolt, which is 36 rows, the column number of each row is along with +.>、/>And changes in the numerical value of (c).
S102: for complete magnetocardiogram signal period samplePerforming integral R peak alignment operation to obtain an aligned cardiac magnetic signal set；
Wherein, byThe R peak of the first row of sequences is taken as a reference, and left and right translation is carried out on the rest 35 rows of sequences, so that the R peak of the rest 35 rows of sequences is aligned with the R peak of the first row of sequences; taking forward with R peak position->Column data, get back->The column data and the column data of R peak together form an aligned cardiac magnetic signal set +.>In millivolts, is oneA two-dimensional matrix of size.
S103: will align the set of magnetocardiogram signalsPerforming integral normalization operation to generate normalized magnetocardiogram signal set；
Wherein,and->Is->Maximum and minimum of>For normalizing the set of magnetocardiographic signals, is a +.>A two-dimensional matrix of size.
Further, in the step S101、/>The selection rule of (1) is->Is>Calculating data point +.1R peak to 2R peak contained in the 1 st R peak based on 3R peaks of the row sequence>Peak 2 to peak 2Data points contained by 3R peaks；/>To take the forward number of the data point of the 2 nd R peak +.>Data points of>To take the data point of the 2 nd R peak to be backward by +.>If->And->The decimal fraction is rounded up. Intercepting according to a formula (3) to obtain a complete magnetocardiogram signal period sample +.>Is>And row data. At the time of interception->Is>In line sequence, according to->、/>Is recalculated->、/>And intercepting until->The 36-line data interception of (2) is completed.
Further, in the step S102、/>The selection rule of (1) is that +.>After the R peaks of the 36-line sequence of (a) are aligned with the R peaks of the first-line sequence, respectively calculating the number of data points contained before the R peaks of each line, and selecting the least number of points as +.>The method comprises the steps of carrying out a first treatment on the surface of the The data points of each row after the R peak are calculated respectively, and the data points with the least points are selected as +.>The method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>、/>Always satisfy、/>。
Further, the step B: for normalized magnetocardiogram signal setDrawing concentric rings to construct a circle graph according to normalized magnetocardiogram signal set +.>Cutting concentric circles in the circle graph into sectors to give sectors and +.>The position relation of the data is based on +.>Performing color filling operation on the data to obtain a magnetocardiogram color space circle diagram>The method comprises the following specific steps:
s104: for the preprocessed 36-channel magnetocardiogram signals, namely normalized magnetocardiogram signal setsDrawing concentric ringsAnd construct a circle +.>；
Wherein,is>Circular ring corresponding to row->Is 1, the inner circle radius is +.>The radius of the outer ring of the circular ring is +.>The method comprises the steps of carrying out a first treatment on the surface of the 36 circles with the width of 1 and different radiuses are added together by taking the circle center as a base point to construct a circle diagram with the radius of 36>。
S105: to be used forIs>Line sequence Length will circle->Is>The rings, i.e. concentric rings +.>Cutting into sector shape to form space round diagram->And give +.>Is>The data of the row sequence corresponds to the space circle diagram +.>Is a sector of the location of (2);
wherein, concentric ringsIs cut into->Fan-shaped->Representation->Is>Line sequence>Data points, start angle of sector, +.>Representation->Is>Line sequence>Data points, the suspension angle of the sector; space round diagram->There are 36 concentric rings, each divided into +.>Sector-shaped, co-operatingA plurality of sectors; />Is>Line sequence>Data points in space circle diagram->The position of (2) is +.>Middle->Concentric rings->Is>Sector-shaped areas.
S106: will beThe matrix values map to a linear mapping to an RGB image matrix and are as +.>Data points and space circle diagrams->Position correspondence pair->Color filling operation is performed on the sector area of (2) and a magnetocardiogram color space circle diagram is drawn>；
The RGB image matrix refers to a 3-channel color image standard matrix, and the numerical value of each data point of the matrix corresponds to one color;is +.>The figure is a color circle image, namely a magnetocardiogram color space circle image, the sector position of which is equal to +.>Matrix data positions are in one-to-one correspondence, and color values thereof represent +.>The data values of the matrix.
Further, the step C: centering magnetic color space circle diagramDrawing an edge magnetocardiogram space circle graph for binarization and noise reduction processing and edge detection of the processed image>Then->Performing Hough transformation linear detection, and calculating the time of the P wave band and the QRS wave band of the magnetocardiogram signal according to the angle value of the linear, wherein the method comprises the following specific steps:
s107: centering magnetic color space circle diagramBinarizing and denoising the binarized image to obtain a binarized magnetocardiogram space circle figure +.>；
Wherein the binarization rule is formula (7), i.e. graphThe middle color value is within the interval->The color in is filled with white, figure->Is filled with black, 255 represents white, 0 represents black, interval +.>Determining a specific value by a doctor;the noise reduction treatment process of the binarized image comprises three steps of hole filling of the binarized image, corrosion of the binarized image and expansion of the binarized image, and the two steps are all image processing methods.
S108: for a pair ofPerforming edge detection to obtain edge magnetocardiogram space circle diagram +.>；
The edge detection operation is an image processing method, which can draw the boundary of the binary image, namely, the black-white boundary of the binary image will become a white line, and finally, an image edge map with black background color and white line is output.
S109: edge-to-edge magnetocardiogram space circle diagramDetecting 4 optimal straight lines by Hough transformation, and visually displaying the 4 optimal straight lines on a magnetocardiogram color space circular graph;
the Hough transform is a straight line detection method, and can output an edge magnetocardiogram space circular graph through the Hough transformThe position of the straight line on the upper surface; in the invention, only P wave band and QRS wave band time are marked, and only 4 optimal straight lines are detected; the visual display means that the 4 lines are arranged in the magnetocardiogram color space round chart according to the position data of the 4 best lines>Drawing and displaying; the included angle area of the 1 st straight line and the 2 nd straight line is the P wave band data area of the magnetocardiogram signal, and the included angle area of the 3 rd straight line and the 4 th straight line is the QRS wave band data area of the magnetocardiogram signal.
S110: according to the detected 4 optimal straight lines and magnetocardiogram color space circular diagramCalculating the time of the P wave band and the QRS wave band of the magnetocardiogram signal according to the angle value of the 0-degree line included angle;
wherein,indicating detected->Optimal straight line and magnetocardiogram color space circle diagram +.>The angle value of the 0-degree line included angle is calculated in the anticlockwise direction, and the value range is +.>；/>Indicating detected->The time value corresponding to the optimal straight line is expressed in seconds; magnetocardiogram color space circle diagram->The 0 degree line of (2) refers to the graph +.>The upper 0 degree angle corresponds to rays; p-band time is +.>QRS band time is +.>。
In yet another aspect, the invention also discloses a computer readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the steps of the method as described above.
In yet another aspect, the invention also discloses a computer device comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of the method as above.
As can be seen from the technical scheme, the magnetocardiogram color space circle map generation and wave band time identification method comprises the steps of obtaining 36-channel magnetocardiogram signal sets based on SQUID probe measurementA window function is designed to intercept to obtain a complete magnetocardiogram signal period sample, R peak alignment operation and integral normalization operation are carried out on intercepted magnetocardiogram signals, pretreatment of data is realized, and a normalized magnetocardiogram signal set is obtained; secondly, a circle diagram is constructed by drawing concentric circles of the preprocessed 36-channel magnetocardiogram signal set, the concentric circles in the circle diagram are subjected to fan-shaped cutting by calculating angle values according to the length of the magnetocardiogram signal sequence, a space circle diagram is generated, meanwhile, the position relation between the normalized magnetocardiogram signal set and the fan shape is given, color filling is carried out on the fan shape at the corresponding position according to the data values of the normalized magnetocardiogram signal set, a magnetocardiogram color space circle diagram is obtained, and conversion of the 36-channel magnetocardiogram signal to the magnetocardiogram color space circle diagram is realized; and finally, binarizing and denoising the magnetocardiogram color space circular graph, performing edge detection on the processed image, performing Hough transformation linear detection on the edge image generated by the edge detection, visually displaying the linear, calculating the P wave band and the QRS wave band time of the magnetocardiogram signal at a linear angle, and realizing the identification of the wave band time. The invention can complete the generation of the magnetocardiogram signal, the magnetocardiogram space circle graph comprises space position information and time change trend in a time period, and can automatically identify Q wave band and QRS wave band in the time change process and calculate wave band time in the magnetocardiogram space circle graph.
In general, the magnetocardiogram color space circle graph generation and wave band time identification method can realize that the space position information and the time change trend of magnetocardiogram signal data are displayed on the same image at the same time, namely, the magnetocardiogram color space circle graph is generated; by adopting the image processing method, the Q wave band and the QRS wave band can be automatically identified in the magnetic color space circular graph, and the visual annotation and the wave band time calculation can be carried out on the image, so that the method has good practical application value.
The magnetocardiogram color space circle diagram of the invention corresponds to the surface of a human body, 36-channel magnetocardiogram signals are obtained based on SQUID probe measurement, the time variation trend of single-channel magnetocardiogram signals is displayed through the sector on a single circle on the circle diagram, and the spatial position up-conversion trend of 36-channel magnetocardiogram signals is displayed through the sector between the circles on the circle diagram. A doctor can observe the space information and the time information of the 36-channel magnetocardiogram signals through the magnetocardiogram color space circular graph, and simultaneously obtain the activity condition of the heart of the patient in one period according to the calculated P wave band and QRS wave band time, so as to help the doctor to judge the health state of the heart of the patient.
Drawings
FIG. 1 is a flow chart of a magnetocardiogram color space circle map generation and band time identification method in an embodiment of the invention;
FIG. 2 is a schematic illustration of a space circle diagram in an embodiment of the invention;
FIG. 3 is a magnetocardiogram color space circle diagram in an embodiment of the present invention;
FIG. 4 is a magnetocardiogram color space circle diagram with angular and concentric circular auxiliary lines in an embodiment of the present invention;
FIG. 5 is a circular view of an edge magnetocardiogram space in an embodiment of the present invention;
fig. 6 is a band visual label diagram in an embodiment of the invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention.
As shown in fig. 1, the magnetocardiogram color space circle map generation and band time identification method according to the present embodiment specifically includes the following steps,
step A: for 36-channel magnetocardiogram signal set obtained based on SQUID probe measurementThe unit is millivolts; />The time is in seconds, a window function is designed to intercept to obtain a complete magnetocardiogram signal period sample, and R peak alignment operation and integral normalization operation are carried out on the intercepted magnetocardiogram signal to obtain a normalized magnetocardiogram signal set->Preprocessing of data is achieved;
s101: 36-channel magnetocardiogram signal set based on SQUID probe measurementExpressed in the form of a sequence, i.e. 36 channel magnetocardiogram signal set +.>Designing a rectangular window function to intercept to obtain a complete magnetocardiogram signal period sample +.>；
Wherein,for the sampling amount +.>Is 36-channel magnetocardiogram signal set->Is a sampling frequency of (a);for the window function +.>、/>Respectively a window function start point and an end point, which correspond to the pair +.>When intercepting each line of (a)>、/>Are all different in value; />Is 36-channel magnetocardiogram signal set sequence, the unit is millivolt, is 36-channel magnetocardiogram signal set +.>In the form of a sequence representation of +.>A two-dimensional matrix of size, each row representing a magnetocardiogram signal sequence of a channel; />For a complete magnetocardiogram signal period sample, the unit is millivolt, which is 36 rows, the column number of each row is along with +.>、/>And changes in the numerical value of (c).
、/>The selection rule of (1) is->Is>Calculating data point +.1R peak to 2R peak contained in the 1 st R peak based on 3R peaks of the row sequence>Data points contained in the 2 nd to 3 rd R peaks +.>；/>To take the forward number of the data point of the 2 nd R peak +.>Data points of>To take the data point of the 2 nd R peak to be backward by +.>If the data points of (1)And->The decimal fraction is rounded up. Intercepting according to a formula (3) to obtain a complete magnetocardiogram signal period sample +.>Is>And row data. At the time of interception->Is>In line sequence, according to->、/>Is recalculated->、/>And intercepting until->The 36-line data interception of (2) is completed.
S102: for complete magnetocardiogram signal period samplePerforming integral R peak alignment operation to obtain an aligned cardiac magnetic signal set；
Wherein, byThe R peak of the first row of sequences is taken as a reference, and left and right translation is carried out on the rest 35 rows of sequences, so that the R peak of the rest 35 rows of sequences is aligned with the R peak of the first row of sequences; taking forward with R peak position->Column data, get back->The column data and the column data of R peak together form an aligned cardiac magnetic signal set +.>In millivolts, is oneA two-dimensional matrix of size. />、/>The selection rule of (1) is that +.>After the R peaks of the 36-line sequence of (a) are aligned with the R peaks of the first-line sequence, respectively calculating the number of data points contained before the R peaks of each line, and selecting the least number of points as +.>The method comprises the steps of carrying out a first treatment on the surface of the The data points of each row after the R peak are calculated respectively, and the data points with the least points are selected as +.>The method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)>、/>Always satisfy、/>。
S103: will align the set of magnetocardiogram signalsPerforming integral normalization operation to generate normalized magnetocardiogram signal set；
Wherein,and->Is->Maximum and minimum of>For normalizing the set of magnetocardiographic signals, is a +.>A two-dimensional matrix of size.
And (B) step (B): for normalized magnetocardiogram signal setDrawing concentric rings to construct a circle graph according to normalized magnetocardiogram signal set +.>Cutting concentric circles in the circle graph into sectors to give sectors and +.>The position relation of the data is based on +.>Performing color filling operation on the data to obtain a magnetocardiogram color space circle diagram>；
S104: for the preprocessed 36-channel magnetocardiogram signals, namely normalized magnetocardiogram signal setsDrawing concentric ringsAnd construct a circle +.>；
Wherein,is>Circular ring corresponding to row->Is 1, the inner circle radius is +.>The radius of the outer ring of the circular ring is +.>The method comprises the steps of carrying out a first treatment on the surface of the 36 circles with the width of 1 and different radiuses are added together by taking the circle center as a base point to construct a circle diagram with the radius of 36>。
S105: to be used forIs>Line sequence Length will circle->Is>The rings, i.e. concentric rings +.>Cutting into sector shape to form space round diagram->And give +.>Is>The data of the row sequence corresponds to the space circle diagram +.>Is a sector of the location of (2);
wherein, concentric ringsIs cut into->Fan-shaped->Representation->Is>Line sequence>Data points, start angle of sector, +.>Representation->Is>Line sequence>Data points, the suspension angle of the sector; space round diagram->There are 36 concentric rings, each divided into +.>Sector-shaped, co-operatingA plurality of sectors; />Is>Line sequence>Data points in space circle diagram->The position of (2) is +.>Middle->Concentric rings->Is>Sector-shaped areas.
As shown in FIG. 2, a space circle diagram is schematically shown, which is formed by overlapping 36 concentric circles of width 1 as shown in FIG. 2, and each concentric circle is divided into sections as shown in FIG. 2Sector of->Is shown, and a line of 0 degrees in the positive X-axis direction is shown.
S106: will beThe matrix values map to a linear mapping to an RGB image matrix and are as +.>Data points and space circle diagrams->Position correspondence pair->Color filling operation is performed on the sector area of (2) and a magnetocardiogram color space circle diagram is drawn>；
The RGB image matrix refers to a 3-channel color image standard matrix, and the numerical value of each data point of the matrix corresponds to one color;is +.>The figure is a color circle image, namely a magnetocardiogram color space circle image, the sector position of which is equal to +.>Matrix data positions are in one-to-one correspondence, and color values thereof represent +.>The data values of the matrix.
As shown in fig. 3, the magnetocardiogram color space circle diagram is a space circle diagram after coloring. When the circle center is observed outwards in a ray form along a certain angle in the figure, the change trend of the 36-channel magnetocardiogram signal at a certain moment, namely the spatial relationship is represented. The time variation trend, namely the time relation of the 36-channel magnetocardiogram signals is represented by the 0-degree line which is observed in the anticlockwise direction in the figure. In the figure, the intensity of different color areas represents the intensity of the heart activity of the human body, the lighter the color is, the gentle the heart activity is, and the darker the color is, the strong the heart activity is.
As shown in FIG. 4, the magnetocardiogram color space circle graph with angle and concentric ring auxiliary lines is marked with angle lines and 36 concentric ring dividing lines, which can effectively and rapidly help doctors to locate magnetocardiogram signal positions of each channel and initially determine angle intervals of each wave band.
Step C: centering magnetic color space circle diagramDrawing an edge magnetocardiogram space circle graph for binarization and noise reduction processing and edge detection of the processed image>Then->Performing Hough transformation linear detection, and calculating the time of the P wave band and the QRS wave band of the magnetocardiogram signal according to the linear angle value;
s107: centering magnetic color space circle diagramBinarizing and denoising the binarized image to obtain a binarized magnetocardiogram space circle figure +.>；
Wherein the binarization rule is formula (7), i.e. graphThe middle color value is within the interval->The color in is filled with white, figure->Is filled with black, 255 represents white, 0 represents black, interval +.>Determining a specific value by a doctor; the noise reduction treatment process of the binarized image comprises three steps of hole filling of the binarized image, corrosion of the binarized image and expansion of the binarized image, and the two steps are all image processing methods.
S108: for a pair ofPerforming edge detection to obtain edge magnetocardiogram space circle diagram +.>；
The edge detection operation is an image processing method, which can draw the boundary of the binary image, namely, the black-white boundary of the binary image will become a white line, and finally, an image edge map with black background color and white line is output.
As shown in fig. 5, the edge magnetocardiogram space circle diagram is shown, the circle diagram boundary and each band boundary are effectively displayed in the diagram, hough transformation detection can be performed according to the diagram, and each band boundary is marked.
S109: edge-to-edge magnetocardiogram space circle diagramDetecting 4 optimal straight lines by Hough transformation, and visually displaying the 4 optimal straight lines on a magnetocardiogram color space circular graph;
the Hough transform is a straight line detection method, and can output an edge magnetocardiogram space circular graph through the Hough transformThe position of the straight line on the upper surface; in the invention, only P wave band and QRS wave band time are marked, and only 4 optimal straight lines are detected; the visual display means that the 4 lines are arranged in the magnetocardiogram color space round chart according to the position data of the 4 best lines>Drawing and displaying; the included angle area of the 1 st straight line and the 2 nd straight line is the P wave band data area of the magnetocardiogram signal, and the included angle area of the 3 rd straight line and the 4 th straight line is the QRS wave band data area of the magnetocardiogram signal.
S110: according to the detected 4 optimal straight lines and magnetocardiogram color space circular diagramCalculating the time of the P wave band and the QRS wave band of the magnetocardiogram signal according to the angle value of the 0-degree line included angle;
wherein,indicating detected->Optimal straight line and magnetocardiogram color space circle diagram +.>The angle value of the 0-degree line included angle is calculated in the anticlockwise direction, and the value range is +.>；/>Indicating detected->The time value corresponding to the optimal straight line is expressed in seconds; magnetocardiogram color space circle diagram->The 0 degree line of (2) refers to the graph +.>The upper 0 degree angle corresponds to rays; p-band time is +.>QRS band time is +.>。
As shown in fig. 6, a band visual label is shown, in which a boundary line, a start angle and a stop angle of a P band are marked, the start angle of the P band is 118.8608 degrees, the stop angle is 137.1949 degrees, and the position of the P band in the magnetocardiogram color space circle diagram is a sector area divided by corresponding straight lines of the two angles. The boundary line, the start angle and the stop angle of the QRS wave band are marked in the figure, the start angle of the QRS wave band is 207.3324 degrees, the stop angle is 265.1894 degrees, and the position of the QRS wave band in the magnetocardiogram color space circular diagram is a divided sector area of the straight line corresponding to the two angles. The P-band time can be calculated to be 0.0326 seconds and the QRS time 0.1029 seconds.
In summary, the magnetocardiogram color space circle graph generation and band time identification method of the embodiment of the invention comprises the following steps: (1) For 36-channel magnetocardiogram signal set obtained based on SQUID probe measurementThe unit is millivolts; />The time is in seconds, a window function is designed to intercept to obtain a complete magnetocardiogram signal period sample, and R peak alignment operation and integral normalization operation are carried out on the intercepted magnetocardiogram signal to obtain a normalized magnetocardiogram signal set->Preprocessing of data is achieved; (2) For normalized cardiac magnetic signal set->Drawing concentric rings to construct a circle graph according to normalized magnetocardiogram signal set +.>Cutting concentric circles in the circle graph into sectors to give sectors and +.>The position relation of the data is based on +.>Performing color filling operation on the data to obtain a magnetocardiogram color space circle diagram>The method comprises the steps of carrying out a first treatment on the surface of the (3) Centered magnetic color space circle diagram->Drawing an edge magnetocardiogram space circle graph for binarization and noise reduction processing and edge detection of the processed image>Then->And carrying out Hough transformation linear detection, and calculating the time of the P wave band and the QRS wave band of the magnetocardiogram signal according to the linear angle value. The invention can generate a magnetocardiogram color space circle graph, and effectively displays the time change trend and the space position up-conversion trend of 36-channel magnetocardiogram signals on the same image. In addition, the P wave band and the QRS wave band are automatically marked by the methods of image processing and straight line detection, and the time is calculated. The magnetocardiogram color space circle map generation and wave band time identification method can help doctors to judge the health state of the heart of the patient, and has good practical application value.
In yet another aspect, the invention also discloses a computer readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the steps of the method as described above.
In yet another aspect, the invention also discloses a computer device comprising a memory and a processor, the memory storing a computer program which, when executed by the processor, causes the processor to perform the steps of the method as above.
In yet another embodiment provided herein, there is also provided a computer program product containing instructions that, when run on a computer, cause the computer to perform the magnetocardiogram color space circular map generation and band time identification method of any of the above embodiments.
It may be understood that the system provided by the embodiment of the present invention corresponds to the method provided by the embodiment of the present invention, and explanation, examples and beneficial effects of the related content may refer to corresponding parts in the above method.
The embodiment of the application also provides an electronic device, which comprises a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface and the memory are communicated with each other through the communication bus,
a memory for storing a computer program;
and the processor is used for realizing the magnetocardiogram color space circle graph generation and band time identification method when executing the program stored in the memory.
The communication bus mentioned by the above electronic device may be a peripheral component interconnect standard (english: peripheral Component Interconnect, abbreviated: PCI) bus or an extended industry standard architecture (english: extended Industry Standard Architecture, abbreviated: EISA) bus, or the like. The communication bus may be classified as an address bus, a data bus, a control bus, or the like.
The communication interface is used for communication between the electronic device and other devices.
The Memory may include random access Memory (Random Access Memory, abbreviated as RAM) or nonvolatile Memory (NVM), such as at least one disk Memory. Optionally, the memory may also be at least one memory device located remotely from the aforementioned processor.
The processor may be a general-purpose processor, including a central processing unit (Central Processing Unit, CPU for short), a network processor (Network Processor, NP for short), etc.; it may also be a digital signal processor (English: digital Signal Processing; DSP; for short), an application specific integrated circuit (English: application Specific Integrated Circuit; ASIC; for short), a Field programmable gate array (English: field-Programmable Gate Array; FPGA; for short), or other programmable logic device, discrete gate or transistor logic device, discrete hardware components.
In the above embodiments, it may be implemented in whole or in part by software, hardware, firmware, or any combination thereof. When implemented in software, may be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When loaded and executed on a computer, produces a flow or function in accordance with embodiments of the present application, in whole or in part. The computer may be a general purpose computer, a special purpose computer, a computer network, or other programmable apparatus. The computer instructions may be stored in or transmitted from one computer-readable storage medium to another, for example, by wired (e.g., coaxial cable, optical fiber, digital Subscriber Line (DSL)), or wireless (e.g., infrared, wireless, microwave, etc.). The computer readable storage medium may be any available medium that can be accessed by a computer or a data storage device such as a server, data center, etc. that contains an integration of one or more available media. The usable medium may be a magnetic medium (e.g., floppy Disk, hard Disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., solid State Disk (SSD)), etc.
It is noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article, or apparatus that comprises the element.
In this specification, each embodiment is described in a related manner, and identical and similar parts of each embodiment are all referred to each other, and each embodiment mainly describes differences from other embodiments. In particular, for system embodiments, since they are substantially similar to method embodiments, the description is relatively simple, as relevant to see a section of the description of method embodiments.
The above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (8)
1. A magnetocardiogram color space circle diagram generation and wave band time identification method is characterized by comprising the following steps,
step A: for 36-channel magnetocardiogram signal set obtained based on SQUID probe measurementThe unit is millivolts; />The time is in seconds, a window function is designed to intercept to obtain a complete magnetocardiogram signal period sample, and R peak alignment operation and integral normalization operation are carried out on the intercepted magnetocardiogram signal to obtain a normalized magnetocardiogram signal set->Preprocessing of data is achieved;
and (B) step (B): for normalized magnetocardiogram signal setDrawing concentric rings to construct a circle graph and normalizing the magnetocardiogram signal setCutting concentric circles in the circle graph into sectors to give sectors and +.>The position relation of the data is based on +.>Performing color filling operation on the data to the corresponding position sector to obtain a magnetocardiogram color space circle map；
Step C: centering magnetic color space circle diagramDrawing an edge magnetocardiogram space circle graph for binarization and noise reduction processing and edge detection of the processed image>Then->Performing Hough transformation linear detection, and calculating the time of the P wave band and the QRS wave band of the magnetocardiogram signal according to the linear angle value;
the step C comprises the following steps of, in particular,
s107: centering magnetic color space circle diagramBinarizing and denoising the binarized image to obtain a binarized magnetocardiogram space circle figure +.>；
Wherein the binarization rule is formula (7), i.e. graphThe middle color value is within the interval->The color in is filled with white, figure->Is filled with black, 255 represents white, 0 represents black, interval +.>Determining a specific value by a doctor;
s108: for a pair ofPerforming edge detection to obtain edge magnetocardiogram space circle diagram +.>；
S109: edge-to-edge magnetocardiogram space circle diagramDetecting 4 optimal straight lines by Hough transformation, and visually displaying the 4 optimal straight lines on a magnetocardiogram color space circular graph;
s110: according to the detected 4 optimal straight lines and magnetocardiogram color space circular diagramCalculating the time of the P wave band and the QRS wave band of the magnetocardiogram signal according to the angle value of the 0-degree line included angle;
wherein,indicating detected->Optimal straight line and magnetocardiogram color space circle diagram +.>The angle value of the 0-degree line included angle is calculated in the anticlockwise direction, and the value range is +.>；/>Indicating detected->The time value corresponding to the optimal straight line is expressed in seconds; magnetocardiogram color space circle diagram->The 0 degree line of (2) refers to the graph +.>The upper 0 degree angle corresponds to rays; p-band time is +.>QRS band time is +.>；
The step A specifically comprises the following steps of,
s101: 36-channel magnetocardiogram signal set based on SQUID probe measurementExpressed in the form of a sequence, i.e. 36 channel magnetocardiogram signal set +.>Designing a rectangular window function to intercept to obtain a complete magnetocardiogram signal period sample +.>；
Wherein,for the sampling amount +.>Is 36-channel magnetocardiogram signal set->Is a sampling frequency of (a); />For the window function +.> 、/>The starting point and the end point of the window function are respectively; />Is 36-channel magnetocardiogram signal set sequence, the unit is millivolt, is 36-channel magnetocardiogram signal set +.>In the form of a sequence representation of +.>A two-dimensional matrix of size, each row representing a magnetocardiogram signal sequence of a channel; />The unit is millivolt for a complete magnetocardiogram signal period sample;
s102: for complete magnetocardiogram signal period samplePerforming whole R peak alignment operation to obtain alignment cardiac magnetic signal set +.>；
Wherein, byThe R peak of the first row of sequences is taken as a reference, and left and right translation is carried out on the rest 35 rows of sequences, so that the R peak of the rest 35 rows of sequences is aligned with the R peak of the first row of sequences; taking forward with R peak position->Column data, get back->The column data and the column data of R peak together form an aligned cardiac magnetic signal set +.>In millivolts, is oneA two-dimensional matrix of size;
s103: will align the set of magnetocardiogram signalsCarrying out integral normalization operation to generate normalized magnetocardiogram signal set +.>；
Wherein,and->Is->Maximum and minimum of>For normalizing the heart magnetic signal set, is oneA two-dimensional matrix of size.
2. The magnetocardiogram color space circle map generation and band time identification method according to claim 1, wherein: in the step S101 、/>The selection rule of (1) is->Is>Calculating data point +.1R peak to 2R peak contained in the 1 st R peak based on 3R peaks of the row sequence>Data points contained in the 2 nd to 3 rd R peaks +.>；/>To take the forward number of the data point of the 2 nd R peak +.>Data points of>To take the data point of the 2 nd R peak to be backward by +.>If->And->Rounding up for decimal places; intercepting according to a formula (3) to obtain a complete magnetocardiogram signal period sample +.>Is>Line data; at the time of interception->Is>In line sequence, according to-> 、/>Is recalculated-> 、/>And intercepting until->The 36-line data interception of (2) is completed.
3. The magnetocardiogram color space circle map generation and band time identification method according to claim 2, characterized in that: in the step S102 、/>The selection rule of (1) is that +.>After the R peaks of the 36-line sequence of (a) are aligned with the R peaks of the first-line sequence, respectively calculating the number of data points contained before the R peaks of each line, and selecting the least number of points as +.>The method comprises the steps of carrying out a first treatment on the surface of the The data points of each row after the R peak are calculated respectively, and the data points with the least points are selected as +.>The method comprises the steps of carrying out a first treatment on the surface of the Wherein (1)> 、/>Always satisfy->、/>。
4. The magnetocardiogram color space circle map generation and band time identification method according to claim 1, wherein: the step B comprises the following steps of, in particular,
s104: for the preprocessed 36-channel magnetocardiogram signals, namely normalized magnetocardiogram signal setsDrawing concentric circle +.>And construct a circle +.>；
Wherein,is>Circular ring corresponding to row->Is 1, the inner circle radius is +.>The radius of the outer ring of the circular ring is +.>The method comprises the steps of carrying out a first treatment on the surface of the 36 circles with the width of 1 and different radiuses are added together by taking the circle center as a base point to construct a circle diagram with the radius of 36>；
S105: to be used forIs>Line sequence Length will circle->Is>The rings, i.e. concentric rings +.>Cutting into sector shape to form space round diagram->And give +.>Is>The data of the row sequence corresponds to the space circle diagram +.>Is a sector of the location of (2);
wherein, concentric ringsIs cut into->Fan-shaped->Representation->Is>Line sequence>Data points, start angle of sector, +.>Representation->Is>Line sequence>Data points, the suspension angle of the sector; space round diagram->There are 36 concentric rings, each divided into +.>Sector-shaped, co-operatingA plurality of sectors; />Is>Line sequence>Data points in space circle diagram->The position of (2) is +.>Middle->Concentric rings->Is>A plurality of sector-shaped regions;
s106: will beThe matrix values map to a linear mapping to an RGB image matrix and are as +.>Data points and space circle diagramsPosition correspondence pair->Color filling operation is performed on the sector area of (2) and a magnetocardiogram color space circle diagram is drawn>；
The RGB image matrix refers to a 3-channel color image standard matrix, and the numerical value of each data point of the matrix corresponds to one color;is +.>The figure is a color circle image, namely a magnetocardiogram color space circle image, the sector position of which is equal to +.>Matrix data positions are in one-to-one correspondence, and color values thereof represent +.>The data values of the matrix.
5. The magnetocardiogram color space circle map generation and band time identification method according to claim 1, wherein: in step S107, the noise reduction processing of the binarized image includes three steps, i.e., hole filling of the binarized image, erosion of the binarized image, and expansion of the binarized image, which are all image processing methods.
6. The magnetocardiogram color space circle map generation and band time identification method according to claim 1, wherein:
in step S108, the edge detection operation draws the boundary of the binary image, i.e. the black-white boundary of the binary image will become a white line, and finally outputs an image edge map with the black line as the background color.
7. The magnetocardiogram color space circle map generation and band time identification method according to claim 1, wherein:
in step S109, the edge magnetocardiogram space circular map is output through Hough transformThe position of the straight line on the upper surface; only the P wave band and the QRS wave band time are marked, and only 4 optimal straight lines are detected; the visual display means that the 4 lines are arranged in the magnetocardiogram color space round chart according to the position data of the 4 best lines>Drawing and displaying; the included angle area of the 1 st straight line and the 2 nd straight line is the P wave band data area of the magnetocardiogram signal, and the included angle area of the 3 rd straight line and the 4 th straight line is the QRS wave band data area of the magnetocardiogram signal.
8. A computer readable storage medium storing a computer program which, when executed by a processor, causes the processor to perform the steps of the method of any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311234687.0A CN116965822B (en) | 2023-09-25 | 2023-09-25 | Magnetocardiogram color space circle diagram generation and wave band time identification method and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311234687.0A CN116965822B (en) | 2023-09-25 | 2023-09-25 | Magnetocardiogram color space circle diagram generation and wave band time identification method and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116965822A CN116965822A (en) | 2023-10-31 |
| CN116965822B true CN116965822B (en) | 2023-12-29 |
Family
ID=88485336
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311234687.0A Active CN116965822B (en) | 2023-09-25 | 2023-09-25 | Magnetocardiogram color space circle diagram generation and wave band time identification method and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116965822B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117274117B (en) * | 2023-11-23 | 2024-02-02 | 合肥工业大学 | Frequency domain pseudo-color enhanced magnetocardiogram signal characteristic image generation method and storage medium |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6187032B1 (en) * | 1997-10-30 | 2001-02-13 | Kabushiki Kaisha Toshiba | Measurement of intracardiac electrophysiological phenomena |
| JP2005342081A (en) * | 2004-06-01 | 2005-12-15 | Japan Science & Technology Agency | Cardiac magnetic field diagnosis apparatus and method of operating the same |
| CN101862183A (en) * | 2010-06-13 | 2010-10-20 | 北京大学 | Method for quickly determining magnetocardiogram measurement range based on magnetocardiogram characteristics |
| WO2012173583A1 (en) * | 2011-06-16 | 2012-12-20 | Chaykovskyy Lllya Anatoliiovych | Method and device for evaluation of myocardial ischemia based on current density maps |
| KR101525485B1 (en) * | 2014-03-21 | 2015-06-03 | 한국표준과학연구원 | Three-dimensional Reconstruction of a Cardiac Outline by Magnetocardiography |
| CN107374610A (en) * | 2017-07-12 | 2017-11-24 | 中国科学院上海微系统与信息技术研究所 | Magnetocardiogram generation method and generation system |
| EP3308703A1 (en) * | 2016-10-11 | 2018-04-18 | Biomagnetik Park GmbH | Magnetocardiographic method and magnetocardiographic system |
| CN108802540A (en) * | 2018-07-19 | 2018-11-13 | 山东大学 | Method for visualizing suitable for power distribution network synchronized phasor real-time state monitoring |
| CN113779841A (en) * | 2021-09-14 | 2021-12-10 | 哈尔滨理工大学 | Stator transposition winding transient circulating current calculation method based on three-dimensional electromagnetic coupling network |
| CN115337016A (en) * | 2022-08-01 | 2022-11-15 | 苏州卡迪默克医疗器械有限公司 | Method and device for evaluating myocardial ischemia based on current density vector diagram |
| CN116712056A (en) * | 2023-08-07 | 2023-09-08 | 合肥工业大学 | Characteristic image generation and identification method, equipment and storage medium for electrocardiogram data |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7805179B2 (en) * | 2007-06-08 | 2010-09-28 | Herng-Er Horng | Method of examining dynamic cardiac electromagnetic activity and detection of cardiac functions using results thereof |
-
2023
- 2023-09-25 CN CN202311234687.0A patent/CN116965822B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6187032B1 (en) * | 1997-10-30 | 2001-02-13 | Kabushiki Kaisha Toshiba | Measurement of intracardiac electrophysiological phenomena |
| JP2005342081A (en) * | 2004-06-01 | 2005-12-15 | Japan Science & Technology Agency | Cardiac magnetic field diagnosis apparatus and method of operating the same |
| CN101862183A (en) * | 2010-06-13 | 2010-10-20 | 北京大学 | Method for quickly determining magnetocardiogram measurement range based on magnetocardiogram characteristics |
| WO2012173583A1 (en) * | 2011-06-16 | 2012-12-20 | Chaykovskyy Lllya Anatoliiovych | Method and device for evaluation of myocardial ischemia based on current density maps |
| KR101525485B1 (en) * | 2014-03-21 | 2015-06-03 | 한국표준과학연구원 | Three-dimensional Reconstruction of a Cardiac Outline by Magnetocardiography |
| EP3308703A1 (en) * | 2016-10-11 | 2018-04-18 | Biomagnetik Park GmbH | Magnetocardiographic method and magnetocardiographic system |
| CN107374610A (en) * | 2017-07-12 | 2017-11-24 | 中国科学院上海微系统与信息技术研究所 | Magnetocardiogram generation method and generation system |
| CN108802540A (en) * | 2018-07-19 | 2018-11-13 | 山东大学 | Method for visualizing suitable for power distribution network synchronized phasor real-time state monitoring |
| CN113779841A (en) * | 2021-09-14 | 2021-12-10 | 哈尔滨理工大学 | Stator transposition winding transient circulating current calculation method based on three-dimensional electromagnetic coupling network |
| CN115337016A (en) * | 2022-08-01 | 2022-11-15 | 苏州卡迪默克医疗器械有限公司 | Method and device for evaluating myocardial ischemia based on current density vector diagram |
| CN116712056A (en) * | 2023-08-07 | 2023-09-08 | 合肥工业大学 | Characteristic image generation and identification method, equipment and storage medium for electrocardiogram data |
Non-Patent Citations (1)
| Title |
|---|
| Multi-channel high-Tc SQUID system for bio-applications;Y. Hatsukade;Solid State Phenomena;424-427 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116965822A (en) | 2023-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116965822B (en) | Magnetocardiogram color space circle diagram generation and wave band time identification method and storage medium | |
| US20090099468A1 (en) | Automated Processing of Electrophysiological Data | |
| CA2912476C (en) | Electrocardiogram analysis | |
| CN113066061B (en) | MRA-based aneurysm detection method, system, terminal and medium | |
| EP4042445A1 (en) | Systems and methods for reduced lead electrocardiogram diagnosis using deep neural networks and rule-based systems | |
| CN109009084B (en) | QRS wave group calibration method, device, equipment and medium for multi-lead electrocardiosignal | |
| CN108577825B (en) | Myocardial ischemia lesion position positioning method and system, storage medium and terminal | |
| WO2013111031A1 (en) | Method and system for cardiac ischemia detection | |
| WO2021071646A1 (en) | Systems and methods for electrocardiogram diagnosis using deep neural networks and rule-based systems | |
| Arnold et al. | Discriminating electrocardiographic responses to His-bundle pacing using machine learning | |
| Fortune et al. | Digitizing ECG image: a new method and open-source software code | |
| US20230222668A1 (en) | Image processing apparatus, image processing method, and recording medium | |
| Kim et al. | Identification of atrial fibrillation with single-lead mobile ECG during normal sinus rhythm using deep learning | |
| CN115177234B (en) | High-density sensor, high-density detection device, data processing method and imaging method | |
| CN109431488B (en) | Method for identifying typical data in electrocardiogram dynamic real-time analysis data | |
| CN114403812B (en) | Auxiliary analysis method, device and system for brain injury condition and storage medium | |
| CN110916649B (en) | Processing device, processing method and detection device for long-range electrocardiogram scatter diagram | |
| CN110464339B (en) | Electrocardiogram analysis method and device | |
| Agrawal et al. | A machine learning approach for classifying healthy and infarcted patients using heart rate variabilities derived vector magnitude | |
| US20200170624A1 (en) | Diagnostic apparatus and diagnostic method | |
| Saadi et al. | Low-cost ecg monitoring system with classification using deep learning | |
| CN112386264B (en) | Method and device for carrying out graphical display on ST-section data in multi-lead data | |
| CN116342497B (en) | Three-dimensional mapping method and system for inner wall of human body cavity | |
| Azman et al. | Variability in the atrial flutter vectorcardiographic loops and non-invasive localization of circuits | |
| CN117274413B (en) | EIT-based conductivity image reconstruction method, system and equipment |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |