X-ray examination apparatus with means for noise reduction

Abstract: Noise reduction means for use in an x-ray examination apparatus is provided for performing weighted temporally averaging in dependence on an amount of motion in parts of successive x-ray images. Further noise reduction is performed by combining temporally averaging with spatially filtering and motion detection. This reduces both noise breakthrough and the occurrence of trailers. In particular noise breakthrough is appropriately reduced by bi-temporal filtering. Threshold-values for discriminating between noise and motion are computed on the basis of images generated by the x-ray detector.

Description: The invention pertains to an x-ray examination apparatus for producing x-ray images and comprising means for noise reduction. The invention also relates to means for noise reduction suitable for use in an x-ray examination apparatus. An x-ray examination apparatus of said kind is described in the US Patent Specification US 4 367 490. In that document an x-ray examination apparatus is described which produces a time-sequence of x-ray images, i.e. it produces a sequence of successive x-ray images and converts said x-ray images into video signals having signal amplitudes in correspondence with brightness of picture elements (pixels) of the x-ray images. An x-ray examination apparatus known from the cited reference comprises means for noise reduction by recursively temporally averaging of pixel values. An important noise source in x-ray imaging is shot noise that is to be attributed to the quantised nature of the x-ray photons. By recursive temporally averaging images random noise can be suppressed, because uncorrelated information contained in an image is averaged out. However, when, in a medical x-ray examination motion is involved of e.g. the wall of a beating heart, or a catheter being moved by a physician, recursive temporally averaging leads to loss of image quality. In particular moving edges in the displayed images would have long trailers behind them. In an x-ray examination apparatus as described in the cited reference noise reduction is omitted for pixels which are involved in motion of a part of an object that is imaged. As a consequence noise-breakthrough occurs in parts of an image involving motion. I.e. noise is displayed in parts involved in motion in a displayed image. EP-A-0 289 152 relates to a noise reduction apparatus for processing an X-ray images, with a recursive filter combining succesive images and a means for mixing the outputs of the recursive filter and a 2D smoothing filter, the mixing ratio being dependent on motion. It is inter alia an object of the invention to provide an x-ray examination apparatus having means for applying noise reduction in stationary parts of an image produced by the x-ray examination apparatus as well as in parts of the image involved in motion. This is achieved by the features of claim 1. The dependent claims relate to embodiments. By the x-ray examination apparatus a sequence of successive images is produced. Noise reduction of said images is achieved by averaging. Motion detection is performed by comparing differences between pixel-values of corresponding pixels in successive images. Trailers and noise-tails in the displayed image are reduced by reducing the contribution to the averaging of preceding images regarding pixels involved in motion. Differences of pixel-values exceeding a predetermined threshold value are attributed to relevant pixels being involved in motion. Simultaneously, noise-breakthrough is reduced since some averaging is carried out even for pixels involved in motion. Careful balancing of averaging and omitting averaging achieves reduction of both trailers, noise-tails and noise-breakthrough. A preferred embodiment of an x-ray examination apparatus in accordance with the invention is characterised in that said means for weighted temporally averaging a sequence of successive images is arranged for weighted temporally averaging pixel-values of corresponding pixels of images from a sequence of successive images, and said motion detection means comprising means arranged for generating weight factors for said weighted temporally averaging in dependence on differences between pixel-values of corresponding pixels of images of said sequence. An advantage of performing noise reduction by processing pixel values is that comparatively simple circuit means can be employed for processing signal amplitudes in correspondence to pixel-values. A further preferred embodiment of an x-ray examination apparatus in accordance with the invention is characterised in that said noise reduction means is arranged for correcting a predetermined image in said sequence of successive image by averaging said predetermined image with images from said sequence preceding said predetermined image and with images from said sequence succeeding said predetermined image. In this preferred embodiment said noise reduction means comprises memory means for storing a time-sequence of images. An actual image being produced is an element substantially at the end of said time-sequence. An actual image being displayed is generated by weighted temporally averaging corresponding pixel-values of images of said time-sequence. Because in the weighted temporally averaging images that precede and images that succeed a predetermined element in a time-sequence are employed, such averaging may be denoted as bi-temporally averaging. A further advantage of bi-temporally averaging is that noise-breakthrough is avoided in places in a displayed image whereto a moving object has moved as well as a noise-tail is avoided in places in a displayed image wherefrom a moving object has moved away. The motion detector means comprises subtraction means for computing differences between pixel-values of images being successive elements of said time-sequence. Should a such difference exceed a predetermined threshold value, then a former one of pixel-values involved in computing said difference is left out of averaging the time-sequence. This provides the advantage that on the one hand trailers are avoided in the displayed image, but on the other hand a considerable noise reduction is sustained because of averaging over the remainder of the time-sequence. A further preferred embodiment of an x-ray examination apparatus in accordance with the invention is characterised in that said means for weighted temporally averaging is arranged for recursive temporally averaging. Noise reduction of an x-ray image by way of recursive averaging is performed by repeatedly for each pixel-value of an image to be corrected, adding to a pixel-value of the x-ray image at issue, pixel-values of corresponding pixels in preceding images, multiplied by a weight factor. I.e. noise reduction is performed by recursive averaging consisting of geometric summation of corresponding pixel-values pertaining to a time-sequence of image preceding an image at issue to an optional extent. A further preferred embodiment of an x-ray examination apparatus in accordance with the invention is characterised in that said noise reduction means comprises means for spatially filtering images and comprises computing means for computing a convex combination, in dependence of values of the weight factors generated by the motion detection means, of a weighted temporally averaged image and a corresponding spatially filtered image. When a substantial majority of pixel-values of images in the time-sequence have differences exceeding said threshold value then so few, possibly none, images in the time-sequence are involved in averaging that noise-breakthrough can occur. To prevent noise-breakthrough, spatial filtering is performed in dependence on the weight factors being determined by the motion detection means. Spatial filtering can be carried out in the form of spatial averaging, blurring weighted median filtering or median hybrid filtering. Pixel-values pertaining to a blurred image are computed by spatially filtering pixel-values of an image arising as a result of the weighted temporally averaging of the time-sequence. Subsequently there are computed convex combinations of pixel-values of the blurred image and pixel-values of the image arising as a result of the weighted averaging of the time-sequence. A convex combination of two images having pixel-values, I A and I B , say, is defined as an image having pixel-values I c = αI A + βI B , where α + β = 1 . As a magnitude of weight factors produced by the motion detection means decreases, then a contribution of pixel-values pertaining to the blurred image in said convex combination is increased. It is noted that independently combining temporal and spatial filtering per se for achieving noise reduction is known from the reference A.K. Katsaggelos et al., SPIE Visual Communications and Image Processing Vol. 1199 (1989)61-70. A further preferred embodiment of an x-ray examination apparatus in accordance with the invention is characterised in that said noise reduction means further comprises threshold-determination means for computing a threshold-value for controlling weighted temporally averaging on the basis of pixel-values in portions of images generated by the x-ray detector. In order to provide an accurate discrimination between noise and motion on the basis of differences of pixel-values of pixels pertaining to successive images in a time-sequence of images a threshold value must be determined. When differences of pixel-values exceed said threshold-value, then the difference at issue is attributed to motion. In order to rightly discriminate between noise and motion the threshold-value is preferably determined as a function of parameters describing the image at issue, e.g. brightness, contrast and artificial time-correlations caused by 'lag' of the imaging system. By way of said threshold-determination means an x-ray examination apparatus in accordance with the invention provides for discriminating between noise and motion in an adaptive way so as to avoid the motion detection means mistaking noise for motion or vice versa . A noise reduction means suitable for use in a x-ray examination apparatus preferably comprises means for weighted temporally averaging a set of images and motion detection means for generating weight factors for use in said weighted temporally averaging in dependence on differences between elements of said set of images. These and other aspects of the invention will become apparent from and elucidated with reference to the embodiments described hereinafter, with reference to the accompanying drawings. Figure 1 shows a block scheme of an embodiment of an x-ray examination apparatus in accordance with the invention. Figure 2 shows a block scheme of an embodiment of an x-ray examination apparatus in accordance with the invention notably comprising means for recursively averaging and means for spatially filtering. Figure 1 shows a block scheme of an embodiment of an x-ray examination apparatus in accordance with the invention, notably comprising means for averaging a time-sequence and means for spatially averaging. An x-ray source 1 is provided for generating a beam of x-rays 2 which is directed to an object 2, e.g. a patient. The x-rays are attenuated by the object and the attenuated x-radiation is collected by a x-ray image intensifier 3 having an input screen 4 and an output screen 5. An x-ray image collected on said input screen is converted into a visible image on the output screen. By means of an optical arrangement 6 said visible image is imaged onto a video camera 7, whereby the visible image is converted into an electrical video signal. The video signal is fed into noise reduction circuit means of which successive frame delay means are a part. In the figure frame delay means 10,11,12 and 13 are shown. Noise reduction for an image produced in the form of video signals is carried out as follows. First a temporally averaging is performed involving weight factors depending on presence of motion in the image. Further noise reduction is performed, notably in parts of the image involving motion and where temporal motion is consequently reduced, by forming a spatially filtered image. A spatially filtered image and a temporally averaged image are combined into a convex combination that is finally output from the noise reduction means. As a consequence video signals pertaining to a set of images in the form of a time sequence of successive images are available at positions indicated by I 1 , I 2 , I 3 , I 4 and I 5 , respectively. The embodiment shown in Figure 1 is discussed with the help of a time-sequence of five successive images; it is to be noted that the particular number of five successive images is not at all essential for the operation of an embodiment in accordance with the invention. Motion detection means are provided comprising subtraction means for computing differences of pixel-values pertaining to pairs of successive images. By means of subtraction means 14, 15, 16 and 17, differences of pixel-values pertaining to pairs of successive images are computed. Said differences are each fed into respective look-up tables, that are denoted by reference numerals 18, 19, 20 and 21, being a further part of the motion detection means. The respective look-up tables produce weight factors denoted as k 1 , k 2 , k 4 and k 5 , having values that decrease as the differences computed by subtraction means 14, 15, 16 and 17 respectively have increasing magnitudes. By way of a comparison means 50, the respective differences computed by subtraction means 14,15, 16 and 17 are compared to a threshold value provided by a threshold-determination means 51. Whenever a relevant difference as computed by anyone of the subtraction means 14, 15, 16 or 17 has a magnitude exceeding the threshold value, then a relevant weight factor from the set {k 1 , k 2 , k 4 , k 5 } is set to zero or to a very small value, less than ¼, say. A weighted average of pixel-values for a time-sequence of successive images is now computed from the time-sequence of images I 1-5 in the following manner. Pixel-values pertaining to image I n are multiplied by weight factors k n , for n = 1, 2, 4, 5, respectively, by way of respective multiplication means 23, 24, 25 and 26. Pixel-value of a predetermined image, here chosen to be the image I 3 are passed to an addition means 22, in other words these pixel-values are multiplied by unity. The results of the multiplication by said weight factors of said pixel-values are passed to the addition means 22. The result of addition by way of addition means 22 of pixel-values of image I 3 with weighted pixel-values of images I 1,2,4,5 is a weighted sum of pixel values. This weighted sum is further passed to a further multiplication means 27 to be converted in a normalised weighted average. To that end further addition means 28 are provided for computation of the sum k 1 + k 2 + k 4 + k 5 . To compute a required value of a normalisation constant there is provided a computation means 29 for computing a multiplicative normalisation factor N = 1/(1 + k 1 + k 2 + k 4 + k 5 ) . By way of multiplication means 27 the weighted sum of pixel-values produced by addition means 28 is multiplied by the normalisation factor produced by computation means 29, so as to provide normalised weighted averages of pixel-values, thereby forming a weighted temporally averaged video signal. To reduce noise-breakthrough, spatial filtering, e.g. denoted as 'blurring' hereinafter, is performed in dependence on the weight factors being determined by the motion detection means.The sum k 1 + k 2 + k 4 + k 5 as computed by the addition means 28 is further employed for performing blurring. To that end the value of the sum k 1 + k 2 + k 4 + k 5 is supplied to a look-up table 40 that provides a value for a weight factor b. The look-up table 40 is arranged such that the value of b increases as the value of the sum k 1 + k 2 + k 4 + k 5 decreases. A further computing means 41 is arranged for computing a value of 1-b. The weighted averaged video signal produced to by computation means 29 is on the one hand passed through a low-pass filter 42 so as to produce a spatially filtered component of the weighted averaged video signal. To produce a convex combination of the weighted temporally averaged video signal and the low-spatial-frequency component of the weighted temporally averaged video signal, the weighted temporally averaged video signal is multiplied by the factor 1-b by means of multiplication means 43 and the low-spatial-frequency component of the weighted temporally averaged video signal is multiplied by the factor b, by way of multiplication means 44. Finally, the required convex combination is computed by adding the results produced by multiplication means 43 and 44, by means of addition means 45. Thus the circuit means formed by the look-up table 40, computing means 41, multiplication means 43 and 44, low-pass filter 42 and addition means 45 constitute means for spatially filtering the weighted temporally averaged video signal in dependence of the result produced by the motion detection means, in that a convex combination is produced of the normalised weighted averaged video signal and a low-spatial frequency component of the actual video signal. The convex combination is formed by the addition means 45 and supplied to and output buffer circuit means 46. Figure 2 shows a block scheme of a embodiment of a x-ray examination apparatus in accordance with the invention notably comprising means for recursively averaging and means for spatially filtering. Circuit means 60 for recursively averaging are shown within the dashed box. A video signal S(t) pertaining to an actual image at time t is generated by the video camera 7 and a video signal S(t-τ) pertaining to a preceding image is produced by delay means 61. The delayed signal pertaining to the preceding image is subtracted from the actual signal by means of a subtraction means and the result of the subtraction is supplied to a look-up table 63 that produces a value for a recursion factor k , said value being decreasing when the magnitude of the result of the subtraction increases. Subsequently, a multiplication means 64 is provided for multiplying the result of the subtraction produced by subtraction means 62 by the recursion factor k. I.e. there is computed the value of k[S(t)-S(t-τ)] . Finally, recursive temporally averaging is completed by adding the video signal pertaining to the delayed image to the result of the multiplication performed by multiplication means 64. Viz. addition means 65 are provided for computing a value of S(t-τ)+k[S(t)-S(t-τ)]=kS(t)+(1-k)S(t-τ) , so that the addition means 65 produces at every stage a convex combination of an actual image and a delayed image. Therefore, by recursively temporally averaging of a actual image with its past is performed. For pixels involved in motion recursive averaging should preferably be omitted so as to avoid trailers in the image. To that end the subtraction means 62 operates also as a motion detection means. The difference of a actual video signal and a delayed video signal is attributed to motion whenever said difference has a magnitude exceeding a threshold-value. The threshold-value is computed by a threshold-determination means 70 using pixel values of the actual image. The difference of pixel-values of the actual video signal and the corresponding pixel-values of the delayed signal as computed by subtraction means 62 are compared to the threshold-value provided by the threshold-determination means 70, by employing a comparator 71. In case a difference at issue exceeds the threshold-value then the value of the recursion factor k is set to zero, so that recursive averaging is then omitted. For avoiding noise breakthrough whenever pixels involved in motion are not recursively averaged, further circuit means 80 are provided for performing spatial filtering. The actual video signal is passed through a low-pass filter 81 for producing a low-spatial-frequency component of the video signal of the actual image. The difference of the signals pertaining to the actual image and the recursively averaged image is supplied to a look-up table 82, which produces a value for a weight factor b . By way of computing means 83 a value for 1-b is computed. Subsequently,,the recursively averaged signal as produced by the addition means 65 is multiplied by means of a multiplication means 84 by the value of 1-b . Together with the low-spatial-frequency component of the actual video signal, the value of b is supplied to a further multiplication means 85, which multiplies the low-spatial-frequency component of the actual video signal by the weight factor b . Finally, an addition means 86 computes, employing the results produced by the multiplications means, 84 and 85, a convex combination of the recursively averaged signal and the low-spatial-frequency component of the actual signal; said convex combination is supplied to an output buffer circuit means 90. Threshold-determination means such as indicated in Figure 1 and in Figure 2 by reference numerals 51 and 70, respectively, comprises means for dividing a central region of interest of an image into portions of preferably 32 × 32 pixels, means for determining for each of said portions a maximum of differences between pixel-values of corresponding pixels of successive images and means for determining a minimum of said maxima of differences. A division of a central region of interest into portions of 32 × 32 pixels has empirically been found to yield reliable results for the threshold-value. if said portions are chosen to be larger, the probability of not having any portion not involving motion increases to an unacceptable level, if on the other hand, said portions are chosen to be smaller the statistical variation of the maxima for each of the portions increases so that the probability of finding a threshold-value that is too low increases to an unacceptable level. The minimal maximum pixel-difference per portion that is computed by the threshold determination means is employed as a threshold-value for discriminating between noise and motion. In stead of employing a maximum of differences between pixel values of successive images one can as an alternative employ other statistical parameters of the distribution of differences of pixel-values of successive images, e.g. one may use an average or a standard deviation. Further refinement of determining a threshold-value is achieved by employing a simple function of such a statistical parameter. For many sequences of medical x-ray images it has been found that noise reduction can be improved by employing a coefficient times said minimal-maximal pixel difference. An optimum value of said coefficient can be found empirically. Another refinement is related to distinguishing noise from motion in a situation when the patient is moved during the x-ray examination. In such a situation preferable the maximal maximum pixel-difference per portion as computed, and a fraction of said maximal maximum is employed as a threshold value. During motion of the patient the threshold value has a high value so that temporal filtering is improved. Upon terminating motion of the patient the threshold value decreases rapidly and temporal filtering for removing noise is restored.

Claims: An x-ray examination apparatus for generating a sequence of x-ray images and comprising motion detection means (14,15,16,17,18,19,20,21) for deriving weight factors from differences between pixel-values of successive x-ray images, said weight factors having a non-zero sum, characterised in that the x-ray examination apparatus is provided with noise reduction means (22,23,24,25,26) arranged to derive a bi-temporal average from pixel-values of a predetermined x-ray image, pixel-values of an x- ray image preceding said predetermined x-ray image, and pixel-values of an x-ray image succeeding said predetermined image, using said weight factors and wherein said weight factors have values that decrease as differences between pixel-values of successive x-ray images increase. An x-ray examination apparatus as claimed in Claim 1, characterised in that the sum of said weight factors equals unity. An x-ray examination apparatus as claimed in Claim 1, characterised in that the noise reduction means (22,27,43,44,45) is arranged to derive a low-frequency component of the bi-temporal average and a convex combination of the bi-temporal average and the low-frequency component of the bi-temporal average, which convex combination depends on the weight factors. An x-ray examination apparatus as claimed in Claim 1 or 2, characterised in that the noise reduction means is arranged to compute a threshold-value, to compare differences between pixel-values of successive x-ray images to the threshold-value (50) and to control derivation of the bi-temporal average on the basis of a difference between the threshold-value and said difference pixel-values of successive x-ray images.

Inventors: DE JONGE DIRK JAN ANDRIES, STEGEHUIS HERMAN, DE JONGE, DIRK JAN ANDRIES, STEGEHUIS, HERMAN
Applicants: KONINKL PHILIPS ELECTRONICS NV, KONINKLIJKE PHILIPS ELECTRONICS N.V.
Application year: 1993
Main classification: H04N5
IPC codes: H04N5