IP Library › Granted Patent US 12,614,339
Granted Patent B2
US 12,614,339 · App. 19/023,894 · Granted Apr 28, 2026

Three-dimensional CT imaging method and apparatus

Inventors: Ping Chen (Shanxi, CN); Zhiqing Wei (Shanxi, CN); Sukai Wang (Shanxi, CN); Jiaotong Wei (Shanxi, CN); Xiaojie Zhao (Shanxi, CN); Yangxu Wu (Shanxi, CN)
Assignee: NORTH UNIVERSITY OF CHINA
G06T15/08G06T3/40G06V10/44G16H30/20G16H30/40
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,614,339
App. No.
19/023,894
Granted
Apr 28, 2026
Kind
B2
Abstract

The disclosure provides a three-dimensional CT imaging method and apparatus. The method includes: collecting two DR images with a perpendicular relationship for an object to be imaged; inputting the two DR images into a preset three-dimensional volume reconstruction model to obtain a three-dimensional volume of the object outputted by the three-dimensional volume reconstruction model; and slicing the three-dimensional volume to obtain three-dimensional CT imaging of the object, wherein the preset three-dimensional volume reconstruction model is trained using multiple training samples, each training sample comprising: a DR image group composed of two DR images with a perpendicular relationship, a pseudo-sinogram corresponding to the DR image group, and a three-dimensional volume corresponding to the DR image group.

Claims (86)

1 . A three-dimensional computed tomography (CT) imaging method based on X-ray dual projection, comprising:

collecting two digital radiography (DR) images with a perpendicular relationship for an object to be imaged;

obtain a three-dimensional volume of the object to be imaged corresponding to the two DR images with a perpendicular relationship based on a preset three-dimensional volume reconstruction model; and

slicing the three-dimensional volume to obtain three-dimensional CT imaging of the object to be imaged,

wherein, generating the preset three-dimensional volume reconstruction model comprises:

obtaining a training data set; wherein the training data set is composed of DR images of reference objects collected circumferentially;

training to-be-trained three-dimensional volume reconstruction model using DR images with a perpendicular relationship in the training data set to obtain the preset three-dimensional volume reconstruction model;

wherein, when obtaining the training data set, the method further comprises:

obtaining a reference three-dimensional volume and a reference pseudo-sinogram of each of the reference objects based on the training data set; wherein the reference pseudo-sinogram is obtained by stacking collected DR images of a same reference object into a three-dimensional matrix according to a collection order, summing sinograms in each layer column by column, and normalizing a sum result;

dividing the DR images of each reference object in the training data set into a plurality of DR image groups, and each of the DR image groups comprises two DR images with a perpendicular relationship;

wherein the training to-be-trained three-dimensional volume reconstruction model using DR images with a perpendicular relationship in the training data set to obtain the preset three-dimensional volume reconstruction model comprises:

inputting the plurality of DR image groups into the to-be-trained three-dimensional volume reconstruction model, training the to-be-trained three-dimensional volume reconstruction model, and respectively obtaining a training pseudo-sinogram and a training three-dimensional volume;

calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram and the reference three-dimensional volume according to a training order;

terminating the training and obtaining the preset three-dimensional volume reconstruction model in a case that any calculated loss value is less than a target loss value, or the number of training times reaches a target number of times.

2 . The method according to claim 1 wherein the to-be-trained three-dimensional 1 volume reconstruction model comprises an enhancement network model and a reconstruction network model;

the enhancement network model is a model based on an encoder-decoder network architecture, and configured to input a DR image group with a perpendicular relationship and output a corresponding training pseudo-sinogram; and

the reconstruction network model is an autoencoder network model for two-dimensional to three-dimensional volume mapping and configured to input the training pseudo-sinogram and output a corresponding training three-dimensional volume.

3 . The method according to claim 2 , wherein the enhancement network model to input a DR image group with a perpendicular relationship, and output a corresponding training pseudo-sinogram, comprises:

performing pattern and feature

extraction on the inputted DR image group, and performing a cascade operation on a results of the extraction;

performing latent space feature encoding on features after the cascade operation to obtain encoded latent space features; and

decoding the encoded latent space features to obtain the training pseudo-sinogram.

4 . The method according to claim 2 , wherein the reconstruction network model to input the training pseudo-sinogram, and output a corresponding training three-dimensional volume, comprises:

performing feature capture and enhancement on the input training pseudo-sinogram;

encoding a result of feature enhancement, and continuously reducing the spatial resolution on the basis of keeping the number of channels unchanged; and

decoding a result of the encoding, and continuously restoring the spatial resolution to a spatial resolution of the reference three-dimensional volume to obtain the training three-dimensional volume.

5 . The method according to claim 1 , wherein the calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume according to the training order, comprises:

calculating the loss value using a preset loss function based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume,

wherein the preset loss function is a function of a weighted summation of a first loss function and a second loss function; the first loss function is configured to calculate a mean square error loss between a training pseudo-sinogram and a reference pseudo-sinogram, and the second loss function is configured to calculate a mean square error loss of a training three-dimensional volume and a reference three-dimensional volume.

6 . The method according to claim 2 , wherein the calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume according to the training order, comprises:

calculating the loss value using a preset loss function based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume,

wherein the preset loss function is a function of a weighted summation of a first loss function and a second loss function; the first loss function is configured to calculate a mean square error loss between a training pseudo-sinogram and a reference pseudo-sinogram, and the second loss function is configured to calculate a mean square error loss of a training three-dimensional volume and a reference three-dimensional volume.

7 . The method according to claim 3 , wherein the calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume according to the training order, comprises;

calculating the loss value using a preset loss function based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume,

wherein the preset loss function is a function of a weighted summation of a first loss function and a second loss function; the first loss function is configured to calculate a mean square error loss between a training pseudo-sinogram and a reference pseudo-sinogram, and the second loss function is configured to calculate a mean square error loss of a training three-dimensional volume and a reference three-dimensional volume.

8 . The method according to claim 4 , wherein the calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume according to the training order comprises:

calculating the loss value using a preset loss function based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume,

wherein the preset loss function is a function of a weighted summation of a first loss function and a second loss function: the first loss function is configured to calculate a mean square error loss between a training pseudo-sinogram and a reference pseudo-sinogram, and the second loss function is configured to calculate a mean square error loss of a training three-dimensional volume and a reference three-dimensional volume.

9 . A three-dimensional computer tomography (CT) imaging device based on X-ray dual projection, comprising:

a storage unit, configured to store a preset three-dimensional volume reconstruction model; generating the preset three-dimensional volume reconstruction model comprises: obtaining a training data set; wherein the training data set is composed of digital radiography (DR) images of reference objects collected circumferentially; training to-be-trained three-dimensional volume reconstruction model using DR images with a perpendicular relationship in the training data set to obtain the preset three-dimensional volume reconstruction model;

an acquisition unit, configured to collect two DR images with a perpendicular relationship for an object to be imaged;

an obtaining unit, configured to obtain a three-dimensional volume of the object to be imaged corresponding to the two DR images with a perpendicular relationship based on the preset three-dimensional volume reconstruction model; and

an imaging unit, configured to slice the three-dimensional volumes to obtain three-dimensional CT imaging of the object to be imaged,

wherein, generating the three-dimensional volume reconstruction model stored in the storage unit comprises:

when obtaining the training data set, obtaining a reference three-dimensional volume and a reference pseudo-sinogram of each of the reference objects based on the training data set; wherein the reference pseudo-sinogram is obtained by stacking collected DR images of a same reference object into a three-dimensional matrix according to a collection order, summing sinograms in each layer column by column, and normalizing a sum result;

dividing the DR images of each reference object in the training data set into a plurality of DR image groups, and each of the DR image groups comprises two DR images with a perpendicular relationship;

inputting the plurality of DR image groups into the to-be-trained three-dimensional volume reconstruction model, training the to-be-trained three-dimensional volume reconstruction model, and respectively obtaining a training pseudo-sinogram and a training three-dimensional volume;

calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram and the reference three-dimensional volume according to a training sequence;

terminating the training and obtaining the preset three-dimensional volume reconstruction model in a case that any calculated loss value is less than a target loss value, or the number of training times reaches a target number of times.

10 . An electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable by the processor, the processor is configured to execute the computer program to implement a three-dimensional computed tomography (CT) imaging method based on X-ray dual projection, the method comprises:

collecting two digital radiography (DR) images with a perpendicular relationship for an object to be imaged;

obtain a three-dimensional volume of the object to be imaged corresponding to the two DR images with a perpendicular relationship based on a preset three-dimensional volume reconstruction model; and

slicing the three-dimensional volume to obtain three-dimensional CT imaging of the object to be imaged,

wherein, generating the preset three-dimensional volume reconstruction model comprises:

obtaining a training data set; wherein the training data set is composed of DR images of reference objects collected circumferentially;

training to-be-trained three-dimensional volume reconstruction model using DR images with a perpendicular relationship in the training data set to obtain the preset three-dimensional volume reconstruction model;

wherein, when obtaining the training data set, the method further comprises:

obtaining a reference three-dimensional volume and a reference pseudo-sinogram of each of the reference objects based on the training data set; wherein the reference pseudo-sinogram is obtained by stacking collected DR images of a same reference object into a three-dimensional matrix according to a collection order, summing sinograms in each layer column by column, and normalizing a sum result;

dividing the DR images of each reference object in the training data set into a plurality of DR image groups, and each of the DR image groups comprises two DR images with a perpendicular relationship;

wherein the training to-be-trained three-dimensional volume reconstruction model using DR images with a perpendicular relationship in the training data set to obtain the preset three-dimensional volume reconstruction model comprises:

inputting the plurality of DR image groups into the to-be-trained three-dimensional volume reconstruction model, training the to-be-trained three-dimensional volume reconstruction model, and respectively obtaining a training pseudo-sinogram and a training three-dimensional volume;

calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram and the reference three-dimensional volume according to a training order;

terminating the training and obtaining the preset three-dimensional volume reconstruction model in a case that any calculated loss value is less than a target loss value, or the number of training times reaches a target number of times.

11 . The electronic device according to claim 10 , wherein the to-be-trained three-dimensional volume reconstruction model comprises an enhancement network model and a reconstruction network model;

the enhancement network model is a model based on an encoder-decoder network architecture, and configured to input a DR image group with a perpendicular relationship and output a corresponding training pseudo-sinogram; and

the reconstruction network model is an autoencoder network model for two-dimensional to three-dimensional volume mapping, and configured to input the training pseudo-sinogram and output a corresponding training three-dimensional volume.

12 . The electronic device according to claim 11 , wherein the enhancement network model to input a DR image group with a perpendicular relationship, and output a corresponding training pseudo-sinogram, comprises:

performing pattern and feature extraction on the input DR image group, and performing a cascade operation on a result of the extraction;

performing latent space feature encoding on features after the cascade operation to obtain encoded latent space features; and

decoding the encoded latent space features to obtain the training pseudo-sinogram.

13 . The electronic device according to claim 11 , wherein the reconstruction network model to input the training pseudo-sinogram and output a corresponding training three-dimensional volume, comprises:

performing feature capture and enhancement on the input training pseudo-sinogram;

encoding a result of feature enhancement, and continuously reducing the spatial resolution on the basis of keeping the number of channels unchanged; and

decoding a result of the encoding, and continuously restoring the spatial resolution to a spatial resolution of the reference three-dimensional volume to obtain the training three-dimensional volume.

14 . The electronic device according to claim 10 , wherein the calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume according to the training order, comprises:

calculating the loss value using a preset loss function based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume,

wherein the preset loss function is a function of a weighted summation of a first loss function and a second loss function; the first loss function is configured to calculate a mean square error loss between a training pseudo-sinogram and a reference pseudo-sinogram, and the second loss function is configured to calculate a mean square error loss of a training three-dimensional volume and a reference three-dimensional volume.

15 . The electronic device according to claim 11 , wherein the calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume according to the training order, comprises;

calculating the loss value using a preset loss function based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume,

wherein the preset loss function is a function of a weighted summation of a first loss function and a second loss function; the first loss function is configured to calculate a mean square error loss between a training pseudo-sinogram and a reference pseudo-sinogram, and the second loss function is configured to calculate a mean square error loss of a training three-dimensional volume and a reference three-dimensional volume.

16 . The electronic device according to claim 12 , wherein the calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume according to the training order, comprises:

calculating the loss value using a preset loss function based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume,

wherein the preset loss function is a function of a weighted summation of a first loss function and a second loss function; the first loss function is configured to calculate a mean square error loss between a training pseudo-sinogram and a reference pseudo-sinogram, and the second loss function is configured to calculate a mean square error loss of a training three-dimensional volume and a reference three-dimensional volume.

17 . The electronic device according to claim 13 , wherein the calculating a loss value based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume according to the training order, comprises:

calculating the loss value using a preset loss function based on the training pseudo-sinogram, the training three-dimensional volume, the reference pseudo-sinogram, and the reference three-dimensional volume,

wherein the preset loss function is a function of a weighted summation of a first loss function and a second loss function; the first loss function is configured to calculate a mean square error loss between a training pseudo-sinogram and a reference pseudo-sinogram, and the second loss function is configured to calculate a mean square error loss of a training three-dimensional volume and a reference three-dimensional volume.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jan 21, 2025
From: CHEN, PING; WEI, ZHIQING; WANG, SUKAI; WEI, JIAOTONG; ZHAO, XIAOJIE; WU, YANGXU
To: NORTH UNIVERSITY OF CHINA
Reel/Frame 069942/0084 →
Priority Claims (1)
CN 202410831945.1 · Jun 25, 2024 · national
Continuity (1)
Related Publication 20250391100A1 · Dec 25, 2025
References Cited (17)
US 12016717B2 · Guo et al. · 2024 [cited by applicant]
US 20200305806A1 · Tang et al. · 2020 [cited by applicant]
US 20210074036A1 · Fuchs · 2021 [cited by examiner]
US 20210192810A1 · Paysan · 2021 [cited by examiner]
US 20220405990A1 · Chen · 2022 [cited by examiner]
US 20230011759A1 · Lu · 2023 [cited by examiner]
US 20230326596A1 · Lu et al. · 2023 [cited by applicant]
US 20250157033A1 · Golden · 2025 [cited by examiner]
CN 109745062A · 2019 [cited by applicant]
CN 114004929A · 2022 [cited by applicant]
CN 119784931A · 2025 [cited by applicant]
WO WO2022127496A1 · 2022 [cited by applicant]
Xueqin Sun et al., Reconstruction method suitable for fast CT imaging, Optics Express, published May 6, 2024, Taiyuan, CN (16 pages). [cited by applicant]
Niu Xiaowei, Research on CT image reconstruction and artifact correction algorithm, published Jan. 15, 2024 (61 pages). [cited by applicant]
Chinese Office Action for Chinese Application No. 202410831945.1 dated May 6, 2025 (32 pages). [cited by applicant]
Danil Kazimirov et al., Adaptive automated sinogram normalization for ring artifacts suppression in CT, Optics Express, May 2024, pp. 17606-17643, vol. 32, No. 10. [cited by applicant]
Qian Yufei, in-depth Learning Based Sparse Angle CT Reconstruction Study, Chinese Excellent Master's Theses Journal of Medical and Health Science and Technology, Feb. 15, 2021, E060-356, No. 02. [cited by applicant]