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 PDF

Info

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
Application number
PCT/CN2019/111028
Other languages
French (fr)
Chinese (zh)
Inventor
葛永帅
梁栋
朱炯滔
刘新
郑海荣
Original Assignee
中国科学院深圳先进技术研究院
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 中国科学院深圳先进技术研究院 filed Critical 中国科学院深圳先进技术研究院
Publication of WO2021062885A1 publication Critical patent/WO2021062885A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [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

Disclosed are a CT imaging method and apparatus, a storage medium, and a medical imaging system. The CT imaging method comprises: acquiring projection data of a target object, the projection data being collected when the target object is detected on the basis of a preset dose of X-ray photons, and the preset dose being smaller than a CT standard dose (S110); and sending the projection data to a pre-trained imaging model, and determining an output image of the imaging model as a CT image of the target object, wherein the imaging model is trained by a CT image at the standard dose and sample projection data which contains noise. The imaging model comprises 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 on projection data of a projection domain, the second sub-network module is used to transform the processed projection data so as to generate CT domain data, and the third sub-network module is used to perform noise reduction on the CT domain data so as to generate a CT image of the target object (S120).

Description

CT成像方法、装置、存储介质及医学成像系统CT imaging method, device, storage medium and medical imaging system
本申请要求在2019年09月30日提交中国专利局、申请号为201910941729.1的中国专利申请的优先权,该申请的全部内容通过引用结合在本申请中。This application claims the priority of the Chinese patent application filed with the Chinese Patent Office with application number 201910941729.1 on September 30, 2019. The entire content of this application is incorporated into this application by reference.
技术领域Technical field
本申请实施例涉及医学成像技术,例如涉及一种电子计算机断层扫描(Computed Tomography,CT)成像方法、装置、存储介质及医学成像系统。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.
背景技术Background technique
基于X射线的CT成像是辅助诊断的常用方式,X摄像具有放射性,会增加被检测对象的患癌风险。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.
为了减少被检对象所受的辐射剂量,低剂量的CT检测成为CT成像领域的研究方向之一。低剂量的CT重建算法是基于迭代重建技术开发的,将CT图像重建模拟为一种数学优化问题,迭代重建算法是在估计的图像上通过前向投影模拟一个综合投影,该过程中尽量模拟真实CT系统中X射线光子穿过被检测对象到达探测器的过程,将综合投影与探测器采集的真实测量值进行比较,并根据两者之间的差值确定下一次更新,以根据下一次更新对当前估计得到的图像进行校正。In order to reduce the radiation dose received by 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 process of X-ray photons passing through the detected object to the detector in the CT system. 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.
但是上述迭代重建算法存在重建时间长的问题,尤其是复杂的计算无法实现实时重建。However, the above iterative reconstruction algorithm has the problem of long reconstruction time, especially complex calculations cannot achieve real-time reconstruction.
发明内容Summary of the invention
本申请提供一种CT成像方法、装置、存储介质及医学成像系统,以实现提高CT图像重建效率和质量。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.
本申请实施例提供了一种CT成像方法,包括:The embodiment of the present application provides a 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.
本申请实施例还提供了一种CT成像装置,该装置包括:An embodiment of the present application also provides a CT imaging device, which includes:
投影数据获取模块,设置为获取目标对象的投影数据,其中,所述投影数据为基于预设剂量的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.
本申请实施例还提供了一种计算机可读存储介质,存储有计算机程序,该程序被处理器执行时实现如本申请任一实施例提供的CT成像方法。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.
本申请实施例还提供了一种医学成像系统,包括医学成像设备和计算机设备,所述计算机设备包括存储器、一个或多个处理器及存储在存储器上并可在处理器上运行的计算机程序,该程序被处理器执行时实现如本申请任一实施例提供的CT成像方法。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.
附图说明Description of the drawings
图1是本申请实施例一提供的一种CT成像方法的流程示意图;FIG. 1 is a schematic flowchart of a CT imaging method provided in Embodiment 1 of the present application;
图2是本申请实施例一提供的一种成像模型的结构示意图;2 is a schematic structural diagram of an imaging model provided in Embodiment 1 of the present application;
图3是本申请实施例一提供的另一种成像模型的结构示意图;3 is a schematic structural diagram of another imaging model provided in Embodiment 1 of the present application;
图4是本申请实施例二提供的另一种CT成像方法的流程示意图;4 is a schematic flowchart of another CT imaging method provided in Embodiment 2 of the present application;
图5是本申请实施例三提供的一种CT成像装置的结构示意图;FIG. 5 is a schematic structural diagram of a CT imaging apparatus according to Embodiment 3 of the present application;
图6是本申请实施例四提供的一种医学成像系统的结构示意图。FIG. 6 is a schematic structural diagram of a medical imaging system provided by Embodiment 4 of the present application.
具体实施方式Detailed ways
下面结合附图和实施例对本申请进行说明。本文所描述的具体实施例仅仅用于解释本申请,而非对本申请的限定。为了便于描述,附图中仅示出了与本申请相关的部分而非全部结构。The application will be described below with reference to the drawings and embodiments. The specific embodiments described herein are only used to explain the application, but not to limit the application. For ease of description, the drawings only show a part of the structure related to the present application instead of all of the structure.
实施例一Example one
图1为本申请实施例一提供的一种CT成像方法的流程示意图,本实施例可适用于基于成像模型快速生成高质量的低剂量CT图像的情况,该方法可以由本申请实施例提供的CT成像装置来执行,该CT成像方法包括S110-S120。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.
S110、获取目标对象的投影数据,其中,所述投影数据为基于预设剂量的X射线光子对所述目标对象进行检测时采集得到,所述预设剂量小于CT标准剂量。S110. Obtain 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.
S120、将所述投影数据发送至预先训练的成像模型,将所述成像模型的输出图像确定为所述目标对象的CT图像。其中,所述成像模型由标准剂量的CT图像和包含噪声的样本投影数据训练得到,所述成像模型包括第一子网络模块、第二子网络模块和第三子网络模块,所述第一子网络模块用于对投影域的投影数据进行降噪处理,所述第二子网络模块用于对处理后的投影数据进行变换,生成CT域数据,所述第三子网络模块用于对所述CT域数据进行降噪处理,生成所述目标对象的CT图像。S120. 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. Wherein, 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, and 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.
在本实施例中,目标对象可以是人或者动物,通过低剂量的X射线光子对目标对象进行检测,可降低检测过程中X射线对目标对象的放射损伤,降低目标对象的患癌风险。本实施例中,成像模型为预先训练的,具有对低剂量投影数据进行重建的功能和对低剂量投影数据中的噪声进行降噪处理的功能,将通过低剂量X射线光子对目标对象进行检测,采集得到的投影数据,输入至成像模型,可生成目标对象的降噪后的高质量CT图像。由于成像模型是预先训练的,可自动对输入的投影数据进行处理,输出目标对象的CT图像,无需对每一组投影数据进行迭代运算,降低了CT图像重建的计算量和难度,提高了CT图像重建效率,同时,成像模型对投影数据进行降噪处理,提高了低剂量检测生成图像的清晰度。In this embodiment, 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. In this embodiment, 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. There is no need to iteratively calculate each set of projection data, which reduces the calculation and difficulty of CT image reconstruction, and improves CT The image reconstruction efficiency, meanwhile, the imaging model performs noise reduction processing on the projection data, which improves the clarity of the image generated by low-dose detection.
本实施例中,成像模型包括第一子网络模块、第二子网络模块和第三子网络模块,示例性的,参见图2,图2是本申请实施例提供的一种成像模型的结构示意图。第一子网络模块用于对输入的投影图像进行卷积滤波处理,投影域为所有投影图像构成的域。第一子网络模块可以是包括依次连接的第一预设数量的卷积模块,所述第一预设数量的卷积模块中每个卷积模块包括卷积层和激活函数层,第一子网络模块可对投影数据进行降噪和滤波处理。一实施例中,第一预设数量可以是5、6、10等。一实施例中,所述第一子网络模块中的卷积层的卷积核为a×b,a和b均为大于或等于1的正整数,b>a。示例性的,第一子网络模块中的卷积层的卷积核可以是1×30、3×30、5×30、3×33、5×33等,通过采用b>a形式的卷积核可同时达到对投影数据的降噪和滤波功能。示例性 的,激活函数层可以是带泄露线性整流函数(Leaky Rectified Linear Unit,Leaky ReLU)。In this embodiment, the imaging model includes a first sub-network module, a second sub-network module, and a third sub-network module. For example, refer to FIG. 2, which 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. In an embodiment, 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. Exemplarily, 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. Exemplarily, the activation function layer may be a Leaky Rectified Linear Unit (Leaky ReLU).
第二子卷积模型用于将投影域的投影数据转换为CT域数据,CT域为所有CT图像构成的域。一实施例中,所述第二子网络模块根据如下公式对处理后的投影数据进行变换:
Figure PCTCN2019111028-appb-000001
该公式中,f(x,y)为第二子网络模块输出的CT域图像,x和y分别为所述CT域图像中的横坐标和纵坐标,p(r,θ)为输入第二子网络模块的投影数据,r为投影数据与原点的距离,δ为投影X射线,θ为投影角度。第二子网络模块根据上述公式对第一子网络模块输入的投影数据进行重建,输出与期望尺寸相同的CT域图像,该CT域图像满足CT图像属性。
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. In an embodiment, the second sub-network module transforms the processed projection data according to the following formula:
Figure PCTCN2019111028-appb-000001
In this 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, and 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, and θ 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.
第三子网络模块用于对CT域图像进行卷积滤波操作,同时对得到的CT域图像进行图像还原处理和降噪处理。一实施例中,对图像进行还原可以是包括依次连接的第二预设数量的卷积模块,所述第二预设数量的卷积模块中每个卷积模块包括卷积层和激活函数层。一实施例中,第二预设数量可以是与第一预设数量相同或不同,第二预设数量例如可以是5、6或者10等,第二预设数量和第一预设数量可以是分别根据功能精度确定。示例性的,激活函数层可以是Leaky ReLU,第三子网络模块中的卷积层的卷积核为m×m,所述m为大于或等于1的正整数,例如,第三子网络模块中的卷积层的卷积核可以是3×3、5×5或者7×7等。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. In an embodiment, 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. . In an embodiment, 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. Exemplarily, the activation function layer may be Leaky ReLU, the convolution kernel of the convolution layer in the third sub-network module is m×m, and the m is a positive integer greater than or equal to 1, for example, the third sub-network module The convolution kernel of the convolution layer in can be 3×3, 5×5, or 7×7.
示例性的,参见图3,图3是本申请实施例一提供的另一种成像模型的结构示意图。图3中的第一子网络模块包括6个卷积模块,卷积模块中卷积层的卷积核为3×30,步长为1,第三子网络模块包括6个卷积模块,卷积模块中卷积层的卷积核为3×3,步长为1。图3中输入投影数据的尺寸为900×848,输出CT图像的尺寸为512×512,上述图像尺寸仅为图3示例,输入投影数据的尺寸可根据成像模型的输入层确定,在输入投影数据之前,若采集的投影数据的尺寸与成像模块输入层的处理尺寸不一致时,预先调节该采集的投影数据,使该采集的投影数据符合成像模块输入层的处理尺寸。Exemplarily, refer to FIG. 3, which 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.
在上述实施例的基础上,所述第一子网络模块的输入端与所述第一子网络模块的输出端短接,所述第三子网络模块的输入端与所述第三子网络模块的输出端短接。通过将第一子网络模块和第三子网络模块的网络层进行短接,其一,对于任一子网络模块,通过短接将输入信息传输至输出端,避免在处理过程中的信息丢失,其二,在训练过程中,避免反向训练过程中,梯度下降导致的训练失败,提高了训练效率和质量。On the basis of the foregoing embodiment, 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. By short-circuiting the network layers of the first sub-network module and the third sub-network module, firstly, for any sub-network module, the input information is transmitted to the output terminal through the short connection, so as to avoid the loss of information in the processing process, Second, in the training process, it avoids the training failure caused by the gradient descent during the reverse training process, which improves the training efficiency and quality.
本实施例提供的技术方案,通过设置三个子网络模块,分别对投影域数据 和CT域数据进行去噪处理,实现在图像重建过程中的双重去噪,提高了去噪效果,使得在低剂量X射线光子进行检测时,得到高质量的CT图像。同时三个子网络模块相互连接,彼此独立,可根据成像需求对任一模块进行更新和置换,提高了成像模型的可配置性和适用性。In the technical solution provided by this embodiment, by setting three sub-network modules, 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. At the same time, 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.
实施例二Example two
图4是本申请实施例二提供的一种CT成像方法的流程示意图,在上述实施例的基础上,提供了成像模型的训练方法,该CT成像方法包括S210-S250。Fig. 4 is a schematic flow chart of a CT imaging method provided in the second embodiment of the present application. On the basis of the foregoing embodiment, a training method of an imaging model is provided. The CT imaging method includes S210-S250.
S210、建立初始成像模型。S210. Establish an initial imaging model.
S220、将包含噪声的样本投影数据输入至初始成像模型中,得到重建图像。S220: Input the sample projection data containing noise into the initial imaging model to obtain a reconstructed image.
S230、根据所述重建图像和所述包含噪声的样本投影数据对应的标准剂量的CT图像确定损失函数,根据所述损失函数调节所述初始成像模型中的网络参数,生成所述成像模型。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.
S240、获取目标对象的投影数据,其中,所述投影数据为基于预设剂量的X射线光子对所述目标对象进行检测时采集得到,所述预设剂量小于CT标准剂量。S240. Obtain 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.
S250、将所述投影数据发送至预先训练的成像模型,将所述成像模型的输出图像确定为所述目标对象的CT图像。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.
在本实施例中,基于样本数据对初始成像模型进行训练,以得到具有图像重建功能和去噪功能成像模型。一实施例中,样本数据中包括包含噪声的低剂量投影数据和标准剂量的CT图像。包含噪声的低剂量投影数据可以是通过低剂量的X射线光子对样本对象进行检测时采集的,还可以是基于标准剂量的CT图像进行加噪处理得到的。In this embodiment, the initial imaging model is trained based on sample data to obtain an imaging model with image reconstruction function and denoising function. In an embodiment, 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.
一实施例中,在将包含噪声的样本投影数据输入至初始成像模型中之前,获取包含噪声的样本投影数据可以是:获取标准剂量的CT图像,对所述标准剂量的CT图像进行预处理,得到标准投影数据;对所述标准投影数据增加噪声数据,生成包含噪声的样本投影数据。一实施例中,对标准剂量的CT图像进行预处理可以是对标准剂量的CT图像进行Radon操作,生成标准投影数据,对标准投影数据中添加噪声,例如可以是泊松噪声,得到包含噪声的样本投影数据。本实施例中,标准剂量的CT图像可以是临床检测得到的CT图像,通过对标准投影数据添加噪声的方式,生成样本投影数据,避免样本生成过程中对样本对象进行X射线扫描导致样本对象的患癌风险。一实施例中,对标准投影数据添加不同基本的噪声,生成不同的样本投影数据,增加样本中噪声的多样性,以 及增大样本数量。一实施例中,所述对所述标准投影数据增加噪声,生成包含噪声的样本投影数据,包括:设置至少两个噪声等级的泊松函数;基于所述至少两个噪声等级的泊松函数对所述标准投影数据增加噪声,生成所述标准剂量的CT图像对应的至少两个包含噪声的样本投影数据。一实施例中,通过Radon操作得到标准剂量的CT图像的正投影数据,即正弦图像(sino图像),通过如下方式生成包含噪声的样本投影数据,sino_ps=(poissrnd(exp(-sino).*I 1,[del_num,angle_num])+1)./I 1,sino_ps为包含噪声的样本投影数据,del_num为探测器数量,angle_num为投影总角度数,I 1=I 0×k,k为噪声强度系数,I 0为基础光子数,I 0=1×10 6。一实施例中,k≤1,例如可以是1、0.1、0.2、0.5或者0.05等,对k值不做限定,通过k值调节添加的噪声等级,k越小噪声越强。本实施例中,通过在标准投影数据添加不同等级的噪声,形成具有不同等级噪声的样本投影数据。基于不同等级噪声的样本投影数据对初始成像模型进行降噪训练,提高成像模型对不同量级的噪声的降噪效果。 In an embodiment, before inputting the noise-containing sample projection data into the initial imaging model, 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. In one embodiment, 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. In this embodiment, 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. In one embodiment, 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. In an embodiment, 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. In one embodiment, the orthographic projection data of the CT image of the standard dose, namely the sine image (sino image), is obtained through the Radon operation, and the sample projection data containing noise is generated in the following manner, sino_ps=(poissrnd(exp(-sino).* I 1 ,[del_num,angle_num])+1)./I 1 , 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 is noise Intensity coefficient, I 0 is the basic photon number, I 0 =1×10 6 . In an embodiment, k≤1, for example, 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. In this embodiment, by adding different levels of noise to the standard projection data, 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.
实施例三Example three
图5是本申请实施例三提供的一种CT成像装置的结构示意图,该装置包括:投影数据获取模块310,设置为获取目标对象的投影数据,其中,所述投影数据为基于预设剂量的X射线光子对所述目标对象进行检测时采集得到,所述预设剂量小于CT标准剂量;CT图像重建模块320,设置为将所述投影数据发送至预先训练的成像模型,将所述成像模型的输出图像确定为所述目标对象的CT图像,其中,所述成像模型由标准剂量的CT图像和包含噪声的样本投影数据训练得到,所述成像模型包括第一子网络模块、第二子网络模块和第三子网络模块,所述第一子网络模块用于对投影域的投影数据进行降噪处理,所述第二子网络模块用于对处理后的投影数据进行变换,生成CT域数据,所述第三子网络模块用于对所述CT域数据进行降噪处理,生成所述目标对象的CT图像。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, and 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.
一实施例中,所述第一子网络模块包括依次连接的第一预设数量的卷积模块,所述第一预设数量的卷积模块中每个卷积模块包括卷积层和激活函数层;所述第三子网络模块包括依次连接的第二预设数量的卷积模块,所述第二预设数量的卷积模块中每个卷积模块包括卷积层和激活函数层。In an embodiment, 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.
一实施例中,所述第一子网络模块中的卷积层的卷积核为a×b,其中,a和b均为大于或等于1的正整数,b>a;所述第三子网络模块中的卷积层的卷积核为m×m,其中,所述m为大于或等于1的正整数。In an embodiment, 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.
一实施例中,所述第一子网络模块的输入端与所述第一子网络模块的输出端短接,所述第三子网络模块的输入端与所述第三子网络模块的输出端短接。In an embodiment, 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.
一实施例中,所述第二子网络模块根据如下公式对处理后的投影数据进行变换:In an embodiment, the second sub-network module transforms the processed projection data according to the following formula:
Figure PCTCN2019111028-appb-000002
Figure PCTCN2019111028-appb-000002
其中,f(x,y)为第二子网络模块输出的CT域图像,x和y分别为所述CT域图像中的横坐标和纵坐标,p(r,θ)为输入第二子网络模块的处理后的投影数据,r为处理后的投影数据与原点的距离,θ为投影角度。Among them, 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 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.
一实施例中,所述装置还包括:初始成像模型建立模块,设置为建立初始成像模型;重建图像生成模块,设置为将包含噪声的样本投影数据输入至初始成像模型中,得到重建图像;损失函数确定模块,设置为根据所述重建图像和所述包含噪声的样本投影数据对应的标准剂量的CT图像确定损失函数;成像模型确定模块,设置为根据所述损失函数调节所述初始成像模型中的网络参数,生成所述成像模型。In one embodiment, 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.
一实施例中,所述装置还包括:标准投影数据获取模块,设置为在将包含噪声的样本投影数据输入至初始成像模型中之前,获取标准剂量的CT图像,对所述标准剂量的CT图像进行预处理,得到标准投影数据;样本投影数据获取模块,设置为对所述标准投影数据增加噪声数据,生成包含噪声的样本投影数据。In an embodiment, 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.
一实施例中,所述样本投影数据获取模块是设置为:设置至少两个噪声等级的泊松函数;基于所述至少两个噪声等级的泊松函数对所述标准投影数据增加噪声,生成所述标准剂量的CT图像对应的至少两个包含噪声的样本投影数据。In an embodiment, 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.
本申请实施例提供的CT成像装置可执行本申请任意实施例所提供的CT成像方法,具备执行CT成像方法相应的功能模块和有益效果。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.
实施例四Example four
图6是本申请实施例四提供的一种医学成像系统的结构示意图,图6示出了适于用来实现本申请实施方式的示例性医学成像系统的框图,图6显示的医学成像系统仅仅是一个示例,不应对本申请实施例的功能和使用范围带来任何限制。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.
医学成像系统包括医学成像设备500和计算机600。The medical imaging system includes a medical imaging device 500 and a computer 600.
计算机600可以被设置为实现实施本申请一些实施例中披露的特定方法和装置。本实施例中的特定装置利用功能框图展示了一个包含显示模块的硬件平台。在一些实施例中,计算机600可以通过硬件设备、软件程序、固件以及它 们的组合来实现本申请一些实施例的实施。在一些实施例中,计算机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. In some embodiments, the computer 600 may implement some embodiments of the present application through hardware devices, software programs, firmware, and combinations thereof. In some embodiments, the computer 600 may be a general purpose computer or a special purpose computer.
如图6所示,计算机600可以包括内部通信总线601,处理器(processor)602,只读存储器(Read-Only Memory,ROM)603,随机存取存储器(Random Access Memory,RAM)604,通信端口605,输入/输出组件606,硬盘607,以及用户界面608。内部通信总线601可以实现计算机600组件间的数据通信。处理器602可以进行判断和发出提示。在一些实施例中,处理器602可以由一个或多个处理器组成。通信端口605可以实现计算机600与其他部件(图中未示出)例如:外接设备、图像采集设备、数据库、外部存储以及图像处理工作站等之间进行数据通信。在一些实施例中,计算机600可以通过通信端口605从网络发送和接受信息及数据。输入/输出组件606支持计算机600与其他部件之间的输入/输出数据流。用户界面608可以实现计算机600和用户之间的交互和信息交换。计算机600还可以包括不同形式的程序储存单元以及数据储存单元,例如硬盘607,ROM603,RAM604,能够存储计算机处理和/或通信使用的多种数据文件,以及处理器602所执行的可能的程序指令。As shown in FIG. 6, 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. In some embodiments, 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. In some embodiments, 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 .
所述处理器602执行程序时可设置为执行一种CT成像方法,所述方法包括:获取目标对象的投影数据,其中,所述投影数据为基于预设剂量的X射线光子对所述目标对象进行检测时采集得到,所述预设剂量小于CT标准剂量;将所述投影数据发送至预先训练的成像模型,将所述成像模型的输出图像确定为所述目标对象的CT图像,其中,所述成像模型由标准剂量的CT图像和包含噪声的样本投影数据训练得到,所述成像模型包括第一子网络模块、第二子网络模块和第三子网络模块,所述第一子网络模块用于对投影域的投影数据进行降噪处理,所述第二子网络模块用于对处理后的投影数据进行变换,生成CT域数据,所述第三子网络模块用于对所述CT域数据进行降噪处理,生成所述目标对象的CT图像。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.
本申请使用了特定词语来描述本申请的实施例。如“一个实施例”、“一实施例”、和/或“一些实施例”意指与本申请至少一个实施例相关的一个特征、结构或特点。因此,本说明书中在不同位置两次或多次提及的“一实施例”或“一个实施例”或“一替代性实施例”并不一定是指同一实施例。此外,本申请的一个或多个实施例中的一些特征、结构或特点可以进行适当的组合。This application uses specific words to describe the embodiments of this application. For example, "one embodiment", "an embodiment", and/or "some embodiments" mean a feature, structure, or characteristic related to at least one embodiment of the present application. Therefore, "an embodiment" or "an embodiment" or "an alternative embodiment" mentioned twice or more in different positions in this specification do not necessarily refer to the same embodiment. In addition, some features, structures, or characteristics in one or more embodiments of the present application can be appropriately combined.
本申请的一个或多个方面可以通过多个种类或情况进行说明和描述,包括任何新的和有用的工序、机器、产品或物质的组合,或对他们的任何新的和有用的改进。相应地,本申请的一个或多个方面可以完全由硬件执行、可以完全由软件(包括固件、常驻软件、微码等)执行、也可以由硬件和软件组合执行。 以上硬件或软件均可被称为“数据块”、“模块”、“子网络模块”、“引擎”、“单元”、“子单元”、“组件”或“系统”。此外,本申请的一个或多个方面可能表现为位于一个或多个计算机可读介质中的计算机产品,该产品包括计算机可读程序编码。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. Correspondingly, 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". In addition, 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.
实施例五Example five
本申请实施例五提供了一种计算机可读存储介质,存储有计算机程序,该程序被处理器执行时实现如本申请所有申请实施例提供的CT成像方法,所述方法包括:获取目标对象的投影数据,其中,所述投影数据为基于预设剂量的X射线光子对所述目标对象进行检测时采集得到,所述预设剂量小于CT标准剂量;将所述投影数据发送至预先训练的成像模型,将所述成像模型的输出图像确定为所述目标对象的CT图像,其中,所述成像模型由标准剂量的CT图像和包含噪声的样本投影数据训练得到,所述成像模型包括第一子网络模块、第二子网络模块和第三子网络模块,所述第一子网络模块用于对投影域的投影数据进行降噪处理,所述第二子网络模块用于对处理后的投影数据进行变换,生成CT域数据,所述第三子网络模块用于对所述CT域数据进行降噪处理,生成所述目标对象的CT图像。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 perform noise reduction processing on the CT domain data to generate a CT image of the target object.
计算机可读信号介质可能包含一个内含有计算机程序编码的传播数据信号,例如在基带上或作为载波的一部分。该传播信号可能有多种表现形式,包括电磁形式、光形式等等、或合适的组合形式。计算机可读信号介质可以是除计算机可读存储介质之外的任何计算机可读介质,该介质可以通过连接至一个指令执行系统、装置或设备以实现通讯、传播或传输供使用的程序。位于计算机可读信号介质上的程序编码可以通过任何合适的介质进行传播,包括无线电、电缆、光纤电缆、射频信号、或类似介质、或任何上述介质的组合。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.
本申请一个或多个部分操作所需的计算机程序编码可以用任意一种或多种程序语言编写,包括面向对象编程语言如Java、Scala、Smalltalk、Eiffel、JADE、Emerald、C++、C#、VB.NET、Python等,常规程序化编程语言如C语言、Visual Basic、Fortran 2003、Perl、COBOL 2002、PHP、ABAP,动态编程语言如Python、Ruby和Groovy,或其他编程语言等。该程序编码可以完全在用户计算机上运行、或作为独立的软件包在用户计算机上运行、或部分在用户计算机上运行部分在远程计算机运行、或完全在远程计算机或服务器上运行。在后种情况下,远程计算机可以通过任何网络形式与用户计算机连接,比如局域网(Local Area Network,LAN)或广域网(Wide Area Network,WAN),或连接至外部计算机(例如通过因特网),或在云计算环境中,或作为服务使用如软件即服务 (Software-as-a-Service,SaaS)。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. In the latter case, 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).
除非权利要求中说明,本申请所述处理元素和序列的顺序、数字字母的使用、或其他名称的使用,并非用于限定本申请流程和方法的顺序。尽管上述披露中通过多种示例讨论了一些申请实施例,该类细节仅起到说明的目的。例如,虽然以上所描述的系统组件可以通过硬件设备实现,但是也可以只通过软件的解决方案得以实现,如在服务器或移动设备上安装所描述的系统。Unless stated in the claims, the order of processing elements and sequences, the use of numbers and letters, or the use of other names in this application are not used to limit the order of the procedures and methods of this application. Although some application embodiments are discussed in the foregoing disclosure through various examples, such details are only for illustrative purposes. For example, although the system components described above can be implemented by hardware devices, they can also be implemented only by software solutions, such as installing the described system on a server or mobile device.
为了简化本申请披露的表述,从而帮助对一个或多个申请实施例的理解,本文对本申请实施例的描述中,有时会将多种特征归并至一个实施例、附图或对该多种特征的描述中。一些实施例中使用了描述成分、属性数量的数字,此类用于实施例描述的数字,在一些示例中使用了修饰词“大约”、“近似”或“大体上”来修饰。除非另外说明,“大约”、“近似”或“大体上”表明所述数字允许有±20%的变化。相应地,在一些实施例中,使用的数值参数均为近似值,该近似值根据个别实施例所需特点可以发生改变。在一些实施例中,数值参数应考虑规定的有效数位并采用一般位数保留的方法。尽管本申请一些实施例中用于确认范围广度的数值域和参数为近似值,在具体实施例中,此类数值的设定在可行范围内尽可能精确。In order to simplify the expressions disclosed in this application, and to help the understanding of one or more application embodiments, in the description of the embodiments of this application, sometimes multiple features are combined into one embodiment, drawings, or the multiple features. In the description. In some embodiments, numbers describing the quantities of ingredients and attributes are used. Such numbers used in the description of the embodiments are modified by the modifiers "about", "approximately" or "substantially" in some examples. Unless otherwise stated, "approximately", "approximately" or "substantially" indicates that the number is allowed to vary by ±20%. Correspondingly, in some embodiments, the numerical parameters used are approximate values, and the approximate values can be changed according to the required characteristics of individual embodiments. In some embodiments, the numerical parameter should consider the prescribed effective digits and adopt the method of general digit retention. Although 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.

Claims (11)

  1. 一种电子计算机断层扫描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.
  2. 根据权利要求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.
  3. 根据权利要求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.
  4. 根据权利要求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.
  5. 根据权利要求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:
    Figure PCTCN2019111028-appb-100001
    Figure PCTCN2019111028-appb-100001
    其中,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.
  6. 根据权利要求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.
  7. 根据权利要求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.
  8. 根据权利要求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.
  9. 一种电子计算机断层扫描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.
  10. 一种计算机可读存储介质,存储有计算机程序,所述程序被处理器执行时实现如权利要求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.
  11. 一种医学成像系统,包括医学成像设备和计算机设备,所述计算机设备包括存储器、至少一个处理器及存储在所述存储器上并可在所述至少一个处理器上运行的计算机程序,所述至少一个处理器执行所述程序时可设置为执行 权利要求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.
PCT/CN2019/111028 2019-09-30 2019-10-14 Ct imaging method and apparatus, storage medium, and medical imaging system WO2021062885A1 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (12)

* Cited by examiner, † Cited by third party
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