Magnetocardiogram multipole track analysis method, device and equipment and readable storage medium

Abstract

The invention provides a magnetocardiogram multipole track analysis method, device and equipment and a readable storage medium. The method comprises the following steps: acquiring a multi-pole track diagram based on an original magnetocardiogram set; counting a plurality of track branches in the multi-pole track diagram, and dividing the plurality of track branches into a positive main branch, a negative main branch and a secondary branch; based on the positive main branch, the negative main branch and the secondary branch, determining a multi-pole sub-track score corresponding to the multi-pole sub-track graph; the trajectory deviation degree of the original magnetocardiogram set and a preset magnetocardiogram set is determined on the basis of the multi-pole trajectory score, a trajectory analysis result of the original magnetocardiogram set is obtained on the basis of the trajectory deviation degree, and by extracting trajectory parameters of different branches, cardiac electrophysiological change details behind different wavebands can be effectively described; according to the method, accurate information extraction of time-space continuity, mutability and energy distribution change of the multipole track is realized, information in the magnetocardiogram can be comprehensively described, and the accuracy of a magnetocardiogram analysis result is improved.

Description

Technical Field

[0001] This application belongs to the field of data technology, and relates to a method for analyzing the trajectories of multipoles in magnetocardiogram, in particular to a method, device, equipment and readable storage medium for analyzing the trajectories of multipoles in magnetocardiogram. Background Art

[0002] Magnetocardiogram is a method for detecting the magnetic field of the heart using superconducting quantum interference devices or other advanced magnetic sensor technologies. It has significant advantages such as non-contact, high spatio-temporal resolution, and high sensitivity to deep myocardial electrical activities. Existing magnetocardiogram analysis methods mainly rely on the single dipole model or only generally analyze the overall multipole degree. However, due to the weakness of the magnetocardiogram signal itself and the complexity of the heart electromagnetic field distribution, the above methods have limitations in decoupling the multipole characteristics of the magnetic field distribution, seriously affecting the accuracy of the analysis results, especially in capturing the continuity, direction mutation, and energy distribution change of the multipole trajectories. These multipole characteristics may contain important cardiac electrophysiological information. Therefore, how to improve the accuracy of the analysis results has become a technical problem that needs to be solved urgently by relevant technical personnel. Summary of the Invention

[0003] This application provides a method, device, equipment and readable storage medium for analyzing the trajectories of multipoles in magnetocardiogram, which is used to solve the technical problem of lacking accurate analysis results in the prior art when analyzing magnetocardiogram.

[0004] In a first aspect, an embodiment of this application provides a method for analyzing the trajectories of multipoles in magnetocardiogram. The method includes: obtaining a multipole trajectory map based on an original magnetocardiogram set; determining multiple trajectory branches in the multipole trajectory map, and dividing the multiple trajectory branches into a positive main branch, a negative main branch, and secondary branches; determining a multipole trajectory score corresponding to the multipole trajectory map based on the positive main branch, the negative main branch, and the secondary branches; determining a trajectory deviation degree between the original magnetocardiogram set and a preset magnetocardiogram set based on the multipole trajectory score, and obtaining a trajectory analysis result of the original magnetocardiogram set based on the trajectory deviation degree.

[0005] In an implementation manner of the first aspect, determining the multipole trajectory score corresponding to the multipole trajectory diagram based on the positive main branch, the negative main branch, and the secondary branch includes: determining a first trajectory parameter of the positive main branch based on the main branch information corresponding to the positive main branch, where the main branch information corresponding to the positive main branch includes the magnetic field value of the positive pole on the positive main branch and the position information of the positive pole; determining a second trajectory parameter of the negative main branch based on the main branch information corresponding to the negative main branch; matching the positive poles on the positive main branch and the negative poles on the negative main branch within the same frame to obtain main branch dipoles, and determining a third trajectory parameter of the main branch dipoles; determining a fourth trajectory parameter of the secondary branch based on the secondary branch information corresponding to the secondary branch; and determining the multipole trajectory score based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter.

[0006] In an implementation manner of the first aspect, the first trajectory parameter includes an overall trajectory parameter and a segmented trajectory parameter. Determining the first trajectory parameter of the positive main branch based on the main branch information corresponding to the positive main branch includes: obtaining the magnetic field value of the positive poles on the positive main branch that are greater than a preset pole threshold and the position information of the positive poles; based on the magnetic field value and the position information, determining the overall trajectory parameter of the positive main branch, where the overall trajectory parameter includes an effective length, a main pole length, and a non-boundary rate; successively determining the planar distance between adjacent positive poles on the positive main branch, and determining a main segment of the positive main branch where the planar distance is less than a preset offset distance; determining a head segment of the positive main branch located in front of the main segment; and determining a tail segment of the positive main branch located behind the main segment.

[0007] Respectively determining the segmented trajectory parameters of the main segment, the head, and the tail, where the segmented trajectory parameters include a trajectory path, a trajectory span, and a trajectory elastic curvature.

[0008] In an implementation manner of the first aspect, determining the overall trajectory parameter of the positive main branch based on the magnetic field value and the position information includes: determining the effective length corresponding to the positive main branch based on the magnetic field values of all the positive poles greater than a preset pole threshold; obtaining the main pole length corresponding to the positive main branch based on the ratio of the magnetic field value equal to the maximum magnetic field value; and obtaining the non-boundary rate corresponding to the positive main branch based on the proportion of the planar coordinates of the positive poles on the positive main branch that are not within the preset range of the preset matrix boundary.

[0009] In one implementation of the first aspect, determining the third trajectory parameter of the main branch dipole includes: determining the main branch of the main branch dipole as the intersection of the positive poles on the positive main branch and the negative poles on the negative main branch; sequentially determining the planar distances between adjacent frames of the positive poles on the positive main branch, and determining a section of the positive main branch with a planar distance less than a preset offset distance as the main section of the positive main branch; sequentially determining the planar distances between adjacent frames of the negative poles on the negative main branch, and determining a section of the negative main branch with a planar distance less than a preset offset distance as the main section of the negative main branch; determining the intersection of the main section of the positive main branch and the main section of the negative main branch as the main section of the main branch dipole; taking the difference between the main branch of the main branch dipole and the main section of the main branch dipole to obtain the secondary section of the main branch dipole; based on the main branch information, sequentially obtaining the third trajectory parameters of the main branch, the main section, and the secondary section of the main branch dipole.

[0010] In one implementation of the first aspect, the fourth trajectory parameter includes the number of effective branches, the branch size, the total trajectory path, the total trajectory span, the total trajectory elastic curvature, and the total trajectory elastic strength. Determining the fourth trajectory parameter of the secondary branch based on the secondary branch information corresponding to the secondary branch includes:

[0011] Determining the non-boundary rate of each secondary branch, removing the secondary branches with a non-boundary rate less than or equal to a preset boundary threshold to obtain the boundary factor of the secondary branch; obtaining the sum of the magnetic field values between the start time and the end time of all secondary branches; dividing the magnetic field value of each secondary branch by the sum of the magnetic field values to obtain the length factor of each secondary branch; summing all the boundary factors to obtain the number of effective branches; sequentially multiplying each boundary factor by the length factor of the secondary branch corresponding to the boundary factor and then summing to obtain the branch size of the secondary branch; determining the trajectory path, the trajectory span, and the trajectory elastic curvature of each secondary branch, and obtaining the trajectory elastic strength based on the trajectory elastic curvature; respectively summing the trajectory paths, the trajectory spans, the trajectory elastic curvatures, and the trajectory elastic strengths of all secondary branches to obtain the total trajectory path, the total trajectory span, the total trajectory elastic curvature, and the total trajectory elastic strength.

[0012] In one implementation of the first aspect, determining the multi-pole trajectory score based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter includes: respectively determining the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the secondary branch trajectory score based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter; obtaining the multi-pole trajectory score based on the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the secondary branch trajectory score.

[0013] An embodiment of the present application provides a method for analyzing the multi - dipole trajectory of a magnetocardiogram. In this method, the first trajectory parameter of the positive main branch, the second trajectory parameter of the negative main branch, the third trajectory parameter of the main - branch dipole, and the fourth trajectory parameter of the secondary branch are respectively obtained, and a multi - dipole trajectory score is determined based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter; based on the multi - dipole trajectory score, the trajectory deviation degree between the original magnetocardiogram set and a preset magnetocardiogram set is determined, and a trajectory analysis result of the original magnetocardiogram set is obtained based on the trajectory deviation degree. By extracting the trajectory parameters corresponding to different branches, the detailed electrophysiological changes of the heart behind different frequency bands can be effectively characterized, which is the information decoupling and detailed characterization of the magnetocardiogram signal in spatial distribution, realizing the accurate information extraction of the temporal - spatial continuity, abruptness, and energy - distribution changes of the multi - dipole trajectory, being able to more comprehensively characterize the information contained in the magnetocardiogram, improving the ability of magnetocardiogram analysis, and greatly improving the accuracy of the magnetocardiogram analysis result.

[0014] In a second aspect, an embodiment of the present application provides a magnetocardiogram multi - dipole trajectory analysis device. The magnetocardiogram multi - dipole trajectory analysis device includes: a multi - dipole trajectory map acquisition module, configured to acquire a multi - dipole trajectory map based on an original magnetocardiogram set; a trajectory branch determination module, configured to determine multiple trajectory branches in the multi - dipole trajectory map and divide the multiple trajectory branches into a positive main branch, a negative main branch, and a secondary branch; a multi - dipole trajectory score determination module, configured to determine a multi - dipole trajectory score corresponding to the multi - dipole trajectory map based on the positive main branch, the negative main branch, and the secondary branch; a trajectory analysis result determination module, configured to determine the trajectory deviation degree between the original magnetocardiogram set and a preset magnetocardiogram set based on the multi - dipole trajectory score, and obtain a trajectory analysis result of the original magnetocardiogram set based on the trajectory deviation degree.

[0015] In a third aspect, an embodiment of the present application provides a computer - readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the method for analyzing the multi - dipole trajectory of a magnetocardiogram according to any one of the first aspects of the embodiments of the present application is implemented.

[0016] In a fourth aspect, an embodiment of the present application provides an electronic device. The electronic device includes a memory and a processor; the memory is used to store a computer program; the processor is used to execute the computer program and implement the method for analyzing the multi - dipole trajectory of a magnetocardiogram according to any one of the first aspects of the embodiments of the present application when executing the computer program. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1A It shows a schematic diagram of an application scenario corresponding to a method for analyzing the multi - dipole trajectory of a magnetocardiogram provided by an embodiment of the present application.

[0018] Figure 1B It shows a flowchart of the magnetocardiogram multipole trajectory analysis method provided by an embodiment of the present application.

[0019] Figure 1C It shows a multipole trajectory diagram provided by an embodiment of the present application.

[0020] Figure 1D It shows all isomagnetic maps of the T band provided by an embodiment of the present application.

[0021] Figure 2 It shows a flowchart for determining the multipole trajectory score in an embodiment of the present application.

[0022] Figure 3 It shows a flowchart for determining the first trajectory parameter in an embodiment of the present application.

[0023] Figure 4 It shows a flowchart for determining the overall trajectory parameter of the positive main branch in an embodiment of the present application.

[0024] Figure 5 It shows a flowchart for determining the third trajectory parameter of the main branch dipole in an embodiment of the present application.

[0025] Figure 6 It shows a flowchart for determining the fourth trajectory parameter of the secondary branch in an embodiment of the present application.

[0026] Figure 7A It shows a flowchart for determining the multipole trajectory score in an embodiment of the present application.

[0027] Figure 7B It shows another multipole trajectory diagram provided by an embodiment of the present application.

[0028] Figure 7C It shows a planar visualization diagram of the positive / negative main branch provided by an embodiment of the present application.

[0029] Figure 8 It shows a schematic diagram of the magnetocardiogram multipole trajectory analysis device provided by an embodiment of the present application.

[0030] Figure 9 It shows a schematic diagram of the structure of an electronic device in an embodiment of the present application.

[0031] Description of component labels

[0032] Steps S11 to S14

[0033] Steps S21 to S25

[0034] Steps S31 to S36

[0035] Steps S4 to S43

[0036] Steps S51 to S56

[0037] Steps S61 to S67

[0038] Steps S71 to S72

[0039] 80 Magnetocardiogram Multipole Trajectory Analysis Device

[0040] 81 Multipole Trajectory Map Acquisition Module

[0041] 82 Trajectory Branch Statistics Module

[0042] 83 Multipole Trajectory Score Determination Module

[0043] 84 Trajectory Analysis Result Determination Module

[0044] 90 Electronic Device

[0045] 91 Processor

[0046] 92 Non-volatile Storage Medium

[0047] 93 System Bus

[0048] 94 Internal Memory

[0049] 95 Network Interface Detailed Implementation Manner

[0050] The following uses specific specific examples to illustrate the implementation manners of the present application. Those skilled in the art can easily understand other advantages and effects of the present application from the content disclosed in this specification. The present application can also be implemented or applied through other different specific implementation manners. Various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of the present application. It should be noted that, without conflict, the following embodiments and the features in the embodiments can be combined with each other.

[0051] It should be noted that the drawings provided in the following embodiments only illustrate the basic concept of the present application in a schematic manner. Therefore, only the components related to the present application are shown in the drawings, rather than being drawn according to the number, shape, and size of the components in actual implementation. The form, quantity, and ratio of each component in actual implementation can be arbitrarily changed, and the component layout pattern may also be more complex.

[0052] In the prior art, during the analysis of magnetocardiogram multipole trajectories, there is a lack of an accurate magnetocardiogram multipole trajectory analysis method.

[0053] At least for the above problems, an embodiment of the present application provides a method for analyzing the multi - dipole trajectory of a magnetocardiogram. The method for analyzing the multi - dipole trajectory of a magnetocardiogram can obtain a multi - dipole trajectory map based on an original magnetocardiogram dataset; count multiple trajectory branches in the multi - dipole trajectory map, and divide the multiple trajectory branches into a positive main branch, a negative main branch, and secondary branches; determine a first trajectory parameter of the positive main branch based on the main branch information corresponding to the positive main branch, where the main branch information corresponding to the positive main branch includes the magnetic field value of the positive dipole on the positive main branch and the position information of the positive dipole; determine a second trajectory parameter of the negative main branch based on the main branch information corresponding to the negative main branch; match the positive dipole on the positive main branch and the negative dipole on the negative main branch within the same frame to obtain a main - branch dipole, and determine a third trajectory parameter of the main - branch dipole; determine a fourth trajectory parameter of the secondary branch based on the secondary - branch information corresponding to the secondary branch; determine a multi - dipole trajectory score based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter; determine the trajectory deviation degree between the original magnetocardiogram dataset and a preset magnetocardiogram dataset based on the multi - dipole trajectory score, and obtain the trajectory analysis result of the original magnetocardiogram dataset based on the trajectory deviation degree, which can solve the technical problem of the lack of accurate multi - dipole trajectory analysis in the prior art.

[0054] Figure 1A Shown is a schematic diagram of an application scenario corresponding to the method for analyzing the multi - dipole trajectory of a magnetocardiogram provided by an embodiment of the present application. As Figure 1A shown, this application scenario includes a multi - channel magnetocardiograph and an electronic device, and the multi - channel magnetocardiograph and the electronic device are communicatively connected. Among them, the multi - channel magnetocardiograph is used to collect an original magnetocardiogram dataset and send the collected original magnetocardiogram dataset to the electronic device. After receiving the original magnetocardiogram dataset sent by the multi - channel magnetocardiograph, the electronic device processes the data of the original magnetocardiogram dataset and draws a multi - dipole trajectory map of the original magnetocardiogram dataset, and respectively obtains a first trajectory parameter of the positive main branch, a second trajectory parameter of the negative main branch, a third trajectory parameter of the main - branch dipole, and a fourth trajectory parameter of the secondary branch based on the multi - dipole trajectory map, and determines a multi - dipole trajectory score based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter; determines the trajectory deviation degree between the original magnetocardiogram dataset and a preset magnetocardiogram dataset based on the multi - dipole trajectory score.

[0055] Next, the technical solutions in the embodiments of the present application will be described in detail with reference to the accompanying drawings in the embodiments of the present application.

[0056] Figure 1B Shown is a flowchart of the method for analyzing the multi - dipole trajectory of a magnetocardiogram provided by an embodiment of the present application. As Figure 1BAs shown, the method for analyzing the multi-pole trajectory of magnetocardiogram provided by the embodiment of the present application includes the following steps S11 to S14.

[0057] S11. Obtain a multi-pole trajectory map based on the original magnetocardiogram atlas.

[0058] In some embodiments, the obtaining of the multi-pole trajectory map based on the original magnetocardiogram atlas includes: obtaining the original magnetocardiogram atlas, determining the multi-pole information of each frame of at least one band according to the original magnetocardiogram atlas, wherein the multi-pole information is greater than a preset pole threshold; and drawing a multi-pole trajectory map based on the multi-pole information of each frame.

[0059] Exemplarily, the original magnetocardiogram atlas includes QR band, RS band, T band, etc. In actual applications, any suitable at least one band can be selected according to requirements, and the present application does not limit this.

[0060] Among them, the multi-pole includes a plurality of positive poles and a plurality of negative poles, and the multi-pole information includes the value of the pole and the position of the pole.

[0061] Exemplarily, the preset pole threshold can be 0.3. When determining the multi-pole information of each frame of at least one band according to the original magnetocardiogram atlas, the multi-pole information indicates that the value of the positive pole and the absolute value of the negative pole are both greater than 0.3.

[0062] It should be noted that the specific value of the preset pole threshold listed above is only for illustrative purposes. In actual applications, other any suitable preset pole thresholds can also be selected according to specific application requirements, and the present application does not limit this.

[0063] Exemplarily, the multi-pole trajectory map can refer to Figure 1C , where Figure 1C the black lines in are the positive main branch and the negative main branch respectively, and the remaining lines are secondary branches.

[0064] Exemplarily, Figure 1D shows all the isomagnetic maps of the T band provided by an embodiment of the present application. Each isomagnetic map contains multi-pole information (one or more positive poles, one or more negative poles).

[0065] It should be noted that the numbers and letters in the lower left corner of each of the above Figure 1D figures represent the labels of the pictures.

[0066] S12. Count the multiple trajectory branches in the multi-pole trajectory map, and divide the multiple trajectory branches into a positive main branch, a negative main branch, and secondary branches.

[0067] Among them, the positive main branch represents the branch with the largest number of positive poles, and all in the positive main branch are positive poles. The negative main branch represents the branch with the largest number of negative poles, and all in the negative main branch are negative poles. The secondary branch is the branch other than the positive main branch and the negative main branch, and the secondary branch includes a positive secondary branch and a negative secondary branch.

[0068] S13. Based on the positive main branch, the negative main branch and the secondary branch, determine the multi-pole trajectory score corresponding to the multi-pole trajectory map.

[0069] S14. Based on the multi-pole trajectory score, determine the trajectory deviation degree between the original magnetocardiogram atlas and the preset magnetocardiogram atlas, and based on the trajectory deviation degree, obtain the trajectory analysis result of the original magnetocardiogram atlas.

[0070] Exemplarily, the multi-pole trajectory score represents the trajectory deviation degree between the original magnetocardiogram atlas and the preset magnetocardiogram atlas. Among them, the larger the multi-pole trajectory score, the greater the trajectory deviation degree between the original magnetocardiogram atlas and the preset magnetocardiogram atlas.

[0071] Exemplarily, if the multi-pole trajectory score is in the range of 0 - 5, it indicates that the trajectory deviation degree between the original magnetocardiogram atlas and the preset magnetocardiogram atlas is small; if the multi-pole trajectory score is in the range of 5 - 10, it indicates that the trajectory deviation degree between the original magnetocardiogram atlas and the preset magnetocardiogram atlas is large. If the deviation degree is low, it can be considered that the original magnetocardiogram is relatively similar to the reference atlas. Therefore, the analysis result of the reference atlas can be used to analyze the original magnetocardiogram to obtain the analysis result of the original magnetocardiogram. Among them, the reference atlas is a magnetocardiogram with an analysis result set in advance. For example, taking the magnetocardiogram of a coronary heart disease patient as the reference map, if the trajectory deviation degree between the original magnetocardiogram and this reference map is low, it can be considered that the cardiac electrical activity reflected by the original magnetocardiogram may be abnormal and there may be a possibility of coronary heart disease.

[0072] Among them, the multi-pole trajectory scores in multiple frequency bands can also be calculated, and based on the sum of the multi-pole trajectory scores in multiple frequency bands, judge the trajectory deviation degree between the multi-pole trajectory of the magnetocardiogram and the reference atlas.

[0073] An embodiment of the present application provides a method for analyzing the trajectory of a magnetocardiogram multipole. In this method, by drawing a multipole trajectory diagram and statistically obtaining the positive main branch, negative main branch, and secondary branch, based on the positive main branch, negative main branch, and secondary branch, the multipole trajectory score corresponding to the multipole trajectory diagram is determined. Based on the multipole trajectory score, the trajectory deviation degree between the original magnetocardiogram set and the preset magnetocardiogram set is determined, and based on the trajectory deviation degree, the trajectory analysis result of the original magnetocardiogram set is obtained. It effectively depicts the details of the cardiac electrophysiological changes behind different wave branches, realizes the accurate information extraction of the time-space continuity, direction mutation, and energy distribution change of the multipole trajectory, can more accurately extract the information contained in the magnetocardiogram, and improves the accuracy of the magnetocardiogram analysis result.

[0074] Figure 2 It is shown as a flowchart for determining the multipole trajectory score in an embodiment of the present application. As Figure 2 shown, the process of determining the first trajectory parameter in the embodiment of the present application includes the following steps S21 to S25.

[0075] S21, determining the first trajectory parameter of the positive main branch based on the main branch information corresponding to the positive main branch.

[0076] Wherein, the main branch information corresponding to the positive main branch includes the magnetic field value of the positive pole on the positive main branch and the position information of the positive pole.

[0077] S22, determining the second trajectory parameter of the negative main branch based on the main branch information corresponding to the negative main branch.

[0078] Wherein, the main branch information corresponding to the negative main branch includes the magnetic field value of the negative pole on the negative main branch and the position information of the negative pole.

[0079] S23, matching the positive poles on the positive main branch and the negative poles on the negative main branch within the same frame to obtain the main branch dipole, and determining the third trajectory parameter of the main branch dipole.

[0080] Exemplarily, if the positive poles on the positive main branch and the negative poles on the negative main branch are both within the frame interval of 500 - 525, then the positive poles on the positive main branch and the negative poles on the negative main branch within the frame interval of 500 - 525 are matched as the main branch dipole. Conversely, if a certain positive pole on the positive main branch is within the frame interval of 500 - 525 and a certain negative pole on the main branch is within the frame interval of 55 - 550, then at this time, the positive pole and the negative pole are not within the same frame interval, and the positive pole and the negative pole cannot be matched. After matching the positive poles on the positive main branch and the negative poles on the negative main branch within the same frame to obtain the main branch dipole, determine the third trajectory parameter of the main branch dipole.

[0081] S24. Determine the fourth trajectory parameter of the sub-branch based on the sub-branch information corresponding to the sub-branch.

[0082] S25. Determine the multipole trajectory score based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter.

[0083] Among them, the multipole trajectory score includes the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the sub-branch trajectory score.

[0084] Specifically, obtain the positive pole trajectory score based on the first trajectory parameter, obtain the negative pole trajectory score based on the second trajectory parameter, obtain the dipole trajectory score based on the third trajectory parameter, and obtain the sub-branch trajectory score based on the fourth trajectory parameter.

[0085] Exemplarily, taking the T band as an example, the expression of the multipole trajectory score is:

[0086] TTtjsc = m1 * mpsc + m2 * mnsc + m3 * mdsc + m4 * srsc

[0087] Among them, mi represents a preset coefficient, where i = 1, 2, 3, 4, mpsc represents the positive pole trajectory score, mnsc represents the negative pole trajectory score, mdsc represents the dipole trajectory score, and srsc represents the sub-branch trajectory score.

[0088] The embodiment of the present application provides a method for determining the multipole trajectory score. This method determines the first trajectory parameter corresponding to the positive main branch, the second trajectory parameter corresponding to the sub-branch, the fourth trajectory parameter corresponding to the sub-branch, and matches the positive poles on the positive main branch and the negative poles on the negative main branch to obtain the main branch dipole, and determines the third trajectory parameter of the main branch dipole, and obtains an accurate multipole trajectory score based on the above four trajectory parameters. By combining the four trajectory parameters, this multipole trajectory scoring method can comprehensively reflect the dynamic characteristics of the multipole trajectory, overcome the problem of one-sidedness that may be caused by single-parameter analysis, and greatly improve the accuracy of magnetocardiogram multipole trajectory analysis.

[0089] Figure 3 It is shown as the flowchart for determining the first trajectory parameter in an embodiment of the present application. As Figure 3 shown, the process of determining the first trajectory parameter in the embodiment of the present application includes the following steps S31 to S36.

[0090] S31. Obtain the magnetic field value and the position information of the positive poles on the positive main branch that are greater than the preset pole threshold.

[0091] Exemplarily, the preset dipole threshold is any suitable value such as 0.3, 0.5, etc., and the present application does not limit this.

[0092] Exemplarily, the magnetic field value of the positive dipole can be detected based on a gaussmeter, a superconducting quantum interference device, and a fluxgate magnetometer.

[0093] It should be noted that the instruments listed above for detecting the magnetic field value of the positive dipole are only for illustrative purposes. In actual applications, other any suitable instruments or devices can also be selected according to specific application requirements to detect the magnetic field value of the positive dipole, and the present application does not limit this.

[0094] Among them, the position information of the positive dipole can be coordinate position information.

[0095] Exemplarily, the position information of the positive dipole can be obtained based on a multi - dipole trajectory diagram.

[0096] Exemplarily, the positive dipole can also be imaged based on an electron microscope or an X - ray imager. Through image analysis software, combined with the magnetic field measurement results, the specific positions of the positive dipoles greater than the preset threshold are located, and the coordinate information of the positive dipoles is recorded.

[0097] S32. Based on the magnetic field value and the position information, determine the overall trajectory parameters of the positive main branch.

[0098] Among them, the overall trajectory parameters include effective length, main dipole length, and non - boundary rate.

[0099] S33. Sequentially determine the planar distance between adjacent frames of positive dipoles on the positive main branch, and determine a section on the positive main branch where the planar distance is less than the preset offset distance as the main section.

[0100] Specifically, count the planar distance between the front and rear frames of positive dipoles on the positive main branch, and take the largest stable section in the middle of the positive main branch as the main section (m).

[0101] S34. Determine a section on the positive main branch in front of the main section as the head.

[0102] S35. Determine a section on the positive main branch behind the main section as the tail.

[0103] S36. Respectively determine the sectional trajectory parameters of the entire section, the main section, the head, and the tail. The sectional trajectory parameters include trajectory path, trajectory span, and trajectory elastic curvature.

[0104] Among them, the entire section represents the entire positive main branch.

[0105] Specifically, sum up the planar distances of all positive dipoles on the entire section, the sectional section, the head, and the tail respectively to obtain the trajectory path (tjpa).

[0106] Specifically, the farthest two-point distances of all the positive poles on the whole section, segmented section, head, and tail are respectively determined as the trajectory span (tjsp).

[0107] Specifically, the trajectory path / trajectory span of the whole section, segmented section, head, and tail are respectively obtained to get the trajectory elastic curvature (tjec).

[0108] It should be noted that the specific steps of the method for determining the first trajectory parameter of the negative poles on the negative main branch are similar to the specific steps of the method for determining the first trajectory parameter of the positive poles on the positive main branch as described above, and thus will not be elaborated herein in this application.

[0109] The embodiments of this application provide a method for determining the first trajectory parameter. In this method, the first trajectory parameter includes the overall trajectory parameter and the segmented trajectory parameter. The overall trajectory parameter is determined through the magnetic field values and position information of the positive poles on the positive main branch. The planar distances between adjacent frames of the positive poles on the positive main branch are sequentially determined, and a section on the positive main branch where the planar distance is less than the preset offset distance is determined as the main section; a section on the positive main branch in front of the main section is determined as the head; a section on the positive main branch behind the main section is determined as the tail; the segmented trajectory parameters of the main section, head, and tail are respectively determined, and the segmented trajectory parameters include the trajectory path, trajectory span, and trajectory elastic curvature. By segmenting the positive main branch, the complex overall structure can be decomposed into simpler and more consistent local units (such as the head, main section, and tail). This method for segmenting and determining the first trajectory parameter enables the numerical values of the calculated trajectory path, trajectory span, and trajectory elastic curvature to more accurately reflect the characteristics or performance of different sections in the positive main branch, facilitating in-depth analysis and understanding of the structures and characteristics of different sections in the positive main branch, and providing an accurate data basis for subsequent trajectory analysis of the original magnetocardiogram.

[0110] Figure 4 Shown is a flowchart of determining the overall trajectory parameter of the positive main branch in an embodiment of this application. As Figure 4 shown, the process of determining the overall trajectory parameter of the positive main branch in the embodiments of this application includes the following steps S41 to S43.

[0111] S41, based on the magnetic field values of all the positive poles greater than the preset pole threshold, determine the effective length corresponding to the positive main branch.

[0112] Specifically, based on the normalization method, the magnetic field values of all the positive poles are sequentially divided by the maximum magnetic field value, so that all the positive poles on the positive main branch are limited between 0 and 1.

[0113] Exemplarily, the preset dipole threshold can be set to 0.3 - 1 or 0.5 - 1. If the preset dipole threshold is set to 0.3 - 1, then the magnetic field values of the positive dipoles on the positive main branch that are within 0.3 - 1 are screened, and the effective length corresponding to the positive main branch is determined based on all the positive dipoles within 0.3 - 1.

[0114] Exemplarily, if the magnetic field values of the positive dipoles on the positive main branch are [1, 1, 1, 1, 0.7, 1] at a total of 6 points, its effective length is 5.7 / 6 * 100 = 95.

[0115] Similarly, based on the normalization method, the negative magnetic field values of all the negative dipoles are successively divided by the absolute value of the maximum negative magnetic field value, so that all the negative dipoles on the negative main branch are limited between -1 and 0. At this time, the preset dipole threshold can be set to -1 to -0.3 or -1 to -0.5, and the effective length corresponding to the negative main branch is determined based on all the negative dipoles within the preset dipole threshold.

[0116] It should be noted that the various numerical values listed above are only for illustrative purposes. In actual applications, any appropriate numerical value can be selected according to specific application requirements, and the present application does not limit this.

[0117] Specifically, the magnetic field values of all the positive dipoles on the positive main branch are summed and then averaged, and the averaged magnetic field value is rounded to obtain the effective length (melg) of the positive main branch. For example, rounding the averaged magnetic field value can be multiplying the averaged magnetic field value by 100, 1000, etc.

[0118] It should be noted that multiplying the averaged magnetic field value by 100 or 1000 above is only for illustrative purposes. In actual applications, other arbitrary appropriate numerical values can also be selected to round the averaged magnetic field value, and the present application does not limit this.

[0119] S42. Based on the proportion of the positive dipole with the maximum magnetic field value among all the positive dipoles on the positive main branch that are greater than the preset dipole threshold, the main dipole length corresponding to the positive main branch is obtained.

[0120] Among them, the positive dipole with the maximum magnetic field value refers to the positive dipole with a magnetic field value of 1.

[0121] Specifically, determine the proportion of the positive dipoles with a magnetic field value of 1 on the positive main branch among all the positive dipoles on the positive main branch that are greater than the preset dipole threshold, and multiply this proportion by 100 or 1000, etc., to obtain the main dipole length (mmlg).

[0122] Exemplarily, if there are 6 positive poles on a positive main branch, and the magnetic field values of the 6 positive poles are [1, 1, 1, 1, 0.7, 1] respectively, the length of its main pole is 5 / 6 * 100 = 83.

[0123] Another exemplarily, if there are 6 positive poles on a positive main branch, and the magnetic field values of the 6 positive poles are [1, 1, 1, 1, 0.7, 1] at a total of 6 points, the length of its main pole is 5 / 6 * 1000 = 830.

[0124] It should be noted that the length of the main pole should be an integer.

[0125] It should be noted that during the process of taking the length of the main pole as an integer, this ratio can be multiplied by 100 or 1000 or any other appropriate value, and the present application does not limit this.

[0126] S43, based on the proportion of the planar coordinates of the positive poles on the positive main branch that are not within the preset range of the preset matrix boundary, obtain the non-boundary rate corresponding to the positive main branch.

[0127] Exemplarily, based on the proportion of the planar coordinates of the poles on the main branch that are not within 5 of the boundary of the 100 * 100 matrix, determine the non-boundary rate (mnbd).

[0128] Another exemplarily, based on the proportion of the planar coordinates of the poles on the main branch that are not within 2 of the boundary of the 10 * 10 matrix, determine the non-boundary rate (mnbd).

[0129] It should be noted that the specific values such as the 100 * 100 matrix, 10 * 10 matrix, within 5 of the matrix boundary, and within 2 of the matrix boundary listed above are only used for exemplary illustration. In actual applications, other arbitrarily appropriate matrices and the corresponding values of the matrix boundary can be selected according to specific application requirements, and the present application does not limit this.

[0130] It should be noted that the method for determining the overall trajectory parameters of the negative main branch is similar to the steps of the method for determining the overall trajectory parameters of the positive main branch as described above, and the present application will not elaborate on this.

[0131] An embodiment of the present application provides a method for determining the overall trajectory parameters of a positive main branch. In this method, based on the magnetic field values of all the positive poles greater than a preset pole threshold, the effective length corresponding to the positive main branch is determined; based on the ratio of the magnetic field value equal to the maximum magnetic field value, the main pole length corresponding to the positive main branch is obtained; based on the proportion of the planar coordinates of the positive poles on the positive main branch that are not within the preset range of the preset matrix boundary, the non-boundary rate corresponding to the positive main branch is obtained. The effective length, main pole length, and non-boundary rate of the positive main branch can be obtained simply and intuitively without a complex calculation process, greatly improving the calculation efficiency and accuracy of the effective length, main pole length, and non-boundary rate of the positive main branch, and the effective length, main pole length, and non-boundary rate of the positive main branch accurately reflect the structural characteristics of the positive main branch.

[0132] Figure 5 It shows a flowchart for determining the third trajectory parameter of the main branch dipole in an embodiment of the present application. As Figure 4 shown, the process of determining the third trajectory parameter of the main branch dipole in the embodiment of the present application includes the following steps S51 to S56.

[0133] S51, determine the intersection part of the positive poles on the positive main branch and the negative poles on the negative main branch as the main branch of the main branch dipole.

[0134] Specifically, the frame interval of each positive pole greater than the preset pole threshold on the positive main branch is determined in sequence, and the frame interval of each negative pole greater than the preset pole threshold on the negative main branch is determined in sequence. Among them, the absolute value of the magnetic field value of each negative pole is compared with the preset pole threshold, and the positive poles and negative poles located in the same frame interval are matched as the main branch dipole, and all the main branch dipoles form the main branch (br) of the main branch dipole.

[0135] Exemplarily, the frame interval of each positive pole greater than the preset pole threshold on the positive main branch and the frame interval of each negative pole greater than the preset pole threshold on the negative main branch can be obtained through the t-axis in the multi-pole trajectory diagram.

[0136] S52, determine the planar distance between the positive poles of adjacent frames on the positive main branch in sequence, and determine a section on the positive main branch where the planar distance is less than the preset offset distance as the main section of the positive main branch.

[0137] Exemplarily, the preset offset distance can be set to any appropriate value such as 4, 6, etc.

[0138] Exemplarily, the position information of the positive pole includes coordinate information. Based on the Pythagorean theorem, the coordinate information of the positive poles of adjacent frames on the positive main branch is calculated to obtain the planar distance between the positive poles of adjacent frames on the positive main branch.

[0139] Exemplarily, if the planar distance between adjacent positive poles from a certain point A to a certain point B on the positive main branch is less than 4, but the planar distance between adjacent positive poles in the adjacent frames after point B is a value greater than or equal to 4, then the section from A to B on the positive main branch is determined as the main section.

[0140] It should be noted that in the process of determining the main section based on the above method, there is no restriction on the length of the main section, as long as the planar distance between adjacent positive poles on the positive main branch meets the requirement that the planar distance is less than the preset offset distance.

[0141] S53. Sequentially determine the planar distances between adjacent negative poles on the negative main branch, and determine a section on the negative main branch where the planar distance is less than the preset offset distance as the main section of the negative main branch.

[0142] S54. Determine the intersection part of the main section of the positive main branch and the main section of the negative main branch as the main section of the main branch dipole.

[0143] Specifically, the frame interval corresponding to the main section of the positive main branch can be read based on the multi-pole trajectory diagram, and the frame interval corresponding to the main section of the negative main branch can be determined based on the multi-pole trajectory diagram. The intersection part of the frame interval corresponding to the main section of the positive main branch and the frame interval corresponding to the main section of the negative main branch is determined as the main section (sg) of the main branch dipole.

[0144] S55. Take the difference between the main branch of the main branch dipole and the main section of the main branch dipole to obtain the secondary section of the main branch dipole.

[0145] Specifically, subtract the time point of the main section from the time point of the main branch, and the remaining frames are the secondary section (df) of the main branch dipole.

[0146] S56. Based on the main branch information, sequentially obtain the third trajectory parameters of the main branch, main section, and secondary section of the main branch dipole.

[0147] Specifically, the third trajectory parameters include 13 kinds of trajectory parameter values, specifically including effective length (melg), non-boundary rate (mnbd), overall rotation angle (tta), quadrant where it is located (q), amplitude (areaf), pole moment (aread), maximum dipole moment (ma), minimum value (mi), mean value (av), mean square deviation (ms), trajectory path (tjpa), trajectory span (tjpa), and trajectory elastic curvature (tjec).

[0148] Exemplarily, the effective lengths of all the positive main branches and all the negative main branches are determined respectively; the effective lengths of all the positive main branches and all the negative main branches are averaged to obtain the effective length of the main branch dipole; wherein, the method or steps for determining the effective length of the positive main branch have been described in detail in S31 above, and the present application will not repeat them here.

[0149] Exemplarily, the non-boundary rates of all the positive main branches and all the negative main branches are determined respectively; the non-boundary rates of all the positive main branches and all the negative main branches are averaged to obtain the non-boundary rate of the main branch dipole, wherein the specific steps for determining the non-boundary rate of the positive main branch have been described in detail in step S33 above, and the present application will not repeat them here.

[0150] Exemplarily, the positive and negative pointing angle represents the angular vector of the positive pole pointing to the negative pole in the main branch dipole on the plane, that is, an angular value. The amplitude difference represents the magnetic field value of the positive pole minus the magnetic field value of the negative pole among the positive and negative poles. The dipole moment represents the distance between the positive and negative poles. The position of the center point of the pole can be obtained by taking the average of the sum of the coordinate information of the positive and negative poles in the main branch dipole.

[0151] Exemplarily, within the entire time interval, the angular value of the overall rotation of the positive and negative pointing angles is obtained. The difference is taken between the front and back positive and negative pointing angles, and then the sum of the differences of different positive and negative pointing angles is calculated as a whole, which is the overall rotation angle. And record the specific quadrants where the positive and negative pointing angles fall during the rotation. The specific quadrants where the positive and negative pointing angles fall during the rotation can be directly read from the multi-pole trajectory diagram, where Figure 1D only the multi-pole trajectory diagram of the first quadrant is listed.

[0152] Exemplarily, the positive and negative pointing angles of the main branch, main segment, and secondary segment of the main branch dipole are respectively adjusted to the same starting point (the positive poles coincide), and the area swept during the rotation is calculated with the normalized amplitude difference / dipole moment as the vector length respectively.

[0153] Exemplarily, based on the sequence of all dipole moments on the main branch of the main branch dipole, the maximum value, minimum value, mean value, and mean square deviation of the dipole moment of the main branch of the main branch dipole are obtained.

[0154] Similarly, based on the sequence of all dipole moments on the main segment of the main branch dipole, the maximum value, minimum value, mean value, and mean square deviation of the dipole moment of the main segment of the main branch dipole are obtained.

[0155] Similarly, based on the sequence of all dipole moments on the secondary segment of the main branch dipole, the maximum value, minimum value, mean value, and mean square deviation of the dipole moment of the secondary segment of the main branch dipole are obtained.

[0156] It should be noted that the specific methods for determining the trajectory path, trajectory span, and trajectory elastic curvature in the third trajectory parameter are similar to the steps for determining the trajectory path, trajectory span, and trajectory elastic curvature in S36 above, and are not elaborated in this application.

[0157] The embodiment of the present application provides a method for determining the third trajectory parameter. In this method, the main branch of the main branch dipole is determined based on the positive pole on the positive main branch and the negative pole on the negative main branch; and the main segment of the positive main branch is determined based on the relationship between the planar distance of the positive poles of adjacent frames on the positive main branch and the preset offset distance. Similarly, the main segment of the negative main branch is determined based on the relationship between the planar distance of the negative poles of adjacent frames on the negative main branch and the preset offset distance, and the main segment of the branch dipole and the secondary segment of the main branch dipole are determined, and the third trajectory parameters of the main branch of the main branch dipole, the main segment of the main branch dipole, and the secondary segment of the main branch dipole are determined in sequence based on the main branch information, where each segment includes 13 third trajectory parameters. The 13 third trajectory parameters provide rich data, which helps to more comprehensively understand the dynamic changes and data characteristics of the main branch dipole.

[0158] Figure 6 Shown is a flowchart for determining the fourth trajectory parameter of the secondary branch in an embodiment of the present application. Among them, the fourth trajectory parameter includes the effective number of branches, branch size, total trajectory path, total trajectory span, total trajectory elastic curvature, and total trajectory elastic strength. As Figure 6 shown, the process of determining the fourth trajectory parameter of the secondary branch in the embodiment of the present application includes the following steps S61 to S67.

[0159] S61, determine the non-boundary rate of each of the secondary branches, and remove the secondary branches with a non-boundary rate less than or equal to the preset boundary threshold to obtain the boundary factor of the secondary branches.

[0160] Among them, the secondary branches include positive secondary branches and negative secondary branches, and the determination of the fourth trajectory parameter of the secondary branches is calculated separately for the positive secondary branches and the negative secondary branches.

[0161] Exemplarily, the preset boundary threshold can be set to 0.1, 0.2, etc.

[0162] It should be noted that the specific values of the above preset boundary threshold are only for exemplary illustration, and in actual applications, other any appropriate preset boundary threshold can also be selected according to specific application requirements, and this application does not limit this.

[0163] S62, obtain the total magnetic field value between the start time and the end time of all secondary branches.

[0164] Exemplarily, taking the positive secondary branch as an example, the start time and end time of all secondary branches are denoted as t1 and t2 respectively, and the sum of all magnetic field values between t1 and t2 is determined as the total magnetic field value (amf).

[0165] S63. Divide the magnetic field value of each said secondary branch by the total magnetic field value to obtain the length factor of each said secondary branch.

[0166] S64. Sum all the said boundary factors to obtain the effective number of branches.

[0167] Among them, the larger the effective number of branches, the closer it is to the boundary and the more likely it is to be noise.

[0168] S65. Multiply each said boundary factor by the length factor of the corresponding said secondary branch in sequence and then sum them to obtain the branch size of the said secondary branch.

[0169] S66. Determine the trajectory path, trajectory span, and trajectory elastic curvature of each said secondary branch, and obtain the trajectory elastic strength based on the trajectory elastic curvature.

[0170] It should be noted that the specific method or steps for determining the trajectory path, trajectory span, and trajectory elastic curvature of the secondary branch are similar to those in S36 for determining the trajectory path, trajectory span, and trajectory elastic curvature above, and this application will not elaborate on them here.

[0171] Exemplarily, the trajectory elastic strength represents the square of the trajectory elastic curvature.

[0172] S67. Sum the trajectory path, trajectory span, trajectory elastic curvature, and trajectory elastic strength of all the said secondary branches respectively to obtain the total trajectory path, total trajectory span, total trajectory elastic curvature, and total trajectory elastic strength.

[0173] It should be noted that the method for determining the fourth trajectory parameter of the secondary branch is applicable to determining the fourth trajectory parameter of the positive secondary branch and the negative secondary branch.

[0174] The embodiment of the present application provides a method for determining the fourth trajectory parameter of a secondary branch. In this method, the fourth trajectory parameter includes the effective number of branches, branch size, total trajectory path, total trajectory span, total trajectory elastic curvature, and total trajectory elastic strength. By sequentially determining the above 6 fourth trajectory parameters, a detailed description of the structural characteristics of the secondary branch is realized. The above 6 fourth trajectory parameters reveal the geometric shape and spatial distribution of the secondary branch, providing an accurate and rich data basis for subsequent trajectory analysis of the original magnetocardiogram.

[0175] Figure 7A Shown is a flowchart for determining the multi-pole trajectory score in an embodiment of the present application. AsFigure 7A As shown in Figure 7A , the process of determining the multipole trajectory score in the embodiments of the present application includes the following steps S71 to S72.

[0176] S71. Based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter, determine the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the sub-branch trajectory score respectively.

[0177] Exemplarily, before determining the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the sub-branch trajectory score, a data set of 1000 samples can be first constructed and labeled according to the analysis results of the preset magnetocardiogram atlas. The magnetocardiogram data sets of each sample are obtained respectively, the multipole trajectories of different bands are plotted, and the multipole trajectory branch data of different bands are obtained; then, using the analysis methods in the above steps, the trajectory parameters of the main branch, the trajectory parameters of the main branch dipole, and the trajectory parameters of the sub-branch are extracted, and according to the distribution of different trajectory parameters of all samples, the value range of the weight parameters for calculating each trajectory score is delimited. Thereafter, the 4 kinds of trajectory scores of the magnetocardiogram to be analyzed and the trajectory deviation degree between the magnetocardiogram to be analyzed and the preset magnetocardiogram atlas (T-band multipole trajectory score / sum of multipole trajectory scores of multiple bands) can be calculated according to the different trajectory parameters of the magnetocardiogram to be analyzed.

[0178] Specifically, the expression for determining the positive pole trajectory score is:

[0179] mpsc = max(a11×p_melg, a12×p_mmlg, a13×p_mnbd×100)

[0180] + [min(a21×p_atjpa÷100, a22×p_atjsp÷10) + a23×p_atjec÷10]

[0181] + [min(a31×p_mtjpa÷100, a32×p_mtjsp÷10) + a33×p_mtjec÷10]

[0182] + [min(a41×p_htjpa÷100, a42×p_htjsp÷10) + a43×p_htjec÷10]

[0183] + [min(a51×p_rtjpa÷100, a52×p_rtjsp÷10) + a53×p_rtjec÷10]

[0184] where aij is a preset coefficient (i = 1 to 5, j = 1 to 5), and some parameters are multiplied by 100 or divided by 100 or divided by 10, which is determined by the interpolation dimension of the isomagnetic map (taking the 100×100 isomagnetic map as an example here, p_melg represents the effective length of the positive pole, p_mmlg represents the main pole length of the positive pole, p_mnbd represents the non-boundary rate of the positive pole, p_atjpa represents the trajectory path of the main segment of the positive pole, p_atjsp represents the trajectory span of the main segment of the positive pole, p_atjec represents the trajectory elastic curvature of the main segment of the positive pole, p_mtjpa represents the trajectory path of the main segment of the positive pole, p_mtjsp represents the trajectory span of the main segment of the positive pole, p_mtjec represents the trajectory elastic curvature of the main segment of the positive pole, p_htjpa represents the trajectory path of the head of the positive pole, p_htjsp represents the trajectory span of the head of the positive pole, p_htjec represents the trajectory elastic curvature of the head of the positive pole, p_rtjpa represents the trajectory path of the tail of the positive pole, p_rtjsp represents the trajectory span of the tail of the positive pole, p_rtjec represents the trajectory elastic curvature of the tail of the positive pole).

[0185] Specifically, the expression for determining the trajectory score of the negative pole is as follows:

[0186] mnsc = max(b11×n_melg, b12×n_mmlg, b13×n_mnbd×100)

[0187] + [min(b21×n_atjpa÷100, b22×n_atjsp÷10) + b23×n_atjec÷10]

[0188] + [min(b31×n_mtjpa÷100, b32×n_mtjsp÷10) + b33×n_mtjec÷10]

[0189] + [min(b41×n_htjpa÷100, b42×n_htjsp÷10) + b43×n_htjec÷10]

[0190] + [min(b51×n_rtjpa÷100, b52×n_rtjsp÷10) + b53×n_rtjec÷10]

[0191] Among them, bij is a preset coefficient (i = 1 to 5, j = 1 to 5), n_melg represents the effective length of the positive pole, n_mmlg represents the main pole length of the positive pole, n_mnbd represents the non-boundary rate of the positive pole, n_atjna represents the trajectory path of the main segment of the positive pole, n_atjsn represents the trajectory span of the main segment of the positive pole, n_atjec represents the trajectory elastic curvature of the main segment of the positive pole, n_mtjna represents the trajectory path of the main segment of the positive pole, n_mtjsn represents the trajectory span of the main segment of the positive pole, n_mtjec represents the trajectory elastic curvature of the main segment of the positive pole, n_htjna represents the trajectory path of the head of the positive pole, n_htjsn represents the trajectory span of the head of the positive pole, n_htjec represents the trajectory elastic curvature of the head of the positive pole, n_rtjna represents the trajectory path of the tail of the positive pole, n_rtjsn represents the trajectory span of the tail of the positive pole, n_rtjec represents the trajectory elastic curvature of the tail of the positive pole.

[0192] Specifically, the expression for determining the dipole trajectory score is:

[0193] mdsc = [max(c10×br_melg, c11×br_mnbd×100) + (c12×br_tta + c13×br_q)

[0194] + 10×(c14×br_areaf + c15×br_aread)

[0195] + min(c16×br_didma÷50, c17×br_didmi÷50, c18×br_didav÷50, c19×br_didms)

[0196] + min(f11×br_ctjpa÷100, f12×br_ctjsp÷10) + f13×br_ctjec÷10]

[0197] + [max(c20×sg_melg, c21×sg_mnbd×100) + (c22×sg_tta + c23×sg_q)

[0198] + 10×(c24×sg_areaf + c25×sg_aread)

[0199] + min(c26×sg_didma÷50, c27×sg_didmi÷50, c28×sg_didav÷50, c29×sg_didms)

[0200] +min(f21×sg_ctjpa÷100,f22×sg_ctjsp÷10)+f23×sg_ctjec÷10

[0201] +[max(c30×df_melg,c31×df_mnbd×100)+(c32×df_tta+c33×df_q)

[0202] +10×(c34×df_areaf+c35×df_aread)

[0203] +min(c36×df_didma÷50,c37×df_didmi÷50,c38×df_didav÷50,c39×df_didms)

[0204] +min(f31×df_ctjpa÷100,f32×df_ctjsp÷10)+f33×df_ctjec÷10

[0205] Wherein, cij and fij are preset coefficients, and c13, c23, and c33 are classification activation coefficients of quadrants (taking the T band as an example, the coefficients are different when the quadrant is 1 / 12 / 14, and the coefficients are the same for other values, which is determined by the parameter distribution of different bands); br represents the main branch, sg represents the main segment, did represents the dipole, c represents the center point of the dipole, and df represents the secondary segment.

[0206] Specifically, the expression for determining the score of the secondary branch trajectory is:

[0207] srsc = [max(d11×p_snbd×10,d12×p_selg×10)

[0208] +min(d13×p_sapa÷100,d14×p_sasp÷10,d15×p_saec÷10,d16×p_saep÷100)]

[0209] +[max(d21×n_snbd×10,d22×n_selg×10)

[0210] +min(d23×n_sapa÷100,d24×n_sasp÷10,d25×n_saec÷10,d26×n_saep÷100)]

[0211] Wherein, dij are preset coefficients.

[0212] S72. Obtain the multipole trajectory score based on the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the secondary branch trajectory score.

[0213] Exemplarily, taking the T band as an example, the expression for obtaining the multipole trajectory score based on the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the sub-branch trajectory score is:

[0214] TTtjsc = m1 × mpsc + m2 × mnsc + m3 × mdsc + m4 × srsc

[0215] where mi is a preset coefficient, where i = 1 to 4.

[0216] It should be noted that the sum of the multipole trajectory scores of multiple bands (such as the QR band, the RS band, and the T band) can also be calculated, and the trajectory deviation degree between the magnetocardiogram multipole trajectory and the reference atlas can be judged based on the sum of the multipole trajectory scores of the calculated multiple bands (such as the QR band, the RS band, and the T band). The specific number of calculation bands can be determined according to specific application requirements, and this application does not limit this.

[0217] In addition, according to actual needs, only one or more trajectory scores can be used to calculate the trajectory deviation degree between it and the preset magnetocardiogram atlas, that is, one or at most 3 of the preset coefficients in the above expression of the multipole trajectory score are set to 0.

[0218] The embodiment of the present application provides a method for determining the multipole trajectory score. In this method, the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the sub-branch trajectory score are respectively determined by the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter calculated through the above respective steps. The positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the sub-branch trajectory score are input into the expression of the multipole trajectory score to obtain the multipole trajectory score. Among them, the processes of determining the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score, and the sub-branch trajectory score and the process of determining the multipole trajectory score are all calculated through specific formulas, which greatly improves the accuracy of the calculation results. The formula-based calculation eliminates the possibility of subjective judgment and human bias. When it is necessary to explain a certain score result, it is possible to trace back to the specific formula and parameters to analyze which factors have a major impact on the score, which helps to understand the system behavior and has strong interpretability; the formula can be easily integrated into software or algorithms to realize the automation of the scoring process without manual evaluation one by one, which greatly improves the efficiency of processing a large amount of data.

[0219] The protection scope of the magnetocardiogram multipole trajectory analysis method in the embodiments of the present application is not limited to the execution order of the steps listed in this embodiment. Any solution achieved by adding or subtracting steps of the prior art and replacing steps according to the principle of the present application is included in the protection scope of the present application.

[0220] In addition, in a magnetocardiogram multipole trajectory analysis method provided by an embodiment of the present application, several characteristic parameters can also be extracted for parameter analysis, trajectory analysis, and score calculation. For example, specifically extract the trajectory path, trajectory span, and trajectory elastic curvature parameters of the head segment of the positive pole in the T band, draw the multipole trajectory of this segment, and design a new score calculation method to study the characteristic laws of the information of this segment.

[0221] Please refer to Figure 7B , Figure 7B which shows another multipole trajectory diagram provided by an embodiment of the present application. Figure 7B It shows the 3D image of the multipole trajectory in the T band in an embodiment (corresponding to the positions of the multipoles in the isomagnetic map above Figure 1D . The t-axis is the time point, and the x / y axes correspond to the rows and columns of the isomagnetic map. In Figure 7B , the branch in the lower left corner represents the positive secondary branch, the branch at the top represents the negative secondary branch, the longest branch in the middle represents the positive main branch, and the branch on the rightmost side represents the negative main branch.

[0222] Please refer to Figure 7C , Figure 7C which shows the plane visualization diagram of the positive / negative main branches provided by an embodiment of the present application.

[0223] The embodiments of the present application also provide a magnetocardiogram multipole trajectory analysis device. The magnetocardiogram multipole trajectory analysis device can implement the magnetocardiogram multipole trajectory analysis method of the present application. However, the implementation devices of the magnetocardiogram multipole trajectory analysis method of the present application include, but are not limited to, the structure of the magnetocardiogram multipole trajectory analysis device listed in this embodiment. Any structural deformation and replacement of the prior art made according to the principle of the present application are included in the protection scope of the present application.

[0224] As Figure 8 shown, in an embodiment, the magnetocardiogram multipole trajectory analysis device 80 of the present application includes a multipole trajectory diagram acquisition module 81, a trajectory branch statistics module 82, a multipole trajectory score determination module 83, and a trajectory analysis result determination module 84.

[0225] The multipole trajectory map acquisition module 81 is used to obtain a multipole trajectory map based on the original magnetocardiogram atlas; the trajectory branch statistics module 82 is used to count multiple trajectory branches in the multipole trajectory map and divide the multiple trajectory branches into a positive main branch, a negative main branch, and secondary branches; the multipole trajectory score determination module 83 is used to determine a multipole trajectory score corresponding to the multipole trajectory map based on the positive main branch, the negative main branch, and the secondary branches; the trajectory analysis result determination module 84 is used to determine a trajectory deviation degree between the original magnetocardiogram atlas and a preset magnetocardiogram atlas based on the multipole trajectory score, and obtain a trajectory analysis result of the original magnetocardiogram based on the trajectory deviation degree.

[0226] Among them, the structure and principle of the multipole trajectory map acquisition module 81, the trajectory branch statistics module 82, the multipole trajectory score determination module 83, and the trajectory analysis result determination module 84 correspond one by one to the steps in the above magnetocardiogram multipole trajectory analysis method, so they will not be elaborated here.

[0227] In several embodiments provided in the present application, it should be understood that the disclosed device or method can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of modules / units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple modules or units can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed mutual coupling or direct coupling or communication connection can be through some interfaces. The indirect coupling or communication connection of the device or module or unit can be in an electrical, mechanical or other form.

[0228] The module / unit described as a separated component may or may not be physically separated, and the component displayed as a module / unit may or may not be a physical module, that is, it may be located in one place, or may be distributed to multiple network units. Some or all of the modules / units can be selected according to actual needs to achieve the purpose of the embodiments of the present application. For example, in each embodiment of the present application, each functional module / unit can be integrated in a processing module, or each module / unit can exist physically alone, or two or more modules / units can be integrated in one module / unit.

[0229] Those of ordinary skill in the art should also be able to further realize that the units and steps of each example described in combination with the embodiments disclosed herein can be implemented by electronic hardware, computer software, or a combination of the two. To clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described according to functions in the above description. Whether these functions are executed in a hardware or software manner depends on the specific application and design constraints of the technical solution. Professional technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of this application.

[0230] The embodiments of this application also provide a computer-readable storage medium. Those of ordinary skill in the art can understand that all or part of the steps in the methods of the above embodiments can be completed by instructing a processor through a program. The program can be stored in a computer-readable storage medium. The storage medium is a non-transitory medium, such as random access memory, read-only memory, flash memory, hard disk, solid-state drive, magnetic tape, floppy disk, optical disc, and any combination thereof. The above storage medium can be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more integrated available media. The available medium can be a magnetic medium (such as a floppy disk, hard disk, magnetic tape), an optical medium (such as a digital video disc (DVD)), or a semiconductor medium (such as a solid-state disk (SSD)), etc.

[0231] The embodiments of this application also provide an electronic device. Figure 9 Shown is a schematic structural diagram of an electronic device 90 in an embodiment of this application. The feature vector extraction method provided by the embodiments of this application can be applied to Figure 9 the shown electronic device 90, but is not limited thereto. As Figure 9 shown, the electronic device 90 includes a processor 91, a memory, a system bus 93, and a network interface 95. Among them, the memory may include a non-volatile storage medium 92 and an internal memory 94.

[0232] The non-volatile storage medium 92 can store an operating system and a computer program. The computer program includes program instructions. When the program instructions are executed, the processor can be made to execute any one of the feature vector extraction methods provided by the embodiments of this application.

[0233] The processor is used to provide computing and control capabilities to support the operation of the entire computer device.

[0234] The internal memory 94 provides an environment for the operation of a computer program in a non-volatile storage medium. When the computer program is executed by the processor, the processor can be caused to execute any one of the feature vector extraction methods provided by the embodiments of the present application.

[0235] The network interface 95 is used for network communication, such as sending assigned tasks and the like. Those skilled in the art can understand that Figure 1A the structure shown in is only a block diagram of some structures related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine some components, or have different component arrangements.

[0236] It should be understood that the processor 91 may be a central processing unit (CPU), and the processor may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field programmable gate arrays (FPGAs) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. Among them, the general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc.

[0237] The electronic device 90 of the embodiments of the present application may be included in terminal devices such as tablet computers, laptop computers, mobile phones, supercomputers, smart wearable devices, etc., and may also be applied to databases, servers, and service response systems based on terminal artificial intelligence. The embodiments of the present application do not impose any restrictions on the specific type of the electronic device.

[0238] For example, the electronic device can be a station (STAION, ST) in a WLAN, a cellular phone, a cordless phone, a Session Initiation Protocol (SIP) phone, a Wireless Local Loop (WLL) station, a handheld device with wireless communication capabilities, a computing device, or other processing devices connected to a wireless modem, a computer, a laptop, a handheld communication device, a handheld computing device, and / or other devices for communicating on a wireless system, as well as next-generation communication systems, such as mobile terminals in a 5G network, mobile terminals in a future evolved Public Land Mobile Network (PLMN), or mobile terminals in a future evolved Non-terrestrial Network (NTN).

[0239] By way of example and not limitation, when the electronic device is a wearable device, the wearable device can also be a general term for devices that are intelligently designed for daily wear using wearable technology and developed into wearable devices, such as gloves, watches, etc. configured with a near-field communication module. A wearable device is a portable device that is either directly worn on the body or integrated into the user's clothes or accessories. By attaching to the user's body and using a pre-bound electronic card, it can perform operations such as payment and authentication. A wearable device is not just a hardware device, but more importantly, it realizes powerful functions through software support, data interaction, and cloud interaction. Broadly speaking, wearable intelligent devices include those with full functions and large sizes that can achieve complete or partial functions without relying on a smartphone, such as smart watches or smart glasses, as well as those that only focus on a certain type of application function and need to cooperate with other devices such as smartphones, such as various smart watches and smart bracelets with a display screen.

[0240] The descriptions of the processes or structures corresponding to the above respective figures each have their own focuses. For parts not detailed in a certain process or structure, reference can be made to the relevant descriptions of other processes or structures.

[0241] The above embodiments merely illustrate the principles and effects of the present application and are not intended to limit the present application. Any person familiar with this technology can modify or change the above embodiments without departing from the spirit and scope of the present application. Therefore, all equivalent modifications or changes made by those with ordinary knowledge in the technical field without departing from the spirit and technical idea disclosed in the present application should still be covered by the claims of the present application.

Claims

1. A method for analyzing the trajectory of a magnetocardiogram multipole, characterized in that, The method includes: Obtaining a multipole trajectory map based on an original magnetocardiogram atlas; Determining multiple trajectory branches in the multipole trajectory map, and dividing the multiple trajectory branches into a positive main branch, a negative main branch, and a secondary branch; Determining a multipole trajectory score corresponding to the multipole trajectory map based on the positive main branch, the negative main branch, and the secondary branch; Determining a trajectory deviation degree between the original magnetocardiogram atlas and a preset magnetocardiogram atlas based on the multipole trajectory score, and obtaining a trajectory analysis result of the original magnetocardiogram atlas based on the trajectory deviation degree.

2. The method for analyzing the multi-pole trajectory of magnetocardiogram according to claim 1, wherein: The determining the multipole trajectory score corresponding to the multipole trajectory map based on the positive main branch, the negative main branch, and the secondary branch includes: Determining a first trajectory parameter of the positive main branch based on the main branch information corresponding to the positive main branch, where the main branch information corresponding to the positive main branch includes the magnetic field value of the positive pole on the positive main branch and the position information of the positive pole; Determining a second trajectory parameter of the negative main branch based on the main branch information corresponding to the negative main branch; Matching the positive poles on the positive main branch and the negative poles on the negative main branch within the same frame to obtain a main branch dipole, and determining a third trajectory parameter of the main branch dipole; Determining a fourth trajectory parameter of the secondary branch based on the secondary branch information corresponding to the secondary branch; Determining the multipole trajectory score based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter.

3. The method for analyzing the multi-pole trajectory of magnetocardiogram according to claim 1, characterized in that: The first trajectory parameter includes an overall trajectory parameter and a segmented trajectory parameter. The determining the first trajectory parameter of the positive main branch based on the main branch information corresponding to the positive main branch includes: Obtaining the magnetic field value of the positive poles on the positive main branch that are greater than a preset pole threshold and the position information of the positive poles; Based on the magnetic field value and the position information, determining the overall trajectory parameter of the positive main branch, where the overall trajectory parameter includes an effective length, a main pole length, and a non-boundary rate; Sequentially determining the planar distance between adjacent positive poles on the positive main branch, and determining a main segment of the positive main branch where the planar distance is less than a preset offset distance; Determining a head segment of the positive main branch located in front of the main segment; Determining a tail segment of the positive main branch located behind the main segment; Respectively determining the segmented trajectory parameters of the main segment, the head segment, and the tail segment, where the segmented trajectory parameters include a trajectory path, a trajectory span, and a trajectory elastic curvature.

4. The method for analyzing the multi-pole trajectory of a magnetocardiogram according to claim 3, characterized in that: The determining the overall trajectory parameter of the positive main branch based on the magnetic field value and the position information includes: Based on the magnetic field values of all the positive poles greater than the preset pole threshold, determining the effective length corresponding to the positive main branch; Based on the proportion of the magnetic field value equal to the maximum magnetic field value, obtaining the main pole length corresponding to the positive main branch; Based on the proportion of the planar coordinates of the positive poles on the positive main branch that are not within a preset range of the preset matrix boundary, obtaining the non-boundary rate corresponding to the positive main branch.

5. The method for analyzing the multi-pole trajectory of magnetocardiogram according to claim 1, wherein: The determining the third trajectory parameter of the main branch dipole includes: The intersection of the positive pole on the positive main branch and the negative pole on the negative main branch is determined as the main branch of the main branch dipole; sequentially determining the plane distances of the positive poles of the adjacent frames on the positive main branch, and determining a section on the positive main branch where the plane distance is less than a preset offset distance as a main section of the positive main branch; sequentially determining the plane distances of the negative poles of the adjacent frames on the negative main branch, and determining a section on the negative main branch where the plane distance is less than a preset offset distance as a main section of the negative main branch; Determine the intersection of the main segment of the positive main branch and the main segment of the negative main branch as the main segment of the main branch dipole; Taking a difference between a main branch of the main branch dipole and a main segment of the main branch dipole to obtain a secondary segment of the main branch dipole; Based on the main branch information, third trajectory parameters of the main branch, the main segment and the secondary segment of the main branch dipole are obtained in sequence.

6. The method for analyzing the multi - dipole trajectory of magnetocardiogram according to claim 1, wherein: The fourth trajectory parameters include an effective number of branches, a branch size, a total trajectory path, a total trajectory span, a total trajectory elastic curvature, and a total trajectory elastic strength. The fourth trajectory parameters of the secondary branch are determined based on the secondary branch information corresponding to the secondary branch, including: Determine the non-boundary rate of each of the sub-branches, remove the sub-branches whose non-boundary rate is less than or equal to a preset boundary threshold, and obtain the boundary factor of the sub-branch; Obtain the sum of the magnetic field values ​​between the start time and the end time of all sub-branches; Dividing the magnetic field value of each of the secondary branches by the total magnetic field values ​​to obtain a length factor of each of the secondary branches; Sum all the boundary factors to obtain the effective number of branches; Multiplying each of the boundary factors and the length factor of the secondary branch corresponding to the boundary factor in turn and then summing the products to obtain the branch size of the secondary branch; Determining the trajectory path, trajectory span, and trajectory elastic curvature of each of the sub-branches, and obtaining the trajectory elastic strength based on the trajectory elastic curvature; The trajectory paths, the trajectory spans, the trajectory elastic curvatures, and the trajectory elastic strengths of all the sub-branches are summed up respectively to obtain a total trajectory path, a total trajectory span, a total trajectory elastic curvature, and a total trajectory elastic strength.

7. The magnetocardiogram multipole trajectory analysis method according to claim 1, wherein: The determining of the multipole trajectory score based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter and the fourth trajectory parameter comprises: Determine a positive pole trajectory score, a negative pole trajectory score, a dipole trajectory score, and a secondary branch trajectory score based on the first trajectory parameter, the second trajectory parameter, the third trajectory parameter, and the fourth trajectory parameter, respectively; The multipole trajectory score is obtained based on the positive pole trajectory score, the negative pole trajectory score, the dipole trajectory score and the secondary branch trajectory score.

8. A magnetocardiogram multipole trajectory analysis device, characterized in that, The device comprises: A multipole trajectory map acquisition module, used for acquiring a multipole trajectory map based on an original magnetocardiographic atlas; A trajectory branch determination module, used to determine a plurality of trajectory branches in the multipole trajectory diagram, and divide the plurality of trajectory branches into a positive main branch, a negative main branch and a secondary branch; A multipole trajectory scoring determination module, configured to determine a multipole trajectory score corresponding to the multipole trajectory map based on the positive main branch, the negative main branch, and the secondary branch; A trajectory analysis result determination module, configured to determine a trajectory deviation degree between the original magnetocardiogram atlas and a preset magnetocardiogram atlas based on the multipole trajectory score, and obtain a trajectory analysis result of the original magnetocardiogram atlas based on the trajectory deviation degree.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the magnetocardiogram multipole trajectory analysis method according to any one of claims 1 to 7.

10. An electronic device, characterized in that, The electronic device includes: A memory storing a computer program; A processor communicatively connected to the memory, and when calling the computer program, executes the magnetocardiogram multipole trajectory analysis method according to any one of claims 1 to 7.

Images