WO2019047545A1 - 低剂量成像方法及装置 - Google Patents
低剂量成像方法及装置 Download PDFInfo
- Publication number
- WO2019047545A1 WO2019047545A1 PCT/CN2018/085628 CN2018085628W WO2019047545A1 WO 2019047545 A1 WO2019047545 A1 WO 2019047545A1 CN 2018085628 W CN2018085628 W CN 2018085628W WO 2019047545 A1 WO2019047545 A1 WO 2019047545A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image
- projection data
- displaying
- iterative reconstruction
- iteration
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 claims abstract description 105
- 230000008569 process Effects 0.000 claims description 67
- 238000007499 fusion processing Methods 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 15
- 230000004927 fusion Effects 0.000 claims description 5
- 238000004904 shortening Methods 0.000 abstract description 4
- 230000000875 corresponding effect Effects 0.000 description 16
- 210000000988 bone and bone Anatomy 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 6
- 210000004872 soft tissue Anatomy 0.000 description 6
- 238000002591 computed tomography Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- 230000036541 health Effects 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 238000013441 quality evaluation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000001959 radiotherapy Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007408 cone-beam computed tomography Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 210000001835 viscera Anatomy 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/17—Circuit arrangements not adapted to a particular type of detector
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/424—Iterative
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/444—Low dose acquisition or reduction of radiation dose
Definitions
- the present invention relates to the field of image processing, and in particular, to a low dose imaging method and apparatus.
- X-ray computed tomography (CT) technology for diagnostics, and cone beam (CB) CT technology for radiotherapy or surgical guidance have been widely used in clinical medical imaging diagnosis and radiotherapy, but Excessive doses of X-ray radiation during CT and CBCT scans can have an impact on patient health.
- the imaging radiation dose to the imaged object can be reduced while the image quality is guaranteed, and then the data generated after the X-ray passes through the target area of the patient can be captured and used as projection data.
- a target image for clinical treatment is obtained, which is also referred to as a low-dose imaging mode.
- the low-dose imaging mode one frame of projection data is acquired for each irradiation.
- the display of the target image needs to wait for the projection data of all the frames to be collected, and the time taken to display the target image is long.
- the embodiment of the invention provides a low-dose imaging method and device, which can solve the problem that the display target image needs to wait for all the frame projection data to be collected after the target image is collected, and the time taken to display the target image is long.
- the technical solution is as follows:
- a low dose imaging method comprising:
- the performing, by the second processing based on the preset number of projection data, generating the second image and displaying the second image includes:
- i is 1, 2, ..., m
- the m is the total number of iterative reconstruction operations in the iterative reconstruction process
- the m is an integer greater than or equal to 1.
- the displaying the i-th iteration to reconstruct the image comprises:
- the performing, by the second processing based on the preset number of projection data, generating the second image and displaying the second image includes:
- performing image fusion processing on the ith iteration reconstructed image and the first image to generate the ith fused image including:
- the displaying the i-th iteration to reconstruct the image comprises:
- the i-th iteration reconstructed image and the progress information from the mth iteration reconstructed image are displayed.
- the displaying the second image includes:
- a low dose imaging device comprising:
- a first processing module configured to perform first processing on the collected projection data, generate a first image, and display the first image, before the data amount of the collected projection data reaches a preset data amount;
- a second processing module configured to perform a second process based on the preset number of projection data, generate a second image, and display the second image when the data amount of the collected projection data reaches a preset data amount.
- the second processing module is configured to:
- i is 1, 2, ..., m
- the m is the total number of iterative reconstruction operations in the iterative reconstruction process
- the m is an integer greater than or equal to 1.
- the second processing module is configured to:
- the second processing module is configured to:
- the second processing module is configured to:
- the second processing module is configured to:
- the second processing module is configured to:
- a low dose imaging apparatus comprising a memory, a processor, and a computer program stored in the memory and operative on the processor, the processor executing the computer program as The steps of the method of the first aspect.
- a computer readable storage medium storing a computer program, the computer program being executed by a processor to implement the steps of the method of the first aspect.
- a computer program product in a fifth aspect, storing instructions that, when run on a computer, cause the computer to perform the low dose imaging method of the first aspect.
- a chip comprising programmable logic circuitry and/or program instructions for implementing the low dose imaging method of the first aspect when the chip is in operation.
- the low-dose imaging method and device continuously collect projection data, and perform the first processing on the collected projection data to generate a first image before the collected data amount of the projection data reaches a preset data amount. Displaying a first image, and when the amount of data of the collected projection data reaches a preset data amount, performing a second process based on the preset number of projection data, generating a second image, and displaying the second image, compared to related technologies, It does not need to wait for the projection data of all frames to be displayed after the image is captured, which shortens the time taken to display the image.
- FIG. 1 is a flow chart of a low dose imaging method according to an embodiment of the present invention.
- FIG. 2 is a flowchart of generating a first image and displaying a first image according to an embodiment of the present invention
- FIG. 3 is a flowchart of generating an analytical reconstructed image according to an embodiment of the present invention.
- FIG. 4 is a flowchart of generating a second image and displaying a second image according to an embodiment of the present invention
- FIG. 5 is a flowchart of generating a second image and displaying a second image according to an embodiment of the present invention
- FIG. 6 is a flowchart of generating a fused image according to an embodiment of the present invention.
- FIG. 7 is a schematic structural diagram of a low dose imaging apparatus according to an embodiment of the present invention.
- FIG. 8 is a schematic structural diagram of another low dose imaging apparatus according to an embodiment of the present invention.
- image reconstruction technology refers to using certain energy waves, such as X-rays, positrons, ultrasounds, etc., to illuminate a patient's target region from various directions, and capture data generated after the energy wave passes through the target region of the patient, and Using this as projection data, the projection data is calculated by a specific algorithm to reconstruct a tomographic image including the target region of the patient.
- X-ray CT technology is used to make X-rays pass through human tissues (such as human internal organs) from different directions, and then capture the data generated by X-rays passing through human tissues, and use the data as projection data.
- a human tomographic image is reconstructed based on the projection data. Excessive doses of X-ray radiation can affect the health of patients. In order to reduce the impact on the health of patients, a low dose imaging method relative to conventional dose imaging is proposed.
- the analytical reconstruction method is used to reconstruct the human tomographic image.
- the analytical reconstruction method is built on the continuous signal model, which is sensitive to noise and requires complete projection data.
- the analytical reconstruction method may be a Fourier transform method or a filtered back projection method; in the low dose imaging mode, an iterative reconstruction method is used to reconstruct a human tomographic image, and an iterative reconstruction method is established on a discrete signal model.
- the iterative reconstruction method can reconstruct a high-quality human tomographic image with low signal-to-noise ratio (the signal-to-noise ratio is often low in low-dose imaging mode) and incomplete projection data.
- the iterative reconstruction method may be an Algebra Reconstruction Technique (ART), an Ordered Subsets Expectation Maximization (OSEM), or a Total Variation Algebra Reconstruction Technique (Total Variation Algebra Reconstruction Technique, TV-ART) or Maximum A Posteriori reconstuction (MAP-ICD).
- ART Algebra Reconstruction Technique
- OSEM Ordered Subsets Expectation Maximization
- MAP-ICD Maximum A Posteriori reconstuction
- the quality of the acquired projection data per frame is poor.
- it is usually obtained after the projection data of all the frames is acquired. Iterative reconstruction processing is performed on the projection data of all the frames, thereby generating a target image and displaying it.
- the process can only perform clinical treatment through the final target image.
- the display target image needs to wait for the projection data of all the frames to be collected, and the time taken to display the target image is longer. The staff was wasting time on unnecessary waiting.
- the collected projection data before the collected data volume of the projection data reaches a preset data amount, the collected projection data is first processed, the first image is generated, and the first image is displayed.
- the second processing is performed based on the preset number of projection data, the second image is generated, and the second image is displayed.
- Embodiments of the present invention provide a low-dose imaging method for a low-dose imaging device having a display function in an imaging system
- the imaging system may further include an imaging source (such as a bulb), the imaging source emits an energy wave, and The energy wave passes through the body tissue from different directions, and the imager (such as the detector plate) captures the data generated by the energy wave passing through the human tissue, and uses the data as projection data.
- the energy wave may be an X-ray, a positron beam, an ultrasonic wave, or the like, and the type of the energy wave is not limited in the embodiment of the present invention.
- the low-dose imaging device collects the above projection data and processes and displays the projection data by using a low-dose imaging method. As shown in FIG. 1, the low-dose imaging method includes:
- Step 101 Continuously collect projection data.
- Step 102 Before the collected data volume of the projection data reaches a preset data amount, perform first processing on the collected projection data, generate a first image, and display the first image.
- Step 103 When the data amount of the collected projection data reaches a preset data amount, perform a second process based on the preset number of projection data, generate a second image, and display the second image.
- the above low dose imaging device may be a computer, a server or the like.
- the low-dose imaging method can continuously collect projection data, and perform the first processing on the collected projection data before the data amount of the collected projection data reaches a preset data amount.
- the first image is displayed and the first image is displayed, and when the data amount of the collected projection data reaches a preset data amount, the second processing is performed based on the preset number of projection data, and the second image is generated and the second image is displayed.
- the first image and the second image may be multiple types of images.
- the first image may be an analytical reconstructed image.
- the second image can be an iteratively reconstructed image.
- the second image may also be an image generated by performing image fusion processing on the first image and the iterative reconstructed image.
- the embodiment of the present invention does not limit the types of the first image and the second image.
- the method generates and displays the first image and the second image in the process of continuously acquiring the projection data, and does not need to wait for the projection data of all the frames to be displayed after the image is captured, shortening the time taken for displaying the image, and providing rich clinical treatment.
- the reference data avoids the staff wasting time on unnecessary waiting.
- performing the first processing on the collected projection data in step 102, generating the first image, and displaying the first image may include:
- Step 1021 Perform an analytical reconstruction process on the collected projection data to generate an analytical reconstructed image.
- Step 1022 Display the reconstructed image.
- the low-dose imaging device can analyze and reconstruct the collected 3 frames of projection data before the collected projection data reaches 10 frames, for example, the acquired projection data is 3 frames. Processing, generating an analytical reconstructed image and displaying the analytical reconstructed image.
- the low-dose imaging device may directly perform the analytical reconstruction processing on the collected projection data;
- the low-dose imaging device may first perform denoising processing on the collected projection data, and then perform analytical reconstruction processing on the collected projection data. Therefore, optionally, in this step, as shown in FIG. 3, performing an analytical reconstruction process on the collected projection data to generate an analytical reconstructed image may include:
- Step 1021a Perform denoising processing on the collected projection data to obtain processed projection data.
- Step 1021b Perform analysis and reconstruction processing on the processed projection data to generate an analytical reconstructed image.
- performing the analytical reconstruction processing on the processed projection data may include performing an analytical reconstruction process on the processed projection data by using a Fourier transform method or a filtered back projection method.
- the low-dose imaging device can display the analytical reconstructed image through the display, thereby providing reference data for clinical treatment, and the staff can perform preliminary clinical treatment tasks based on the analytical reconstructed image, such as coarse matching in image guidance.
- the quasi-task avoids the staff wasting time on unnecessary waiting.
- step 103 when the data amount of the collected projection data reaches a preset data amount, the second processing is performed based on the preset number of projection data, and the second image may be generated and displayed in a plurality of manners.
- the second process may be an iterative reconstruction process, and the generated and displayed second image is an iteratively reconstructed image; on the other hand, the second process may include an iterative reconstruction process and an image fusion process, and the generated and displayed second image includes iterations. Reconstruct images and fused images. Step 103 will be described below by taking these two aspects as an example.
- performing a second process based on the preset number of projection data in step 103, generating a second image, and displaying the second image may include:
- Step 1031 Perform an i-th iterative reconstruction operation in an iterative reconstruction process on the preset number of projection data to generate an i-th iterative reconstructed image.
- the low-dose imaging device may generate the i-th iterative reconstructed image by using an ART method, an OSEM method, a TV-ART method, or a MAP-ICD method.
- Step 1032 Display the i-th iteration to reconstruct the image.
- i is 1, 2, ..., m
- m is the total number of iterative reconstruction operations in the iterative reconstruction process, and m is an integer greater than or equal to 1.
- the total number m of iterative reconstruction operations in the iterative reconstruction process is equal to 3, and the low-dose imaging device performs the first iteration of the 10 frame projection data in the iterative reconstruction process when 10 frames of projection data is acquired.
- Iterative reconstruction operation generating the first iteration reconstructed image, and displaying the first iteration reconstructed image; performing the second iterative reconstruction operation in the iterative reconstruction process on the 10-frame projection data, generating the second iteration reconstructed image, and displaying The second iteration reconstructs the image; performs a third iteration reconstruction operation in the iterative reconstruction process on the 10-frame projection data, generates a third iteration reconstructed image, and displays the third iteration reconstructed image.
- the low dose imaging device can display iteratively reconstructed images to provide reference data for clinical treatment, and the staff can perform more detailed clinical treatment tasks based on iteratively reconstructed images.
- displaying the i-th iteration to reconstruct the image in step 1032 may include:
- the i-th iteration reconstructed image and the progress information of the reconstructed image from the mth iteration are displayed.
- the progress information of the reconstructed image from the mth iteration may be the duration from the generation time of the i-th iterative reconstructed image to the mth iteration of the reconstructed image.
- the low-dose imaging device acquires 10 frames of projection data, it performs the first iterative reconstruction operation in the iterative reconstruction process on the 10-frame projection data, generates the first iteration reconstructed image, and displays the reconstructed image from the third iteration.
- the low-dose imaging device displays the i-th iterative reconstructed image and the progress information of the reconstructed image from the mth iteration, so that the staff can timely understand the generation progress of the iterative reconstructed image, and facilitate the execution of the corresponding clinical treatment. task.
- the second processing is performed based on the preset number of projection data in the step 103, and the second image is generated and the second image is displayed, which may include:
- Step 1034 Perform image fusion processing on the i-th iterative reconstructed image and the first image to generate an i-th merged image and display the i-th merged image.
- Step 1035 Perform an mth iterative reconstruction operation on the projection data of the preset data amount in an iterative reconstruction process, generate an mth iteration reconstructed image, and display the mth iteration reconstructed image.
- the low-dose imaging device collects 10 frames of projection data, it performs the first iterative reconstruction operation in the iterative reconstruction process of the 10-frame projection data, generates the first iteration reconstructed image, and reconstructs the image with the first iteration and the first image.
- the image is subjected to image fusion processing to generate a first fused image and display the first fused image; the second iterative reconstruction operation in the iterative reconstruction process of the 10-frame projection data is generated, and the second iteration reconstructed image is generated, and the second image is reconstructed.
- the iterative reconstruction image and the first image are subjected to image fusion processing to generate a second fused image and display the second fused image.
- the low-dose imaging device performs a third iterative reconstruction operation in the iterative reconstruction process on the 10-frame projection data, generates a third iteration reconstructed image, and displays the third iteration reconstructed image.
- the first image used to perform image fusion processing may be an analytical reconstructed image. Since the noise of the reconstructed image is large, the quality of the edge information is high, and the noise of the iteratively reconstructed image is small, and the quality of the edge information is low, in the embodiment of the present invention, the low-dose imaging device can reconstruct the image from the ith iteration.
- the image fusion processing is performed on the analytically reconstructed image generated in step 1021, and an ith fusion image having less noise and high edge information quality is generated.
- the low-dose imaging device performs an analytical reconstruction process on the three-frame projection data to generate an analytical reconstructed image J.
- the low-dose imaging device performs the first iterative reconstruction operation on the 6-frame projection data, generates the first iteration reconstructed image D1, and performs image fusion processing on the iterative reconstructed image D1 and the analytical reconstructed image J.
- the first fused image B1 is generated and the first fused image B1 is displayed; the low-dose imaging device performs a second iterative reconstruction operation on the 6-frame projection data, generates a second iteration reconstructed image D2, and reconstructs the image from the second iteration.
- D2 performs image fusion processing with the analysis reconstructed image J, generates a second fused image B2, and displays the second fused image B2.
- the low-dose imaging device performs a third iteration reconstruction operation on the 6-frame projection data, generates a third iteration reconstructed image D3, and displays the third iteration reconstructed image D3. In this way, the low-dose imaging device displays a total of 2 fused images: B1 and B2, and an iterative reconstruction image: D3.
- displaying the ith fusion image in step 1034 may include:
- the i-th merged image and the progress information of the reconstructed image from the mth iteration are displayed.
- the progress information of the reconstructed image from the mth iteration may be the duration from the generation time of the i-th fused image to the generation time of the m-th iteration reconstructed image.
- the low-dose imaging device displays the i-th fusion image and the progress information of the reconstructed image from the m-th iteration, so that the staff can know the progress of the fusion image generation in time, and facilitate the execution of the corresponding clinical treatment task.
- the i-th iterative reconstructed image is image-fused with the first image to generate an i-th merged image, which may include:
- Step 1034a Determine, according to the number of times i of the iterative reconstruction operation, the weight of the pixel value (or pixel) in the i-th iterative reconstructed image and the weight of the pixel value in the first image.
- Step 1034b Perform a fusion process on the i-th iterative reconstructed image and the first image according to the weight of the pixel value in the image reconstructed in the i-th iteration and the weight of the pixel value in the first image to generate the i-th merged image.
- the sum of the weight of one pixel in the i-th iterative reconstructed image and the weight of the pixel fused with the pixel in the first image is 1, and the weight of the pixel value in the reconstructed image of the i-th iteration and the iterative reconstruction operation are
- the magnitude of the number of times i is positively correlated
- the weight of the pixel values in the first image is inversely related to the magnitude of the number of iterations of the iterative reconstruction operation. That is, the greater the number i of iterative reconstruction operations, the greater the weight of the pixel values in the reconstructed image of the i-th iteration, and the smaller the weight of the pixel values in the first image.
- the pixels in the image may include pixels corresponding to the soft tissue and pixels corresponding to the bone tissue, but the weight of the pixel value corresponding to the soft tissue in the image increases or decreases the rate and/or mode of the weight of the pixel value corresponding to the bone tissue. It can be different.
- the low-dose imaging device may increase the weight of the pixel value in the corresponding i-th iterative reconstructed image, including the weight of the pixel value of the pixel corresponding to the soft tissue.
- the weight of the pixel value of the pixel corresponding to the bone tissue but the weight and/or mode of the weight of the pixel value of the pixel corresponding to the soft tissue is different, and here, in particular, the pixel value of the pixel corresponding to the bone tissue is increased.
- the weight of the pixel value in the image reconstructed in the first iteration is q1, and the weight of the pixel value in the first image is p1; when the number i of the iterative reconstruction operation is equal to 2, The weight of the pixel value in the reconstructed image in the second iteration is q2, and the weight of the pixel value in the first image is p2, then q1 ⁇ q2, p1>p2.
- the displaying the second image in step 103 may include:
- the image quality evaluation method may be an objective evaluation method of digital image quality.
- the type of digital image quality objective evaluation method may be a full reference type (Full Reference, FR), a partial reference type (Reduced Reference, RR) or a no reference type (No Reference, NR).
- FR refers to the fact that the original image is known, and the quality of the current image is evaluated based on the original image.
- NR refers to predicting the quality of an overall image based on image local features with discriminative power in the current image without the original image.
- RR is between FR and NR, which refers to the use of partial information of the original image to evaluate the quality of the current image.
- the original image is the first image in the embodiment of the present invention
- the current image is the second image in the embodiment of the present invention.
- the first image may be an analytically reconstructed image
- the second image may include an iteratively reconstructed image
- the second image may also include a fused image.
- the low-dose imaging method can continuously collect projection data, and perform the first processing on the collected projection data before the data amount of the collected projection data reaches a preset data amount.
- the first image displays the first image, and when the data amount of the collected projection data reaches a preset data amount, performs a second process based on the preset number of projection data, generates a second image, and displays the second image, and further
- the progress information of the second image and the quality information of the second image can be displayed.
- the embodiment of the present invention provides a low-dose imaging device having a display function and disposed in an imaging system. As shown in FIG. 7, the device 700 includes:
- the acquisition module 710 is configured to continuously collect projection data.
- the first processing module 720 is configured to perform first processing on the collected projection data, generate a first image, and display the first image, before the data amount of the collected projection data reaches a preset data amount.
- the acquisition module continuously collects projection data
- the first processing module performs the collected projection data before the collected data amount of the projection data reaches a preset data amount.
- a first process the first image is generated and the first image is displayed
- the second processing module performs a second process based on the preset number of projection data to generate a second image when the data amount of the collected projection data reaches a preset data amount
- displaying the second image compared with the related art, there is no need to wait for the projection data of all the frames to be displayed after the image is captured, which shortens the time taken to display the image.
- the first processing module 720 is configured to:
- the second processing module 730 is configured to:
- i is 1, 2, ..., m
- m is the total number of iterative reconstruction operations in the iterative reconstruction process, and m is an integer greater than or equal to 1.
- the second processing module 730 is configured to:
- the i-th iteration reconstructed image and the progress information of the reconstructed image from the mth iteration are displayed.
- the second processing module 730 is configured to:
- the second processing module 730 is configured to:
- the i-th iterative reconstructed image and the first image are subjected to fusion processing according to the weight of the pixel value in the image reconstructed in the i-th iteration and the weight of the pixel value in the first image to generate an i-th merged image.
- the second processing module 730 is configured to:
- the i-th merged image and the progress information of the reconstructed image from the mth iteration are displayed.
- the second processing module 730 is configured to:
- the acquisition module continuously collects projection data
- the first processing module performs the collected projection data before the collected data amount of the projection data reaches a preset data amount.
- the first image is generated and the first image is displayed
- the second processing module performs a second process based on the preset number of projection data to generate a second image when the data amount of the collected projection data reaches a preset data amount
- displaying the second image compared with the related art, there is no need to wait for the projection data of all the frames to be displayed after the image is captured, which shortens the time taken to display the image.
- the embodiment of the invention further provides a low dose imaging device, as shown in FIG. 8, comprising:
- the low-dose imaging device can continuously collect projection data, and perform the first processing on the collected projection data before the data amount of the collected projection data reaches a preset data amount.
- the first image displays the first image, and when the data amount of the collected projection data reaches a preset data amount, performs a second process based on the preset number of projection data, generates a second image, and displays the second image, and further
- the progress information of the second image and the quality information of the second image can be displayed.
- the embodiment of the present invention further provides a computer readable storage medium, which is a non-volatile readable storage medium, where the computer readable storage medium stores a computer program, and the computer program is implemented by a processor The steps of the low dose imaging method provided by the above embodiments.
- Embodiments of the present invention also provide a computer program product having instructions stored therein that, when run on a computer, cause the computer to implement the steps of the low dose imaging method provided by the above embodiments.
- Embodiments of the present invention also provide a chip that includes programmable logic circuitry and/or program instructions for implementing the steps of the low dose imaging method provided by the above embodiments when the chip is in operation.
- a person skilled in the art may understand that all or part of the steps of implementing the above embodiments may be completed by hardware, or may be instructed by a program to execute related hardware, and the program may be stored in a computer readable storage medium.
- the storage medium mentioned may be a read only memory, a magnetic disk or an optical disk or the like.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种低剂量成像方法及装置,属于图像处理领域。该方法包括:持续采集投影数据;在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示所述第一图像;当采集到的投影数据的数据量达到预设数据量时,基于所述预设数量的投影数据进行第二处理,生成第二图像并显示所述第二图像。本发明解决了相关技术中显示目标图像需要等待所有帧的投影数据被采集完后进行,显示目标图像所用的时间较长的问题,达到了无需等待所有帧的投影数据被采集完后显示图像,缩短了显示图像所用的时间的效果,用于低剂量成像。
Description
本发明涉及图像处理领域,特别涉及一种低剂量成像方法及装置。
应用于诊断领域的X射线计算机断层扫描(Computed Tomography,CT)技术,以及应用于放疗或手术引导的锥形束(Cone Beam,CB)CT技术已经广泛应用于临床医学影像诊断、放疗中,但是CT与CBCT扫描过程中过高的X射线辐射剂量会对患者的健康造成影响。为了减小对患者的健康造成的影响,可以在保证成像质量的同时,降低对被成像对象的成像辐射剂量,然后捕获X射线穿过患者靶区后产生的数据,将其作为投影数据,再基于该投影数据得到用于进行临床治疗的目标图像,该方式也称作低剂量成像方式,在低剂量成像方式中,每进行一次辐射,就会获取到一帧投影数据。
然而由于X射线辐射剂量较低,所以获取的每帧投影数据的质量较差,为了保证最终得到的目标图像的质量,通常会在所有帧的投影数据都获取完毕后再对所有帧的投影数据进行迭代重建处理,进而生成用于进行临床治疗的目标图像,并对其进行显示。
上述方式中显示目标图像需要等待所有帧的投影数据被采集完后进行,显示目标图像所用的时间较长。
发明内容
本发明实施例提供了一种低剂量成像方法及装置,可以解决相关技术中显示目标图像需要等待所有帧的投影数据被采集完后进行,显示目标图像所用的时间较长的问题。所述技术方案如下:
第一方面,提供了一种低剂量成像方法,所述方法包括:
持续采集投影数据;
在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示所述第一图像;
当采集到的投影数据的数据量达到预设数据量时,基于所述预设数量的投影数据进行第二处理,生成第二图像并显示所述第二图像。
可选的,所述基于所述预设数量的投影数据进行第二处理,生成第二图像并显示所述第二图像,包括:
对所述预设数量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像;
显示所述第i次迭代重建图像;
其中,所述i为1,2,…,m,所述m为所述迭代重建处理中迭代重建操作的总次数,且所述m为大于等于1的整数。
可选的,所述显示所述第i次迭代重建图像,包括:
显示所述第i次迭代重建图像以及距生成所述第m次迭代重建图像的进度信息。
可选的,所述基于所述预设数量的投影数据进行第二处理,生成第二图像并显示所述第二图像,包括:
对所述预设数据量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像;
将所述第i次迭代重建图像与所述第一图像进行图像融合处理,生成所述第i次融合图像并显示所述第i次融合图像,其中,所述i为1,2,…,m-1,所述m为所述迭代重建处理中迭代重建操作的总次数,且所述m为大于等于2的整数;
对所述预设数据量的投影数据进行迭代重建处理中的第m次迭代重建操作,生成第m次迭代重建图像并显示第m次迭代重建图像。
可选的,将所述第i次迭代重建图像与所述第一图像进行图像融合处理,生成所述第i次融合图像,包括:
基于所述迭代重建操作的次数i,确定所述第i次迭代重建图像中像素值的权重以及所述第一图像中像素值的权重;
根据所述第i次迭代重建图像中像素值的权重和所述第一图像中像素值的权重,将所述第i次迭代重建图像和所述第一图像进行融合处理,生成所述第i次融合图像。
可选的,所述显示所述第i次迭代重建图像,包括:
显示所述第i次迭代重建图像以及距生成所述第m次迭代重建图像的进度 信息。
可选的,所述显示所述第二图像,包括:
显示所述第二图像以及所述第二图像的质量信息,所述质量信息用于表征所述第二图像相较于所述第一图像的图像质量。
第二方面,提供了一种低剂量成像装置,所述装置包括:
采集模块,用于持续采集投影数据;
第一处理模块,用于在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示所述第一图像;
第二处理模块,用于在采集到的投影数据的数据量达到预设数据量时,基于所述预设数量的投影数据进行第二处理,生成第二图像并显示所述第二图像。
可选的,所述第二处理模块,用于:
对所述预设数量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像;
显示所述第i次迭代重建图像;
其中,所述i为1,2,…,m,所述m为所述迭代重建处理中迭代重建操作的总次数,且所述m为大于等于1的整数。
可选的,所述第二处理模块,用于:
显示所述第i次迭代重建图像以及距生成所述第m次迭代重建图像的进度信息。
可选的,所述第二处理模块,用于:
对所述预设数据量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像;
将所述第i次迭代重建图像与所述第一图像进行图像融合处理,生成所述第i次融合图像并显示所述第i次融合图像,其中,所述i为1,2,…,m-1,所述m为所述迭代重建处理中迭代重建操作的总次数,且所述m为大于等于2的整数;
对所述预设数据量的投影数据进行迭代重建处理中的第m次迭代重建操作,生成第m次迭代重建图像并显示第m次迭代重建图像。
可选的,所述第二处理模块,用于:
基于所述迭代重建操作的次数i,确定所述第i次迭代重建图像中像素值的 权重以及所述第一图像中像素值的权重;
根据所述第i次迭代重建图像中像素值的权重和所述第一图像中像素值的权重,将所述第i次迭代重建图像和所述第一图像进行融合处理,生成所述第i次融合图像。
可选的,所述第二处理模块,用于:
显示所述第i次融合图像以及距生成所述第m次迭代重建图像的进度信息。
可选的,所述第二处理模块,用于:
显示所述第二图像以及所述第二图像的质量信息,所述质量信息用于表征所述第二图像相较于所述第一图像的图像质量。
第三方面,提供了一种低剂量成像装置,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如第一方面所述方法的步骤。
第四方面,提供了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如第一方面所述方法的步骤。
第五方面,提供了一种计算机程序产品,所述计算机程序产品中存储有指令,当其在计算机上运行时,使得计算机执行如第一方面所述的低剂量成像方法。
第六方面,提供了一种芯片,所述芯片包括可编程逻辑电路和/或程序指令,当所述芯片运行时用于实现如第一方面所述的低剂量成像方法。
本发明实施例提供的技术方案带来的有益效果是:
本发明实施例提供的低剂量成像方法及装置,持续采集投影数据,在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示第一图像,并在采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示第二图像,相较于相关技术,无需等待所有帧的投影数据被采集完后显示图像,缩短了显示图像所用的时间。
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所 需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种低剂量成像方法的流程图;
图2是本发明实施例提供的生成第一图像并显示第一图像的流程图;
图3是本发明实施例提供的一种生成解析重建图像的流程图;
图4是本发明实施例提供的一种生成第二图像并显示第二图像的流程图;
图5是本发明实施例提供的另一种生成第二图像并显示第二图像的流程图;
图6是本发明实施例提供的一种生成融合图像的流程图;
图7是本发明实施例提供的一种低剂量成像装置的结构示意图;
图8是本发明实施例提供的另一种低剂量成像装置的结构示意图。
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述,显然,所描述的实施例仅仅是本发明一部份实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
相关技术中,图像重建技术指的是采用某种能量波,如X射线、正电子射线、超声波等,从各个不同方向照射患者靶区,捕获能量波穿过患者靶区后产生的数据,并将其作为投影数据,再通过特定的算法对该投影数据进行计算,重建出包括患者靶区的断层图像。在图像处理领域中,X射线CT技术用于使X射线从不同方向穿过人体组织(比如人体内脏),然后捕获X射线穿过人体组织后产生的数据,并将该数据作为投影数据,再基于该投影数据重建出人体断层图像。由于过高的X射线辐射剂量会对患者的健康造成影响。为了减小对患者的健康造成的影响,提出了相对于常规剂量成像方式的低剂量成像方式。
通常,在常规剂量成像方式中,采用解析重建法来重建人体断层图像,解析重建法建立在连续信号模型上,对噪声比较敏感且要求投影数据完备。示例的,解析重建法可以为Fourier(傅立叶)变换法或滤波反投影法等;在低剂量成像方式中,采用迭代重建法来重建人体断层图像,迭代重建法建立在离散信号模型上,相较于解析重建法,迭代重建法在低信噪比(低剂量成像方式下的 信噪比常常较低)以及投影数据不完备的情况下,能重建出质量较高的人体断层图像。示例的,迭代重建法可以为代数重建法(Algebra Reconstruction Technique,ART)、有序子集最大期望值法(Ordered Subsets Expectation Maximization,OSEM)、最小化图像总变差迭代法(Total Variation Algebra Reconstruction Technique,TV-ART)或最大后验概率迭代梯度下降法(Maximum A Posteriori reconstuction,MAP-ICD)等。
在低剂量成像方式中,由于X射线辐射剂量较低,所以获取的每帧投影数据的质量较差,为了保证最终得到的目标图像的质量,通常会在所有帧的投影数据都获取完毕后再对所有帧的投影数据进行迭代重建处理,进而生成目标图像,并对其进行显示。但是该过程在等待较长时间后,仅能够通过最终的目标图像进行临床治疗,显示目标图像需要等待所有帧的投影数据被采集完后进行,显示目标图像所用的时间较长。工作人员将时间都浪费在了无谓的等待上面。
而在本发明实施例中,能够在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示该第一图像,当采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示该第二图像,相较于相关技术,无需等待所有帧的投影数据被采集完后显示图像,缩短了显示图像所用的时间,为临床治疗提供了丰富的参考数据,避免了工作人员将时间都浪费在了无谓的等待上面。
本发明实施例提供了一种低剂量成像方法,用于成像系统中的具有显示功能的低剂量成像装置,该成像系统还可以包括成像源(如球管),成像源发出能量波,并使能量波从不同方向穿过人体组织,成像器(如探测器平板)捕获能量波穿过人体组织后产生的数据,并将该数据作为投影数据。示例的,该能量波可以为X射线、正电子射线、超声波等,本发明实施例对该能量波的类型不作限定。低剂量成像装置采集上述投影数据,并利用低剂量成像方法对这些投影数据进行处理并显示,如图1所示,该低剂量成像方法包括:
步骤101、持续采集投影数据。
步骤102、在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示该第一图像。
步骤103、当采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示该第二图像。
上述低剂量成像装置可以为计算机、服务器等装置。
综上所述,本发明实施例提供的低剂量成像方法,能够持续采集投影数据,在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示第一图像,并在采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示第二图像,相较于相关技术,无需等待所有帧的投影数据被采集完后显示图像,缩短了显示图像所用的时间。在本发明实施例中,第一图像和第二图像可以为多种类型的图像。示例的,第一图像可以为解析重建图像。示例的,第二图像可以为迭代重建图像。第二图像也可以是对第一图像和迭代重建图像进行图像融合处理后生成的图像,本发明实施例对第一图像和第二图像的类型不做限定。该方法在持续采集投影数据的过程中,生成并显示第一图像和第二图像,无需等待所有帧的投影数据被采集完后显示图像,缩短了显示图像所用的时间,为临床治疗提供了丰富的参考数据,避免了工作人员将时间都浪费在了无谓的等待上面。
可选的,如图2所示,步骤102中对采集到的投影数据进行第一处理,生成第一图像并显示该第一图像,可以包括:
步骤1021、对采集到的投影数据进行解析重建处理,以生成解析重建图像。
步骤1022、显示解析重建图像。
假设预设数据量为10帧,低剂量成像装置在采集到的投影数据的数据量达到10帧之前,比如采集到的投影数据为3帧,那么可以对采集到的3帧投影数据进行解析重建处理,生成解析重建图像并显示该解析重建图像。
可选的,对采集到的投影数据进行解析重建处理时,当不需要对采集到的投影数据进行去噪处理时,低剂量成像装置可以直接对采集到的投影数据进行解析重建处理;而当需要对采集到的投影数据进行去噪处理时,低剂量成像装置可以先对采集到的投影数据进行去噪处理,然后对采集到的投影数据进行解析重建处理。所以,可选的,在本步骤,如图3所示,对采集到的投影数据进行解析重建处理,以生成解析重建图像,可以包括:
步骤1021a、对采集到的投影数据进行去噪处理,得到处理后的投影数据。
步骤1021b、对处理后的投影数据进行解析重建处理,以生成解析重建图 像。
示例的,对处理后的投影数据进行解析重建处理,可以包括:采用Fourier变换法或滤波反投影法对处理后的投影数据进行解析重建处理。
低剂量成像装置在生成解析重建图像后,可以通过显示器显示该解析重建图像,进而为临床治疗提供参考数据,工作人员可以基于该解析重建图像执行初步的临床治疗任务,比如图像引导中的粗配准任务,避免了工作人员将时间都浪费在了无谓的等待上面。
在步骤103中,当采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示该第二图像可以有多种方式。一方面,第二处理可以为迭代重建处理,生成并显示的第二图像为迭代重建图像;另一方面,第二处理可以包括迭代重建处理和图像融合处理,生成并显示的第二图像包括迭代重建图像和融合图像。下面以这两个方面为例对步骤103进行说明。
可选的,一方面,如图4所示,步骤103中基于预设数量的投影数据进行第二处理,生成第二图像并显示第二图像,可以包括:
步骤1031、对预设数量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像。
示例的,在本步骤中,低剂量成像装置可以采用ART法、OSEM法、TV-ART法或MAP-ICD法等方式,生成第i次迭代重建图像。
步骤1032、显示第i次迭代重建图像。
其中,i为1,2,…,m,m为迭代重建处理中迭代重建操作的总次数,且m为大于等于1的整数。
假设预设数据量为10帧,迭代重建处理中迭代重建操作的总次数m等于3,低剂量成像装置在采集到10帧投影数据时,对10帧投影数据进行迭代重建处理中的第1次迭代重建操作,生成第1次迭代重建图像,并显示该第1次迭代重建图像;对10帧投影数据进行迭代重建处理中的第2次迭代重建操作,生成第2次迭代重建图像,并显示该第2次迭代重建图像;对10帧投影数据进行迭代重建处理中的第3次迭代重建操作,生成第3次迭代重建图像,并显示该第3次迭代重建图像。
在本发明实施例中,低剂量成像装置可以显示迭代重建图像,为临床治疗提供参考数据,工作人员可以基于迭代重建图像执行更加细化的临床治疗任 务。
其中,在一种可实现方式中,步骤1032中显示第i次迭代重建图像,可以包括:
显示第i次迭代重建图像以及距生成第m次迭代重建图像的进度信息。
示例的,距生成第m次迭代重建图像的进度信息可以是从第i次迭代重建图像的生成时刻开始到第m次迭代重建图像的生成时刻的时长。
假设预设数据量为10帧,迭代重建处理中迭代重建操作的总次数m等于3那么,i依次为1,2,3。低剂量成像装置在采集到10帧投影数据时,对10帧投影数据进行迭代重建处理中的第1次迭代重建操作,生成第1次迭代重建图像,以及显示距生成第3次迭代重建图像的进度信息;对10帧投影数据进行迭代重建处理中的第2次迭代重建操作,生成第2次迭代重建图像,显示该第2次迭代重建图像,以及显示距生成第3次迭代重建图像的进度信息;对10帧投影数据进行迭代重建处理中的第3次迭代重建操作,生成第3次迭代重建图像,显示该第3次迭代重建图像,以及显示距生成第3次迭代重建图像的进度信息。
在本发明实施例中,低剂量成像装置显示第i次迭代重建图像以及距生成第m次迭代重建图像的进度信息,使得工作人员能够及时了解迭代重建图像的生成进度,便于执行对应的临床治疗任务。
另一方面,如图5所示,步骤103中基于预设数量的投影数据进行第二处理,生成第二图像并显示第二图像,可以包括:
步骤1033、对预设数据量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像。
步骤1034、将第i次迭代重建图像与第一图像进行图像融合处理,生成第i次融合图像并显示第i次融合图像。
其中,i为1,2,…,m-1,m为迭代重建处理中迭代重建操作的总次数,且m为大于等于2的整数。
步骤1035、对预设数据量的投影数据进行迭代重建处理中的第m次迭代重建操作,生成第m次迭代重建图像并显示第m次迭代重建图像。
假设预设数据量为10帧,迭代重建处理中迭代重建操作的总次数m等于3,那么i依次为1,2。低剂量成像装置在采集到10帧投影数据时,对10帧投影数据进行迭代重建处理中的第1次迭代重建操作,生成第1次迭代重建图 像,并将第1次迭代重建图像与第一图像进行图像融合处理,生成第1次融合图像并显示第1次融合图像;对10帧投影数据进行迭代重建处理中的第2次迭代重建操作,生成第2次迭代重建图像,并将第2次迭代重建图像与第一图像进行图像融合处理,生成第2次融合图像并显示第2次融合图像。低剂量成像装置再对10帧投影数据进行迭代重建处理中的第3次迭代重建操作,生成第3次迭代重建图像并显示第3次迭代重建图像。
在本发明实施例中,低剂量成像装置可以显示融合图像和第m次迭代重建图像,进而为临床治疗提供参考数据,工作人员可以基于融合图像和第m次迭代重建图像执行相应的临床治疗任务。
示例的,用于进行图像融合处理的第一图像可以为解析重建图像。由于解析重建图像的噪声较大,边缘信息质量较高,而迭代重建图像的噪声较小,边缘信息质量较低,所以在本发明实施例中,低剂量成像装置可以将第i次迭代重建图像与步骤1021中生成的解析重建图像进行图像融合处理,进而生成噪声较小,边缘信息质量较高的第i次融合图像。
假设预设数据量为6帧,迭代重建处理中迭代重建操作的总次数m等于3。那么i依次为1,2。在采集到3帧投影数据时,低剂量成像装置对3帧投影数据进行解析重建处理,生成解析重建图像J。当采集到6帧投影数据时,低剂量成像装置对6帧投影数据进行第1次迭代重建操作,生成第1次迭代重建图像D1,将迭代重建图像D1与解析重建图像J进行图像融合处理,生成第1次融合图像B1并显示第1次融合图像B1;低剂量成像装置对6帧投影数据进行第2次迭代重建操作,生成第2次迭代重建图像D2,并将第2次迭代重建图像D2与解析重建图像J进行图像融合处理,生成第2次融合图像B2并显示第2次融合图像B2。低剂量成像装置再对6帧投影数据进行第3次迭代重建操作,生成第3次迭代重建图像D3并显示第3次迭代重建图像D3。这样一来,低剂量成像装置共显示2个融合图像:B1和B2,1个迭代重建图像:D3。
其中,在一种可实现方式中,步骤1034中显示第i次融合图像,可以包括:
显示第i次融合图像以及距生成第m次迭代重建图像的进度信息。
示例的,距生成第m次迭代重建图像的进度信息可以是从第i次融合图像的生成时刻开始到第m次迭代重建图像的生成时刻的时长。
假设预设数据量为10帧,迭代重建处理中迭代重建操作的总次数m等于3那么,i依次为1,2。低剂量成像装置在采集到10帧投影数据时,生成第1 次融合图像,并显示第1次融合图像,以及显示第1次融合图像以及距生成第3次迭代重建图像的进度信息;低剂量成像装置生成第2次融合图像,并显示第2次融合图像,以及显示第2次融合图像以及距生成第3次迭代重建图像的进度信息。
在本发明实施例中,低剂量成像装置显示第i次融合图像以及距生成第m次迭代重建图像的进度信息,使得工作人员能够及时了解融合图像的生成进度,便于执行对应的临床治疗任务。
可选的,如图6所示,步骤1034中将第i次迭代重建图像与第一图像进行图像融合处理,生成第i次融合图像,可以包括:
步骤1034a、基于迭代重建操作的次数i,确定第i次迭代重建图像中像素值(或像素)的权重以及第一图像中像素值的权重。
步骤1034b、根据第i次迭代重建图像中像素值的权重和第一图像中像素值的权重,将第i次迭代重建图像和第一图像进行融合处理,生成第i次融合图像。
其中,第i次迭代重建图像中的一像素的权重与第一图像中与该像素融合的像素的权重之和为1,另外,第i次迭代重建图像中像素值的权重与迭代重建操作的次数i的大小正相关,第一图像中像素值的权重与迭代重建操作的次数i的大小负相关。也即是,迭代重建操作的次数i越大,第i次迭代重建图像中像素值的权重越大,而第一图像中像素值的权重越小。通常,图像中的像素可以包括软组织对应的像素和骨组织对应的像素,但图像中的软组织对应的像素值的权重与骨组织对应的像素值的权重增大或者减小的速率和/或模式可以不一样。
在本发明实施例中,当迭代重建操作的次数i增大时,低剂量成像装置可以增大对应的第i次迭代重建图像中像素值的权重,包括软组织对应的像素的像素值的权重和骨组织对应的像素的像素值的权重,但骨组织与软组织对应的像素的像素值的权重增大的速率和/或模式不一样,这里,尤其是增大骨组织对应的像素的像素值的权重;同时减小第一图像中像素值的权重,包括软组织对应的像素的像素值的权重和骨组织对应的像素的像素值的权重,但骨组织与软组织对应的像素的像素值的权重减小的速率和/或模式不一样。当迭代重建操作的次数i增大至预设值时,第i次融合图像中的像素值等于第i次迭代重建图像中的像素值,示例的,预设值可以为m-1。
示例的,当迭代重建操作的次数i等于1时,第1次迭代重建图像中像素值的权重为q1,第一图像中像素值的权重为p1;当迭代重建操作的次数i等于2时,第2次迭代重建图像中像素值的权重为q2,第一图像中像素值的权重为p2,则有q1<q2,p1>p2。
可选的,步骤103中显示第二图像,可以包括:
显示第二图像以及第二图像的质量信息,该质量信息用于表征第二图像相较于第一图像的图像质量。
在本发明实施例中,低剂量成像装置除了可以显示第二图像外,还可以采用图像质量评价方式确定第二图像的质量信息,并显示第二图像的质量信息,这样一来,工作人员能够及时了解第二图像的质量,进而确定该第二图像可以用于临床治疗中的何种临床治疗任务,避免了工作人员将时间都浪费在了无谓的等待上面。
可选的,图像质量评价方式可以为数字图像质量客观评价方式。示例的,数字图像质量客观评价方式的类型可以为全参考类型(Full Reference,FR)、部分参考类型(Reduced Reference,RR)或无参考类型(No Reference,NR)。其中,FR指的是原始图像是已知的,基于原始图像来评价当前图像的质量。NR指的是在没有原始图像的情况下,基于当前图像中具有辨别力的图像局部特征预测整体图像的质量。RR介于FR和NR之间,其指的是利用原始图像的部分信息来评价当前图像的质量。其中,原始图像即为本发明实施例中的第一图像,当前图像即为本发明实施例中的第二图像。示例的,第一图像可以为解析重建图像,第二图像可以包括迭代重建图像,第二图像还可以包括融合图像。
此外,图像质量评价方式还可以为主观实验评价方式。主观实验评价方式指的是在一定的条件(图像源、显示器和观看条件等)下,给观看者同时提供两幅图像,该两幅图像即为本发明实施例中的第二图像和第一图像。观看者基于第二图像和第一图像得到大量的评分数据,并对大量的评分数据进行统计,进而得到第二图像的质量信息。示例的,评分数据可以包括均值、标准差等数据。其中,在主观实验评价方式中,第二图像的质量信息可以有两种表达形式,一种表达形式为绝对评分表达形式,即表示第二图像的绝对质量;另一种表达形式为差值表达形式,即表示第二图像和第一图像的评价成绩的绝对差。
需要说明的是,本发明实施例提供的低剂量成像方法的步骤的先后顺序可以进行适当调整,该低剂量成像方法的步骤也可以根据情况进行相应增减, 任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化的方法,都应涵盖在本发明的保护范围之内,因此不再赘述。
综上所述,本发明实施例提供的低剂量成像方法,能够持续采集投影数据,在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示第一图像,并在采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示第二图像,且还能够显示第二图像的进度信息,以及第二图像的质量信息,相较于相关技术,无需等待所有帧的投影数据被采集完后显示图像,缩短了显示图像所用的时间,为临床治疗提供参考数据。
本发明实施例提供了一种低剂量成像装置,该低剂量成像装置具有显示功能,设置于成像系统中,如图7所示,该装置700包括:
采集模块710,用于持续采集投影数据。
第一处理模块720,用于在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示第一图像。
第二处理模块730,用于在采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示第二图像。
综上所述,本发明实施例提供的低剂量成像装置,采集模块持续采集投影数据,第一处理模块在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示第一图像,第二处理模块在采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示第二图像,相较于相关技术,无需等待所有帧的投影数据被采集完后显示图像,缩短了显示图像所用的时间。可选的,第一处理模块720,用于:
对采集到的投影数据进行解析重建处理,以生成解析重建图像;
显示解析重建图像。
可选的,第二处理模块730,用于:
对预设数量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像;
显示第i次迭代重建图像;
其中,i为1,2,…,m,m为迭代重建处理中迭代重建操作的总次数, 且m为大于等于1的整数。
可选的,第二处理模块730,用于:
显示第i次迭代重建图像以及距生成第m次迭代重建图像的进度信息。
可选的,第二处理模块730,用于:
对预设数据量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像;
将第i次迭代重建图像与第一图像进行图像融合处理,生成第i次融合图像并显示第i次融合图像,其中,i为1,2,…,m-1,m为迭代重建处理中迭代重建操作的总次数,且m为大于等于2的整数;
对预设数据量的投影数据进行迭代重建处理中的第m次迭代重建操作,生成第m次迭代重建图像并显示第m次迭代重建图像。
可选的,第二处理模块730,用于:
基于迭代重建操作的次数i,确定第i次迭代重建图像中像素值的权重以及第一图像中像素值的权重;
根据第i次迭代重建图像中像素值的权重和第一图像中像素值的权重,将第i次迭代重建图像和第一图像进行融合处理,生成第i次融合图像。
可选的,第二处理模块730,用于:
显示第i次融合图像以及距生成第m次迭代重建图像的进度信息。
可选的,第二处理模块730,用于:
显示第二图像以及第二图像的质量信息,该质量信息用于表征第二图像相较于第一图像的图像质量。综上所述,本发明实施例提供的低剂量成像装置,采集模块持续采集投影数据,第一处理模块在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示第一图像,第二处理模块在采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示第二图像,相较于相关技术,无需等待所有帧的投影数据被采集完后显示图像,缩短了显示图像所用的时间。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的装置和模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
本发明实施例还提供了一种低剂量成像装置,如图8所示,包括:
包括存储器801、处理器802以及存储在所述存储器801中并可在所述处理器802上运行的计算机程序8011,所述处理器802执行所述计算机程序8011时实现上述实施例所提供的低剂量成像方法的步骤。
综上所述,本发明实施例提供的低剂量成像装置,能够持续采集投影数据,在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示第一图像,并在采集到的投影数据的数据量达到预设数据量时,基于预设数量的投影数据进行第二处理,生成第二图像并显示第二图像,且还能够显示第二图像的进度信息,以及第二图像的质量信息,相较于相关技术,无需等待所有帧的投影数据被采集完后显示图像,缩短了显示图像所用的时间,为临床治疗提供参考数据。
本发明实施例还提供了一种计算机可读存储介质,该存储介质为非易失性可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现上述实施例所提供的低剂量成像方法的步骤。
本发明实施例还提供了一种计算机程序产品,所述计算机程序产品中存储有指令,当其在计算机上运行时,使得计算机实现上述实施例所提供的低剂量成像方法的步骤。
本发明实施例还提供了一种芯片,所述芯片包括可编程逻辑电路和/或程序指令,当所述芯片运行时用于实现上述实施例所提供的低剂量成像方法的步骤。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
- 一种低剂量成像方法,其特征在于,所述方法包括:持续采集投影数据;在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示所述第一图像;当采集到的投影数据的数据量达到预设数据量时,基于所述预设数量的投影数据进行第二处理,生成第二图像并显示所述第二图像。
- 根据权利要求1所述的方法,其特征在于,所述基于所述预设数量的投影数据进行第二处理,生成第二图像并显示所述第二图像,包括:对所述预设数量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像;显示所述第i次迭代重建图像;其中,所述i为1,2,…,m,所述m为所述迭代重建处理中迭代重建操作的总次数,且所述m为大于等于1的整数。
- 根据权利要求2所述的方法,其特征在于,所述显示所述第i次迭代重建图像,包括:显示所述第i次迭代重建图像以及距生成所述第m次迭代重建图像的进度信息。
- 根据权利要求1所述的方法,其特征在于,所述基于所述预设数量的投影数据进行第二处理,生成第二图像并显示所述第二图像,包括:对所述预设数据量的投影数据进行迭代重建处理中的第i次迭代重建操作,生成第i次迭代重建图像;将所述第i次迭代重建图像与所述第一图像进行图像融合处理,生成所述第i次融合图像并显示所述第i次融合图像,其中,所述i为1,2,…,m-1,所述m为所述迭代重建处理中迭代重建操作的总次数,且所述m为大于等于2的整数;对所述预设数据量的投影数据进行迭代重建处理中的第m次迭代重建操 作,生成第m次迭代重建图像并显示第m次迭代重建图像。
- 根据权利要求4所述的方法,其特征在于,将所述第i次迭代重建图像与所述第一图像进行图像融合处理,生成所述第i次融合图像,包括:基于所述迭代重建操作的次数i,确定所述第i次迭代重建图像中像素值的权重以及所述第一图像中像素值的权重;根据所述第i次迭代重建图像中像素值的权重和所述第一图像中像素值的权重,将所述第i次迭代重建图像和所述第一图像进行融合处理,生成所述第i次融合图像。
- 根据权利要求5所述的方法,其特征在于,所述显示所述第i次融合图像,包括:显示所述第i次融合图像以及距生成所述第m次迭代重建图像的进度信息。
- 根据权利要求1所述的方法,其特征在于,所述显示所述第二图像,包括:显示所述第二图像以及所述第二图像的质量信息,所述质量信息用于表征所述第二图像相较于所述第一图像的图像质量。
- 一种低剂量成像装置,其特征在于,所述装置包括:采集模块,用于持续采集投影数据;第一处理模块,用于在采集到的投影数据的数据量达到预设数据量之前,对采集到的投影数据进行第一处理,生成第一图像并显示所述第一图像;第二处理模块,用于在采集到的投影数据的数据量达到预设数据量时,基于所述预设数量的投影数据进行第二处理,生成第二图像并显示所述第二图像。
- 一种低剂量成像装置,包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至7任一项所述方法的步骤。
- 一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程 序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至7任一项所述方法的步骤。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201880002818.3A CN110730977B (zh) | 2018-05-04 | 2018-05-04 | 低剂量成像方法及装置 |
PCT/CN2018/085628 WO2019047545A1 (zh) | 2018-05-04 | 2018-05-04 | 低剂量成像方法及装置 |
US17/053,061 US20210233293A1 (en) | 2018-05-04 | 2018-05-04 | Low-dose imaging method and apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2018/085628 WO2019047545A1 (zh) | 2018-05-04 | 2018-05-04 | 低剂量成像方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2019047545A1 true WO2019047545A1 (zh) | 2019-03-14 |
Family
ID=65634749
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2018/085628 WO2019047545A1 (zh) | 2018-05-04 | 2018-05-04 | 低剂量成像方法及装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20210233293A1 (zh) |
CN (1) | CN110730977B (zh) |
WO (1) | WO2019047545A1 (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115778412B (zh) * | 2023-02-09 | 2023-04-28 | 之江实验室 | X光声成像中造影剂剂量的优化方法、装置及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103180879A (zh) * | 2010-10-26 | 2013-06-26 | 皇家飞利浦电子股份有限公司 | 用于从投影数据对对象进行混合重建的设备和方法 |
CN103810735A (zh) * | 2014-02-28 | 2014-05-21 | 南方医科大学 | 一种低剂量x射线ct图像统计迭代重建方法 |
CN103810734A (zh) * | 2014-02-28 | 2014-05-21 | 南方医科大学 | 一种低剂量x射线ct投影数据恢复方法 |
US20160210762A1 (en) * | 2015-01-15 | 2016-07-21 | General Electric Company | Method for tomographic reconstruction |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3006722B2 (ja) * | 1991-02-22 | 2000-02-07 | 株式会社東芝 | コンピュータ断層撮影装置 |
DE19812118A1 (de) * | 1998-03-19 | 1999-09-23 | Coin Inst Fuer Computer Integr | Verfahren zur Erhöhung der Auflösung eines Bildes |
EP1802235A1 (en) * | 2004-10-11 | 2007-07-04 | Philips Intellectual Property & Standards GmbH | Imaging system for the generation of high-quality x-ray projections |
US10918305B2 (en) * | 2010-12-21 | 2021-02-16 | Klinikum Mannheim Gmbh Universitätsklinikum Medizinische | Method and system for 4D radiological intervention guidance (4D-cath) |
CN102334979B (zh) * | 2011-08-03 | 2012-12-12 | 中国科学院自动化研究所 | 一种基于迭代收缩的双模态融合断层成像方法 |
CN103366389A (zh) * | 2013-04-27 | 2013-10-23 | 中国人民解放军北京军区总医院 | Ct图像重建方法 |
CN104021582A (zh) * | 2014-05-28 | 2014-09-03 | 山东大学 | 一种ct迭代图像重建方法 |
CN106373165B (zh) * | 2016-08-31 | 2017-11-28 | 广州华端科技有限公司 | 断层合成图像重建方法和系统 |
-
2018
- 2018-05-04 WO PCT/CN2018/085628 patent/WO2019047545A1/zh active Application Filing
- 2018-05-04 US US17/053,061 patent/US20210233293A1/en not_active Abandoned
- 2018-05-04 CN CN201880002818.3A patent/CN110730977B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103180879A (zh) * | 2010-10-26 | 2013-06-26 | 皇家飞利浦电子股份有限公司 | 用于从投影数据对对象进行混合重建的设备和方法 |
CN103810735A (zh) * | 2014-02-28 | 2014-05-21 | 南方医科大学 | 一种低剂量x射线ct图像统计迭代重建方法 |
CN103810734A (zh) * | 2014-02-28 | 2014-05-21 | 南方医科大学 | 一种低剂量x射线ct投影数据恢复方法 |
US20160210762A1 (en) * | 2015-01-15 | 2016-07-21 | General Electric Company | Method for tomographic reconstruction |
Also Published As
Publication number | Publication date |
---|---|
CN110730977B (zh) | 2024-03-29 |
CN110730977A (zh) | 2020-01-24 |
US20210233293A1 (en) | 2021-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6999576B2 (ja) | 空間とスペクトル情報に基づく複数エネルギーのct画像におけるノイズ制御のためのシステムと方法 | |
US8107695B2 (en) | Methods and systems for assessing patient movement in diagnostic imaging | |
McCollough et al. | Achieving routine submillisievert CT scanning: report from the summit on management of radiation dose in CT | |
US9002089B2 (en) | Image data registration for dynamic perfusion CT | |
JP6442243B2 (ja) | コンピュータ断層撮影画像再構成のアーチファクト低減の方法および装置 | |
US12073492B2 (en) | Method and system for generating attenuation map from SPECT emission data | |
JP6662880B2 (ja) | 放射線放出撮像システム、記憶媒体及び撮像方法 | |
JP2010500901A (ja) | エネルギー感受性コンピュータ断層撮影における動き補償 | |
Allie et al. | Pitfalls and artifacts using the D-SPECT dedicated cardiac camera | |
JP6416572B2 (ja) | コンピュータ断層撮影(ct)画像再構成におけるアーチファクトを低減するための方法および装置 | |
CN113116370A (zh) | 用于正电子发射断层摄影中的运动检测的方法和系统 | |
US20200226800A1 (en) | Tomographic imaging apparatus and method of generating tomographic image | |
KR101946576B1 (ko) | 의료 영상 장치 및 의료 영상 처리 방법 | |
US12070348B2 (en) | Methods and systems for computed tomography | |
KR20190079371A (ko) | 조영제를 주입하여 컴퓨터 단층 촬영하는 방법 및 장치 | |
KR102399792B1 (ko) | 인공지능 기반 Hounsfield unit (HU) 정규화 및 디노이징 (Denoising)을 통한 전처리 장치 및 방법 | |
WO2024174984A1 (zh) | 能谱ct成像的方法、装置和ct扫描成像系统 | |
JP4997575B2 (ja) | 減弱係数マップ作成装置、減弱係数マップ作成方法およびプログラム | |
US7569828B2 (en) | Application-driven optimization of acquisition and reconstruction of SPECT/PET projection data | |
US20030048937A1 (en) | Method of processing visual imagery from a medical imaging device | |
CN110215226B (zh) | 图像衰减校正方法、装置、计算机设备和存储介质 | |
JP4879472B2 (ja) | 脳血流定量解析プログラム、記録媒体および脳血流画像データ処理方法 | |
WO2019047545A1 (zh) | 低剂量成像方法及装置 | |
EP4148680A1 (en) | Attenuation correction-based weighting for tomographic inconsistency detection | |
CN110811665A (zh) | Pet图像衰减校正方法、装置、计算机设备和存储介质 |
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: 18853513 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: 18853513 Country of ref document: EP Kind code of ref document: A1 |