WO2021062885A1 - Ct imaging method and apparatus, storage medium, and medical imaging system - Google Patents
Ct imaging method and apparatus, storage medium, and medical imaging system Download PDFInfo
- Publication number
- WO2021062885A1 WO2021062885A1 PCT/CN2019/111028 CN2019111028W WO2021062885A1 WO 2021062885 A1 WO2021062885 A1 WO 2021062885A1 CN 2019111028 W CN2019111028 W CN 2019111028W WO 2021062885 A1 WO2021062885 A1 WO 2021062885A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- projection data
- sub
- network module
- image
- noise
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 79
- 238000002059 diagnostic imaging Methods 0.000 title claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 45
- 238000013170 computed tomography imaging Methods 0.000 claims abstract description 26
- 230000009467 reduction Effects 0.000 claims abstract description 21
- 238000012545 processing Methods 0.000 claims description 28
- 230000006870 function Effects 0.000 claims description 27
- 238000012549 training Methods 0.000 claims description 14
- 230000004913 activation Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 7
- 238000002591 computed tomography Methods 0.000 claims description 4
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000003325 tomography Methods 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 238000004891 communication Methods 0.000 description 9
- 238000001514 detection method Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 229910052704 radon Inorganic materials 0.000 description 2
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013501 data transformation Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000010977 jade Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003607 modifier Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
Definitions
- the embodiments of the present application relate to medical imaging technology, for example, to a computer tomography (Computed Tomography, CT) imaging method, device, storage medium, and medical imaging system.
- a computer tomography Computer tomography (Computed Tomography, CT) imaging method
- device storage medium
- medical imaging system for example, to a computer tomography (Computed Tomography, CT) imaging method, device, storage medium, and medical imaging system.
- CT Computer tomography
- X-ray-based CT imaging is a common way to assist diagnosis.
- X-ray imaging is radioactive and will increase the risk of cancer in the subject.
- low-dose CT detection has become one of the research directions in the field of CT imaging.
- the low-dose CT reconstruction algorithm is developed based on iterative reconstruction technology, which simulates CT image reconstruction as a mathematical optimization problem.
- the iterative reconstruction algorithm simulates a comprehensive projection on the estimated image through forward projection, and the process simulates the reality as much as possible.
- the integrated projection is compared with the real measurement value collected by the detector, and the next update is determined based on the difference between the two. Correct the current estimated image.
- the present application provides a CT imaging method, device, storage medium, and medical imaging system to improve the efficiency and quality of CT image reconstruction.
- the embodiment of the present application provides a CT imaging method, including:
- the projection data is sent to a pre-trained imaging model, and the output image of the imaging model is determined as the CT image of the target object, wherein the imaging model is composed of a standard dose CT image and sample projection data containing noise
- the imaging model includes a first sub-network module, a second sub-network module, and a third sub-network module.
- the first sub-network module is used to perform noise reduction processing on projection data in the projection domain
- the second The sub-network module is used to transform the processed projection data to generate CT domain data
- the third sub-network module is used to perform noise reduction processing on the CT domain data to generate a CT image of the target object.
- An embodiment of the present application also provides a CT imaging device, which includes:
- a projection data acquisition module configured to acquire projection data of a target object, wherein the projection data is acquired when the target object is detected based on a preset dose of X-ray photons, and the preset dose is less than the CT standard dose;
- the CT image reconstruction module is configured to send the projection data to a pre-trained imaging model, and determine the output image of the imaging model as the CT image of the target object, wherein the imaging model consists of a standard dose of CT image And the sample projection data containing noise is obtained by training, the imaging model includes a first sub-network module, a second sub-network module and a third sub-network module, the first sub-network module is used to reduce the projection data of the projection domain Noise processing, the second sub-network module is used to transform the processed projection data to generate CT domain data, and the third sub-network module is used to perform noise reduction processing on the CT domain data to generate the target CT image of the subject.
- the embodiment of the present application also provides a computer-readable storage medium that stores a computer program, and when the program is executed by a processor, the CT imaging method as provided in any embodiment of the present application is implemented.
- the embodiment of the present application also provides a medical imaging system, which includes a medical imaging device and a computer device.
- the computer device includes a memory, one or more processors, and a computer program stored in the memory and running on the processor, When the program is executed by the processor, the CT imaging method as provided in any embodiment of the present application is realized.
- FIG. 1 is a schematic flowchart of a CT imaging method provided in Embodiment 1 of the present application;
- FIG. 2 is a schematic structural diagram of an imaging model provided in Embodiment 1 of the present application.
- FIG. 3 is a schematic structural diagram of another imaging model provided in Embodiment 1 of the present application.
- FIG. 4 is a schematic flowchart of another CT imaging method provided in Embodiment 2 of the present application.
- FIG. 5 is a schematic structural diagram of a CT imaging apparatus according to Embodiment 3 of the present application.
- FIG. 6 is a schematic structural diagram of a medical imaging system provided by Embodiment 4 of the present application.
- Fig. 1 is a schematic flow chart of a CT imaging method provided in the first embodiment of this application. This embodiment can be applied to the situation of rapidly generating high-quality low-dose CT images based on an imaging model.
- the CT imaging method includes S110-S120.
- the imaging model is obtained by training a standard dose CT image and sample projection data containing noise.
- the imaging model includes a first sub-network module, a second sub-network module, and a third sub-network module.
- the network module is used to denoise the projection data in the projection domain
- the second sub-network module is used to transform the processed projection data to generate CT domain data
- the third sub-network module is used to The CT domain data undergoes noise reduction processing to generate a CT image of the target object.
- the target object may be a human or an animal, and the detection of the target object by low-dose X-ray photons can reduce radiation damage to the target object by X-rays during the detection process, and reduce the risk of cancer of the target object.
- the imaging model is pre-trained and has the function of reconstructing low-dose projection data and the function of reducing noise in the low-dose projection data.
- the target object will be detected by low-dose X-ray photons.
- the collected projection data is input to the imaging model to generate a high-quality CT image of the target object after noise reduction. Since the imaging model is pre-trained, it can automatically process the input projection data and output the CT image of the target object.
- the imaging model performs noise reduction processing on the projection data, which improves the clarity of the image generated by low-dose detection.
- the imaging model includes a first sub-network module, a second sub-network module, and a third sub-network module.
- FIG. 2 is a schematic structural diagram of an imaging model provided by an embodiment of the present application.
- the first sub-network module is used to perform convolution filtering processing on the input projection image, and the projection domain is a domain formed by all the projection images.
- the first sub-network module may include a first preset number of convolution modules connected in sequence. Each convolution module in the first preset number of convolution modules includes a convolution layer and an activation function layer.
- the network module can perform noise reduction and filtering processing on the projection data. In an embodiment, the first preset number may be 5, 6, 10, and so on.
- the convolution kernel of the convolution layer in the first sub-network module is a ⁇ b, a and b are both positive integers greater than or equal to 1, and b>a.
- the convolution kernel of the convolution layer in the first sub-network module may be 1 ⁇ 30, 3 ⁇ 30, 5 ⁇ 30, 3 ⁇ 33, 5 ⁇ 33, etc., by adopting the convolution of the form b>a
- the verification can simultaneously achieve the noise reduction and filtering functions of the projection data.
- the activation function layer may be a Leaky Rectified Linear Unit (Leaky ReLU).
- the second subconvolution model is used to convert projection data in the projection domain into CT domain data, which is a domain formed by all CT images.
- the second sub-network module transforms the processed projection data according to the following formula:
- f(x, y) is the CT domain image output by the second sub-network module
- x and y are the abscissa and ordinate in the CT domain image, respectively
- p(r, ⁇ ) is the second input
- the projection data of the sub-network module r is the distance between the projection data and the origin
- ⁇ is the projection X-ray
- ⁇ is the projection angle.
- the second sub-network module reconstructs the projection data input by the first sub-network module according to the above formula, and outputs a CT domain image of the same size as the expected size, and the CT domain image satisfies the CT image attributes.
- the third sub-network module is used to perform convolution filtering operations on the CT domain image, and simultaneously perform image restoration processing and noise reduction processing on the obtained CT domain image.
- restoring the image may include a second preset number of convolution modules connected in sequence, and each convolution module of the second preset number of convolution modules includes a convolution layer and an activation function layer.
- the second preset number may be the same as or different from the first preset number, the second preset number may be, for example, 5, 6, or 10, etc., and the second preset number and the first preset number may be Respectively determined according to the functional accuracy.
- the activation function layer may be Leaky ReLU
- the convolution kernel of the convolution layer in the third sub-network module is m ⁇ m
- the m is a positive integer greater than or equal to 1
- the convolution kernel of the convolution layer in can be 3 ⁇ 3, 5 ⁇ 5, or 7 ⁇ 7.
- FIG. 3 is a schematic structural diagram of another imaging model provided in Embodiment 1 of the present application.
- the first sub-network module in Figure 3 includes 6 convolution modules.
- the convolution kernel of the convolution layer in the convolution module is 3 ⁇ 30 with a step size of 1.
- the third sub-network module includes 6 convolution modules.
- the convolution kernel of the convolution layer in the product module is 3 ⁇ 3, and the step size is 1.
- the size of the input projection data in Figure 3 is 900 ⁇ 848, and the size of the output CT image is 512 ⁇ 512.
- the above image size is only an example of Figure 3.
- the size of the input projection data can be determined according to the input layer of the imaging model. Previously, if the size of the collected projection data was inconsistent with the processing size of the input layer of the imaging module, the collected projection data was adjusted in advance to make the collected projection data meet the processing size of the input layer of the imaging module.
- the input end of the first sub-network module is short-circuited with the output end of the first sub-network module, and the input end of the third sub-network module is connected to the third sub-network module
- the output terminal is shorted.
- the projection domain data and the CT domain data are respectively denoised, so as to realize double denoising in the image reconstruction process, improve the denoising effect, and make When X-ray photons are detected, high-quality CT images are obtained.
- the three sub-network modules are connected to each other and are independent of each other. Any module can be updated and replaced according to imaging requirements, which improves the configurability and applicability of the imaging model.
- Fig. 4 is a schematic flow chart of a CT imaging method provided in the second embodiment of the present application.
- a training method of an imaging model is provided.
- the CT imaging method includes S210-S250.
- S220 Input the sample projection data containing noise into the initial imaging model to obtain a reconstructed image.
- S230 Determine a loss function according to the reconstructed image and the CT image of the standard dose corresponding to the noise-containing sample projection data, and adjust the network parameters in the initial imaging model according to the loss function to generate the imaging model.
- S250 Send the projection data to a pre-trained imaging model, and determine an output image of the imaging model as a CT image of the target object.
- the initial imaging model is trained based on sample data to obtain an imaging model with image reconstruction function and denoising function.
- the sample data includes low-dose projection data containing noise and CT images of standard doses.
- the low-dose projection data containing noise can be collected when the sample object is detected by low-dose X-ray photons, or can be obtained by adding noise based on the standard-dose CT image.
- acquiring the noise-containing sample projection data may be: acquiring a standard dose CT image, and preprocessing the standard dose CT image, Obtain standard projection data; add noise data to the standard projection data to generate sample projection data containing noise.
- the preprocessing of the standard dose CT image may be to perform Radon operation on the standard dose CT image, generate standard projection data, and add noise to the standard projection data, such as Poisson noise, to obtain noise-containing Sample projection data.
- the standard-dose CT image may be a CT image obtained by clinical testing.
- the sample projection data is generated by adding noise to the standard projection data, so as to avoid the X-ray scanning of the sample object during the sample generation process, which will cause the sample object to be damaged. Risk of cancer.
- different basic noises are added to the standard projection data, different sample projection data are generated, the variety of noise in the sample is increased, and the number of samples is increased.
- the adding noise to the standard projection data to generate sample projection data containing noise includes: setting a Poisson function of at least two noise levels; and a pair of Poisson functions based on the at least two noise levels Noise is added to the standard projection data, and at least two noise-containing sample projection data corresponding to the CT image of the standard dose are generated.
- the orthographic projection data of the CT image of the standard dose namely the sine image (sino image)
- sino_ps is the sample projection data containing noise
- del_num is the number of detectors
- angle_num is the total number of projection angles
- I 1 I 0 ⁇ k
- k noise Intensity coefficient
- I 0 is the basic photon number
- I 0 1 ⁇ 10 6 .
- k ⁇ 1 can be 1, 0.1, 0.2, 0.5, or 0.05, etc.
- the value of k is not limited, and the added noise level is adjusted by the value of k. The smaller the k, the stronger the noise.
- sample projection data with different levels of noise are formed. Based on the sample projection data of different levels of noise, the initial imaging model is trained to reduce noise, so as to improve the noise reduction effect of the imaging model on different levels of noise.
- FIG. 5 is a schematic structural diagram of a CT imaging apparatus provided in Embodiment 3 of the present application.
- the apparatus includes: a projection data acquisition module 310 configured to acquire projection data of a target object, wherein the projection data is based on a preset dose X-ray photons are collected when the target object is detected, and the preset dose is less than the CT standard dose; the CT image reconstruction module 320 is configured to send the projection data to the pre-trained imaging model, and the imaging model
- the output image of is determined to be the CT image of the target object, wherein the imaging model is obtained by training a standard dose of CT image and sample projection data containing noise, and the imaging model includes a first sub-network module and a second sub-network Module and a third sub-network module.
- the first sub-network module is used for noise reduction processing on projection data in the projection domain
- the second sub-network module is used for transforming the processed projection data to generate CT domain data
- the third sub-network module is configured to perform noise reduction processing on the CT domain data to generate a CT image of the target object.
- the first sub-network module includes a first preset number of convolution modules connected in sequence, and each convolution module in the first preset number of convolution modules includes a convolution layer and an activation function Layer;
- the third sub-network module includes a second preset number of convolution modules connected in sequence, and each convolution module in the second preset number of convolution modules includes a convolution layer and an activation function layer.
- the convolution kernel of the convolution layer in the first sub-network module is a ⁇ b, where a and b are both positive integers greater than or equal to 1, and b>a;
- the third sub-network module The convolution kernel of the convolution layer in the network module is m ⁇ m, where m is a positive integer greater than or equal to 1.
- the input end of the first sub-network module is short-circuited with the output end of the first sub-network module, and the input end of the third sub-network module is connected to the output end of the third sub-network module Short.
- the second sub-network module transforms the processed projection data according to the following formula:
- f(x, y) is the CT domain image output by the second sub-network module
- x and y are the abscissa and ordinate in the CT domain image, respectively
- p(r, ⁇ ) is the input to the second sub-network
- the processed projection data of the module r is the distance between the processed projection data and the origin, and ⁇ is the projection angle.
- the device further includes: an initial imaging model establishment module configured to establish an initial imaging model; a reconstructed image generation module configured to input noise-containing sample projection data into the initial imaging model to obtain a reconstructed image; loss The function determining module is configured to determine a loss function according to the reconstructed image and the standard dose CT image corresponding to the sample projection data containing noise; the imaging model determining module is configured to adjust the initial imaging model according to the loss function Of the network parameters to generate the imaging model.
- an initial imaging model establishment module configured to establish an initial imaging model
- a reconstructed image generation module configured to input noise-containing sample projection data into the initial imaging model to obtain a reconstructed image
- loss The function determining module is configured to determine a loss function according to the reconstructed image and the standard dose CT image corresponding to the sample projection data containing noise
- the imaging model determining module is configured to adjust the initial imaging model according to the loss function Of the network parameters to generate the imaging model.
- the device further includes: a standard projection data acquisition module configured to acquire a standard dose CT image before inputting noise-containing sample projection data into the initial imaging model, and compare the standard dose CT image Preprocessing is performed to obtain standard projection data; the sample projection data acquisition module is configured to add noise data to the standard projection data to generate sample projection data containing noise.
- a standard projection data acquisition module configured to acquire a standard dose CT image before inputting noise-containing sample projection data into the initial imaging model, and compare the standard dose CT image Preprocessing is performed to obtain standard projection data
- the sample projection data acquisition module is configured to add noise data to the standard projection data to generate sample projection data containing noise.
- the sample projection data acquisition module is configured to: set a Poisson function of at least two noise levels; add noise to the standard projection data based on the Poisson function of the at least two noise levels to generate At least two noise-containing sample projection data corresponding to the CT image of the standard dose.
- the CT imaging device provided in the embodiment of the present application can execute the CT imaging method provided in any embodiment of the present application, and has the corresponding functional modules and beneficial effects for executing the CT imaging method.
- FIG. 6 is a schematic structural diagram of a medical imaging system provided in the fourth embodiment of the present application.
- FIG. 6 shows a block diagram of an exemplary medical imaging system suitable for implementing the embodiments of the present application.
- the medical imaging system shown in FIG. 6 is only It is an example, and should not bring any limitation to the function and scope of use of the embodiments of the present application.
- the medical imaging system includes a medical imaging device 500 and a computer 600.
- the computer 600 may be configured to implement specific methods and devices disclosed in some embodiments of the present application.
- the specific device in this embodiment uses a functional block diagram to show a hardware platform including a display module.
- the computer 600 may implement some embodiments of the present application through hardware devices, software programs, firmware, and combinations thereof.
- the computer 600 may be a general purpose computer or a special purpose computer.
- the computer 600 may include an internal communication bus 601, a processor (processor) 602, a read-only memory (Read-Only Memory, ROM) 603, a random access memory (Random Access Memory, RAM) 604, and a communication port 605, input/output component 606, hard disk 607, and user interface 608.
- the internal communication bus 601 can implement data communication between the components of the computer 600.
- the processor 602 can make a judgment and issue a prompt.
- the processor 602 may be composed of one or more processors.
- the communication port 605 can implement data communication between the computer 600 and other components (not shown in the figure), such as external devices, image acquisition devices, databases, external storage, and image processing workstations.
- the computer 600 can send and receive information and data from the network through the communication port 605.
- the input/output component 606 supports the input/output data flow between the computer 600 and other components.
- the user interface 608 can implement interaction and information exchange between the computer 600 and the user.
- the computer 600 may also include different forms of program storage units and data storage units, such as hard disk 607, ROM603, RAM604, capable of storing various data files used for computer processing and/or communication, and possible program instructions executed by the processor 602 .
- the processor 602 can be configured to execute a CT imaging method when executing a program, the method includes: acquiring projection data of a target object, wherein the projection data is based on a preset dose of X-ray photons on the target object It is acquired during detection, the preset dose is less than the CT standard dose; the projection data is sent to a pre-trained imaging model, and the output image of the imaging model is determined as the CT image of the target object, wherein The imaging model is obtained by training standard dose CT images and sample projection data containing noise.
- the imaging model includes a first sub-network module, a second sub-network module, and a third sub-network module.
- the first sub-network module uses To perform noise reduction processing on projection data in the projection domain, the second sub-network module is used to transform the processed projection data to generate CT domain data, and the third sub-network module is used to perform processing on the CT domain data. Perform noise reduction processing to generate a CT image of the target object.
- One or more aspects of this application can be illustrated and described in multiple categories or situations, including any new and useful process, machine, product, or combination of substances, or any new and useful improvement to them.
- one or more aspects of the present application may be completely executed by hardware, may be completely executed by software (including firmware, resident software, microcode, etc.), or may be executed by a combination of hardware and software.
- the above hardware or software can be called “data block”, “module”, “sub-network module”, “engine”, “unit”, “sub-unit”, “component” or “system”.
- one or more aspects of the present application may be embodied as a computer product located in one or more computer-readable media, and the product includes computer-readable program codes.
- the fifth embodiment of the present application provides a computer-readable storage medium that stores a computer program that, when executed by a processor, implements the CT imaging method as provided in all application embodiments of this application, and the method includes: acquiring a target object Projection data, wherein the projection data is collected when the target object is detected based on a preset dose of X-ray photons, the preset dose is less than the CT standard dose; the projection data is sent to the pre-trained imaging Model, the output image of the imaging model is determined as the CT image of the target object, wherein the imaging model is obtained by training a CT image of a standard dose and sample projection data containing noise, and the imaging model includes a first sub A network module, a second sub-network module, and a third sub-network module, the first sub-network module is used to perform noise reduction processing on projection data in the projection domain, and the second sub-network module is used to perform processing on the processed projection data Transformation is performed to generate CT domain data, and the third sub-network module is used to
- the computer-readable signal medium may include a propagated data signal containing a computer program code, for example on a baseband or as part of a carrier wave.
- the propagated signal may have multiple manifestations, including electromagnetic forms, optical forms, etc., or suitable combinations.
- the computer-readable signal medium may be any computer-readable medium other than a computer-readable storage medium, and the medium may be connected to an instruction execution system, apparatus, or device to realize communication, propagation, or transmission of a program for use.
- the program code located on the computer-readable signal medium can be propagated through any suitable medium, including radio, cable, fiber optic cable, radio frequency signal, or similar medium, or any combination of the foregoing medium.
- the computer program codes required for the operation of one or more parts of this application can be written in any one or more programming languages, including object-oriented programming languages such as Java, Scala, Smalltalk, Eiffel, JADE, Emerald, C++, C#, VB. NET, Python, etc., conventional programming languages such as C language, Visual Basic, Fortran 2003, Perl, COBOL 2002, PHP, ABAP, dynamic programming languages such as Python, Ruby and Groovy, or other programming languages.
- the program code can be run entirely on the user's computer, or run as an independent software package on the user's computer, or partly run on the user's computer and partly run on a remote computer, or run entirely on the remote computer or server.
- the remote computer can be connected to the user's computer through any network, such as a local area network (LAN) or a wide area network (WAN), or connected to an external computer (for example, via the Internet), or In the cloud computing environment, or used as a service such as Software-as-a-Service (SaaS).
- LAN local area network
- WAN wide area network
- SaaS Software-as-a-Service
- the numerical parameter should consider the prescribed effective digits and adopt the method of general digit retention.
- the numerical ranges and parameters used to confirm the breadth of the range in some embodiments of the present application are approximate values, in specific embodiments, the setting of such numerical values is as accurate as possible within the feasible range.
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Optics & Photonics (AREA)
- Heart & Thoracic Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
Claims (11)
- 一种电子计算机断层扫描CT成像方法,包括:An electronic computed tomography CT imaging method, including:获取目标对象的投影数据,其中,所述投影数据为基于预设剂量的X射线光子对所述目标对象进行检测时采集得到,所述预设剂量小于CT标准剂量;Acquiring projection data of the target object, where the projection data is collected when detecting the target object based on a preset dose of X-ray photons, and the preset dose is less than the CT standard dose;将所述投影数据发送至预先训练的成像模型,将所述成像模型的输出图像确定为所述目标对象的CT图像,其中,所述成像模型由标准剂量的CT图像和包含噪声的样本投影数据训练得到,所述成像模型包括第一子网络模块、第二子网络模块和第三子网络模块,所述第一子网络模块用于对投影域的投影数据进行降噪处理,所述第二子网络模块用于对处理后的投影数据进行变换,生成CT域数据,所述第三子网络模块用于对所述CT域数据进行降噪处理,生成所述目标对象的CT图像。The projection data is sent to a pre-trained imaging model, and the output image of the imaging model is determined as the CT image of the target object, wherein the imaging model is composed of a standard dose CT image and sample projection data containing noise After training, the imaging model includes a first sub-network module, a second sub-network module, and a third sub-network module. The first sub-network module is used to perform noise reduction processing on projection data in the projection domain, and the second The sub-network module is used to transform the processed projection data to generate CT domain data, and the third sub-network module is used to perform noise reduction processing on the CT domain data to generate a CT image of the target object.
- 根据权利要求1所述的方法,其中,所述第一子网络模块包括依次连接的第一预设数量的卷积模块,所述第一预设数量的卷积模块中每个卷积模块包括卷积层和激活函数层;The method according to claim 1, wherein the first sub-network module comprises a first preset number of convolution modules connected in sequence, and each convolution module of the first preset number of convolution modules includes Convolutional layer and activation function layer;所述第三子网络模块包括依次连接的第二预设数量的卷积模块,所述第二预设数量的卷积模块中每个卷积模块包括卷积层和激活函数层。The third sub-network module includes a second preset number of convolution modules connected in sequence, and each convolution module in the second preset number of convolution modules includes a convolution layer and an activation function layer.
- 根据权利要求2所述的方法,其中,所述第一子网络模块中的卷积层的卷积核为a×b,其中,a和b均为大于或等于1的正整数,b>a;The method according to claim 2, wherein the convolution kernel of the convolution layer in the first sub-network module is a×b, wherein a and b are both positive integers greater than or equal to 1, and b>a ;所述第三子网络模块中的卷积层的卷积核为m×m,其中,m为大于或等于1的正整数。The convolution kernel of the convolution layer in the third sub-network module is m×m, where m is a positive integer greater than or equal to 1.
- 根据权利要求1至3任一项所述的方法,其中,所述第一子网络模块的输入端与所述第一子网络模块的输出端短接,所述第三子网络模块的输入端与所述第三子网络模块的输出端短接。The method according to any one of claims 1 to 3, wherein the input end of the first sub-network module is short-circuited with the output end of the first sub-network module, and the input end of the third sub-network module is Short-circuit with the output terminal of the third sub-network module.
- 根据权利要求1至4任一项所述的方法,其中,所述第二子网络模块根据如下公式对所述处理后的投影数据进行变换:The method according to any one of claims 1 to 4, wherein the second sub-network module transforms the processed projection data according to the following formula:其中,f(x,y)为所述第二子网络模块输出的所述CT域图像,x和y分别为所述CT域图像中的横坐标和纵坐标,p(r,θ)为输入所述第二子网络模块的所述处理后的投影数据,r为所述处理后的投影数据与原点的距离,δ为投影X射线,θ为投影角度。Wherein, f(x, y) is the CT domain image output by the second sub-network module, x and y are the abscissa and ordinate in the CT domain image, respectively, and p(r, θ) is the input For the processed projection data of the second sub-network module, r is the distance between the processed projection data and the origin, δ is the projection X-ray, and θ is the projection angle.
- 根据权利要求1至5任一项所述的方法,其中,所述成像模型的训练方法,包括:The method according to any one of claims 1 to 5, wherein the training method of the imaging model comprises:建立初始成像模型;Establish an initial imaging model;将所述包含噪声的样本投影数据输入至所述初始成像模型中,得到重建图像;Input the noise-containing sample projection data into the initial imaging model to obtain a reconstructed image;根据所述重建图像和所述包含噪声的样本投影数据对应的标准剂量的CT图像确定损失函数,根据所述损失函数调节所述初始成像模型中的网络参数,生成所述成像模型。A loss function is determined according to the reconstructed image and the CT image of the standard dose corresponding to the sample projection data containing noise, and the network parameters in the initial imaging model are adjusted according to the loss function to generate the imaging model.
- 根据权利要求6所述的方法,在所述将所述包含噪声的样本投影数据输入至所述初始成像模型中之前,还包括:The method according to claim 6, before said inputting said noise-containing sample projection data into said initial imaging model, further comprising:获取所述标准剂量的CT图像,对所述标准剂量的CT图像进行预处理,得到标准投影数据;Acquiring the CT image of the standard dose, and preprocessing the CT image of the standard dose to obtain standard projection data;对所述标准投影数据增加噪声数据,生成所述包含噪声的样本投影数据。Noise data is added to the standard projection data to generate the noise-containing sample projection data.
- 根据权利要求7所述的方法,其中,所述对所述标准投影数据增加噪声,生成所述包含噪声的样本投影数据,包括:8. The method according to claim 7, wherein said adding noise to said standard projection data to generate said noise-containing sample projection data comprises:设置至少两个噪声等级的泊松函数;Set the Poisson function of at least two noise levels;基于所述至少两个噪声等级的泊松函数对所述标准投影数据增加噪声,生成所述标准剂量的CT图像对应的至少两个包含噪声的样本投影数据。Noise is added to the standard projection data based on the Poisson function of the at least two noise levels, and at least two noise-containing sample projection data corresponding to the CT image of the standard dose are generated.
- 一种电子计算机断层扫描CT成像装置,包括:An electronic computed tomography CT imaging device, including:投影数据获取模块,设置为获取目标对象的投影数据,其中,所述投影数据为基于预设剂量的X射线光子对所述目标对象进行检测时采集得到,所述预设剂量小于CT标准剂量;A projection data acquisition module configured to acquire projection data of a target object, wherein the projection data is acquired when the target object is detected based on a preset dose of X-ray photons, and the preset dose is less than the CT standard dose;CT图像重建模块,设置为将所述投影数据发送至预先训练的成像模型,将所述成像模型的输出图像确定为所述目标对象的CT图像,其中,所述成像模型由标准剂量的CT图像和包含噪声的样本投影数据训练得到,所述成像模型包括第一子网络模块、第二子网络模块和第三子网络模块,所述第一子网络模块用于对投影域的投影数据进行降噪处理,所述第二子网络模块用于对处理后的投影数据进行变换,生成CT域数据,所述第三子网络模块用于对所述CT域数据进行降噪处理,生成所述目标对象的CT图像。The CT image reconstruction module is configured to send the projection data to a pre-trained imaging model, and determine the output image of the imaging model as the CT image of the target object, wherein the imaging model consists of a standard dose of CT image And the sample projection data containing noise is obtained by training, the imaging model includes a first sub-network module, a second sub-network module and a third sub-network module, the first sub-network module is used to reduce the projection data of the projection domain Noise processing, the second sub-network module is used to transform the processed projection data to generate CT domain data, and the third sub-network module is used to perform noise reduction processing on the CT domain data to generate the target CT image of the subject.
- 一种计算机可读存储介质,存储有计算机程序,所述程序被处理器执行时实现如权利要求1-8中任一所述的电子计算机断层扫描CT成像方法。A computer-readable storage medium storing a computer program, which when executed by a processor, realizes the computerized tomography CT imaging method according to any one of claims 1-8.
- 一种医学成像系统,包括医学成像设备和计算机设备,所述计算机设备包括存储器、至少一个处理器及存储在所述存储器上并可在所述至少一个处理器上运行的计算机程序,所述至少一个处理器执行所述程序时可设置为执行 权利要求1-8任一所述的电子计算机断层扫描CT成像方法。A medical imaging system includes a medical imaging device and a computer device. The computer device includes a memory, at least one processor, and a computer program stored in the memory and running on the at least one processor. When a processor executes the program, it can be set to execute the computerized tomography CT imaging method of any one of claims 1-8.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910941729.1 | 2019-09-30 | ||
CN201910941729.1A CN112581554B (en) | 2019-09-30 | 2019-09-30 | CT imaging method, device, storage equipment and medical imaging system |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2021062885A1 true WO2021062885A1 (en) | 2021-04-08 |
Family
ID=75116282
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2019/111028 WO2021062885A1 (en) | 2019-09-30 | 2019-10-14 | Ct imaging method and apparatus, storage medium, and medical imaging system |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN112581554B (en) |
WO (1) | WO2021062885A1 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023279316A1 (en) * | 2021-07-08 | 2023-01-12 | 深圳高性能医疗器械国家研究院有限公司 | Pet reconstruction method based on denoising score matching network |
CN114972118B (en) * | 2022-06-30 | 2023-04-28 | 抖音视界有限公司 | Noise reduction method and device for inspection image, readable medium and electronic equipment |
CN116593504B (en) * | 2023-07-17 | 2023-10-03 | 中国科学院深圳先进技术研究院 | CT imaging method, device, equipment and storage medium |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090262996A1 (en) * | 2008-04-17 | 2009-10-22 | Samsonov Alexey A | Method for image reconstruction employing sparsity-constrained iterative correction |
CN102737392A (en) * | 2012-06-07 | 2012-10-17 | 南方医科大学 | Non-partial regularization prior reconstruction method for low-dosage X-ray captive test (CT) image |
CN102906791A (en) * | 2010-05-19 | 2013-01-30 | 威斯康星校友研究基金会 | Method for radiation dose reduction using prior image constrained image reconstruction |
CN105194692A (en) * | 2010-07-12 | 2015-12-30 | 通用电气医疗集团股份有限公司 | X-ray imaging at low contrast agent concentrations and/or low dose radiation |
CN106844524A (en) * | 2016-12-29 | 2017-06-13 | 北京工业大学 | A kind of medical image search method converted based on deep learning and Radon |
WO2018200493A1 (en) * | 2017-04-25 | 2018-11-01 | The Board Of Trustees Of The Leland Stanford Junior University | Dose reduction for medical imaging using deep convolutional neural networks |
CN109166161A (en) * | 2018-07-04 | 2019-01-08 | 东南大学 | A kind of low-dose CT image processing system inhibiting convolutional neural networks based on noise artifacts |
CN109272472A (en) * | 2018-10-15 | 2019-01-25 | 天津大学 | Noise and artifact eliminating method towards medical power spectrum CT image |
US20190076102A1 (en) * | 2017-09-13 | 2019-03-14 | LiteRay Medical, LLC | Systems and methods for ultra low dose ct fluoroscopy |
CN109509235A (en) * | 2018-11-12 | 2019-03-22 | 深圳先进技术研究院 | Method for reconstructing, device, equipment and the storage medium of CT image |
CN109613462A (en) * | 2018-11-21 | 2019-04-12 | 河海大学 | A kind of scaling method of CT imaging |
CN109717886A (en) * | 2017-10-30 | 2019-05-07 | 上海交通大学 | A kind of CT scanning method of low radiation dose |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106373163B (en) * | 2016-08-29 | 2019-05-28 | 东南大学 | A kind of low-dose CT imaging method indicated based on three-dimensional projection's distinctive feature |
WO2019128660A1 (en) * | 2017-12-29 | 2019-07-04 | 清华大学 | Method and device for training neural network, image processing method and device and storage medium |
CN110211194A (en) * | 2019-05-21 | 2019-09-06 | 武汉理工大学 | A method of sparse angular CT imaging artefacts are removed based on deep learning |
-
2019
- 2019-09-30 CN CN201910941729.1A patent/CN112581554B/en active Active
- 2019-10-14 WO PCT/CN2019/111028 patent/WO2021062885A1/en active Application Filing
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090262996A1 (en) * | 2008-04-17 | 2009-10-22 | Samsonov Alexey A | Method for image reconstruction employing sparsity-constrained iterative correction |
CN102906791A (en) * | 2010-05-19 | 2013-01-30 | 威斯康星校友研究基金会 | Method for radiation dose reduction using prior image constrained image reconstruction |
CN105194692A (en) * | 2010-07-12 | 2015-12-30 | 通用电气医疗集团股份有限公司 | X-ray imaging at low contrast agent concentrations and/or low dose radiation |
CN102737392A (en) * | 2012-06-07 | 2012-10-17 | 南方医科大学 | Non-partial regularization prior reconstruction method for low-dosage X-ray captive test (CT) image |
CN106844524A (en) * | 2016-12-29 | 2017-06-13 | 北京工业大学 | A kind of medical image search method converted based on deep learning and Radon |
WO2018200493A1 (en) * | 2017-04-25 | 2018-11-01 | The Board Of Trustees Of The Leland Stanford Junior University | Dose reduction for medical imaging using deep convolutional neural networks |
US20190076102A1 (en) * | 2017-09-13 | 2019-03-14 | LiteRay Medical, LLC | Systems and methods for ultra low dose ct fluoroscopy |
CN109717886A (en) * | 2017-10-30 | 2019-05-07 | 上海交通大学 | A kind of CT scanning method of low radiation dose |
CN109166161A (en) * | 2018-07-04 | 2019-01-08 | 东南大学 | A kind of low-dose CT image processing system inhibiting convolutional neural networks based on noise artifacts |
CN109272472A (en) * | 2018-10-15 | 2019-01-25 | 天津大学 | Noise and artifact eliminating method towards medical power spectrum CT image |
CN109509235A (en) * | 2018-11-12 | 2019-03-22 | 深圳先进技术研究院 | Method for reconstructing, device, equipment and the storage medium of CT image |
CN109613462A (en) * | 2018-11-21 | 2019-04-12 | 河海大学 | A kind of scaling method of CT imaging |
Also Published As
Publication number | Publication date |
---|---|
CN112581554A (en) | 2021-03-30 |
CN112581554B (en) | 2024-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107133996B (en) | Method for generating an attenuation map for PET data reconstruction and PET/CT system | |
CN107516330B (en) | Model generation method, image processing method and medical imaging equipment | |
US20210074036A1 (en) | Deep encoder-decoder models for reconstructing biomedical images | |
WO2021062885A1 (en) | Ct imaging method and apparatus, storage medium, and medical imaging system | |
CN110809782A (en) | Attenuation correction system and method | |
US11790578B2 (en) | System and method for computed tomography | |
JP6746676B2 (en) | Image processing apparatus, image processing method, and program | |
CN111540025B (en) | Predicting images for image processing | |
WO2017031982A1 (en) | Methods, systems, and media for noise reduction in computed tomography images | |
JP2021013725A (en) | Medical apparatus | |
JP6222813B2 (en) | X-ray computed tomography apparatus, image processing apparatus and image processing method | |
WO2013097390A1 (en) | Attenuation correction method and device for image in pet system | |
US10314558B2 (en) | Image processing apparatus, image processing system, image processing method, and recording medium | |
JPWO2017104700A1 (en) | Image processing apparatus and image processing method | |
CN107928690B (en) | Method for determining proportion of main rays, image de-scattering method, system and storage medium | |
KR20170087320A (en) | Tomography apparatus and method for reconstructing a tomography image thereof | |
US20200240934A1 (en) | Tomography apparatus and controlling method for the same | |
CN113096210A (en) | Image reconstruction method and device, electronic equipment and storage medium | |
CN111862255A (en) | Regularization image reconstruction method, system, readable storage medium and device | |
CN107019496A (en) | Electron density information acquisition method, device and equipment | |
WO2005076038A1 (en) | Method and apparatus for providing high-resolution reconstruction of an observed object | |
EP4190243A1 (en) | Image processing device, image processing method, learning device, learning method, and program | |
CN113643394A (en) | Scattering correction method, device, computer equipment and storage medium | |
CN114511666A (en) | Model generation method, image reconstruction method, device, equipment and medium | |
CN102499705B (en) | Method and system for eliminating ring artifact in tomography |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 19947983 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 19947983 Country of ref document: EP Kind code of ref document: A1 |
|
32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 17.10.2022) |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 19947983 Country of ref document: EP Kind code of ref document: A1 |