Patent Grant
9743910
Embodiments described herein relate generally to ultrasonic diagnostic apparatus, an ultrasonic image processing apparatus, a medical image diagnostic apparatus, and a medical image processing apparatus to be used for the purpose of performing appropriate tone correction when generating an image from a signal obtained by ultrasonic scanning using an ultrasonic probe and displaying the image.
An ultrasonic diagnostic apparatus is designed to apply ultrasonic pulses generated from vibration elements provided on an ultrasonic probe into an object and acquire biological information by receiving reflected ultrasonic waves caused by acoustic impedance differences in the tissue of the object through the vibration elements. This apparatus can display image data in real time by simple operation of bringing the ultrasonic probe into contact with the body surface. The apparatus allows to observe a moving object such as a heart and is therefore widely used for morphological diagnosis and functional diagnosis of the circulatory region and various kinds of organs. To improve the viewability of an ultrasonic image acquired by such an ultrasonic diagnostic apparatus, it is necessary to appropriately set and adjust the intensity (tone) of brightness and the like in accordance with the object.
In general, according to one embodiment, the amount of noise and the like unwanted for diagnosis increases depending on a structure or a region (organ) in a living body. This may make it difficult to do adjustment to optimize the contrast of an ultrasonic image or perform appropriate adjustment within the range adjustable by the user himself/herself. However, the conventional ultrasonic diagnostic apparatus has poor operability and, for example, in the above-described cases, cannot appropriately assist in setting/adjusting the tone of an image.
The present invention has been made in consideration of the above-described situations, and has as its object to provide an ultrasonic diagnostic apparatus, an ultrasonic image processing apparatus, a medical image diagnostic apparatus, and a medical image processing apparatus that can appropriately set/adjust the tone of an image as compared to the conventional ultrasonic diagnostic apparatus.
An ultrasonic diagnostic apparatus according to an embodiment, an ultrasonic diagnostic apparatus comprising: an image data acquisition unit configured to transmit an ultrasonic wave into a predetermined region including a diagnosis target of an object, receive a reflected wave from the predetermined region, and acquire ultrasonic image data based on the reflected wave; and a correction unit configured to execute tone correction of the ultrasonic image data, wherein in the tone correction, the correction unit calculates a histogram associated with brightness of the image data, and calculates a brightness distribution range corresponding to the diagnosis target and a tone correction function using the histogram.
According to the present invention, it is possible to implement an ultrasonic diagnostic apparatus, an ultrasonic image processing apparatus, a medical image diagnostic apparatus, and a medical image processing apparatus that can appropriately set/adjust the tone of an image as compared to the conventional ultrasonic diagnostic apparatus.
The embodiments will now be described with reference to the accompanying drawings. Note that the same reference numerals in the following description denote constituent elements having almost the same functions and arrangements, and a repetitive description will be made only when required.
The ultrasonic probe 12 is a device (probe) which transmits ultrasonic waves to an object and receives reflected waves from the object based on the transmitted ultrasonic waves. The ultrasonic probe 12 has, on its distal end, an array of a plurality of piezoelectric transducers, a matching layer, a backing member, and the like. In the ultrasonic probe 12, each of the piezoelectric transducers transmits an ultrasonic wave in a desired direction in a scan region based on a driving signal from the ultrasonic transmission unit 21 and converts a reflected wave from the object into an electrical signal. The matching layer is an intermediate layer which is provided for the piezoelectric transducers to make ultrasonic energy efficiently propagate. The backing member prevents ultrasonic waves from propagating backward from the piezoelectric transducers. When the ultrasonic probe 12 transmits an ultrasonic wave to an object P, the transmitted ultrasonic wave is sequentially reflected by a discontinuity surface of acoustic impedance of internal body tissue, and is received as an echo signal by the ultrasonic probe 12. The amplitude of this echo signal depends on an acoustic impedance difference on the discontinuity surface by which the echo signal is reflected. The echo produced when a transmitted ultrasonic pulse is reflected by the surface of a moving blood flow is subjected to a frequency shift depending on the velocity component of the moving body in the ultrasonic transmission/reception direction due to the Doppler effect.
Note that to acquire volume data, for example, a two-dimensional array probe (a probe having a plurality of ultrasonic transducers arrayed in a two-dimensional matrix) or a mechanical 4D probe (a probe which can execute ultrasonic scanning while mechanically swinging a piezoelectric transducer array in a direction perpendicular to the array direction) may be adopted as the ultrasonic probe 12. However, the ultrasonic probe to be used is not limited to these examples. For example, it is possible to use a one-dimensional array probe as the ultrasonic probe 12 and acquire volume data by performing ultrasonic scanning while manually swinging the probe.
The input device 13 is connected to an apparatus body 11 and includes various types of switches, buttons, a trackball, a mouse, and a keyboard which are used to input, to the apparatus body 11, various types of instructions, conditions, an instruction to set a region of interest (ROI), various types of image quality condition setting instructions, and the like from an operator. The input device 13 also has a start instruction, an instruction to select a target ultrasonic image, and the like for the tone optimization function (to be described later).
The monitor 14 displays morphological information and blood flow information in a living body as images based on video signals from the image processing unit 28.
The ultrasonic transmission unit 21 includes a trigger generation circuit, a delay circuit, and a pulser circuit (none of which are shown). The trigger generation circuit repetitively generates trigger pulses for the formation of transmission ultrasonic waves at a predetermined rate frequency fr Hz (period: 1/fr sec). The delay circuit gives each trigger pulse a delay time necessary to focus an ultrasonic wave into a beam and determine transmission directivity for each channel. The pulser circuit applies a driving pulse to the probe 12 at the timing based on this trigger pulse.
The ultrasonic transmission unit 21 has a function of instantly changing a transmission frequency, a transmission driving voltage, or the like to execute a predetermined scan sequence in accordance with an instruction from the control processor 28. In particular, the function of changing a transmission driving voltage is implemented by a linear amplifier type transmission circuit capable of instantly switching its value or a mechanism of electrically switching a plurality of power supply units.
The ultrasonic reception unit 22 includes an amplifier circuit, an A/D converter, and an adder (none of which are shown). The amplifier circuit amplifies an echo signal received via the probe 12 for each channel. The A/D converter gives each amplified echo signal the delay time required to determine reception directivity and perform reception dynamic focusing. The adder then performs addition processing. This addition processing will enhance a reflection component from a direction corresponding to the reception directivity of the echo signal to form a composite beam for ultrasonic transmission/reception in accordance with the reception directivity and transmission directivity.
The B-mode processing unit 23 receives an echo signal from the ultrasonic reception unit 22, and performs logarithmic amplification, envelope detection processing, and the like for the signal to generate data whose signal intensity is expressed by a brightness level.
The blood flow detection unit 24 extracts a blood flow signal from the echo signal received from the reception unit 22, and generates blood flow data. In general, CFM (Color Flow Mapping) is used for blood flow extraction. In this case, the blood flow detection unit 24 analyzes a blood flow signal to obtain an average velocity, variance, power, and the like as blood flow data at multiple points.
The RAW data memory 25 generates B-mode RAW data and blood flow RAW data for each frame using B-mode data received from the B-mode processing unit 23 and blood flow data received from the blood flow detection unit 24, respectively. The RAW data memory 25 also generates volume data from the RAW data by executing RAW/voxel conversion, as needed.
The image processing unit 28 executes scan conversion processing for the RAW data received from the RAW data memory 25. The image processing unit 28 also performs predetermined image processing such as volume rendering, multi planar reconstruction (MPR), and maximum intensity projection (MIP) for the volume data received from the RAW data memory 25. Note that for the purpose of reducing noise and improving image concatenation, it is possible to perform spatial smoothing by inserting a two-dimensional filter after the image processing unit 28.
The display processing unit 28 executes various kinds of processes associated with a dynamic range, brightness, contrast, γ curve correction, RGB conversion, and the like for various kinds of image data generated/processed by the image processing unit 28. The tone optimization unit 280 included in the display processing unit 28 executes processing corresponding to the tone optimization function (to be described later) under the control of the control processor 29.
The control processor 29 has a function as an information processing apparatus (computer), and controls the operation of this ultrasonic diagnostic apparatus. The control processor 29 reads out a dedicated program for implementing the tone optimization function (to be described later) from the storage unit 30, expands the program on the memory of its own, and executes computation/control and the like associated with various kinds of processes.
The storage unit 30 stores a dedicated program for implementing the tone optimization function (to be described later), diagnosis information (patient ID, findings by doctors, and the like), a diagnostic protocol, transmission/reception conditions, a program for implementing a speckle removal function, a body mark generation program, and other data. The storage unit is also used to store images in the RAW data memory, as needed. It is possible to transfer data in the storage unit 30 to an external peripheral device via the interface unit 31.
The interface unit 31 is an interface associated with the input device 13, a network, and a new external storage device (not shown). The interface unit 31 can transfer data such as ultrasonic images, analysis results, and the like obtained by this apparatus to another apparatus via a network.
The tone optimization function of the ultrasonic diagnostic apparatus 1 will be described next. This function calculates feature points from a histogram generated using individual acquired ultrasonic image data and the first and second deviations of the histogram, and sets a brightness distribution range corresponding to a diagnosis target. The function also sets control conditions using the calculated feature points, and calculates a tone correction function for each image using the set control conditions.
Note that for the sake of a concrete description, tone optimization processing is assumed to be executed for an ultrasonic image acquired by ultrasonic scanning of a two-dimensional region including a cardiac muscle in the B mode. When the ultrasonic image processing apparatus implements the tone optimization processing, for example, ultrasonic image data generated in step S2 is stored in advance, and processing from step S3 is executed. In this case, the ultrasonic image processing apparatus comprises the arrangement indicated by the broken line in
[Patient Information: Reception of Transmission/Reception Conditions as Inputs: Step S1]
The operator inputs patient information and selects transmission/reception conditions (the size of a region to be scanned, a focal position, a transmission voltage, and the like), an imaging mode for ultrasonic scanning on a predetermined region of an object, a scan sequence, and the like via the input device 13 (step S1). The apparatus automatically stores the input and selected various kinds of information and conditions in the storage unit 30.
[Ultrasonic Scanning: Generation of Ultrasonic Image: Steps S2 and S3]
The ultrasonic probe 12 is brought into contact with a desired position on an object surface. Ultrasonic scanning in the B mode is executed for a two-dimensional region including a diagnosis region (the heart, in this case) set as a region to be scanned (step S2). Echo signals acquired by the ultrasonic scanning in the B mode are sequentially sent to the B-mode processing unit 23 via the ultrasonic reception unit 22. The B-mode processing unit 23 executes logarithmic amplification, envelope detection processing, and the like to generate a plurality of B-mode data. The RAW data memory 25 generates B-mode RAW data using the plurality of B-mode data received from the B-mode processing unit 23. The image processing unit 26 executes scan conversion for the generated B-mode RAW data, thereby generating an ultrasonic image for each frame (step S3).
[Tone Optimization Processing: Step S4]
Next, the tone optimization unit 280 executes tone optimization processing as shown in
The tone optimization unit 280 calculates the differential (or difference) of the histogram and then executes smoothing processing, thereby obtaining a result as shown in
Next, the tone optimization unit 280 calculates the feature points of the histogram using the calculated second order differential of the histogram (step S45). In this embodiment, as examples of the feature points of the histogram, four points are employed, that is, the negative minimum (maximal) value of the second order differential of the histogram (corresponding to the positive peak position (maximum value or maximal value) of the histogram), the second zero point of the second order differential of the histogram (the inflection point in the tail region of the histogram: corresponding to the negative peak position of the differential of the histogram), the inflection point in the tail region of the histogram, and the third zero point of the second order differential of the histogram (for example, corresponding to the branch point between the signal from the cardiac muscle region and the signal from another region).
The tone optimization unit 280 sets a brightness distribution range using the calculated feature points (step S46). For example, the tone optimization unit 280 sets a brightness distribution range L by setting a negative minimum value (maximal value) P1 of the second order differential of the histogram as the lower limit and a third zero point P4 of the second order differential of the histogram as the upper limit, as shown in
The tone optimization unit 280 calculates a tone characteristic function using the calculated feature points, set control conditions, and the brightness distribution range L (step S47). For example, the tone optimization unit 280 calculates the tone characteristic function by following control conditions (1) to (5).
(1) The slope of the smoothed cumulative value at the negative minimum value (maximal value) P1 of the second order differential of the histogram is normalized to obtain the slope of the tone characteristic function.
(2) The input minimum value is associated with the output minimum value (a point B1 in the lower graph of
(3) The upper limit of the noise level is set at the inflection point (the second zero point of the second order differential) and associated with a desired output value (a point B2 in the lower graph of
(4) The object region (the cardiac muscle region, in this case) is set to the third zero point from the noise level condition (inflection point), and the output value at the third zero point is associated with a value obtained by a desired slope (a point B3 in the lower graph of
(5) The input maximum value is associated with the output maximum value (a point B4 in the lower graph of
The tone optimization unit 280 calculates, for example, a tone correction function C1 in the lower graph of
[Display of Ultrasonic Image: Step S5]
Next, the monitor 14 displays, in a desired format, the ultrasonic image that has undergone the tone optimization processing (step S5).
Note that the contents of the tone optimization processing are not limited to those described above. Modifications of the tone optimization processing will be explained below.
For example, as shown in
The tone characteristic function may be calculated by adding another control condition to the above-described conditions (1) to (5). For example, it is possible to add the following condition aiming at suppressing the slope of the high brightness region.
(6) An input value that exists from the third zero point to the input maximum value (for example, a tone of 256) is selected and associated with a desired output value (for example, a point B5 in the lower graph of
Note that the control condition (6) may parallelly be added to the control conditions (1) to (5) when calculating the tone characteristic function. The tone characteristic function may newly be calculated by adding the control condition (6) as needed in consideration of the result of the tone characteristic function calculated using the control conditions (1) to (5). If the tone characteristic at the position corresponding to the added control condition is less than, for example, 255, correction is preferably done to make the tone characteristic function monotonically increase so as to suppress the fluctuations in the characteristic.
Tone expansion processing may be performed as needed. Whether to do tone expansion processing is determined in, for example, the following way. If the cumulative sum of histograms of the input image reaches a predetermined amount before the end processing region of smoothing, it is determined that the dynamic range (DR) or tone characteristic at the time of image acquisition is not appropriate. Hence, the tone characteristic is expanded, as shown in
The brightness distribution range can arbitrarily be changed. For example, the above-described brightness distribution range L (lower limit: the negative minimum value of the second order differential of the histogram, upper limit: third zero point) in
The tone optimization processing is also applicable to a case in which a moving image is acquired and displayed.
In the histogram analysis image selection processing (step S40), one of a plurality of acquired ultrasonic images (for example, ultrasonic images of a heart for one or a plurality of heartbeats) is selected, and the tone optimization processing from step S41 is then executed. As the histogram analysis image, for example, an image corresponding to a preset phase or timing, an image corresponding to the start, middle, or end phase of a predetermined period, an image having a maximum average brightness, an image whose average brightness is average in a plurality of images, or the like can be selected. Which image is to be selected as the histogram analysis image can freely be set. If the diagnosis target is the heart, selecting a totally dark image or an image that has not appropriately visualized the object should be prohibited because the cumulative sum of histograms periodically varies in accordance with the dilatation/contraction.
The tone optimization processing can be executed, for example, every time ultrasonic image data is acquired. If the processing need not always be executed for all images, it may be executed at an arbitrary timing in response to, for example, an instruction input by the user via the input device 13.
In the above-described embodiment, an example has been explained in which the first extreme value, the second and third zero points, and the like of the second order differential equation of the histogram associated with the brightness value of an ultrasonic image are calculated as feature points. However, the tone optimization function is not limited to this example. For example, the ith extreme value and the jth zero point of the nth order differential equation (n, i, and j are arbitrary natural numbers) may be calculated as feature points, and the brightness distribution range and the tone correction curve may be calculated using the obtained feature points.
According to the above-described arrangement, the following effects can be obtained.
According to the ultrasonic diagnostic apparatus, feature points are calculated from a histogram generated using acquired ultrasonic image data and the first order differential, the second order differential, and the like of the histogram, and the brightness distribution range corresponding to the diagnosis target is set. Control conditions are set using the calculated feature points, and the tone correction function is calculated using the set control conditions. This allows to define appropriate brightness distribution range and tone correction function corresponding to the tone characteristic of an individual image. It is consequently possible to quickly and easily optimize tone correction in accordance with the diagnosis target, individual differences, and imaging situation.
An ultrasonic diagnostic apparatus according to the second embodiment will be described next. When gain adjustment is performed for an ultrasonic image that is the target of tone optimization processing, the shape of the histogram associated with the brightness value after the gain adjustment changes. As a result, the cumulative sum of histograms and the shape of the nth order differential equation (n is a natural number) also change. In this embodiment, tone optimization processing when the histogram shape has changed due to gain adjustment will be described.
Assume that a predetermined ultrasonic image is acquired, and a histogram associated with the brightness value of the image before gain adjustment is obtained as shown in
When processing of raising the gain is performed for the current ultrasonic image (the ultrasonic image corresponding to the histogram in
When processing of lowering the gain is performed for the current ultrasonic image, the histogram associated with the brightness value and the first order differential equation, the second order differential equation, and the third order differential equation of the histogram after gain lowering are obtained as shown in
That is, the peak position of a histogram, and the polarity of each extreme value and the position of each zero point in the nth order differential equation of the histogram change upon gain adjustment. Hence, the ith extreme value and the jth zero point (i and j are natural numbers) to be used to calculate feature points in tone optimization processing are preferably appropriately selected individually in accordance with the positions and shapes of the histogram and the nth order differential equation after gain correction.
Note that the present invention is not limited to each embodiment described above, and constituent elements can be modified and embodied in the execution stage within the spirit and scope of the invention. The followings are concrete modifications.
(1) Each function associated with each embodiment can also be implemented by installing programs for executing the corresponding processing in a computer such as a workstation and expanding them on a memory. In this case, the programs which can cause the computer to execute the corresponding techniques can be distributed by being stored in recording media such as magnetic disks ((floppy®) disks, hard disks, and the like), optical disks (CD-ROMs, DVDs, and the like), and semiconductor memories.
(2) In the above embodiments, an example has been described in which tone optimization processing is performed for an ultrasonic image acquired by ultrasonic scanning of a two-dimensional region including a cardiac muscle in the B mode. However, the tone optimization function is not limited to this example. The function is also applicable to, for example, a case in which a part other than the heart is the diagnosis target, a case in which imaging is performed not in the B mode but in a blood flow detection mode (CFM mode or the like), a case in which ultrasonic scanning of a three-dimensional region is performed, and the like. Especially in ultrasonic scanning of a three-dimensional region, suitable ultrasonic images can be acquired by applying the tone optimization processing to the two-dimensional images and MPR images included in volume data.
(3) The above embodiments have exemplified the case in which the tone optimization processing is executed for ultrasonic image data acquired by the ultrasonic diagnostic apparatus. However, the tone optimization processing is applicable not only to ultrasonic image data but also to a medical image acquired by a medical image diagnostic apparatus such as an X-ray diagnostic apparatus, an X-ray computed tomography apparatus, or a magnetic resonance imaging apparatus. The tone optimization processing may be performed ex post facto using a medical image processing apparatus for images acquired by various kinds of medical image diagnostic apparatuses.
Various inventions can be formed by proper combinations of a plurality of constituent elements disclosed in the above embodiments. For example, several constituent elements may be omitted from all the constituent elements in each embodiment. In addition, constituent elements of the different embodiments may be combined as needed.
While certain embodiments have been described, these embodiments have been presented by way of example only, and are not intended to limit the scope of the inventions. Indeed, the novel embodiments described herein may be embodied in a variety of other forms; furthermore, various omissions, substitutions and changes in the form of the embodiments described herein may be made without departing from the spirit of the inventions. The accompanying claims and their equivalents are intended to cover such forms or modifications as would fall within the scope and spirit of the inventions.
| Number | Date | Country | Kind |
|---|---|---|---|
| 2010-240112 | Oct 2010 | JP | national |
This application is a Continuation Application of PCT Application No. PCT/JP2011/074684, filed Oct. 26, 2011 and based upon and claiming the benefit of priority from Japanese Patent Application No. 2010-240112, filed Oct. 26, 2010, the entire contents of all of which are incorporated herein by reference.
| Number | Name | Date | Kind |
|---|---|---|---|
| 4662380 | Riley | May 1987 | A |
| 4868651 | Chou et al. | Sep 1989 | A |
| 5469850 | Iizuka et al. | Nov 1995 | A |
| 5579768 | Klesenski | Dec 1996 | A |
| 7586653 | Pulsifer | Sep 2009 | B2 |
| 7676091 | Zwirn | Mar 2010 | B2 |
| 8526729 | Mitsunaga | Sep 2013 | B2 |
| 9390485 | Nakamura | Jul 2016 | B2 |
| 20030144592 | Jeong et al. | Jul 2003 | A1 |
| 20060078182 | Zwirn | Apr 2006 | A1 |
| 20060079777 | Karasawa | Apr 2006 | A1 |
| 20080130964 | Zwirn | Jun 2008 | A1 |
| 20090279764 | Kaji et al. | Nov 2009 | A1 |
| 20100157161 | Tanahashi | Jun 2010 | A1 |
| Number | Date | Country |
|---|---|---|
| 0 220 631 | May 1987 | EP |
| 0 544 328 | Jun 1993 | EP |
| 2 031 873 | Mar 2009 | EP |
| 2002-209891 | Jul 2002 | JP |
| 2005-103007 | Apr 2005 | JP |
| 2005-169155 | Jun 2005 | JP |
| 2007-014525 | Jan 2007 | JP |
| 2007-117351 | May 2007 | JP |
| Entry |
|---|
| Hazewinkel, Michiel, ed. (2001), “Point of inflection”, Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4, http://www.encyclopediaofmath.org/index.php?title=Point—of—inflection&oldid=14389. |
| G Zwirn et al. “A histogram-based technique for echocardiographic image enhancement.” Computers in Cardiology 2004; 31: 81-84. |
| Extended European Search Report issued Jul. 30, 2015 in Patent Application No. 11836344.9. |
| International Preliminary Report on Patentability issued May 16, 2013 in PCT/JP2011/074684 filed Oct. 26, 2011. |
| Written Opinion issued Dec. 20, 2010 in PCT/JP2011/074684 filed Oct. 26, 2011 submitting English language translation only. |
| International Search Report mailed Dec. 20, 2011, in PCT/JP2011/074684 filed Oct. 26, 2011 (with English Translation). |
| International Written Opinion mailed Dec. 20, 2011, in PCT/JP2011/074684 filed Oct. 26, 2011. |
| Number | Date | Country | |
|---|---|---|---|
| 20130237825 A1 | Sep 2013 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | PCT/JP2011/074684 | Oct 2011 | US |
| Child | 13871463 | US |
