WO2013068887A1 - Adaptive application of metal artifact correction algorithms - Google Patents
Adaptive application of metal artifact correction algorithms Download PDFInfo
- Publication number
- WO2013068887A1 WO2013068887A1 PCT/IB2012/056041 IB2012056041W WO2013068887A1 WO 2013068887 A1 WO2013068887 A1 WO 2013068887A1 IB 2012056041 W IB2012056041 W IB 2012056041W WO 2013068887 A1 WO2013068887 A1 WO 2013068887A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- image
- neighbourhood
- corrector
- artifact
- corrective action
- Prior art date
Links
- 238000012937 correction Methods 0.000 title claims abstract description 34
- 238000004422 calculation algorithm Methods 0.000 title description 19
- 239000002184 metal Substances 0.000 title description 8
- 230000003044 adaptive effect Effects 0.000 title description 4
- 230000009471 action Effects 0.000 claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000012545 processing Methods 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 16
- 230000001447 compensatory effect Effects 0.000 claims description 14
- 230000001419 dependent effect Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 5
- 230000003321 amplification Effects 0.000 claims description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 description 12
- 238000002591 computed tomography Methods 0.000 description 7
- 230000002238 attenuated effect Effects 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000003702 image correction Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 210000000988 bone and bone Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000007943 implant Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 210000004872 soft tissue Anatomy 0.000 description 2
- 230000003936 working memory Effects 0.000 description 2
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000135 prohibitive effect Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000002195 synergetic effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/77—Retouching; Inpainting; Scratch removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- 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/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- 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
- G06T2207/20012—Locally adaptive
-
- 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/20172—Image enhancement details
- G06T2207/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal filtering
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Definitions
- the invention relates to an apparatus for correcting an image for an image artifact, to a method of correcting an image for an image artifact, to a medical image processing system for correcting an image for an image artifact, to a computer program product, and to a computer readable medium.
- MAR Metal artifact reduction
- an apparatus for correcting an image for an artifact According to one aspect there is provided an apparatus for correcting an image for an artifact.
- the apparatus comprises an input unit for receiving an initial image, a processing unit for processing the image to effect the correction and an output unit for outputting a corrected version of the initial image.
- the processing unit comprises a correction sampler configured to forward the initial image to an artifact corrector and to receive therefrom a sample of a corrected image.
- the sample image is the result of a corrective action applied to the initial image by the artifact corrector.
- the processing unit further comprises a comparator configured to compare the initial image and the corrected sample image to so establish a corrector image representing the corrective action of the artifact corrector on the initial image.
- the processing unit further comprises an artifact correction controller configured to adaptively re-apply the respective corrective action in a weighted manner to any one of a plurality of image points in the initial image.
- the weights or weight factor used per image point for the weighting of the corrective action are related to a degree to which image information in a neighbourhood around that image point is compensable by
- the so adaptively re-corrected initial image is then output by the output unit as the corrected image.
- the corrective action is adaptively re-applied by the controller to obtain an improved, final corrected image, thereby ensuring previously present artifacts are removed and creation of new artifacts are avoided.
- the apparatus allows applying an existing artifact corrector implementing a known metal artifact reduction (MAR) algorithm.
- MAR metal artifact reduction
- the apparatus can be used as an "add-on" to existing MAR systems.
- the apparatus is a "meta-image-artifact corrector" because it examines correction results output by existing MAR systems and then improves the correction results by adaptively re-applying the corrective action using the corrector image as a "roadmap".
- the corrector image records the existing MAR system's estimate of corrective action for the initial image.
- the apparatus acts in a "doubly" adaptive manner , namely in a quantitative and a spatial sense: controller determines where, that is at which image points in the initial image, the estimated corrective action is to be applied and how much, that is attenuated or amplified, of the corrective action is to be applied at that locale.
- the apparatus operates in a two phase fashion: in a first, sample run, the artifact corrector is applied globally to the initial image to obtain the sample corrected image and to so gather information on the corrective action of the algorithm on the initial image.
- the corrective action is recorded in the corrector image.
- the corrector image uses the corrector image, the corrective action recorded therein is adaptively and selectively applied in a second, final run, but this time the corrective action is applied only locally in image regions where the "correction image" shows features which are also present in the original, initial image.
- the controller acts neighbourhood-wise. Centred around each image point in a previously defined region of the image plane (which may include the whole of the image plane safe for a border area around the image frame), a sub-region or neighbourhood
- patch is defined. This patch is then mapped into the corrector image to define a corrector image neighbourhood corresponding to the initial image neighbourhood to so obtain for each image point a pair of neighbourhoods.
- the compensatory degree measures the suitability of the image information in the two neighbourhoods to compensate or cancel each other.
- the degree or the corrective weights associated with the degree indicate the extent to which image information is mirrored or is similar across the pairs of neighbourhoods. As such the degree is a property of each of the respective pairs of initial image neighbourhoods and corrector image neighbourhoods. Because the neighbourhoods correspond to each other, either of the neighbourhood can be said to "have" the degree.
- An image artifact is an image feature defined by a local pixel pattern which is expected to show in the correction image also albeit in opposite pixel or voxel intensity so that the artifact can be compensated or annihilated after adaptive application of the corrector image to the initial image.
- MAR algorithms can be formulated in a way that there is an initial image, that suffers from metal artifacts and the algorithm creates a corrector image that is added to the initial image.
- the apparatus may therefore be put to use for MAR algorithms do not follow this form explicitly but they rather create directly a corrected image.
- the difference between the corrected image and the original image is calculated and this difference is used as the corrector image.
- the apparatus carries into effect the idea of applying the correction image "locally" in regions where it compensates artifacts but to avoid the application where it would create new artifacts.
- the apparatus allows distinguishing between these two cases by using the concept of the compensatory degree.
- the weights are computed using an entropy measure for the combined image information formed from the neighbourhood pair per centre initial image point.
- the statistical correlation coefficient between the image information in each of the pairs of neighbourhoods is computed.
- the controller is configured to adjust a previously set default size for the corrector image neighbourhood until the entropy of the image information in the corrector image neighbourhood exceeds a predetermined threshold value.
- Dynamically adjusting the neighbourhood size allows to keep the run-time of the algorithm at bay as this would increase with neighbourhood size.
- Choosing the threshold entropy value allows to balancing size for runtime: the entropy cost function may turn out rather flat for a too small a neighbourhood size due to lack of image information or structure in the small neighbourhood. A too large a neighbourhood however is computationally prohibitive and it further impedes a proper correction of artifact present in a small region.
- the apparatus selects an appropriate neighbourhood size by accounting for image structure in the neighbourhood, where "structure" is measured by the entropy function, where the entropy of an image is preferably defined as the entropy of the normalized histogram of grey values in that neighbourhood. If the correction image is rather flat in a neighbourhood, then the pixel grey value histogram is highly peaked and has high entropy. Using entropy as a structure measure, this means that a large neighbourhood should be selected. On the other hand, if the correction image has fine streaks in a neighbourhood, the histogram shows several peaks and has lower entropy. Consequently, a smaller
- the apparatus may also be put to use with artifact correction algorithms other than MARs.
- the artifact may be caused not necessarily by "metal" parts but by any other highly radiation attenuating part.
- the apparatus may be used for any image artifact corrector whatever the cause for the image artifact or whatever the particular algorithm underlying the existing artifact corrector.
- the controller controls the production of the final corrected image at the remote artifact controller.
- Amplifying corrective action includes maintaining the corrective action as provided by the MAR at a single given image point, the weight equalling at least unity in this later case.
- Attenuation of corrective action includes eliminating the corrective weight at a single given image point, the weight being around naught in this case.
- Image is to be construed broadly as an at least 2-dimensional array, matrix or similar data structure holding numerical data items, each addressable by at least two- dimensional coordinates
- Image information or feature is a particular pixel value pattern given by a particular pixel value distribution across the pixels making up the patch or region in the image plane.
- Figure 1 schematically shows a block diagram of an apparatus for image correction according to one embodiment of the invention
- Figure 2 diagrammatically shows operation of the apparatus of figure 1 according to one embodiment of the present invention
- Figure 3 shows a flow chart of a method of correcting images according to one embodiment of the present invention
- FIG. 1 To the right of Figure 1 there is shown a block diagram of an apparatus for correcting an image for an artifact.
- the apparatus comprises an input unit or interface means 105, a processing unit 110 and an output unit or interface means 150.
- Input 105 is configured to access a data base system 180 to retrieve therefrom an initial digital image having an image artifact.
- Data base 180 may be arranged as a PACS in a medical facility and is holding medical image data such as 2D computed tomography (CT) projection images or
- reconstructed 2D or 3D cross sectional images also commonly referred to as slice images or "slices”.
- the images are in a suitable digital format such as DICOM format.
- the apparatus is arranged to connect via interface means 105 and a suitable computer communication network to data base 180.
- FIG 1 To the left of Fig 1 there is shown an artifact corrector module MAR 190 implementing a known metal artifact reduction (MAR) algorithm.
- the embodiment of the apparatus as diagrammatically shown in Figure 1 is arranged as an "add-on" for an existing MAR 190 system.
- the MAR 190 may have access to additional data related to the image such as the projection data on which the initial image is based on, or a priori information about the highly attenuating part such as CAD information of orthopaedic implants.
- the apparatus may include a native MAR as a component of processing unit 110.
- the apparatus' processing unit 110 comprises a correction sampler 120, a comparator 130 and an artifact correction control 140.
- the apparatus components 105, 150, 120, 130 and 140 are running as software routines on processing unit 110.
- a distributed architecture of the apparatus where the components are connected in a suitable communication network is also contemplated in alternative embodiments.
- the components may also be arranged as dedicated FPGAs or hardwired standalone chips.
- the components may be programmed in a suitable scientific computing platform such as Matlab® or Simulink® and may then be translated into for, example, C++ or C routines maintained in a library and linked when called on by processing unit 110.
- the apparatus for correcting an image for an artifact connects via input unit 105 to data base 180 and retrieves therefrom an initial image having an image artifact.
- the image artifact may have been the result of a high density part such as a metal part residing in an object whilst CT images are taken of the object.
- the image may be a slice image reconstructed from a set of projection data acquired with a CT system.
- the object of interest may be a part of a human or animal body and the high density parts may include a metallic part or other high density part such as a bone structure embedded in patient's soft tissue.
- Image artifact such as streaks, distortions and shades may present in the acquired initial image due to the part exercising a higher attenuation on the radiation used in the CT run than the average attenuation of the surrounding soft tissue.
- correction sampler 120 After reading-in the initial image for image correction, correction sampler 120 forwards initial image to MAR 190 and requests MAR 190 to process the initial image to output a MAR corrected sample image which is then received back at the correction sampler 120.
- Comparator 130 compares the two images and produces a corrector image recording the corrective action of MAR 190 on the initial image.
- Corrector image represents the corrective action per image point ij of the
- Corrector image is then forwarded to artifact correction controller 140 to directly obtain the final corrected image by using the corrector image in combination with the initial image.
- Controller 140 then re-applies the corrective action to the initial image as per corrector image but this time uses weights to modify the corrective action per image point in the initial image.
- the individual weights or weight factors are computed by controller 140 as a parameter.
- a patch or neighbourhood is defined by the controller.
- controller 140 computes a numerical weight. Each weight measures the degree to which image information in the initial neighbourhood is compensable by the image information in the corrector image
- Controller 140 looks up in the corrector image the corrective action that was previously estimated by the MAR 190 to be applied to any given point.
- the estimated corrective action is then weighted by the respective weight and the so re-weighted corrective action is then applied by controller 140 to the respective image point in the initial image. Proceeding in this fashion for each or a selection of image points in the initial image, the final corrected image is build-up.
- the weights or parameters are computed by controller 140 according to different principles and each embodiment will be explained in more detail below under the heading "Operation”.
- Controller 140 then passes the final corrected image to output interface means 150.
- the final corrected image can then be dispatched across the communication network to data base 180 or can be forwarded to an image renderer and rendered for view on a screen.
- controller 140 instructs MAR unit to re- correct the initial image but this time the correction operation at MAR 190 is controlled by artifact correction controller 140 using the corrector image as explained above.
- controller 140 interfaces during runtime with MARs 190's algorithm and instructs the weights to be used in MAR 190's foreign algorithm.
- controller 140 is equipped with suitably programmed API's to effect the interaction and the image correction algorithm at MAR 190 may have to be suitably adapted to carry into effect the control function of controller 140 at MAR 190.
- MAR 190 may output the final corrected image and so release same for further post-processing to database 180 or otherwise. Operation
- Initial image is formed as a matrix including rows i and columns j.
- Matrix entry at row i and column j is a numerical value representing a grey value level of a pixel/voxel element.
- Each artifact in initial image is formed as distinctive image features defined by region of those pixel and voxel elements. The following reference will be made only to pixel elements but it is understood that the following applies equally to voxel elements.
- the sample correction image Gy produced at MAR 190 and requested by correction sampler 120 is formed of the same number of rows and columns as the image but has in general different pixel value entries because of the corrective action experienced at MAR 190.
- Comparator 130 is configured to generate the corrector image Cy based on initial image and sample image ⁇ ⁇ - .
- the comparator 130 forms the pixel-wise difference between the initial image and sample image ⁇ ⁇ - resulting in the corrector image Cy having the same number of rows and columns as initial image or sample image ⁇ ⁇ - . This difference value is then representative for the corrective action at that image point Corrector image Cy along with initial image is then forwarded for processing to controller 140.
- Controller 140 is configured to use corrector image ⁇ 3 ⁇ 4 ⁇ and to adaptively apply a repeated correction to initial image .
- Controller 140 is configured to apply the correction image c3 ⁇ 4 to initial image only locally, that is, application is restricted to regions where the corrective action of the corrector image Cy would compensate an artifact but application of corrected action is avoided or sufficiently attenuated in other regions of initial image where new artifacts would be created.
- controller 140 allows distinguishing between those two cases which will be explained in more detail below.
- Fig 2 affords a diagrammatical explanation of the operation of controller 140.
- initial image patch ⁇ may include an artifact diagrammatically shown as a dark ellipse.
- MAR 190 for dark ellipse artifact, correctly estimated the corrective action because corrector image ⁇ 3 ⁇ 4 ⁇ when applied to the image points in patch Q mn would compensate dark ellipse artifact completely.
- Image information in patch ⁇ _turn has therefore a high compensatory degree.
- Controller 140 therefore computes a parameter A larger unity to at least maintain or even amplify the corrective action as recorded by corrector image ⁇ 3 ⁇ 4 ⁇ for image points in that path patch ⁇ cultural tract .
- Patch ⁇ . ⁇ shows the opposite scenario.
- MAR 190 incorrectly estimated the corrective action because it would, if applied, introduce a dark rectangular artifact (shown to the lower right in Fig 3) where there is no artifact at all as shown in "clean" path ⁇ « ⁇ the lower left in Fig 3.
- Controller 140 therefore computes for patch ⁇ « a very low compensatory degree, because the structures in both patches do not match.
- Control parameter A is therefore attenuative with a value less than unity and close to naught to so annihilate or eliminate the MAR 190 proposed corrective action for image points in patch ⁇ « .
- controller 140 reads in initial image ⁇ defines for each point i, j neighbourhood called a "patch") around that point i,j.
- the point may be referred to as the centre point of that patch, each patch having at least one such centre or "seed" point.
- Neighbourhood ⁇ # defines a sub set in the image plain and may be given by a rectangular region with n x n pixel height and width. In an embodiment an 11 x 1 1 pixel square is chosen as the default size for neighbourhood ⁇ # .
- neighbourhood ⁇ # is a circle having a specified radius r around each of the image points Controller 140 then maps this neighbourhood ⁇ ,, ⁇ a corresponding neighbourhood ⁇ # in corrector image ⁇ 3 ⁇ 4. According to one embodiment this is done by using the same pixel co-ordinates as defined for the corresponding neighbourhood in corrector image ⁇ 3 ⁇ 4.
- the neighbourhoods are not defined for each and every pixel image pixel point in initial image but are restricted to a region which is likely to include the artifact.
- This region of interest in the image plane can be established for example by use of a suitably programmed segmentor.
- the image points i, j as centre points for the respective neighbourhoods ⁇ # are chosen sufficiently far away from the border of image to ensure that the neighbourhoods are well defined and would not extend beyond the current frame.
- controller 140 calculates the weights as the statistical correlation coefficient between pixel values in the patch ⁇ , ⁇ ⁇ initial image and the corresponding patch ⁇ in corrector image Cy . Calculation for statistical correlation coefficient is according to the following formula:
- ⁇ 1 ⁇ and c l ⁇ denote the average pixel value in initial image patch and corrector image patch, respectively.
- the pixel values of the final corrected image are then calculated according to:
- the range of values of are windowed and transformed into a selectable range by a clipping- function A.
- a partly sinusoidal function is defined: t ⁇ t 0
- Cut-off parameter to varies between minus 1 and zero and defines a cut-off point at which the correlation coefficients t tj are considered to be "negative enough" to warrant unchanged that is un-attenuated application of the corrective action by corrector image Cy at that point
- correlation between the patches is made more robust by tracking with a segmentor, regions that are representative of bones or other high density matter. Those regions are then excluded from the respective patches prior to the calculation of the correlation coefficients in respect of that patch.
- controller 140 is configured to calculate the weights according to the entropy with respect to each pair of patches.
- the weight for each patch pair is determined by the following formula: min ⁇ , H ⁇ ⁇ .. + A' c a ..)
- Entropy function H can be calculated by establishing a pixel level intensity histogram for each patch as indicated in formula (1) above. Intervals ("bins") of pixel values are defined for each patch and the number of pixels having a value falling in any one of the bins are recorded.
- a linear combination of the restricted images is formed in the vector space of matrices.
- the scalar measures how much of the corrective action of corrector image ⁇ 3 ⁇ 4 ⁇ should added in patch Q mn , thereby adding as little image information as possible, that is, to minimize the entropy of the linearly combined image.
- the bins for the histogram are chosen at increments of 10 Hounsfield Units (HU) and the default neighbourhood size is l lxl 1 pixels.
- the apparatus allows the user to configure bin levels and the default neighbourhood size.
- the corrective action at the respective neighbourhood centre point ij is then weighted by according to its degree as measured by the entropy of its neighbourhood pair as the image structure measure:
- the parameter calculated as the minimum entropy in the unit interval will assume a value 1 or close to 1 if the corrector image contains only artifact structures that are already present in the initial image .
- Weight Aij will assume a value of 0 or close to 0 if the corrector image ⁇ 3 ⁇ 4 ⁇ contains only new image information (a new artifact) that is not present in the initial image ⁇
- Improbable values are discarded and reset to a most probable value, for instance by determining according to:
- Aij min A , ⁇ ( ⁇ ⁇ .. + A' 3 ⁇ 4.,) - p(A'),
- p is an a priori probability function (also called “a prior") for the weighting factor, for example a Gaussian distribution function with expectation value of 1.
- corrector 140 can further be configured to use the calculated entropy when selecting the size, for example, edge width or diameter of the patches.
- the patch is selected based on the information content in that patch of the corrector image ⁇ 3 ⁇ 4 ⁇ .
- controller 140 is configured to choose a smaller than average patch in that region. According to one embodiment patches are therefore adaptively chosen according to the entropy of the region enclosed by the patch. According to one embodiment a default patch size is chosen and scaled according to the entropy calculated in that patch. Controller 140 uses a configurable average entropy value to affect a scaling with respect to that average value.
- controller 140 implements other image structure measures than entropy such as variance of image values or mean absolute differences from the mean.
- a flow chart is shown of a method of correcting an image for an image artifact.
- step S305 the initial image is received.
- step S310 the initial image is forwarded to an existing artifact corrector to receive therefrom a sample of a corrected image
- the sample corrected image is the result of an estimated corrective action applied to the initial image by the artifact corrector algorithm implemented by the artifact corrector.
- step S320 the initial image and the corrected sample image are compared.
- the result of the comparison is then recorded as a corrector image representing the corrective action of the artifact corrector per initial image point ij in the initial image.
- a default neighbourhood size is adjusted (S330) until the entropy of the image information in the previously set default sized corrector image neighbourhood exceeds a predetermined threshold value.
- the threshold and the default size are both configurable and the adjustment means in general to enlarge the neighbourhood to secure an adequate amount of entropy.
- a reasonable amount of entropy may for example be
- step S340 the respective corrective action is then adaptively re-applied in a weighted manner to any one of a plurality of image points in the initial image.
- the weights used per image point for the weighting of the corrective action relate to the degree to which image information in a neighbourhood around that image point is compensable by corresponding image information in a corresponding neighbourhood in the corrector image.
- the plurality may include all image points in the initial image (other than a band of image points along the border or the image frame) or a selection of image points that has been established to most likely include the image artifacts.
- a segmentation step may be used to so suitable prune down the image plane to a region of interest.
- the image plane may be trimmed to exclude image background areas.
- the weights are computed to minimize per neighbourhood the entropy of a combination of image information from both, the respective initial image neighbourhood and the corresponding corrector image neighbourhood.
- the weights are computed per neighbourhood as a statistical correlation coefficient between image information in the respective initial image neighbourhood and the corresponding image information in the corresponding corrector image neighbourhood.
- a computer program or a computer program element is provided that is characterized by being adapted to execute the method steps of the method according to one of the preceding embodiments, on an appropriate system.
- the computer program element might therefore be stored on a computer unit, which might also be part of an embodiment of the present invention.
- This computing unit may be adapted to perform or induce a performing of the steps of the method described above. Moreover, it may be adapted to operate the components of the above-described apparatus.
- the computing unit can be adapted to operate automatically and/or to execute the orders of a user.
- a computer program may be loaded into a working memory of a data processor. The data processor may thus be equipped to carry out the method of the invention.
- This exemplary embodiment of the invention covers both, a computer program that right from the beginning uses the invention and a computer program that by means of an up-date turns an existing program into a program that uses the invention.
- the computer program element might be able to provide all necessary steps to fulfil the procedure of an exemplary embodiment of the method as described above.
- a computer readable medium such as a CD-ROM
- the computer readable medium has a computer program element stored on it which computer program element is described by the preceding section.
- a computer program may be stored and/or distributed on a suitable medium, such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the internet or other wired or wireless telecommunication systems.
- a suitable medium such as an optical storage medium or a solid-state medium supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the internet or other wired or wireless telecommunication systems.
- the computer program may also be presented over a network like the World Wide Web and can be downloaded into the working memory of a data processor from such a network.
- a medium for making a computer program element available for downloading is provided, which computer program element is arranged to perform a method according to one of the previously described embodiments of the invention.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
Claims
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
IN3500CHN2014 IN2014CN03500A (en) | 2011-11-08 | 2012-10-31 | |
BR112014010843A BR112014010843A8 (en) | 2011-11-08 | 2012-10-31 | apparatus and method of correcting images for an artifact; medical image processing system for correcting an image for an image artifact; computer program element for controlling an apparatus; and computer readable media |
US14/354,685 US10089720B2 (en) | 2011-11-08 | 2012-10-31 | Adaptive application of metal artifact correction algorithms |
RU2014123282/08A RU2014123282A (en) | 2011-11-08 | 2012-10-31 | ADAPTIVE APPLICATION OF METAL ARTIFACT CORRECTION ALGORITHMS |
JP2014539451A JP6140716B2 (en) | 2011-11-08 | 2012-10-31 | Apparatus, method, system, and storage medium for correcting artifact |
CN201280054611.3A CN103918004B (en) | 2011-11-08 | 2012-10-31 | The self adaptation application of metal artifacts reduction algorithm |
EP12797985.4A EP2754122B1 (en) | 2011-11-08 | 2012-10-31 | Adaptive application of metal artifact correction algorithms |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201161556859P | 2011-11-08 | 2011-11-08 | |
US61/556,859 | 2011-11-08 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013068887A1 true WO2013068887A1 (en) | 2013-05-16 |
Family
ID=47297339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/IB2012/056041 WO2013068887A1 (en) | 2011-11-08 | 2012-10-31 | Adaptive application of metal artifact correction algorithms |
Country Status (8)
Country | Link |
---|---|
US (1) | US10089720B2 (en) |
EP (1) | EP2754122B1 (en) |
JP (1) | JP6140716B2 (en) |
CN (1) | CN103918004B (en) |
BR (1) | BR112014010843A8 (en) |
IN (1) | IN2014CN03500A (en) |
RU (1) | RU2014123282A (en) |
WO (1) | WO2013068887A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016064497A1 (en) * | 2014-10-21 | 2016-04-28 | General Electric Company | Methods and systems for normalizing contrast across multiple acquisitions |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3362987B1 (en) | 2015-10-14 | 2021-09-22 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image correction |
CN105225208B (en) * | 2015-10-14 | 2017-10-03 | 上海联影医疗科技有限公司 | A kind of computer tomography metal artifacts reduction method and device |
CN106651986B (en) * | 2016-01-21 | 2021-05-18 | 上海联影医疗科技股份有限公司 | Computer tomography artifact correction method |
US10013780B2 (en) * | 2016-02-29 | 2018-07-03 | General Electric Company | Systems and methods for artifact removal for computed tomography imaging |
CN106296615B (en) * | 2016-08-16 | 2017-09-29 | 广州华端科技有限公司 | The method and system of metal artifacts is corrected in CT images |
EP3398515B1 (en) * | 2017-05-03 | 2020-06-24 | Siemens Healthcare GmbH | Adaptive method for generating of ct image data with reduced artefacts, as well as image reconstruction unit and corresponding computer program product. |
US11151703B2 (en) | 2019-09-12 | 2021-10-19 | International Business Machines Corporation | Artifact removal in medical imaging |
CN113936068A (en) * | 2020-07-14 | 2022-01-14 | 上海联影医疗科技股份有限公司 | Artifact correction method, artifact correction device and storage medium |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006082563A1 (en) * | 2005-02-03 | 2006-08-10 | Koninklijke Philips Electronics N.V. | Radial adaptive filter for metal artifact correction |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5243664A (en) * | 1991-09-16 | 1993-09-07 | Picker International, Inc. | Post-processing technique for reducing metallic clip artifacts in CT images |
US5933471A (en) | 1997-02-26 | 1999-08-03 | International Business Machines Corporation | System and method for reducing reconstruction artifacts in computed tomography images |
US7001243B1 (en) | 2003-06-27 | 2006-02-21 | Lam Research Corporation | Neural network control of chemical mechanical planarization |
US7925074B2 (en) * | 2006-10-16 | 2011-04-12 | Teradyne, Inc. | Adaptive background propagation method and device therefor |
JP5405048B2 (en) * | 2007-07-10 | 2014-02-05 | 株式会社東芝 | X-ray computed tomography apparatus, reconstruction processing apparatus, and image processing apparatus |
EP2014232B1 (en) * | 2007-07-10 | 2021-06-30 | Toshiba Medical Systems Corporation | X-ray computed tomography apparatus, reconstruction processing apparatus and image processing apparatus |
ATE511161T1 (en) * | 2007-08-31 | 2011-06-15 | Koninkl Philips Electronics Nv | UNCERTAINTY MAP FOR SEGMENTATION IN THE PRESENCE OF METAL ARTIFACTS |
WO2010016425A1 (en) * | 2008-08-07 | 2010-02-11 | 株式会社 日立メディコ | X-ray ct image formation method and x-ray ct device using same |
US8244016B2 (en) | 2009-07-20 | 2012-08-14 | Wisconsin Alumni Research Foundation | Method for suppressing streak artifacts in images produced with an x-ray imaging system |
JP2011115404A (en) * | 2009-12-03 | 2011-06-16 | Canon Inc | X-ray image combining apparatus and x-ray image combining method |
CN101777177A (en) * | 2009-12-29 | 2010-07-14 | 上海维宏电子科技有限公司 | Attenuation filter-based metal artifact removing mixed reconstruction method for CT images |
US8768027B2 (en) | 2010-02-23 | 2014-07-01 | Carestream Health, Inc. | Method and system for cone beam computed tomography high density object artifact reduction |
-
2012
- 2012-10-31 JP JP2014539451A patent/JP6140716B2/en active Active
- 2012-10-31 US US14/354,685 patent/US10089720B2/en active Active
- 2012-10-31 WO PCT/IB2012/056041 patent/WO2013068887A1/en active Application Filing
- 2012-10-31 RU RU2014123282/08A patent/RU2014123282A/en not_active Application Discontinuation
- 2012-10-31 CN CN201280054611.3A patent/CN103918004B/en active Active
- 2012-10-31 BR BR112014010843A patent/BR112014010843A8/en not_active IP Right Cessation
- 2012-10-31 IN IN3500CHN2014 patent/IN2014CN03500A/en unknown
- 2012-10-31 EP EP12797985.4A patent/EP2754122B1/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006082563A1 (en) * | 2005-02-03 | 2006-08-10 | Koninklijke Philips Electronics N.V. | Radial adaptive filter for metal artifact correction |
US7991243B2 (en) | 2005-02-03 | 2011-08-02 | Koninklijke Philips Electronics N.V. | Radial adaptive filter for metal artifact correction |
Non-Patent Citations (1)
Title |
---|
WATZKE AND W A KALENDER O: "A pragmatic approach to metal artifact reduction in CT: merging of metal artifact reduced images", EUROPEAN RADIOLOGY, SPRINGER INTERNATIONAL, BERLIN, DE, vol. 14, no. 5, 1 January 2004 (2004-01-01), pages 849 - 856, XP007906226, ISSN: 0938-7994, DOI: 10.1007/S00330-004-2263-Y * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016064497A1 (en) * | 2014-10-21 | 2016-04-28 | General Electric Company | Methods and systems for normalizing contrast across multiple acquisitions |
US9619889B2 (en) | 2014-10-21 | 2017-04-11 | General Electric Company | Methods and systems for normalizing contrast across multiple acquisitions |
KR20170073640A (en) * | 2014-10-21 | 2017-06-28 | 제네럴 일렉트릭 컴퍼니 | Methods and systems for normalizing contrast across multiple acquisitions |
KR102493193B1 (en) | 2014-10-21 | 2023-01-27 | 제네럴 일렉트릭 컴퍼니 | Methods and systems for normalizing contrast across multiple acquisitions |
Also Published As
Publication number | Publication date |
---|---|
JP2014532506A (en) | 2014-12-08 |
BR112014010843A2 (en) | 2017-06-13 |
US20140286558A1 (en) | 2014-09-25 |
RU2014123282A (en) | 2015-12-20 |
BR112014010843A8 (en) | 2017-06-20 |
EP2754122A1 (en) | 2014-07-16 |
US10089720B2 (en) | 2018-10-02 |
JP6140716B2 (en) | 2017-05-31 |
CN103918004A (en) | 2014-07-09 |
EP2754122B1 (en) | 2019-01-23 |
CN103918004B (en) | 2016-12-14 |
IN2014CN03500A (en) | 2015-10-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10089720B2 (en) | Adaptive application of metal artifact correction algorithms | |
US9275454B2 (en) | Motion compensated second pass metal artifact correction for CT slice images | |
CN106683144B (en) | Image iterative reconstruction method and device | |
JP5462865B2 (en) | Use of non-attenuated corrected PET emission images to compensate for imperfect anatomical images | |
US10282823B2 (en) | Simulating dose increase by noise model based multi scale noise reduction | |
JP5080986B2 (en) | Diagnostic imaging system and diagnostic imaging method | |
US20060291711A1 (en) | Imaging chain for digital tomosynthesis on a flat panel detector | |
US20150228092A1 (en) | Digital breast tomosynthesis reconstruction using adaptive voxel grid | |
KR20120116940A (en) | Enhanced image data/dose reduction | |
US8433154B2 (en) | Enhanced contrast for scatter compensation in X-ray imaging | |
JP2007014759A (en) | Method for reducing image-based artifact in pet/ct imaging | |
JP2007520300A (en) | Reduction of overall image artifacts caused by highly attenuated objects in CT distributing voxel tissue classes | |
FR2862789A1 (en) | Tomosynthesis image data processing for medical diagnostics, involves identifying non-uniform weighting factors for back-projection processing of image data, and accessing and back-projecting data by application of factors | |
US20190156560A1 (en) | Method for enhanced display of image slices from 3-d volume image | |
KR102472464B1 (en) | Image Processing Method and Image Processing Device using the same | |
US7680352B2 (en) | Processing method, image processing system and computer program | |
KR102297972B1 (en) | Low Dose Cone Beam Computed Tomography Imaging System Using Total Variation Denoising Technique | |
CN111091514B (en) | Oral cavity CBCT image denoising method and system | |
JP2005021456A (en) | Device, method and program of image processing for radiographic image | |
JP2010509652A (en) | Image signal noise reduction | |
KR102473009B1 (en) | Low Dose Cone Beam Computed Tomography Imaging System Using Non-Local Total Variation Denoising Technique | |
Wiegert et al. | Improved CT image quality using a new fully physical imaging chain | |
CN117252943A (en) | Metal artifact removing method based on probability diffusion model (DDMP) | |
JP2023540128A (en) | Determination of target object type and location | |
CN115553798A (en) | Hardening correction method, system and storage medium |
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: 12797985 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2012797985 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14354685 Country of ref document: US |
|
ENP | Entry into the national phase |
Ref document number: 2014539451 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2014123282 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: 112014010843 Country of ref document: BR |
|
ENP | Entry into the national phase |
Ref document number: 112014010843 Country of ref document: BR Kind code of ref document: A2 Effective date: 20140506 |