CA2276634A1 - Spatially-selective edge enhancement for discrete pixel images - Google Patents
Spatially-selective edge enhancement for discrete pixel images Download PDFInfo
- Publication number
- CA2276634A1 CA2276634A1 CA002276634A CA2276634A CA2276634A1 CA 2276634 A1 CA2276634 A1 CA 2276634A1 CA 002276634 A CA002276634 A CA 002276634A CA 2276634 A CA2276634 A CA 2276634A CA 2276634 A1 CA2276634 A1 CA 2276634A1
- Authority
- CA
- Canada
- Prior art keywords
- image
- images
- frequency
- gain
- frequency level
- 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
- 230000035945 sensitivity Effects 0.000 claims abstract description 9
- 238000000034 method Methods 0.000 claims description 36
- 238000001914 filtration Methods 0.000 claims description 14
- 230000005855 radiation Effects 0.000 claims description 11
- 238000003384 imaging method Methods 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 9
- 230000002708 enhancing effect Effects 0.000 claims description 4
- 238000004846 x-ray emission Methods 0.000 claims 1
- 238000000354 decomposition reaction Methods 0.000 abstract description 20
- 230000006870 function Effects 0.000 description 10
- 239000002131 composite material Substances 0.000 description 4
- 230000000750 progressive effect Effects 0.000 description 4
- 238000012886 linear function Methods 0.000 description 2
- WURBVZBTWMNKQT-UHFFFAOYSA-N 1-(4-chlorophenoxy)-3,3-dimethyl-1-(1,2,4-triazol-1-yl)butan-2-one Chemical compound C1=NC=NN1C(C(=O)C(C)(C)C)OC1=CC=C(Cl)C=C1 WURBVZBTWMNKQT-UHFFFAOYSA-N 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G06T5/73—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Abstract
A discrete pixel image is enhanced to bring out particular features of interest such as edges. The enhancement includes decomposition of the original processed image data, enhancement by application of gain images to the decomposed images, and reconstruction of the enhanced image. The decomposition proceeds through a series of low-pass filters to arrive at decomposed images of progressively lower spatial frequencies. These decomposed images are then multiplied by gain images during the enhancement phase. Gain images for at least the higher spatial frequency level images are derived from lower spatial frequency level images. These spatial frequency-based gain images may be based upon operator inputs, including a spatial sensitivity function and an edge enhancement value.
Lower spatial frequency level image may be processed by application of predetermined gain values. The reconstruction sequence recombines the decomposed images after application of the gains to arrive at an enhanced image of the same dimensions as the original image.
Lower spatial frequency level image may be processed by application of predetermined gain values. The reconstruction sequence recombines the decomposed images after application of the gains to arrive at an enhanced image of the same dimensions as the original image.
Description
Spatially-Selective Edge Enhancement For Discrete Pixel Images Technical Field The present invention relates'. to digital imaging techniques and, more particularly, to techniques for enhancement of discrete pixel images, such as those produced in medical im~iging systems.
Background Art A wide range of systems have keen proposed and are currently in use for producing discrete pixel images. In such systems, data is gathered representative of characteristics of a large number of picture elements or pixels arranged in an array or matrix. For example, in digital radiography systems, signals are produced which are representative of a level or intensity of radiation received within each pixel region during an exposure. The signals are processed and filtered to provide consistent and meaningfizl information over the matrix of pixels. lFollowing such processing, the data is used to reconstruct a useful composite; image made up of the pixels. The particular filtering and processing of the signals may be adapted to various modalities and features of interest; so as to provide a user with meaningfial information in their reconstructed image. In many medical applications, an attending physician or radiologist will consult the composite image for identification of internal features within a subject as defined by edges, textural regions, contrasted regions, and. so forth.
Algorithms have been developed for digitally processing discrete pixel image data to enhance diagnostic portions of the image while suppressing noise. For example, in one; known method pixel data is filtered through progressive low pass filtering steps. The original image data is thus decomposed into a sequence of images having known frequency bands.
Gain values are applied to the resulting decomposed images for enhancement of image features, such as edges. Additional filtering, contrast equalization, and gradation steps may be employed for further enhancement of the image.
While such techniques provide; useful mechanisms for certain types of image enhancement, they are not mithout drawbacks. For example, gains . applied to decomposed images can result in inadvertent enhancement of noise present in the discrete pixel data. Such noise, when enhanced, renders the reconstructed image difficult to interpret, and may produce visual artifacts which reduce the utility of the reconstructed image, such as by rendering features of interest difficult t:o discern or to distinguish from non relevant information.
There is a need, therefore, for an improved technique for discrete pixel image enhancement. In particular, there is a need for a technique which is capable of enhancing features of interest in such images, such as edges of structures, without enhancement of noise or other artifacts. There is also a particular need 'for an image enhancement technique which can be employed in existing systems, and which can be adjusted by users for various situations, depending upon the; type of subject being imaged, the features of interest, and so forth.
Disclosure of the Invention The present invention provides a technique for enhancing digital pixel images designed to respond to these needs. The technique affords enhancement of particular features in the image as defined by the digital pixel data, while avoiding enhancement of noise and other artifacts. The technique applies feature enhancement gain factors which are based upon the level of an enhanced image in a decomposition pyramid. Gains employed for different levels in the decomposed images can be adjusted, with lower gains being employed for decomposed images at lower spatial frequency levels. The technique may be applied to images acquired through a range of imaging modalities. The technique is particularly well suited to enhancement of images acquired through digital X-ray systems.
In accordance with certain .aspects of the technique, processed acquired pixel data for an image is filtered to decompose the image into a series of spatial frequency level images. The series of images includes images at progressively lower spatial frequency levels. A gain image for at least one of the spatial frequency imap;es is derived from image data from a lower spatial frequency level. This gain is applied to the higher spatial frequency image to obtain an enhanced image for that spatial frequency level. Such gains may be derived for and applied to a series of the decomposed spatial frequency level innages. Such spatial frequency-based gain images may be employed for spa~,tial frequency level images at which valuable information, such as edge ir,~formation, is expected to be found.
1 S Other lower spatial frequency level images may be processed by applying gains equal to unity or gains higher th~,an or less than unity. The resulting images are then recombined to arrive at an enhanced image. The decomposition, enhancement, and reconstruction phases of the technique preferably produce an enhanced image of equal dimensions to the input image.
Brief Description of the Drawings Fig. 1 is a diagrammatical representation of a digital X-ray imaging system for producing discrete pixel images of a subject of interest;
Fig. 2 is plan view of an exemplary discrete pixel image produced via the system of Fig. l;
Fig. 3 is a diagrammatic~il representation of an image decomposition sequence for separating the image of Fig. 2 into decomposed images referenced by desired frequency bands; and Fig. 4 is a block diagram illustrating steps in exemplary control logic for enhancing the image of Fig. 2 based upon the decomposed or recomposed images analyzed as illusb~ated in Fig. 3.
Best Mode for Carrying Out the Invention Turning now to the drawings, and referring first to Figure 1, an imaging system 10 is illustrated diagrammatically in the form of a digital x-ray system. System 10 includes a controlled source of radiation 12 capable of emitting a stream 14 of radiation upon command. The stream of radiation is directed toward a subject of interest 16. A portion 18 of the radiation passes through the subject and impacts a detector 20. The detector, which includes an array of discrete picture elements or pixels generates signals representative of the amount or intensity of the radiation received within the bounds of each pixel region, and transmits the signals to data acquisition circuitry 22. Circuitry 22 processes the acquired signals, such as by appropriate filtering, dynamic range adjustment, and so forth, outputting a string of data representative of the radiation received by detector 20.
The digital signals output by a data acquisition circuitry 22 are transmitted to an image processing unut 24 where the signals are further processed to produce data representative. of features of interest within subject 16. Thus, the signals from data acquisition circuitry 22 are transmitted first to interface circuitry 26 within unit 24. This interface circuitry may perform additional filtering functions, before forwarding the signals to signal processing circuitry 28. As is described in greater detail below, signal processing circuitry 28 performs a series of computations based upon the acquired and filtered signals. In partiacular, signal processing circuit 28 filters and enhances the pixel data so as to render the image defined by the collection of pixels more readily understandable or useful to an attending physician or radiologist.
Background Art A wide range of systems have keen proposed and are currently in use for producing discrete pixel images. In such systems, data is gathered representative of characteristics of a large number of picture elements or pixels arranged in an array or matrix. For example, in digital radiography systems, signals are produced which are representative of a level or intensity of radiation received within each pixel region during an exposure. The signals are processed and filtered to provide consistent and meaningfizl information over the matrix of pixels. lFollowing such processing, the data is used to reconstruct a useful composite; image made up of the pixels. The particular filtering and processing of the signals may be adapted to various modalities and features of interest; so as to provide a user with meaningfial information in their reconstructed image. In many medical applications, an attending physician or radiologist will consult the composite image for identification of internal features within a subject as defined by edges, textural regions, contrasted regions, and. so forth.
Algorithms have been developed for digitally processing discrete pixel image data to enhance diagnostic portions of the image while suppressing noise. For example, in one; known method pixel data is filtered through progressive low pass filtering steps. The original image data is thus decomposed into a sequence of images having known frequency bands.
Gain values are applied to the resulting decomposed images for enhancement of image features, such as edges. Additional filtering, contrast equalization, and gradation steps may be employed for further enhancement of the image.
While such techniques provide; useful mechanisms for certain types of image enhancement, they are not mithout drawbacks. For example, gains . applied to decomposed images can result in inadvertent enhancement of noise present in the discrete pixel data. Such noise, when enhanced, renders the reconstructed image difficult to interpret, and may produce visual artifacts which reduce the utility of the reconstructed image, such as by rendering features of interest difficult t:o discern or to distinguish from non relevant information.
There is a need, therefore, for an improved technique for discrete pixel image enhancement. In particular, there is a need for a technique which is capable of enhancing features of interest in such images, such as edges of structures, without enhancement of noise or other artifacts. There is also a particular need 'for an image enhancement technique which can be employed in existing systems, and which can be adjusted by users for various situations, depending upon the; type of subject being imaged, the features of interest, and so forth.
Disclosure of the Invention The present invention provides a technique for enhancing digital pixel images designed to respond to these needs. The technique affords enhancement of particular features in the image as defined by the digital pixel data, while avoiding enhancement of noise and other artifacts. The technique applies feature enhancement gain factors which are based upon the level of an enhanced image in a decomposition pyramid. Gains employed for different levels in the decomposed images can be adjusted, with lower gains being employed for decomposed images at lower spatial frequency levels. The technique may be applied to images acquired through a range of imaging modalities. The technique is particularly well suited to enhancement of images acquired through digital X-ray systems.
In accordance with certain .aspects of the technique, processed acquired pixel data for an image is filtered to decompose the image into a series of spatial frequency level images. The series of images includes images at progressively lower spatial frequency levels. A gain image for at least one of the spatial frequency imap;es is derived from image data from a lower spatial frequency level. This gain is applied to the higher spatial frequency image to obtain an enhanced image for that spatial frequency level. Such gains may be derived for and applied to a series of the decomposed spatial frequency level innages. Such spatial frequency-based gain images may be employed for spa~,tial frequency level images at which valuable information, such as edge ir,~formation, is expected to be found.
1 S Other lower spatial frequency level images may be processed by applying gains equal to unity or gains higher th~,an or less than unity. The resulting images are then recombined to arrive at an enhanced image. The decomposition, enhancement, and reconstruction phases of the technique preferably produce an enhanced image of equal dimensions to the input image.
Brief Description of the Drawings Fig. 1 is a diagrammatical representation of a digital X-ray imaging system for producing discrete pixel images of a subject of interest;
Fig. 2 is plan view of an exemplary discrete pixel image produced via the system of Fig. l;
Fig. 3 is a diagrammatic~il representation of an image decomposition sequence for separating the image of Fig. 2 into decomposed images referenced by desired frequency bands; and Fig. 4 is a block diagram illustrating steps in exemplary control logic for enhancing the image of Fig. 2 based upon the decomposed or recomposed images analyzed as illusb~ated in Fig. 3.
Best Mode for Carrying Out the Invention Turning now to the drawings, and referring first to Figure 1, an imaging system 10 is illustrated diagrammatically in the form of a digital x-ray system. System 10 includes a controlled source of radiation 12 capable of emitting a stream 14 of radiation upon command. The stream of radiation is directed toward a subject of interest 16. A portion 18 of the radiation passes through the subject and impacts a detector 20. The detector, which includes an array of discrete picture elements or pixels generates signals representative of the amount or intensity of the radiation received within the bounds of each pixel region, and transmits the signals to data acquisition circuitry 22. Circuitry 22 processes the acquired signals, such as by appropriate filtering, dynamic range adjustment, and so forth, outputting a string of data representative of the radiation received by detector 20.
The digital signals output by a data acquisition circuitry 22 are transmitted to an image processing unut 24 where the signals are further processed to produce data representative. of features of interest within subject 16. Thus, the signals from data acquisition circuitry 22 are transmitted first to interface circuitry 26 within unit 24. This interface circuitry may perform additional filtering functions, before forwarding the signals to signal processing circuitry 28. As is described in greater detail below, signal processing circuitry 28 performs a series of computations based upon the acquired and filtered signals. In partiacular, signal processing circuit 28 filters and enhances the pixel data so as to render the image defined by the collection of pixels more readily understandable or useful to an attending physician or radiologist.
Signals processed by circuitry 28, typically in the form of a series of digital values, are stored in a memory circuit 30 within unit 24. Memory circuit 30, which may include any suitable type of electronic memory, such as random access memory, read-only memory, dynamic random access memory, disk storage, and so forth, also preferably stores programming code implemented by signal processing circuit 28 for executing its preconfigured functions. Signal processing circuitry 28 may, in turn, include any suitable general purpose or application-specific microprocessor, CPU, or signal processor.
In addition to the image enhancement functions carried out by signal processing circuitry 28, such circuitr;~ is also available for outputting or displaying composite images defined by acquired and processed pixel data on an output device as indicated by reference numeral 32. Such output or display devices may include conventional computer monitors, printers, and so forth. An input terminal 34 is also provided for communicating input parameters, image types, examination configures, and so forth to signal processing circuitry 28. Input termvial 34 may include a conventional keyboard, serial mouse, or any other suutable interface device. Terminal 34 may also include circuitry associated directly with a terminal screen for permitting selection or input of configuration parameters by touching appropriate locations in an on-screen display. Signal processing circuit 28 is further coupled to radiation source 12 for controlling the generation and transmission of radiation stream 14 as innage data is desired.
Figure 2 represents an exemplary image matrix 24 defined by data acquired through system 10. As illustrated in Figure 2, the matrix permits reconstruction of a composite image 38 in which particular features of interest 40 are defined by varying characteristics of adjacent pixels 42, typically gray level or color. Pixels 42 are arranged in a series of rows 44 and columns 46. The resulting image matrix thus includes a large number of such rows and columns defining an imal;e width 48 and an image height 50.
In addition to the image enhancement functions carried out by signal processing circuitry 28, such circuitr;~ is also available for outputting or displaying composite images defined by acquired and processed pixel data on an output device as indicated by reference numeral 32. Such output or display devices may include conventional computer monitors, printers, and so forth. An input terminal 34 is also provided for communicating input parameters, image types, examination configures, and so forth to signal processing circuitry 28. Input termvial 34 may include a conventional keyboard, serial mouse, or any other suutable interface device. Terminal 34 may also include circuitry associated directly with a terminal screen for permitting selection or input of configuration parameters by touching appropriate locations in an on-screen display. Signal processing circuit 28 is further coupled to radiation source 12 for controlling the generation and transmission of radiation stream 14 as innage data is desired.
Figure 2 represents an exemplary image matrix 24 defined by data acquired through system 10. As illustrated in Figure 2, the matrix permits reconstruction of a composite image 38 in which particular features of interest 40 are defined by varying characteristics of adjacent pixels 42, typically gray level or color. Pixels 42 are arranged in a series of rows 44 and columns 46. The resulting image matrix thus includes a large number of such rows and columns defining an imal;e width 48 and an image height 50.
The particular dimensions of matrix 3~i may depend upon such factors as the particular feature being imaged, the rE;solution desired in the image, and so forth. Moreover, system 10 may be~ configured to generate images of a number of different matrix dimension;., such as 2,048 X 2,048 pixels, 1,024 X 1,024 pixels, and so forth.
System 10 preferably filters the acquired signals to enhance certain types of features in the image. Such enhancement is performed in such a way as to selectively emphasize the features, while avoiding amplification or enhancement of noise or other artifact; in the image. Figure 3 illustrates a diagrarnmatical representation of exemplary fimctional signal processing circuitry employed for such enhancement. It should be noted that, while the elements illustrated in Figure 3 and described below are represented as discrete components for processing image data, in practice, such components may be defined by appropriate programming code executed by signal processing circuitry 28.
As shown in Figure 3, a filtering sequence 52 for spatially-selective edge enhancement of a discrete pixel i~r~age includes a decomposition phase, designated generally by the reference muneral 54, an enhancement phase 56, and a reconstruction phase 58. Decomposition phase 54 involves progressive decomposition of the image data into a series of spatial frequency images of progressively lower spatial frequency. As used herein, the tenor spatial frequency relates generally to the number of individual image components within a defined image area. Through the progressive decomposition of the image, filtering sequence 52 generates a series of decomposed images having different :spatial frequencies. Each of these images is then enhanced by application of a gain image during enhancement phase 56. The resulting images are then recombined during reconstruction phase 58 to complete an enhanced image. The individual processing of the different spatial frequency levels of the decomposed images is represented in Figure 3 by reference numerals 60, 62, fi4, 66 and 68, with 68 being a final level of decomposition and processing;. As noted, each progressive step or level produces an intermediate image having a progressively lower spatial frequency, as noted at reference numerals 70, 72, 74, 76, and 78.
As shown in Figure 3, processing blocks for each spatial frequency level in decomposition phase 54 include a decimator (LP~L) 80 and an interpolator (LPT) 82, represented by low-pass filters (only the filters of the highest spatial frequency level 60 being labeled in Figure 3). These filters progressively combine data representing neighboring pixels of the immediately preceding image level in .a manner generally known in the art.
The filters output a first, filtered imal;e which is employed in subsequent decomposition levels, as well as a second image which is subtracted from the input image for that spatial frequency level at a summer 84. Thus, at a highest level 60, the original processed image data 86 is fed to the decimating low pass filter 80, as well as to summer 84. The original image data is filtered at blocks 80 and 82, with the output of block 80 being passed on as input data to level 62. In this nranner, the original processed image data is decomposed into the progressively lower spatial frequency images.
In a present embodiment, each subsequent step or level in the decomposition phase produces an image which is half' the size of the input image for that level. However, other reduction rates rnay be employed. Moreover, while the decomposition phase filtering may terminate at any desired level, in a present embodiment, such filtering is carried out until a single pixel ( 1 X 1 ) image is attained at lowest level 68.
The filtering sequence 52 then proceeds by applying (e.g., multiplying) each of the spatial frequency images by a gain image in enhancement phase 56. The gain images, designated in Figure 3 by the numerals 88, 90, 92, 94, and 96, are assigned values which enhance specific features of the image, such as edges. Nforeover, certain of the gain images are generated from data from lower spatial frequency levels. Thus, gain image 88, designated G1, is generated from image 72 resulting from the decomposition phase at level 62. Similarly, gain image 90, designated G2, is generated from image 74 resulting from processing at level 64. Any suitable number of levels in the filtering sequence 52 may employ such lower frequency-based gain images. For e:Kample, it has been found that in an image having original dimensions of 2,048 X 2,048 pixels, such lower spatial frequency-based gain images produce satisfactory edge enhancement when used at the first three levels of decomposition, such as levels 60, 62, and 64 in Figure 3. Other lower levels may also employ this technique, or may be multiplied by a gain image having values of unity, or other values, either greater or less than unity.
Following application of the gain images at enhancement phase 56, the resulting enhanced images are recombined in reconstruction phase 58.
As shown in Figure 3, each step or level in the reconstruction phase includes interpolation of output data from an immediately preceding (i.e., lower) level in a block 98 (see level 60 in Figure 3), The resulting image is added to the enhanced image in a summer 100 at each level. The resulting enhanced and reconstructed image data is then forw~~rded to the next (i.e., higher) level, until a reconstructed enhanced image 102 of the original dimensions is obtained.
It should be noted that, in particular applications, the image data need not be fully decomposed or reconstructed to the lowest levels. For example, where a series of lower spatial frequency images are to have a similar gain applied in phase 56, decomposition may be stopped at that level. The gain is then applied to the image at the terminal level, and reconstruction may be begin. Similarly, while as noted above, the higher level gains are derived from lower spatial frequency images resulting from decomposition, these could also be obtained based upon reconstructed images from a lower spatial frequency level (i.e., following application of the gain images during enhancement phase 56).
System 10 preferably filters the acquired signals to enhance certain types of features in the image. Such enhancement is performed in such a way as to selectively emphasize the features, while avoiding amplification or enhancement of noise or other artifact; in the image. Figure 3 illustrates a diagrarnmatical representation of exemplary fimctional signal processing circuitry employed for such enhancement. It should be noted that, while the elements illustrated in Figure 3 and described below are represented as discrete components for processing image data, in practice, such components may be defined by appropriate programming code executed by signal processing circuitry 28.
As shown in Figure 3, a filtering sequence 52 for spatially-selective edge enhancement of a discrete pixel i~r~age includes a decomposition phase, designated generally by the reference muneral 54, an enhancement phase 56, and a reconstruction phase 58. Decomposition phase 54 involves progressive decomposition of the image data into a series of spatial frequency images of progressively lower spatial frequency. As used herein, the tenor spatial frequency relates generally to the number of individual image components within a defined image area. Through the progressive decomposition of the image, filtering sequence 52 generates a series of decomposed images having different :spatial frequencies. Each of these images is then enhanced by application of a gain image during enhancement phase 56. The resulting images are then recombined during reconstruction phase 58 to complete an enhanced image. The individual processing of the different spatial frequency levels of the decomposed images is represented in Figure 3 by reference numerals 60, 62, fi4, 66 and 68, with 68 being a final level of decomposition and processing;. As noted, each progressive step or level produces an intermediate image having a progressively lower spatial frequency, as noted at reference numerals 70, 72, 74, 76, and 78.
As shown in Figure 3, processing blocks for each spatial frequency level in decomposition phase 54 include a decimator (LP~L) 80 and an interpolator (LPT) 82, represented by low-pass filters (only the filters of the highest spatial frequency level 60 being labeled in Figure 3). These filters progressively combine data representing neighboring pixels of the immediately preceding image level in .a manner generally known in the art.
The filters output a first, filtered imal;e which is employed in subsequent decomposition levels, as well as a second image which is subtracted from the input image for that spatial frequency level at a summer 84. Thus, at a highest level 60, the original processed image data 86 is fed to the decimating low pass filter 80, as well as to summer 84. The original image data is filtered at blocks 80 and 82, with the output of block 80 being passed on as input data to level 62. In this nranner, the original processed image data is decomposed into the progressively lower spatial frequency images.
In a present embodiment, each subsequent step or level in the decomposition phase produces an image which is half' the size of the input image for that level. However, other reduction rates rnay be employed. Moreover, while the decomposition phase filtering may terminate at any desired level, in a present embodiment, such filtering is carried out until a single pixel ( 1 X 1 ) image is attained at lowest level 68.
The filtering sequence 52 then proceeds by applying (e.g., multiplying) each of the spatial frequency images by a gain image in enhancement phase 56. The gain images, designated in Figure 3 by the numerals 88, 90, 92, 94, and 96, are assigned values which enhance specific features of the image, such as edges. Nforeover, certain of the gain images are generated from data from lower spatial frequency levels. Thus, gain image 88, designated G1, is generated from image 72 resulting from the decomposition phase at level 62. Similarly, gain image 90, designated G2, is generated from image 74 resulting from processing at level 64. Any suitable number of levels in the filtering sequence 52 may employ such lower frequency-based gain images. For e:Kample, it has been found that in an image having original dimensions of 2,048 X 2,048 pixels, such lower spatial frequency-based gain images produce satisfactory edge enhancement when used at the first three levels of decomposition, such as levels 60, 62, and 64 in Figure 3. Other lower levels may also employ this technique, or may be multiplied by a gain image having values of unity, or other values, either greater or less than unity.
Following application of the gain images at enhancement phase 56, the resulting enhanced images are recombined in reconstruction phase 58.
As shown in Figure 3, each step or level in the reconstruction phase includes interpolation of output data from an immediately preceding (i.e., lower) level in a block 98 (see level 60 in Figure 3), The resulting image is added to the enhanced image in a summer 100 at each level. The resulting enhanced and reconstructed image data is then forw~~rded to the next (i.e., higher) level, until a reconstructed enhanced image 102 of the original dimensions is obtained.
It should be noted that, in particular applications, the image data need not be fully decomposed or reconstructed to the lowest levels. For example, where a series of lower spatial frequency images are to have a similar gain applied in phase 56, decomposition may be stopped at that level. The gain is then applied to the image at the terminal level, and reconstruction may be begin. Similarly, while as noted above, the higher level gains are derived from lower spatial frequency images resulting from decomposition, these could also be obtained based upon reconstructed images from a lower spatial frequency level (i.e., following application of the gain images during enhancement phase 56).
In the present embodiment, the lower spatial frequency-based gains are derived as follows. Following .decomposition of the original filtered image data, signal processing circuit 28 applies an algorithm which may be expressed:
G; (x,y) = max ( 1.0, E; ~k S(O~ (x,y))), where G; (x,y) is the value for each pixel in the gain image of level i to be determined, E; is a user-defined edge strength value greater than 1.0, S is a spatial sensitivity function which may be set by the user, and O~ (x,y) is the value of each pixel in a lower spatial frequency level image j. It should be noted, thus, that the foregoing algorithm effectively maintains a value of unity in the derived gain images when the resulting computed value is less than unity. Also, the algorithm permits the user to effectively define a desired edge strength value, tending to increase or decrease the relative sharpness of the edges and the reconstructed enhanced image. Similarly, a spatial sensitivity function may be selected by the user to regulate the amount of influence in the gain calculation exerted by values of the lower spatial frequency image. The spatial sensitivity function could include linear functions, piecewise linear functions a sigmoidal function, and so forth. In a present embodiment, the values of S resulting from the function S may be set anywhere between 0 and 1. Moreover, the edge strength value or the spatial sensitivity function, or both, may be unique to specific spatial frequency levels (i.e., unique to individual levels 60, 62, and 64 in Figure 3). Also, the particular level at which such spatial frequency-based gain images is derived may be a parameter which is user configurable. Thus, the user may select spatial frequency-based enhancement for one, two, three, or more levels, depending upon the enhancement desired and the frequency levels at which useful information or noise is anticipated.
G; (x,y) = max ( 1.0, E; ~k S(O~ (x,y))), where G; (x,y) is the value for each pixel in the gain image of level i to be determined, E; is a user-defined edge strength value greater than 1.0, S is a spatial sensitivity function which may be set by the user, and O~ (x,y) is the value of each pixel in a lower spatial frequency level image j. It should be noted, thus, that the foregoing algorithm effectively maintains a value of unity in the derived gain images when the resulting computed value is less than unity. Also, the algorithm permits the user to effectively define a desired edge strength value, tending to increase or decrease the relative sharpness of the edges and the reconstructed enhanced image. Similarly, a spatial sensitivity function may be selected by the user to regulate the amount of influence in the gain calculation exerted by values of the lower spatial frequency image. The spatial sensitivity function could include linear functions, piecewise linear functions a sigmoidal function, and so forth. In a present embodiment, the values of S resulting from the function S may be set anywhere between 0 and 1. Moreover, the edge strength value or the spatial sensitivity function, or both, may be unique to specific spatial frequency levels (i.e., unique to individual levels 60, 62, and 64 in Figure 3). Also, the particular level at which such spatial frequency-based gain images is derived may be a parameter which is user configurable. Thus, the user may select spatial frequency-based enhancement for one, two, three, or more levels, depending upon the enhancement desired and the frequency levels at which useful information or noise is anticipated.
Figure 4 illustrates steps in exennplary control logic for carrying out the foregoing processing sequence. As ;>hown in Figure 4, the control logic, designated by reference numeral 104, begins with acquisition of the pixel data at step 106, as described above. ~Che data is then filtered at step 108 through data acquisition circuitry 22, interface circuitry 26, and signal processing circuitry 28, such as to adjust the dynamic range of the input data, and so forth. At step 110, the image is decomposed through a level of spatial frequency as indicated above in phase 'i4 of Figure 3. At step 112, signal processing circuit 28 determines whether the image has been decomposed through the desired number of spatial frequency levels, and if not, returns to step 110 for fi,uther decomposition. Once the image has been decomposed through the desired spatial frequency level, the process advances to step 114.
At step 114, signal processing circuit 28 reads any configuration parameters which have been input by tt~e user or stored in memory circuit 30. As noted above, these may include wl edge strength parameter, a spatial sensitivity function, and so forth. In addition, such parameters may include predetermined gains to be applied to specific spatial frequency images, particularly to the lower level images. Flased upon the decomposed images and upon the parameters read at step 114, signal processing circuit 28 generates gain images (such as image;~ 88, 90, and 92 in Figure 3) as indicated at step 116. As noted above, this step preferably includes application of an algorithm whereby tt~e gain image for levels at which particular features of interest are to be enhanced are generated based upon one or more lower spatial frequency level image. At step 118, these gain images are applied to the decomposed image at the desired spatial frequency levels, such as at the upper levels where valuable edge data would be expected to be found. Subsequently, the process continues as indicated at step 120, by reconstruction of the enhanced images through the processing described above with respect to Figure 3.
At step 114, signal processing circuit 28 reads any configuration parameters which have been input by tt~e user or stored in memory circuit 30. As noted above, these may include wl edge strength parameter, a spatial sensitivity function, and so forth. In addition, such parameters may include predetermined gains to be applied to specific spatial frequency images, particularly to the lower level images. Flased upon the decomposed images and upon the parameters read at step 114, signal processing circuit 28 generates gain images (such as image;~ 88, 90, and 92 in Figure 3) as indicated at step 116. As noted above, this step preferably includes application of an algorithm whereby tt~e gain image for levels at which particular features of interest are to be enhanced are generated based upon one or more lower spatial frequency level image. At step 118, these gain images are applied to the decomposed image at the desired spatial frequency levels, such as at the upper levels where valuable edge data would be expected to be found. Subsequently, the process continues as indicated at step 120, by reconstruction of the enhanced images through the processing described above with respect to Figure 3.
Claims (21)
1. A method for enhancing an image defined by a matrix of discrete pixels, the method comprising the steps of:
generating data representative of the pixels defining the image;
at least partially decomposing the image into a plurality of spatial frequency images including at least a first frequency image and a second frequency image;
generating gain data for the first frequency image based upon the second frequency image;
generating an enhanced first frequency image by applying the gain data to the first frequency image; and reconstructing an enhanced image based upon the enhanced first frequency image and the second frequency image.
generating data representative of the pixels defining the image;
at least partially decomposing the image into a plurality of spatial frequency images including at least a first frequency image and a second frequency image;
generating gain data for the first frequency image based upon the second frequency image;
generating an enhanced first frequency image by applying the gain data to the first frequency image; and reconstructing an enhanced image based upon the enhanced first frequency image and the second frequency image.
2. The method of claim 1, wherein the image is decomposed by successive low pass filtering to obtain a plurality of spatial frequency images of successively smaller matrix size.
3. The method of claim 2, wherein a plurality of greater frequency images are generated, gain data are generated for each greater frequency image, and gain data for each greater frequency image is applied to the respective greater frequency image to generate a plurality of enhanced higher frequency images.
4. The method of claim 3, wherein the gain data for each respective greater frequency image is generated from at least one frequency image at a spatial frequency lower than the respective higher frequency image.
5. The method of claim 1, comprising the step of applying at least one gain to the second frequency image prior to reconstruction of the enhanced image.
6. The method of claim 1, wherein the gain data is at least partially based upon an operator selected spatial sensitivity function.
7. The method of claim 1, wherein a number of frequency levels into which the image is decomposed for generation of the gain data is operator adjustable.
8. A method for enhancement of edges, in a discrete pixel image, the method comprising the steps of:
acquiring a set of pixel data representative of the image;
filtering the pixel data into a series of spatial frequency level images, the series including at least one higher frequency level image and one lower frequency level image;
deriving a gain image from the lower frequency level image;
combining the gain image and the higher frequency level image to obtain an enhanced higher frequency level image; and reconstructing an edge enhanced image based upon the enhanced higher frequency level image and the lower frequency level image.
acquiring a set of pixel data representative of the image;
filtering the pixel data into a series of spatial frequency level images, the series including at least one higher frequency level image and one lower frequency level image;
deriving a gain image from the lower frequency level image;
combining the gain image and the higher frequency level image to obtain an enhanced higher frequency level image; and reconstructing an edge enhanced image based upon the enhanced higher frequency level image and the lower frequency level image.
9. The method of claim 8, wherein the set of pixel data is acquired via a digital x-ray imaging system.
10. The method of claim 8, wherein the gain image is derived from the lower frequency level image by applying values representative of the lower frequency level image into a sensitivity function.
11. The method of claim 10, wherein the gain image is derived from the lower frequency level image by applying values representative of
12 the lower frequency level image into a sensitivity function and multiplying by an edge enhancement value.
12. The method of claim 8, wherein a plurality of higher frequency level images are included in the series, and wherein dissimilar gain images are applied to each higher frequency level image to obtain respective enhanced higher frequency level images.
12. The method of claim 8, wherein a plurality of higher frequency level images are included in the series, and wherein dissimilar gain images are applied to each higher frequency level image to obtain respective enhanced higher frequency level images.
13. The method of claim 12, wherein a first higher frequency level image is derived from the image and successive higher frequency level images are derived therefrom by successive low pass filtering.
14. The method of claim 13, wherein the number of higher frequency level images derived is an operator adjustable parameter.
15. The method of claim 8, comprising the futher step of applying a gain to the lower frequency level image prior to reconstruction of the edge enhanced image.
16. The method of claim 8, wherein the lower frequency level image is generated by successively decomposing an initial lower frequency level image into a series of lower frequency level images, and reconstructing the series of lower frequency level images into a processed lower frequency level image, and wherein the gain image is derived from the processed lower frequency level image.
17. The method of claim 16, wherein dissimilar gains are applied to each lower frequency level image prior to reconstructing the series of lower frequency level images into a processed lower frequency level image.
18. A discrete pixel imaging system comprising:
a data acquisition unit configured to generate pixel values representative of pixels in a discrete pixel image; and a signal processing circuit configured to receive the pixel values, to decompose the image into a series of spatial frequency images, to derive a gain image for a first of the spatial frequency images from a second of the spatial frequency images, to apply the gain image to the first frequency image to obtain an edge enhanced image, and to reconstruct a filtered image based upon the edge enhanced image and the second frequency image.
a data acquisition unit configured to generate pixel values representative of pixels in a discrete pixel image; and a signal processing circuit configured to receive the pixel values, to decompose the image into a series of spatial frequency images, to derive a gain image for a first of the spatial frequency images from a second of the spatial frequency images, to apply the gain image to the first frequency image to obtain an edge enhanced image, and to reconstruct a filtered image based upon the edge enhanced image and the second frequency image.
19. The imaging system of claim 18, wherein the data acquisition unit includes an x-ray radiation source and a discrete pixel detector unit for receiving a portion of an x-ray emission from the x-ray radiation source.
20. The imaging system of claim 18, further including an operator input device for receiving operator selectable parameter data, and wherein the signal processing circuit is configured to derive the gain image based at least partially on parameter data input via the input device.
21. The imaging system of claim 18, wherein the signal processing circuit is configured to decompose the image into a plurality of spatial frequency images, and to derive dissimilar gain images for at least two of the spatial frequency images from respective lower spatial frequency images.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/113,653 US6175658B1 (en) | 1998-07-10 | 1998-07-10 | Spatially-selective edge enhancement for discrete pixel images |
US09/113,653 | 1998-07-10 |
Publications (1)
Publication Number | Publication Date |
---|---|
CA2276634A1 true CA2276634A1 (en) | 2000-01-10 |
Family
ID=22350739
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002276634A Abandoned CA2276634A1 (en) | 1998-07-10 | 1999-06-25 | Spatially-selective edge enhancement for discrete pixel images |
Country Status (4)
Country | Link |
---|---|
US (1) | US6175658B1 (en) |
EP (1) | EP0971316A3 (en) |
JP (1) | JP2000059631A (en) |
CA (1) | CA2276634A1 (en) |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10017551C2 (en) * | 2000-04-08 | 2002-10-24 | Carl Zeiss Vision Gmbh | Process for cyclic, interactive image analysis and computer system and computer program for executing the process |
US6956975B2 (en) | 2001-04-02 | 2005-10-18 | Eastman Kodak Company | Method for improving breast cancer diagnosis using mountain-view and contrast-enhancement presentation of mammography |
US20020196965A1 (en) * | 2001-06-22 | 2002-12-26 | Wallace Edward S. | Image transformation and analysis system and method |
DE60202588T2 (en) * | 2002-02-22 | 2006-01-05 | Agfa-Gevaert N.V. | Method for noise reduction |
US7167598B2 (en) | 2002-02-22 | 2007-01-23 | Agfa-Gevaert N.V. | Noise reduction method |
US7362915B2 (en) | 2002-02-22 | 2008-04-22 | Agfa Healthcare, Nv | Method for enhancing the contrast of an image |
DE60214967T2 (en) * | 2002-02-22 | 2007-07-26 | Agfa-Gevaert | Method for improving the contrast of an image |
JP2004005364A (en) * | 2002-04-03 | 2004-01-08 | Fuji Photo Film Co Ltd | Similar image retrieval system |
US6662088B1 (en) | 2002-06-28 | 2003-12-09 | General Electric Company | Methods and systems for inspecting aircraft fuselage frames |
US6901135B2 (en) * | 2003-08-28 | 2005-05-31 | Bio-Imaging Research, Inc. | System for extending the dynamic gain of an X-ray detector |
DE10343496B4 (en) * | 2003-09-19 | 2015-08-06 | Siemens Aktiengesellschaft | Correction of an X-ray image taken by a digital X-ray detector and calibration of the X-ray detector |
US20050243334A1 (en) * | 2004-04-28 | 2005-11-03 | Konica Minolta Medical & Graphic, Inc. | Image processing method, image processing apparatus and image processing program |
JP4724124B2 (en) * | 2004-10-08 | 2011-07-13 | パナソニック株式会社 | Image processing apparatus and image processing program |
DE102005028892A1 (en) * | 2005-06-22 | 2006-12-28 | Siemens Ag | Processing a two-dimensional initial image by dismantling initial image into partial images and residual image, assigning one of the partial images as pilot image and assigning pilot frequency, and determining and summing weighting factors |
US7653230B2 (en) * | 2006-02-21 | 2010-01-26 | General Electric Company | Methods and systems for image reconstruction using low noise kernel |
CN101416218A (en) * | 2006-03-29 | 2009-04-22 | 皇家飞利浦电子股份有限公司 | Method and system for improving visual quality of an image signal |
JP5467958B2 (en) * | 2010-07-16 | 2014-04-09 | 富士フイルム株式会社 | Radiation image processing apparatus and method, and program |
CN104166974B (en) * | 2013-08-01 | 2015-05-13 | 上海联影医疗科技有限公司 | CT locating film image enhancing method and CT locating film image enhancing device |
JP5918198B2 (en) * | 2013-11-26 | 2016-05-18 | 日立アロカメディカル株式会社 | Ultrasonic diagnostic equipment |
US10362288B2 (en) * | 2015-01-09 | 2019-07-23 | Sony Corporation | Method and system for improving detail information in digital images |
CN109801247A (en) * | 2017-11-16 | 2019-05-24 | 通用电气公司 | Medical imaging procedure, computer program and computer storage medium |
JP2023027474A (en) * | 2021-08-17 | 2023-03-02 | キヤノンメディカルシステムズ株式会社 | X-ray diagnostic apparatus, medical image processing apparatus, and program |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2555430B2 (en) | 1988-11-07 | 1996-11-20 | 大日本スクリーン製造 株式会社 | Image edge enhancement method |
US5022091A (en) | 1990-02-28 | 1991-06-04 | Hughes Aircraft Company | Image processing technique |
EP0527525B1 (en) | 1991-08-14 | 1996-10-02 | Agfa-Gevaert N.V. | Method and apparatus for contrast enhancement |
DE69214229T2 (en) | 1991-08-14 | 1997-04-30 | Agfa Gevaert Nv | Method and device for improving the contrast of images |
DE69331719T2 (en) | 1992-06-19 | 2002-10-24 | Agfa Gevaert Nv | Method and device for noise suppression |
EP0610604A1 (en) | 1993-02-11 | 1994-08-17 | Agfa-Gevaert N.V. | Multiple processing of radiographic images based on a pyramidal image decomposition |
EP0610603B1 (en) | 1993-02-11 | 1999-09-08 | Agfa-Gevaert N.V. | Fast interactive off-line processing method for radiographic images |
JPH0922460A (en) * | 1995-02-09 | 1997-01-21 | Fuji Photo Film Co Ltd | Image processing method and device therefor |
GB2311432B (en) * | 1996-03-20 | 2000-05-03 | Sony Uk Ltd | Method and apparatus for processing an input image |
EP0923760B1 (en) | 1997-06-06 | 2005-11-16 | Koninklijke Philips Electronics N.V. | Noise reduction in an image |
-
1998
- 1998-07-10 US US09/113,653 patent/US6175658B1/en not_active Expired - Fee Related
-
1999
- 1999-06-25 CA CA002276634A patent/CA2276634A1/en not_active Abandoned
- 1999-07-08 EP EP99305432A patent/EP0971316A3/en not_active Withdrawn
- 1999-07-08 JP JP11193840A patent/JP2000059631A/en not_active Withdrawn
Also Published As
Publication number | Publication date |
---|---|
JP2000059631A (en) | 2000-02-25 |
EP0971316A3 (en) | 2000-01-19 |
US6175658B1 (en) | 2001-01-16 |
EP0971316A2 (en) | 2000-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6175658B1 (en) | Spatially-selective edge enhancement for discrete pixel images | |
US6782137B1 (en) | Digital image display improvement system and method | |
US5467404A (en) | Method and apparatus for contrast enhancement | |
EP0527525B1 (en) | Method and apparatus for contrast enhancement | |
US7079700B2 (en) | Image processing apparatus, image processing method, storage medium, and program | |
US20020054713A1 (en) | Image processing apparatus and its method, program, and storage medium | |
US20070086671A1 (en) | Image processing apparatus | |
US7149358B2 (en) | Method and system for improving contrast using multi-resolution contrast based dynamic range management | |
EP2904580B1 (en) | Real-time image processing for optimizing sub-images views | |
US7155044B2 (en) | Multiscale gradation processing method | |
JP2005160107A (en) | Resolution adaptive image filtering system and method | |
US7298383B2 (en) | Method and user interface for modifying at least one of contrast and density of pixels of a processed image | |
US5757952A (en) | Method of displaying (a) radiographic image (s) | |
US7626592B1 (en) | Multiscale gradation processing method | |
US20050110802A1 (en) | Method and apparatus for segmentation-based image operations | |
US7298373B2 (en) | Method and user interface for modifying at least one of contrast and density of pixels of a processed image | |
JP2002304622A (en) | Device and method for image processing, storage medium, and program | |
JP2008142178A (en) | Radiological image processing apparatus | |
JP2002183727A (en) | Image processing device | |
EP1341125A2 (en) | Multiscale gradation processing method | |
WO2008056305A1 (en) | Noise reduction of an image signal | |
JPH10105701A (en) | Method and device for radio graph emphasis processing | |
CN106780400A (en) | A kind of image processing method and device | |
JP2003259125A (en) | Multiscale gradation processing method | |
EP1486914B1 (en) | Method and user interface for modifying at least one of contrast and density of pixels of a processed image. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
FZDE | Discontinued |