US20150093002A1 - Method of compensating ultrasound image - Google Patents
Method of compensating ultrasound image Download PDFInfo
- Publication number
- US20150093002A1 US20150093002A1 US14/164,588 US201414164588A US2015093002A1 US 20150093002 A1 US20150093002 A1 US 20150093002A1 US 201414164588 A US201414164588 A US 201414164588A US 2015093002 A1 US2015093002 A1 US 2015093002A1
- Authority
- US
- United States
- Prior art keywords
- brightness
- compensation
- image
- values
- ultrasound image
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000002604 ultrasonography Methods 0.000 title claims abstract description 55
- 239000000463 material Substances 0.000 claims abstract description 23
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 238000003759 clinical diagnosis Methods 0.000 description 9
- 238000012285 ultrasound imaging Methods 0.000 description 9
- 230000007423 decrease Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/94—Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
-
- G06T5/008—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52033—Gain control of receivers
-
- G06T7/0081—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
-
- 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/10132—Ultrasound 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/20172—Image enhancement details
Definitions
- the present invention relates to a method of compensating ultrasound image, and more particularly to a method which demarcates the main material region on the ultrasound image and then compensates the ultrasound image and enhances image contrast.
- acoustic detection technology such as ultrasound imaging technology
- ultrasound imaging technology has the advantages of low price, non-invasive, no danger of radiation, real-time imaging, mm level spatial resolution, portability, and blood detectable, and thus is widely used in clinical diagnosis of various departments.
- the ultrasound imaging technology creates images based on scattering and reflection of sound wave. Because strength of echo of the sound wave produced by the probe decreases attending with the increasing of medium depth, brightness and effective dynamic range of the resulted ultrasound image would be declined accordingly. In addition, decline of brightness of the ultrasound images for various organs or positions of human body would be different, which influences reliability of clinical diagnosis.
- the clinical diagnosis nowadays usually applies a predetermined linear brightness compensation to the ultrasound image and then adjust the brightness gain and dynamic range based on the experience of the operator.
- the process wastes lots of time and is quite inconvenient.
- ultrasound imaging applied to various organs or positions of human body would have different brightness decline, which influences reliability of clinical diagnosis.
- the ultrasound image also needs the operators to adjust brightness gain and dynamic range based on their experience, which derives the issue of low efficiency and operational inconvenience.
- the provided method demarcates the main material regions on the ultrasound image and then compensates the ultrasound image based on the depth and enhances image contrast.
- reliability of clinical diagnosis and operational convenience can be improved.
- a method of compensating ultrasound image comprises the steps of: (a) receiving an ultrasound image with a plurality of first pixel brightness values; (b) demarcating a plurality of main material regions of the ultrasound image to generate a first attenuation curve of the main material regions; (c) executing a full compensation process to generate a full compensation image with a plurality of second pixel brightness values corresponding to the first pixels according to the first attenuation curve; and (d) generating a brightness comparison table including a plurality of brightness zones corresponding to a plurality of distinct brightness compensation values according to the first pixel brightness values and the second pixel brightness values.
- the present invention further comprises the step of: (e) executing a linear sum process to the ultrasound image and the full compensation image generate a compensation image with a second attenuation curve according to the brightness compensation values; (f) generating a compensation curve according to the first attenuation curve and the second attenuation curve; (g) accessing a depth, a first brightness enhance value, and a second brightness enhance value of the compensation curve; and (h) executing a space comparison process to the compensation image to generate a better(modified) compensation image according to the depth, the first brightness enhance value, and the second brightness enhance value.
- the linear sum process of the above mentioned step (e) is executed by using a 2D weight mask.
- the brightness compensation values are ranged between 0 to 1
- the space comparison process of step (g) is executed by using an adaptive 2D gamma filter.
- the depth is defined as d
- the first brightness enhance value is defined as r1
- the second brightness enhance value is defined as r2
- the threshold value rth is a predetermined value in brightness.
- the contrast enhance value ⁇ is further identified by using the corresponding first brightness value r1 and the second brightness value r2, and also the threshold value rth.
- the method of compensating ultrasound image provided in the present invention is capable of compensating brightness along depth direction automatically, thus, the unwanted influence of brightness decline because of medium depth can be effectively reduced to enhance reliability of clinical diagnosis.
- the method provided in the present invention is capable to adaptively compensate brightness automatically, which should be an effective approach to resolve the problem of inconvenience and image variation of transitional ultrasound imaging system in which the compensation curve needs to be adjusted manually.
- the operation of ultrasound imaging in clinical diagnosis can be more efficiency.
- FIG. 1 and FIG. 1A are flow charts showing the method of compensating ultrasound image in accordance with a preferred embodiment of the present invention.
- FIG. 2 is a schematic view showing the ultrasound image in accordance with a preferred embodiment of the present invention.
- FIG. 2A is a schematic view showing the main material region in accordance with a preferred embodiment of the present invention.
- FIG. 3 is a schematic view showing the first attenuation curve, the second attenuation curve, and the third attenuation curve in accordance with a preferred embodiment of the present invention.
- FIG. 4 is a schematic view showing a brightness comparison table in accordance with a preferred embodiment of the present invention.
- FIG. 5 is a schematic view showing the comparison of the first attenuation curve and the second attenuation curve.
- FIG. 6 is a schematic view showing the compensation image in accordance with a preferred embodiment of the present invention.
- FIG. 7 is a diagram showing the first attenuation curve and the second attenuation curve in accordance with a preferred embodiment of the present invention.
- FIG. 7A is a diagram showing the compensation curve in accordance with a preferred embodiment of the present invention.
- FIG. 8 is a schematic view showing the better (modified) compensation image in accordance with a preferred embodiment of the present invention.
- FIG. 1 and FIG. 1A are flow charts showing the method of compensating ultrasound image in accordance with a preferred embodiment of the present invention
- FIG. 2 is a schematic view showing the ultrasound image in accordance with a preferred embodiment of the present invention
- FIG. 2A is a schematic view showing the main material region in accordance with a preferred embodiment of the present invention
- FIG. 3 is a schematic view showing the first attenuation curve, the second attenuation curve, and the third attenuation curve in accordance with a preferred embodiment of the present invention
- FIG. 4 is a schematic view showing a brightness comparison table in accordance with a preferred embodiment of the present invention
- FIG. 5 is a schematic view showing the comparison of the first attenuation curve and the second attenuation curve
- FIG. 5 is a schematic view showing the comparison of the first attenuation curve and the second attenuation curve
- FIG. 6 is a schematic view showing the compensation image in accordance with a preferred embodiment of the present invention
- FIG. 7 is a diagram showing the first attenuation curve and the second attenuation curve in accordance with a preferred embodiment of the present invention
- FIG. 7A is a diagram showing the compensation curve in accordance with a preferred embodiment of the present invention
- FIG. 8 is a schematic view showing the better (modified) compensation image in accordance with a preferred embodiment of the present invention.
- the method of compensating ultrasound image provided in the preferred embodiment of the present invention can be applied to ultrasound image of all the imaging medium, such as the organs of human body or blood flows.
- the method comprises the steps as below:
- Step S 101 receiving an ultrasound image with a plurality of first pixel brightness values
- Step S 102 demarcating a plurality of main material regions of the ultrasound image to generate a first attenuation curve of the main material regions;
- Step S 103 executing a full compensation process to generate a full compensation image with a plurality of second pixel brightness values
- Step S 104 generating a brightness comparison table including a plurality of brightness zones corresponding to a plurality of distinct brightness compensation values according to the first pixel brightness values and the second pixel brightness values;
- Step S 105 executing a linear sum process to the ultrasound image and the full compensation image to generate a compensation image with a second attenuation curve according to the brightness compensation values;
- Step S 106 generating a compensation curve according to the first attenuation curve and the second attenuation curve;
- Step S 107 accessing a depth, a first brightness enhance value, and a second brightness enhance value of the compensation curve.
- Step S 108 executing a space comparison process to the compensation image to generate a better (modified) compensation image according to the depth, the first brightness enhance value, and the second brightness enhance value.
- step S 101 is carried out to receive an ultrasound image including a plurality of first pixel brightness values.
- the ultrasound image 100 received in this step is an original ultrasound image as shown in FIG. 2 , which is the original image based on the echo of the ultrasound wave radiated to the medium such as tissue of human organ.
- the ultrasound image is composed of a plurality of first pixels (not labeled in the figure) corresponding to a plurality of first pixel brightness values.
- step S 102 is carried out to demarcate a plurality of main material regions of the ultrasound image so as to generate a first attenuation curve.
- the regions of main material are demarcated as the main material regions A, and the main material regions A are regarded as portion of the main tissue.
- determining the regions along the depth direction from the main material regions based on energy distribution and average energy value.
- the portions with similar condition are determined as the same tissue. Theoretically, where there has sudden change of brightness can be determined as of different tissues, and the continuation of main material region A is stopped thereby. However, if it is determined to be of the same tissue, the main material region A keeps extending downward.
- an attenuation curve along the depth direction of the main material region A is formed as the first attenuation curve 300 as shown in FIG. 3 .
- step S 103 is carried out to execute a full compensation process to the first attenuation curve so as to generate a full compensation image with a plurality of second pixel brightness values.
- a full compensation process is executed to the first attenuation curve to have the first attenuation curve reach the same brightness.
- the full compensation process is well known for those skilled in the art and thus is not repeated here.
- the full compensation image 300 as shown in FIG. 4 is generated.
- the full compensation image 300 includes a plurality of second pixels (not labeled in the figure) corresponding to a plurality of second pixel brightness values for generating the full compensation curve 400 as shown in FIG. 3 .
- step S 104 is carried out to generate a brightness comparison table including a plurality of brightness zones corresponding to a plurality of distinct brightness compensation values according to the first pixel brightness values and the second pixel brightness values.
- SNR signal to noise ratio
- step S 104 is carried out to generate a brightness comparison table including a plurality of brightness zones corresponding to a plurality of distinct brightness compensation values according to the first pixel brightness values and the second pixel brightness values.
- SNR signal to noise ratio
- the present step is carried out prior to the compensation step to generate the brightness comparison table as shown in FIG. 5 based on the first pixel brightness values included in the ultrasound image 100 and the second pixel brightness values included in the full compensation image 300 .
- the brightness comparison table is of logarithmic scale.
- the horizontal axis value corresponds to the logarithm value of brightness before compensation, which is the logarithm value of the first pixel brightness value.
- the vertical axis value corresponds to the logarithm value of brightness after compensation, which is the logarithm value of the second pixel brightness value.
- the brightness comparison table is divided into a plurality of brightness zones B.
- the grouping process of brightness zones B can be executed automatically according to different brightness areas using algorithm.
- each brightness zones B is with respect to a distinct brightness compensation value.
- the brightness compensation value defined as ⁇ is ranged between 0 to 1. The greater brightness compensation value implies that the output is more close to the high quality image. The smaller brightness compensation value implies that the output is more close to the image with even brightness.
- step S 105 a linear sum process is carried out to the ultrasound image and the full compensation image so as to generate a compensation image with a second attenuation curve according to the brightness compensation values.
- various tissues included in the ultrasound image 100 received in step S 101 which might have different attenuation values during the ultrasound imaging process.
- the main material regions A the black region in the figure
- the demarcated regions and the un-demarcated regions might have different attenuation values in fact. It is unreasonable to use the attenuation values of the demarcated regions as those for the un-demarcated regions.
- step S 104 is executed prior to step S 105 to reduce the unwanted influence of wrong attenuation compensation used in the un-demarcated regions.
- a linear sum process is executed to the ultrasound image 100 and the full compensation image 300 according to the brightness compensation values of each of the brightness zones B so as to generate a compensation image 500 as show in FIG. 6 , which having a second attenuation curve 600 as shown in FIG. 3 .
- the linear sum process is executed by using an algorithm of 2D weight mask.
- the compensation image 500 includes a plurality of third pixels (not labeled in this figure) corresponding to a plurality of third pixel brightness values. The third pixel brightness values should be with respect to the first pixel brightness values and the second pixel brightness values because the compensation image 500 is generated based on the ultrasound image 100 and the full compensation image 300 using the linear sum process.
- step S 104 and step S 105 are introduced to compensate the drawback that the un-demarcated tissues/regions are missed in the attenuation calculation of step S 102 and step S 103 so as to reduce the influence effectively.
- step S 106 is carried out to generate a compensation curve according to the first attenuation curve and the second attenuation curve.
- the first attenuation curve 200 and the second attenuation curve 600 are modified to generate the first attenuation average curve 700 and the second attenuation average curve 800 as shown in FIG. 7 by mapping the vertical axis from brightness to average value.
- the difference of the two curves 700 and 800 are used to generate the compensation curve 900 as shown in FIG. 7A .
- step S 107 is executed to access a depth, a first brightness enhance value, and a second brightness enhance value of the compensation curve 900 .
- this step fetches a point on the compensation curve 900 , which includes the depth, the first brightness enhance value, and the second brightness enhance value.
- the compensation point C is fetched with the depth as d, the first brightness enhance value as r1, and the second brightness enhance value as r2.
- the second brightness enhance value r2 is the average value of brightness enhancement of the whole image.
- the first attenuation average curve 700 shows brightness attenuation of each layer of the ultrasound image 100
- the second attenuation average curve 800 shows brightness attenuation of each layer of the ultrasound image after compensation.
- the second attenuation average curve 800 has greater compensation in compared with the first attenuation average curve at the deeper location.
- the compensation curve 900 of FIG. 7A the compensation changes attending with the increasing of depth. For a given depth value (depth d of the compensation point C), contrast of the compensation image 500 should be enhanced as the brightness is greater than that corresponding to the depth value.
- the step S 107 the needed contrast enhancement with respect to the depth.
- the step S 108 can be carried out to execute a space comparison process to the compensation image to generate a better (modified) compensation image according to the depth, the first brightness enhance value, and the second brightness enhance value.
- the space comparison process is carried out by using an adaptive 2D gamma filter so as to find out a contrast enhance value ⁇ based on a threshold value rth, the depth d, the first brightness enhance value r1, and the second brightness enhance value r2.
- the threshold value rth is a predetermined value in brightness.
- the contrast enhance value ⁇ is identified by using the corresponding first brightness value r1 and the second brightness value r2, and also the threshold value rth. Then, the contrast enhance value is used to generate the modified compensation image 1000 as shown in FIG. 8 to complete the compensation process.
- the method of compensating ultrasound image provided in the present invention is capable of compensating brightness along depth direction automatically, thus, the unwanted influence of brightness decline because of medium depth can be effectively reduced to enhance reliability of clinical diagnosis.
- the method provided in the present invention is capable to adaptively compensate brightness automatically, which should be an effective approach to resolve the problem of inconvenience and image variation of transitional ultrasound imaging system in which the compensation curve needs to be adjusted manually.
- the operation of ultrasound imaging in clinical diagnosis can be more efficiency, which would be quite useful for the industry.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
- The present invention relates to a method of compensating ultrasound image, and more particularly to a method which demarcates the main material region on the ultrasound image and then compensates the ultrasound image and enhances image contrast.
- Attending with the technological development and progress, acoustic detection technology, such as ultrasound imaging technology, has been widely used in the modern diagnostic procedure. In compared with other medical imaging systems being used in clinical medicine, such as X-ray, CT, MRI or nuclear medicine imaging, the ultrasound imaging technology has the advantages of low price, non-invasive, no danger of radiation, real-time imaging, mm level spatial resolution, portability, and blood detectable, and thus is widely used in clinical diagnosis of various departments.
- The ultrasound imaging technology creates images based on scattering and reflection of sound wave. Because strength of echo of the sound wave produced by the probe decreases attending with the increasing of medium depth, brightness and effective dynamic range of the resulted ultrasound image would be declined accordingly. In addition, decline of brightness of the ultrasound images for various organs or positions of human body would be different, which influences reliability of clinical diagnosis.
- In order to access the actual image without attenuation and correct diagnosis information, the clinical diagnosis nowadays usually applies a predetermined linear brightness compensation to the ultrasound image and then adjust the brightness gain and dynamic range based on the experience of the operator. The process wastes lots of time and is quite inconvenient.
- In view of the traditional technology, ultrasound imaging applied to various organs or positions of human body would have different brightness decline, which influences reliability of clinical diagnosis. In addition, the ultrasound image also needs the operators to adjust brightness gain and dynamic range based on their experience, which derives the issue of low efficiency and operational inconvenience.
- Accordingly, it is a main object of the present invention to provide a method of compensating ultrasound image to compensate the ultrasound image automatically. The provided method demarcates the main material regions on the ultrasound image and then compensates the ultrasound image based on the depth and enhances image contrast. Thus, reliability of clinical diagnosis and operational convenience can be improved.
- Based on the above mentioned object, a method of compensating ultrasound image is provided in accordance with an embodiment of the present invention. The method comprises the steps of: (a) receiving an ultrasound image with a plurality of first pixel brightness values; (b) demarcating a plurality of main material regions of the ultrasound image to generate a first attenuation curve of the main material regions; (c) executing a full compensation process to generate a full compensation image with a plurality of second pixel brightness values corresponding to the first pixels according to the first attenuation curve; and (d) generating a brightness comparison table including a plurality of brightness zones corresponding to a plurality of distinct brightness compensation values according to the first pixel brightness values and the second pixel brightness values.
- After executing step (d), the present invention further comprises the step of: (e) executing a linear sum process to the ultrasound image and the full compensation image generate a compensation image with a second attenuation curve according to the brightness compensation values; (f) generating a compensation curve according to the first attenuation curve and the second attenuation curve; (g) accessing a depth, a first brightness enhance value, and a second brightness enhance value of the compensation curve; and (h) executing a space comparison process to the compensation image to generate a better(modified) compensation image according to the depth, the first brightness enhance value, and the second brightness enhance value.
- In accordance with a preferred embodiment of the present invention, the linear sum process of the above mentioned step (e) is executed by using a 2D weight mask. In addition, the compensation image further comprises a plurality of third pixel brightness values corresponding to the first pixel brightness values and the second pixel brightness values, the brightness compensation values are defined as a, the first pixel brightness values are defined as P1, the second pixel brightness values are defined as P2, the third pixel brightness values are defined as P3, and the linear sum process has P3=α*P1+(1−α)P2. Moreover, the brightness compensation values are ranged between 0 to 1, and the space comparison process of step (g) is executed by using an adaptive 2D gamma filter.
- In accordance with a preferred embodiment of the present invention, the depth is defined as d, the first brightness enhance value is defined as r1, the second brightness enhance value is defined as r2, and a contrast enhance value γ is further identified by using the space comparison process to generate the modified compensation image according to a threshold value rth, the depth d, the first brightness enhance value r1, and the second brightness enhance value r2 such that γ=(r1−rth)/r2. The threshold value rth is a predetermined value in brightness. When depth is greater than the predetermined depth d, the contrast enhance value γ is further identified by using the corresponding first brightness value r1 and the second brightness value r2, and also the threshold value rth.
- The method of compensating ultrasound image provided in the present invention is capable of compensating brightness along depth direction automatically, thus, the unwanted influence of brightness decline because of medium depth can be effectively reduced to enhance reliability of clinical diagnosis. In addition, the method provided in the present invention is capable to adaptively compensate brightness automatically, which should be an effective approach to resolve the problem of inconvenience and image variation of transitional ultrasound imaging system in which the compensation curve needs to be adjusted manually. Thus, the operation of ultrasound imaging in clinical diagnosis can be more efficiency.
- The embodiments adopted in the present invention would be further discussed by using the flowing paragraph and the figures for a better understanding.
-
FIG. 1 andFIG. 1A are flow charts showing the method of compensating ultrasound image in accordance with a preferred embodiment of the present invention. -
FIG. 2 is a schematic view showing the ultrasound image in accordance with a preferred embodiment of the present invention. -
FIG. 2A is a schematic view showing the main material region in accordance with a preferred embodiment of the present invention. -
FIG. 3 is a schematic view showing the first attenuation curve, the second attenuation curve, and the third attenuation curve in accordance with a preferred embodiment of the present invention. -
FIG. 4 is a schematic view showing a brightness comparison table in accordance with a preferred embodiment of the present invention. -
FIG. 5 is a schematic view showing the comparison of the first attenuation curve and the second attenuation curve. -
FIG. 6 is a schematic view showing the compensation image in accordance with a preferred embodiment of the present invention. -
FIG. 7 is a diagram showing the first attenuation curve and the second attenuation curve in accordance with a preferred embodiment of the present invention. -
FIG. 7A is a diagram showing the compensation curve in accordance with a preferred embodiment of the present invention. -
FIG. 8 is a schematic view showing the better (modified) compensation image in accordance with a preferred embodiment of the present invention. - There are various embodiments of the method of compensating ultrasound image in accordance with the present invention, which are not repeated hereby. The preferred embodiments are mentioned in the following paragraph as an example. It should be understood by those skilled in the art that the preferred embodiments disclosed in the following paragraph are merely an example instead of restricting the scope of the invention itself.
-
FIG. 1 andFIG. 1A are flow charts showing the method of compensating ultrasound image in accordance with a preferred embodiment of the present invention,FIG. 2 is a schematic view showing the ultrasound image in accordance with a preferred embodiment of the present invention,FIG. 2A is a schematic view showing the main material region in accordance with a preferred embodiment of the present invention,FIG. 3 is a schematic view showing the first attenuation curve, the second attenuation curve, and the third attenuation curve in accordance with a preferred embodiment of the present invention,FIG. 4 is a schematic view showing a brightness comparison table in accordance with a preferred embodiment of the present invention,FIG. 5 is a schematic view showing the comparison of the first attenuation curve and the second attenuation curve,FIG. 6 is a schematic view showing the compensation image in accordance with a preferred embodiment of the present invention,FIG. 7 is a diagram showing the first attenuation curve and the second attenuation curve in accordance with a preferred embodiment of the present invention,FIG. 7A is a diagram showing the compensation curve in accordance with a preferred embodiment of the present invention, andFIG. 8 is a schematic view showing the better (modified) compensation image in accordance with a preferred embodiment of the present invention. - As shown, the method of compensating ultrasound image provided in the preferred embodiment of the present invention can be applied to ultrasound image of all the imaging medium, such as the organs of human body or blood flows. The method comprises the steps as below:
- Step S101: receiving an ultrasound image with a plurality of first pixel brightness values;
- Step S102: demarcating a plurality of main material regions of the ultrasound image to generate a first attenuation curve of the main material regions;
- Step S103: executing a full compensation process to generate a full compensation image with a plurality of second pixel brightness values;
- Step S104: generating a brightness comparison table including a plurality of brightness zones corresponding to a plurality of distinct brightness compensation values according to the first pixel brightness values and the second pixel brightness values;
- Step S105: executing a linear sum process to the ultrasound image and the full compensation image to generate a compensation image with a second attenuation curve according to the brightness compensation values;
- Step S106: generating a compensation curve according to the first attenuation curve and the second attenuation curve;
- Step S107: accessing a depth, a first brightness enhance value, and a second brightness enhance value of the compensation curve; and
- Step S108: executing a space comparison process to the compensation image to generate a better (modified) compensation image according to the depth, the first brightness enhance value, and the second brightness enhance value.
- After the process starts, step S101 is carried out to receive an ultrasound image including a plurality of first pixel brightness values. The
ultrasound image 100 received in this step is an original ultrasound image as shown inFIG. 2 , which is the original image based on the echo of the ultrasound wave radiated to the medium such as tissue of human organ. The ultrasound image is composed of a plurality of first pixels (not labeled in the figure) corresponding to a plurality of first pixel brightness values. - After completion of step S101, step S102 is carried out to demarcate a plurality of main material regions of the ultrasound image so as to generate a first attenuation curve. Concretely speaking, as shown in
FIG. 2 andFIG. 2A , the regions of main material are demarcated as the main material regions A, and the main material regions A are regarded as portion of the main tissue. Thereafter, determining the regions along the depth direction from the main material regions based on energy distribution and average energy value. The portions with similar condition are determined as the same tissue. Theoretically, where there has sudden change of brightness can be determined as of different tissues, and the continuation of main material region A is stopped thereby. However, if it is determined to be of the same tissue, the main material region A keeps extending downward. After demarcate the main material regions A, an attenuation curve along the depth direction of the main material region A is formed as thefirst attenuation curve 300 as shown inFIG. 3 . - After completion of step S102, step S103 is carried out to execute a full compensation process to the first attenuation curve so as to generate a full compensation image with a plurality of second pixel brightness values. Concretely speaking, in this step, a full compensation process is executed to the first attenuation curve to have the first attenuation curve reach the same brightness. The full compensation process is well known for those skilled in the art and thus is not repeated here. After the full compensation process, the
full compensation image 300 as shown inFIG. 4 is generated. Thefull compensation image 300 includes a plurality of second pixels (not labeled in the figure) corresponding to a plurality of second pixel brightness values for generating thefull compensation curve 400 as shown inFIG. 3 . - After completion of step S103, step S104 is carried out to generate a brightness comparison table including a plurality of brightness zones corresponding to a plurality of distinct brightness compensation values according to the first pixel brightness values and the second pixel brightness values. Wherein, signal to noise ratio (SNR) of
ultrasound image 100 is greater thefull compensation image 300, which implies better image quality. In contrast, brightness distribution of thefull compensation image 300 which undergoes a brightness enhancing process is more uniform than theultrasound image 100. Thus, to enhance compensation accuracy, the present step is carried out prior to the compensation step to generate the brightness comparison table as shown inFIG. 5 based on the first pixel brightness values included in theultrasound image 100 and the second pixel brightness values included in thefull compensation image 300. The brightness comparison table is of logarithmic scale. The horizontal axis value corresponds to the logarithm value of brightness before compensation, which is the logarithm value of the first pixel brightness value. The vertical axis value corresponds to the logarithm value of brightness after compensation, which is the logarithm value of the second pixel brightness value. - In addition, the brightness comparison table is divided into a plurality of brightness zones B. The grouping process of brightness zones B can be executed automatically according to different brightness areas using algorithm. Moreover, each brightness zones B is with respect to a distinct brightness compensation value. In the preferred embodiment of the present invention, the brightness compensation value defined as α is ranged between 0 to 1. The greater brightness compensation value implies that the output is more close to the high quality image. The smaller brightness compensation value implies that the output is more close to the image with even brightness.
- After completion of step S104, in step S105, a linear sum process is carried out to the ultrasound image and the full compensation image so as to generate a compensation image with a second attenuation curve according to the brightness compensation values. Concretely speaking, there might be various tissues included in the
ultrasound image 100 received in step S101, which might have different attenuation values during the ultrasound imaging process. Thus, although the main material regions A (the black region in the figure) are demarcated in step S102 for attenuation determination, the demarcated regions and the un-demarcated regions might have different attenuation values in fact. It is unreasonable to use the attenuation values of the demarcated regions as those for the un-demarcated regions. Thus, according to a preferred embodiment of the present invention, in order to make image brightness of different tissues in the compensated image reasonable, step S104 is executed prior to step S105 to reduce the unwanted influence of wrong attenuation compensation used in the un-demarcated regions. - In the present step S105, a linear sum process is executed to the
ultrasound image 100 and thefull compensation image 300 according to the brightness compensation values of each of the brightness zones B so as to generate acompensation image 500 as show inFIG. 6 , which having asecond attenuation curve 600 as shown inFIG. 3 . In detail, as a preferred embodiment, the linear sum process is executed by using an algorithm of 2D weight mask. In addition, thecompensation image 500 includes a plurality of third pixels (not labeled in this figure) corresponding to a plurality of third pixel brightness values. The third pixel brightness values should be with respect to the first pixel brightness values and the second pixel brightness values because thecompensation image 500 is generated based on theultrasound image 100 and thefull compensation image 300 using the linear sum process. - Take the first pixel brightness values as P1, the second pixel brightness values as P2, the third pixel brightness values as P3. The third pixel brightness values P3 will follow the function P3=α*P1+(1−α)P2 based on the linear sum process of a preferred embodiment of the present invention. Advantageously, because SNR becomes worse when doing brightness compensation, the un-demarcated tissues/regions defined in step S102, which have an average brightness different from that of the demarcated main material regions A, will have such problem and should be adequately compensated. The step S104 and step S105 are introduced to compensate the drawback that the un-demarcated tissues/regions are missed in the attenuation calculation of step S102 and step S103 so as to reduce the influence effectively.
- After completion of step S105, the step S106 is carried out to generate a compensation curve according to the first attenuation curve and the second attenuation curve. In this step, the
first attenuation curve 200 and thesecond attenuation curve 600 are modified to generate the first attenuationaverage curve 700 and the second attenuationaverage curve 800 as shown inFIG. 7 by mapping the vertical axis from brightness to average value. Then, the difference of the twocurves compensation curve 900 as shown inFIG. 7A . - After the completion of step S106, the step S107 is executed to access a depth, a first brightness enhance value, and a second brightness enhance value of the
compensation curve 900. Concretely speaking, this step fetches a point on thecompensation curve 900, which includes the depth, the first brightness enhance value, and the second brightness enhance value. As shown, the compensation point C is fetched with the depth as d, the first brightness enhance value as r1, and the second brightness enhance value as r2. - When the brightness compensation value is 1, which means no compensation, the first brightness enhance value should be 0, and when the brightness compensation value is not 1, which means the weight of the full compensation image is included, the brightness has certain enhancement and the enhance value, the first brightness enhance value r1, would be changed according attending with the depth. The second brightness enhance value r2 is the average value of brightness enhancement of the whole image. As shown in
FIG. 7 , the first attenuationaverage curve 700 shows brightness attenuation of each layer of theultrasound image 100, the second attenuationaverage curve 800 shows brightness attenuation of each layer of the ultrasound image after compensation. As shown, the second attenuationaverage curve 800 has greater compensation in compared with the first attenuation average curve at the deeper location. In addition, referring to thecompensation curve 900 ofFIG. 7A , the compensation changes attending with the increasing of depth. For a given depth value (depth d of the compensation point C), contrast of thecompensation image 500 should be enhanced as the brightness is greater than that corresponding to the depth value. - That is, the step S107 the needed contrast enhancement with respect to the depth. Then, the step S108 can be carried out to execute a space comparison process to the compensation image to generate a better (modified) compensation image according to the depth, the first brightness enhance value, and the second brightness enhance value. Concretely speaking, as a preferred embodiment, the space comparison process is carried out by using an adaptive 2D gamma filter so as to find out a contrast enhance value γ based on a threshold value rth, the depth d, the first brightness enhance value r1, and the second brightness enhance value r2. The threshold value rth is a predetermined value in brightness. When the depth is greater than the predetermined depth d, the contrast enhance value γ is identified by using the corresponding first brightness value r1 and the second brightness value r2, and also the threshold value rth. Then, the contrast enhance value is used to generate the modified
compensation image 1000 as shown inFIG. 8 to complete the compensation process. - Moreover, it should be noted that, according to a preferred embodiment of the present invention, the contrast enhance value γ, which follows the function: γ=(r1−rth)/r2, indicates the filter parameter along the depth direction. Because the brightness enhancement of certain regions of the compensated image will result in a feeling of reducing contrast, the step S107 adopts the adaptive 2D gamma filter to generate linear contrast enhance value γ to deeper location with brightness enhancement greater than average so as to have the contrast enhancement corresponding with the brightness enhancement, meanwhile, this step also considers the issue that shallow layers should not have too much contrast adjustment so as to generate a better compensation image which is easier to analyze.
- In conclusion, the method of compensating ultrasound image provided in the present invention is capable of compensating brightness along depth direction automatically, thus, the unwanted influence of brightness decline because of medium depth can be effectively reduced to enhance reliability of clinical diagnosis. In addition, the method provided in the present invention is capable to adaptively compensate brightness automatically, which should be an effective approach to resolve the problem of inconvenience and image variation of transitional ultrasound imaging system in which the compensation curve needs to be adjusted manually. Thus, the operation of ultrasound imaging in clinical diagnosis can be more efficiency, which would be quite useful for the industry.
- The detail description of the aforementioned preferred embodiments is for clarifying the feature and the spirit of the present invention. The present invention should not be limited by any of the exemplary embodiments described herein, but should be defined only in accordance with the following claims and their equivalents. Specifically, those skilled in the art should appreciate that they can readily use the disclosed conception and specific embodiments as a basis for designing or modifying other structures for carrying out the same purposes of the present invention without departing from the scope of the invention as defined by the appended claims.
Claims (6)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW102135058 | 2013-09-27 | ||
TW102135058A | 2013-09-27 | ||
TW102135058A TWI485420B (en) | 2013-09-27 | 2013-09-27 | A method of compensating ultrasound image |
Publications (2)
Publication Number | Publication Date |
---|---|
US20150093002A1 true US20150093002A1 (en) | 2015-04-02 |
US9008403B1 US9008403B1 (en) | 2015-04-14 |
Family
ID=52740238
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/164,588 Active US9008403B1 (en) | 2013-09-27 | 2014-01-27 | Method of compensating ultrasound image |
Country Status (2)
Country | Link |
---|---|
US (1) | US9008403B1 (en) |
TW (1) | TWI485420B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108366784A (en) * | 2015-11-30 | 2018-08-03 | 奥林巴斯株式会社 | The working procedure of ultrasound observation apparatus, the working method of ultrasound observation apparatus and ultrasound observation apparatus |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112763573B (en) * | 2020-12-31 | 2024-02-06 | 绵阳同成智能装备股份有限公司 | Banner curve interference compensation method and system based on paper detection scanning device |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5579768A (en) * | 1995-03-21 | 1996-12-03 | Acuson Corporation | Automatic gain compensation in an ultrasound imaging system |
US6159153A (en) * | 1998-12-31 | 2000-12-12 | Duke University | Methods and systems for ultrasound scanning using spatially and spectrally separated transmit ultrasound beams |
US6413218B1 (en) * | 2000-02-10 | 2002-07-02 | Acuson Corporation | Medical diagnostic ultrasound imaging system and method for determining an acoustic output parameter of a transmitted ultrasonic beam |
US6572547B2 (en) * | 2001-07-31 | 2003-06-03 | Koninklijke Philips Electronics N.V. | Transesophageal and transnasal, transesophageal ultrasound imaging systems |
US6592520B1 (en) * | 2001-07-31 | 2003-07-15 | Koninklijke Philips Electronics N.V. | Intravascular ultrasound imaging apparatus and method |
US20080249419A1 (en) * | 2006-12-12 | 2008-10-09 | Sekins K Michael | Time-of-flight triangulation based methods of device spatial registration for multiple-transducer therapeutic ultrasound systems |
US7604594B2 (en) * | 2004-08-04 | 2009-10-20 | General Electric Company | Method and system of controlling ultrasound systems |
US20100305443A1 (en) * | 2007-08-31 | 2010-12-02 | Stewart Gavin Bartlett | Apparatus and method for medical scanning |
US20140088428A1 (en) * | 2012-09-24 | 2014-03-27 | Samsung Medison Co., Ltd. | Ultrasound apparatus and information providing method of the ultrasound apparatus |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6743174B2 (en) * | 2002-04-01 | 2004-06-01 | Koninklijke Philips Electronics N.V. | Ultrasonic diagnostic imaging system with automatically controlled contrast and brightness |
DE602006018229D1 (en) * | 2005-08-24 | 2010-12-30 | Medison Co Ltd | Apparatus and method for processing an ultrasound image |
KR100748858B1 (en) * | 2005-11-24 | 2007-08-13 | 주식회사 메디슨 | Image processing system and method for improving quality of images |
EP2453257A3 (en) * | 2010-11-10 | 2012-08-22 | Samsung Medison Co., Ltd. | Enhancing quality of ultrasound image in ultrasound system |
US8837792B2 (en) * | 2011-06-02 | 2014-09-16 | International Business Machines Corporation | Ultrasound image processing |
-
2013
- 2013-09-27 TW TW102135058A patent/TWI485420B/en active
-
2014
- 2014-01-27 US US14/164,588 patent/US9008403B1/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5579768A (en) * | 1995-03-21 | 1996-12-03 | Acuson Corporation | Automatic gain compensation in an ultrasound imaging system |
US6159153A (en) * | 1998-12-31 | 2000-12-12 | Duke University | Methods and systems for ultrasound scanning using spatially and spectrally separated transmit ultrasound beams |
US6413218B1 (en) * | 2000-02-10 | 2002-07-02 | Acuson Corporation | Medical diagnostic ultrasound imaging system and method for determining an acoustic output parameter of a transmitted ultrasonic beam |
US6572547B2 (en) * | 2001-07-31 | 2003-06-03 | Koninklijke Philips Electronics N.V. | Transesophageal and transnasal, transesophageal ultrasound imaging systems |
US6592520B1 (en) * | 2001-07-31 | 2003-07-15 | Koninklijke Philips Electronics N.V. | Intravascular ultrasound imaging apparatus and method |
US7604594B2 (en) * | 2004-08-04 | 2009-10-20 | General Electric Company | Method and system of controlling ultrasound systems |
US20080249419A1 (en) * | 2006-12-12 | 2008-10-09 | Sekins K Michael | Time-of-flight triangulation based methods of device spatial registration for multiple-transducer therapeutic ultrasound systems |
US20100305443A1 (en) * | 2007-08-31 | 2010-12-02 | Stewart Gavin Bartlett | Apparatus and method for medical scanning |
US20140088428A1 (en) * | 2012-09-24 | 2014-03-27 | Samsung Medison Co., Ltd. | Ultrasound apparatus and information providing method of the ultrasound apparatus |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108366784A (en) * | 2015-11-30 | 2018-08-03 | 奥林巴斯株式会社 | The working procedure of ultrasound observation apparatus, the working method of ultrasound observation apparatus and ultrasound observation apparatus |
Also Published As
Publication number | Publication date |
---|---|
TW201512686A (en) | 2015-04-01 |
US9008403B1 (en) | 2015-04-14 |
TWI485420B (en) | 2015-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Koonsanit et al. | Image enhancement on digital x-ray images using N-CLAHE | |
US10235766B2 (en) | Radiographic image analysis device and method, and storage medium having stored therein program | |
US8374412B2 (en) | Data processing apparatus, medical diagnostic apparatus, data processing method and medical diagnostic method | |
US8433154B2 (en) | Enhanced contrast for scatter compensation in X-ray imaging | |
JP2007144181A (en) | Image processing system and method | |
US10194881B2 (en) | Radiographic image processing device, method, and recording medium | |
US20120108973A1 (en) | Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus | |
US10052080B2 (en) | Medical imaging apparatus and control method thereof | |
US20090171206A1 (en) | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing method | |
WO2012060318A1 (en) | Ultrasonic diagnostic apparatus and ultrasonic image processing apparatus | |
JP2002527834A (en) | Imaging method | |
CN104091309B (en) | Balanced display method and system for flat-plate X-ray image | |
Carpentier et al. | Enhancing your acoustics: ultrasound image optimization of breast lesions | |
WO2019097796A1 (en) | Medical image processing device and medical image processing method | |
US9008403B1 (en) | Method of compensating ultrasound image | |
JP2003153891A (en) | Method for suppressing artifact in computer tomographic raw data, computer tomographic apparatus and computer software product | |
US10646181B2 (en) | Breast type identification device, method, and program | |
Cunitz et al. | Quantification of renal stone contrast with ultrasound in human subjects | |
EP3029634B1 (en) | Echogenicity quantification method and calibration method for ultrasonic device using echogenicity index | |
JP2009268726A (en) | Image display, mammographic apparatus and image display program | |
US8891840B2 (en) | Dynamic steered spatial compounding in ultrasound imaging | |
CN117437144A (en) | Method and system for image denoising | |
US20220319072A1 (en) | Medical image processing apparatus and medical image processing method | |
US10758747B2 (en) | Method and device for determining the stopping power for proton therapy | |
US20210068792A1 (en) | Methods and systems for motion corrected wide-band pulse inversion ultrasonic imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: NATIONAL TAIWAN UNIVERSITY, TAIWAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:LI, PAI-CHI;WEI, YU-MING;SIGNING DATES FROM 20140112 TO 20140114;REEL/FRAME:032126/0202 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2551); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 4 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2552); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 8 |