The present invention belongs to the field of PET (positron emission tomography) imaging, which involves a dual-tracer PET separation method based on Deep Belief Networks (DBN) aiming at the tracer pairs labelled with the same isotope.
Positron emission tomography (PET) is a kind of functional medical imaging modality used to detect the physiochemical activity in living bodies. The technology is often used to measure the physiological indexes (i.e. glucose metabolism, hypoxia, blood flow, and hypoxia) of the certain parts. The technology of imaging different functions in different parts in the living body mainly relies on various radioactive tracers. Researchers replaced one or more atoms in a chemical compound (i.e. Glucose, protein, nucleic acid) by a radionuclide (^{11}C, ^{13}N, ^{15}O, ^{18}F is most common) to make radioactive tracers. By using bolus and continuous infusion protocols, the tracers will concentrate in a certain part in the living body that referred as Regions of Interest (ROIs), where it decays in a certain rate. During the decay period, it emits a positron, which travels a short distance before annihilating with an electron. This annihilation produces two high-energy (511 keV) γ photons propagating in nearly opposite directions. The two γ photons will be detected by detectors outside of the body, then the reconstruction algorithm will be performed to recover a radioactive concentration distribution image.
However, for a specific part of the body, multi-faceted, multi-angle detection can provide more information, and depicting the physiological and functional status of the tumor from different aspects also helps to improve the diagnostic accuracy, so the dual-tracer PET imaging is very necessary. In order to save scan time, reduce patient suffering and improve diagnostic efficiency, single-scan dual-tracer PET dynamic imaging has become a key technology to be solved. Since an annihilation reaction produces a pair of 511 ekV γ photons, it cannot be from the energy point of view to distinguish between the two tracers. At present, there are two main types of mainstream dual-tracer PET separation methods: (1) separation of reconstructed dual-tracer PET images; (2) integration of separation algorithm into the separation process of dual-tracer PET, and direct reconstruction of the concentration distribution of the two radiotracers. The first type of algorithm is essentially a signal separation problem, and the solution is relatively extensive. The second type of algorithm requires complex reconstruction algorithms to support, and the practicality is not strong. At present, the researchers mainly focus on the first type of algorithm.
The separation operation based on the mixed dual-tracer PET image mainly depends on the difference of the tracer half-life and kinetic parameters, combined with the compartment model to solve, and most of them need to be interlaced. Two tracers were injected to provide partial non-mixed single tracer information to facilitate modeling and separation, which lengthens the scanning process for the entire PET.
The dual tracers labeled with the same radioactive isotope play important roles for clinical use. For example, [^{18}F]FDG-[^{18}F]FLT is used for detecting cell proliferation, [^{11}C]FMZ-[^{11}C]DTBZ is used for measuring neurotransmitter system, [^{62}Cu]ATSM-[^{62}Cu]PTSM is used for observing blood flow and hypoxia. However, when dual tracers are labeled with the same isotope, the distinguishability of the two signals is diminished, and clinical requirements require simultaneous injection of two tracers to reduce scan time, which renders most algorithms ineffective. The separation of PET signals with the same isotope labeling has become a difficult problem to solve.
According to all described above, the present invention provides a DBN based separation method of a mixture of dual-tracer single-acquisition PET signals labelled with the same isotope. With the help of information exaction ability of deep learning, two individual PET volumetric images can be separated from a mixture of dual-tracer PET images.
A DBN based separation method of a mixture of dual-tracer single-acquisition PET signals labelled with the same isotope comprises the following steps:
Furthermore, the method used to extract TACs from dynamic PET image sequences X^{Dual}, X^{I }and X^{II }in step (5) can be formulated as:
X
^{dual}=[TAC_{1}^{dual }TAC_{2}^{dual }. . . TAC_{n}^{dual}]^{T }
X
^{I}=[TAC_{1}^{I }TAC_{2}^{I }. . . TAC_{n}^{I}]^{T }
X
^{II}=[TAC_{1}^{II }TAC_{2}^{II }. . . TAC_{n}^{II}]^{T }
wherein TAC_{1}^{dual}˜TAC_{n}^{dual }are the TACs of the 1st-n-th pixels in the dynamic PET images sequence X^{Dual}, TAC_{1}^{I}˜TAC_{n}^{I }are the TACs of the 1st-n-th pixels in the dynamic PET images sequence X^{I}, TAC_{1}^{II}˜TAC_{n}^{II }are the TACs of the 1st-n-th pixels in the dynamic PET images sequence X^{II}, n represents the total pixels of the PET image, and T represents transposition.
Furthermore, the training process of DBN in the described step (5) can be detailed as follows:
5.1 initializing a DBN framework consisting of an input layer, hidden layers and an output layer, wherein the hidden layers are composed of three stacked Restricted Boltzmann Machines (RBMs);
5.2 initializing the parameters in the described DBN, which includes the number of nodes in hidden layers, the offset vector and the weight vector between layers, learning rate, activation function and the maximum number of iterations;
5.3 pre-training the stacked RBMs in the hidden layers;
5.4 transmitting the parameters obtained from pre-training to initialized DBN, and then inputting the TAC in the dynamic PET image sequence X^{dual }one by one into the DBN, calculating the error function L between the output result and the corresponding ground truth and applying the gradient descent method to continuously update the parameters of the whole network until the error function L converges or the maximum number of iteration is reached then a trained PET reconstruction model can be acquired.
Furthermore, in step 5.3, the restricted Boltzmann machine in the hidden layer is pre-trained, that is, each restricted Boltzmann machine is composed of a display layer and a hidden layer, and by the contrast divergence algorithm, the weight of each layer will be updated until the hidden layer can accurately represent the features of the display layer and can reverse the display layer. Furthermore, the loss function L described in step 5.4 is as follows:
L=∥−TAC_{j}^{I}∥_{2}^{2}+∥−TAC_{j}^{II}∥_{2}^{2}+ζ∥(+)−(+∥_{2}^{2 }
wherein TAC_{j}^{I }is the TAC of the j-th pixel in the dynamic PET image sequence X^{I}; TAC_{j}^{II }is the TAC of the j-th pixel in the dynamic PET image sequence X^{II}; TAC_{j}^{I }and TAC_{j}^{II }represent the two output results correspond to those in X^{I }and X^{II }that obtained after the TAX of the j-th pixel in the dynamic PET image sequence X^{dual }are substituted into the DBN; j is a natural number and 1≤j≤n, n is the total number of the PET images, ∥ ∥_{2 }the L2 norm; ζ is a user defined constant.
Furthermore, in the step (6), the TAC of the j-th pixel of X^{dual }in the test set is input into the dual-tracer PET reconstruction model, and the output of the pixel is [
X
_{test}
^{I}=[
X
_{test}
^{II}=[
The present invention described above can achieve the separation of a mixture of dual-tracer single-acquisition PET signals labelled with the same isotope by the trained DBN framework. A point-to-point mapping relationship between a mixed dual-tracer PET image and two single-tracer PET images was learned by inputting training data and ground truth into the constructed neural network. It is worth mentioning that this network is universal. The training set contains multiple sets of tracer combinations labeled with the same isotope, and the results show the good performance of this model on dual-tracer separation.
In summary, the present invention includes a universal framework constructed by deep belief network to establish a mapping relationship for the separation of dual tracers. By the powerful feature extraction capabilities of it, the signal separation of dual-tracer PET with the same isotope can be achieved well
In order to more specifically describe the present invention, the detailed instructions are provided in conjunction with the attached figures and following specific embodiments:
As
(1) Preparation of Training Data
Extracting pixel based time activity curves (TACs) from the reconstructed images X^{dual}, Y^{I}, Y^{II}, the process can be described as follows:
X
^{dual}=[x_{1},x_{2},x_{3}, . . . ,x_{N}]^{T},x_{i}=[x_{i}^{1},x_{i}^{2}, . . . ,x_{i}^{M}]^{T }
Y
^{I}=[(y^{1})_{1},(y^{1})_{2}, . . . ,(y^{1})_{N}]^{T},(y^{1})_{i}=[(y^{1})_{i}^{1},(y^{1})_{i}^{2}, . . . ,(y^{1})_{i}^{M}]^{T }
Y
^{II}=[(y^{2})_{1},(y^{2})_{2}, . . . ,(y^{2})_{N}]^{T},(y^{2})_{i}=[(y^{2})_{i}^{1},(y^{2})_{i}^{2}, . . . ,(y^{2})_{i}^{M}]^{T }
wherein x_{i}^{j }represents the radioactive concentration of pixel i at j-th frame, (y^{1})_{i}^{j }and (y^{2})_{i}^{j }represent the radiotracer concentration of pixel point i of Tracer I and Tracer II at j-th frame, respectively, and N is the total number of pixels of the resulting PET image, M is the total number of frames acquired by dynamic PET.
(2) Preparation of Training Sets and Test Set Data.
70% of the TAC data set X^{dual }is extracted as the training set X_{train}^{dual}, and the remaining 30% is used as the testing set X_{test}^{dual}, and Y^{I }and Y^{II }corresponding to the training set and the test set are respectively connected in series as the labels of the training set and the ground truth of the testing set, which can be shown as follows:
(3) Constructing a deep belief network for the signal separation of the dual-tracer PET with the same isotope; as shown in
(4) The training set is input into this network for training. The training process is as follows:
The error function L is as follows:
L=∥
−Y
_{train}
^{I}∥_{2}^{2}+∥−Y_{train}^{II}∥_{2}^{2}−γ(∥−Y_{train}^{I}∥_{2}^{2}+∥−Y_{train}^{II}∥_{2}^{2})
wherein, the first two reflect the error between the ground truth and the predicted value, while the latter two are the difference between the two tracer signals; γ is a defined constant used to adjust the proportion of the latter two items in the loss function.
Loss(_{j}^{I},_{j}^{II})=∥_{j}^{I}−TAC_{j}^{I}∥_{2}^{2}+∥_{j}^{II}−TAC_{j}^{II}∥_{2}^{2}+ξ∥_{j}^{I}+_{j}^{II})−(TAC_{j}^{I}+TAC_{j}^{II})∥_{2}^{2 }
wherein _{j}^{I }and _{j}^{II }and are the prediction value of TAC_{j}^{I }and TAC_{j}^{II }by the DBN respectively, and ξ is the weigh coefficient.
(5) Inputting the TACs of the test set data X_{test}^{dual }into the trained neural network to obtain the separation signal of the dual-tracer PET labeled with the same isotope.
Next, we validate the present invention by simulated experiment.
(1) Phantom Selection
There are three different dual tracers in training datasets, and every group is equipped with a different phantom with different Regions of interest (ROIs), which indicates different biochemical environment.
(2) The Simulation of PET Concentration Distribution
The parallel compartment model was used to simulate the motion of dual tracers, and then the stable dual-tracer concentration distribution can be acquired by solving the dynamic ordinary differential equations (ODEs). The single compartment model based on kinetic parameters was used to simulate the motion of a single tracer in vivo, and then the stable single-tracer concentration distribution can be acquired by solving the dynamic ordinary differential equations (ODEs).
(3) The Simulation of PET Scanning Process
Computer Monte Carlo simulations were used to perform the dynamic dual-tracer PET scans with the help of software GATE. All simulations are based on the geometry of the full 3D whole body PET scanner SHR74000, designed by Nippon Kanematsu Photonics Co., Ltd. The PET scanner is designed with 6 rings and each ring has 48 detector blocks. Each detector block consists of a 16×16 array of lutetium yttrium orthosilicate (LYSO) crystals. The ring diameter of the scanner is 826 mm. When three groups of dual tracers and a single tracer concentration distribution map was input the Monte Carlo system, corresponding sinogram data were acquired.
(4) Reconstruction Process
The sinogram was reconstructed using the classical ML-EM reconstruction algorithm to obtain the concentration distribution of the simulated radiotracer pairs.
(5) Acquisition of TAC Curve;
TACs Based on pixel were obtained by recombining the concentration distribution matrix of the three groups of mixed radioactive concentration tracers in the body.
(6) Training Process;
the 70% TACs of three dual-tracer groups ([^{8}F] FDG+[^{18}F] FLT, [^{11}C]FMZ+[^{11}C] DTBZ and [^{62}Cu] ATSM+[^{62}Cu] PTSM) was inputted to the DBN as training data to do the pre-training. The TAC curve of a single tracer serves as a label to provide feedback for fine-tuning the entire network.
(7) Testing
The remaining 30% was used to evaluate the validity of the network.
Comparing the predicted image with the simulated ground truth, it can be found that the constructed deep belief network can separate the dual-tracer PET signal with the same isotope label well. This confirms the effectiveness of the deep belief network in feature extraction and signal separation, and also demonstrates that the method of the invention is effective in processing PET signals labeled with the same isotope.
The description of the specific instances is to facilitate ordinary technicians in the technical field to understand and apply the present invention. It is obvious that a person familiar with the technology in this field can easily modify the specific implementation mentioned above and apply the general principles described here into other instances without the creative Labor. Therefore, the present invention is not limited to the above instances. According to the disclosure of the present invention, the improvement and modification of the present invention will be all within the protection scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
201810507876.3 | Aug 2018 | CN | national |