Cardiac magnetic dipole center track parameter determination system

Abstract

The invention provides a magnetocardiogram dipole center track parameter determination system. The system comprises an acquisition device, a generation device, a determination device and a calculation device. Wherein the acquisition device is used for acquiring a magnetocardiogram data set of a target wave band in a target magnetocardiogram data set; the generating device is used for generating a target isomagnetic map of the target wave band based on the magnetocardiogram data set of the target wave band; the determining device is used for determining a target dipole center pixel point set of a target wave band based on the dipole position of each frame of isomagnetic map in the target isomagnetic map; and the calculation device is used for calculating the center track parameters of the target dipole based on the center pixel point set of the target dipole and a pre-created standard center track distribution map. Therefore, the change rule of the center of the magnetocardiogram dipole can be accurately and intuitively described, and the distribution states of the dipoles in different wave bands can be accurately quantified through the calculated center track parameters of the dipoles.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of magnetocardiogram, in particular to a magnetocardiogram dipole center trajectory parameter determination system. BACKGROUND

[0002] Magnetocardiogram is a non-invasive, non-radiation and non-contact detection technology of cardiac electromagnetic function. It can record the weak magnetic field generated by the cardiac electrical activity by using a super-high sensitivity magnetic field detection array. The biggest feature of magnetocardiogram is the dipole formed by the positive and negative magnetic poles. The relative position of the dipole changes constantly in the whole cycle of cardiac electrical activity, and presents different rules in different wave bands, such as the rotation of the dipole in the QRS wave band and the stable state in the T wave band. However, it is difficult to directly analyze the state of the dipole, and it is also difficult to directly quantify the change rule and phenomenon of the dipole center. SUMMARY

[0003] Therefore, the present application aims to provide a magnetocardiogram dipole center trajectory parameter determination system and method, which can accurately and intuitively describe the change rule of the magnetocardiogram dipole center, and can accurately quantify the dipole distribution state in different wave bands through the dipole center trajectory parameter.

[0004] The present application provides a magnetocardiogram dipole center trajectory parameter determination system, which comprises an acquisition device, a generation device, a determination device and a calculation device.

[0005] The acquisition device is configured to acquire a magnetocardiogram data set of a target wave band. The target wave band comprises at least one of a QRS wave band, a T wave band and a P wave band.

[0006] The generation device is configured to generate a target isomagnetic map of the target wave band based on the magnetocardiogram data set of the target wave band.

[0007] The determination device is configured to determine a target dipole center pixel point set of the target wave band based on the dipole position of each frame of the target isomagnetic map.

[0008] The calculation device is configured to calculate a target dipole center trajectory parameter based on the target dipole center pixel point set and a pre-created standard center trajectory distribution atlas. The standard center trajectory distribution atlas is generated based on a plurality of sample magnetocardiogram data sets. The target dipole center trajectory parameter is used to represent the distribution of the target dipole center pixel point set in each pixel position level of the standard center trajectory distribution atlas.

[0009] Further, the determining apparatus is specifically configured to:

[0010] For any frame of isodynamic map, based on the magnetic field strength of the frame of isodynamic map, determine the position coordinates of the positive pole and the position coordinates of the negative pole in the frame of isodynamic map;

[0011] According to the position coordinates of the positive pole and the position coordinates of the negative pole, determine the center point position coordinates of the dipole;

[0012] Based on the center point position coordinates of the dipole determined by each frame of isodynamic map, determine the target dipole center pixel point set.

[0013] Further, the system further comprises a creating apparatus, and the creating apparatus comprises a first obtaining module, a first determining module, a second determining module and a third determining module.

[0014] The first obtaining module is configured to obtain a magnetocardiogram data set of a target wave band in a plurality of sample magnetocardiogram data sets.

[0015] The first determining module is configured to, for any sample magnetocardiogram data set, determine a sample dipole center pixel point set based on the magnetocardiogram data set of the target wave band in the sample magnetocardiogram data set.

[0016] The second determining module is configured to determine a density distribution map composed of each sample dipole center pixel point based on a plurality of sample dipole center pixel point sets.

[0017] The third determining module is configured to determine the standard center trajectory distribution map based on the density distribution map.

[0018] Further, the third determining module is specifically configured to:

[0019] Determine the abnormal probability of each pixel position based on the total number of sample dipole center pixel points included in each pixel position in the density distribution map and the label information of each sample dipole center pixel point.

[0020] Generate a probability distribution map based on the abnormal probability of each pixel position in the density distribution map.

[0021] Determine the position level of each pixel position based on the probability distribution map.

[0022] Determine the standard center trajectory distribution map based on the position level of each pixel position.

[0023] Further, the computing device comprises a fourth determining module and a computing module;

[0024] The fourth determining module is configured to determine the number of dipole center pixel points of each pixel position level in the standard center trajectory distribution map.

[0025] The computing module is configured to calculate the target dipole center trajectory parameter based on the number of dipole center pixel points included in each pixel position level and a preset dipole center trajectory parameter formula.

[0026] Further, the preset dipole center trajectory parameter formula is:

[0027]

[0028] Wherein, dctj is the dipole center trajectory parameter, nc1 is the number of dipole center pixel points in the pixel position level 1, nc2 is the number of dipole center pixel points in the pixel position level 2, …, nc n is the number of dipole center pixel points in the pixel position level n, a1, a2, …, a n are the weighting coefficients of nc1, nc2, …, nc n respectively.

[0029] Further, the acquisition device comprises a second acquisition module and a fifth determining module;

[0030] The second acquisition module is configured to acquire the target magnetocardiogram data set collected by the multi-channel magnetocardiograph.

[0031] The fifth determining module is configured to determine the magnetocardiogram data set of the target wave band based on the target magnetocardiogram data set.

[0032] Further, the number of channels of the multi-channel magnetocardiograph is one of 4 channels, 9 channels, 36 channels and 64 channels.

[0033] Further, when the generation device is used to generate the target isomagnetogram of the target wave band based on the magnetocardiogram data set of the target wave band, the generation device is specifically configured to:

[0034] Generate the target isomagnetogram of the target wave band based on the magnetocardiogram data set of the target wave band and the position of each channel in the multi-channel magnetocardiograph.

[0035] Further, the third determining module is specifically configured to determine the anomaly probability of each pixel position according to the following steps:

[0036] For any pixel position, the total number of sample dipole center pixel points included in the pixel position and the target number of sample dipole center pixel points of an abnormal type are acquired;

[0037] A ratio between the target number and the total number is determined as an abnormal probability of the pixel position.

[0038] An embodiment of the present application provides a heart magnetic dipole center trajectory parameter determination system, which determines a target dipole center pixel point set of the target wave band through a magnetogram data set of the target wave band, and calculates a target dipole center trajectory parameter based on the target dipole center pixel point set and a pre-created standard center trajectory distribution atlas. Compared with related technologies which cannot directly quantify the dipole distribution state of different wave bands, the present application can accurately quantify the dipole distribution state of different wave bands through the dipole center trajectory parameter.

[0039] In order to make the above objectives, characteristics and advantages of the present application more apparent, clear and easy to understand, the following preferred embodiments are specifically described below, and the accompanying drawings are referred to for detailed description. BRIEF DESCRIPTION OF DRAWINGS

[0040] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings needed to be used in the embodiments. It should be understood that the following drawings only show some embodiments of the present application, and therefore should not be regarded as a limitation to the scope. For those skilled in the art, other related drawings can also be obtained without creative labor on the basis of these drawings.

[0041] Figure 1 Fig. 1 shows a structure schematic diagram of a heart magnetic dipole center trajectory parameter determination system provided by an embodiment of the present application;

[0042] Figure 2 Fig. 2 shows a structure schematic diagram of a heart magnetic dipole center trajectory parameter determination system provided by an embodiment of the present application;

[0043] Figure 3 Fig. 3 shows a structure schematic diagram of a heart magnetic dipole center trajectory parameter determination system provided by an embodiment of the present application;

[0044] Figure 4 Fig. 4 shows a structure schematic diagram of a heart magnetic dipole center trajectory parameter determination system provided by an embodiment of the present application;

[0045] Figure 5 Fig. 5 shows a density distribution diagram provided by an embodiment of the present application;

[0046] Figure 6A standard central track distribution map provided by the embodiment of the present application is shown. DETAILED DESCRIPTION

[0047] To make the objectives, technical solutions, and superiorities of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described below in connection with the drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. The components of the embodiments of the present application described and shown in the drawings herein can be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of the present application provided in the drawings is not intended to limit the scope of the claimed present application, but only represents selected embodiments of the present application. Based on the embodiments of the present application, every other embodiment obtained by a person skilled in the art without creative work falls within the scope of the present application.

[0048] First, the application scenarios applicable to the present application are introduced. The present application can be applied to the wave band data processing scenario in magnetocardiogram.

[0049] The biggest feature of magnetocardiogram is the dipole formed by the positive and negative magnetic poles. The relative position of the dipole changes constantly in the whole cycle of cardiac electrical activity, and presents different rules in different wave bands, such as the dipole rotation in QRS wave band and the stable state in T wave band. However, it is very difficult to analyze the state of the dipole directly, and it is also very difficult to quantify the change rule and phenomenon of the dipole center directly.

[0050] Please refer to Figures 1-4 , Figure 1 FIG. 1 is a structural schematic diagram of a magnetocardiogram dipole center track parameter determination system provided by the embodiment of the present application. Figure 2 FIG. 2 is another structural schematic diagram of a magnetocardiogram dipole center track parameter determination system provided by the embodiment of the present application. Figure 3 FIG. 3 is a third structural schematic diagram of a magnetocardiogram dipole center track parameter determination system provided by the embodiment of the present application. Figure 4 FIG. 4 is a fourth structural schematic diagram of a magnetocardiogram dipole center track parameter determination system provided by the embodiment of the present application. As shown in FIG. 4, the magnetocardiogram dipole center track parameter determination system 100 provided by the embodiment of the present application comprises an acquisition device 110, a generation device 120, a determination device 130, and a calculation device 140. Figure 1

[0051] The acquisition device 110 is configured to acquire a magnetocardiogram data set of a target wave band. The target wave band comprises at least one of a QRS wave band, a T wave band, and a P wave band.

[0052] ​Here, the magnetocardiogram dataset of the target wave band in the entire cycle target magnetocardiogram dataset is obtained, and the target wave band includes at least one of a QRS wave band, a T wave band, a P wave band, and the like. The QRS wave band represents electrical activity generated when the ventricle depolarizes and activates, the T wave band represents the process of cardiac repolarization, that is, the process of the ventricle from depolarization and activation to the state of returning to rest, which is usually the process of ventricular diastole, and the P wave band represents atrial electrical activity.

[0053] Here, the magnetocardiogram dataset of the target wave band can be a magnetocardiogram dataset obtained after filtering and noise reduction processing, which can ensure that the processed dataset is a high-quality dataset with almost no noise influence, ensure the reliability of the obtained magnetocardiogram dataset of the target wave band, and further improve the accuracy of the target dipole center trajectory parameters finally calculated.

[0054] The generation apparatus 120 is configured to generate a target isomagnetogram of the target wave band based on the magnetocardiogram dataset of the target wave band.

[0055] The determination apparatus 130 is configured to determine a target dipole center pixel point set of the target wave band based on the dipole positions of each frame of isomagnetogram in the target isomagnetogram.

[0056] In specific implementation, the target isomagnetogram of the target wave band can be generated based on the magnetocardiogram dataset of the target wave band after filtering and noise reduction processing. The target isomagnetogram is composed of multiple frames of isomagnetogram, and each frame of isomagnetogram has a corresponding dipole position. Based on the dipole positions of each frame of isomagnetogram in the target isomagnetogram, that is, based on the dipole positions corresponding to all frames of isomagnetogram in the target isomagnetogram, a target dipole center pixel point set under the target wave band can be determined. Here, each dipole position includes a positive dipole position and a negative dipole position.

[0057] The calculation apparatus 140 is configured to calculate target dipole center trajectory parameters based on the target dipole center pixel point set and a pre-created standard center trajectory distribution atlas. The standard center trajectory distribution atlas is generated based on multiple sample magnetocardiogram datasets. The target dipole center trajectory parameters are used to represent the distribution of the dipole center pixel point set in each pixel position level in the standard center trajectory distribution atlas.

[0058] In specific implementation, the target dipole center trajectory parameters can be calculated according to the obtained target dipole center pixel point set and the pre-created standard center trajectory distribution atlas.

[0059] Here, the pixel position level is the pixel position to which a level is assigned. It can be understood that in the standard central trajectory distribution atlas, a plurality of pixel positions can be divided, and each pixel position corresponds to a level. The levels corresponding to each pixel position can have the same number of levels. Here, by way of example, in the standard central trajectory distribution atlas, there are 100 pixel positions, of which 20 pixel positions are level 1, 40 pixel positions are level 2, 25 pixel positions are level 3, and 15 pixel positions are level 4.

[0060] Here, it should be noted that the target dipole central trajectory parameter calculated can quantify the distribution state of the dipole, i.e., the numerical value of the dipole central trajectory parameter can intuitively represent the distribution of the dipole central pixel point set in each pixel position level in the standard central trajectory distribution atlas. That is, the numerical value of the dipole central trajectory parameter can determine whether the number of dipole central pixels distributed in each pixel position level is relatively large or small. Further, it can also represent whether the number of the overall dipole distribution in each pixel position level is relatively large or small. By way of example, the larger the numerical value of the dipole central trajectory parameter, the more dipole central pixels in the pixel position level in the standard central trajectory distribution atlas.

[0061] The present application considers that the related art cannot directly quantify the dipole state of different wave bands. Therefore, the present application determines a target dipole central pixel point set of a target wave band based on a magnetogram data set of the target wave band, calculates a target dipole central trajectory parameter based on the target dipole central pixel point set and a pre-created standard central trajectory distribution atlas. In this way, the variation of the magnetocardiographic dipole center can be accurately and intuitively described, and the dipole distribution state of different wave bands can be accurately quantified by the calculated dipole central trajectory parameter.

[0062] Further, as shown in Figure 1 When determining the target dipole central pixel point set of the target wave band based on the positions of the dipoles in each frame of magnetogram in the target magnetogram, the determining apparatus 130 is specifically configured to:

[0063] For any frame of magnetogram, based on the magnetic field intensity of the frame of magnetogram, determine the position coordinates of the positive dipole and the position coordinates of the negative dipole in the frame of magnetogram;

[0064] According to the position coordinates of the positive dipole and the position coordinates of the negative dipole, determine the center point position coordinates of the dipole;

[0065] Based on the center point position coordinates of the dipole determined from each frame of magnetogram, determine the target dipole central pixel point set.

[0066] Here, for any frame of isogam, there are multiple sizes of magnetic field strength in any frame of isogam, the maximum magnetic field strength position coordinates are taken as the positive pole position coordinates of the dipole, the minimum magnetic field strength position coordinates are taken as the negative pole position coordinates of the dipole, and the center point position coordinates between the positive pole position coordinates and the negative pole position coordinates can be determined according to the positive pole position coordinates and the negative pole position coordinates, which are the center point position coordinates of the dipole. After determining the center point position coordinates of the dipole of each frame of isogam, the center point position coordinates of the dipole of each frame of isogam are combined, and a target dipole center pixel point set can be determined. Here, the center point position coordinates of a dipole can be composed of the center pixel point position coordinates of 1-4 dipoles.

[0067] Here, it should be noted that the center point position coordinates of the dipole are the position coordinates of the middle point of the straight line connecting the positive pole position coordinates of the dipole and the negative pole position coordinates of the dipole. The position coordinates are composed of horizontal coordinates and vertical coordinates, here, the horizontal coordinates of the center point position coordinates of a dipole are obtained by taking the arithmetic mean of the horizontal coordinates of the positive pole position coordinates and the horizontal coordinates of the negative pole position coordinates, and similarly, the vertical coordinates of the center point position coordinates of a dipole are obtained by taking the arithmetic mean of the vertical coordinates of the positive pole position coordinates and the vertical coordinates of the negative pole position coordinates, so that the center point position coordinates of the dipole can be obtained. The center point position coordinates of a dipole can represent the position of the center point of the dipole.

[0068] Here, exemplary, for each frame magnetogram of the QRS wave band (for example, a 100*100 tensor), assuming that there are three frames of magnetograms in the target magnetogram, the three frames of magnetograms determine the center point position coordinates of three dipoles as P0=(56, 72), P1=(48, 65), and P2=(67, 77), so that a target dipole center pixel point set corresponding to the QRS wave band is B={P0, P1, P2}={(56, 72), (48, 65), (67, 77)}. Of course, if the center point position coordinates of the dipole determined by the frame magnetogram are not integers, that is, in this special case, the position coordinates contain decimals, such as P0=(56.5, 72.5), at this time, the integer part of the coordinates can be obtained, that is, the center pixel point position coordinates of all dipoles covering the center point position coordinates of the dipole, that is, the center point position coordinates of the dipole P0={(56, 72), (56, 73), (57, 72), (57, 73)}. In this way, a target dipole center pixel point set B={P0, P1, P2}={(56, 72), (56, 73), (57, 72), (57, 73), (48, 65), (67, 77)}. Here, a target dipole center pixel point set is composed of multiple dipole center point position coordinates, which can also be understood as a target dipole center pixel point set composed of multiple target dipole center pixel point position coordinates representing the center point of the dipole. At the same time, if the center point position coordinates of a dipole are integers, the center pixel point position coordinates of the dipole are the center point position coordinates of the dipole, and if the center point position coordinates of a dipole contain decimals, the center pixel point position coordinates of two or four dipoles represent the center point position coordinates of the dipole. Of course, any one of the two or four center pixel point position coordinates can describe or represent a dipole. Here, it should be noted that the manner of determining a target dipole center pixel point set based on the center point position coordinates of each dipole determined by each frame magnetogram includes but is not limited to the above manner.

[0069] Further, as shown in Figure 2 The system further comprises a creating device 150, and the creating device 150 comprises a first acquiring module 151, a first determining module 152, a second determining module 153, and a third determining module 154.

[0070] The first acquiring module 151 is configured to acquire magnetogram data sets of a target wave band in a plurality of sample magnetocardiogram data sets.

[0071] The first determining module 152 is configured to determine a sample dipole center pixel point set based on the magnetogram data sets of the target wave band in the sample magnetocardiogram data set for any sample magnetocardiogram data set.

[0072] The second determining module 153 is configured to determine and generate a density distribution map composed of the central pixel points of each sample dipole based on the plurality of sample dipole central pixel point sets;

[0073] The third determination module 154 is configured to determine the standard center trajectory distribution map based on the density distribution map.

[0074] Here, each sample corresponds to a sample magnetocardiogram dataset, and each sample magnetocardiogram dataset corresponds to a magnetocardiogram dataset with a target band. A magnetocardiogram dataset of the target band is obtained from multiple sample magnetocardiogram datasets. For any sample magnetocardiogram dataset, a sample isomagnetic map of the target band can be obtained based on the magnetocardiogram dataset of the target band in the sample magnetocardiogram dataset. Then, based on the dipole positions of each frame of the isomagnetic map in the sample isomagnetic map, the center point positions of the sample dipoles corresponding to each frame of the isomagnetic map of the target band are determined. The coordinates of the center pixels of duplicate dipoles are then removed to obtain a sample dipole center pixel set for the target band of the sample. Based on multiple dipole center pixel sets of multiple samples, according to the center pixel position coordinates of the sample dipoles, all the sample dipole center pixel points in the dipole center pixel set corresponding to the multiple samples are added to the corresponding pixel positions in the tensor of preset size according to the center pixel position coordinates of each sample dipole. Then, a density distribution map consisting of the center pixel points of the dipoles of each sample can be generated. It can also be understood as generating a density distribution map consisting of the center points of the dipoles of each sample. Then, based on the density distribution map, the standard center trajectory distribution map is determined. Here, for example, the preset size tensor can be a 100*100 tensor.

[0075] Among them, Figure 5 The density distribution map shown is formed by accumulating the center points of all sample dipoles in the multiple dipole center pixel points corresponding to the multiple samples under the target band according to the center pixel position coordinates of each sample dipole to the corresponding pixel position in the tensor of the preset size. For example, if there are 1000 samples, there are 1000 sample dipole center pixel points. Then, all the sample dipole centers in the 1000 sample dipole center pixel points under the QRS band are accumulated to a 100*100 tensor to form a density distribution map.

[0076] It should be noted that, for a sample dipole center pixel set in a target wave band of a sample, the center pixel position coordinates of the dipoles in the sample dipole center pixel set are usually not repeated; for multiple sample dipole center pixel sets, the center pixel position coordinates of the dipoles in different sample dipole center pixel sets can be repeated.

[0077] Further, as shown in Figure 2 The third determining module 154 is specifically configured to:

[0078] determine the abnormal probability of each pixel position based on the total number of sample dipole center pixel points included in each pixel position in the density distribution map and the label information of each sample dipole center pixel point;

[0079] generate a probability distribution map based on the abnormal probability of each pixel position in the density distribution map;

[0080] determine the position level of each pixel position based on the probability distribution map;

[0081] determine the standard center trajectory distribution atlas based on the position level of each pixel position.

[0082] Here, for any pixel position in the density distribution map, the total number of sample dipole center pixel points included in the pixel position in the density distribution map is counted, and the label information carried by each sample dipole center pixel point included in the pixel position is viewed, and then the abnormal probability of each pixel position is calculated. After determining the abnormal probability of each pixel position, the probability distribution map can be generated. Then, according to the pre-specified abnormal probability and level corresponding rule, the abnormal probability of each pixel position in the probability distribution map is assigned to the position level of each pixel position. After determining the corresponding position level of each pixel position in the probability distribution map, the standard center trajectory distribution atlas is obtained.

[0083] Here, the label information of each sample is known, so the label information of the corresponding sample dipole center pixel set in each sample is known, and therefore the label information of all sample dipole center pixel points in the sample dipole center pixel set is known. The label information can represent an abnormal type or a normal type. For example, the label information can be represented by 0 or 1, for example, 0 represents an abnormal type and 1 represents a normal type. If the label of a sample is 0, the label information of all sample dipole center pixel points in the sample dipole center pixel set is 0.

[0084] The pixel position on the density distribution map is generally of four types, i.e., a pixel position only with abnormal objects, a pixel position only with normal objects, a pixel position containing both normal objects and abnormal objects, and a pixel position without normal and abnormal objects (blank pixel position). Here, the pixel position only with abnormal objects is that the label information of each sample dipole center pixel point in the pixel position is all of abnormal type; the pixel position only with normal objects is that the label information of each sample dipole center pixel point in the pixel position is all of normal type; the pixel position containing both normal and abnormal objects is that the label information of each sample dipole center pixel point in the pixel position is both of normal type and abnormal type; and the pixel position without normal and abnormal objects (blank pixel position) is that the pixel position has no sample dipole center pixel point.

[0085] Here, it should be noted that since the blank pixel position has no sample dipole center pixel point, the abnormal probability value cannot be calculated, and therefore the blank pixel position can be directly recorded as an abnormal probability of 0.

[0086] Here, by way of example, the pre-specified abnormal probability and grade corresponding rule can be that an abnormal probability in the range of [0, 0.2] corresponds to a grade of 1, an abnormal probability in the range of (0.2, 0.4] corresponds to a grade of 2, an abnormal probability in the range of (0.4, 0.6] corresponds to a grade of 3, an abnormal probability in the range of (0.6, 0.8] corresponds to a grade of 4, and an abnormal probability in the range of (0.8, 1] corresponds to a grade of 5. Then, according to the above corresponding rule, the abnormal probability determined in each pixel position in the probability distribution map can be assigned a position grade of each pixel position, and different grades of pixel positions can be represented by different colors, for example, the colors corresponding to grades from 1 to 5 are from light to dark, so that a standard center trajectory distribution map can be obtained. Further, from the standard center trajectory distribution map, it can be seen that the deeper the color of a pixel position, the greater the grade and the greater the abnormal probability. Here, the pre-specified abnormal probability and grade corresponding rule can also be extended, i.e., the grade is divided into n grades, and the abnormal probability corresponding to the n grades covers the probability range of 0-1, and the abnormal probability from low to high corresponds to grades 1, 2, 3, …, n in turn.

[0087] Here, as Figure 6The standard center trajectory distribution atlas shown can represent the distribution of the dipole center pixel points, that is, can represent the overall distribution of the dipoles. For example, the pixel position with a darker color (the greater the anomaly probability) represents a pixel position where the abnormal objects are more concentrated, that is, represents that the dipole center points representing the abnormal objects are more concentrated at the pixel position, and can also be understood as that the dipoles of the abnormal objects are more concentrated at the pixel position. The color distribution of the standard center trajectory distribution atlas has a certain rule, and the larger the number of samples, the more obvious the rule of the standard center trajectory distribution atlas (that is, the fixed color position is darker). Here, an object can be a sample, and for an abnormal object, the label information of the sample can be represented as an abnormal type.

[0088] Further, as shown in Figure 3 The computing device 140 includes a fourth determination module 141 and a calculation module 142.

[0089] The fourth determination module 141 is configured to determine the number of dipole center pixel points of the target dipole center pixel points that are concentrated in the pixel position level of the standard center trajectory distribution atlas.

[0090] The calculation module 142 is configured to calculate the target dipole center trajectory parameter based on the number of dipole center points included in each pixel position level and a preset dipole center trajectory parameter formula.

[0091] The number of target dipole center pixel points of the target dipole center pixel points that are concentrated in the pixel position level of the standard center trajectory distribution atlas is determined. According to the number of dipole center points included in each level and the preset dipole center trajectory parameter formula, the target dipole center trajectory parameter can be calculated.

[0092] Further, the preset dipole center trajectory parameter formula is:

[0093]

[0094] wherein dctj is the dipole center trajectory parameter, nc1 is the number of dipole center pixel points in the pixel position level 1, nc2 is the number of dipole center pixel points in the pixel position level 2, …, nc n is the number of dipole center pixel points in the pixel position level n, a1, a2, …, a n are the weighting coefficients of nc1, nc2, …, nc n , respectively.

[0095] The preset dipole center trajectory parameter formula can be determined based on the number of levels and the level weighting. Here, a1+a n= 0, a2+a n-1 = 0, a3+a n-2 = 0, that is, all the pairs of weighting coefficients whose indices add up to n+1 add up to 0, and if n is odd, then Here, for example, if n is 5, then

[0096] Here, for the convenience of understanding the preset dipole center trajectory parameter formula, for example, assuming that the number of levels is 5, then based on the number of levels, the preset dipole center trajectory parameter formula can be determined to be specifically:

[0097]

[0098] wherein dctj is a dipole center trajectory parameter, nc1 is the number of dipole center pixel points in the pixel position level 1, nc2 is the number of dipole center pixel points in the pixel position level 2, nc3 is the number of dipole center pixel points in the pixel position level 3, nc4 is the number of dipole center pixel points in the pixel position level 4, and nc5 is the number of dipole center pixel points in the pixel position level 5.

[0099] Here, since nc3 in the preset dipole center trajectory parameter formula is the number of dipole center pixel points in the pixel position level 3, and the abnormal probability corresponding to the pixel position level 3 is usually in the range of about 0.5, for example, (0.4, 0.6], which belongs to the “middle ground”, that is, the abnormal probability is neither too large nor too small, so no weighting is performed, and therefore the weighting parameter in front of nc3 is 0. When (nc1, nc2, nc3, nc4, nc5) = (2, 5, 0, 7, 5), it means that the number of dipole center pixel points in the pixel position level 1 is 2, the number of dipole center pixel points in the pixel position level 2 is 5, the number of dipole center pixel points in the pixel position level 3 is 0, the number of dipole center pixel points in the pixel position level 4 is 7, and the number of dipole center pixel points in the pixel position level 5 is 5. Here, it can be easily known that the sum of the number of dipole center pixel points in the pixel position level 4 and the number of dipole center pixel points in the pixel position level 5 is 5 more than the sum of the number of dipole center pixel points in the pixel position level 1 and the number of dipole center pixel points in the pixel position level 2. Then, substituting (nc1, nc2, nc3, nc4, nc5) = (2, 5, 0, 7, 5) into the above preset dipole center trajectory parameter formula, we can get When (nc1, nc2, nc3, nc4, nc5) = (1, 2, 0, 9, 9), here, it can be easily known that the sum of the number of the dipole center pixels in the pixel position level of the 4th level and the number of the dipole center pixels in the pixel position level of the 5th level is 15 more than the sum of the number of the dipole center pixels in the pixel position level of the 1st level and the number of the dipole center pixels in the pixel position level of the 2nd level. It can be easily known that Thus, the value of the dipole center trajectory parameter obtained by the above calculation is 0.842 and the value of the dipole center trajectory parameter is 2.19, 2.19 is greater than 0.842, the greater the value of the dipole center trajectory parameter, the relatively greater the difference between the sum of the number of the dipole center pixels in the pixel position level of the 4th level and the number of the dipole center pixels in the pixel position level of the 5th level and the sum of the number of the dipole center pixels in the pixel position level of the 1st level and the number of the dipole center pixels in the pixel position level of the 2nd level. That is, the greater the value of the dipole center trajectory parameter, the relatively more the number of the dipole center pixels in the pixel position level of the 4th level and the 5th level compared to the number of the dipole center pixels in the pixel position level of the 1st level and the 2nd level which are concentrated in the standard center trajectory distribution atlas. This also indicates that the target dipole center pixels are more concentrated in the 4th level and the 5th level, or the overall can also reflect that the target dipole is more concentrated in the 4th level and the 5th level. Similarly, it can be easily verified that the smaller the value of the dipole center trajectory parameter, the relatively more the number of the dipole center pixels in the pixel position level of the 1st level and the 2nd level compared to the number of the dipole center pixels in the pixel position level of the 4th level and the 5th level which are concentrated in the standard center trajectory distribution atlas. This also indicates that the target dipole center pixels are more concentrated in the 1st level and the 2nd level, and further, the overall can also reflect that the target dipole is more concentrated in the 1st level and the 2nd level.

[0100] Here, for example, the value of the calculated dipole center trajectory parameter is the result of quantifying the dipole distribution state, that is, the value of the dipole center trajectory parameter can intuitively and accurately determine whether the target dipole center pixels are more concentrated in the 4th level and the 5th level or more concentrated in the 1st level and the 2nd level. It can also be that the number of dipoles in the 4th level and the 5th level is relatively more or less than the number of dipoles in the 1st level and the 2nd level.

[0101] Here, when the preset dipole center trajectory parameter formula is The dipole center trajectory parameter value range is between [-4, 4], which can be derived according to the dipole center trajectory parameter formula.

[0102] Further, such as Figure 4 As shown, the acquisition device 110 includes a second acquisition module 111 and a fifth determination module 112:

[0103] The second acquisition module 111 is configured to acquire a target magnetocardiogram dataset collected by a multi-channel magnetocardiograph;

[0104] The fifth determining module 112 is configured to determine a target band magnetocardiogram dataset based on the target magnetocardiogram dataset.

[0105] Here, a target magnetocardiogram dataset for the entire cycle acquired by a multi-channel magnetocardiograph is obtained. The target magnetocardiogram dataset is then segmented into individual bands to determine the target band's magnetocardiogram dataset. Segmentation of the individual bands can be performed to separate the QRS and T bands, as well as other desired bands. Each band has its own corresponding magnetocardiogram dataset.

[0106] Furthermore, the number of channels of the multi-channel magnetocardiograph is one of 4 channels, 9 channels, 36 channels, and 64 channels.

[0107] Here, multi-channel magnetocardiographs include, but are not limited to, superconducting magnetocardiographs and atomic magnetometer magnetocardiographs. Each type of multi-channel magnetocardiograph has a corresponding number of channels. The number of channels of a multi-channel magnetocardiograph can be any number of channels, such as 4 channels, 9 channels, 36 channels, or 64 channels.

[0108] Further, such as Figure 4 As shown, when the generating device 120 is used to generate the target isomagnetic map of the target band based on the magnetocardiogram dataset of the target band, the generating device 120 is specifically used to:

[0109] A target isomagnetic map of the target band is generated based on the magnetocardiogram dataset of the target band and the positions of each channel in the multi-channel magnetocardiograph.

[0110] Here, a target isomagnetic map of the target band is generated based on the magnetocardiogram dataset of the target band, the magnetic field intensity of each channel in the multi-channel magnetocardiograph, and the position of each channel in the multi-channel magnetocardiograph.

[0111] Further, such as Figure 4 As shown, the third determination module 154 is further configured to determine the abnormality probability of each pixel position according to the following steps:

[0112] For any pixel position, the total number of sample dipole center pixel points included in the pixel position and the target number of sample dipole center pixel points of an abnormal type are obtained;

[0113] The ratio between the target number and the total number is determined as the abnormal probability of the pixel position.

[0114] Here, in the density distribution map, for any pixel position, first, the total number of sample dipole center pixel points in the pixel position is obtained, and then the target number of sample dipole center pixel points of an abnormal type is obtained, and then the target number of sample dipole center pixel points of an abnormal type is obtained. The total number of sample dipole center pixel points is obtained, and the abnormal probability corresponding to the pixel position is obtained.

[0115] Here, for example, for any pixel position, the total number of sample dipole center pixel points in the pixel position is 10, and the target number of sample dipole center pixel points of an abnormal type is 8, and the abnormal probability corresponding to the pixel position is 0.8.

[0116] The application provides a magnetocardiogram dipole center trajectory parameter determination system, the system comprises an acquisition device, a generation device, a determination device, a calculation device; wherein the acquisition device is used for acquiring a target wave band magnetocardiogram data set; the target wave band comprises at least one of QRS wave band, T wave band and P wave band; the generation device is used for generating a target isomagnetogram of the target wave band based on the target wave band magnetocardiogram data set; the determination device is used for determining a target dipole center pixel point set of the target wave band based on the dipole position of each frame isomagnetogram in the target isomagnetogram; the calculation device is used for calculating a target dipole center trajectory parameter based on the target dipole center pixel point set and a pre-created standard center trajectory distribution map; wherein the standard center trajectory distribution map is generated by a plurality of sample magnetocardiogram data sets; the target dipole center trajectory parameter is used for representing the distribution of the dipole center pixel point set in each pixel position level in the standard center trajectory distribution map. In this way, the change rule of the magnetocardiogram dipole center can be accurately and intuitively described, and the dipole center trajectory parameter obtained by calculation can accurately quantify the dipole distribution state of different wave bands.

[0117] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of the system, device and unit described above can refer to the corresponding process in the foregoing method embodiments, which will not be described here.

[0118] In several embodiments provided in the present application, it should be understood that the disclosed system, device and method can be implemented in other manners. The described device embodiments are merely schematic, for example, the division of the units is only a logical function division, and there can be another division manner in actual implementation, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the displayed or discussed mutual couplings or direct couplings or communication connections between different units, or the among different units, can be indirect couplings or communication connections through some interfaces, devices or units, and can be in electrical, mechanical or other forms.

[0119] The units described as separated components can or can not be physically separated, and the components displayed as units can or can not be physical units, i.e., can be located in one place, or can be distributed on a plurality of network units. In actual implementation, some or all of the units can be selected according to actual needs to achieve the purposes of the embodiments of the present application.

[0120] In addition, each functional unit in the embodiments of the present application can be integrated in one processing unit, or each unit can exist physically as a separate unit, or two or more units can be integrated in one unit.

[0121] If the functions are implemented in the form of software function units and sold or used as independent products, they can be stored in a non-volatile computer readable storage medium executable by a processor. Based on this understanding, the technical solutions of the present application essentially or the part that contributes to the prior art or part of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes a plurality of instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in the embodiments of the present application. The aforementioned storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk, and various media that can store program codes.

[0122] Finally, it should be noted that the above-described embodiments are merely specific embodiments of the present application, which are used to illustrate the technical solutions of the present application, but not to limit the same. The protection scope of the present application is not limited thereto. Although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those skilled in the art that any skilled person in the art can still modify or easily think of changes to the technical solutions recorded in the foregoing embodiments, or make equivalent replacements to some of the technical features, within the technical scope disclosed by the present application. The modifications, changes or replacements do not make the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A cardiac magnetic dipole center trajectory parameter determination system, characterized by, The system comprises an acquisition device, a generation device, a determination device, a calculation device; wherein The acquisition device is configured to acquire a magnetocardiogram data set of a target wave band; the target wave band comprises at least one of a QRS wave band, a T wave band and a P wave band; The generation device is configured to generate a target isomagnetic map of the target wave band based on the magnetocardiogram data set of the target wave band; The determination device is configured to determine a target dipole center pixel point set of the target wave band based on a dipole position of each frame of the target isomagnetic map; The calculation device is configured to calculate a target dipole center trajectory parameter based on the target dipole center pixel point set and a pre-created standard center trajectory distribution atlas; the standard center trajectory distribution atlas is generated based on a plurality of sample magnetocardiogram data sets; and the target dipole center trajectory parameter is used to represent a distribution of the dipole center pixel point set in each pixel position level in the standard center trajectory distribution atlas.

2. The magnetocardiographic dipole center trajectory parameter determination system of claim 1, wherein, When the determination device is used to determine the target dipole center pixel point set of the target wave band based on a dipole position of each frame of the target isomagnetic map, the determination device is specifically configured to: For any frame of the isomagnetic map, determine the position coordinates of the positive pole and the position coordinates of the negative pole in the frame of the isomagnetic map based on the magnetic field intensity of the frame of the isomagnetic map; Determine the center point position coordinates of the dipole according to the position coordinates of the positive pole and the position coordinates of the negative pole; Determine the target dipole center pixel point set based on the center point position coordinates of the dipole determined in each frame of the isomagnetic map.

3. The magnetometric dipole center trajectory parameter determination system of claim 1, wherein, The system further comprises a creation device, and the creation device comprises a first acquisition module, a first determination module, a second determination module and a third determination module; The first acquisition module is configured to acquire a magnetocardiogram data set of a target wave band in a plurality of sample magnetocardiogram data sets; The first determination module is configured to determine a sample dipole center pixel point set based on the magnetocardiogram data set of the target wave band in the sample magnetocardiogram data set for any sample magnetocardiogram data set; The second determination module is configured to determine a density distribution map formed by each sample dipole center pixel point based on a plurality of sample dipole center pixel point sets; The third determination module is configured to determine the standard center trajectory distribution atlas based on the density distribution map.

4. The magnetometric dipole center trajectory parameter determination system of claim 3, wherein, When the third determination module is used to determine the standard center trajectory distribution atlas based on the density distribution map, the third determination module is specifically configured to: Determine an abnormal probability of each pixel position based on a total number of sample dipole center pixel points included in each pixel position in the density distribution map and label information of each sample dipole center pixel point; Generate a probability distribution map based on the abnormal probability of each pixel position in the density distribution map; Determine a position level of each pixel position based on the probability distribution map; Determine the standard center trajectory distribution atlas based on the position level of each pixel position.

5. The magnetometric dipole center trajectory parameter determination system of claim 1, wherein, The calculation device comprises a fourth determination module and a calculation module; The fourth determination module is configured to determine the number of dipole center pixel points of each pixel position level in the standard center trajectory distribution map. The calculation module is configured to calculate the target dipole center trajectory parameter based on the number of dipole center pixel points of each pixel position level and a preset dipole center trajectory parameter formula.

6. The magnetometric dipole center trajectory parameter determination system of claim 5, wherein, The preset dipole center trajectory parameter formula is as follows: Where dctj is the dipole center trajectory parameter, nc1 is the number of dipole center pixels in the pixel position level 1, nc2 is the number of dipole center pixels in the pixel position level 2, ..., nc n is the number of dipole center pixels in the pixel position level n, a1, a2, ..., a n They are nc1, nc2, ..., nc n The weighting coefficient of .

7. The magnetometric dipole center trajectory parameter determination system of claim 1, wherein, The acquisition device includes a second acquisition module and a fifth determination module. The second acquisition module is configured to acquire a target magnetocardiogram data set collected by the multi-channel magnetocardiograph. The fifth determination module is configured to determine a magnetocardiogram data set of a target wave band based on the target magnetocardiogram data set.

8. The magnetometric dipole center trajectory parameter determination system of claim 7, wherein, The number of channels of the multi-channel magnetocardiograph is one of 4 channels, 9 channels, 36 channels, and 64 channels.

9. The magnetometric dipole center trajectory parameter determination system of claim 7, wherein, When the generation device is used to generate a target isomagnetogram of the target wave band based on the magnetocardiogram data set of the target wave band, the generation device is specifically configured to: generate the target isomagnetogram of the target wave band based on the magnetocardiogram data set of the target wave band and the positions of each channel in the multi-channel magnetocardiograph.

10. The magnetometric dipole center trajectory parameter determination system of claim 4, wherein, The third determination module is specifically configured to determine the anomaly probability of each pixel position according to the following steps: For any pixel position, the total number of sample dipole center pixel points included in the pixel position and the target number of sample dipole center pixel points with abnormal types of label information are acquired; the ratio between the target number and the total number is determined as the anomaly probability of the pixel position.

Images