CA2385631A1 - A method for generating gray scale transfer functions for use in displaying a digital radiogram - Google Patents
A method for generating gray scale transfer functions for use in displaying a digital radiogram Download PDFInfo
- Publication number
- CA2385631A1 CA2385631A1 CA002385631A CA2385631A CA2385631A1 CA 2385631 A1 CA2385631 A1 CA 2385631A1 CA 002385631 A CA002385631 A CA 002385631A CA 2385631 A CA2385631 A CA 2385631A CA 2385631 A1 CA2385631 A1 CA 2385631A1
- Authority
- CA
- Canada
- Prior art keywords
- curve
- values
- representing
- displaying
- contrast
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000012546 transfer Methods 0.000 title claims abstract description 11
- 230000006870 function Effects 0.000 title description 18
- 238000004590 computer program Methods 0.000 claims 2
- 230000000007 visual effect Effects 0.000 description 11
- 230000005855 radiation Effects 0.000 description 8
- 238000004422 calculation algorithm Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 238000002601 radiography Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 208000010392 Bone Fractures Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013481 data capture Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G06T5/92—
Abstract
A method for generating on demand a gray scale transfer curve for use in displaying an image comprised of plurality of digital values representing pixel gray scale levels with a desired contrast and brightness (Fig. 3). This is achieved by using a single sigmoid shaped function Y=f(X,a,b,c) representing an H & D curve of a photographic film. The Y axis represents digital output values to be used in displaying an image (Fig.2), while the X
axis represents input digital values representing captured image information (Fig. 2). "a" is a number representing boundary conditions for X and Y, "b" is a first parameter controlling the location of the curve along the X axis and "c" is a second parameter controlling the slope of the curve slope. This function is used to derive a Look-Up-Table (LUT) representing a continuous GST
curve having the desired contrast and brightness.
axis represents input digital values representing captured image information (Fig. 2). "a" is a number representing boundary conditions for X and Y, "b" is a first parameter controlling the location of the curve along the X axis and "c" is a second parameter controlling the slope of the curve slope. This function is used to derive a Look-Up-Table (LUT) representing a continuous GST
curve having the desired contrast and brightness.
Description
W~ 01/39104 CA 02385631 2002-03-22 PCT/US00/30938 TITLE
A METHOD FOR GENERATING GRAY SCALE TRANSFER FUNCTIONS FOR
USE IN DISPLAYING A DIGITAL RADIOGRAM
FIELD OF THE INVENTION
This invention relates to a method for displaying a radiographic image composed of a plurality of digital values, and more particularly to a method for generating a look up table (LUT) representing a desired gray scale transfer (GST) function for displaying a visually enhanced radiogram.
BACKGROUND OF THE INVENTION
There exists significant activity in the development of digital x-ray image data capture systems. In such systems direct conversion to an electrical signal of the incident radiation is obtained using a plurality of sensors (also known as pixels) in an array. The sensor output is almost invariably immediately converted to a digital signal by an analog-to-digital converter as known in the art and further processed and stored in a databank for use in the eventual display of the data as a radiograph.
United States patent Number 5,313,066 issued to Lee et al. (hereinafter the '063 patent) and United States patent 5,315,101 issued to Hughes et al. describe typical such sensor arrays and their contents are incorporated herein by reference. Even though several different technologies are being utilized, the output data are quite similar.
A major advantage of digital data detection systems is the wide dynamic range of signal capture. Display media, such as radiographic film or cathode ray tube (CRT) displays, on the other hand, have a substantially more limited dynamic range. A
typical digital x-ray capture system can have a useful dynamic detection range of greater than a 1,000:1. However, the typical currently available display media are generally limited to a dynamic range of less than 100:1. There is, therefore, need to WO 01/39104 CA 02385631 2002-03-22 pCT~jS00/30938 determine and select the optimal limited range of useful data for diagnostic display, and then properly display such range on the available display medium.
This problem, which reduces to a need for a method whereby the exposure sensor output is mapped onto the display transfer function of the display device, has been addressed by the art in numerous ways. Typically the sensor output is digitized, and a histogram of the frequency of occurrence of digital values representing detected exposure is constructed. Histogram analysis is often used to determine the relevant portions of the data, that is data carrying significant diagnostic information, and the digital values in this range are mapped onto the display transfer function, usually using a look-up table (LUT), as is well known in the art. These steps are rather fundamental and are well known to the person skilled in this art. What is significant and the subject of continuing research is the manner in which the digital values are processed and mapped to display the maximum amount of relevant information within the limitations of the available display media. This in turn requires the generation of the transfer function to create an optimal radiograph wherein the features of interest are distinguishable from background features and wherein maximum visual contrast is applied to the region of interest.
The method most often used in actual practice for mapping the data derived from the detector, is to store a plurality of GST curves representing empirically derived GST shapes in the form of LUTs in a memory , and to select one of the stored GSTs to reproduce a given image. If the desired region of interest in the reproduction proves too dark, or has insufficient contrast, another GST curve is chosen and the output observed again and so on, until an acceptable image is produced. Such method, however, is limited by the number of curves stored, and does not allow for continuous changing of the contrast and brightness of the displayed image.
United States patent number 4,641, 267 issued Feb. 3 1987 to Asai et al. shows an early attempt to overcome some of the above limitations. Asai begins by generating first WO 01/39104 CA 02385631 2002-03-22 pCT/[JS00/30938 a set of reference tone correction curves. In displaying an image, Asai selects one of the reference curves and a point Y( x0) on the curve corresponding to a desired point on the Y (output value) axis. By rotating the curve around this point Asai can change the contrast of the output image.
United States patent number 5,946,407 issued August 31, 1999 to Bamberger et al. is an improvement on the method taught by Asai and is believed to represent the current state in developing GSTs. Bamberger teaches creating the GST as a combination of two curves, by using two algorithms to describe two distinct portions of a GST. The point where the two curves connect is set at the 50 % of the Y
(output value) axis. The slope of the two curves around the connecting point determines the contrast of the displayed image in the area of most interest, while moving the connecting point along the X (input value) changes the brightness (or film density if the display medium is film) of the areas of interest. The term brightness as used hereinafter is to be considered as referring to the gray scale gradations of an image in general, and to include optical density when the display medium is hard copy.
While this method provides good results, it is computationally intensive.
Furthermore, in certain applications it is desirable to calculate the point of maximum slope of the GST (maximum contrast) which is usually done by obtaining the second derivative of the curve. The discontinuity due to the use of two algorithms to form the GST introduces complications in this calculation. It is desirable, therefore to have a method for quickly developing accurate GSTs on demand, given certain input parameters representing desirable contrast and brightness characteristics with optimum visual appearance similar to the visual appearance obtained as a result of the typical H&D response of photographic materials used in capturing radiograms using the traditional method of exposing the photographic material to X-ray radiation.
WO 01/39104 CA 02385631 2002-03-22 pCT~JS00/30938 SUMMARY OF THE INVENTION
It is an object of this invention to provide a method for generating on demand a gray scale transfer curve for use in displaying an image comprised of a plurality of digital values representing pixel gray scale levels with a desired contrast and brightness. This is achieved by using a single sigmoid shaped function Y=f(X, a, b, c) wherein the function resembles an H&D curve of a photographic film. The Y axis represents digital output values to be used in displaying an image, while the X axis represents input digital values representing captured image information.
"a" is a number representing boundary conditions for X and Y, " b" is a first parameter controlling the location of the curve along the X axis and "c" is a second parameter controlling the slope of the curve slope. This function is used to derive a Look-Up-Table (LUT) representing a continuous GST curve having the desired contrast and brightness.
The above described method permits the derivation of families of GST curves and corresponding LUTs representing different contrasts and brightness, by varying the values of the parameters "b" and "c" .
Particularly useful is the function Y = a/(1+b*e~~x). Parameter "c" is dependent on the maximum input digital value possible. Experience has shown that in a system where the input values have been normalized and range between 0 and 100, good results are obtained with values for "c" selected between 0.02% and 0.2% of the maximum scale value for X, and preferably between 0.06% and 0.10%. Parameter "b" is given by:
Xmax contrast 0 = ~ b and the value of "a" may be calculated by solving the above given function for "a"
using:
a = Y° * ( 1+ b*e~xu ).
WO 01/39104 cA 02385631 2002-03-22 pCT/[JS00/30938 The input digital values for the X and Y axis may be normalized preferably between 0-100, in which case "a" is derived for Y~=100 and X~=100.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 shows a typical H&D curve for a screen-film combination used in traditional radiography.
Figure 2 shows a family of GST curves obtained through the process of this invention illustrating the change of contrast as a function of parameter "c".
Figure 3 shows a family of GST curves obtained through the process of this invention illustrating the change in output values and thus brightness levels (or film density if the display medium is film) for a given region of input values as a function of parameter "b".
DETAIL DESCRIPTION OF THE INVENTION
The invention will next be described with reference to the figures, wherein similar number indicate similar parts in all figures.
Figure 1 shows two typical H&D curves for two hypothetical photographic films used in capturing X-ray radiation. What these curves show are the exposed and developed film density as a result of radiation exposure. The Y axis represents the film density and the X axis the logarithm of the exposure. Exposure is the product of the incident radiation intensity times the duration of incidence of the radiation onto the film. The exposure may be the result of irradiating the film with x-rays only or it may be the result of exposing the film through a film intensifying screen combination.
Different films have different H&D curves. The slope of the curve at any point indicates the contrast characteristic of the particular film or film screen combination at a WO 01/39104 CA 02385631 2002-03-22 pCT/jJS00/30938 given exposure level. The shift of the curve to the right or left relative to the Y axis is related to the reproduced density, (or brightness) of portions of the exposed areas, and is, in photographic film, a function of the film screen system sensitivity.
The medical profession has been trained over the years to evaluate the results of S radiographic examinations in the form of radiograms, that is films exposed to x-ray radiation, that display areas of different radiation absorption in accordance to the film's H&D curve. Photographic film happens to have a response wherein density as a function of the logarithm of exposure closely matches the eye's response and thus produces desirable visual appearance images.
Electronic image capture is almost always a linear process of extremely wide range. As stated in the introduction above, in traditional radiography the exposure latitude is of the order of 100 to 1 or less, while in electronic radiography the exposure latitude is of the order of 1,000 to 1 or more, and has introduced a set of display problems unique to this field.
In an effort to display an electronic radiogram in such manner as to simulate the visual appearance of a traditional radiogram while controlling contrast and brightness it is essential to process the raw data captured by the electronic device and following digital conversion. Such processing typically includes: (a) determining the exposure range and image densities that carry significant information for diagnostic purposes, (b) compressing the extensive data obtained from the detector usually from 14 bits to 12 or 8 bits and (c) mapping the data through an input - output transformation so that the output data used in displaying the image will display an image of desirable visual characteristics with optimum contrast at optimum brightness, or when the display medium is hard copy such as film, density.
The determination of the data values representing a significant brightness or density range, and the compression of the data values are separate subjects not directly related to understanding the present invention. For the present invention purposes the imaging data values have already undergone such processing and have been reduced to the input data values that will be displayed, typically 12 bits representing image values from 0 to 4095 ( 4096 digital values), or 8 bits representing image values varying from 0 to 255 (256 digital values). United states patent 5,633,511 and U.S. patent 5,946,407 are two reference among others, that disclose methods for performing steps (a) and (b) .
Once the data has reached the stage where it is available in format for display as a plurality of digital values, it is next mapped onto a gray scale transfer (GST) curve through the use of a look up table (LUT), and the new values are sent to a display device for image display. The digital values prior to processing through the LUT are referred to hereinafter as the input digital values while the values resulting from the LUT operation on the input values are referred to as the output digital values. The digital output values are used to drive a display device which can be a CRT
monitor in which case the output values represent brightness levels of the displayed image on the monitor, or the display device may be a hard copy output device with the output being a photosensitive material, in which case the output values represent optical density (gray scale) of the output image. Typically a higher output value represents a darker image in CRT displays and a higher density image (darker) image in hard copy display.
The present invention proposes the use of a single algorithm to simulate the generally "S" shaped H&D curve of the typical film screen combination, and use this algorithm to construct GST curves having different contrast and brightness characteristics. While there are a number of possible algorithms that produce a generally "S" shaped curve, it has been found that a function Y= f (a, b, c, X) where "a" has a value dependent on selected boundary conditions, and where "b" and "c" are parameters that control the slope and relative position along the X-axis of the curve represented by the function can be used to quickly derive any one of an almost infinite number of "S" type curves representing different displayed image contrast and brightness, or density. The ability to vary the curves is provided ;~y the two parameters WO 01/39104 CA 02385631 2002-03-22 PCT/jJS00/30938 "b" and "c", which may be varied within predetermined limits to provide families of "S" type curves that generate images exhibiting different contrast characteristics and/or brightness (or density).
More particularly according to a preferred embodiment of this invention the following algorithm satisfies the above requirements:
(1) Y = a/(1+b*e~'X)~
The slope of this curve, represents the contrast of an image produced using the output values obtained from input values. Selecting the point of maximum contrast along the X axis defines "b" and is given by:
A METHOD FOR GENERATING GRAY SCALE TRANSFER FUNCTIONS FOR
USE IN DISPLAYING A DIGITAL RADIOGRAM
FIELD OF THE INVENTION
This invention relates to a method for displaying a radiographic image composed of a plurality of digital values, and more particularly to a method for generating a look up table (LUT) representing a desired gray scale transfer (GST) function for displaying a visually enhanced radiogram.
BACKGROUND OF THE INVENTION
There exists significant activity in the development of digital x-ray image data capture systems. In such systems direct conversion to an electrical signal of the incident radiation is obtained using a plurality of sensors (also known as pixels) in an array. The sensor output is almost invariably immediately converted to a digital signal by an analog-to-digital converter as known in the art and further processed and stored in a databank for use in the eventual display of the data as a radiograph.
United States patent Number 5,313,066 issued to Lee et al. (hereinafter the '063 patent) and United States patent 5,315,101 issued to Hughes et al. describe typical such sensor arrays and their contents are incorporated herein by reference. Even though several different technologies are being utilized, the output data are quite similar.
A major advantage of digital data detection systems is the wide dynamic range of signal capture. Display media, such as radiographic film or cathode ray tube (CRT) displays, on the other hand, have a substantially more limited dynamic range. A
typical digital x-ray capture system can have a useful dynamic detection range of greater than a 1,000:1. However, the typical currently available display media are generally limited to a dynamic range of less than 100:1. There is, therefore, need to WO 01/39104 CA 02385631 2002-03-22 pCT~jS00/30938 determine and select the optimal limited range of useful data for diagnostic display, and then properly display such range on the available display medium.
This problem, which reduces to a need for a method whereby the exposure sensor output is mapped onto the display transfer function of the display device, has been addressed by the art in numerous ways. Typically the sensor output is digitized, and a histogram of the frequency of occurrence of digital values representing detected exposure is constructed. Histogram analysis is often used to determine the relevant portions of the data, that is data carrying significant diagnostic information, and the digital values in this range are mapped onto the display transfer function, usually using a look-up table (LUT), as is well known in the art. These steps are rather fundamental and are well known to the person skilled in this art. What is significant and the subject of continuing research is the manner in which the digital values are processed and mapped to display the maximum amount of relevant information within the limitations of the available display media. This in turn requires the generation of the transfer function to create an optimal radiograph wherein the features of interest are distinguishable from background features and wherein maximum visual contrast is applied to the region of interest.
The method most often used in actual practice for mapping the data derived from the detector, is to store a plurality of GST curves representing empirically derived GST shapes in the form of LUTs in a memory , and to select one of the stored GSTs to reproduce a given image. If the desired region of interest in the reproduction proves too dark, or has insufficient contrast, another GST curve is chosen and the output observed again and so on, until an acceptable image is produced. Such method, however, is limited by the number of curves stored, and does not allow for continuous changing of the contrast and brightness of the displayed image.
United States patent number 4,641, 267 issued Feb. 3 1987 to Asai et al. shows an early attempt to overcome some of the above limitations. Asai begins by generating first WO 01/39104 CA 02385631 2002-03-22 pCT/[JS00/30938 a set of reference tone correction curves. In displaying an image, Asai selects one of the reference curves and a point Y( x0) on the curve corresponding to a desired point on the Y (output value) axis. By rotating the curve around this point Asai can change the contrast of the output image.
United States patent number 5,946,407 issued August 31, 1999 to Bamberger et al. is an improvement on the method taught by Asai and is believed to represent the current state in developing GSTs. Bamberger teaches creating the GST as a combination of two curves, by using two algorithms to describe two distinct portions of a GST. The point where the two curves connect is set at the 50 % of the Y
(output value) axis. The slope of the two curves around the connecting point determines the contrast of the displayed image in the area of most interest, while moving the connecting point along the X (input value) changes the brightness (or film density if the display medium is film) of the areas of interest. The term brightness as used hereinafter is to be considered as referring to the gray scale gradations of an image in general, and to include optical density when the display medium is hard copy.
While this method provides good results, it is computationally intensive.
Furthermore, in certain applications it is desirable to calculate the point of maximum slope of the GST (maximum contrast) which is usually done by obtaining the second derivative of the curve. The discontinuity due to the use of two algorithms to form the GST introduces complications in this calculation. It is desirable, therefore to have a method for quickly developing accurate GSTs on demand, given certain input parameters representing desirable contrast and brightness characteristics with optimum visual appearance similar to the visual appearance obtained as a result of the typical H&D response of photographic materials used in capturing radiograms using the traditional method of exposing the photographic material to X-ray radiation.
WO 01/39104 CA 02385631 2002-03-22 pCT~JS00/30938 SUMMARY OF THE INVENTION
It is an object of this invention to provide a method for generating on demand a gray scale transfer curve for use in displaying an image comprised of a plurality of digital values representing pixel gray scale levels with a desired contrast and brightness. This is achieved by using a single sigmoid shaped function Y=f(X, a, b, c) wherein the function resembles an H&D curve of a photographic film. The Y axis represents digital output values to be used in displaying an image, while the X axis represents input digital values representing captured image information.
"a" is a number representing boundary conditions for X and Y, " b" is a first parameter controlling the location of the curve along the X axis and "c" is a second parameter controlling the slope of the curve slope. This function is used to derive a Look-Up-Table (LUT) representing a continuous GST curve having the desired contrast and brightness.
The above described method permits the derivation of families of GST curves and corresponding LUTs representing different contrasts and brightness, by varying the values of the parameters "b" and "c" .
Particularly useful is the function Y = a/(1+b*e~~x). Parameter "c" is dependent on the maximum input digital value possible. Experience has shown that in a system where the input values have been normalized and range between 0 and 100, good results are obtained with values for "c" selected between 0.02% and 0.2% of the maximum scale value for X, and preferably between 0.06% and 0.10%. Parameter "b" is given by:
Xmax contrast 0 = ~ b and the value of "a" may be calculated by solving the above given function for "a"
using:
a = Y° * ( 1+ b*e~xu ).
WO 01/39104 cA 02385631 2002-03-22 pCT/[JS00/30938 The input digital values for the X and Y axis may be normalized preferably between 0-100, in which case "a" is derived for Y~=100 and X~=100.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 shows a typical H&D curve for a screen-film combination used in traditional radiography.
Figure 2 shows a family of GST curves obtained through the process of this invention illustrating the change of contrast as a function of parameter "c".
Figure 3 shows a family of GST curves obtained through the process of this invention illustrating the change in output values and thus brightness levels (or film density if the display medium is film) for a given region of input values as a function of parameter "b".
DETAIL DESCRIPTION OF THE INVENTION
The invention will next be described with reference to the figures, wherein similar number indicate similar parts in all figures.
Figure 1 shows two typical H&D curves for two hypothetical photographic films used in capturing X-ray radiation. What these curves show are the exposed and developed film density as a result of radiation exposure. The Y axis represents the film density and the X axis the logarithm of the exposure. Exposure is the product of the incident radiation intensity times the duration of incidence of the radiation onto the film. The exposure may be the result of irradiating the film with x-rays only or it may be the result of exposing the film through a film intensifying screen combination.
Different films have different H&D curves. The slope of the curve at any point indicates the contrast characteristic of the particular film or film screen combination at a WO 01/39104 CA 02385631 2002-03-22 pCT/jJS00/30938 given exposure level. The shift of the curve to the right or left relative to the Y axis is related to the reproduced density, (or brightness) of portions of the exposed areas, and is, in photographic film, a function of the film screen system sensitivity.
The medical profession has been trained over the years to evaluate the results of S radiographic examinations in the form of radiograms, that is films exposed to x-ray radiation, that display areas of different radiation absorption in accordance to the film's H&D curve. Photographic film happens to have a response wherein density as a function of the logarithm of exposure closely matches the eye's response and thus produces desirable visual appearance images.
Electronic image capture is almost always a linear process of extremely wide range. As stated in the introduction above, in traditional radiography the exposure latitude is of the order of 100 to 1 or less, while in electronic radiography the exposure latitude is of the order of 1,000 to 1 or more, and has introduced a set of display problems unique to this field.
In an effort to display an electronic radiogram in such manner as to simulate the visual appearance of a traditional radiogram while controlling contrast and brightness it is essential to process the raw data captured by the electronic device and following digital conversion. Such processing typically includes: (a) determining the exposure range and image densities that carry significant information for diagnostic purposes, (b) compressing the extensive data obtained from the detector usually from 14 bits to 12 or 8 bits and (c) mapping the data through an input - output transformation so that the output data used in displaying the image will display an image of desirable visual characteristics with optimum contrast at optimum brightness, or when the display medium is hard copy such as film, density.
The determination of the data values representing a significant brightness or density range, and the compression of the data values are separate subjects not directly related to understanding the present invention. For the present invention purposes the imaging data values have already undergone such processing and have been reduced to the input data values that will be displayed, typically 12 bits representing image values from 0 to 4095 ( 4096 digital values), or 8 bits representing image values varying from 0 to 255 (256 digital values). United states patent 5,633,511 and U.S. patent 5,946,407 are two reference among others, that disclose methods for performing steps (a) and (b) .
Once the data has reached the stage where it is available in format for display as a plurality of digital values, it is next mapped onto a gray scale transfer (GST) curve through the use of a look up table (LUT), and the new values are sent to a display device for image display. The digital values prior to processing through the LUT are referred to hereinafter as the input digital values while the values resulting from the LUT operation on the input values are referred to as the output digital values. The digital output values are used to drive a display device which can be a CRT
monitor in which case the output values represent brightness levels of the displayed image on the monitor, or the display device may be a hard copy output device with the output being a photosensitive material, in which case the output values represent optical density (gray scale) of the output image. Typically a higher output value represents a darker image in CRT displays and a higher density image (darker) image in hard copy display.
The present invention proposes the use of a single algorithm to simulate the generally "S" shaped H&D curve of the typical film screen combination, and use this algorithm to construct GST curves having different contrast and brightness characteristics. While there are a number of possible algorithms that produce a generally "S" shaped curve, it has been found that a function Y= f (a, b, c, X) where "a" has a value dependent on selected boundary conditions, and where "b" and "c" are parameters that control the slope and relative position along the X-axis of the curve represented by the function can be used to quickly derive any one of an almost infinite number of "S" type curves representing different displayed image contrast and brightness, or density. The ability to vary the curves is provided ;~y the two parameters WO 01/39104 CA 02385631 2002-03-22 PCT/jJS00/30938 "b" and "c", which may be varied within predetermined limits to provide families of "S" type curves that generate images exhibiting different contrast characteristics and/or brightness (or density).
More particularly according to a preferred embodiment of this invention the following algorithm satisfies the above requirements:
(1) Y = a/(1+b*e~'X)~
The slope of this curve, represents the contrast of an image produced using the output values obtained from input values. Selecting the point of maximum contrast along the X axis defines "b" and is given by:
(2) Xmaxconcrast= In (b) / c Two parameters, (c) and (b), control the slope and location of maximum slope along the X-axis of the resulting GST curve generated using equation ( 1 ).
Once parameter "c" is selected, parameter "b" may be calculated using equation (2).
Parameter "c" is related to the scaling of the X-axis, and the maximum possible digital value, i.e. maximum scale value, of the input digital values on the X-axis, Xmax value "c" has been found by experience to be preferably between 0.01 % and 0.2% of Xmax value . When the input digital values are normalized, that is the X axis value range is between 0 and 100, preferred values for "c" are between 0.01 and 0.20 with the most preferred range being between 0.06 and 0.1. Where the X-axis values are from 0-255, "c" values range, i.e. 0.05 to 0.5, and so on. Values outside these ranges may be also be used, however experience has shown that most applications will obtain optimum visual results using values within the given range.
Depending on where the maximum contrast point is placed on the X axis, a different value for parameter "b" is obtained. Once parameters "c" and "b"
have been obtained, "a" is derived by solving equation (1) for desired boundary conditions. From ( 1 ) we obtain (3) a = I'~ ( 1 + b * e-~x~ ) The boundary conditions depend on the digital input and output value range. Xu and Y~ are the X and Y values representing the upper limits of possible input and output digital values. In an 8 bit system the input and output boundary conditions for X" and Y
are 255 and 255. Preferably both the input and output values are normalized digital values ranging from 0-100, in which case X~ and Y~ are equal to 100. In such case, a=100(1+b*evoo~) Figure 2 shows a family of GTS curves obtained using equation (1) with parameter "c" varied between c = 0.06, 0.08, and 0.10. In all cases XmaX
contrast = 45. As is seen from the figure the slope of the curve changes resulting in more contrast in the displayed image, increasing as "c" increases. Figure 3 on the other hand shows another family of curves where the parameter "c" controlling contrast is held constant while parameter "b" has been changed. In this case the point of maximum contrast has been shifted progressively to the right along the X-axis, Xmas contrast being 45, 60 and 75 respectively. Equation (2) yields the following values for "b": 28.79, 121.51 and 403.
The result of this curve shift is a change in the output value for a given input value. For example a digital value of 60 will result in an output value of 28, 54 and 82 respectively. A value of 28 will be represented as much brighter than one of 82 in the final image. Depending on the particular portion of the data one wants to emphasize any one of these curves may provide better visual results. For instance if the information sought to be observed is in the range of between 60 and 80 input digital normalized values, then curve III is probably the best suited to display the information since this data range is displayed over a range of 30 to 80 output normalized values. On the other hand if the desired input range is between 30 and 60 in the input axis, curve I is best.
Once parameter "c" is selected, parameter "b" may be calculated using equation (2).
Parameter "c" is related to the scaling of the X-axis, and the maximum possible digital value, i.e. maximum scale value, of the input digital values on the X-axis, Xmax value "c" has been found by experience to be preferably between 0.01 % and 0.2% of Xmax value . When the input digital values are normalized, that is the X axis value range is between 0 and 100, preferred values for "c" are between 0.01 and 0.20 with the most preferred range being between 0.06 and 0.1. Where the X-axis values are from 0-255, "c" values range, i.e. 0.05 to 0.5, and so on. Values outside these ranges may be also be used, however experience has shown that most applications will obtain optimum visual results using values within the given range.
Depending on where the maximum contrast point is placed on the X axis, a different value for parameter "b" is obtained. Once parameters "c" and "b"
have been obtained, "a" is derived by solving equation (1) for desired boundary conditions. From ( 1 ) we obtain (3) a = I'~ ( 1 + b * e-~x~ ) The boundary conditions depend on the digital input and output value range. Xu and Y~ are the X and Y values representing the upper limits of possible input and output digital values. In an 8 bit system the input and output boundary conditions for X" and Y
are 255 and 255. Preferably both the input and output values are normalized digital values ranging from 0-100, in which case X~ and Y~ are equal to 100. In such case, a=100(1+b*evoo~) Figure 2 shows a family of GTS curves obtained using equation (1) with parameter "c" varied between c = 0.06, 0.08, and 0.10. In all cases XmaX
contrast = 45. As is seen from the figure the slope of the curve changes resulting in more contrast in the displayed image, increasing as "c" increases. Figure 3 on the other hand shows another family of curves where the parameter "c" controlling contrast is held constant while parameter "b" has been changed. In this case the point of maximum contrast has been shifted progressively to the right along the X-axis, Xmas contrast being 45, 60 and 75 respectively. Equation (2) yields the following values for "b": 28.79, 121.51 and 403.
The result of this curve shift is a change in the output value for a given input value. For example a digital value of 60 will result in an output value of 28, 54 and 82 respectively. A value of 28 will be represented as much brighter than one of 82 in the final image. Depending on the particular portion of the data one wants to emphasize any one of these curves may provide better visual results. For instance if the information sought to be observed is in the range of between 60 and 80 input digital normalized values, then curve III is probably the best suited to display the information since this data range is displayed over a range of 30 to 80 output normalized values. On the other hand if the desired input range is between 30 and 60 in the input axis, curve I is best.
4 CA 02385631 2002-03-22 PCT/iJS00/30938 The equations are typically used as follows. Following capture of a radiogram and conversion to a digital set of values, the raw data is processed to compensate for individual pixel gain lack of uniformity and to replace dead or bad pixels.
This is well known technology and needs no further discussion. See for example currently pending patent application serial number 09/255,772 filed Feb 23, 1999 by Chuande Liu and assigned to the assignee of the present invention, as well as United States patent Number 4,602,291 issued to Temes on Jul. 22, 1986 and United States patent Number 5,164,993 issued to Capozzi et al. on Nov.l7, 1992 as examples of the technology used for this purpose.
Again typically the data after this processing stage is compressed from 14 bits to 12 bits and may undergo further processing as described in co-pending application filed contemporaneously herewith in the name of Gregory Powell and assigned to the assignee of the present application, or as described in the aforementioned United States patent Number 5,633,511 issued to Lee et al.
Following this and any other processing the image data is now ready for display.
At this time, before the data is forwarded to the display device, it is mapped onto a GST
curve. This curve is generated on demand, based on default values for the parameters "c", "b" and "a" . The point of maximum contrast for use as Xmax contrast, usually varies between 30 and 80 and is preferably about 4~.
Default parameters have been developed based on experimental results for different types of radiograms, i.e. chest lateral, chest frontal, GI, bone fractures etc.
Typical values for parameters were given above in this description. The default parameters provide a beginning point from which one may, after observing the displayed image, select parameters and modify the GST to produce a visual appearance best suited to the particular examination needs. For example, if the image appears generally flat, parameter "c" may be increased to provide better contrast. Similarly if a particular region of special interest in the image appears dark, Xmax concr~t ~d p~~eter "b" can be altered to shift the GST to the left or to the right changing the brightness of that region.
WO 01/39104 CA 02385631 2002-03-22 pCT/[JS00/30938 The GST may be applied to the input values on the fly, so to speak, by calculating output (Y) values for each input (X) value of the digital image, or a LUT
table may be created prior to displaying the image and stored in memory for later use during the image display. In either case there is never a need to store a large number of LUTs representing different GSTs in the system, and to attempt to optimize the visual appearance of the displayed radiogram by selecting one of the stored GSTs.
This invention provides an infinite number of GSTs any one of which can be rapidly generated and tailored to produce an image with proper contrast, brightness (or density), and generally a visual appearance that simulates accurately the appearance of similar radiograms obtained through traditional radiography using a film screen combination.
The invention is best implemented through a computer programmed to derive the GST curve and/or a corresponding LUT representing the function Y = f (X, a, b, c). Such program may be embodied in a computer readable medium such as magnetic tape or disc, or in a CD or any other medium usable for programming a computer. Programming a computer to solve for values of Y given a function such as equation 1 is well known in the art. Such program will be structured to accept input values for "c" arid Xmas ~o~n~t from a keyboard or other input device, or may retrieve such values from a stored data base containing preselected values for "c" and Xmax co~ttast for different types of radiograms. The program will next calculate "b" and "a" by solving equations (2) and (3). Values for Y~ and X~ may be input either by the operator or, preferably be stored in memory in the computer as part of an initialization process.
Once the two equations have been solved and the values of "a", "b" and "c"
obtained, the program will use these parameter values to solve equation ( 1 ) for each input value. This may be done by providing the program with the ability to receive input 2~ values and calculate an output value each time an input value is received, or to pre-calculate all output values for the full range of input, X axis, values and form a LUT
which is stored in a memory and used to supply output values for each input value received.
Such program may be a stand alone program, or, preferably, may be part of a larger program used to process the data received from a radiation detector and to supply such data for display to a display medium.
Those skilled in the art having the benefit of the above description can make numerous modifications particularly as to the actual numerical values used in the examples given above. Any such modifications are to be construed as encompassed within the scope of the invention as claimed herein below.
This is well known technology and needs no further discussion. See for example currently pending patent application serial number 09/255,772 filed Feb 23, 1999 by Chuande Liu and assigned to the assignee of the present invention, as well as United States patent Number 4,602,291 issued to Temes on Jul. 22, 1986 and United States patent Number 5,164,993 issued to Capozzi et al. on Nov.l7, 1992 as examples of the technology used for this purpose.
Again typically the data after this processing stage is compressed from 14 bits to 12 bits and may undergo further processing as described in co-pending application filed contemporaneously herewith in the name of Gregory Powell and assigned to the assignee of the present application, or as described in the aforementioned United States patent Number 5,633,511 issued to Lee et al.
Following this and any other processing the image data is now ready for display.
At this time, before the data is forwarded to the display device, it is mapped onto a GST
curve. This curve is generated on demand, based on default values for the parameters "c", "b" and "a" . The point of maximum contrast for use as Xmax contrast, usually varies between 30 and 80 and is preferably about 4~.
Default parameters have been developed based on experimental results for different types of radiograms, i.e. chest lateral, chest frontal, GI, bone fractures etc.
Typical values for parameters were given above in this description. The default parameters provide a beginning point from which one may, after observing the displayed image, select parameters and modify the GST to produce a visual appearance best suited to the particular examination needs. For example, if the image appears generally flat, parameter "c" may be increased to provide better contrast. Similarly if a particular region of special interest in the image appears dark, Xmax concr~t ~d p~~eter "b" can be altered to shift the GST to the left or to the right changing the brightness of that region.
WO 01/39104 CA 02385631 2002-03-22 pCT/[JS00/30938 The GST may be applied to the input values on the fly, so to speak, by calculating output (Y) values for each input (X) value of the digital image, or a LUT
table may be created prior to displaying the image and stored in memory for later use during the image display. In either case there is never a need to store a large number of LUTs representing different GSTs in the system, and to attempt to optimize the visual appearance of the displayed radiogram by selecting one of the stored GSTs.
This invention provides an infinite number of GSTs any one of which can be rapidly generated and tailored to produce an image with proper contrast, brightness (or density), and generally a visual appearance that simulates accurately the appearance of similar radiograms obtained through traditional radiography using a film screen combination.
The invention is best implemented through a computer programmed to derive the GST curve and/or a corresponding LUT representing the function Y = f (X, a, b, c). Such program may be embodied in a computer readable medium such as magnetic tape or disc, or in a CD or any other medium usable for programming a computer. Programming a computer to solve for values of Y given a function such as equation 1 is well known in the art. Such program will be structured to accept input values for "c" arid Xmas ~o~n~t from a keyboard or other input device, or may retrieve such values from a stored data base containing preselected values for "c" and Xmax co~ttast for different types of radiograms. The program will next calculate "b" and "a" by solving equations (2) and (3). Values for Y~ and X~ may be input either by the operator or, preferably be stored in memory in the computer as part of an initialization process.
Once the two equations have been solved and the values of "a", "b" and "c"
obtained, the program will use these parameter values to solve equation ( 1 ) for each input value. This may be done by providing the program with the ability to receive input 2~ values and calculate an output value each time an input value is received, or to pre-calculate all output values for the full range of input, X axis, values and form a LUT
which is stored in a memory and used to supply output values for each input value received.
Such program may be a stand alone program, or, preferably, may be part of a larger program used to process the data received from a radiation detector and to supply such data for display to a display medium.
Those skilled in the art having the benefit of the above description can make numerous modifications particularly as to the actual numerical values used in the examples given above. Any such modifications are to be construed as encompassed within the scope of the invention as claimed herein below.
Claims (8)
1. A method for generating on demand a gray scale transfer (GST) curve for use in deriving a Look-Up-Table (LUT) for displaying an image comprised of a plurality of digital values representing pixel densities with a desired contrast and brightness, the method comprising using a single sigmoid function Y=f(X, a, b, c) wherein Y
represents output values used in displaying an image, X represents input digital values representing captured image information, "b" is a first parameter controlling the location of the curve along the X axis, "c" is a second parameter controlling the curve slope and "a" is determined from boundary conditions.
represents output values used in displaying an image, X represents input digital values representing captured image information, "b" is a first parameter controlling the location of the curve along the X axis, "c" is a second parameter controlling the curve slope and "a" is determined from boundary conditions.
2. The method according to claim 1 wherein the parameters "b" and "c" are varied to produce a family of LUTs representing a family of GST curves having desired different contrast and brightness.
3. The method according to claim 1 wherein the function F=f(X, a, b, c) is:
Y= a/(1+b*e-cX) where b = inv In (X max contrast*c);"c" is between .01% to 0.2% of X max value and:
a = Y u*(1+b*e-cXu).
Y= a/(1+b*e-cX) where b = inv In (X max contrast*c);"c" is between .01% to 0.2% of X max value and:
a = Y u*(1+b*e-cXu).
4. The method according to claim 3 wherein the input digital values for the X
axis are normalized input digital values.
axis are normalized input digital values.
5. The method according to claim 4 wherein the output digital values from the Y
axis are also normalized output digital values.
axis are also normalized output digital values.
6. A program embodied in a machine readable medium for programming a control computer to generate on demand a gray scale transfer (GST) curve for use in deriving a Look-Up-Table (LUT) for displaying an image comprised of a plurality of digital values representing pixel densities with a desired contrast and brightness, using a single sigmoid function Y=f(X, a, b, c) wherein Y represents output values used in displaying an image, X represents input digital values representing captured image information, "b"
is a first parameter controlling the location of the curve along the X axis, "c" is a second parameter controlling the curve slope and "a" is a number representing boundary conditions.
is a first parameter controlling the location of the curve along the X axis, "c" is a second parameter controlling the curve slope and "a" is a number representing boundary conditions.
7. The computer program according to claim 6 wherein the function F=f(X, a, b, c) is:
Y=a/(1+b*e-cX) where b = inv In (X max contrast*c) ; "c" is between .01% t0 0.2% of X max value and:
a=Y u*(1+b*e-cXu).
Y=a/(1+b*e-cX) where b = inv In (X max contrast*c) ; "c" is between .01% t0 0.2% of X max value and:
a=Y u*(1+b*e-cXu).
8. The computer program according to claim 7 wherein the output values are digital output values representing a radiogram.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US09/449,075 US6370265B1 (en) | 1999-11-24 | 1999-11-24 | Method for generating gray scale transfer functions for use in displaying a digital radiogram |
| US09/449,075 | 1999-11-24 | ||
| PCT/US2000/030938 WO2001039104A1 (en) | 1999-11-24 | 2000-11-10 | A method for generating gray scale transfer functions for use in displaying a digital radiogram |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2385631A1 true CA2385631A1 (en) | 2001-05-31 |
Family
ID=23782777
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002385631A Abandoned CA2385631A1 (en) | 1999-11-24 | 2000-11-10 | A method for generating gray scale transfer functions for use in displaying a digital radiogram |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US6370265B1 (en) |
| EP (1) | EP1236161A4 (en) |
| JP (1) | JP2003515825A (en) |
| CA (1) | CA2385631A1 (en) |
| WO (1) | WO2001039104A1 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117078676A (en) * | 2023-10-16 | 2023-11-17 | 青岛盛鑫达机械有限公司 | Breaking hammer shell visual detection method based on image processing |
Families Citing this family (26)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE3721074A1 (en) * | 1986-12-24 | 1988-07-07 | Licentia Gmbh | CIRCUIT ARRANGEMENT FOR SHIFTING VECTOR SEARCH IN DIGITAL IMAGE ANALYSIS |
| JP2002366952A (en) * | 2001-04-05 | 2002-12-20 | Konica Corp | Method, device and program for displaying medical image and recording medium |
| US7046834B2 (en) * | 2001-05-09 | 2006-05-16 | Electronics And Telecommunications Research Institute | Method for measuring bone mineral density by using X-ray image |
| US6956977B2 (en) * | 2001-08-08 | 2005-10-18 | General Electric Company | Methods for improving contrast based dynamic range management |
| US7102789B2 (en) * | 2001-09-17 | 2006-09-05 | Sharp Laboratories Of America, Inc. | Method for rendering an image comprising multi-level pixels |
| US7113649B2 (en) * | 2002-06-24 | 2006-09-26 | Eastman Kodak Company | Enhancing the tonal characteristics of digital images |
| US7158686B2 (en) * | 2002-09-19 | 2007-01-02 | Eastman Kodak Company | Enhancing the tonal characteristics of digital images using inflection points in a tone scale function |
| US7043090B2 (en) * | 2002-10-02 | 2006-05-09 | Eastman Kodak Company | Enhancing the tonal characteristics of digital images using expansive and compressive tone scale functions |
| JP4323770B2 (en) * | 2002-10-18 | 2009-09-02 | キヤノン株式会社 | Image processing apparatus, image processing method, program, and recording medium |
| US7116838B2 (en) * | 2002-10-25 | 2006-10-03 | Eastman Kodak Company | Enhancing the tonal and spatial characteristics of digital images using selective spatial filters |
| JP2004212598A (en) * | 2002-12-27 | 2004-07-29 | Sharp Corp | Converting device, correcting circuit, driving device, display device, inspecting device, and display method |
| JP2004336516A (en) * | 2003-05-09 | 2004-11-25 | Konica Minolta Photo Imaging Inc | Image processing method, image processor, and image recording apparatus |
| FR2858087B1 (en) * | 2003-07-25 | 2006-01-21 | Eastman Kodak Co | METHOD FOR DIGITIGLY SIMULATING AN IMAGE SUPPORT RENDER |
| JP2005348237A (en) * | 2004-06-04 | 2005-12-15 | Brother Ind Ltd | Image density adjusting apparatus and image reader with the same |
| US20060034502A1 (en) * | 2004-08-10 | 2006-02-16 | Konica Minolta Medical & Graphic, Inc. | Image processing system and image processing method |
| JP4372747B2 (en) * | 2005-01-25 | 2009-11-25 | シャープ株式会社 | Luminance level conversion device, luminance level conversion method, solid-state imaging device, luminance level conversion program, and recording medium |
| US8004538B2 (en) * | 2007-09-21 | 2011-08-23 | Global Oled Technology Llc | Preferential tone scale for electronic displays |
| US20090263015A1 (en) * | 2008-04-17 | 2009-10-22 | Guoyi Fu | Method And Apparatus For Correcting Underexposed Digital Images |
| CN102484736A (en) * | 2009-07-07 | 2012-05-30 | 汤姆森特许公司 | Method and system for brightness correction for three-dimensional (3D) projection |
| JP5868119B2 (en) * | 2010-12-09 | 2016-02-24 | キヤノン株式会社 | Image processing apparatus, radiation imaging system, image processing method, and recording medium |
| CN102298772A (en) * | 2011-09-07 | 2011-12-28 | 谭洪舟 | Method for eliminating pole of liquid crystal display (LCD) motion image inverse model |
| US8988552B2 (en) | 2011-09-26 | 2015-03-24 | Dolby Laboratories Licensing Corporation | Image formats and related methods and apparatuses |
| KR101768435B1 (en) | 2011-12-06 | 2017-08-17 | 돌비 레버러토리즈 라이쎈싱 코오포레이션 | Device and method of improving the perceptual luminance nonlinearity-based image data exchange across different display capabilities |
| US10242650B2 (en) | 2011-12-06 | 2019-03-26 | Dolby Laboratories Licensing Corporation | Perceptual luminance nonlinearity-based image data exchange across different display capabilities |
| US9911036B2 (en) * | 2012-11-14 | 2018-03-06 | Tascent, Inc. | Focusing method for optically capturing an iris image |
| US11193898B1 (en) * | 2020-06-01 | 2021-12-07 | American Science And Engineering, Inc. | Systems and methods for controlling image contrast in an X-ray system |
Family Cites Families (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5983149A (en) | 1982-11-04 | 1984-05-14 | Fuji Photo Film Co Ltd | Forming method of curve for correcting gradation |
| US4602291A (en) | 1984-05-09 | 1986-07-22 | Xerox Corporation | Pixel non-uniformity correction system |
| JP2703010B2 (en) * | 1988-12-23 | 1998-01-26 | 株式会社日立製作所 | Neural net signal processing processor |
| JP2906281B2 (en) * | 1990-09-05 | 1999-06-14 | セイコーインスツルメンツ株式会社 | Optical pattern recognition device |
| US5220398A (en) * | 1990-09-28 | 1993-06-15 | Massachusetts Institute Of Technology | Analog VLSI microchip for object position and orientation |
| US5164993A (en) | 1991-11-25 | 1992-11-17 | Eastman Kodak Company | Method and apparatus for automatic tonescale generation in digital radiographic images |
| GB9202693D0 (en) | 1992-02-08 | 1992-03-25 | Philips Electronics Uk Ltd | A method of manufacturing a large area active matrix array |
| US5313066A (en) | 1992-05-20 | 1994-05-17 | E. I. Du Pont De Nemours And Company | Electronic method and apparatus for acquiring an X-ray image |
| IL119767A (en) | 1993-08-13 | 1998-02-08 | Sophis View Tech Ltd | System and method for diagnosis of living tissue diseases |
| JPH07248841A (en) * | 1994-03-09 | 1995-09-26 | Mitsubishi Electric Corp | Nonlinear function generator and format converter |
| JPH0836644A (en) * | 1994-07-26 | 1996-02-06 | Matsushita Electric Works Ltd | Image recognition method by neural net |
| WO1996009717A1 (en) * | 1994-09-19 | 1996-03-28 | Apple Computer, Inc. | Generation of tone reproduction curves using psychophysical data |
| JP3557673B2 (en) * | 1994-12-08 | 2004-08-25 | 富士ゼロックス株式会社 | Image processing device |
| US5662113A (en) * | 1995-06-30 | 1997-09-02 | Siemens Medical Systems, Inc | Edge enhancement system for ultrasound images |
| US5588435A (en) * | 1995-11-22 | 1996-12-31 | Siemens Medical Systems, Inc. | System and method for automatic measurement of body structures |
| US5633511A (en) | 1995-12-22 | 1997-05-27 | Eastman Kodak Company | Automatic tone scale adjustment using image activity measures |
| KR100251967B1 (en) * | 1998-02-28 | 2000-04-15 | 윤종용 | Scan format converter |
-
1999
- 1999-11-24 US US09/449,075 patent/US6370265B1/en not_active Expired - Fee Related
-
2000
- 2000-11-10 JP JP2001540696A patent/JP2003515825A/en active Pending
- 2000-11-10 WO PCT/US2000/030938 patent/WO2001039104A1/en not_active Application Discontinuation
- 2000-11-10 CA CA002385631A patent/CA2385631A1/en not_active Abandoned
- 2000-11-10 EP EP00978500A patent/EP1236161A4/en not_active Withdrawn
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117078676A (en) * | 2023-10-16 | 2023-11-17 | 青岛盛鑫达机械有限公司 | Breaking hammer shell visual detection method based on image processing |
| CN117078676B (en) * | 2023-10-16 | 2024-01-12 | 青岛盛鑫达机械有限公司 | Breaking hammer shell visual detection method based on image processing |
Also Published As
| Publication number | Publication date |
|---|---|
| EP1236161A4 (en) | 2005-05-11 |
| JP2003515825A (en) | 2003-05-07 |
| EP1236161A1 (en) | 2002-09-04 |
| US6370265B1 (en) | 2002-04-09 |
| WO2001039104A1 (en) | 2001-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6370265B1 (en) | Method for generating gray scale transfer functions for use in displaying a digital radiogram | |
| US4317179A (en) | Method and apparatus for processing a radiographic image | |
| JP2551552B2 (en) | A method of displaying or recording an image in which details of the original image are emphasized | |
| US5357549A (en) | Method of dynamic range compression of an X-ray image and apparatus effectuating the method | |
| US5493622A (en) | Radiation image processing method which increase and decreases a frequency region of the radiation image | |
| JP2000132674A (en) | Picture processing method | |
| JP2005522237A (en) | X-ray image enhancement | |
| JP3538244B2 (en) | How to find saturated pixels in displaying radiographic images. | |
| US6215900B1 (en) | Image data processing for digital x-ray detectors | |
| JP3201853B2 (en) | Signal-to-density mapping method with controlled contrast | |
| JP3599376B2 (en) | Method for monitoring the sensitivity of a device for reading radiation images recorded on a photostimulable phosphor screen | |
| US6118892A (en) | Method for automatic detection of region of interest for digital x-ray detectors using a filtered histogram | |
| JP3597272B2 (en) | Abnormal shadow candidate detection method | |
| US6381351B1 (en) | Weighted inverse topography method for digital x-ray image data processing | |
| JPH05200021A (en) | Reproduction method of medical image that provides optimum quality of image for diagnosis | |
| JPH01261787A (en) | Radiograph processor | |
| Angelhed et al. | Optimization and evaluation of an image intensifier TV system for digital chest imaging | |
| USH1627H (en) | Method of and apparatus for image processing using a variable focal spot size | |
| EP0527524B1 (en) | Image quality enhancing method | |
| JPH0572624B2 (en) | ||
| JP3197952B2 (en) | Method for generating an enhanced visible image of an image represented by an original image signal | |
| JPS6262378B2 (en) | ||
| JPH09261471A (en) | Digital image photographing device | |
| JPH04347142A (en) | Radiation image processor |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FZDE | Discontinued | ||
| FZDE | Discontinued |
Effective date: 20061110 |