Apex calibration method, system, medium and electronic device based on magnetocardiogram

Abstract

The present invention provides a method, system, medium and electronic device for apex calibration based on magnetocardiogram, the method comprising the following steps: acquiring target data, the target data comprising magnetocardiogram data and angiography image data; performing gradient calculation based on the magnetocardiogram data to obtain a first apex point; reconstructing a heart model based on the angiography image, thereby obtaining a second apex point based on the heart model; performing coordinate conversion based on the first apex point and the second apex point to complete the calibration of the apex. The method, system, medium and electronic device for apex calibration based on magnetocardiogram of the present invention can solve the problems in the prior art of CT and MCG data calibration deviation caused by manual positioning, inability to use past CT and MCG data for retrospective analysis, and difficulty in batch processing of data calibration; at the same time, large quantities of data can also be automatically processed in batches, which is convenient and efficient, and can assist in the diagnosis of heart diseases and preoperative guidance.

Description

Technical Field

[0001] The present invention belongs to the technical field of image processing, and in particular relates to a method, system, medium and electronic equipment for apex calibration based on magnetocardiogram. Background Art

[0002] Magnetocardio Gram (MCG), as a non-invasive examination technology with no radiation, no contact and high spatial resolution, has become a research hotspot in the diagnosis of heart diseases. Computed Tomography Angiography (CTA) is based on the data acquisition of contrast agents passing through arteries, and has high sensitivity and specificity for the diagnosis of coronary artery lesions. MCG data and CTA data provide different functional and anatomical information respectively. By fusing them in multimodal way, their advantages can be comprehensively utilized to obtain more comprehensive information. However, due to the differences in the acquisition process of MCG data and CTA data and the specificity of the equipment itself, there are spatial misalignment and morphological differences between the two, which in turn affects the accuracy and reliability of source positioning after fusion. Therefore, it is necessary to align medical images of different modalities to achieve accurate spatial correspondence and provide more accurate information for medical diagnosis and treatment.

[0003] At present, there are two common solutions for the calibration of CTA data and MCG data. The first is to determine the spatial coordinates based on the scanning range of the two devices and the subject's suprasternal fossa. This method is called bony positioning. The disadvantage is that the manual bony positioning method may deviate due to different operating methods of the collector or incorrect posture of the subject. The second is to attach markers to the surface of the subject's chest skin and use external markers to calibrate the collected CT images and MCG images. The disadvantage is that this method requires individual operation for each subject and cannot be processed in large quantities. The existing data are generally marker-free CTA and MCG data that have been collected, and this type of data cannot be analyzed retrospectively. Summary of the invention

[0004] In view of the above-mentioned shortcomings of the prior art, an object of the present invention is to provide a method, system, medium and electronic device for apex calibration based on magnetocardiogram to solve the above-mentioned problems.

[0005] In a first aspect, the present invention provides a method for apex calibration based on magnetocardiogram, the method comprising the following steps:

[0006] Acquiring target data, wherein the target data includes magnetocardiogram data and contrast image data;

[0007] Performing gradient calculation based on the magnetocardiogram data to obtain a first apex point;

[0008] Reconstructing a heart model based on the angiography image, thereby obtaining a second apex point based on the heart model;

[0009] Coordinate transformation is performed based on the first apex point and the second apex point to complete the calibration operation of the apex.

[0010] In a possible implementation of the present application, performing gradient calculation based on the magnetocardiogram data to obtain the first apex point specifically includes:

[0011] Performing gradient calculation based on the magnetocardiogram data to obtain a point source position corresponding to each frame point source;

[0012] A lower right corner point within a target time period is identified based on the point source position to obtain the first apex point.

[0013] In a possible implementation of the present application, the method further includes: before calculating the gradient of the magnetocardiogram data, performing noise reduction preprocessing on the magnetocardiogram data.

[0014] In a possible implementation of the present application, the step of identifying the lower right corner point within the target time period based on the point source position to obtain the first apex point specifically includes:

[0015] Measuring magnetic field strength on the plane of the magnetic map;

[0016] Calculate the gradient based on the magnetic field strength and solve the target intensity point source corresponding to each moment;

[0017] The point with the largest horizontal coordinate and the smallest vertical coordinate in the target period is selected as the lower right corner point, so as to obtain the first apex point based on the target intensity point source corresponding to the lower right corner point. The specific coordinates are (MCG x ,MCG y ,MCG z ).

[0018] In a possible implementation of the present application, reconstructing a heart model based on the angiography image, thereby obtaining a second apex point based on the heart model, specifically includes:

[0019] Reconstructing a heart model based on the angiography image, wherein the plane corresponding to the horizontal and vertical coordinate axes of the heart model is the actual frontal plane;

[0020] Based on the upper left corner of the model as the origin, the distance between each point on the heart model and the origin is calculated;

[0021] Based on the point with the largest distance, the second apex point is obtained, and the specific coordinates are (CTA x ,CTA y,CTA z ).

[0022] In a possible implementation of the present application, the coordinate transformation based on the first apex point and the second apex point to complete the calibration operation of the apex specifically includes:

[0023] The coordinates of the first apex point and the second apex point are aligned and translated to obtain the target apex point of the calibrated heart model. The calibration formula is as follows:

[0024] heart′ x =heart x -CTA x +MCG x ;

[0025] heart′ y =heart y -CTA y +MCG y ;

[0026] Among them, heart x and heart y is the point coordinate of the heart model before calibration, heart′ x and heart′ y The coordinates of the points after calibration of the heart model.

[0027] In a possible implementation of the present application, the method further includes: before calculating the second apex point, segmenting and gridding the heart model.

[0028] In a second aspect, the present invention provides a system for apex calibration based on magnetocardiogram, the system comprising:

[0029] An acquisition module, used for acquiring target data, wherein the target data includes magnetocardiogram data and contrast image data;

[0030] A calculation module, configured to perform gradient calculation based on the magnetocardiogram data to obtain a first apex point;

[0031] A reconstruction module, used for reconstructing a heart model based on the angiography image, so as to obtain a second apex point based on the heart model;

[0032] The calibration module is used to perform coordinate conversion based on the first apex point and the second apex point to complete the calibration operation of the apex.

[0033] In a third aspect, the present invention provides an electronic device, the electronic device comprising: a processor and a memory;

[0034] The memory is used to store computer programs;

[0035] The processor is used to execute the computer program stored in the memory, so that the electronic device performs the above-mentioned magnetocardiogram-based apex calibration method.

[0036] In a fourth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by an electronic device, implements the above-mentioned magnetocardiogram-based apex calibration method.

[0037] As described above, the magnetocardiogram-based apex calibration method, system, medium and electronic device described in the present invention have the following beneficial effects: they can solve the problems in the prior art of CT and MCG data calibration deviation caused by manual positioning, inability to use past CT and MCG data for retrospective analysis, and difficulty in batch processing of data calibration; at the same time, large quantities of data can also be automatically processed in batches, which is convenient and efficient, and can assist in the diagnosis of heart disease and preoperative guidance. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] Figure 1 It is a flowchart showing the steps of a method for apex calibration based on magnetocardiography in one embodiment of the present invention;

[0039] Figure 2 It is a schematic diagram showing the process of apex calibration of the magnetocardiogram-based apex calibration method in one embodiment of the present invention;

[0040] Figure 3 It is a flowchart showing the steps of a method for apex calibration based on magnetocardiography in one embodiment of the present invention;

[0041] Figure 4 It is a flowchart showing the steps of a method for apex calibration based on magnetocardiography in one embodiment of the present invention;

[0042] Figure 5 It is a schematic diagram showing the bone positioning of the apex calibration method based on magnetocardiography in one embodiment of the present invention;

[0043] Figure 6 It is a schematic structural diagram of a magnetocardiogram-based apex calibration system in one embodiment of the present invention;

[0044] Figure 7 It is a schematic structural diagram of an electronic device of the present invention in one embodiment. DETAILED DESCRIPTION

[0045] The following describes the embodiments of the present invention by specific examples, and those skilled in the art can easily understand other advantages and effects of the present invention from the contents disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and the details in this specification can also be modified or changed in various ways based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that the following embodiments and features in the embodiments can be combined with each other without conflict.

[0046] It should be noted that the illustrations provided in the following embodiments are only schematic illustrations of the basic concept of the present invention, and thus the drawings only show components related to the present invention rather than being drawn according to the number, shape and size of components in actual implementation. In actual implementation, the type, quantity and proportion of each component may be changed arbitrarily, and the component layout may also be more complicated.

[0047] The present invention provides a method for calibrating the apex of the heart based on a magnetocardiogram, and the principle is that the electrical activity of the heart is closely related to the depolarization process in different bands. Specifically, after the ventricle starts to depolarize, a ventricular septum depolarization vector from left to right is first generated, and then the excitation is transmitted to the apex, and the ventricular wall near the apex begins to be excited through the left and right bundle branches and Purkinje fibers, and then the depolarization surface is depolarized on the ventricular wall through the apex, and the depolarization surface of the left ventricle is more obviously larger than that of the right ventricle, and then most of the right ventricle ends depolarization, while the left ventricle still maintains a relatively large depolarization surface and begins to point to the upper left side, and finally most of the ventricles are depolarized, and only a small part of the myocardium at the rear bottom of the left ventricle or the myocardium at the pulmonary artery root (vertebral body) of the right ventricle is still in depolarization. According to the above theory, the point source will pass through the apex during the QRS band reconstruction process, and the point with the largest gradient intensity of the magnetocardiogram at each moment can be regarded as the location of the point source activity. Therefore, combining the above two points, the gradient intensity point located in the lower right corner of the entire QRS band can be aligned with the apex in the CTA reconstructed heart model to complete the apex calibration.

[0048] The technical solutions in the embodiments of the present invention will be described in detail below in conjunction with the accompanying drawings in the embodiments of the present invention.

[0049] Specifically, see Figure 1 In one embodiment of the invention, the apex calibration method based on magnetocardiogram of the present invention comprises the following steps:

[0050] Step S102, acquiring target data, wherein the target data includes magnetocardiogram data and contrast image data;

[0051] Step S104, performing gradient calculation based on the magnetocardiogram data to obtain a first apex point;

[0052] Step S106, reconstructing a heart model based on the angiography image, thereby obtaining a second apex point based on the heart model;

[0053] Step S108: performing coordinate transformation based on the first apex point and the second apex point to complete the calibration operation of the apex.

[0054] It should be noted that, in this embodiment, Figure 2 As shown, it is a flow chart of apex calibration, wherein the target data is first acquired, and correspondingly, the target data includes magnetocardiogram data (MCG data, Magneto Cardio Gram) and angiography image data (CTA data, Computed Tomography Angiography), wherein the processing operation steps corresponding to each data are different, specifically, data denoising (noise reduction) preprocessing is specifically performed for the magnetocardiogram data, and heart model segmentation and gridding are specifically performed for the angiography image data, so as to obtain the first apex point and the second apex point respectively, and finally coordinate transformation is performed based on the first apex point and the second apex point to complete the calibration operation of the apex, wherein the coordinate transformation specifically includes translation alignment, and the gradient of the magnetocardiogram is used to find the corresponding apex position.

[0055] Furthermore, in one embodiment of the invention, the gradient calculation based on the magnetocardiogram data to obtain the first apex point specifically includes: performing a gradient calculation based on the magnetocardiogram data to obtain a point source position corresponding to a point source in each frame; and identifying a lower right corner point within a target time period based on the point source position to obtain the first apex point.

[0056] It should be noted that, in this embodiment, specifically, Figure 3 As shown, it is shown as a step flow chart, wherein, based on the point source position, identifying the lower right corner point in the target time period to obtain the first apex point specifically includes the following steps:

[0057] Step S302, measuring the magnetic field strength on the magnetic map plane;

[0058] Step S304, calculating the gradient based on the magnetic field strength and solving the target intensity point source corresponding to each moment;

[0059] Step S306 , selecting a point with the largest horizontal coordinate and the smallest vertical coordinate within the target time period as the lower right corner point, thereby obtaining the first apex point based on the target intensity point source corresponding to the lower right corner point.

[0060] It should be noted that, in this embodiment, before calculating the gradient of the magnetocardiogram data, the magnetocardiogram data is preprocessed for noise reduction. The gradient is calculated based on the magnetic map, and the point with the largest point source intensity on the magnetic map at each moment can be obtained. In the QR period, the point with the point source position at the lower right is considered to be the point passing through the apex. The coordinates of this point are aligned with the actual apex point to complete the calibration of the apex.

[0061] Specifically, assume that the magnetic field at any point on the plane is given by the coordinates (r x ,r y ,r z ) is generated by the magnetic dipole, defining the square of the magnetic field gradient The maximum value corresponds to the point on the plane where the point source intensity is the largest. Obviously, the spatial derivatives of the gradient q along the x-axis and y-axis at this point become "0" at the same time, and the corresponding formula is:

[0062]

[0063] Among them, B z is the magnetic induction intensity, x is the horizontal coordinate, y is the vertical coordinate, and the magnetic field gradient on the measurement plane The point with the maximum value is located above the point source, and its coordinates are the same as the dipole coordinates (r x ,r y ) overlap, solve the maximum intensity point source corresponding to each moment according to the gradient and record the coordinates. Therefore, for the three-dimensional heart model, within the QR period, the point with the maximum x coordinate and the minimum y coordinate of the point source can be obtained as the first apex point, and the specific coordinates are (MCG x ,MCG y ,MCG z ).

[0064] Furthermore, if Figure 4 As shown, it is shown as a step flow chart, wherein reconstructing a heart model based on the angiography image, thereby obtaining a second apex point based on the heart model, specifically includes the following steps:

[0065] Step S402, reconstructing a heart model based on the angiography image, wherein the plane corresponding to the horizontal and vertical coordinate axes of the heart model is the actual frontal plane;

[0066] Step S404, based on the upper left corner of the model as the origin, calculating the distance between each point on the heart model and the origin;

[0067] Step S406: obtaining the second apex point based on the point with the largest distance.

[0068] It should be noted that, in this embodiment, before calculating the second apex point, the heart model is segmented and gridded. Specifically, according to the heart model reconstructed from the CTA data, the plane corresponding to the horizontal and vertical coordinate axes (XY) of the heart model is the actual frontal plane, wherein the point at the upper left corner of the heart model is the origin, and thus the distance between each point on the heart model and the origin is calculated, wherein the point with the largest (farthest) distance from the origin is the actual apex on the heart model, i.e., the second apex point, and the specific coordinates are (CTA x ,CTA y ,CTA z ).

[0069] Furthermore, in one embodiment of the invention, the coordinate conversion based on the first apex point and the second apex point to complete the calibration operation of the apex specifically includes:

[0070] The coordinates of the first apex point and the second apex point are aligned and translated to obtain the target apex point of the calibrated heart model. The calibration formula is as follows:

[0071] heart′ x =heart x -CTA x +MCG x ;

[0072] heart′ y =heart y -CTA y +MCG y ;

[0073] Among them, heart x and heart y is the point coordinate of the heart model before calibration, heart′ x and heart′ y The coordinates of the points after calibration of the heart model.

[0074] It should be noted that, in this embodiment, after obtaining the first apex point and the second apex point, the coordinates of the two apexes are aligned and translated to complete the apex calibration based on the magnetocardiogram. The calibration formula is as follows: x =heart x -CTA x +MCG x ; and heart′ y =heart y -CTA y +MCG y Among them, heart x and heart yis the point coordinate of the heart model before calibration, heart′ x and heart′ y The coordinates of the points after calibration of the heart model.

[0075] Furthermore, in a specific application, the upper body of the subject is scanned by CT, and the obtained CTA image is preprocessed, the heart model is segmented and meshed to generate a personalized heart mesh model, and the upper body of the subject is scanned by a magnetocardiograph, wherein the upper boundary of the scanning range is "4 cm" below the suprasternal fossa, and the right boundary is "8 cm" to the right of the suprasternal fossa. According to the heart model reconstructed by CTA data, the point on the heart farthest from the origin is calculated as the actual apex (second apex point) on the heart model, and the coordinates are (CTA x ,CTA y ,CTA z ); and according to the collected MCG data, the gradient of the magnetic map of each frame is calculated, and the maximum intensity point source corresponding to each moment is solved according to the gradient and the coordinates are recorded. In the QR period, the point located in the lower right corner is considered to be the first apex point, and the coordinates are (MCG x ,MCG y ,MCG z ), thereby performing coordinate transformation of the heart model according to the coordinates of the two apex points and completing the calibration of the two data.

[0076] In addition, in addition to the apex calibration of this patent described in the above embodiments, this embodiment describes a comparative verification experiment, in which two subjects were selected, wherein the first is calibration using markers: three radiopaque acrylic markers and coil markers are attached to the chest surface of the subject, and the upper body of the subject with the markers is scanned by CT, and the obtained CTA image is segmented and gridded to generate a personalized heart grid model, and the MCG data containing the coil markers is scanned, and the heart grid model obtained by the MCG and CTA scans is physically calibrated using external markers; the second is calibration using clinical bone positioning: such as Figure 5 As shown in the figure, it is a schematic diagram of bony positioning. MCG data is obtained through magnetocardiograph scanning. The suprasternal fossa of the subject is used as the anatomical landmark for the alignment of the two examination images. The suprasternal fossa can be directly observed in the chest CT reconstructed image, and the relative position of the heart can be clearly seen. During the MCG acquisition process, laser positioning is used to ensure that the scanning range of all subjects maintains the same relative position to the suprasternal fossa. The upper boundary of the scanning range is b cm below the suprasternal fossa, and the right boundary is a cm to the right of the suprasternal fossa.

[0077] Finally, the multimodal fusion data calibrated by the above three methods were used to solve the inverse problem of the magnetic field using the volume conductor model, restore the position of the point source of each frame of the QRS band, and connect them into a trajectory. The shape of the QRS trajectory was compared, and the same fixed point in the heart was selected. The coordinates of the point under the three methods were recorded and compared. The results showed that the two subjects used almost the same apex calibration method and marker calibration, indicating that the apex calibration based on this patent can effectively calibrate the two data without the assistance of external markers, which is conducive to the retrospective study of historical CTA and MCG data. In addition, it can also automatically process batches of data at the same time, which is convenient and efficient, and solves the problems of CT and MCG data calibration deviation caused by manual positioning in the prior art, the inability to use past CT and MCG data for retrospective analysis, and the difficulty in batch processing of data calibration.

[0078] Furthermore, an embodiment of the present application also provides an apex calibration system based on a magnetocardiogram, which can implement the apex calibration method based on a magnetocardiogram described in the present application. However, the implementation device of the apex calibration method based on a magnetocardiogram described in the present application includes but is not limited to the structure of the apex calibration system based on a magnetocardiogram listed in the present embodiment. All structural deformations and replacements of the prior art made according to the principles of the present application are included in the protection scope of the present application.

[0079] See also Figure 6 In one embodiment, this embodiment provides a magnetocardiogram-based apex calibration system 60, the system comprising:

[0080] An acquisition module 61 is used to acquire target data, wherein the target data includes magnetocardiogram data and contrast image data;

[0081] A calculation module 62, configured to perform gradient calculation based on the magnetocardiogram data to obtain a first apex point;

[0082] A reconstruction module 63, configured to reconstruct a heart model based on the angiography image, thereby obtaining a second apex point based on the heart model;

[0083] The calibration module 64 is used to perform coordinate conversion based on the first apex point and the second apex point to complete the calibration operation of the apex.

[0084] It should be noted that, in this embodiment, the acquisition module 61 is specifically used to acquire target data when applied. Accordingly, the target data specifically includes magnetocardiogram data (MCG data, Magneto Cardio Gram) and contrast image data (CTA data, Computed Tomography Angiography). Among them, the magnetocardiogram data can be obtained by scanning the upper body of the subject with a magnetocardiograph, and the contrast image data can be obtained by scanning the upper body of the subject with a CT scan.

[0085] Further, in this embodiment, the calculation module 62 is specifically used to perform gradient calculation based on the magnetocardiogram data to obtain the first apex point, which includes measuring the magnetic field strength on the magnetocardiogram plane, calculating the gradient based on the magnetic field strength and solving the target intensity point source corresponding to each moment, and finally selecting the point with the largest horizontal axis coordinate and the smallest vertical axis coordinate as the lower right corner point within the target time period, so as to obtain the first apex point based on the corresponding target intensity point source of the lower right corner point, and the specific coordinates are (MCG x ,MCG y ,MCG z ).

[0086] Further, in this embodiment, the reconstruction module 63 is specifically used to reconstruct the heart model based on the contrast image, so as to obtain the second apex point based on the heart model, which includes reconstructing the heart model based on the contrast image, wherein the plane corresponding to the horizontal and vertical axes of the heart model is the actual frontal plane, and based on the upper left corner of the model as the origin, the distance between each point on the heart model and the origin is calculated, and finally the second apex point is obtained based on the point with the largest distance, and the specific coordinates are (CTA x ,CTA y ,CTA z ).

[0087] Further, in this embodiment, the calibration module 64 is used to perform coordinate conversion based on the first apex point and the second apex point to complete the calibration of the apex. After obtaining the first apex point and the second apex point, the coordinates of the two apexes are aligned and translated to complete the apex calibration based on the magnetocardiogram. The calibration formula is as follows: x =heart x -CTA x +MCG x ; and heart′ y =heart y -CTA y +MCG y Among them, heart x and hearty is the point coordinate of the heart model before calibration, heart′ x and heart′ y The coordinates of the points after calibration of the heart model.

[0088] Furthermore, since the specific implementation of this embodiment corresponds to the aforementioned method embodiment, the same details will not be repeated here, and those skilled in the art should also understand that Figure 6 The division of the modules in the embodiment is merely a division of logical functions, and in actual implementation, all or part of them can be integrated into one or more physical entities, and these modules can be implemented entirely in the form of software called by processing elements, or entirely in the form of hardware, or some modules can be implemented in the form of software called by processing elements, and some modules can be implemented in the form of hardware.

[0089] In the several embodiments provided by the present invention, it should be understood that the disclosed system, device or method can be implemented in other ways. For example, the device embodiments described above are only schematic, for example, the division of modules / units is only a logical function division, and there may be other division methods in actual implementation, such as multiple modules or units can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or modules or units, which can be electrical, mechanical or other forms.

[0090] The modules / units described as separate components may or may not be physically separated, and the components displayed as modules / units may or may not be physical modules, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules / units may be selected according to actual needs to achieve the purpose of the embodiments of the present invention. For example, the functional modules / units in the various embodiments of the present invention may be integrated into one processing module, or each module / unit may exist physically separately, or two or more modules / units may be integrated into one module / unit.

[0091] Those of ordinary skill in the art should further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in terms of function in the above description. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0092] The embodiment of the present invention also provides a computer-readable storage medium. A person of ordinary skill in the art can understand that all or part of the steps in the method for implementing the above embodiment can be completed by instructing the processor through a program, and the program can be stored in a computer-readable storage medium, and the storage medium is a non-transitory medium, such as a random access memory, a read-only memory, a flash memory, a hard disk, a solid-state hard disk, a magnetic tape, a floppy disk, an optical disc, and any combination thereof. The above storage medium can be any available medium that can be accessed by a computer or a data storage device such as a server or a data center that includes one or more available media integrated. The available medium can be a magnetic medium (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a digital video disc (DVD)), or a semiconductor medium (e.g., a solid-state disk (SSD)), etc.

[0093] An embodiment of the present invention further provides an electronic device, which includes a processor and a memory.

[0094] The memory is used to store computer programs.

[0095] The memory includes: ROM, RAM, disk, USB flash drive, memory card or CD and other media that can store program codes.

[0096] The processor is connected to the memory and is used to execute the computer program stored in the memory so as to enable the electronic device to perform the above-mentioned magnetocardiogram-based apex calibration method.

[0097] Preferably, the processor can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components.

[0098] like Figure 7 As shown, the electronic device of the present invention is in the form of a general computing device. The components of the electronic device may include but are not limited to: one or more processors or processing units 71, a memory 72, and a bus 73 connecting different system components (including the memory 72 and the processing unit 71).

[0099] Bus 73 represents one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, a processor or a local bus using any of a variety of bus architectures. By way of example, these architectures include, but are not limited to, an Industry Standard Architecture (ISA) bus, a Micro Channel Architecture (MAC) bus, an Enhanced ISA bus, a Video Electronics Standards Association (VESA) local bus, and a Peripheral Component Interconnect (PCI) bus.

[0100] Electronic devices typically include a variety of computer system readable media. These media can be any available media that can be accessed by the electronic device, including volatile and non-volatile media, removable and non-removable media.

[0101] The memory 72 may include computer system readable media in the form of volatile memory, such as random access memory (RAM) 721 and / or cache memory 722. The electronic device may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, the storage system 723 may be used to read and write non-removable, non-volatile magnetic media ( Figure 7 not shown, usually called a "hard drive"). Although Figure 7Not shown, a disk drive for reading and writing to a removable non-volatile disk (e.g., a "floppy disk"), and an optical disk drive for reading and writing to a removable non-volatile optical disk (e.g., a CD-ROM, a DVD-ROM, or other optical media) may be provided. In these cases, each drive may be connected to the bus 73 via one or more data medium interfaces. The memory 72 may include at least one program product having a set (e.g., at least one) of program modules that are configured to perform the functions of the various embodiments of the present invention.

[0102] A program / utility 724 having a set (at least one) of program modules 7241 may be stored, for example, in the memory 72, such program modules 7241 including but not limited to an operating system, one or more application programs, other program modules, and program data, each of which or some combination may include an implementation of a network environment. The program modules 7241 generally perform the functions and / or methods of the embodiments described herein.

[0103] The electronic device may also communicate with one or more external devices (e.g., keyboards, pointing devices, displays, etc.), one or more devices that enable a user to interact with the electronic device, and / or any device that enables the electronic device to communicate with one or more other computing devices (e.g., network cards, modems, etc.). Such communication may be performed via input / output (I / O) interface 74. Furthermore, the electronic device may also communicate with one or more networks (e.g., local area networks (LANs), wide area networks (WANs), and / or public networks, such as the Internet) via network adapter 75. Figure 7 As shown, the network adapter 75 communicates with other modules of the electronic device via the bus 73. It should be understood that, although not shown in the figure, other hardware and / or software modules can be used in conjunction with the electronic device, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.

[0104] The above embodiments are merely illustrative of the principles and effects of the present invention, and are not intended to limit the present invention. Anyone familiar with the art may modify or alter the above embodiments without departing from the spirit and scope of the present invention. Therefore, all equivalent modifications or alterations made by a person of ordinary skill in the art without departing from the spirit and technical concept disclosed by the present invention shall still be covered by the claims of the present invention.

Claims

1. A method for apex calibration based on magnetocardiogram, characterized in that: include: Acquiring target data, wherein the target data includes magnetocardiogram data and contrast image data; Performing gradient calculation based on the magnetocardiogram data to obtain a first apex point, specifically comprising performing gradient calculation based on the magnetocardiogram data to obtain a point source position corresponding to each frame point source; identifying a lower right corner point within a target time period based on the point source position to obtain the first apex point; Reconstructing a heart model based on the angiography image, thereby obtaining a second apex point based on the heart model; Coordinate transformation is performed based on the first apex point and the second apex point to complete the calibration operation of the apex.

2. The method for apex calibration based on magnetocardiogram according to claim 1, characterized in that: The method further includes: before calculating the gradient of the magnetocardiogram data, performing noise reduction preprocessing on the magnetocardiogram data.

3. The method for apex calibration based on magnetocardiogram according to claim 1, characterized in that: The step of identifying the lower right corner point within the target time period based on the point source position to obtain the first apex point specifically includes: Measuring magnetic field strength on the plane of the magnetic map; Calculate the gradient based on the magnetic field strength and solve the target intensity point source corresponding to each moment; The point with the largest horizontal coordinate and the smallest vertical coordinate in the target period is selected as the lower right corner point, so as to obtain the first apex point based on the target intensity point source corresponding to the lower right corner point. The specific coordinates are (MCG x , MCG y , MCG z ).

4. The method for apex calibration based on magnetocardiogram according to claim 3, characterized in that: The reconstructing the heart model based on the angiography image, thereby obtaining a second apex point based on the heart model, specifically includes: Reconstructing a heart model based on the angiography image, wherein the plane corresponding to the horizontal and vertical coordinate axes of the heart model is the actual frontal plane; Based on the upper left corner of the model as the origin, the distance between each point on the heart model and the origin is calculated; Based on the point with the largest distance, the second apex point is obtained, and the specific coordinates are (CTA x , CTA y , CTA z ).

5. The method for apex calibration based on magnetocardiogram according to claim 4, characterized in that: The coordinate conversion based on the first apex point and the second apex point to complete the calibration operation of the apex specifically includes: The coordinates of the first apex point and the second apex point are aligned and translated to obtain the target apex point of the calibrated heart model. The calibration formula is as follows: heart ′ x =heart x -CTA x +MCG x ; heart ′ y =heart y -CTA y +MCG y ; Among them, heart x and heart y is the point coordinate of the heart model before calibration, heart ′ x and heart ′ y The coordinates of the points after calibration of the heart model.

6. The method for apex calibration based on magnetocardiogram according to claim 4, characterized in that: The method further includes: before calculating the second apex point, segmenting and meshing the heart model.

7. A magnetocardiogram-based apex calibration system, characterized in that: include: An acquisition module, used for acquiring target data, wherein the target data includes magnetocardiogram data and contrast image data; A calculation module, configured to perform gradient calculation based on the magnetocardiogram data to obtain a first apex point, wherein the calculation specifically includes performing gradient calculation based on the magnetocardiogram data to obtain a point source position corresponding to a point source of each frame; and identifying a lower right corner point within a target time period based on the point source position to obtain the first apex point; A reconstruction module, used for reconstructing a heart model based on the angiography image, so as to obtain a second apex point based on the heart model; The calibration module is used to perform coordinate conversion based on the first apex point and the second apex point to complete the calibration operation of the apex.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the apex calibration method based on magnetocardiogram according to any one of claims 1 to 6 is implemented.

9. An electronic device, characterized in that: The electronic device comprises: a processor and a memory; wherein the memory is used to store a computer program, and the processor is used to execute the computer program stored in the memory, so that the electronic device performs the apex calibration method based on magnetocardiogram as claimed in any one of claims 1 to 6.

Images