CN103714560B - A kind of image rebuilding method based on Katsevich algorithm and system - Google Patents

A kind of image rebuilding method based on Katsevich algorithm and system Download PDF

Info

Publication number
CN103714560B
CN103714560B CN201310734168.0A CN201310734168A CN103714560B CN 103714560 B CN103714560 B CN 103714560B CN 201310734168 A CN201310734168 A CN 201310734168A CN 103714560 B CN103714560 B CN 103714560B
Authority
CN
China
Prior art keywords
data
projection
image
interpolation
described step
Prior art date
Legal status (The legal status 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 status listed.)
Expired - Fee Related
Application number
CN201310734168.0A
Other languages
Chinese (zh)
Other versions
CN103714560A (en
Inventor
张岩
纪东梅
高伦
景丽
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Graduate School Harbin Institute of Technology
Original Assignee
Shenzhen Graduate School Harbin Institute of Technology
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 Shenzhen Graduate School Harbin Institute of Technology filed Critical Shenzhen Graduate School Harbin Institute of Technology
Priority to CN201310734168.0A priority Critical patent/CN103714560B/en
Publication of CN103714560A publication Critical patent/CN103714560A/en
Application granted granted Critical
Publication of CN103714560B publication Critical patent/CN103714560B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Apparatus For Radiation Diagnosis (AREA)
  • Image Processing (AREA)

Abstract

The present invention proposes a kind of image rebuilding method based on Katsevich algorithm, by the error of each operation module in Katsevich algorithm for reconstructing is optimized, the impact on rebuilding image of the error artifact can be effectively reduced, and the information of original true CT image can be kept.Prove through repeatedly simulation practice, the method of the present invention is simple to operate, can effectively reduce the error artifact in CT image, the reconstruction image obtained by the method for the present invention can not only reduce the metal artifacts impact on image, and can keep the detailed information of area-of-interest.

Description

A kind of image rebuilding method based on Katsevich algorithm and system
Technical field
The invention belongs to the image processing field of medical image, be specifically related to a kind of based on Katsevich algorithm pattern picture weight The method and system built.
Background technology
In the reconstruction image of exact reconstruction algorithm katsevich, it may appear that various difform artifacts, these artifacts Have impact on the accuracy of algorithm imaging.Wherein owing to the data for projection during rebuilding is carried out artifact produced by discretization just It it is the most important a kind of source of artifact.Directly result in that to change artifact Producing reason mainly Katsevich algorithm itself be a kind of Exact reconstruction algorithm, but in actual application or the data for projection that lacks is discrete, follow-up operation is all to dispersion number According to having carried out corresponding approximation operation, thus result in final imaging containing substantial amounts of artifact.High-quality in order to obtain CT rebuilds image, needs the error artifact that the discrete operations as far as possible reducing data for projection is brought, the method pretreatment generally used The method optimized, cardinal principle is to using different derivative algorithms that data for projection is carried out derivation module in image reconstruction process It is optimized, thus reaches to reduce the effect of artifact.
Rebuilding, from exact reconstruction algorithm Katsevich, the mechanism that image artifacts produces, multiple CT artifact pretreatment disappears Except method is suggested in succession.The method of pretreatment mainly includes using a kind of higher Method of Seeking Derivative of precision to carry out derivation etc..Use Single method carries out derivation, and the result being easily caused derivation can be lost, thus can produce new error, so weight to image Build effect improved the most obvious.
Summary of the invention
In order to solve the deficiencies in the prior art, the present invention provide a kind of simple to operate, reconstructed image quality is high and rebuilds Fireballing image rebuilding method based on Katsevich algorithm and system.
The present invention is achieved through the following technical solutions:
A kind of image rebuilding method based on Katsevich algorithm, comprises the following steps:
(1) data for projection rebuild required for image is obtained;
(2) described data for projection is carried out mixing derivation based on line-of-sight course and Richardson extrapolation, particularly as follows: to projection Each row, column of data, frame start most and the derivation of last at least two data for projection use line-of-sight course derivation;Right Data for projection in other uses Richardson extrapolation to carry out derivation, additionally carries out this module at fixed point in concrete operations Reason;
(3) data for projection processing described step (2) carries out length weight correction;
(4) interpolation operation before the data for projection processing described step (3) is filtered, sits data for projection The conversion of mark system, in this step, the married operation closest to interpolation Yu cubic convolution interpolation that data for projection is used, additionally exist Concrete operations carry out fixed point process to this module;
(5) data for projection processing described step (4) carries out Hilbert transform, then carries out adding respective window function Operation;
(6) data for projection processing described step (5) carries out interpolation, and data for projection carries out coordinate system change again Change, recover the coordinate system of data for projection, additionally in concrete operations, this module is carried out fixed point process;
(7) data for projection processing described step (6) carries out back projection's generation reconstruction image, in this step to throwing When shadow data carry out cumulative, the method using mixing difference to choose data for projection, additionally in concrete operations, this module is carried out Fixed point processes.
Further, when the data for projection obtained contains noise, described step (1) also includes obtained projection number According to being filtered process, to reduce noise in subsequent operation to rebuilding the impact of artifact in image.
Further, when described step (2) carries out derivation to data for projection, the method for employing is the method for mixing derivation, That is: data acquisition difference coefficient, line-of-sight course or the five-spot derivation to edge, for middle data acquisition center difference coefficient, Li Chasen Extrapolation or fitting process derivation.
Further, described step (5) uses fast Fourier transform FFT to complete Hilbert transform.Have only to into FFT computing of row, and operation result is temporarily stored for subsequent calculations.
Further, described step (6) uses interpolation method to be the hybrid interpolation closest to interpolation Yu cubic convolution interpolation Method.
Further, the method that the mixing difference in described step (7) chooses data for projection, particularly as follows: when coordinate (u*, When w*) being positioned at the outermost of detecting plate, closest difference is used to choose the data for projection participating in integration;When coordinate (u*, w*) position Except time on the intersection point of the rib in upper step and rib on detecting plate, use bilinear interpolation method and choose data for projection;When Coordinate (u*, w*) is positioned on detecting plate except, time on the rib in above-mentioned two steps, the method using closest interpolation chooses throwing Shadow data;When coordinate (u*, w*) is positioned on detecting plate the region except above-mentioned three steps, the method for closest interpolation is used to choose Data for projection.
Further, when described step (7) obtains integrating range [sb, st] required during back projection, need completing Solving of nonlinear equation, uses interative computation based on two way classification or Newton iteration method or Secant Method when solving.
Further, the method carrying out fixed point process in described step (2), (4), (6), (7), particularly as follows: use Cordic algorithm surmounts function linearisation in solving to implement, and realizes using in pipeline organization 14 for cordic algorithm Level flowing water progression, it is possible to meet the required precision of medical image, it addition, in fixed point, surface analysis devises based on affine The algorithm of form realizes;Precision analysis devises algorithm based on quantization error propagation model and simulation self adaptation annealing algorithm phase It is implemented in combination with;Meanwhile, according to different constraints, respectively obtaining range of error is 2-4, tri-groups of precision analysis of 2-8,2-13 Result.
Another aspect of the present invention, it is provided that a kind of image reconstruction side based on Katsevich algorithm performing the present invention The system of method, it is characterised in that: this system includes that CT machine, image reconstruction work station, high-definition image display device, the X of CT machine penetrate Detecting plate in ray tube, CT machine;Its working method is: the relative invariant position of the X-ray tube in CT machine and detecting plate, they Rotating along scanning gantry 360 ° and be irradiated reconstruction object, the data for projection obtained is by PCI-E bus on a frame-by-frame basis It is transferred in image reconstruction work station, image reconstruction work station utilizes the image based on Katsevich algorithm of the present invention Method for reconstructing carries out the most quickly rebuilding, and the reconstruction image obtained is shown by high-definition image display device.
The invention has the beneficial effects as follows: the present invention is by carrying out the error of modules in Katsevich algorithm for reconstructing Optimize, it is possible to be effectively reduced the impact on rebuilding image of the error artifact, correlation module has been carried out fixed point process simultaneously and has added The fast speed rebuild, and the information of original true CT image can be kept.Prove through repeatedly simulation practice, the side of the present invention Method is simple to operate, it is possible to effectively reduce the error artifact in CT image, by the reconstruction image of the method acquisition of the present invention not only The metal artifacts impact on image can be reduced, and the detailed information of area-of-interest can be kept.
Accompanying drawing explanation
Fig. 1 is the image rebuilding method flow chart based on Katsevich algorithm of the present invention;
Fig. 2 is the method flow diagram that fixed point processes;
Fig. 3 is the detailed annotation figure of mixing derivation operations;
Fig. 4 is the interpolation schematic diagram of individual frames data for projection;
The middle hybrid interpolation schematic diagram of Tu5Shi back projection operation;
Fig. 6 is carried out the system cross sectional representation in the Y direction of the inventive method;
Fig. 7 (a) is the Shepp-Logan body mould of standard, is used as the ideal body mould of CT image reconstruction;
Fig. 7 (b) is the image not carrying out the Katsevich method for reconstructing reconstruction that invention optimizes;
Fig. 7 (c) is by the image rebuild based on the rear Katsevich method for reconstructing that the present invention optimizes;
In figure, the implication of each sign is: 1-CT machine, 2-image reconstruction work station, 3-high-definition image display device, 4-CT machine X-ray bulb, detecting plate in 5-CT machine.
Detailed description of the invention
The present invention is further described for explanation and detailed description of the invention below in conjunction with the accompanying drawings.
Image rebuilding method based on Katsevich algorithm as shown in Figure 1, in turn includes the following steps:
(1) obtain CT machine systematic parameter and rebuild the data for projection required for image.When the data for projection obtained contains During noise, need to be filtered obtained data for projection processing, to reduce noise in subsequent operation to rebuilding in image The impact of artifact.
(2) data for projection obtaining step (1) carries out mixing derivation based on line-of-sight course and Richardson extrapolation, specifically Embodiment as shown in Figure 3, the later last at least two started most of each row, column of data for projection, frame is thrown The derivation of shadow data (MARG) uses line-of-sight course derivation;Data for projection for other uses Richardson extrapolation to ask Lead.The main purpose using line-of-sight course is that it can with the derivative obtaining MARG of degree of precision, thus be avoided throwing The most obtainable situation of shadow data derivative;And the purpose using Richardson extrapolation to carry out derivation be utilize Richardson extrapolation with The advantage of the change derivation precision change extrapolation number of times requirement as the case may be can use different extrapolation situations Reason looks into Senn process, meets the required precision of derivation to greatest extent.Be can be obtained by and original projection number by operation above According to the data for projection after the derivation that amount is equal, on the one hand improve the derivation precision of data to greatest extent, on the other hand eliminate After derivation, the reconstruction efficiency lost to improve algorithm for reconstructing of data, has carried out pinpointing the optimisation strategy of streamline to this module This module is optimized.
(3) data for projection processing step (2) carries out length weight correction;
(4) interpolation operation before the data for projection processing step (3) is filtered, carries out coordinate system to data for projection Conversion, in this step, the married operation closest to interpolation Yu cubic convolution interpolation using data for projection, such as accompanying drawing 4 institute Show, be i.e. closest to interpolation to edge data acquisition, and cubic convolution interpolation is used for non-edge data.Due to Precision closest to interpolation is the lowest so that rebuild image effect very poor, but simple to operate, bilinear interpolation can produce low pass The effect of wave filter and to make to rebuild image fuzzyyer, and cubic convolution interpolation can both improve precision and also be able to avoid simultaneously Obscuring of image, therefore improves the computational accuracy of interpolating module, and cubic convolution interpolation makes computing quantitative change big simultaneously.In order to carry The reconstruction efficiency of high algorithm for reconstructing, this module is optimized by the optimisation strategy that this module has carried out pinpointing streamline, can So that this module area of whole hardware circuit when carrying out hardware and realizing is less, faster, power consumption is lower for speed.
(5) data processing step (4) carry out Hilbert transform, and usual data for projection is to be come by convolution operation Completing Hilbert transform, in the present invention, the fast Fourier transform (FFT) of employing completes Hilbert transform, complete After having become Hilbert transform, the data for projection processed is added respective window function operation.
(6) data for projection processing step (5) carries out interpolation, and data for projection carries out coordinate system transformation again, Recover the coordinate system of data for projection, remain the married operation closest to interpolation Yu cubic convolution interpolation what this step used. Module is optimized by the streamline simultaneously using fixed point.
(7) data for projection producing step (6) carries out back projection's generation reconstruction image, in this step to data for projection When carrying out cumulative, the present invention uses and a kind of mixes the method that difference chooses data for projection, as shown in Figure 5: work as coordinate When (u*, w*) is positioned at the outermost of detecting plate, the present invention uses closest difference to choose the data for projection participating in integration;Work as coordinate When (u*, w*) is positioned on detecting plate on the intersection point except the rib in upper step and rib, uses bilinear interpolation method and choose throwing Shadow data;When coordinate (u*, w*) is positioned on detecting plate except time on the rib in above-mentioned two steps, use the side of closest interpolation Method chooses data for projection;When coordinate (u*, w*) is positioned on detecting plate the region except above-mentioned three steps, use closest interpolation Method chooses data for projection, so that precision is higher when choosing data for projection and being updated reconstruction point.In addition to add The fast speed rebuild, the present invention has carried out fixed point streamline optimization to this module, so that rebuild in hgher efficiency, power consumption is more Low.
When step (2) carries out derivation to the data for projection acquired in step (1), the method using mixing derivation: to edge Data can use difference coefficient, line-of-sight course or five-spot etc., can use, for middle data, the Method of Seeking Derivative that precision is higher, The such as method of center difference coefficient, Richardson extrapolation or fitting derivative so that the data volume after derivation and derivation it Before holding consistent, it is to avoid the loss of data and error therefore.
When step (5) carries out FFT computing to the data in step (4), owing to needing repeatedly to use the FFT of Hilbert core Operation result, in order to accelerate computing in process of reconstruction, it is only necessary to it is carried out once-through operation, and is deposited by its operation result temporarily Storage is for subsequent calculations.
Fixed point process method as shown in Figure 2, particularly as follows: use cordic algorithm solve to implement in surmount Function linearization, realizes using in pipeline organization 14 grades of flowing water progression for cordic algorithm, it is possible to meet medical image Required precision, it addition, in fixed point, surface analysis devises algorithm based on affine form and realizes;Precision analysis devises Algorithm based on quantization error propagation model and simulation self adaptation annealing algorithm combine realization;Meanwhile, according to different constraints Condition, respectively obtaining range of error is 2-4, the result of tri-groups of precision analysis of 2-8,2-13
As shown in Figure 6, being carried out the system of the image rebuilding method based on Katsevich algorithm of the present invention, this is System includes the detecting plate in CT machine, image reconstruction work station, high-definition image display device, the X-ray bulb of CT machine, CT machine;Its Working method is: the relative invariant position of the X-ray tube in CT machine and detecting plate, and they rotate counterweight along scanning gantry 360 ° Building object to be irradiated, the data for projection obtained on a frame-by-frame basis is transferred in image reconstruction work station by PCI-E bus, The image rebuilding method based on Katsevich algorithm utilizing the present invention in image reconstruction work station carries out the most quickly rebuilding, The reconstruction image obtained is shown by high-definition image display device.
As shown in Figure 7, Fig. 7-a is the Shepp-Logan body mould of standard, is used as the ideal body mould of CT image reconstruction, Fig. 7-b is the image not carrying out the Katsevich method for reconstructing reconstruction that invention optimizes, and Fig. 7-c is by based on optimization of the present invention Rear Katsevich method for reconstructing rebuild image, by the mean square deviation of image and rectangular histogram are carried out right, after optimization figure Significantly declining as mean square deviation has had, the most histogrammic high-frequency noise is also significantly improved.
Above content is to combine concrete preferred implementation further description made for the present invention, it is impossible to assert Being embodied as of the present invention is confined to these explanations.For general technical staff of the technical field of the invention, On the premise of present inventive concept, it is also possible to make some simple deduction or replace, all should be considered as belonging to the present invention's Protection domain.

Claims (10)

1. an image rebuilding method based on Katsevich algorithm, it is characterised in that said method comprising the steps of:
(1) data for projection rebuild required for image is obtained;
(2) described data for projection is carried out mixing derivation based on line-of-sight course and Richardson extrapolation, particularly as follows: to data for projection Each row, column, frame start most and the derivation of last at least two data for projection use line-of-sight course derivation;For it His data for projection uses Richardson extrapolation to carry out derivation, uses fixed point to process this step in concrete operations;
(3) data for projection processing described step (2) carries out length weight correction;
(4) interpolation operation before the data for projection processing described step (3) is filtered, carries out coordinate system to data for projection Conversion, in this step, the married operation closest to interpolation Yu cubic convolution interpolation using data for projection, in concrete operations In to this step use fixed point process;
(5) data for projection processing described step (4) carries out Hilbert transform, then carries out adding respective window function behaviour Make;
(6) data for projection processing described step (5) carries out interpolation, and data for projection carries out coordinate system transformation again, Recover the coordinate system of data for projection, use fixed point to process this step in concrete operations;
(7) data for projection processing described step (6) carries out back projection's generation reconstruction image, in this step to projection number According to when carrying out cumulative, use and mix the method that difference chooses data for projection, in concrete operations, this step is used at fixed point Reason.
Method the most according to claim 1, it is characterised in that: when the data for projection obtained contains noise, described step (1) also include being filtered obtained data for projection processing, to reduce noise in subsequent operation to rebuilding puppet in image The impact of shadow.
Method the most according to claim 1, it is characterised in that: when described step (2) carries out derivation to data for projection, use Method be the method for mixing derivation, it may be assumed that data acquisition difference coefficient, line-of-sight course or the five-spot derivation to edge, for middle Data acquisition center difference coefficient, Richardson extrapolation or fitting process derivation.
Method the most according to claim 1, it is characterised in that: described step (5) uses fast Fourier transform FFT to come Become Hilbert transform.
Method the most according to claim 4, it is characterised in that: have only to carry out a FFT computing, and operation result is faced Time storage for subsequent calculations.
Method the most according to claim 1, it is characterised in that: described step (6) uses interpolation method to be closest to interpolation Hybrid interpolation method with cubic convolution interpolation.
Method the most according to claim 1, it is characterised in that: the mixing difference in described step (7) chooses data for projection Method, particularly as follows: when coordinate (u*, w*) is positioned at the outermost of detecting plate, uses closest difference to choose and participate in integration Data for projection;When coordinate (u*, w*) is positioned on detecting plate on the intersection point except the rib in upper step and rib, use bilinearity Interpolation method chooses data for projection;When coordinate (u*, w*) is positioned on detecting plate except time on the rib in above-mentioned two steps, use The method of closest interpolation chooses data for projection;When coordinate (u*, w*) is positioned on detecting plate the region except above-mentioned three steps, adopt Data for projection is chosen by the method for closest interpolation.
Method the most according to claim 1, it is characterised in that: described step obtains integration required during back projection in (7) Time interval [sb, st], need to solve completing nonlinear equation, when solving, use interative computation based on two way classification or newton Iterative method or Secant Method.
Method the most according to claim 1, it is characterised in that: described step (2), (4), (6), (7) carry out fixed point The method processed, particularly as follows: use cordic algorithm to surmount function linearisation, for cordic algorithm in solving to implement Realize pipeline organization uses 14 grades of flowing water progression, it is possible to meet the required precision of medical image, it addition, in fixed point, Surface analysis devises algorithm based on affine form and realizes;Precision analysis devises algorithm based on quantization error propagation model Combine realization with simulation self adaptation annealing algorithm;Meanwhile, according to different constraints, respectively obtaining range of error is 2-4, The result of tri-groups of precision analysis of 2-8,2-13.
10. the image rebuilding method based on Katsevich algorithm performed as described in any one of claim 1-8 is System, it is characterised in that: this system include CT machine, image reconstruction work station, high-definition image display device, the X-ray bulb of CT machine, Detecting plate in CT machine;Its working method is: the relative invariant position of the X-ray bulb in CT machine and detecting plate, they along Scanning gantry 360 ° rotates and is irradiated reconstruction object, and the data for projection obtained on a frame-by-frame basis is transmitted by PCI-E bus In image reconstruction work station, utilize in image reconstruction work station as described in any one of claim 1-8 based on The image rebuilding method of Katsevich algorithm carries out accurately quickly rebuilding, and the reconstruction image obtained is shown by high-definition image and sets For showing.
CN201310734168.0A 2013-12-27 2013-12-27 A kind of image rebuilding method based on Katsevich algorithm and system Expired - Fee Related CN103714560B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310734168.0A CN103714560B (en) 2013-12-27 2013-12-27 A kind of image rebuilding method based on Katsevich algorithm and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310734168.0A CN103714560B (en) 2013-12-27 2013-12-27 A kind of image rebuilding method based on Katsevich algorithm and system

Publications (2)

Publication Number Publication Date
CN103714560A CN103714560A (en) 2014-04-09
CN103714560B true CN103714560B (en) 2016-10-05

Family

ID=50407501

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310734168.0A Expired - Fee Related CN103714560B (en) 2013-12-27 2013-12-27 A kind of image rebuilding method based on Katsevich algorithm and system

Country Status (1)

Country Link
CN (1) CN103714560B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105809723B (en) * 2016-03-09 2018-10-19 哈尔滨工业大学深圳研究生院 CBCT method for reconstructing and system
CN106056646B (en) * 2016-05-24 2018-02-23 广州华端科技有限公司 CT image rebuilding methods and system

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1788684A (en) * 2004-09-29 2006-06-21 株式会社东芝 Image reconstruction method and x-ray computed tomography apparatus

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7539281B2 (en) * 2004-10-15 2009-05-26 General Electric Company Methods and apparatus for reconstruction in helical cone beam volumetric CT

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1788684A (en) * 2004-09-29 2006-06-21 株式会社东芝 Image reconstruction method and x-ray computed tomography apparatus

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Studies on implementation of the Katsevichalgorithm for spiral cone-beam CT;Hengyong Yua et.al.;《Journal of X-Ray Science and Technology》;20041231;第12卷(第2期);97-116 *
圆弧轨迹Katsevich类型的锥束CT重建;王瑜等;《光电子激光》;20080731;第19卷(第7期);967-971 *
螺旋锥束CT三维图像重建算法中求导模块的设计与实现;贾喜元;《中国优秀硕士学位论文全文数据库_信息科技辑》;20111215;I138-1510 *
螺旋锥束工业CT的Katsevich精确重建算法研究;牛小明;《中国优秀硕士学位论文全文数据库_信息科技辑》;20090615;I138-768 *

Also Published As

Publication number Publication date
CN103714560A (en) 2014-04-09

Similar Documents

Publication Publication Date Title
Kurz et al. CBCT correction using a cycle-consistent generative adversarial network and unpaired training to enable photon and proton dose calculation
Yu et al. Deep sinogram completion with image prior for metal artifact reduction in CT images
Landry et al. Comparing Unet training with three different datasets to correct CBCT images for prostate radiotherapy dose calculations
CN103279929B (en) A kind of prediction of the CT image metallic traces based on integral cosine and artifact minimizing technology
CN102314698B (en) Total variation minimization dosage CT (computed tomography) reconstruction method based on Alpha divergence constraint
CN102945328B (en) Based on the x-ray imaging image simulation method of GPU concurrent operation
CN112102213B (en) Low-dose CT image processing method, scanning system and computer storage medium
CN104268914B (en) Reestablishing method of 4D-CT (Four Dimensional-Computed Tomography) different time phase sequence image
CN103295234B (en) Based on the medical image segmentation system and method for deformation surface model
CN103810734B (en) A kind of low dose X-ray CT data for projection restoration methods
CN103247061A (en) Augmented lagrangian iterative reconstruction method of X-ray image and CI image
CN106127825B (en) A kind of X ray CT image rebuilding method based on broad sense punishment weighted least-squares
WO2024051018A1 (en) Pet parameter image enhancement method and apparatus, device, and storage medium
CN105931280A (en) GPU-based fast three-dimensional CT (computed tomography) iterative reconstruction method
CN103793890A (en) Method for recovering and processing energy spectrum CT images
CN103714560B (en) A kind of image rebuilding method based on Katsevich algorithm and system
Dillon et al. Evaluating reconstruction algorithms for respiratory motion guided acquisition
Zhang et al. 3D remote healthcare for noisy CT images in the Internet of things using edge computing
CN103226815A (en) Low dose CT image filtering method
CN102024267A (en) Low-dose computed tomography (CT) image processing method based on wavelet space directional filtering
CN103340602B (en) Method for optimizing optimal viewing angle of multiple branch interesting blood vessel section
Li et al. MARGANVAC: metal artifact reduction method based on generative adversarial network with variable constraints
CN102243766A (en) New method for accelerating reconstruction of computerized tomography (CT) image
Yang et al. Cycle-consistent learning-based hybrid iterative reconstruction for whole-body PET imaging
Sun et al. ACnerf: enhancement of neural radiance field by alignment and correction of pose to reconstruct new views from a single x-ray

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20161005

Termination date: 20201227