CN115880440B - Magnetic particle three-dimensional reconstruction imaging method based on generation countermeasure network - Google Patents

Magnetic particle three-dimensional reconstruction imaging method based on generation countermeasure network Download PDF

Info

Publication number
CN115880440B
CN115880440B CN202310047876.0A CN202310047876A CN115880440B CN 115880440 B CN115880440 B CN 115880440B CN 202310047876 A CN202310047876 A CN 202310047876A CN 115880440 B CN115880440 B CN 115880440B
Authority
CN
China
Prior art keywords
image
dimensional
mpi
images
representing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202310047876.0A
Other languages
Chinese (zh)
Other versions
CN115880440A (en
Inventor
田捷
张利文
杨帆
申钰松
尚亚欣
惠辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Automation of Chinese Academy of Science
Original Assignee
Institute of Automation of Chinese Academy of Science
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN202310047876.0A priority Critical patent/CN115880440B/en
Publication of CN115880440A publication Critical patent/CN115880440A/en
Application granted granted Critical
Publication of CN115880440B publication Critical patent/CN115880440B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

The invention belongs to the technical field of medical imaging, in particular relates to a magnetic particle three-dimensional reconstruction imaging method, system and device based on a generation countermeasure network, and aims to solve the problems of long reconstruction time and low imaging resolution of the existing MPI high-resolution image. The method comprises the following steps: acquiring a multi-view sparse two-dimensional MPI image of an imitation body of an object to be imaged and reconstructed, and performing downsampling to obtain a downsampled sparse two-dimensional MPI image; inputting the downsampled sparse two-dimensional MPI image into a trained generation model of a neighborhood point average diffusion search self-attention generation countermeasure network to obtain a dense two-dimensional MPI image; reconstructing each dense two-dimensional MPI image through a filtering back projection reconstruction algorithm to finally obtain a three-dimensional MPI image of the object to be imaged; the invention accelerates the reconstruction speed of the MPI image, improves the imaging resolution, and ensures that the magnetic particle imaging equipment has a larger application prospect in the medical field.

Description

Magnetic particle three-dimensional reconstruction imaging method based on generation countermeasure network
Technical Field
The invention belongs to the technical field of medical imaging, and particularly relates to a magnetic particle three-dimensional reconstruction imaging method, system and device based on a generation countermeasure network.
Background
Magnetic particle imaging is a novel imaging technology at present, and has the imaging advantages of no radiation, high sensitivity and the like. The current magnetic particle imaging equipment has wide application prospect in clinical research, including early detection of focus, drug targeting treatment and other leading directions. In the existing magnetic particle imaging devices, how to quickly realize the scanning of an object and the reconstruction of a high-resolution three-dimensional image is a challenging problem. At present, part of equipment performs fault scanning on an object by rotating an equipment coil, scans two-dimensional images of the object at different angles by physically rotating an imaged coil, and then performs filtered back projection reconstruction on the scanned images. However, the imaging technology needs to perform coil rotation movement, and the service life of equipment can be reduced due to too high rotation speed. Second, if the time resolution of the scan is to be improved, multiple scans are required, which is time consuming and laborious. Therefore, how to improve the reconstruction speed and the imaging resolution of imaging under the condition of considering the construction cost of hardware equipment is a current challenge. Based on the method, the invention provides a magnetic particle three-dimensional reconstruction imaging method based on generation of an antagonism network.
Disclosure of Invention
In order to solve the problems in the prior art, namely, the problems of long reconstruction time and low imaging resolution of the existing MPI three-dimensional high-resolution image, the invention provides a magnetic particle three-dimensional reconstruction imaging method based on generation of a countermeasure network, which comprises the following steps:
step S100, acquiring a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged, and performing downsampling to obtain a downsampled sparse two-dimensional MPI image; the multi-view sparse two-dimensional MPI images are a plurality of MPI projection images which are collected in sequence according to set angles in a rotating mode;
step S200, inputting the downsampled sparse two-dimensional MPI image into a trained generation model of a self-attention generation countermeasure network by average diffusion search of neighborhood points, and obtaining a dense two-dimensional MPI image;
step S300, reconstructing each dense two-dimensional MPI image by a filtering back projection reconstruction algorithm to finally obtain a three-dimensional MPI image of the object to be imaged;
the neighborhood point average diffusion search self-attention generation countermeasure network comprises a generation model and a discrimination model;
the generating model comprises five diffusion searching attention mechanism convolution modules; the diffusion search attention mechanism convolution module is constructed based on a neighborhood point average diffusion convolution sub-network, a self-attention fusion network and an activation function layer which are connected in sequence; the activation function adopted by the activation function layer in the first four diffusion searching attention mechanism convolution modules is a leakage ReLU, and the activation function adopted by the activation function layer in the fifth diffusion searching attention mechanism convolution module is a hyperbolic tangent function;
The neighborhood point average diffusion convolution sub-network is configured to sequentially encode the downsampled sparse two-dimensional MPI image through an encoding function to obtain a sequentially encoded image; performing regional convolution operation on each sequence coding image to obtain particle concentration information, angle information and association degree information among the images of each sequence coding image;
the self-attention fusion network is configured to multiply the particle concentration information, the angle information and the association degree information among the images of each sequence coding image with the corresponding weight matrix for linear transformation, fuse the images through a multi-head attention mechanism, and up-sample the fused images to obtain a dense two-dimensional MPI image;
the discrimination model comprises five convolutional neural network modules; the convolutional neural network module is constructed based on a convolutional layer, a regularization layer and an activation function layer which are sequentially connected; the activation function adopted by the activation function layer of the first four convolutional neural network modules is a leak ReLU, and the activation function adopted by the activation function layer of the fifth convolutional neural network module is a softmax.
In some preferred embodiments, the training process of the neighborhood point average spread search self-attention generation countermeasure network is as follows:
Step A100, acquiring a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged through MPI imaging equipment, and performing downsampling to obtain a downsampled sparse two-dimensional MPI image; obtaining a training data set based on the downsampled sparse two-dimensional MPI image and a corresponding truth value label thereof;
step A200, inputting the downsampled sparse two-dimensional MPI image into a pre-constructed neighborhood point average diffusion search self-attention generation countermeasure network generation model to obtain a dense two-dimensional MPI image;
step A300, inputting the dense two-dimensional MPI image and the truth value label corresponding to the dense two-dimensional MPI image into the neighborhood point average diffusion search self-attention generation countermeasure network discrimination model, and acquiring discrimination results of the dense two-dimensional MPI image;
step A400, based on the discrimination result, combining each dense two-dimensional MPI image, each sparse two-dimensional MPI image and the corresponding truth value labels thereof, calculating a total loss value through a pre-constructed loss function, and updating the network parameters of the neighborhood point average diffusion search self-attention generation countermeasure network generation model and discrimination model;
and step A500, training the generator network and the discriminator network of the neighborhood point average spread search self-attention generation countermeasure network in a circulating way until the trained neighborhood point average spread search self-attention generation countermeasure network is obtained.
In some preferred embodiments, the number of angles of formation is set to
Figure SMS_1
At an angle of
Figure SMS_2
At the time of the first
Figure SMS_3
Zhang Xishu two-dimensional MPI image has a coding function of
Figure SMS_4
Figure SMS_5
In some preferred embodiments, the sequential encoded images are subjected to a region convolution operation by: processing the concentration values of the pixel points in the sequential coding image, and performing operation by adopting a 3 multiplied by 3 convolution kernel after processing;
the method for processing the concentration value of the pixel point in the sequential coding image comprises the following steps:
Figure SMS_6
wherein ,
Figure SMS_9
representing the first of said sequential encoded pictures
Figure SMS_12
Line
1
Figure SMS_15
The pixel points of the column are arranged,
Figure SMS_8
representing a point neighborhood
Figure SMS_11
The average value of all the values within the range,
Figure SMS_13
representation of
Figure SMS_16
The minimum of all values within the neighborhood of points,
Figure SMS_7
representation of
Figure SMS_10
The maximum of all values within the neighborhood of points,
Figure SMS_14
representation of
Figure SMS_17
The variable values in the point neighborhood consist of points that are adjacent laterally, longitudinally and diagonally.
In some preferred embodiments, the particle concentration information, the angle information, and the correlation degree information between images of each sequentially encoded image are multiplied by a corresponding weight matrix to perform linear transformation, and the method includes:
Figure SMS_18
Figure SMS_19
wherein ,
Figure SMS_37
representing self-attention mechanism, vector
Figure SMS_40
The information on the concentration of the particles is indicated,
Figure SMS_43
the method comprises the steps of carrying out a first treatment on the surface of the Vector quantity
Figure SMS_23
The information of the angle is represented by a set of angles,
Figure SMS_27
the method comprises the steps of carrying out a first treatment on the surface of the Vector quantity
Figure SMS_30
Information indicating the degree of association between images,
Figure SMS_34
Figure SMS_32
representing the deformation of a matrix into a one-dimensional column vector
Figure SMS_36
Representing sequentially encoded images of the input neighborhood point average diffusion convolution sub-network;
Figure SMS_39
Figure SMS_42
Figure SMS_35
respectively representing weight matrices corresponding to particle concentration information, angle information and correlation degree information between images of sequentially encoded images,
Figure SMS_38
representing sequential encoded images
Figure SMS_41
Neurons in all output layers
Figure SMS_44
The sum of the output results from 1 to N,
Figure SMS_21
representing sequential encoded images
Figure SMS_24
Through the first
Figure SMS_28
The output values of the individual neurons are then,
Figure SMS_31
representing sequential encoded images
Figure SMS_20
Through the first
Figure SMS_26
The output values of the individual neurons are then,
Figure SMS_29
representing sequential encoded images
Figure SMS_33
Through the first
Figure SMS_22
The output value of the individual neurons is the total output value
Figure SMS_25
Is a specific gravity of (c).
In some preferred embodiments, the particle concentration information, the angle information and the association degree information between the images of each sequentially encoded image are fused through a multi-head attention mechanism, and the method comprises the following steps:
Figure SMS_45
Figure SMS_46
wherein ,
Figure SMS_47
is a function of the splice-in,
Figure SMS_48
represent the first
Figure SMS_49
The outputs of the individual attention machine sub-modules,
Figure SMS_50
a total weight matrix representing the output of the multi-headed attentiveness mechanism,
Figure SMS_51
and information obtained by fusing the particle concentration information, the angle information and the association degree information between the images of the sequential encoded images.
In some preferred embodiments, the total loss value is calculated by a pre-constructed loss function by:
Figure SMS_52
Figure SMS_53
Figure SMS_54
Figure SMS_55
wherein ,
Figure SMS_58
as a function of the total loss,
Figure SMS_62
representing a pixel loss function and an contrast loss function,
Figure SMS_65
representing an angular sequential ordering optimization function,
Figure SMS_57
representing the error of the pixel information and,
Figure SMS_63
representing the loss of the generative model,
Figure SMS_67
a sparse two-dimensional MPI image is represented,
Figure SMS_70
representation generation of the first
Figure SMS_59
Zhang Miji a two-dimensional MPI image,
Figure SMS_60
representation generation of the first
Figure SMS_64
Zhang Miji a two-dimensional MPI image,
Figure SMS_69
represent the first
Figure SMS_72
Zhang Xishu truth labels corresponding to two-dimensional MPI images,
Figure SMS_77
a coefficient representing the loss of a pixel,
Figure SMS_78
a coefficient representing the countering loss is provided,
Figure SMS_81
indicating the probability of whether the discrimination network discriminates true,
Figure SMS_74
the output of the generative model is represented,
Figure SMS_76
representing the function of the angle difference,
Figure SMS_79
as a function of the hyperbolic tangent,
Figure SMS_80
representing the predicted value of the angle at the corresponding input,
Figure SMS_56
Figure SMS_61
representing arbitrary two dense two-dimensional MPI images
Figure SMS_66
An angle judgment function is provided, wherein the angle judgment function is provided with a plurality of angle judgment functions,
Figure SMS_71
and
Figure SMS_68
respectively represent the first
Figure SMS_73
The angle labels corresponding to the densely-generated two-dimensional MPI images,
Figure SMS_75
parameters indicating the richness settings for the information contained in the different mimics.
In a second aspect of the invention, a magnetic particle three-dimensional reconstruction imaging system based on generation of a countermeasure network is presented, the system comprising: an image acquisition module 100, an image reconstruction module 200 and a three-dimensional reconstruction module 300;
The image acquisition module 100 is configured to acquire a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged, and perform downsampling to obtain a downsampled sparse two-dimensional MPI image; the multi-view sparse two-dimensional MPI images are a plurality of MPI projection images which are collected in sequence according to set angles in a rotating mode;
the image reconstruction module 200 is configured to input the downsampled sparse two-dimensional MPI image into a trained generation model of a neighborhood point average diffusion search self-attention generation countermeasure network to obtain a dense two-dimensional MPI image;
the three-dimensional reconstruction module 300 is configured to reconstruct each dense two-dimensional MPI image through a filtered back projection reconstruction algorithm, and finally obtain a three-dimensional MPI image of the reconstructed object to be imaged;
the neighborhood point average diffusion search self-attention generation countermeasure network comprises a generation model and a discrimination model;
the generating model comprises five diffusion searching attention mechanism convolution modules; the diffusion search attention mechanism convolution module is constructed based on a neighborhood point average diffusion convolution sub-network, a self-attention fusion network and an activation function layer which are connected in sequence; the activation function adopted by the activation function layer in the first four diffusion searching attention mechanism convolution modules is a leakage ReLU, and the activation function adopted by the activation function layer in the fifth diffusion searching attention mechanism convolution module is a hyperbolic tangent function;
The neighborhood point average diffusion convolution sub-network is configured to sequentially encode the downsampled sparse two-dimensional MPI image through an encoding function to obtain a sequentially encoded image; performing regional convolution operation on each sequence coding image to obtain particle concentration information, angle information and association degree information among the images of each sequence coding image;
the self-attention fusion network is configured to multiply the particle concentration information, the angle information and the association degree information among the images of each sequence coding image with the corresponding weight matrix for linear transformation, fuse the images through a multi-head attention mechanism, and up-sample the fused images to obtain a dense two-dimensional MPI image;
the discrimination model comprises five convolutional neural network modules; the convolutional neural network module is constructed based on a convolutional layer, a regularization layer and an activation function layer which are sequentially connected; the activation function adopted by the activation function layer of the first four convolutional neural network modules is a leak ReLU, and the activation function adopted by the activation function layer of the fifth convolutional neural network module is a softmax.
In a third aspect of the invention, a storage device is provided in which a plurality of programs are stored, said programs being adapted to be loaded and executed by a processor for implementing a magnetic particle three-dimensional reconstruction imaging method based on generating a countermeasure network as described above.
In a fourth aspect of the present invention, a processing arrangement is provided, comprising a processor, a storage device; a processor adapted to execute each program; a storage device adapted to store a plurality of programs; the program is adapted to be loaded and executed by a processor for carrying out one of the above-mentioned magnetic particle three-dimensional reconstruction imaging methods based on generating a countermeasure network.
The invention has the beneficial effects that:
the invention improves the reconstruction speed and imaging resolution of MPI imaging;
the invention integrates the average diffusion search self-attention mechanism of the neighborhood points into the generation countermeasure network, optimizes by utilizing the angle sequence ordering optimization function, accelerates the reconstruction speed of MPI images, improves the reconstruction quality of the MPI images, overcomes the hardware defects of the current equipment, and ensures that the magnetic particle imaging equipment has a larger application prospect in the medical field.
Drawings
Other features, objects and advantages of the present application will become more apparent upon reading of the detailed description of non-limiting embodiments, made with reference to the following drawings, in which:
FIG. 1 is a flow diagram of a method for three-dimensional reconstruction imaging of magnetic particles based on generation of a countermeasure network in accordance with an embodiment of the present invention;
FIG. 2 is a schematic diagram of a training flow of a magnetic particle three-dimensional reconstruction imaging method based on generating a countermeasure network in accordance with an embodiment of the present invention;
FIG. 3 is a schematic diagram of a network module of an attention mechanism of one embodiment of the present invention;
FIG. 4 is a schematic diagram of a network structure of a method for three-dimensional reconstruction imaging of magnetic particles based on generation of a countermeasure network in accordance with an embodiment of the present invention;
FIG. 5 is a schematic diagram of a frame of a magnetic particle three-dimensional reconstruction imaging system based on generation of a countermeasure network in accordance with an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The present application is described in further detail below with reference to the drawings and examples. It is to be understood that the specific embodiments described herein are merely illustrative of the invention and are not limiting of the invention. It should be noted that, for convenience of description, only the portions related to the present invention are shown in the drawings.
It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be combined with each other. The present application will be described in detail below with reference to the accompanying drawings in conjunction with embodiments.
The invention provides a magnetic particle three-dimensional reconstruction imaging method based on a generation countermeasure network, which is shown in figure 1 and comprises the following steps:
step S100, acquiring a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged, and performing downsampling to obtain a downsampled sparse two-dimensional MPI image; the multi-view sparse two-dimensional MPI images are a plurality of MPI projection images which are collected in sequence according to set angles in a rotating mode;
step S200, inputting the downsampled sparse two-dimensional MPI image into a trained generation model of a self-attention generation countermeasure network by average diffusion search of neighborhood points, and obtaining a dense two-dimensional MPI image;
step S300, reconstructing each dense two-dimensional MPI image by a filtering back projection reconstruction algorithm to finally obtain a three-dimensional MPI image of the object to be imaged;
the neighborhood point average diffusion search self-attention generation countermeasure network comprises a generation model and a discrimination model;
the generating model comprises five diffusion searching attention mechanism convolution modules; the diffusion search attention mechanism convolution module is constructed based on a neighborhood point average diffusion convolution sub-network, a self-attention fusion network and an activation function layer which are connected in sequence; the activation function adopted by the activation function layer in the first four diffusion searching attention mechanism convolution modules is a leakage ReLU, and the activation function adopted by the activation function layer in the fifth diffusion searching attention mechanism convolution module is a hyperbolic tangent function;
The neighborhood point average diffusion convolution sub-network is configured to sequentially encode the downsampled sparse two-dimensional MPI image through an encoding function to obtain a sequentially encoded image; performing regional convolution operation on each sequence coding image to obtain particle concentration information, angle information and association degree information among the images of each sequence coding image;
the self-attention fusion network is configured to multiply the particle concentration information, the angle information and the association degree information among the images of each sequence coding image with the corresponding weight matrix for linear transformation, fuse the images through a multi-head attention mechanism, and up-sample the fused images to obtain a dense two-dimensional MPI image;
the discrimination model comprises five convolutional neural network modules; the convolutional neural network module is constructed based on a convolutional layer, a regularization layer and an activation function layer which are sequentially connected; the activation function adopted by the activation function layer of the first four convolutional neural network modules is a leak ReLU, and the activation function adopted by the activation function layer of the fifth convolutional neural network module is a softmax.
In order to more clearly describe the three-dimensional reconstruction imaging method of magnetic particles based on the generation of the countermeasure network, each step in the embodiment of the method of the present invention is described in detail below with reference to the accompanying drawings.
In the existing commercial equipment for realizing magnetic particle (Magnetic particle Imaging, MPI) imaging, in order to image the distribution of particles in an object, the object needs to be scanned, but the current physical scanning mode based on coil rotation has long time for reconstructing an image, and multi-angle tomography is limited by the rotation speed of the coil, so that the hardware cost of the equipment is greatly increased, and the service life of the equipment is shortened. In order to effectively reduce reconstruction time and save hardware cost, the invention designs a Multi-angle three-dimensional magnetic particle imaging (Multi-view Three Dimensional Magnetic Particle Imaging, MV-3D-MPI) method and system based on a generated discrimination model. The invention provides a neighborhood point average diffusion search self-attention generation countermeasure network model, which aims to realize rapid reconstruction of an MPI high-resolution image and improve imaging resolution.
In the following embodiments, a training process of generating an countermeasure network by average spread search self-attention of a neighborhood point is described first, and then a testing process of generating an countermeasure network by average spread search self-attention of a neighborhood point is described in detail.
1. The training process of the neighborhood point average diffusion search self-attention generation countermeasure network comprises the following specific steps:
Step A100, acquiring a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged through MPI imaging equipment, and performing downsampling to obtain a downsampled sparse two-dimensional MPI image; obtaining a training data set based on the downsampled sparse two-dimensional MPI image and a corresponding truth value label thereof;
in the embodiment, acquiring multi-view dense MPI projection images by an own MPI device, and taking 15000 images as a training set; for a simulation, an image is preferably acquired every 5 ° (in other embodiments, the images can be acquired at set angle intervals according to actual requirements), and a total of 36-angle projection MPI images are acquired; the pixel size is changed to 12×12 by any conventional downsampling method, and the downsampled image is used as an input object of the network.
Step A200, inputting the downsampled sparse two-dimensional MPI image into a pre-constructed neighborhood point average diffusion search self-attention generation countermeasure network generation model to obtain a dense two-dimensional MPI image;
in this embodiment, the neighborhood point average diffusion search self-attention generation countermeasure network includes a generation model and a discrimination model;
the generating model comprises five diffusion searching attention mechanism convolution modules; the diffusion search attention mechanism convolution module is constructed based on a neighborhood point average diffusion convolution sub-network, a self-attention fusion network and an activation function layer which are connected in sequence; the activation function adopted by the activation function layer in the first four diffusion searching attention mechanism convolution modules is a leakage ReLU, and the activation function adopted by the activation function layer in the fifth diffusion searching attention mechanism convolution module is a hyperbolic tangent function;
The neighborhood point average diffusion convolution sub-network is configured to sequentially encode the downsampled sparse two-dimensional MPI image through an encoding function to obtain a sequentially encoded image; performing regional convolution operation on each sequence coding image to obtain particle concentration information, angle information and association degree information among the images of each sequence coding image;
the self-attention fusion network is configured to multiply the particle concentration information, the angle information and the association degree information among the images of each sequence coding image with the corresponding weight matrix for linear transformation, fuse the images through a multi-head attention mechanism, and up-sample the fused images to obtain a dense two-dimensional MPI image;
in the invention, because of the sequential relation of angles existing in the input, in order to enable the model to learn the relevance of adjacent angles and the global relevance of different angles, the input layer of the network adopts sequential coding, the input of the set model comprises 6 images of different angles, one image is taken every 30 degrees in [0 degree, 180 degrees ], and the angle is
Figure SMS_82
At the time of the first
Figure SMS_83
Zhang Xishu encoding function of two-dimensional MPI image
Figure SMS_84
The method comprises the following steps:
Figure SMS_85
wherein ,
Figure SMS_86
the sine and cosine functions in the expression can ensure that the constructed model learns the relative position relation between different images.
Aiming at the characteristics of an MPI image, the neighborhood point diffusion search value is mainly used for more fully modeling the relation between a certain pixel value of the MPI image and the neighborhood thereof, and as the magnetic particle response of a magnetic field free region can be influenced by the particles beside, the concentration value of imaging at part of points is related to the concentration value of other peripheral pixel points. Performing regional convolution operation on the sequence coding image, wherein the method comprises the following steps: processing the concentration values of the pixel points in the sequential coding images, and preferably adopting a 3 multiplied by 3 convolution kernel to operate after the processing;
the method for processing the concentration value of the pixel point in the sequential coding image comprises the following steps:
Figure SMS_87
to be used for
Figure SMS_88
The method for processing the concentration value by taking the point as an example comprises the following steps:
Figure SMS_89
wherein ,
Figure SMS_92
representing the first of said sequential encoded pictures
Figure SMS_93
Line
1
Figure SMS_97
The pixel points of the column are arranged,
Figure SMS_91
representing a point neighborhood
Figure SMS_95
The average value of all the values within the range,
Figure SMS_98
representation of
Figure SMS_100
The minimum of all values within the neighborhood of points,
Figure SMS_90
representation of
Figure SMS_94
The maximum of all values within the neighborhood of points,
Figure SMS_96
representation of
Figure SMS_99
The variable values in the point neighborhood consist of points that are adjacent laterally, longitudinally and diagonally.
The invention provides a self-attention fusion network module for fusing key information extracted from images by different sub-networks in a network model generation part. Setting the input of each sub-network as the preprocessed image
Figure SMS_101
Processing input by three sub-networks to obtain different types of information, multiplying the obtained three types of different information with corresponding weight matrixes respectively for linear transformation, wherein the corresponding weight matrixes are respectively
Figure SMS_102
Figure SMS_103
Figure SMS_104
Then fusing different types of information through a self-attention fusion network module, multiplying the particle concentration information, the angle information and the association degree information among the images of each sequence coding image with a corresponding weight matrix for linear transformation, wherein the method comprises the following steps:
Figure SMS_105
Figure SMS_106
wherein ,
Figure SMS_115
representing self-attention mechanism, vector
Figure SMS_119
The information on the concentration of the particles is indicated,
Figure SMS_123
the method comprises the steps of carrying out a first treatment on the surface of the Vector quantity
Figure SMS_109
The information of the angle is represented by a set of angles,
Figure SMS_114
the method comprises the steps of carrying out a first treatment on the surface of the Vector quantity
Figure SMS_117
Information indicating the degree of association between images,
Figure SMS_120
Figure SMS_110
representing the deformation of a matrix into a one-dimensional column vector
Figure SMS_111
Representing sequentially encoded images of the input neighborhood point average diffusion convolution sub-network;
Figure SMS_118
Figure SMS_122
Figure SMS_125
respectively representing weight matrices corresponding to particle concentration information, angle information and correlation degree information between images of sequentially encoded images,
Figure SMS_126
representing sequential encoded images
Figure SMS_128
Neurons in all output layers
Figure SMS_130
The sum of the output results from 1 to N,
Figure SMS_124
representing sequential encoded images
Figure SMS_127
Through the first
Figure SMS_129
The output values of the individual neurons are then,
Figure SMS_131
representing sequential encoded images
Figure SMS_107
Through the first
Figure SMS_113
The output values of the individual neurons are then,
Figure SMS_116
representing sequential encoded images
Figure SMS_121
Through the first
Figure SMS_108
The output value of the individual neurons is the total output value
Figure SMS_112
Is a specific gravity of (c).
The particle concentration information, the angle information and the association degree information among the images of each sequence coding image are fused through a multi-head attention mechanism, and the method comprises the following steps:
Figure SMS_132
Figure SMS_133
Figure SMS_134
wherein ,
Figure SMS_136
is a function of the splice-in,
Figure SMS_138
represent the first
Figure SMS_140
The outputs of the individual attention machine sub-modules,
Figure SMS_137
a total weight matrix representing the output of the multi-headed attentiveness mechanism,
Figure SMS_139
information representing the combination of particle concentration information, angle information and correlation information between images of sequentially encoded images, and the output of each attention sub-module is spliced and then passed through a processing unit
Figure SMS_141
Multiplying to perform linear transformation, and finally obtaining output of the generated model through up-sampling
Figure SMS_142
Figure SMS_135
Representing the upsampling function.
Step A300, inputting the dense two-dimensional MPI image and the truth value label corresponding to the dense two-dimensional MPI image into the neighborhood point average diffusion search self-attention generation countermeasure network discrimination model, and acquiring discrimination results of the dense two-dimensional MPI image;
the discrimination model comprises five convolutional neural network modules; the convolutional neural network module is constructed based on a convolutional layer, a regularization layer and an activation function layer which are sequentially connected; the activation function adopted by the activation function layer of the first four convolutional neural network modules is a leak ReLU, and the activation function adopted by the activation function layer of the fifth convolutional neural network module is a softmax.
The convolution layer carries out convolution operation on input data of the judging model, regularization is carried out through the regularization layer, and finally the activation function layer is input to carry out function operation.
Step A400, based on the discrimination result, combining each dense two-dimensional MPI image, each sparse two-dimensional MPI image and the corresponding truth value labels thereof, calculating a total loss value through a pre-constructed loss function, and updating the network parameters of the neighborhood point average diffusion search self-attention generation countermeasure network generation model and discrimination model;
in this embodiment, the total loss value is calculated by a pre-constructed loss function, and the method is as follows:
Figure SMS_143
Figure SMS_144
Figure SMS_145
Figure SMS_146
wherein ,
Figure SMS_164
as a function of the total loss,
Figure SMS_166
representing a pixel loss function and an contrast loss function,
Figure SMS_168
representing an angular sequential ordering optimization function,
Figure SMS_148
representing the error of the pixel information and,
Figure SMS_152
representing the loss of the generative model,
Figure SMS_157
a sparse two-dimensional MPI image is represented,
Figure SMS_160
representation generation of the first
Figure SMS_149
Zhang Miji a two-dimensional MPI image,
Figure SMS_151
representation generation of the first
Figure SMS_154
Zhang Miji a two-dimensional MPI image,
Figure SMS_161
represent the first
Figure SMS_167
Zhang Xishu truth labels corresponding to two-dimensional MPI images,
Figure SMS_169
a coefficient representing the loss of a pixel,
Figure SMS_173
representing the coefficient of countermeasures against losses, in the present embodiment
Figure SMS_175
And
Figure SMS_163
it is preferably set to 1 and,
Figure SMS_165
indicating the probability of whether the discrimination network discriminates true,
Figure SMS_170
The output of the generative model is represented,
Figure SMS_172
representing the function of the angle difference,
Figure SMS_147
as a function of the hyperbolic tangent,
Figure SMS_153
representing the predicted value of the angle at the corresponding input,
Figure SMS_156
Figure SMS_159
representing arbitrary two dense two-dimensional MPI images
Figure SMS_150
The angle judgment function takes a value of 1 if the judgment value is true, takes a value of 0 if the judgment value is not true,
Figure SMS_155
and
Figure SMS_158
respectively represent the first
Figure SMS_162
An angle label corresponding to the generated dense two-dimensional MPI image;
Figure SMS_171
parameters indicating the richness settings of the information contained in the different mimics, in the present invention, if the collected mimics are letter mimics, for the letter mimics
Figure SMS_174
Preferably set to 1, for more complex Chinese character mimics
Figure SMS_176
Preferably set to 2, for a simulation dedicated to evaluating resolution, since the details are rich,
Figure SMS_177
preferably set to 3. And when the image pairs are taken, taking the imitation as a unit, supervising the image pairs taken by each imitation image, and punishing error sequencing. The three sets of imitation body fills preferably have an inner diameter of 0.5mm containing 10 mg/ml Perimag solution. The field of view (FOV) of the three groups of simulations is preferably 6cm ×6× 6cm ×6cm (letters), 6cm ×6cm ×10× 10 cm (spiral, acquisition time of one sparse projection is 3 minutes), 6cm×6cm×cm (container, acquisition time of one sparse projection is 2 minutes), respectively.
Step A500, training a generator network and a discriminator network of the neighborhood point average spread search self-attention generation countermeasure network in a circulating manner until the trained neighborhood point average spread search self-attention generation countermeasure network is obtained;
2. test procedure for generating countermeasure network by neighborhood point average diffusion search self-attention
Step S100, acquiring a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged, and performing downsampling to obtain a downsampled sparse two-dimensional MPI image; the multi-view sparse two-dimensional MPI images are a plurality of MPI projection images which are collected in sequence according to set angles in a rotating mode;
in this embodiment, 200 instance images acquired by the MPI apparatus are taken as a test set.
Step S200, inputting the downsampled sparse two-dimensional MPI image into a trained generation model of a self-attention generation countermeasure network by average diffusion search of neighborhood points, and obtaining a dense two-dimensional MPI image;
in this embodiment, the self-attention generation countermeasure network is searched based on the average diffusion of the trained neighborhood points, and preferably, the sparse two-dimensional MPI image is generated from any angle of [0 °,5 °,10 °,15 °,20 °,25 °, 30 °,35 °,40 °,45 °,50 °,55 °,60 °,65 °, …,170 °,175 ° ], and the pixel size of the generated dense two-dimensional MPI image is 48×48.
Step S300, reconstructing each dense two-dimensional MPI image by a filtering back projection reconstruction algorithm to finally obtain a three-dimensional MPI image of the object to be imaged;
a magnetic particle three-dimensional reconstruction imaging system based on generation of a countermeasure network according to a second embodiment of the present invention, as shown in fig. 4, includes: an image acquisition module 100, an image reconstruction module 200 and a three-dimensional reconstruction module 300;
the image acquisition module 100 is configured to acquire a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged, and perform downsampling to obtain a downsampled sparse two-dimensional MPI image; the multi-view sparse two-dimensional MPI images are a plurality of MPI projection images which are collected in sequence according to set angles in a rotating mode;
the image reconstruction module 200 is configured to input the downsampled sparse two-dimensional MPI image into a trained generation model of a neighborhood point average diffusion search self-attention generation countermeasure network to obtain a dense two-dimensional MPI image;
the three-dimensional reconstruction module 300 is configured to reconstruct each dense two-dimensional MPI image through a filtered back projection reconstruction algorithm, and finally obtain a three-dimensional MPI image of the reconstructed object to be imaged;
the neighborhood point average diffusion search self-attention generation countermeasure network comprises a generation model and a discrimination model;
The generating model comprises five diffusion searching attention mechanism convolution modules; the diffusion search attention mechanism convolution module is constructed based on a neighborhood point average diffusion convolution sub-network, a self-attention fusion network and an activation function layer which are connected in sequence; the activation function adopted by the activation function layer in the first four diffusion searching attention mechanism convolution modules is a leakage ReLU, and the activation function adopted by the activation function layer in the fifth diffusion searching attention mechanism convolution module is a hyperbolic tangent function;
the neighborhood point average diffusion convolution sub-network is configured to sequentially encode the downsampled sparse two-dimensional MPI image through an encoding function to obtain a sequentially encoded image; performing regional convolution operation on each sequence coding image to obtain particle concentration information, angle information and association degree information among the images of each sequence coding image;
the self-attention fusion network is configured to multiply the particle concentration information, the angle information and the association degree information among the images of each sequence coding image with the corresponding weight matrix for linear transformation, fuse the images through a multi-head attention mechanism, and up-sample the fused images to obtain a dense two-dimensional MPI image;
The discrimination model comprises five convolutional neural network modules; the convolutional neural network module is constructed based on a convolutional layer, a regularization layer and an activation function layer which are sequentially connected; the activation function adopted by the activation function layer of the first four convolutional neural network modules is LeakyReLU, and the activation function adopted by the activation function layer of the fifth convolutional neural network module is softmax.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working processes and related descriptions of the above-described system may refer to corresponding processes in the foregoing method embodiments, which are not repeated herein.
It should be noted that, in the magnetic particle three-dimensional reconstruction imaging system provided in the foregoing embodiment, only the division of the foregoing functional modules is illustrated, and in practical application, the foregoing functional allocation may be performed by different functional modules according to needs, that is, the modules or steps in the foregoing embodiment of the present invention are further decomposed or combined, for example, the modules in the foregoing embodiment may be combined into one module, or may be further split into a plurality of sub-modules to complete all or part of the functions described above. The names of the modules and steps related to the embodiments of the present invention are merely for distinguishing the respective modules or steps, and are not to be construed as unduly limiting the present invention.
A storage device of a third embodiment of the present invention stores therein a plurality of programs adapted to be loaded by a processor and to implement a magnetic particle three-dimensional reconstruction imaging method based on generation of a countermeasure network as described above.
A processing device according to a fourth embodiment of the present invention includes a processor, a storage device; a processor adapted to execute each program; a storage device adapted to store a plurality of programs; the program is adapted to be loaded and executed by a processor for carrying out one of the above-mentioned magnetic particle three-dimensional reconstruction imaging methods based on generating a countermeasure network.
It will be clear to those skilled in the art that, for convenience and brevity of description, the specific working process of the storage device and the processing device and the related description of the foregoing description may refer to the corresponding process in the foregoing method example, which is not repeated herein.
Those of skill in the art will appreciate that the various illustrative modules, method steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both, and that the program(s) corresponding to the software modules, method steps, may be embodied in Random Access Memory (RAM), memory, read Only Memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art. To clearly illustrate this interchangeability of electronic hardware and software, various illustrative components and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as electronic hardware or software depends upon the particular application and design constraints imposed on the solution. Those skilled in the art may implement the described functionality using different approaches for each particular application, but such implementation is not intended to be limiting.
The terms "first," "second," "third," and the like, are used for distinguishing between similar objects and not for describing a particular sequential or chronological order.
Thus far, the technical solution of the present invention has been described in connection with the preferred embodiments shown in the drawings, but it is easily understood by those skilled in the art that the scope of protection of the present invention is not limited to these specific embodiments. Equivalent modifications and substitutions for related technical features may be made by those skilled in the art without departing from the principles of the present invention, and such modifications and substitutions will be within the scope of the present invention.

Claims (10)

1. A magnetic particle three-dimensional reconstruction imaging method based on generation of a countermeasure network, the method comprising the steps of:
step S100, acquiring a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged, and performing downsampling to obtain a downsampled sparse two-dimensional MPI image; the multi-view sparse two-dimensional MPI images are a plurality of MPI projection images which are collected in sequence according to set angles in a rotating mode;
step S200, inputting the downsampled sparse two-dimensional MPI image into a trained generation model of a self-attention generation countermeasure network by average diffusion search of neighborhood points, and obtaining a dense two-dimensional MPI image;
Step S300, reconstructing each dense two-dimensional MPI image by a filtering back projection reconstruction algorithm to finally obtain a three-dimensional MPI image of the object to be imaged;
the neighborhood point average diffusion search self-attention generation countermeasure network comprises a generation model and a discrimination model;
the generating model comprises five diffusion searching attention mechanism convolution modules; the diffusion search attention mechanism convolution module is constructed based on a neighborhood point average diffusion convolution sub-network, a self-attention fusion network and an activation function layer which are connected in sequence; the activation function adopted by the activation function layer in the first four diffusion searching attention mechanism convolution modules is a leakage ReLU, and the activation function adopted by the activation function layer in the fifth diffusion searching attention mechanism convolution module is a hyperbolic tangent function;
the neighborhood point average diffusion convolution sub-network is configured to sequentially encode the downsampled sparse two-dimensional MPI image through an encoding function to obtain a sequentially encoded image; performing regional convolution operation on each sequence coding image to obtain particle concentration information, angle information and association degree information among the images of each sequence coding image;
the self-attention fusion network is configured to multiply the particle concentration information, the angle information and the association degree information among the images of each sequence coding image with the corresponding weight matrix for linear transformation, fuse the images through a multi-head attention mechanism, and up-sample the fused images to obtain a dense two-dimensional MPI image;
The discrimination model comprises five convolutional neural network modules; the convolutional neural network module is constructed based on a convolutional layer, a regularization layer and an activation function layer which are sequentially connected; the activation function adopted by the activation function layer of the first four convolutional neural network modules is a leak ReLU, and the activation function adopted by the activation function layer of the fifth convolutional neural network module is a softmax.
2. The method for three-dimensional reconstruction imaging of magnetic particles based on generation of an countermeasure network according to claim 1, wherein the training process of generating the countermeasure network by average spread search self-attention of the neighborhood points is as follows:
step A100, acquiring a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged through MPI imaging equipment, and performing downsampling to obtain a downsampled sparse two-dimensional MPI image; obtaining a training data set based on the downsampled sparse two-dimensional MPI image and a corresponding truth value label thereof;
step A200, inputting the downsampled sparse two-dimensional MPI image into a pre-constructed neighborhood point average diffusion search self-attention generation countermeasure network generation model to obtain a dense two-dimensional MPI image;
step A300, inputting the dense two-dimensional MPI image and the truth value label corresponding to the dense two-dimensional MPI image into the neighborhood point average diffusion search self-attention generation countermeasure network discrimination model, and acquiring discrimination results of the dense two-dimensional MPI image;
Step A400, based on the discrimination result, combining each dense two-dimensional MPI image, each sparse two-dimensional MPI image and the corresponding truth value labels thereof, calculating a total loss value through a pre-constructed loss function, and updating the network parameters of the neighborhood point average diffusion search self-attention generation countermeasure network generation model and discrimination model;
and step A500, training the generator network and the discriminator network of the neighborhood point average spread search self-attention generation countermeasure network in a circulating way until the trained neighborhood point average spread search self-attention generation countermeasure network is obtained.
3. The method for three-dimensional reconstruction imaging of magnetic particles based on generation of countermeasure network according to claim 1, wherein the number of generation angles is set to be
Figure QLYQS_1
At an angle of +>
Figure QLYQS_2
At the time->
Figure QLYQS_3
Zhang Xishu two-dimensional MPI imageIs a function of the code of (a)
Figure QLYQS_4
The method comprises the following steps: />
Figure QLYQS_5
4. The method for three-dimensional reconstruction imaging of magnetic particles based on generation of countermeasure network according to claim 1, wherein the sequentially encoded images are subjected to a region convolution operation, which comprises the following steps: processing the concentration values of the pixel points in the sequential coding image, and performing operation by adopting a 3 multiplied by 3 convolution kernel after processing;
The method for processing the concentration value of the pixel point in the sequential coding image comprises the following steps:
Figure QLYQS_6
wherein ,
Figure QLYQS_8
representing the +.o in the sequentially encoded pictures>
Figure QLYQS_12
Line->
Figure QLYQS_15
Pixel points of column->
Figure QLYQS_9
Representing the neighborhood of points->
Figure QLYQS_10
Mean value of all values in>
Figure QLYQS_13
Representation->
Figure QLYQS_16
Minimum of all values in the point neighborhood, +.>
Figure QLYQS_7
Representation->
Figure QLYQS_11
Maximum value of all values in the point neighborhood, +.>
Figure QLYQS_14
Representation->
Figure QLYQS_17
Variable values in the point neighborhood.
5. The method for three-dimensional reconstruction imaging of magnetic particles based on generation of countermeasure network according to claim 1, wherein the particle concentration information, angle information, and correlation degree information between images of each sequentially encoded image are multiplied by corresponding weight matrices to perform linear transformation, the method is as follows:
Figure QLYQS_18
Figure QLYQS_19
wherein ,
Figure QLYQS_36
representing the self-attention mechanism, vector->
Figure QLYQS_40
Information indicating particle concentration>
Figure QLYQS_43
The method comprises the steps of carrying out a first treatment on the surface of the Vector->
Figure QLYQS_21
The information of the angle is represented by a set of angles,/>
Figure QLYQS_25
the method comprises the steps of carrying out a first treatment on the surface of the Vector->
Figure QLYQS_29
Information representing the degree of association between images, +.>
Figure QLYQS_34
,/>
Figure QLYQS_31
Representing the deformation of the matrix into a one-dimensional column vector, vector +.>
Figure QLYQS_35
Representing sequentially encoded images of the input neighborhood point average diffusion convolution sub-network; />
Figure QLYQS_38
、/>
Figure QLYQS_41
、/>
Figure QLYQS_37
Respectively representing weight matrices corresponding to particle concentration information, angle information and correlation degree information between images of sequentially encoded images, ++>
Figure QLYQS_39
Representing each sequentially encoded image +. >
Figure QLYQS_42
Neurons in all output layers +.>
Figure QLYQS_44
Sum of output results from 1 to N, < ->
Figure QLYQS_23
Representing each sequentially encoded image +.>
Figure QLYQS_26
Through->
Figure QLYQS_28
Output values of individual neurons, < >>
Figure QLYQS_32
Representing each sequentially encoded image +.>
Figure QLYQS_20
Through->
Figure QLYQS_27
Output values of individual neurons, < >>
Figure QLYQS_30
Representing each sequentially encoded image +.>
Figure QLYQS_33
Through->
Figure QLYQS_22
The output value of the individual neurons is the total output value +.>
Figure QLYQS_24
Is a specific gravity of (c).
6. The method for three-dimensional reconstruction imaging of magnetic particles based on generation of countermeasure network according to claim 5, wherein the particle concentration information, the angle information and the association degree information between images of each sequentially encoded image are fused by a multi-head attention mechanism, the method comprising:
Figure QLYQS_45
Figure QLYQS_46
wherein ,
Figure QLYQS_47
is a splicing function->
Figure QLYQS_48
Indicate->
Figure QLYQS_49
Output of the attention machine submodule, +.>
Figure QLYQS_50
Total weight matrix representing the output of the multi-head attention mechanism,/->
Figure QLYQS_51
And information obtained by fusing the particle concentration information, the angle information and the association degree information between the images of the sequential encoded images.
7. A method of three-dimensional reconstruction imaging of magnetic particles based on generation of a countermeasure network according to claim 2, characterized in that the total loss value is calculated by a pre-constructed loss function by:
Figure QLYQS_52
Figure QLYQS_53
Figure QLYQS_54
Figure QLYQS_55
wherein ,
Figure QLYQS_59
for the total loss function- >
Figure QLYQS_63
Representing pixel loss function and contrast loss function, < ->
Figure QLYQS_67
Representing an angular sequential ordering optimization function,/->
Figure QLYQS_58
Representing pixel information error,/>
Figure QLYQS_60
Representing the loss of the generative model->
Figure QLYQS_65
A sparse two-dimensional MPI image is represented,
Figure QLYQS_69
representing the generated->
Figure QLYQS_57
Zhang Miji two-dimensional MPI image,>
Figure QLYQS_62
representing the generated->
Figure QLYQS_64
Zhang Miji two-dimensional MPI image,>
Figure QLYQS_68
indicate->
Figure QLYQS_71
Zhang Xishu truth label corresponding to two-dimensional MPI image, < ->
Figure QLYQS_74
Coefficients representing pixel loss, < >>
Figure QLYQS_77
Coefficient representing countermeasures against losses->
Figure QLYQS_80
Represents the probability of judging whether or not it is true by the judging network, < >>
Figure QLYQS_72
Representing the output of the generative model,/>
Figure QLYQS_75
Representing the angle difference function, ++>
Figure QLYQS_78
As hyperbolic tangent function, +.>
Figure QLYQS_81
Representing the predicted value of the angle at the corresponding input,
Figure QLYQS_56
,/>
Figure QLYQS_61
representing any two dense two-dimensional MPI images +.>
Figure QLYQS_66
Angle judging function (F)>
Figure QLYQS_70
and />
Figure QLYQS_73
Respectively represent +.>
Figure QLYQS_76
Angle label corresponding to dense two-dimensional MPI image generated by sheet,/->
Figure QLYQS_79
Parameters indicating the richness settings for the information contained in the different mimics.
8. A magnetic particle three-dimensional reconstruction imaging system based on a generation countermeasure network, the system comprising: an image acquisition module 100, an image reconstruction module 200 and a three-dimensional reconstruction module 300;
the image acquisition module 100 is configured to acquire a multi-view sparse two-dimensional MPI image of an imitation body of a reconstructed object to be imaged, and perform downsampling to obtain a downsampled sparse two-dimensional MPI image; the multi-view sparse two-dimensional MPI images are a plurality of MPI projection images which are collected in sequence according to set angles in a rotating mode;
The image reconstruction module 200 is configured to input the downsampled sparse two-dimensional MPI image into a trained generation model of a neighborhood point average diffusion search self-attention generation countermeasure network to obtain a dense two-dimensional MPI image;
the three-dimensional reconstruction module 300 is configured to reconstruct each dense two-dimensional MPI image through a filtered back projection reconstruction algorithm, and finally obtain a three-dimensional MPI image of the reconstructed object to be imaged;
the neighborhood point average diffusion search self-attention generation countermeasure network comprises a generation model and a discrimination model;
the generating model comprises five diffusion searching attention mechanism convolution modules; the diffusion search attention mechanism convolution module is constructed based on a neighborhood point average diffusion convolution sub-network, a self-attention fusion network and an activation function layer which are connected in sequence; the activation function adopted by the activation function layer in the first four diffusion searching attention mechanism convolution modules is a leakage ReLU, and the activation function adopted by the activation function layer in the fifth diffusion searching attention mechanism convolution module is a hyperbolic tangent function;
the neighborhood point average diffusion convolution sub-network is configured to sequentially encode the downsampled sparse two-dimensional MPI image through an encoding function to obtain a sequentially encoded image; performing regional convolution operation on each sequence coding image to obtain particle concentration information, angle information and association degree information among the images of each sequence coding image;
The self-attention fusion network is configured to multiply the particle concentration information, the angle information and the association degree information among the images of each sequence coding image with the corresponding weight matrix for linear transformation, fuse the images through a multi-head attention mechanism, and up-sample the fused images to obtain a dense two-dimensional MPI image;
the discrimination model comprises five convolutional neural network modules; the convolutional neural network module is constructed based on a convolutional layer, a regularization layer and an activation function layer which are sequentially connected; the activation function adopted by the activation function layer of the first four convolutional neural network modules is a leak ReLU, and the activation function adopted by the activation function layer of the fifth convolutional neural network module is a softmax.
9. A storage device in which a plurality of programs are stored, characterized in that the programs are adapted to be loaded and executed by a processor to implement a method of generating a countermeasure network-based magnetic particle three-dimensional reconstruction imaging as claimed in any of claims 1 to 7.
10. A processing device, comprising a processor and a storage device; a processor adapted to execute each program; a storage device adapted to store a plurality of programs; characterized in that the program is adapted to be loaded and executed by a processor for implementing a method for three-dimensional reconstruction imaging of magnetic particles based on generation of a countermeasure network as claimed in any of claims 1 to 7.
CN202310047876.0A 2023-01-31 2023-01-31 Magnetic particle three-dimensional reconstruction imaging method based on generation countermeasure network Active CN115880440B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310047876.0A CN115880440B (en) 2023-01-31 2023-01-31 Magnetic particle three-dimensional reconstruction imaging method based on generation countermeasure network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310047876.0A CN115880440B (en) 2023-01-31 2023-01-31 Magnetic particle three-dimensional reconstruction imaging method based on generation countermeasure network

Publications (2)

Publication Number Publication Date
CN115880440A CN115880440A (en) 2023-03-31
CN115880440B true CN115880440B (en) 2023-04-28

Family

ID=85758588

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310047876.0A Active CN115880440B (en) 2023-01-31 2023-01-31 Magnetic particle three-dimensional reconstruction imaging method based on generation countermeasure network

Country Status (1)

Country Link
CN (1) CN115880440B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116563412B (en) * 2023-06-26 2023-10-20 中国科学院自动化研究所 MPI image reconstruction method, system and equipment based on sparse system matrix
CN116797739B (en) * 2023-07-07 2024-09-06 深圳市魔数智擎人工智能有限公司 Real tree interactive modeling method based on sparse image

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021017372A1 (en) * 2019-08-01 2021-02-04 中国科学院深圳先进技术研究院 Medical image segmentation method and system based on generative adversarial network, and electronic equipment
CN113808234A (en) * 2021-11-08 2021-12-17 北京航空航天大学 Rapid magnetic particle imaging reconstruction method based on undersampling
CN113850883A (en) * 2021-10-14 2021-12-28 北京航空航天大学 Magnetic particle imaging reconstruction method based on attention mechanism
CN114792287A (en) * 2022-03-25 2022-07-26 南京航空航天大学 Medical ultrasonic image super-resolution reconstruction method based on multi-image fusion
CN115409945A (en) * 2022-09-01 2022-11-29 中国科学院自动化研究所 Three-dimensional magnetic particle imaging system and method for carrying out quantitative analysis by fusing imaging parameters
CN115526946A (en) * 2022-10-13 2022-12-27 中国科学院自动化研究所 Method, system and equipment for denoising magnetic particle imaging image based on feature fusion

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110420026B (en) * 2019-07-15 2020-05-19 中国科学院自动化研究所 Magnetic particle imaging three-dimensional reconstruction method, system and device based on FFL

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021017372A1 (en) * 2019-08-01 2021-02-04 中国科学院深圳先进技术研究院 Medical image segmentation method and system based on generative adversarial network, and electronic equipment
CN113850883A (en) * 2021-10-14 2021-12-28 北京航空航天大学 Magnetic particle imaging reconstruction method based on attention mechanism
CN113808234A (en) * 2021-11-08 2021-12-17 北京航空航天大学 Rapid magnetic particle imaging reconstruction method based on undersampling
CN114792287A (en) * 2022-03-25 2022-07-26 南京航空航天大学 Medical ultrasonic image super-resolution reconstruction method based on multi-image fusion
CN115409945A (en) * 2022-09-01 2022-11-29 中国科学院自动化研究所 Three-dimensional magnetic particle imaging system and method for carrying out quantitative analysis by fusing imaging parameters
CN115526946A (en) * 2022-10-13 2022-12-27 中国科学院自动化研究所 Method, system and equipment for denoising magnetic particle imaging image based on feature fusion

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Jiwon Kim 等.Accurate Image Super-Resolution Using Very Deep Convolutional Networks.IEEE:COMPUTE` .\R SOCIETY.2016,正文第1646-1654页. *
基于改进超分辨率生成对抗网络的图像重建算法;查体博 等;基于改进超分辨率生成对抗网络的图像重建算法;第58卷(第8期);08100050-1-11 *

Also Published As

Publication number Publication date
CN115880440A (en) 2023-03-31

Similar Documents

Publication Publication Date Title
CN115880440B (en) Magnetic particle three-dimensional reconstruction imaging method based on generation countermeasure network
Shen et al. Patient-specific reconstruction of volumetric computed tomography images from a single projection view via deep learning
EP3506209B1 (en) Image processing method, image processing device and storage medium
CN112085829A (en) Spiral CT image reconstruction method and equipment based on neural network and storage medium
Nguyen et al. Evaluating transformer-based semantic segmentation networks for pathological image segmentation
Liu et al. Isotropic reconstruction of electron tomograms with deep learning
CN112884819B (en) Image registration and neural network training method, device and equipment
Lee et al. High-throughput computer method for 3D neuronal structure reconstruction from the image stack of the Drosophila brain and its applications
Tang et al. Point scene understanding via disentangled instance mesh reconstruction
CN113628296B (en) Magnetic particle imaging reconstruction method from time-frequency domain signal to two-dimensional image
CN114170146A (en) Image processing method, image processing device, electronic equipment and computer readable storage medium
CN117152441B (en) Biological image instance segmentation method based on cross-scale decoding
Powell et al. Learning structural heterogeneity from cryo-electron sub-tomograms with tomoDRGN
CN114913262B (en) Nuclear magnetic resonance imaging method and system with combined optimization of sampling mode and reconstruction algorithm
Kim et al. A convolutional neural network-based anthropomorphic model observer for signal detection in breast CT images without human-labeled data
Zheng et al. Identifying sets of favorable projections for few-view low-dose cone-beam CT scanning
Venkatesan et al. Nodule detection with convolutional neural network using apache spark and GPU frameworks
CN112233227B (en) CT projection graph generation method and device
Lou et al. Predicting radiologists' gaze with computational saliency models in mammogram reading
Yan et al. Vaccine-style-net: Point cloud completion in implicit continuous function space
CN116188615A (en) Sparse angle CT reconstruction method based on sine domain and image domain
Guo et al. Noise-resilient deep learning for integrated circuit tomography
Guo et al. Iterative Image Reconstruction for Limited‐Angle CT Using Optimized Initial Image
US20230360366A1 (en) Visual Explanation of Classification
Borrelli et al. Deep learning for accelerating Radon inversion in single-cells tomographic phase imaging flow cytometry

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant