WO2013088294A1 - Automatic determination of regularization factor for iterative image reconstruction with regularization and/or image de-noising - Google Patents
Automatic determination of regularization factor for iterative image reconstruction with regularization and/or image de-noising Download PDFInfo
- Publication number
- WO2013088294A1 WO2013088294A1 PCT/IB2012/056929 IB2012056929W WO2013088294A1 WO 2013088294 A1 WO2013088294 A1 WO 2013088294A1 IB 2012056929 W IB2012056929 W IB 2012056929W WO 2013088294 A1 WO2013088294 A1 WO 2013088294A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image
- quality
- iteration
- regularization
- quality metric
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20004—Adaptive image processing
-
- 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
Definitions
- the following generally relates to iterative image reconstruction with regularization and/or image de-noising and finds particular application to computed tomography (CT), and is also amenable to other imaging modalities.
- CT computed tomography
- a CT scanner includes an x-ray tube mounted on a rotatable gantry that rotates around an examination region about a z-axis.
- a detector array subtends an angular arc opposite the examination region from the x-ray tube.
- the x-ray tube emits radiation that traverses the examination region.
- the detector array detects radiation that traverses the examination region and generates projection data indicative thereof.
- a reconstructor reconstructs, using an iterative or non- iterative reconstruction algorithm, the projection data, generating images indicative of the examination region.
- the noise level of the images depends on various imaging parameters and can be reduced through de-noising algorithms.
- Cost(f) F(f) + ⁇ ⁇ R ⁇ f) , where / represents an image, F(f) represents the data fidelity term, R(f) represents the regularization term, and ⁇ is a free parameter and represents the regularization factor that controls the strength of the regularization.
- the data fidelity term compares the forward projected image with the measured projections, optionally including a noise model on the projection.
- the data fidelity term can be the difference with the initial, noisy image.
- a noise model can be included.
- the regularization factor ⁇ is manually determined.
- ⁇ tends to be image and dataset dependent. Therefore, ⁇ has to be manually determined for each different image and/or dataset. This has been achieved by running several iterations with different ⁇ 's (in parallel or series) and selecting the ⁇ that satisfies predetermined criteria.
- this process can be computational and time intensive. For example, given a typical iteration takes two (2) to three (3) hours and ten (10) or more iterations may be ran before reaching a desirable final image, determining ⁇ for a particular image may take twenty (20) or more hours.
- a processing component processes images based on an iterative reconstruction algorithm with regularization and/or de-noising algorithm.
- the processing component includes a set point determiner that determines a quality set point between predetermined lower and upper quality bounds based on a quality variable indicative of an image quality of interest.
- the processing component further includes a comparator that compares, each processing iteration, a quality metric of a current generated image with the quality set point and generates a difference value indicative of a difference between the quality metric and the quality set point.
- the processing component further includes a regularization factor updater that generates an updated regularization factor for a next processing iteration based on a current value of the regularization factor and at least the quality metric in response to the difference value indicating that the quality metric is outside of a predetermined range about the quality set point.
- a method in another aspect, includes processing an initial image during a first processing iteration using an algorithm, which employs regularization, and an initial regularization factor, and generating a first image.
- the method further includes determining a first quality metric of the first image.
- the method further includes comparing the first quality metric of the first image with a quality set point and determining a difference value there between.
- the method further includes generating an updated regularization factor for a subsequent processing iteration of the image in response to the first quality metric being outside of a predetermined range about the quality set point.
- the updated regularization factor is generated based on the initial regularization factor and the first quality metric, and the updated regularization factor is used to process the first image during the subsequent processing iteration and generate a subsequent image.
- a processing component in another aspect, includes a component that processes an image using an iterative algorithm employing regularization and generates an image, wherein the iterative algorithm at least includes a regularization term multiplied by a regularization factor.
- the processing component further includes an analyzer that generates a quality metric after each processing iteration for a respective generated image.
- the processing component further includes a regularization factor determiner that automatically generates an updated regularization factor for a next processing iteration by the component based on a quality metric of a generated image of a current processing iteration and a predetermined quality set point for the generated image.
- the invention may take form in various components and arrangements of components, and in various steps and arrangements of steps.
- the drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention.
- FIGURE 1 schematically illustrates an example imaging system.
- FIGURE 2 schematically illustrates an example reconstructor that can be utilized with the imaging system of FIGURE 1.
- FIGURE 3 schematically illustrates an example image de-noiser that can be utilized with the imaging system of FIGURE 1.
- FIGURE 4 graphically illustrates convergence of an automatically determined regularization factor to a value corresponding to a predetermined quality metric.
- FIGURE 5 graphically illustrates a desired behavior for the convergence of the predetermined quality metric.
- FIGURE 6 illustrates an example method in accordance with at least one of FIGURES 2 or 3.
- the regularization factor ⁇ is automatically determined and the user only has to provide a relative level of regularization (e.g., between 0% (i.e., no regularization, or unregularized) and 100% (i.e., full regularization)), even though the regularization factor ⁇ may vary from image to image and/or data set to data set. In one instance, this mitigates manual determination of a suitable regularization factor ⁇ through trial and error each image and/or data set by a user and thus reduces the amount of time to reach a suitable
- FIGURE 1 illustrates an example imaging system 100, such as a computed tomography (CT) scanner.
- the imaging system 100 includes a generally stationary gantry 102 and a rotating gantry 104, which is rotatably supported by the stationary gantry 102 and rotates around an examination region 106 about a z-axis.
- CT computed tomography
- a radiation source 110 such as an x-ray tube, is rotatably supported by the rotating gantry 104.
- the radiation source 110 rotates with the rotating gantry 104 and emits radiation that traverses the examination region 106.
- a one or two dimensional radiation sensitive detector array 112 subtends an angular arc opposite the radiation source 110 across the examination region 106.
- the detector array 112 detects radiation traversing the examination region 106 and generates projection data indicative thereof.
- a subject support 118 such as a couch, supports an object or subject such as a human or animal subject in the examination region 106.
- a general-purpose computing system or computer serves as an operator console 120.
- the console 120 includes a human readable output device such as a monitor and an input device such as a keyboard, mouse, etc.
- Software resident on the console 120 allows the operator to interact with and/or operate the scanner 100 via a graphical user interface (GUI) or otherwise.
- GUI graphical user interface
- a reconstructor 114 reconstructs the projection data and generates image data indicative thereof.
- the reconstructor 114 employs an iterative reconstruction algorithm with regularization and determines the regularization factor ⁇ for each iteration based on an image quality metric of the image generated for the current iteration and a predetermined regularization level of interest, which can be a default and/or user defined level.
- the regularization factor is determined automatically.
- An optional image de-noiser 116 de-noises images generated by the reconstructor 114 and/or other reconstructor, including a reconstructor employing a conventional filtered-backprojection reconstruction algorithm and/or other regularized or non-regularized based reconstruction algorithm. As described in greater detail below, in one instance, the de-noiser 116, likewise, automatically determines the regularization factor ⁇ for a subsequent iteration based on an image quality metric of the current generated image generated for the current iteration and a predetermined regularization level of interest, which can be a default and/or use define level.
- the reconstructor 114 and the image de-noiser 116 may be part of the same component, for example, a processing component 122, as shown in FIGURE 1, or separate and distinct components. It is to be appreciated that automatically determining the regularization factor ⁇ for either or both of the iterative reconstruction algorithm with regularization or the image de-noising may enhance usability of iterative reconstruction algorithms with regularization and/or image de-noising, as a user does not have to know an absolute value of ⁇ , but rather just provides a relative level of regularization.
- the reconstructor 114 and/or the image de-no iser 116 are implemented via a processor (or microprocessor, controller, etc.) executing one or more computer readable instructions encoded or embedded on computer readable storage medium such as physical memory.
- a processor can be part of the console 120 and/or other computing device such as a dedicated computer, and/or other computing device.
- the processor can also execute at least one computer readable instructions carried by a carrier wave, a signal, or other non-computer readable storage medium such as a transitory medium.
- FIGURE 2 illustrates an example of the reconstructor 114.
- An image updater 202 receives an initial image and processes the image using an iterative reconstruction algorithm with regularization 204.
- An example of such an algorithm is described in Fessler, "Statistical image reconstruction methods for transmission tomography," Handbook of Medical Imaging, Volume 2: Medical Image Processing and Analysis, pages 1.70, SPIE, Bellingham, 2000.
- this reference discusses an approach for determining an estimate (I of the true attenuation map by maximizing a penalized- likelihood object function of the following form: , where
- ⁇ ( ⁇ ) ⁇ ( ⁇ ) - ⁇ ⁇ ⁇ ( ⁇ ) , ⁇ is the objective function, and R( i) is a roughness penalty.
- the regularization factor for the first reconstruction iteration, is set to an initial value of ⁇ 206, which can be a predetermined default value and/or user defined, and can be a zero or a non-zero value.
- the image updater 202 After such processing, the image updater 202 generates an image and likewise processes the generated image during a subsequent reconstruction iteration. In the illustrated example, the image updater 202 iteratively reconstructs a final image using a different updated value of ⁇ each reconstruction iteration.
- An analyzer 208 receives the generated image after each reconstruction iteration and analyzes the received image based on one or more given quality criteria such as noise, artifacts, etc.
- the analyzer 208 generates a quality metric indicative thereof.
- the quality metric can be a single value or multiple values. Examples of quality metrics include, but are not limited to, a global or a local image noise standard deviation or variance. Other quantitative criteria may optionally be used.
- the analyzer 208 outputs the quality metric after each reconstruction iteration.
- a ⁇ determiner 212 receives the quality metric of each generated image after each reconstruction iteration and determines the updated ⁇ for each next reconstruction iteration.
- the ⁇ determiner 212 includes a comparator 214 which compares the received quality metric of each reconstruction iteration with a predetermined quality set point 216 stored in memory 218, determining a value indicative of an error or difference between the received quality metric and the predetermined quality set point 216 each reconstruction iteration.
- the ⁇ determiner 212 conveys a signal to the image updater 202, which terminates the iterative reconstruction, and the image updater 202 outputs a final reconstructed image.
- Other stopping criteria may include the number of iterations, lapse a time duration, etc.
- a ⁇ updater 220 updates a current value of ⁇ 222, and conveys the updated ⁇ to the image updater 202 for the next reconstruction iteration.
- the set point 216 can be variously determined.
- a set point determiner 224 determines the set point 216 based on lower and upper quality bounds 226 and a quality variable 228, which indicates a level between the lower and upper quality bounds 226.
- the upper bound is determined as the quality metric of the initial unregularized image
- the lower bound is determined as the quality metric of a fully (or 100%) regularized initial image, which can be determined by fully regularizing the image or approximated without fully regularizing the initial image.
- the lower and upper quality bounds 226 can be normalized by the upper bound. Other approaches for determining the lower and upper quality bounds 226 are also contemplated herein.
- the quality variable 228, in this embodiment, represents a relative target level (e.g., expressed as a percentage or otherwise) between the lower and upper bounds 226 (e.g., 50%) between the lower and upper bounds 226).
- This relative target level can be determined and set by a user (e.g., a radiologist or other practitioner) or computer.
- FIGURE 3 illustrates an example of the image de-noiser 116.
- a de-noiser 302 receives an initial (e.g., a non- de-noised) image and de- noises the image based on a de-noising algorithm 304.
- the initial image can be generated by the imaging system 100 or other system.
- PCT international
- IB2011/054588 A non-limiting example of a suitable de-noising algorithm with regularization is described in international (PCT) application serial number PCT/IB2011/054588, entitled “Low Dose CT De-Noising,” filed on October 17, 2011, and assigned to Koninklijke Philips Electronics N V (Eindhoven, NL), the entirety of which is incorporated herein by reference.
- the de-noiser 302 uses an initial value of ⁇ 306 with the de-noising algorithm 304 to de-noise the received initial image.
- the initial value of ⁇ 306 can be determined as described herein for the initial value of ⁇ 206.
- the image de-noiser 302 After such processing, the image de-noiser 302 generates a de-noised image.
- the de-noiser 302 may also further de- noise a generated de-noised image, if an image quality of the de-noised image does not satisfy an image quality of interest, as described below. In the illustrated example, the de- noiser 302 iteratively further de-noises the generated de-noised image using a different updated value of ⁇ each iteration.
- the analyzer 208 and the ⁇ determiner 212 are as described in connection with FIGURE 2.
- the analyzer 208 receives the generated image after each de-noising iteration and analyzes the received image based on one or more given quality criteria such as noise, artifacts, etc.
- the analyzer 208 generates a quality metric indicative thereof.
- the quality metric can be a single value or multiple values. Examples of quality metrics include, but are not limited to, a global or a local image noise standard deviation or variance. Other quantitative criteria may optionally be used.
- the analyzer 208 outputs the quality metric after each image de-noising iteration.
- a ⁇ determiner 212 receives the quality metric of each generated image after each de-noising iteration and determines the updated ⁇ for each next de-noising iteration.
- the ⁇ determiner 212 includes a comparator 214 which compares the received quality metric of each de-noising iteration with a predetermined quality set point 216 stored in memory 218, determining a value indicative of an error or difference between the received quality metric and the predetermined quality set point 216 each de-noising iteration.
- the ⁇ determiner 212 conveys a signal to the image de-noiser 302, which terminates de-noising, and the image de-noiser 302 outputs a final de-noised image.
- other stopping criteria may include the number of iterations, lapse of a time duration, etc.
- a ⁇ updater 220 updates a current value of ⁇ 222, and conveys the updated ⁇ to the image de- noiser 302 for the next image de-noising iteration.
- the upper bound is determined as the quality metric of the initial (non- de-noised) image
- the lower bound is determined as the quality metric of a fully (100%) de-noised initial image, which can be determined by fully de-noising the image or approximated without fully de-noise the initial image.
- the lower and upper quality bounds 226 can be normalized by the upper bound. Other approaches for determining the lower and upper quality bounds 226 are also contemplated herein.
- the quality variable 228 can be set based on dose, which allows for dose reduction of the image for a given image quality. For example, where the practitioner is going to set the quality variable 228 at 50%, patient dose can be reduced four times for a given image quality, as the relationship between dose and noise can be determined. In another instance, the final image quality can be improved for a given patient dose. In yet another example, both image quality can be improved and dose can be reduced. Other approaches for determining the set point 216 are also contemplated herein.
- ⁇ (k) ⁇ (k - 1) ⁇ 2 p ⁇ e ⁇ k)+IE ⁇ k)+D e ⁇ k)) , where ⁇ (k) is the updated ⁇ for the next iteration, ⁇ (k-1) is the previous updated ⁇ , P, I and D respectively are proportionality, integrating and derivative constants for PID control, e(k) represents the difference between the quality metric of the current iteration and the set point of the current iteration, E(k) is the summation of e(k) up to the k-th iteration, and Ae(k) represents the change (i.e., e(k) - e(k - V) ) of the difference between the quality metric and the set point between the (k-l)-th to the k-th iteration.
- ⁇ can be updated using P, PI, and/or other control.
- FIGURE 4 graphically illustrates a non-limiting example of updating ⁇ over multiple iterations using a maximum likelihood reconstruction with a Huber penalty where the quality metric is the global noise level of the image.
- a first y-axis 402 represents the quality metric range between the lower and the upper bounds 226, normalized to the upper bound
- a second y-axis 404 represents the value of ⁇
- an x-axis 406 represents the number of reconstruction or de-noising iterations.
- a quality set point 408 is set to 50% between normalized upper and lower bounds (or 0.5).
- a quality metric curve 410 shows that the quality metric converges to the 50% point after approximately 500 iterations.
- a ⁇ curve 412 shows the value of ⁇ converges to and remains fairly constant at approximately 250 once the quality metric reaches the set point. This example is provided for explanatory purposes and is not limiting.
- FIGURE 5 graphically illustrates an example target behavior 500 for the quality metric.
- a y-axis 502 represents the quality range between the lower and the upper bounds 226, normalized to the upper bound
- an x-axis 504 represents the number of reconstruction or de-noising iterations.
- the target behavior 500 includes an approximately linear region 506 (i.e., a linear decay with a fixed slope) until a set point 508 is reached to mitigate possible overshoot of the set point 216.
- Other target behavior is also contemplated herein. This example behavior is provided for explanatory purposes and is not limiting.
- FIGURE 6 illustrates an example method in accordance with the embodiments described herein.
- an image is obtained for processing with an algorithm employing a regularization factor.
- the regularization factor is set to an initial value.
- the image is processed, for example, as described herein via an iterative reconstruction algorithm with regularization or a de-noising algorithm, producing a processed image.
- a quality metric is determined for the processed image.
- the quality metric is compared with a quality set point and a difference value indicative of a difference between the quality metric and the quality set point is determined.
- the quality set point can be determined based on target level between a lower quality bound and an upper quality bound 226.
- the regularization factor is updated as described herein and acts 604-610 are repeated with the updated regularization factor. If the quality metric is within the predetermined range, then at 614, a final image is output.
Abstract
Description
Claims
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP12812393.2A EP2748798B1 (en) | 2011-12-13 | 2012-12-04 | Automatic determination of regularization factor for iterative image reconstruction with regularization and/or image de-noising |
US14/364,715 US9076191B2 (en) | 2011-12-13 | 2012-12-04 | Automatic determination of regularization factor for iterative image reconstruction with regularization and/or image de-noising |
IN4759CHN2014 IN2014CN04759A (en) | 2011-12-13 | 2012-12-04 | |
CN201280061612.0A CN104025156B (en) | 2011-12-13 | 2012-12-04 | Processing component and method used for processing images on the basis of regularization iteration reconstruction algorithm and/or de-noising algorithm |
RU2014128596A RU2014128596A (en) | 2011-12-13 | 2012-12-04 | AUTOMATIC DETERMINATION OF THE REGULARIZATION COEFFICIENT FOR ITERATIVE RECOVERY OF THE IMAGE WITH REGULARIZATION AND / OR NOISE SUPPRESSION OF THE IMAGE |
JP2014546681A JP6141313B2 (en) | 2011-12-13 | 2012-12-04 | Automatic determination of regularization coefficients for iterative image reconstruction using regularization and / or image denoising |
BR112014014093A BR112014014093A2 (en) | 2011-12-13 | 2012-12-04 | processing component, which processes images based on iterative reconstruction algorithm with smoothing and / or noise reduction algorithm and method |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161569835P | 2011-12-13 | 2011-12-13 | |
US61/569,835 | 2011-12-13 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013088294A1 true WO2013088294A1 (en) | 2013-06-20 |
Family
ID=47520208
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/IB2012/056929 WO2013088294A1 (en) | 2011-12-13 | 2012-12-04 | Automatic determination of regularization factor for iterative image reconstruction with regularization and/or image de-noising |
Country Status (8)
Country | Link |
---|---|
US (1) | US9076191B2 (en) |
EP (1) | EP2748798B1 (en) |
JP (1) | JP6141313B2 (en) |
CN (1) | CN104025156B (en) |
BR (1) | BR112014014093A2 (en) |
IN (1) | IN2014CN04759A (en) |
RU (1) | RU2014128596A (en) |
WO (1) | WO2013088294A1 (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015044237A1 (en) | 2013-09-30 | 2015-04-02 | Koninklijke Philips N.V. | Method for local adjustment of regularization parameters for image quality optimization in fully 3d iterative ct reconstruction |
WO2016042466A3 (en) * | 2014-09-15 | 2016-05-12 | Koninklijke Philips N.V. | Iterative image reconstruction with a sharpness driven regularization parameter |
WO2016091594A1 (en) * | 2014-12-10 | 2016-06-16 | Koninklijke Philips N.V. | Statistically weighted regularization in multi-contrast imaging |
WO2019149359A1 (en) | 2018-01-31 | 2019-08-08 | Mitos Gmbh | Method for image reconstruction of an object, in particular based on computed-tomography image reconstruction, and apparatus, system and computer program product for the same |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2923332B1 (en) * | 2012-11-26 | 2019-01-16 | Koninklijke Philips N.V. | Projection data de-noising |
US9406107B2 (en) * | 2013-12-18 | 2016-08-02 | General Electric Company | System and method of computed tomography signal restoration via noise reduction |
US9615803B2 (en) * | 2014-06-16 | 2017-04-11 | General Electric Company | System and method for determining X-ray exposure parameters |
KR102372165B1 (en) * | 2015-01-22 | 2022-03-11 | 삼성전자주식회사 | X-ray imaging apparatus, image processing apparatus and method for image processing |
GB2542764A (en) * | 2015-09-23 | 2017-04-05 | Pathxl Ltd | Image processing method and apparatus for normalisation and artefact correction |
US9911208B2 (en) * | 2016-04-11 | 2018-03-06 | Toshiba Medical Systems Corporation | Apparatus and method of iterative image reconstruction using regularization-parameter control |
CN109844815A (en) * | 2016-08-22 | 2019-06-04 | 皇家飞利浦有限公司 | The image procossing based on feature is carried out using from the characteristic image of different iterative extractions |
US20200175732A1 (en) * | 2017-06-02 | 2020-06-04 | Koninklijke Philips N.V. | Systems and methods to provide confidence values as a measure of quantitative assurance for iteratively reconstructed images in emission tomography |
US10853689B2 (en) * | 2017-10-09 | 2020-12-01 | WebFurther, LLC | Methods for more effectively moderating one or more images and devices thereof |
CN107997780B (en) * | 2018-01-19 | 2020-11-06 | 重庆大学 | Cone beam CT instantaneous scanning device and reconstruction method |
CN112529796B (en) * | 2020-12-02 | 2022-07-05 | 华北理工大学 | CNN medical CT image denoising method based on multi-feature extraction |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61225640A (en) * | 1985-03-30 | 1986-10-07 | Toshiba Corp | Tomographic camera |
US20040034304A1 (en) * | 2001-12-21 | 2004-02-19 | Chikayoshi Sumi | Displacement measurement method and apparatus, strain measurement method and apparatus elasticity and visco-elasticity constants measurement apparatus, and the elasticity and visco-elasticity constants measurement apparatus-based treatment apparatus |
JP3887774B2 (en) * | 2001-12-21 | 2007-02-28 | 親良 炭 | Displacement vector measuring device and strain tensor measuring device |
US7864999B2 (en) | 2005-10-19 | 2011-01-04 | Siemens Medical Solutions Usa, Inc. | Devices systems and methods for processing images |
DE102005052061A1 (en) * | 2005-11-01 | 2007-05-16 | Carl Zeiss Imaging Solutions G | Method and device for image processing |
US7558414B2 (en) * | 2006-09-11 | 2009-07-07 | Case Western Reserve University | Iterative image reconstruction |
US8468189B2 (en) * | 2006-11-08 | 2013-06-18 | Okinawa Institute Of Science And Technology Promotion Corporation | Iterated variational regularization combined with componentwise regularization |
US8233682B2 (en) * | 2007-06-05 | 2012-07-31 | General Electric Company | Methods and systems for improving spatial and temporal resolution of computed images of moving objects |
EP2232446B1 (en) | 2007-12-20 | 2013-04-17 | Wisconsin Alumni Research Foundation | Method for prior image constrained image reconstruction |
US7885371B2 (en) | 2008-08-28 | 2011-02-08 | General Electric Company | Method and system for image reconstruction |
CN102159965B (en) * | 2008-09-17 | 2014-09-24 | 皇家飞利浦电子股份有限公司 | B1-mapping and b1l-shimming for mri |
DE102009014723B4 (en) | 2009-03-25 | 2012-10-25 | Siemens Aktiengesellschaft | Contrast-dependent regularization strength in the iterative reconstruction of CT images |
EP2424438A4 (en) * | 2009-04-30 | 2013-04-03 | Univ California | System and methods for fast implementation of equally-sloped tomography |
DE102010029281A1 (en) | 2010-05-25 | 2011-12-01 | Siemens Aktiengesellschaft | Method and image reconstruction device for the reconstruction of image data |
EP2633494B1 (en) | 2010-10-27 | 2019-12-11 | Koninklijke Philips N.V. | Low dose ct denoising |
US20140003688A1 (en) * | 2011-03-17 | 2014-01-02 | Koninklijke Philips N.V. | Multiple modality cardiac imaging |
CN102231204A (en) * | 2011-06-20 | 2011-11-02 | 南京航空航天大学 | Sequence image self-adaptive regular super resolution reconstruction method |
US20130108010A1 (en) * | 2011-10-31 | 2013-05-02 | Eigenor Oy | Method and apparatus for performing tomographic reconstruction |
US8879811B2 (en) * | 2012-03-28 | 2014-11-04 | Siemens Aktiengesellschaft | Alternating direction of multipliers method for parallel MRI reconstruction |
-
2012
- 2012-12-04 CN CN201280061612.0A patent/CN104025156B/en active Active
- 2012-12-04 JP JP2014546681A patent/JP6141313B2/en active Active
- 2012-12-04 US US14/364,715 patent/US9076191B2/en active Active
- 2012-12-04 RU RU2014128596A patent/RU2014128596A/en not_active Application Discontinuation
- 2012-12-04 BR BR112014014093A patent/BR112014014093A2/en not_active IP Right Cessation
- 2012-12-04 WO PCT/IB2012/056929 patent/WO2013088294A1/en active Application Filing
- 2012-12-04 EP EP12812393.2A patent/EP2748798B1/en active Active
- 2012-12-04 IN IN4759CHN2014 patent/IN2014CN04759A/en unknown
Non-Patent Citations (4)
Title |
---|
ANNA K JEREBKO ET AL: "Regularization Parameter Selection in Maximum a Posteriori Iterative Reconstruction for Digital Breast Tomosynthesis", 16 June 2010, DIGITAL MAMMOGRAPHY : 10TH INTERNATIONAL WORKSHOP, IWDM 2010, GIRONA, SPAIN, JUNE 16-18, 2010 ; PROCEEDINGS, SPRINGER-VERLAG, BERLIN, PAGE(S) 548 - 555, ISBN: 978-3-642-13665-8, XP019144399 * |
FESSLER: "Handbook of Medical Imaging, Volume 2: Medical Image Processing and Analysis", vol. 2, 2000, article "Statistical image reconstruction methods for transmission tomography", pages: 1.70 |
HAIYONG LIAO ET AL: "Selection of regularization parameter in total variation image restoration", JOURNAL OF THE OPTICAL SOCIETY OF AMERICA A, vol. 26, no. 11, 1 November 2009 (2009-11-01), pages 2311 - 2320, XP055058094, ISSN: 1084-7529, DOI: 10.1364/JOSAA.26.002311 * |
HE LINGJUN ET AL: "Adaptive regularized MAP of CT image reconstruction method", COMPUTER ENGINEERING AND APPLICATIONS, vol. 47, no. 28, 1 October 2011 (2011-10-01), pages 198 - 200, XP055056173, ISSN: 1002-8331, DOI: 10.3778/j.issn.1002-8331.2011.28.055 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015044237A1 (en) | 2013-09-30 | 2015-04-02 | Koninklijke Philips N.V. | Method for local adjustment of regularization parameters for image quality optimization in fully 3d iterative ct reconstruction |
US9478049B2 (en) | 2013-09-30 | 2016-10-25 | Koninklijke Philips N.V. | Method for local adjustment of regularization parameters for image quality optimization in fully 3D iterative CT reconstruction |
WO2016042466A3 (en) * | 2014-09-15 | 2016-05-12 | Koninklijke Philips N.V. | Iterative image reconstruction with a sharpness driven regularization parameter |
JP2017521124A (en) * | 2014-09-15 | 2017-08-03 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Iterative image reconstruction using sharpness-driven regularization parameters |
US9959640B2 (en) | 2014-09-15 | 2018-05-01 | Koninklijke Philips N.V. | Iterative image reconstruction with a sharpness driven regularization parameter |
WO2016091594A1 (en) * | 2014-12-10 | 2016-06-16 | Koninklijke Philips N.V. | Statistically weighted regularization in multi-contrast imaging |
CN107004259A (en) * | 2014-12-10 | 2017-08-01 | 皇家飞利浦有限公司 | Statistical weight regularization in many contrast imagings |
US10297049B2 (en) | 2014-12-10 | 2019-05-21 | Koninklijke Philips N.V. | Statistically weighted regularization in multi-contrast imaging |
WO2019149359A1 (en) | 2018-01-31 | 2019-08-08 | Mitos Gmbh | Method for image reconstruction of an object, in particular based on computed-tomography image reconstruction, and apparatus, system and computer program product for the same |
Also Published As
Publication number | Publication date |
---|---|
EP2748798B1 (en) | 2018-04-11 |
JP6141313B2 (en) | 2017-06-07 |
CN104025156B (en) | 2017-05-24 |
EP2748798A1 (en) | 2014-07-02 |
BR112014014093A2 (en) | 2017-06-13 |
US9076191B2 (en) | 2015-07-07 |
CN104025156A (en) | 2014-09-03 |
JP2015500118A (en) | 2015-01-05 |
IN2014CN04759A (en) | 2015-09-18 |
US20140348440A1 (en) | 2014-11-27 |
RU2014128596A (en) | 2016-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9076191B2 (en) | Automatic determination of regularization factor for iterative image reconstruction with regularization and/or image de-noising | |
JP5694357B2 (en) | Improved image data / dose reduction | |
CN106683144B (en) | Image iterative reconstruction method and device | |
US9159122B2 (en) | Image domain de-noising | |
JP6275826B2 (en) | Noise removal reconstruction image data edge improvement | |
US20160225169A1 (en) | Method for local adjustment of regularization parameters for image quality optimization in fully 3d iterative ct reconstruction | |
JP6158910B2 (en) | Iterative image reconstruction by regularization | |
US20160007948A1 (en) | Method and system for reducing localized artifacts in imaging data | |
EP3195265A2 (en) | Iterative image reconstruction with a sharpness driven regularization parameter | |
WO2011161557A1 (en) | Method and system for noise reduction in low dose computed tomography | |
JP6062250B2 (en) | X-ray computed tomography (X-ray CT) using the successive approximation method | |
WO2012073167A1 (en) | Iterative reconstruction algorithm with a constant variance based weighting factor | |
US20140153809A1 (en) | Image generation apparatus |
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: 12812393 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2014546681 Country of ref document: JP Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14364715 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2014128596 Country of ref document: RU Kind code of ref document: A |
|
REG | Reference to national code |
Ref country code: BR Ref legal event code: B01A Ref document number: 112014014093 Country of ref document: BR |
|
ENP | Entry into the national phase |
Ref document number: 112014014093 Country of ref document: BR Kind code of ref document: A2 Effective date: 20140610 |