Method and apparatus for calibration and correction of gray levels in images
A method and apparatus for calibrating an apparatus that acquires a sequence of radiographic images and correcting images of an object under observation. For each image of a sequence acquired by the apparatus and for a given frequency of acquisition of the sequence, the apparatus is calibrated by determining the value of the variation of a mean of gray levels in at least one zone of interest of the current image of at least one calibration device, the variation being determined relative to the mean gray level of the first image of the sequence in each zone of interest. The determination of the variation is reiterated for a series of images sequences acquired using calibration devices resulting in first images of mean gray levels different from one sequence to another. Each image of an image sequence of the object under observation is corrected, comprising zones of observation having different gray levels by subtracting from the current image the variation of one gray level relative to the first image of the object, the subtraction being a function of the gray level considered from each zone of observation.
This application claims the benefit of a priority under 35 USC 119(a)-(d) to French Patent Application No. 03 01797 filed Feb. 14, 2003, the entire contents of which are hereby incorporated by reference.
BACKGROUND OF THE INVENTIONThe present invention relates to a method apparatus for calibrating and correction of gray levels in images. In particular, the present invention is directed to acquiring a sequence of radiographic images and correction images of an object under observation. More particularly, the present invention relates to a method and apparatus for an acquiring sequence of radiographic images and calibration and correction of images of an object under observation by subtracting from each image of the sequence the spurious variation of a gray level between the images of the same sequence. The present invention can be particularly used in the medical field, such as, for example, in mammography and in the detection of cancerous tumors.
A known radiographic apparatus is comprises a console, means for providing a beam of radiation in the direction of means for detection. The means for detecting receives the radiation after passing through an object under observation placed in the observation space arranged in the beam, between on the means for providing the beam of radiation and on the means for detection. The known apparatus also comprises means for processing enabling acquiring and processing a sequence of images of an object sent from the means for detection.
The difference in absorption of the radiation by the different parts of the object under observation enables obtaining information on the composition of the object. The image formed on the means for detection comprises different gray levels, from which information can be derived. Thus, if the object under observation is a human body part, for example, the bones will appear clearly on the image acquired by the means for detection and are distinctly separate from the part formed by the muscles.
Generally, at the time of acquisition of a sequence of radiographic images, there is always a variation in the gray levels found between the successive images of the sequence. In a succession of images acquired the means for detection measures an increase in the variation of the gray level. This phenomenon is due to a persistence or remanence of the radiographic information from one image to the other, which causes the gray level to vary between the images. The variation of the gray level in the sequence depends on the thickness and the composition of the object that is being observed. Thus, for an object observed having a first thickness, will have curve of the variation versus time different than a curve for an object having a second thickness different from the first thickness.
The variation of the gray level is due principally to the trapping of charges in the photodiodes of the means for detection. The variation can also have a number of causes. It can be a question especially of an increase in temperature of the different elements of the apparatus.
In all cases, the variation of the gray levels from one image to another perturbs the measurements acquired by the apparatus. In the case of a large remanence that is, of very considerable variations, the quality and the interpretation of the images acquired may deteriorate considerably. There may be the appearance of a “ghost” or multiple images; that is, superimposing of images acquired previously onto a new image of an object. When the remanence is weaker and does not cause the appearance of ghosts, the measures made are similarly distorted by the spurious variation of the gray level from one image to the other. The variations can be of the same order of magnitude as the dynamics in gray levels of the signal that one wants to detect.
Certain methods enable elimination of this remanence in special applications using special devices. In certain applications, it is possible to determine a law according to which the remanence diminishes over time and to subtract the remanence in the images of the sequence. Consequently, at least one black measurement is made between two acquisition instants corresponding to gray levels supplied by the means for detection in the absence of exposure by radiation. The black measurement enables determining the value of the remanence for that instant. As a result, by virtue of the knowledge of the values of remanence at a given first time the law of decrease of remanence can be deduced. From the law of decrease of remanence between two acquisitions, the value of the remanence at given second time can be deduced of the following acquired image and thus correct the images acquired in a sequence.
This method of correcting sequences of images presumes the use of devices for measuring blacks. Such devices are not always available or accessible on prior art radiographic devices. Consequently, in the majority of cases, it is practically impossible to be able to correct the variation in gray level in a sequence of successive images.
BRIEF DESCRIPTION OF THE INVENTIONAn embodiment of the invention provides a method and an apparatus for correction of gray levels in images. An embodiment of the invention is directed to correcting the remanence in a sequence of radiographic images. The method and the apparatus are able to eliminate the utilization of black measurement devices.
An embodiment of the invention provides a method and an apparatus for calibration a device capable of acquiring a sequence of radiographic images. Calibration is done so as to be able to correct the effects of the variation of gray level in a sequence of radiographic images of an object under observation.
An embodiment of the invention provides a method and an apparatus for calibration and correction of radiographic images applicable to all repeatable variation phenomena; that is, variation phenomena that repeat from one measurement of an acquisition sequence to the next when maintaining identical acquisition phenomena from one measurement of an acquisition sequence to the next.
An embodiment of the invention provides a method for calibrating an apparatus capable of acquiring a sequence of radiographic images and correcting images of an object under observation comprising:
-
- for each image of a sequence acquired by the apparatus and for a given frequency of acquisition of the sequence, the apparatus is calibrated by determining the value of the variation of a mean of gray levels in at least one zone of interest of the current image of at least one calibration device, the variation being determined relative to the mean gray level of the first image of the sequence in each zone of interest;
- the determination of the variation is reiterated for a series of images sequences acquired using calibration devices resulting in first images of mean gray levels different from one sequence to another;
- each image of an image sequence of the object under observation is corrected, comprising zones of observation having different gray levels by subtracting from the current image the variation of one gray level relative to the first image of the object, the subtraction being a function of the gray level considered from each zone of observation.
The invention similarly relates to an apparatus capable of acquiring a sequence of radiographic images implementing a method according to an embodiment of the invention.
BRIEF DESCRIPTION OF THE DRAWINGSThe embodiments of the invention will better understood from the following description, which is purely illustrative and non-limiting, and which should be read with reference to the drawings annexed hereto, wherein:
An example of a known radiographic apparatus is represented schematically in
As noted above, the difference in absorption of the radiation by the different parts of the object under observation enables obtaining information on the composition of the object. In fact, an image formed on the detector 2 comprises different gray levels, from which information can be derived. Thus, if the object under observation is a human body part, for example, the bones will appear clearly on the image acquired by the detector and are distinctly separate from the part formed by the muscles.
As noted above, generally, at the time of acquisition of a sequence of radiographic images, there is always a variation in the gray levels found between the successive images of the sequence. Thus,
As noted above, the variation L of the gray level is due principally to the trapping of charges in the photodiodes of the detector. The variation L can also have a number of causes. It can be a question especially of an increase in temperature of the different elements of the device.
As noted above, in all cases, the variation of the gray levels from one image to another perturbs the measurements acquired by the device. In fact, in the case of a large remanence that is, of very considerable variations, the quality and the interpretation of the images acquired may deteriorate considerably. In fact, there may be the appearance of a “ghost” or multiple images; that is, superimposing of images acquired previously onto a new image of an object. When the remanence is weaker and does not cause the appearance of ghosts, the measures made are similarly distorted by the spurious variation of the gray level from one image to the other. The variations can be of the same order of magnitude as the dynamics in gray levels of the signal that one wants to detect.
As noted above, in fact, in certain applications, it is possible to determine the law according to which the remanence diminishes over time and to subtract said remanence in the images of the sequence.
An embodiment of the invention is a calibration method of an apparatus capable of acquiring a sequence of radiographic images and correction of images of an object under observation enabling correcting the unwanted effects of the gray level variations in a sequence of radiographic images.
The embodiments of the method can be considered according to two approaches. In a first approach, the calibration step is done prior to acquisition of the sequence of images of the object under observation. The second approach allows performance of the calibration at the same time as acquisition of the images of the object under observation.
In all of the figures, similar elements or steps are referenced using identical numbers.
An embodiment of a method for implementing the first approach is represented schematically in
The calibration part is described as follows. At the time of step (a), using the radiographic apparatus a sequence of images #1, #2, . . . #N are acquired. Acquisition of the sequence is done by observing a calibration device 40. The calibration device 40 is positioned in the zone of observation over the detector of the apparatus and covers the greater part of the detector surface.
This sequence is acquired for a given acquisition frequency. For example, the acquisition frequency can be acquisition of an image every 30 seconds or acquisition of an image every 60 seconds. The images can also be acquired at irregular time intervals. The acquisition frequency of the calibration sequence is preferentially the same as the acquisition frequency that is to be used for acquiring the sequence of images of the object under observation.
At the time of acquisition and for each image a mean gray level is determined on a selected homogeneous zone of interest 41, as shown in
A series of acquisition of sequences is then done using calibration devices 40 having different thicknesses. Thus, the mean gray levels are obtained regarding the zone of interest 41 different between each sequence of the series.
As noted above, the mean gray level evolves as a function of time; that is, it varies between successive images of the same sequence. Thus, the average gray levels in the zone of interest 41 of each image of all the sequences of the series are made. Then, for each nth image C(n) of an sequence acquired by the apparatus and for the zone of interest 41, using means for processing in each apparatus, the value of the mean variation of gray level between the current image and the first image is determined. L(n) records this difference. It is equal to:
L(n)={overscore (C(n))}−{overscore (C(n))} (1)
One then proceeds to the step (b) in
L(n)={overscore (C(1))}×[α(n) {overscore (C(1))}+β(n)] (2)
wherein α(n) and β(n) are the coefficients of the linear regression of the curves plotted at the time of step (b) and calculated at the time of a step (c).
Step (c) shows that using the means for processing included in the apparatus, the coefficients of regression α(i) and β(i) corresponding respectively to the director coefficients and at the ordinate to the origin of each line of image i is determined. These coefficients are stored in the means for providing a memory of the apparatus included in the means for processing 8 or in means arranged outside of the apparatus. The operation of step (c) ends the calibration.
Each line depends on the one hand on the priority of the nth image in the sequence as indicated by the superscript (n) and, on the other hand, on the frequency of acquisition of the sequence. In other words, and as can be confirmed in
The set of curves are plotted corresponding to all of the images of the calibration sequence.
The correction part for the image is described as follows.
In order to apply a correction to the pixel of coordinates (i,j) of an uncorrected nth image Y(n), the mean gray level
is measured in a zone of observation R(i,j) centered on the point (i,j) of the first image, for example, zone 43 in
The method for correction is then applied to the image Y(n) by subtracting from the current image the variation of a gray level relative to the first image of the object.
at the time of step (d) and using the calibration data, the value of the variation L(n) can be reached which is a function of
So, it is sufficient to subtract this value from the value of the current gray level.
A measurement is thus made of the mean gray level in a plurality of zones of observation, for example, similarly in the zone of observation 44 in
Mathematically, a corrected nth image
is obtained using the formula:
By applying such a correction to the set of images of the sequence, a sequence is obtained, in which all of the variations of gray level have been suppressed or at least reduced.
In the hereinbefore described equations, the measure of the mean level of gray around the pixel (i,j) has been used. A plurality of alternatives is possible as a function of the gray level being considered.
The value of the medial gray level in each zone of observation can also be determined. Thus, the mean gray level is no longer considered. The medial value of a series is the value situated in the middle of the series of values arranged in ascending or descending order.
The use of the median avoids the affect of the measured aberrant gray levels on the value of gray levels taken into account in the correction step. In effect, the extreme values of the series have no influence on the calculation of the median value. Such aberrant values can be measured in the zone of the object 42 having a strong thickness gradient and around an abrupt change in thickness of the object under observation.
The median gray level instead of mean gray level is considered for the object having such thickness gradients. The calibration steps and the equations (1) and (2) remain the same. Equation (3) is slightly modified, because the median gray level is being applied instead of the mean gray level.
The calculation of the median generally slightly increases the processing time of information done in the processing means of the apparatus.
According to second alternative, the mean gray level can also be replaced by the value of
of gray level of the pixel (i,j). Of course, the value of a zone of interest R(i,j) (43 and/or 44, for example) is no longer averaged and the median value is not calculated.
As for the first alternative, the steps of the method remain the same and only equation (3) is slightly modified, because it is applied to the value
of the gray level at the pixel (i,j).
The second alternative allows, as did the first alternative, having a good estimation of the variation of gray level in pixels situation near a zone in which the thickness of the object under observation rapidly varies. In contrast to the first alternative, it allows a reduction in the processing time. However, determination of the value of the variation of the gray level is less precise, because there is an enhancement of the quantum noise effect. The effects of quantum noise are significantly reduced by averaging or calculation of the median.
In all of the aforesaid embodiments, it is presumed that the variation of the gray level is homogeneous over the entire surface of the detector of the apparatus. The coefficients α(n) and β(n) are determined in fact only in one single zone of interest 41.
An adaptive approach of the first approach is described as follows. The adaptive approach takes into account the inhomogeneity of the variations in gray level over the surface of the detector of the apparatus. It appears that the detectors are not perfect and that there is a disparity of variations of gray level depending on the position on the detector.
The steps of an adaptive approach are represented schematically in
Thus, at the time of step (a) in
is calculated using the equation:
-
- (which is similar to equation (1)), wherein each Rc(k,l),(k,l) εΩ represents a zone of interest (41) centered on the pixel having the coordinates (k,l) and where Ω is the set of pixels on which the zones of interest 41 in
FIG. 5 are centered. Similarly, LΩ(n)/{overscore (C(1) )}is calculated in that it is plotted as a function of {overscore (C(1))}.
- (which is similar to equation (1)), wherein each Rc(k,l),(k,l) εΩ represents a zone of interest (41) centered on the pixel having the coordinates (k,l) and where Ω is the set of pixels on which the zones of interest 41 in
Thus, at the time of step (b) it is not a graph of a set of straight lines corresponding to the different images that is obtained but a set of graphs, each one corresponding to a zone of interest 41. Thus, one has as many graphs as zones of interest on the surface of the calibration device 40, i.e., a set of graphs 411, 412, . . . 41N corresponding to the set of measurements in the N zones of interest 41.
Step (c) approximates the curves obtained by the straight lines that were defined for each graph, the coefficients
henceforth dependent on the zone of interest 41 for which they have been calculated.
At step (d), a sequence of images of an object 42 under observation are acquired. The zones of observation 43 and 44 are determined for which one wishes to effect a correction. Thus the position of the zones of observation is marked relative to the different zones of interest 41.
In step (e) in order to determine the variation of gray level of each zone of observation, the values of
are taken, which have been determined for the corresponding zone of interest 41. Then the correction is applied using an equation similar to equation (3). The corrected gray level
becomes:
(which is similar to formula (3)) wherein
is the gray level before correction and
are the coefficients determined for each zone of interest and applied to the zone of observation centered on the pixel (i,j). In this fashion, a corrected image
is obtained in step (f).
Several alternatives of the adaptive approach are possible. A first possible alternative consists of interpolating the
coefficients of zones of observation arranged outside of the zone of interest. Thus, one could have coefficients and
for any zone of observation. Other alternatives are also possible and utilize the median gray level value or the gray level values on a pixel instead of considering the mean gray level in a zone of observation.
As described above, an embodiment of the method for the first approach considers a calibration before acquiring the image sequence of the object under observation. An embodiment of the method for a second approach does the calibration at the same time as acquisition of the image sequence of the object.
In the second approach of a method according to an embodiment of the invention will now be described the calibration steps take place in the same time as acquisition of the images of the object under observation. The steps of the second approach are represented schematically in
In
In the second approach, the calibration device 40 comprises at least two zones of interest that have mean gray levels that are different from one zone to the other for each image. This is what has been represented in
Thus, at the time of an identical acquisition of a sequence of images of the calibration device 40, there will be a plurality of points of reference enabling generating the plotting according to
Thus, at the time of an identical acquisition of a sequence of images of the calibration device 40, there will be a plurality of points of reference enabling generating the plotting according to
In step (b) the curves of these relative values are plotted with respect to the mean gray level in the first image as a function of the gray level in the first image.
In step (c) the different curves of the variations of gray level are approximated as a function of the gray level of the first image by the straight lines and the coefficients of the representative function are calculated.
In step (e) the calibration of steps (a), (b), (c) are used. The calibration enables calculation of the variation for each image of the object 42 of the sequence. Thus, in (f), a corrected series of images is obtained.
In an understanding of this approach, the following can be recognized. Firstly, the most precise correction of the gray level variations is obtained. The remanence is determined directly using the images of the sequence of the object under observation. It is then reasonably certain that there are no differences in behavior of the apparatus between the calibration sequences and acquisition of the images of the object under observation. Secondly, it is generally not necessary to do an entire series of measurements in order to obtain the calibration curves. The calibration is done directly, at the same time as acquisition of the images of the object under observation. Thus, there is a considerable time savings for the operator of the apparatus.
It is not easy to apply an adaptive approach to the second approach: the presence of the object under observation, calibration of the entirety of the acquisition field cannot be done. Thus, according to a first alternative of the second approach, a spatial model is used in order to take into account any disparities of variation of gray level as a function of the position on the detector. The steps of such a variant are represented schematically in
Steps (a), (b), (c) and (e) remain the same as those of
Thus, in order to apply a correction to a pixel (i,j) of image Y(n), the mean gray level
is measured in a zone of observation 44, for example, centered around the point (i,j). After calibration, the gray level Zi,j(n) of the pixel (i,j) becomes:
-
- (which is similar to formula (3)), wherein
is the gray level before calibration α(n) and β(n) are the linear regression coefficients of approximation of the variation of gray level and λi,j is a measured or modeled gain factor for the pixel (i,j).
- (which is similar to formula (3)), wherein
In order to obtain a measured λi,j gain factor one can, for example, utilize a calibration for a plurality of zones of interest as in the adaptive approach described for the first approach. The different values of λi,j will be thus entered in the processing means of the apparatus and applied at the time of the correction step.
In order to model the λi,j gain factor, one can utilize the inverse of the apparatus gain factor. The Kij gain factor of the apparatus compensates the inhomogeneity of illumination of the detector by the emitter and the inhomogeneity of the response of the photodiodes of the detector.
In this fashion, an image corrected in gain is observed on the detector of the form: Kij (Im+remanence), wherein Im is the image of the object uncorrected in gain and inhomogeneous. The part (Kij Im) is a homogenous image. Thus, the inhomogeneous remanence that is observed on the detector is: (Kij remanence). The inhomogeneity of the remanence observed is due to the Kij gain factor. Thus, in order to find a homogeneous remanence, the λi,j factor is applied to the term (Kij remanence), such that: λi,j=1/Ki,j. It is the homogeneous remanence obtained after this multiplication that can be subtracted from the current image.
According to a second possible alternative of the second approach, one compensates for the fact that few of the values for the plotting of the calibration curve can be acquired. As the calibration data are acquired in the same time as the image observation data, there is only one variation value per image for a given gray level.
In order to increase the number of measurement points for improving the precision of the calculation of the regression coefficients α(i) and β(i) of the approximation, it is possible to combine the information from a plurality of successive acquisition sequences of different objects under observation.
The alternative of
Another alternative of the second approach compensates for the fact that few functional points are available for the plotting of the calibration curves. Actually, due to the fact of the presence of the object under observation in the acquisition field of the detector, calibration devices of larger dimensions cannot be arranged in the field. Consequently, the devices no not have a large number of different zones of attenuation characteristics. The present variant utilizes different calibration devices from one sequence to another in order to augment the number of functional points. The steps of such an alternative are represented schematically in
In the alternatives of
Other variants of the second approach are also possible and they utilize gray level median values or the gray level values on a pixel, instead of considering the mean gray level in a zone of observation.
The introduction of the gain factor λi,j during the correction step can be similarly applied to the first approach.
In the developments of the above, the relative estimated value L(n)/{overscore (C(1) )}in the nth image in function of the mean gray level in the first image {overscore (C(1) )}can be approximated by a straight line. Of course, in the above developments, the visible graphic representation of
The embodiments of the invention is therefore complemented by the following features, taken singly or in any technically possible combination thereof:
-
- a graphical representation is approximated having respectively as its ordinate and abscissa;
- the ratio having in the numerator the mean gray level variation of the current image of the sequence of the calibration device relative to the mean gray level of the first image and in the denominator the mean gray level of the first image, and
- the different mean gray levels of each first image;
- a function for which the characteristics are known;
- the function is a straight line;
- the gray level that is corrected in each zone of each image of the object under observation is the mean gray level;
- the gray level that is corrected in each zone of each image of the object under observation is the median gray level;
- the gray level of at least one selected pixel in each image of the object under observation is corrected;
- the calibration takes place before the acquisition of the sequence of images of the object under observation;
- each mean gray level value of the series of sequences of the calibration is given by the observation of at least one calibration plate of defined thickness constituting each calibration device, its thickness change from one series to the other;
- the average gray level is determined on a plurality of zones of interest simultaneously at the time of calibration;
- the subtraction at the time of the correction step depends on the one hand on the position of the zone of observation relative to each zone of interest;
- the value subtracted from each image of the sequence of images of the object is a function on the one hand of the position of the zone of observation and on the other hand a function of the defined spatial gain;
- the calibration takes place during the acquisition of the sequence of images of the object under observation;
- each calibration device is placed in a field of acquisition of the apparatus comprising also the object under observation;
- each calibration device comprises at least two zones of interest having a mean gray level different from one zone to another for each image;
- the value subtracted from each image of the sequence of images of the object is a function on the one hand of the position of the zone of observation and on the other of the function of the spatial gain of the apparatus;
- the calibration measurements of at least two acquisitions of successive sequences are combined, and
- the calibration measurements of at least two acquisitions of successive sequences having different calibration devices are combined.
One skilled in the art may propose or make various modifications to the way/structure and/or function and/or result and/or steps of the disclosed embodiments without departing from the scope and extant of protection.
Claims
1. A method for calibrating an apparatus capable of acquiring a sequence of radiographic images and correcting images of an object under observation, comprising:
- for each image of a sequence acquired by the apparatus and for a given frequency of acquisition of the sequence, the apparatus is calibrated by determining the value of the variation of a mean of gray levels in at least one zone of interest of the current image of at least one calibration device, the variation being determined relative to the mean gray level of the first image of the sequence in each zone of interest;
- the determination of the variation is reiterated for a series of images sequences acquired using calibration devices resulting in first images of mean gray levels different from one sequence to another; and
- each image of an image sequence of the object under observation is corrected, comprising zones of observation having different gray levels by subtracting from the current image the variation of one gray level relative to the first image of the object, the subtraction being a function of the gray level considered from each zone of observation.
2. The method according to claim 1 wherein a graphic representation is approximated having respective for its ordinate and abscissa:
- the ratio having in the numerator the mean gray level variation of the current image of the sequence of the calibration device relative to the mean gray level of the first image and in the denominator the mean gray level of the first image; and
- the different gray levels of each first image;
- by a function for which the characteristics are known.
3. The method according to claim 2 wherein the function is a straight line.
4. The method according to claim 1 wherein the gray level that is corrected in each zone of each image of the object under observation is the mean gray level.
5. The method according to claim 2 wherein the gray level that is corrected in each zone of each image of the object under observation is the mean gray level.
6. The method according to claim 3 wherein the gray level that is corrected in each zone of each image of the object under observation is the mean gray level.
7. The method according to claim 1 wherein the gray level that is corrected in each zone of each image of the object under observation is the median gray level.
8. The method according to claim 2 wherein the gray level that is corrected in each zone of each image of the object under observation is the median gray level.
9. The method according to claim 3 wherein the gray level that is corrected in each zone of each image of the object under observation is the median gray level.
10. The method according to claim 1 wherein the gray level of at least one pixel chose in each image of the object under observation is corrected.
11. The method according to claim 2 wherein the gray level of at least one pixel chose in each image of the object under observation is corrected.
12. The method according to claim 2 wherein the gray level of at least one pixel chose in each image of the object under observation is corrected.
13. The method according to claim 1 wherein the calibration is done before the acquisition of the image sequence of the object under observation.
14. The method according to claim 2 wherein the calibration is done before the acquisition of the image sequence of the object under observation.
15. The method according to claim 3 wherein the calibration is done before the acquisition of the image sequence of the object under observation.
16. The method according to claim 4 wherein the calibration is done before the acquisition of the image sequence of the object under observation.
17. The method according to claim 7 wherein the calibration is done before the acquisition of the image sequence of the object under observation.
18. The method according to claim 10 wherein the calibration is done before the acquisition of the image sequence of the object under observation.
19. The method according to claim 1 wherein each mean gray level value of the series of sequences of the calibration is given by the observation of at least one calibration plate of a defined thickness comprising each calibration device, its thickness changing from one series to another.
20. The method according to claim 2 wherein each mean gray level value of the series of sequences of the calibration is given by the observation of at least one calibration plate of a defined thickness comprising each calibration device, its thickness changing from one series to another.
21. The method according to claim 3 wherein each mean gray level value of the series of sequences of the calibration is given by the observation of at least one calibration plate of a defined thickness comprising each calibration device, its thickness changing from one series to another.
22. The method according to claim 4 wherein each mean gray level value of the series of sequences of the calibration is given by the observation of at least one calibration plate of a defined thickness comprising each calibration device, its thickness changing from one series to another.
23. The method according to claim 7 wherein each mean gray level value of the series of sequences of the calibration is given by the observation of at least one calibration plate of a defined thickness comprising each calibration device, its thickness changing from one series to another.
24. The method according to claim 10 wherein each mean gray level value of the series of sequences of the calibration is given by the observation of at least one calibration plate of a defined thickness comprising each calibration device, its thickness changing from one series to another.
25. The method according to claim 13 wherein each mean gray level value of the series of sequences of the calibration is given by the observation of at least one calibration plate of a defined thickness comprising each calibration device, its thickness changing from one series to another.
26. The method according to claim 13 wherein the mean gray level is determined using a plurality of zones of interest simultaneously at the time of calibration.
27. The method according to claim 19 wherein the mean gray level is determined using a plurality of zones of interest simultaneously at the time of calibration.
28. The method according to claim 19 wherein the subtraction of the correction step depends on the one hand on the gray level in each zone of observation and on the other hand on the position of the zone of observation relative to each zone of interest.
29. The method according to claim 26 wherein the mean gray level is determined using a plurality of zones of interest simultaneously at the time of calibration.
30. The method according to claim 13 wherein the value subtracted from each image of the image sequence of the object is a function on the one hand of the position of the zone of observation and on the other hand of a defined spatial gain function.
31. The method according to claim 19 wherein the value subtracted from each image of the image sequence of the object is a function on the one hand of the position of the zone of observation and on the other hand of a defined spatial gain function.
32. The method according to claim 1 wherein the calibration is done during the acquisition of the image sequence of the object under observation.
33. The method according to claim 2 wherein the calibration is done during the acquisition of the image sequence of the object under observation.
34. The method according to claim 3 wherein the calibration is done during the acquisition of the image sequence of the object under observation.
35. The method according to claim 4 wherein the calibration is done during the acquisition of the image sequence of the object under observation.
36. The method according to claim 7 wherein the calibration is done during the acquisition of the image sequence of the object under observation.
37. The method according to claim 10 wherein the calibration is done during the acquisition of the image sequence of the object under observation.
38. The method according to claim 1 wherein each calibration device is placed in a field of acquisition of the apparatus also comprising the object under observation.
39. The method according to claim 2 wherein each calibration device is placed in a field of acquisition of the apparatus also comprising the object under observation.
40. The method according to claim 3 wherein each calibration device is placed in a field of acquisition of the apparatus also comprising the object under observation.
41. The method according to claim 4 wherein each calibration device is placed in a field of acquisition of the apparatus also comprising the object under observation.
42. The method according to claim 7 wherein each calibration device is placed in a field of acquisition of the apparatus also comprising the object under observation.
43. The method according to claim 10 wherein each calibration device is placed in a field of acquisition of the apparatus also comprising the object under observation.
44. The method according to claim 32 wherein each calibration device is placed in a field of acquisition of the apparatus also comprising the object under observation.
45. The method according to claim 1 wherein each calibration device comprises at least two zones of interest having a mean gray level different from one zone to another for each image.
46. The method according to claim 2 wherein each calibration device comprises at least two zones of interest having a mean gray level different from one zone to another for each image.
47. The method according to claim 3 wherein each calibration device comprises at least two zones of interest having a mean gray level different from one zone to another for each image.
48. The method according to claim 4 wherein each calibration device comprises at least two zones of interest having a mean gray level different from one zone to another for each image.
49. The method according to claim 7 wherein each calibration device comprises at least two zones of interest having a mean gray level different from one zone to another for each image.
50. The method according to claim 10 wherein each calibration device comprises at least two zones of interest having a mean gray level different from one zone to another for each image.
51. The method according to claim 32 wherein each calibration device comprises at least two zones of interest having a mean gray level different from one zone to another for each image.
52. The method according to claim 38 wherein each calibration device comprises at least two zones of interest having a mean gray level different from one zone to another for each image.
53. The method according to claim 32 wherein the value subtracted from each image of the image sequence of the object is a function on the one hand of the observation zone and on the other hand a function of the spatial gain of the apparatus.
54. The method according to claim 38 wherein the value subtracted from each image of the image sequence of the object is a function on the one hand of the observation zone and on the other hand a function of the spatial gain of the apparatus.
55. The method according to claim 39 wherein the value subtracted from each image of the image sequence of the object is a function on the one hand of the observation zone and on the other hand a function of the spatial gain of the apparatus.
56. The method according to claim 32 wherein the calibration measurements of at least two acquisitions of successive sequences are combined.
57. The method according to claim 38 wherein the calibration measurements of at least two acquisitions of successive sequences are combined.
58. The method according to claim 45 wherein the calibration measurements of at least two acquisitions of successive sequences are combined.
59. The method according to claim 53 wherein the calibration measurements of at least two acquisitions of successive sequences are combined.
60. The method according to claim 32 wherein the calibration measurements of at least two acquisitions of successive sequences are combined.
61. The method according to claim 38 wherein the calibration measurements of at least two acquisitions of successive sequences are combined.
62. The method according to claim 45 wherein the calibration measurements of at least two acquisitions of successive sequences are combined.
63. The method according to claim 53 wherein the calibration measurements of at least two acquisitions of successive sequences are combined.
64. The method according to claim 56 wherein the calibration measurements of at least two acquisitions of successive sequences are combined.
65. An apparatus capable of acquiring a sequence of radiographic images of an object under observation comprising:
- means for each image of a sequence acquired by the apparatus and for a given acquisition frequency, calibrating the apparatus comprising means for determining the value of the variation of a mean of gray levels in at least one chosen zone of interest for the current image of at least one calibration device and determining the variation relative to the mean gray level of the first image of the sequence in each zone of interest;
- means for reiterating the determination of the variation for a series of image sequences acquired using calibration devices resulting in first images of mean gray levels different from one sequence to another;
- means for correcting each image of an image sequence of the object under observation comprising zones of observation having different gray levels and subtracting from the current image the variation of one gray level relative to the first image of the object, and effecting the subtraction as a function of the gray level considered of each zone of observation.
66. The apparatus according to claim 65 comprising:
- means for approximating a graphic representation having, respectively, as ordinate and abscissa, wherein the ratio having in the numerator the mean gray level variation of the current image of the sequence of the calibration device relative to the mean gray level of the first image and in the denominator the mean gray level of the first image, and the different mean gray levels of each first image; and
- means for using a function whose characteristics can be determined by an operator of the apparatus.
67. A computer program comprising code means that when executed on a computer or a means for processing carry out all the steps of claim 1.
68. A computer program on a carrier carrying code that when executed on a computer or a means for processing carry out all the steps of claim 1.
69. A computer apparatus or means for processing for carrying out all the steps of claim 1.
70. A method of operating a means for data processing comprising:
- for each image of a sequence acquired by an apparatus and for a given frequency of acquisition of the sequence, the apparatus is calibrated by determining the value of the variation of a mean of gray levels in at least one zone of interest of the current image of at least one calibration device, the variation being determined relative to the mean gray level of the first image of the sequence in each zone of interest;
- the determination of the variation is reiterated for a series of images sequences acquired using calibration devices resulting in first images of mean gray levels different from one sequence to another; and
- each image of an image sequence of the object under observation is corrected, comprising zones of observation having different gray levels by subtracting from the current image the variation of one gray level relative to the first image of the object, the subtraction being a function of the gray level considered from each zone of observation.
Type: Application
Filed: Feb 10, 2004
Publication Date: Feb 10, 2005
Inventors: Fanny Jeunehomme (Versailles), Serge Muller (Guyancourt), Razvan Iordache (Paris)
Application Number: 10/775,912