Patent US5164993A successfully refreshed from Google Patents.
← Back to Patent Explorer
US5164993A Expired / Lapsed US ⏱️ 0.0 years remaining

Method and apparatus for automatic tonescale generation in digital radiographic images

Assignee: Carestream Health Inc
•
Filing: 1991-11-25
•
Granted: 1992-11-17
•
Expires: 2011-11-25
📄 Original PDF 🌐 Google Patents 🔄 Sync Now

Abstract

Abstract A method and apparatus for automatically and adaptively generating tonescale transformation functions that are robust with respect to imaging systems, exposure conditions, and body parts. The technique uses the histogram of the digital input image, the cumulative distribution function of the histogram, and the entropy of subsections of the histogram to create the final tonescale transformation. Using these three functions, the histogram can be divided into a region of interest, a low-signal foreground region, and a high-signal background region. The tonescale is constructed to be substantially linear over the region of interest, joining smoothly with a nonlinear portion extending from the end of the low-signal foreground region to the start of the region of interest, and another nonlinear portion of the high-signal background region. The substantially linear region of interest and the entire image are subject to certain output density constraints to optimize the diagnostic utility of the final image. The technique depends entirely on the histogram of the input image, and, hence, it is adaptive and robust.

Inventors

  • 👤 John R. Capozzi
  • 👤 Ralph Schaetzing

Bibliographic Details

Application Number —
Jurisdiction US
Data Source google_patents
Last Synchronized 2026-10-24 06:08
CPC Classifications
G; G06; G06T; G06T5/00; G06T5/90; G06T5/94; G06T5/40; G06T7/00; G06T7/10; G06T7/194; G06T2207/00; G06T2207/10; G06T2207/10116; Y; Y10; Y10S; Y10S378/00; Y10S378/901

Claims Analysis

Total Claims: 23
Claim 1
Independent Claim 310 words
A method for generating a tonescale transformation for digital radiographic images comprising the steps of: providing a digital radiographic image; in a digital computer: generating both linear and logarithmic histograms from said digital radiographic image; forming a linear cumulative distribution function (CDF) from said linear histogram; forming a logarithmic CDF from said logarithmic histogram; determining start and end points for each of said linear and logarithmic histograms, respectively, where said start and end points are the first and last gray levels whose histogram values are non zero; determining for each of said linear and logarithmic histograms, respective left points which delineate the border between foreground and region of interest, right points which delineate the border between the region of interest and remaining anatomy and background points which delineate the border between remaining anatomy and true direct x-ray background; constructing a tonescale function, (a) having a linear portion from the left point to the right point in the region of interest, computing the average density in the region of interest, and depending upon a classification of the radiographic image, raising or lowering the median value until it is in a predetermined optimum range for the digital radiographic image classification; (b) having one or more line segments from the right point to the background point, wherein the end point coordinates of the line segment(s) are determined by the local structure of the linear CDF; (c) having a roll off between the start point and left point using a moving average window; and (d) raising or lowering the average density of the entire tonescale depending upon a predetermined average density for a classification of the radiographic image; and either processing the digital radiographic image with the tonescale function to produce a processed radiographic image, or creating a look-up table containing the tonescale function and transmitting the look-up table with the digital radiographic image.
Claim 2
Independent Claim 310 words
A method for generating a tonescale transformation for digital radiographic images comprising the steps of: providing a digital radiographic image; in a digital computer: generating both linear and logarithmic histograms from said digital radiographic image; forming a linear cumulative distribution function (CDF) from said linear histogram; forming a logarithmic CDF from said logarithmic histogram; determining start and end points for each of said linear and logarithmic histograms, respectively, where said start and end points are the first and last gray levels whose histogram values are non zero; determining for each of said linear and logarithmic histograms, respective left points which delineate the border between foreground and region of interest, right points which delineate the border between the region of interest and remaining anatomy and background points which delineate the border between remaining anatomy and true direct x-ray background; constructing a tonescale function, (a) having a linear portion from the left point to the right point in the region of interest, computing the average density in the region of interest, and depending upon a classification of the radiographic image, raising or lowering the median value until it is in a predetermined optimum range for the digital radiographic image classification; (b) having one or more line segments from the right point to the background point, wherein the end point coordinates of the line segment(s) are determined by the local structure of the linear CDF; (c) having a roll off between the start point and left point using a moving average window; and (d) raising or lowering the average density of the entire tonescale depending upon a predetermined average density for a classification of the radiographic image; and either processing the digital radiographic image with the tonescale function to produce a processed radiographic image, or creating a look-up table containing the tonescale function and transmitting the look-up table with the digital radiographic image.
Claim 3
Dependent Claim 97 words
The method of claim 1 wherein said step of determining the left point for each of said linear and logarithmic histograms is calculated by using the CDF of each of the histograms or each gray level from 0 to a predetermined gray level, calculating the ratio of the slopes of the gray levels of a predefined number of gray levels above the current one to a predefined number of gray levels below the current one and defining the left point as the gray level having the maximum slope change in the region below said predetermined gray level.
Claim 4
Dependent Claim 195 words
The method of claim 2 wherein said slope ratio is determined as follows: (a) in the gray level region from 0-n, a window size of n is used to calculate slope changes so that for each gray level, two CDF slopes are computed such that the average CDF slope in the n gray levels above the current gray level is computed, the average CDF slope for the n gray levels below the current gray level is computed and the ratio of the slopes of the gray levels above the current one to the gray levels below the current one is stored in a memory; (b) in the gray level range n-2n, the same calculation is done except that a varying window size is used starting at n at gray level n and increasing to 2n at gray level 2n; (c) between gray level 2n and background point-2n, the same calculation is done with a fixed window of 2n; and (d) the slope change values stored in memory are scanned to find the maximum slope change in the region where the CDF is less than 0.5 and this gray level is defined as the left point.
Claim 5
Dependent Claim 64 words
The method of claim 2 wherein roll-off of the tonescale curve from the start point to the left point is determined by looking at the magnitude of the maximum slope change used to find the left point and as the left point gets closer to 0, decreasing the amount of roll-off, and as the magnitude of the maximum slope change decreases, increasing the roll-off.
Claim 6
Dependent Claim 58 words
The method of claim 1 wherein the right points of said linear and logarithmic histograms are determined using a partial entropy technique wherein an effective width of the region of interest is determined by analyzing the structure of the respective histogram and the right point is determined by the sum of the left point and the effective width.
Claim 7
Dependent Claim 135 words
The method of claim 1 wherein said right point of said linear and logarithmic histograms is determined using the CDF of each histogram such that a CDF slope change ratio is found for each gray level starting at the gray level where the respective CDF has reached 0.5 and using a window of 2n gray levels such that for each gray level, the average CDF slope of the following 2n gray levels is found and the average slope of the 2n preceding gray levels is found and the ratio of the following slope to the preceding slope is stored in memory, continuing to find the ratio for each gray level within two window widths of the initial background point and determining the gray level where the maximum slope change ratio occurs as the right point.
Claim 8
Dependent Claim 54 words
The method of claim 6 wherein said digital image is a 12 bit image and wherein said window is set at 200 so that the ratio for each gray level is determined by taking the average CDF slope of the following 200 gray levels over the average slope of the 200 preceding gray levels.
Claim 9
Dependent Claim 134 words
The method of claim 6, including the further step of finding the right points of said linear and said logarithmic histogram by means of a partial entropy technique which produces an effective width of the region of interest by analyzing the structure of the histogram and the respective right point is determined by the sum of the left point and the effective width, and further including the step of mapping the right point of the logarithmic histogram determined by the CDF slope ratio technique and the right point of the logarithmic histogram determined by the partial entropy technique back into linear space to produce four estimates of the right point and then selecting an optimum right point for tone scale construction using a set of closeness criteria in comparing the various right point estimates.
Claim 10
Dependent Claim 96 words
The method of claim 5 wherein said provided digital radiograph image is a image of an anatomical extremity of a human and wherein the correct right point is chosen by determining the percentage of pixels added to the respective CDF per gray level in the input, and if this number is greater than or equal to 0.0075%, the higher of the two right point estimates is chosen as the correct right point and if the percentage of pixels per gray level is less than 0.0075%, then the lower of the two right point estimates is used.
Claim 11
Dependent Claim 93 words
The method of claim 1 wherein said background point is an initial estimate and a final refinement of the background point is made using the maximum slope change technique such that starting at the right point for each gray level, the average CDF slope of the following x gray levels and the average slope of the proceeding x gray levels are computed and their ratio found and stored in a memory and choosing the gray level in which the slope change is maximum as the final background point for the tone scale function.
Claim 12
Dependent Claim 170 words
The method of claim 1 wherein the step of constructing a tonescale function in the linear portion from the left point to the right point is effected as a graph of output density as a function of gray level and the initial tonescale is a straight line going from left point 0 to right point median where median is initially set to a value of 1.7, and including average density contributed by pixels in the gray level range from left point to right point as the sum of the values of the tonescale at each gray level multiplied by the percentage of pixels at that gray level summed from left to right point, wherein said predetermined optimum range for the digital radiographic image classification is a density value for a digital radiographic image of an anatomical part of a human and said median value of said linear portion of said tonescale function is raised or lowered until it is in a predetermined optical density range for a specified anatomical part.
Claim 13
Dependent Claim 15 words
The method of claim 11 wherein the optical density range for a chest is D=0.77-0.83.
Claim 14
Dependent Claim 15 words
The method of claim 11 wherein the optimum density range for an abdomen is D=0.83-0.90.
Claim 15
Dependent Claim 15 words
The method of claim 11 wherein the optimum density range for an extremity is D=1.03-1.09.
Claim 16
Dependent Claim 115 words
The method of claim 1 wherein said step of constructing a tonescale function between said right point and said background point, includes determining whether the slope of a line from the right point to the background point is larger than the slope from the left point to the right point, and if it is, then the tonescale curve is classified as a roll-under curve and the tonescale curve from the right point to the background point is constructed so that the slope is the same as that in the region of interest, that is the line segment representing the region of interest is extended to the right until it hits the background point gray level.
Claim 17
Dependent Claim 66 words
The method of claim 1 wherein the step of constructing a tonescale function between said right point and said background point includes determining whether the slope from the right point to the background point is less than that of the region of interest and, if it is, the tonescale is classified as a roll-over curve and multiple line segments are used to connect the background point.
Claim 18
Dependent Claim 130 words
The method of claim 16 wherein said line segments are determined as follows. If the slope in the region of interest is higher than the slope between the right point and background point by more than a maximum gain decrement, the background point is moved closer to the right point until this change is equal to the maximum gain decrement and the region between right point and background point is filled in with a set of line segments, each starting where the previous one ended with slopes less than or equal to the segment before with the actual length and slope of each segment depending on the percentage increase in the CDF for fixed percentage increases in the remaining gray levels between the current gray level and the background point.
Claim 19
Dependent Claim 45 words
The method of claim 17 wherein the fixed percentage increase in gray levels used to control the segment slopes is set as 20% for a digital radiographic image of a human chest and at 25% for a digital radiographic image of other human anatomical parts.
Claim 20
Dependent Claim 106 words
The method of claim 18 wherein each newly constructed line segment is temporarily extrapolated towards a maximum gray level to see if it hits the maximum density before or after the background point and if it hits maximum density before reaching the background point, another line segment is added to the tonescale, and if it does not reach the maximum density before the background point, the current line segment is altered to reach the maximum density at the background point and the tonescale construction in this region is terminated so that the tonescale is constant at maximum density from the background point to the end point.
Claim 21
Dependent Claim 77 words
The method of claim 1 wherein the step of constructing a tonescale function to raise or lower the average density of the entire tonescale is dependent upon a predetermined average density of a digital radiographic image of a human anatomical part such that the predetermined density range for a chest radiographic image is D=0.87-0.90, the predetermined density range for an abdomen radiographic image is D=0.93-0.96, and the predetermined density range for an extremity radiographic image is D=1.12-1.15.
Claim 22
Dependent Claim 88 words
The method of claim 1 wherein said method of constructing a tonescale is smoothed to eliminate any slope discontinuities in the region of the tonescale containing line segments, by using a window size approximately one-tenth the number of available gray levels, starting at a half a window size below the right point and proceeding to half a window size away from the maximum gray level and between this point and the maximum gray level, a continuously decreasing window size is used until the maximum gray level is reached.
Claim 23
Dependent Claim 30 words
The method of claim 1 wherein if the input calibration tonescale function is logarithmic, the linear tonescale function is transformed into log space using a logarithmic to linear conversion function.

References Cited (Backward Citations)

US4302672A Cited
US4641268A Cited
US4853967A Cited
US4823394A Cited
US4778745A Cited
US4845558A Cited
US5046012A Cited
US5046118A Cited

USPTO Annuity & Maintenance Windows

Annuity Stage Window Opens Payment Deadline Surcharge Period Until Estimated Fee Status
3.5 Years (Year 4) 1995-11-17 1996-05-18 1996-11-17 $2000.00 Window Closed / Lapsed
7.5 Years (Year 8) 1999-11-17 2000-05-18 2000-11-17 $3760.00 Window Closed / Lapsed
11.5 Years (Year 12) 2003-11-17 2004-05-18 2004-11-17 $7700.00 Window Closed / Lapsed

Assigned Categories & Tags

No categories assigned to this patent yet.

Assign Category Tag