A system for determining the center trajectory parameters of a cardiac magnetic dipole

Abstract

The application provides a heart magnetic dipole center trajectory parameter determination system, which comprises an acquisition device, a generation device, a determination device and a calculation device. The acquisition device is used to acquire heart magnetogram data sets of a target wave band in a target heart magnetogram data set. The generation device is used to generate a target isomagnetogram of the target wave band based on the heart magnetogram data sets of the target wave band. The determination device is used to determine a target dipole center pixel point set of the target wave band based on dipole positions of each frame of the isomagnetogram. The calculation device is used 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. In this way, the change rule of the heart magnetic dipole center can be accurately and intuitively described, and the dipole distribution state of different wave bands can be accurately quantified through the calculated dipole center trajectory parameters.

Description

Technical Field

[0001] This application relates to the field of magnetocardiography technology, and more specifically, to a system for determining the trajectory parameters of the center of a magnetic dipole. Background Technology

[0002] Magnetocardiography (MCG) is a non-invasive, radiation-free, and contactless technique for detecting cardiac electromagnetic function. It utilizes a highly sensitive magnetic field detector array to record the weak magnetic fields generated by the heart's electrical activity. The most distinctive feature of MCG is the dipole formed by the positive and negative magnetic poles. The relative position of these dipoles changes continuously throughout the entire cycle of cardiac electrical activity, exhibiting different patterns in different wavelength bands, such as dipole rotation in the QRS complex and a stable state in the T wave. However, directly analyzing the state of the dipoles is very difficult, and the changing patterns and phenomena at the dipole center are also hard to quantify directly. Summary of the Invention

[0003] In view of this, the purpose of this application is to provide a system and method for determining the trajectory parameters of the center of the magnetic dipole, which can accurately and intuitively describe the changing law of the center of the magnetic dipole, and can accurately quantify the distribution state of the dipole in different bands through the trajectory parameters of the center of the dipole.

[0004] This application provides a system for determining the center trajectory parameters of a magnetic dipole, the system comprising an acquisition device, a generation device, a determination device, and a calculation device; wherein... The acquisition device is used to acquire a magnetocardiogram dataset of a target band; the target band includes at least one band selected from the QRS band, T band, and P band. The generating device is used to generate a target isomagnetic map of the target band based on the magnetocardiogram dataset of the target band. The determining device is used to determine the target dipole center pixel set of the target band based on the dipole positions of each frame of the target isomagnetic map. The computing device is used to calculate the target dipole center trajectory parameters based on the target dipole center pixel set and a pre-created standard center trajectory distribution map; wherein the standard center trajectory distribution map is generated from multiple sample magnetocardiogram datasets; the target dipole center trajectory parameters are used to characterize the distribution of the dipole center pixel set in each pixel position level in the standard center trajectory distribution map.

[0005] Furthermore, when determining the target dipole center pixel set of the target band based on the dipole positions of each frame of the target isomagnetic map, the determining device is specifically used for: For any frame of isomagnetic map, the position coordinates of the positive pole and the position coordinates of the negative pole in the frame of isomagnetic map are determined based on the magnetic field strength of the frame of isomagnetic map; The position coordinates of the center point of the dipole are determined based on the position coordinates of the positive pole and the position coordinates of the negative pole. Based on the coordinates of the center point of the dipole determined by the isomagnetic maps of each frame, the set of center pixels of the target dipole is determined.

[0006] Furthermore, the system also includes a creation device, which includes a first acquisition module, a first determination module, a second determination module, and a third determination module; The first acquisition module is used to acquire the magnetocardiogram dataset of the target band from multiple sample magnetocardiogram datasets; The first determining module is used to determine the set of sample dipole center pixels for any sample magnetocardiogram dataset based on the magnetocardiogram dataset of the target band in the sample magnetocardiogram dataset; The second determining module is used to determine the density distribution map formed by the center pixels of each sample dipole based on the set of center pixels of multiple sample dipoles. The third determining module is used to determine the standard center trajectory distribution map based on the density distribution map.

[0007] Furthermore, when the third determining module is used to determine the standard center trajectory distribution map based on the density distribution map, the third determining module is specifically used for: Based on the total number of sample dipole center pixels included in each pixel location in the density distribution map, and the label information of each sample dipole center pixel, the anomaly probability of each pixel location is determined. A probability distribution map is generated based on the anomaly probability of each pixel position in the density distribution map; Based on the probability distribution map, the position level of each pixel location is determined; The standard center trajectory distribution map is determined based on the location level of each pixel position.

[0008] Furthermore, the computing device includes a fourth determining module and a computing module; The fourth determining module is used to determine the number of dipole center pixels at the pixel position level of the target dipole center pixels concentrated in the standard center trajectory distribution map; The calculation module is used to calculate the target dipole center trajectory parameters based on the number of dipole center pixels included in each pixel position level and the preset dipole center trajectory parameter formula.

[0009] Furthermore, the formula for the preset dipole center trajectory parameters is as follows: ; Where dctj represents the dipole center trajectory parameters. This represents the number of dipole center pixels in pixel location level 1. The number of dipole center pixels in pixel location level 2, ..., This represents the number of dipole center pixels in pixel location level n. , , ..., They are respectively , , ..., The weighting coefficients.

[0010] Furthermore, the acquisition device includes a second acquisition module and a fifth determination module; The second acquisition module is used to acquire the target magnetocardiogram dataset collected by the multi-channel magnetocardiogram instrument; The fifth determining module is used to determine the target band of magnetocardiogram dataset based on the target magnetocardiogram dataset.

[0011] Furthermore, the number of channels in a multi-channel magnetocardiograph is one of 4, 9, 36, or 64 channels.

[0012] Furthermore, when the generating device is used to generate a target isomagnetic map of the target band based on the magnetocardiogram dataset of the target band, the generating device is specifically used for: Based on the magnetocardiogram dataset of the target band and the positions of each channel in the multi-channel magnetocardiogram, a target isomagnetic map of the target band is generated.

[0013] Furthermore, the third determining module is specifically used to determine the anomaly probability of each pixel position according to the following steps: For any pixel location, obtain the total number of sample dipole center pixels included in the pixel location and the target number of sample dipole center pixels with abnormal label information; The ratio between the target quantity and the total quantity is determined as the anomaly probability of the pixel location.

[0014] This application provides a system for determining the center trajectory parameters of a magnetocardiogram dipole. Using a magnetocardiogram dataset of the target band, it determines the set of target dipole center pixels for that band. Based on this set of target dipole center pixels and a pre-created standard center trajectory distribution map, it calculates the target dipole center trajectory parameters. Compared to related technologies that cannot directly quantify the dipole distribution state in different bands, this application can accurately quantify the dipole distribution state in different bands using the dipole center trajectory parameters.

[0015] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description

[0016] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0017] Figure 1 This illustration shows one of the structural schematic diagrams of a magnetic dipole center trajectory parameter determination system provided in an embodiment of this application; Figure 2 This is a second schematic diagram of a system for determining the center trajectory parameters of a magnetic dipole provided in an embodiment of this application; Figure 3 This is shown as a third schematic diagram of a system for determining the center trajectory parameters of a magnetic dipole provided in an embodiment of this application; Figure 4 The fourth schematic diagram shows the structure of a system for determining the center trajectory parameters of a magnetic dipole provided in this application embodiment; Figure 5 A schematic diagram of the density distribution provided in an embodiment of this application is shown; Figure 6 A schematic diagram of the standard center trajectory distribution map provided in the embodiments of this application is shown. Detailed Implementation

[0018] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. The components of the embodiments of this application described and shown in the accompanying drawings can generally be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of this application provided in the accompanying drawings is not intended to limit the scope of the claimed application, but merely represents selected embodiments of this application. Based on the embodiments of this application, every other embodiment obtained by those skilled in the art without inventive effort falls within the scope of protection of this application.

[0019] First, the applicable application scenarios of this application are introduced. This application can be applied to the scenario of band data processing in magnetocardiography.

[0020] The most distinctive feature of magnetocardiography (MCG) lies in the dipoles formed by the positive and negative magnetic poles. The relative positions of these dipoles change continuously throughout the entire cycle of cardiac electrical activity, exhibiting different patterns in different wavebands, such as dipole rotation in the QRS complex and a stable state in the T wave. However, directly analyzing the state of the dipoles is very difficult, and the changing patterns and phenomena of the dipole center are also hard to quantify directly.

[0021] Please see Figures 1-4 , Figure 1 This is one of the structural schematic diagrams of a system for determining the center trajectory parameters of a magnetic dipole provided in this application. Figure 2 This is the second schematic diagram of a system for determining the center trajectory parameters of a magnetic dipole provided in this application. Figure 3 This is the third schematic diagram of a system for determining the center trajectory parameters of a magnetic dipole provided in this application. Figure 4 This is the fourth schematic diagram of a system for determining the center trajectory parameters of a magnetic dipole provided in this application. Figure 1 As shown in the illustration, an embodiment of this application provides a system 100 for determining the center trajectory parameters of a magnetic dipole. The system 100 includes an acquisition device 110, a generation device 120, a determination device 130, and a calculation device 140; wherein, The acquisition device 110 is used to acquire a magnetocardiogram dataset of a target band; the target band includes at least one band selected from the QRS band, T band, and P band. Here, we obtain the target band of the magnetocardiogram dataset from the target magnetocardiogram dataset for the entire cycle. The target band includes at least one of the following bands: QRS complex, T band, and P band. Among them, the QRS complex represents the electrical activity generated during ventricular depolarization and excitation, the T band represents the process of cardiac repolarization, that is, the process of the ventricle returning to a resting state from depolarization and excitation, which is usually the process of ventricular diastole, and the P band represents atrial electrical activity.

[0022] Here, the magnetocardiogram (MCC) dataset for the target band can be obtained after filtering and noise reduction. This ensures that the processed dataset is a high-quality dataset with almost no noise impact, thus guaranteeing the reliability of the obtained MCC dataset for the target band and improving the accuracy of the final calculated target dipole center trajectory parameters.

[0023] The generating device 120 is used to generate a target isomagnetic map of the target band based on the magnetocardiogram dataset of the target band. The determining device 130 is used to determine the target dipole center pixel set of the target band based on the dipole positions of each frame of the target isomagnetic map. In practical implementation, a target isomagnetic map of the target band can be generated based on the magnetocardiogram dataset of the target band after filtering and noise reduction. The target isomagnetic map is composed of multiple frames of isomagnetic maps, each frame corresponding to a dipole position. Based on the dipole positions of each frame of the target isomagnetic map—that is, based on the dipole positions corresponding to all frames of the target isomagnetic map—a set of target dipole center pixels in the target band can be determined. Here, each dipole position includes one positive dipole position and one negative dipole position.

[0024] The computing device 140 is used to calculate the target dipole center trajectory parameters based on the target dipole center pixel set and a pre-created standard center trajectory distribution map; wherein the standard center trajectory distribution map is generated from multiple sample magnetocardiogram datasets; the target dipole center trajectory parameters are used to characterize the distribution of the dipole center pixel set in each pixel position level in the standard center trajectory distribution map.

[0025] In practice, the center trajectory parameters of the target dipole can be calculated based on the obtained set of target dipole center pixels and the pre-created standard center trajectory distribution map.

[0026] Here, pixel position level refers to the pixel position assigned a level. This can be understood as the standard center trajectory distribution map being divided into multiple pixel positions, each with a corresponding level. Each pixel position can have the same number of levels. For example, in a standard center trajectory distribution map with 100 pixel positions, 20 are at level 1, 40 at level 2, 25 at level 3, and 15 at level 4.

[0027] It's important to note that by calculating the target dipole center trajectory parameters, the distribution state of the dipoles can be quantified. Specifically, the magnitude of these parameters directly reflects the distribution of the dipole center pixel set across different pixel positions within the standard center trajectory distribution map. In other words, the magnitude of the dipole center trajectory parameters determines whether the number of dipole center pixels at each pixel position is relatively high or low. Furthermore, it can also indicate whether the overall number of dipoles across different pixel positions is relatively high or low. For example, a larger value for the dipole center trajectory parameters indicates a greater number of dipole center pixels at a given pixel position within the standard center trajectory distribution map.

[0028] This application addresses the limitation of current technologies in directly quantifying the dipole states across different wavelength bands. Therefore, it uses a magnetocardiogram dataset of the target wavelength band to determine the target dipole center pixel set. Based on this target dipole center pixel set and a pre-created standard center trajectory distribution map, it calculates the target dipole center trajectory parameters. This allows for an accurate and intuitive description of the variation patterns of the magnetocardiogram dipole centers, and the calculated dipole center trajectory parameters enable accurate quantification of the dipole distribution states across different wavelength bands.

[0029] Furthermore, such as Figure 1 As shown, when determining the target dipole center pixel set of the target band based on the dipole positions of each frame of the target isomagnetic map, the determining device 130 is specifically used for: For any frame of isomagnetic map, the position coordinates of the positive pole and the position coordinates of the negative pole in the frame of isomagnetic map are determined based on the magnetic field strength of the frame of isomagnetic map; The position coordinates of the center point of the dipole are determined based on the position coordinates of the positive pole and the position coordinates of the negative pole. Based on the coordinates of the center point of the dipole determined by the isomagnetic maps of each frame, the set of center pixels of the target dipole is determined.

[0030] Here, for any given frame of isomagnetic map, which contains multiple magnetic field intensities, the coordinates of the largest magnetic field intensity are taken as the positive pole coordinates of the dipole, and the coordinates of the smallest magnetic field intensity are taken as the negative pole coordinates. Based on the positive and negative pole coordinates, the coordinates of the center point between them can be determined; this center point coordinate is the center point coordinate of the dipole. After determining the center point coordinates of the dipoles in each frame of isomagnetic map, these coordinates are combined to determine a set of target dipole center pixels. Here, the center point coordinates of one dipole can be composed of the center pixel coordinates of 1 to 4 dipoles.

[0031] It should be noted that the coordinates of the center point of a dipole are the coordinates of the midpoint of the straight line connecting the coordinates of the positive and negative poles of the dipole. These coordinates consist of an x-coordinate and a y-coordinate. The x-coordinate of the center point is obtained by taking the arithmetic mean of the x-coordinates of the positive and negative poles. Similarly, the y-coordinate is obtained by taking the arithmetic mean of the y-coordinates of the positive and negative poles. Thus, the coordinates of the center point of the dipole can be obtained. The coordinates of the center point of a dipole characterize the position of its center point.

[0032] Here, for example, for each frame of the isomagnetic map in the QRS band (taking a 100*100 tensor as an example), assuming there are three frames of isomagnetic maps in the target isomagnetic map, the coordinates of the center points of the three dipoles determined by these three frames of isomagnetic maps are P0=(56,72), P1=(48,65), P2=(67,77), so the set of center pixels of a target dipole corresponding to the QRS band is B={P0, P1, P2}={(56,72), (48,65), (67,77)}. Of course, if the center point coordinates of the dipole determined by the frame is not an integer, that is, if the coordinates contain decimals, for example, P0 = (56.5, 72.5), then the coordinates can be rounded up or down to obtain the center pixel coordinates of all dipoles covering the center point coordinates of the dipole. That is, the center point coordinates of the dipole are P0 = {(56, 72), (56, 73), (57, 72), (57, 73)}. Thus, the set of center pixels of a target dipole is B = {P0, P1, P2} = {(56, 72), (56, 73), (57, 72), (57, 73), (48, 65), (67, 77)}. Here, a set of target dipole center pixels consists of the coordinates of the center points of multiple dipoles. This can also be understood as the set of target dipole center pixels consisting of the coordinates of the center pixels of multiple target dipoles representing the center point of the dipole. Furthermore, if the coordinates of the center point of a dipole are integers, then the coordinates of the center pixels of that dipole are the center point coordinates of that dipole. If the coordinates of the center point of a dipole contain decimals, then the center point coordinates of that dipole are represented by the coordinates of the center pixels of 2 or 4 dipoles. Of course, any one of these 2 or 4 center pixel coordinates can describe or represent the corresponding dipole. It should be noted that the methods for determining a set of target dipole center pixels based on the coordinates of the center points of each dipole determined by the isomagnetic maps of each frame include, but are not limited to, the methods described above.

[0033] Furthermore, such as Figure 2 As shown, the system further includes a creation device 150, which comprises a first acquisition module 151, a first determination module 152, a second determination module 153, and a third determination module 154. The first acquisition module 151 is used to acquire the magnetocardiogram dataset of the target band from multiple sample magnetocardiogram datasets; The first determining module 152 is used to determine the set of sample dipole center pixels for any sample magnetocardiogram dataset based on the magnetocardiogram dataset of the target band in the sample magnetocardiogram dataset. The second determining module 153 is used to determine the density distribution map formed by the center pixels of each sample dipole based on the set of center pixels of multiple sample dipoles. The third determining module 154 is used to determine the standard center trajectory distribution map based on the density distribution map.

[0034] Here, one sample corresponds to one sample magnetocardiogram (MCG) dataset, and one sample MCG dataset corresponds to a target band MCG dataset. The target band MCG dataset is obtained from multiple sample MCG datasets. For any sample MCG dataset, based on the target band MCG dataset within that sample MCG dataset, the sample isomagnetic map of the target band can be obtained. Then, based on the dipole positions of each frame of the sample isomagnetic map, the center point positions of the sample dipoles corresponding to each frame of the target band isomagnetic map are determined. After removing the coordinates of the center pixels of duplicate dipoles, a set of sample dipole center pixels for the target band can be obtained. Based on multiple sets of dipole center pixels from multiple samples, and according to the coordinates of the center pixels of each sample dipole, all the sample dipole center pixels in the set are accumulated to the corresponding pixel positions in a tensor of a preset size. This generates a density distribution map composed of the center pixels of the dipoles for each sample, or, in other words, a density distribution map composed of the center points of the dipoles for each sample. Then, based on the density distribution map, a standard center trajectory distribution map is determined. Here, for example, the preset size tensor can be a 100*100 tensor.

[0035] Among them, such as Figure 5 The density distribution map shown is generated by summing all the sample dipole center points from multiple sample dipole center pixel sets corresponding to multiple samples under the target band, according to the coordinates of the center pixel positions of each sample dipole, and then adding them to the corresponding pixel positions in a tensor of a preset size. For example, if there are 1000 samples, there are 1000 sample dipole center pixel sets. Then, summing all the sample dipole center points from these 1000 sample dipole center pixel sets under the QRS band into a 100*100 tensor generates the density distribution map.

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

[0037] Furthermore, such as Figure 2 As shown, when the third determining module 154 is used to determine the standard center trajectory distribution map based on the density distribution map, the third determining module 154 is specifically used for: Based on the total number of sample dipole center pixels included in each pixel location in the density distribution map, and the label information of each sample dipole center pixel, the anomaly probability of each pixel location is determined. A probability distribution map is generated based on the anomaly probability of each pixel position in the density distribution map; Based on the probability distribution map, the position level of each pixel location is determined; The standard center trajectory distribution map is determined based on the location level of each pixel position.

[0038] Here, for any pixel location in the density distribution map, the total number of sample dipole center pixels included at that pixel location is counted. Simultaneously, the label information carried by each sample dipole center pixel at that location is examined to calculate the anomaly probability of each pixel location. After determining the anomaly probability of each pixel location, a probability distribution map can be generated. Then, according to a pre-defined rule corresponding to anomaly probabilities and levels, a positional level is assigned to each pixel location in the probability distribution map based on its determined anomaly probability. After determining the corresponding positional levels for each pixel location in the probability distribution map, the standard center trajectory distribution map is obtained.

[0039] 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 also known. Therefore, the label information of all sample dipole center pixels in the sample dipole center pixel set is also 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, such as 0 representing an abnormal type and 1 representing a normal type. If the label of a sample is 0, then the label information of all sample dipole center pixels in the sample dipole center pixel set is 0.

[0040] In the density distribution map, pixel locations typically fall into four categories: pixel location containing only abnormal objects, pixel location containing only normal objects, pixel location containing both normal and abnormal objects, and pixel location containing neither normal nor abnormal objects (blank pixel location). Specifically: a pixel location containing only abnormal objects means all sample dipole center pixels within that location have abnormal label information; a pixel location containing only normal objects means all sample dipole center pixels within that location have normal label information; a pixel location containing both normal and abnormal objects means the sample dipole center pixels within that location have both normal and abnormal label information; and a pixel location containing neither normal nor abnormal objects (blank pixel location) means there are no sample dipole center pixels within that location.

[0041] It should be noted that since there is no sample dipole center pixel at the blank pixel position, it is impossible to calculate the anomaly probability value. Therefore, the blank pixel position can be directly recorded as an anomaly probability of 0.

[0042] Here, for example, the pre-defined rule for corresponding anomaly probabilities and levels can be as follows: an anomaly probability range of [0, 0.2] corresponds to level 1, an anomaly probability range of (0.2, 0.4] corresponds to level 2, an anomaly probability range of (0.4, 0.6] corresponds to level 3, an anomaly probability range of (0.6, 0.8] corresponds to level 4, and an anomaly probability range of (0.8, 1] corresponds to level 5. Then, according to the above correspondence rule, the anomaly probability determined at each pixel position in the probability distribution map can be assigned to each pixel. The location levels of pixels are determined, and different levels of pixel locations can be represented by different colors. For example, levels 1 to 5 correspond to colors that increase in darkness, thus obtaining a standard center trajectory distribution map. Furthermore, the standard center trajectory distribution map shows that the darker the pixel location, the higher its level and the greater its anomaly probability. Here, the pre-defined rule corresponding to anomaly probabilities and levels can be extended, dividing the levels into n levels, with the anomaly probabilities corresponding to the n levels covering the probability range of 0-1, and the anomaly probabilities corresponding to levels 1, 2, 3, ..., n from low to high.

[0043] Among them, such as Figure 6The standard center trajectory distribution map shown can characterize the distribution of dipole center pixels, that is, it can characterize the overall distribution of dipoles. For example, the darker the color (the higher the probability of anomaly), the more concentrated the distribution of anomalous objects is at that pixel location. This means that the dipole center points representing anomalous objects are more concentrated at that pixel location, or the dipoles of the anomalous object are more concentrated at that pixel location. The color depth distribution of the standard center trajectory distribution map follows a certain pattern; the larger the sample size, the more obvious the pattern becomes (i.e., the more fixed the color depth positions). Here, an object can be considered a sample, and for an anomalous object, this can be represented by the sample's label information indicating an anomalous type.

[0044] Furthermore, such as Figure 3 As shown, the computing device 140 includes a fourth determining module 141 and a computing module 142: The fourth determining module 141 is used to determine the number of dipole center pixels at the pixel position level of the target dipole center pixels concentrated in the standard center trajectory distribution map. The calculation module 142 is used to calculate the target dipole center trajectory parameters based on the number of dipole center points included in each pixel position level and the preset dipole center trajectory parameter formula.

[0045] The number of target dipole center pixels is determined by the pixel position level of the target dipole center pixels concentrated in the standard center trajectory distribution map. Based on the number of dipole center points included in each level and the preset dipole center trajectory parameter formula, the target dipole center trajectory parameters can be calculated.

[0046] Furthermore, the formula for the preset dipole center trajectory parameters is as follows: ; Where dctj represents the dipole center trajectory parameters. This represents the number of dipole center pixels in pixel location level 1. The number of dipole center pixels in pixel location level 2, ..., This represents the number of dipole center pixels in pixel location level n. , , ..., They are respectively , , ..., The weighting coefficients.

[0047] The formula for the preset dipole center trajectory parameters can be determined based on the number of levels and a weighted approach. Here, + =0, + =0, + =0……, that is, the sum of all pairs of weighted coefficients whose indices add up to n+1 is 0. If n is odd, then =0, here, for example, if n is 5, then =a3=0.

[0048] Here, to facilitate understanding of the formula for the preset dipole center trajectory parameters, for example, assuming there are 5 levels, the formula for determining the preset dipole center trajectory parameters based on the number of levels can be specifically as follows: ; Where dctj represents the dipole center trajectory parameters. This represents the number of dipole center pixels in pixel location level 1. This represents the number of dipole center pixels in pixel location level 2. This represents the number of dipole center pixels in pixel location level 3. This represents the number of dipole center pixels in pixel location level 4. This represents the number of dipole center pixels in pixel position level 5.

[0049] Here, due to the pre-defined formula for the dipole center trajectory parameters... Let be the number of dipole center pixels in pixel location level 3. The anomaly probability corresponding to pixel location level 3 is usually around 0.5, for example (0.4, 0.6], which belongs to the "middle zone", that is, the anomaly probability is neither large nor small, so it is not weighted. The preceding weighting parameter is 0. When When, it is represented by the number of dipole center pixels in pixel position level 1. The number of dipole center pixels at pixel position level 2 is 5, the number of dipole center pixels at pixel position level 3 is 0, the number of dipole center pixels at pixel position level 4 is 7, and the number of dipole center pixels at pixel position level 5 is 5. It is easy to see that the sum of the number of dipole center pixels at pixel position level 4 and pixel position level 5 is 5 more than the sum of the number of dipole center pixels at pixel position level 1 and pixel position level 2. Then, [the text abruptly ends here]. Substituting these values ​​into the above formula for the preset dipole center trajectory parameters, we can obtain... =0.842. When At this point, it's easy to see that the sum of the number of dipole center pixels at pixel position level 4 and pixel position level 5 is 15 more than the sum of the number of dipole center pixels at pixel position level 1 and pixel position level 2. This is easily obtained. =2.19. Thus, the calculated values ​​for the dipole center trajectory parameters are 0.842 and 2.19. 2.19 is larger than 0.842. A larger value for the dipole center trajectory parameter indicates a larger difference between the sum of the number of dipole center pixels at pixel position level 4 and level 5, and the sum of the number of dipole center pixels at pixel position level 1 and level 2. In other words, a larger value for the dipole center trajectory parameter means a greater number of dipole center pixels at pixel position levels 4 and 5 compared to the number of target dipole center pixels concentrated in the standard center trajectory distribution map at pixel position levels 1 and 2. This suggests that the target dipole center pixels are more concentrated in levels 4 and 5, or that the overall distribution of the target dipoles is more concentrated in levels 4 and 5. Similarly, it can be easily verified that the smaller the value of the dipole center trajectory parameter, the more dipole center pixels of level 1 and level 2 are concentrated in the standard center trajectory distribution map, compared to the number of dipole center pixels of level 4 and 5. This means that the target dipole center pixels are more concentrated in level 1 and level 2. Furthermore, the overall distribution also shows that the target dipole is more concentrated in level 1 and level 2.

[0050] Here, for example, the calculated dipole center trajectory parameters represent the quantified result of the dipole distribution state. That is, the values ​​of these parameters allow for a direct and accurate determination of whether the target dipole center pixels are more concentrated in levels 4 and 5, or more concentrated in levels 1 and 2. Alternatively, it can indicate whether the number of dipoles in levels 4 and 5 is relatively greater or less than the number of dipoles in levels 1 and 2.

[0051] Here, when the preset dipole center trajectory parameter formula is: The value range of the dipole center trajectory parameter is between [-4, 4]. This value range of the dipole center trajectory parameter can be derived from the dipole center trajectory parameter formula.

[0052] Furthermore, such as Figure 4 As shown, the acquisition device 110 includes a second acquisition module 111 and a fifth determination module 112: The second acquisition module 111 is used to acquire the target magnetocardiogram dataset collected by the multi-channel magnetocardiogram instrument; The fifth determining module 112 is used to determine the target band of magnetocardiogram dataset based on the target magnetocardiogram dataset.

[0053] Here, the target magnetocardiogram (MCC) dataset for the entire cycle is acquired by a multi-channel magnetocardiograph. Then, the target MCC dataset is segmented into various bands to determine the target band's MCC dataset. Segmentation of each band can separate the QRS complex, T band, and other required bands, with each band having its corresponding MCC dataset.

[0054] Furthermore, the number of channels in a multi-channel magnetocardiograph is one of 4, 9, 36, or 64 channels.

[0055] Here, multichannel magnetocardiography (MCG) includes, but is not limited to, superconducting MCG and atomic magnetocardiography. Each type of multichannel MCG has a corresponding number of channels. The number of channels in a multichannel MCG can be any of 4, 9, 36, or 64 channels.

[0056] Furthermore, such as Figure 4 As shown, when the generating device 120 generates a target isomagnetic map of the target band based on the magnetocardiogram dataset of the target band, the generating device 120 is specifically used for: Based on the magnetocardiogram dataset of the target band and the positions of each channel in the multi-channel magnetocardiogram, a target isomagnetic map of the target band is generated.

[0057] Here, based on the magnetocardiogram dataset of the target band, as well as the magnetic field strength of each channel in the multi-channel magnetocardiogram and the position of each channel in the multi-channel magnetocardiogram, a target isomagnetic map of the target band is generated.

[0058] Furthermore, such as Figure 4 As shown, the third determining module 154 is further configured to determine the anomaly probability of each pixel position according to the following steps: For any pixel location, obtain the total number of sample dipole center pixels included in the pixel location and the target number of sample dipole center pixels with abnormal label information; The ratio between the target quantity and the total quantity is determined as the anomaly probability of the pixel location.

[0059] Here, in the density distribution map, for any pixel location, firstly, the total number of pixels in the sample dipole at that pixel location is obtained, and then the target number of sample dipole center pixels with the label information of anomalous type is obtained. Then, the ratio of the obtained target number of sample dipole center pixels with the label information of anomalous type to the total number of sample dipole center pixels can be used to obtain the anomalous probability corresponding to that pixel location.

[0060] Here, for example, for any pixel location, the total number of sample dipole center pixels in that pixel location is 10, and the target number of sample dipole center pixels with the label information of abnormal type is 8. Then the abnormal probability corresponding to that pixel location is 0.8.

[0061] This application provides a system for determining the center trajectory parameters of a magnetocardiogram (MCG) dipole. The system includes an acquisition device, a generation device, a determination device, and a calculation device. The acquisition device acquires a magnetocardiogram (MCG) dataset for a target band. The target band includes at least one of the QRS, T, and P bands. The generation device generates a target isomagnetic map of the target band based on the MCG dataset. The determination device determines the target dipole center pixel set for the target band based on the dipole positions of each frame of the isomagnetic map. The calculation device calculates the target dipole center trajectory parameters based on the target dipole center pixel set and a pre-created standard center trajectory distribution map. The standard center trajectory distribution map is generated from multiple sample MCG datasets. The target dipole center trajectory parameters characterize the distribution of the dipole center pixel set at each pixel position level in the standard center trajectory distribution map. In this way, the variation law of the magnetic dipole center can be accurately and intuitively described, and the dipole distribution state in different bands can be accurately quantified by using the calculated dipole center trajectory parameters.

[0062] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0063] In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. The apparatus embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the shown or discussed mutual couplings, direct couplings, or communication connections may be through some communication interfaces; indirect couplings or communication connections between devices or units may be electrical, mechanical, or other forms.

[0064] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0065] In addition, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0066] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a processor-executable, non-volatile, computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0067] Finally, it should be noted that the above-described embodiments are merely specific implementations of this application, used to illustrate the technical solutions of this application, and not to limit them. The scope of protection of this application is not limited thereto. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments, or make equivalent substitutions for some of the technical features, within the scope of the technology disclosed in this application. Such modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application, and should all be covered within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A system for determining the trajectory parameters of a magnetic dipole center, characterized in that, The system includes an acquisition device, a generation device, a determination device, a creation device, and a calculation device; the creation device includes a first acquisition module, a first determination module, a second determination module, and a third determination module; the calculation device includes a fourth determination module and a calculation module; wherein, The acquisition device is used to acquire a magnetocardiogram dataset of a target band; the target band includes at least one band selected from the QRS band, T band, and P band. The generating device is used to generate a target isomagnetic map of the target band based on the magnetocardiogram dataset of the target band. The determining device is used to determine the target dipole center pixel set of the target band based on the dipole positions of each frame of the target isomagnetic map. The computing device is used to calculate the target dipole center trajectory parameters based on the target dipole center pixel set and a pre-created standard center trajectory distribution map; wherein the standard center trajectory distribution map is generated from multiple sample magnetocardiogram datasets; the target dipole center trajectory parameters are used to characterize the distribution of the dipole center pixel set in each pixel position level in the standard center trajectory distribution map; The first acquisition module is used to acquire the magnetocardiogram dataset of the target band from multiple sample magnetocardiogram datasets; The first determining module is used to determine the set of sample dipole center pixels for any sample magnetocardiogram dataset based on the magnetocardiogram dataset of the target band in the sample magnetocardiogram dataset; The second determining module is used to, based on multiple sample dipole center pixel sets, according to the center pixel position coordinates of each sample dipole, accumulate all the sample dipole center pixels in the multiple sample corresponding dipole center pixel sets to the corresponding pixel positions in a tensor of a preset size, thereby generating a density distribution map composed of the center pixels of each sample dipole. The third determining module is used to determine the anomaly probability of each pixel location based on the total number of sample dipole center pixels included in each pixel location of the density distribution map and the label information of each sample dipole center pixel; for any pixel location, it obtains the total number of sample dipole center pixels included in the pixel location and the target number of sample dipole center pixels with anomaly type label information; the ratio between the target number and the total number is determined as the anomaly probability of the pixel location; a probability distribution map is generated based on the anomaly probabilities of each pixel location in the density distribution map; according to the pre-defined anomaly probability and level correspondence rules, a position level is assigned to each pixel location based on the anomaly probability determined in each pixel location in the probability distribution map; after determining the corresponding position level for each pixel location in the probability distribution map, the standard center trajectory distribution map is obtained. The fourth determining module is used to determine the number of dipole center pixels at the pixel position level of the target dipole center pixels concentrated in the standard center trajectory distribution map; The calculation module is used to calculate based on the number of dipole center pixels included in each pixel position level and a preset dipole center trajectory parameter formula. Calculate the trajectory parameters of the target dipole center; where, These are the parameters of the dipole center trajectory. This represents the number of dipole center pixels in pixel location level 1. The number of dipole center pixels in pixel location level 2, ..., This represents the number of dipole center pixels in pixel location level n. , , ..., They are respectively , , ..., The weighting coefficients.

2. The system for determining the trajectory parameters of the magnetic dipole center according to claim 1, characterized in that, When determining the target dipole center pixel set of the target band based on the dipole positions of each frame of the target isomagnetic map, the determining device is specifically used for: For any frame of isomagnetic map, the position coordinates of the positive pole and the position coordinates of the negative pole in the frame of isomagnetic map are determined based on the magnetic field strength of the frame of isomagnetic map; The position coordinates of the center point of the dipole are determined based on the position coordinates of the positive pole and the position coordinates of the negative pole. Based on the coordinates of the center point of the dipole determined by the isomagnetic maps of each frame, the set of center pixels of the target dipole is determined.

3. The system for determining the trajectory parameters of the magnetic dipole center according to claim 1, characterized in that, The acquisition device includes a second acquisition module and a fifth determination module; The second acquisition module is used to acquire the target magnetocardiogram dataset collected by the multi-channel magnetocardiogram instrument; The fifth determining module is used to determine the target band of magnetocardiogram dataset based on the target magnetocardiogram dataset.

4. The system for determining the trajectory parameters of the magnetic dipole center according to claim 3, characterized in that, The number of channels in a multi-channel magnetocardiograph is one of 4, 9, 36, or 64 channels.

5. The system for determining the trajectory parameters of the magnetic dipole center according to claim 3, characterized in that, When the generating device is used to generate a target isomagnetic map of the target band based on the magnetocardiogram dataset of the target band, the generating device is specifically used for: Based on the magnetocardiogram dataset of the target band and the positions of each channel in the multi-channel magnetocardiogram, a target isomagnetic map of the target band is generated.

Images