CN1595957A - Method for determining automatic detection threshold of bad pixel of medical image - Google Patents

Abstract

The invention is a confirming method for medical image bad dots automatic detection threshold. The invention detects the bad dot with formula alpha(x, y) = (grav(x, y)-avg(x, y))/ avg(x, y) and I median= Median (I original), I (x, y) =| I median - I original|, the automatic checking threshold of bad dot includes: collecting detection image, threshold detection, video detection, and detection threshold definition. It is a process method combining threshold detection and video detection. The method is simple, upgrades the detection speed, and guarantees the high accuracy.

Description

The automatic detection threshold of the bad point of medical image is determined method

Technical field

The inventive method belongs to the medical image processing technical field, definite method of the automatic detection threshold of the bad point of particularly a kind of medical image.

Background technology

Bad point in the digital X ray machine in the image system generally is meant the unusual pixel of gray value performance, is usually expressed as " bright spot " or " dim spot " with respect to background." bright spot " is meant gradation of image obviously greater than the pixel of facing the territory average gray, and " dim spot " is meant that gradation of image is significantly less than the pixel of facing the territory average gray.At present vision-based detection is that bad point detects an a kind of method more accurately, but for the image of 1024*1024, shows at data line that at every turn need carry out 1024 times vision-based detection altogether, too much manual detection number of times is unfavorable for increasing work efficiency.And the method for the bad point of Threshold detection, though can carry out the problem that exists detection threshold to set automatically.Detection threshold is the key point that bad point detects.United States Patent (USP) " Automaticidentification and correction of bad pixels in a large area solid state x-ray detector " (U.S.Pat.No.5,657,400), " Defective pixel detection circuit of a solid state image pick-up device capable ofdetecting defective pixels with low power consumption and high precision, and image pick-updevice having such detection circuit " (U.S.Pat.No.5,854,655), detect though proposed to utilize detection threshold to go bad a little, do not relate to and how to set detection threshold.

Summary of the invention

Problem at the prior art existence, the inventive method provides definite method of the automatic detection threshold of the bad point of a kind of medical image, be a kind of simple and practical automatic detection of bad point and the image processing method of bad point calibration, from mathematical computations and two angles of human eye vision, Threshold detection and the visible detection method of differentiating dead pixel points of images are combined, on automatic detection of thresholding and vision-based detection basis, determine suitable detection threshold, reduced the bad point of medical image probability of false detection.

The definition of the bad point of the medical image that the present invention relates to, a kind of is the definition of adopting Agilent (Agilent Technologies):

$\mathrm{\α}(x,y)=\frac{\mathrm{gray}(x,y)-\mathrm{avg}(x,y)}{\mathrm{avg}(x,y)}--\left(1\right)$

(x is that ((x is that (x is the mean value of the image-region S pixel grey scale at center y), and its computing formula is with pixel y) to avg to pixel for x, gray value y) y) to gray in the formula $\mathrm{avg}(x,y)=\frac{1}{N}\underset{(i,j)\∈S}{\mathrm{\Σ}}\mathrm{gray}(i,j),$ Wherein S is that (x y) is the image-region at center, and (i is that (i, j) gradation of image, N are that the whole zone removal of image pixel grey scale is the number of pixels of two pixels of extreme value to pixel j) to gray with pixel.If | α (x, y) |＞detection threshold, (x y) is bad point then to adjudicate this point.

Another kind is at U.S.Pat.No.5, in 854,655, adopts medium filtering to handle, and bad point of elimination with the image subtraction of original image and Filtering Processing, is judged as bad point for subtracting each other the pixel of result greater than thresholding then.Formula is described as (2) (3).

I _median＝Median(I _originall) (2)

I(x，y)＝|I _median-I _originall| (3)

I in the formula _OriginallBe original image, I _MedianBe median-filtered result, Median () expression medium filtering is handled.If I (x, y)＞detection threshold, (x y) is bad point then to adjudicate this point.

According to above definition, the inventive method determines that the automatic detection threshold of bad point comprises the bad point of collection detected image, thresholding detection, vision-based detection, determines four steps of detection threshold, and its bad some testing process as shown in Figure 1.

Step 1: obtain bad some detected image

The images acquired sequence obtains bad some detected image.

Step 2: Threshold detection

It is as follows to utilize formula (1) to carry out the step of Threshold detection:

1) according to the definition of image size and bad some image is divided into several zonules according to the size of certain pixel; Bad point is meant and faces the territory pixel has notable difference on gray scale pixel around it, when therefore adopting formula (1) to calculate, in order to guarantee accuracy, facing the territory selects unsuitable excessive, consider detection speed simultaneously and image evenly can be divided equally, the image of 1024*1024 is divided into several regions according to the size of 16*16 pixel.

2) at each zone, determine its gray scale maximum, minimum value, the mean value of the pixel grey scale after maximum, the minimum gradation value is removed in calculating; When calculating average gray, must eliminate the influence of extreme point (potential " bad point "): for example to mean value, in the zone of image 16*16, if the gray value of a bad point is 1000, all the other each pixel averages are 10, and then this area pixel gray value is 13.9, like this for 15% detection threshold, not only bad point (gray value is 1000) has reached detection threshold, and other each points also might surpass detection threshold; The reason that causes this result is exactly the influence of extreme point (potential " bad point ") to mean value.Therefore when computation of mean values, at first obtain the maximum and the minimum value of this area grayscale, the mean value of all pixels after calculating removal maximum and the minimum value.

3) utilize formula (1), set less detection threshold, at each zone, Threshold detection is carried out in pointwise;

4) to pixel greater than thresholding, its position of mark.

It is as follows to utilize formula (2), (3) to carry out the step of Threshold detection:

1) carry out medium filtering according to formula (2), the medium filtering parameter is (2n+1) * (2n+1), for the image of 1024*1024, and 1＜n＜511;

2), medium filtering and original image are carried out subtraction process, and the result is taken absolute value according to formula (3);

3) set less detection threshold, carry out Threshold detection;

4) to pixel greater than thresholding, its position of mark.

Step 3: vision-based detection

" vision-based detection " described in the inventive method is to utilize human-eye visual characteristic, adopts the display mode of " the mark mark is relevant ", the small-signal in the detection background.So-called " the mark mark is relevant " is meant the multiframe sequence image that will collect under the situation of not doing any noise reduction process, and every two field picture nominated bank is shown successively side by side, and bad like this point can become bright (secretly) line by a bit.Detect bad point by such processing.

The setting of detection threshold mainly is to determine by practical experience.Whether pixel is bad point, finally needs the observer to pass through to observe judgement.Though the relative mean value of some gray values of pixel points can rise and fall to some extent, as long as by observing the not obvious adjacent pixels gray scale on every side that is different from of its gray scale, such pixel can not adjudicated and is bad point.

Adopt the relevant mode of mark mark to show the corresponding row of bad point in the image sequence of gathering, utilize multirow to show the effect of bringing, can very clearly find bad point.The bad point that Threshold detection is obtained carries out vision-based detection, investigate the relevant result displayed of mark mark and whether have " bright (secretly) line ", Fig. 2 is that the relevant display result of mark mark exists the schematic diagram of " bright line ", as long as there is the bad point of place's Threshold detection " bright (secretly) line " not occur, illustrate that then detection threshold is too small, should improve detection threshold, carry out Threshold detection again, the amplitude that each step-length is adjusted is unsuitable excessive.

Step 4: determine detection threshold

Repeating step two and three is carrying out till vision-based detection still is judged as bad point until the bad point that all Threshold detection obtain, and the thresholding of this moment is rational detection threshold.

The invention provides a kind of processing method that Threshold detection and vision-based detection are combined, both used the automaticity of Threshold detection, taken into account the accuracy of vision-based detection again, its advantage is simple and practical, has guaranteed accuracy simultaneously.

Description of drawings

Fig. 1 is the bad point of the inventive method testing process figure;

Fig. 2 is the display result of " the mark mark is relevant " of bad some corresponding row in the inventive method;

Fig. 3 is for being used for the image that bad point detects in the inventive method;

Fig. 4 amplifies for the part that is used for bad some image that detects in the inventive method;

Fig. 5 is the display result of bad point in the inventive method (500,62) " the mark mark is relevant ";

Fig. 6 is the display result of bad point in the inventive method (492,24) " the mark mark is relevant ";

Fig. 7 is the result behind the bad point calibration in the inventive method;

Fig. 8 is local amplification behind the bad point calibration in the inventive method;

Fig. 9 is the display effect of " the mark mark is relevant " behind the bad point calibration in the inventive method.

Embodiment

In conjunction with the accompanying drawings, be that 1024*1024 is an example with the image size, concrete implementation step of the present invention is as follows:

Utilize formula (1) to detect bad point:

1) obtain detected image:

Bad point is a relative background " bright spot " perhaps " dim spot ": detect " bright spot ", need to gather the darker image of background and detect; On the contrary, detect " dim spot ", the image background of collection should be brighter.Bad point has on gray value than big difference with other pixels that face the territory, and therefore should be noted that 2 points when going bad some detection: 1, the removal of images background fluctuation is to avg (x, influence y); 2, using the higher image of signal to noise ratio to go bad a little detects.

Present embodiment detects at " bright spot ", and " dim spot " detection case is similar with it.

At first close the CCD camera aperture, the removal of images background fluctuation; Gather multiple image, utilize the method for multiple image stack to obtain being used for the higher image of signal to noise ratio that bad point detects.Gather 98 frame image sequence under these conditions, carry out multi-frame mean and handle, the result as shown in Figure 3.By visual observation, can find wherein to have a bad point in the two field picture top.Local amplification effect as shown in Figure 4.

2) set a detection threshold 10%, utilize formula (1) to carry out Threshold detection, judgement pixel (500,62) is a bright spot;

3) the bad point that Threshold detection is obtained carries out vision-based detection, investigate the relevant result displayed of its mark mark as shown in Figure 5, obvious " bright line " that does not have similar Fig. 2 among Fig. 5, though illustrating that this pixel grey scale departs from faces territory mean value, but visual effect can not be adjudicated and is bad point, therefore judge that 10% detection threshold is less than normal, should adjust detection threshold, the amplitude that each step-length is adjusted is unsuitable excessive;

4) repeating step 2), 3), adjust detection threshold, the result is adjudicated by vision, determine that at last detection threshold is 30% comparatively suitable, the bad point that this moment, all Threshold detection obtained still is judged as bad point carrying out vision-based detection, so this thresholding is rational detection threshold.

Utilize formula (2), (3) to detect bad point:

1) obtain detected image:

Method is the same.

2) setting detection threshold is 5, by the medium filtering of formula (2) selection 5*5, carries out the medium filtering of 5*5, image and source images after by formula (3) medium filtering being handled carry out subtraction process, and carry out Threshold detection after the result taken absolute value, judgement pixel (492,24) is a bright spot;

3) the bad point that Threshold detection is obtained carries out vision-based detection, and the relevant display result of its mark mark as shown in Figure 6.Do not have similar Fig. 4 obvious " bright line " among Fig. 6, though this grey scale pixel value is described greater than detection threshold, visual effect can not be adjudicated and is bad point.Illustrate that 5 detection threshold is less than normal, adjust detection threshold,

4) repeating step 2), 3), adjust detection threshold, the result is adjudicated by vision, determine that at last detection threshold is 20 comparatively suitable, the bad point that this moment, all Threshold detection obtained still is judged as bad point carrying out vision-based detection, so this thresholding is rational detection threshold.

Behind definite badly some position, adopt its pixel value at the bad some place of non-bad some pixel weighted sum replacement of neighborhood on every side, carry out bad point calibration.

Result after Fig. 3 detected and proofread and correct locally amplifies the result as shown in Figure 8 as shown in Figure 7, and Fig. 9 is for proofreading and correct the relevant display result of back mark mark.

Comparison diagram 4 and Fig. 8, visible the inventive method detects validity to bad point.Comparison diagram 2 and Fig. 9, Fig. 9 do not exist the tangible similar Fig. 2's of vision " bright line ", the correctness of bearing calibration is described.

Claims (5)

1. the automatic detection threshold of the bad point of medical image is determined method, and the bad point that the inventive method will be referred at first is defined as:

$\mathrm{\α}(x,y)=\frac{\mathrm{gray}(x,y)-\mathrm{avg}(x,y)}{\mathrm{avg}(x,y)}$

(x is that ((x is that (x is the mean value of the image-region pixel grey scale at center y), and computing formula is with pixel y) to avg to pixel for x, gray value y) y) to gray in the formula $\mathrm{avg}(x,y)=\frac{1}{N}\underset{(i,j)\∈S}{\mathrm{\Σ}}\mathrm{gray}(i,j),$ Wherein S is that (x y) is the image-region at center, and (i j) is pixel (i to gray with pixel, j) gradation of image, N are that the whole zone removal of image pixel grey scale is the number of pixels of two pixels of extreme value, if | a (x, y) |＞detection threshold, (x y) is bad point then to adjudicate this point; Perhaps by formula I _Median=Median (I _Originall), I (x, y)=| I _Median-I _Originall| judge I in the formula _OriginallBe original image, I _MedianBe median-filtered result, Median () expression medium filtering is handled, if I (x, y)＞detection threshold, (x y) is bad point then to adjudicate this point; It is characterized in that the automatic detection threshold of its definite bad point comprises collection bad some detected image, thresholding detection, vision-based detection, determines four steps of detection threshold.

2. determine method according to the automatic detection threshold of the bad point of the described medical image of claim 1, it is characterized in that at first closing the CCD camera aperture, the removal of images background fluctuation in the bad point of the described collection detected image step; Gather multiple image, utilize the method for multiple image stack to obtain being used for the higher image of signal to noise ratio that bad point detects.

3. determine method according to the automatic detection threshold of the bad point of the described medical image of claim 1, it is characterized in that described thresholding detects in the step, set a less detection threshold, according to the formula of determining bad point, carry out Threshold detection, to its position of pixel mark greater than thresholding.

4. determine method according to the automatic detection threshold of the bad point of the described medical image of claim 1, it is characterized in that " vision-based detection " described in the inventive method, be to utilize the people to follow visual characteristic, adopt the display mode of " the mark mark is relevant ", the small-signal in the detection background; So-called " the mark mark is relevant " is meant the multiframe sequence image that will collect under the situation of not doing any noise reduction process, and every two field picture nominated bank is shown successively side by side, and bad like this point can become bright (secretly) line by a bit; The bad point that Threshold detection is obtained carries out vision-based detection, investigate the relevant result displayed of mark mark and whether have " bright (secretly) line ", as long as there is the bad point of place's Threshold detection " bright (secretly) line " not occur, illustrate that then detection threshold is too small, should improve detection threshold, again carry out Threshold detection, the amplitude that each step-length is adjusted is unsuitable excessive.

5. determine method according to the automatic detection threshold of the bad point of the described medical image of claim 1, it is characterized in that determining that detection threshold is that the repetition thresholding detects, the vision-based detection step, carrying out till vision-based detection still is judged as bad point until the bad point that all Threshold detection obtain, the thresholding of this moment is rational detection threshold.

Images