WO2007046272A1 - 弾性画像を生成する超音波診断装置 - Google Patents
弾性画像を生成する超音波診断装置 Download PDFInfo
- Publication number
- WO2007046272A1 WO2007046272A1 PCT/JP2006/320221 JP2006320221W WO2007046272A1 WO 2007046272 A1 WO2007046272 A1 WO 2007046272A1 JP 2006320221 W JP2006320221 W JP 2006320221W WO 2007046272 A1 WO2007046272 A1 WO 2007046272A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- histogram
- image
- data
- elasticity
- elastic
- Prior art date
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0048—Detecting, measuring or recording by applying mechanical forces or stimuli
- A61B5/0053—Detecting, measuring or recording by applying mechanical forces or stimuli by applying pressure, e.g. compression, indentation, palpation, grasping, gauging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/42—Details of probe positioning or probe attachment to the patient
- A61B8/4209—Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4411—Device being modular
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
- A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- 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/52079—Constructional features
- G01S7/52084—Constructional features related to particular user interfaces
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/02—Details of sensors specially adapted for in-vivo measurements
- A61B2562/0247—Pressure sensors
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/04—Arrangements of multiple sensors of the same type
- A61B2562/046—Arrangements of multiple sensors of the same type in a matrix array
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61N—ELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
- A61N7/00—Ultrasound therapy
-
- 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/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
-
- 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/52053—Display arrangements
- G01S7/52057—Cathode ray tube displays
- G01S7/52074—Composite displays, e.g. split-screen displays; Combination of multiple images or of images and alphanumeric tabular information
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Definitions
- the present invention relates to other quantitative information related to the elasticity of the tissue in addition to the elasticity image generated by measuring the elasticity (hardness or softness) of the living tissue by applying pressure to the subject.
- the present invention relates to an ultrasonic diagnostic apparatus that can provide accurate diagnosis to a user.
- the ultrasonic transmission / reception surface of an ultrasonic probe is applied to a subject by a technique or a mechanical method. Then, pressure is applied to the subject, the displacement of each part of the living body is obtained using the correlation calculation of the ultrasonic reception signals moving back and forth in time series, and elasticity data relating to the elasticity of the tissue is obtained based on the obtained displacement. It has been proposed to image elasticity data.
- the elasticity data of the tissue includes the displacement of each part of the living body and the strain obtained by spatially differentiating the displacement, as well as the Young's modulus of the living tissue obtained from the stress distribution and strain inside the living body caused by external compression. Representative elastic modulus has been proposed. By measuring such elasticity data and displaying an elasticity image with gradation according to hue and brightness according to the size, elasticity of living tissue can be easily recognized, and cancer tumors It is possible to appropriately distinguish such lesion tissues.
- an elasticity image is an elasticity image with gradation depending on hue and brightness according to the degree of elasticity such as displacement, strain or elastic modulus at each measurement point, the hardness and softness of living tissue. Can be easily recognized
- HIFU high-intensity focused ultrasound
- Patent Document 1 Japanese Patent Laid-Open No. 5-317313
- Patent Document 2 JP 2000-60853 A
- the ultrasonic diagnostic apparatus of the present invention is based on ultrasonic tomographic data measured by applying pressure to the tissue of the subject, and the tomographic site of the subject.
- Elasticity data calculating means for obtaining tissue elasticity data at a plurality of measurement points
- elasticity image generating means for generating an elasticity image of the tomographic site based on the elasticity data
- an elastic distribution based on the elasticity data A histogram evaluation unit that generates a histogram to be represented, and an image display unit that displays at least one of the elasticity image and the image of the histogram are provided.
- a region of interest can be set in the elastic image, and the elastic distribution in the region of interest can be histogrammed.
- the region of interest including the region of interest or the like such as displacement, strain or elastic modulus.
- the histogram is a group of elasticity data such as tissue displacement, strain, or elastic modulus at all measurement points included in the region of interest, and the number of measurement points having the same or similar elasticity data. It is a graph. Therefore, since the examiner can quantitatively and relatively identify the tissue elastic distribution in the region of interest including the region of interest by observing the histogram, it can effectively support the tissue differentiation of the lesion region. . In addition, the therapeutic effect can be accurately determined.
- tissue characteristics such as cancerous tumors that are solid tissues, fibrous glandular species that are soft solid tissues, or cysts that are fluid cystic tissues, and provide clinically useful diagnostic information. it can.
- the region of interest set in the elastic image can be automatically changed following the tissue in the region of interest that is displaced according to the pressure applied to the subject. . According to this, it is possible to correctly recognize how the distribution of the elasticity information of the tissue has changed at an arbitrary time.
- a tomographic image of a tomographic site generated based on ultrasonic tomographic data, an elasticity image, and a histogram can be displayed side by side on the same screen. According to this, since the tissue distribution image appearing in the tomographic image, the elasticity information of each region by the elasticity image, and the auxiliary elasticity information of the histogram can be compared and observed, the accuracy of the tissue discrimination can be further improved.
- the image of the histogram is not limited to the histogram, and a feature amount obtained by deriving the histogram force can be included as auxiliary elasticity information.
- a feature amount obtained by deriving the histogram force can be included as auxiliary elasticity information.
- it includes at least one of the average, standard deviation, and median of elasticity data of the histogram.
- at least one of a function approximating the distribution of the histogram, elasticity data of the peak appearing in the histogram, standard deviation of the peak, and half-width of the peak can be displayed together with the histogram.
- Sarakuko displays a cursor on the histogram image and numerically displays the number of measurement points in one section of the histogram specified by the cursor and the elasticity data in that section. can do.
- multiple cursors are displayed on the histogram image, the number of measurement points in the multiple sections of the histogram sandwiched between the two cursors, the average value of the elasticity data in the multiple sections, and the multiple sections At least one of the standard deviation of the elasticity data and the half width of the peak can be displayed numerically.
- the tomographic image of the tomographic site generated based on the ultrasonic tomographic data, the elastic image, and the histogram are displayed side by side on the same screen, and a plurality of force sols are displayed on the image of the histogram. And at least one of the number of measurement points in the plurality of sections of elasticity data sandwiched between the two cursors, the average value of the elasticity data in the sections, and the standard deviation of the inertial data in the sections.
- the positions of a plurality of measurement points belonging to the section between the two cursors can be displayed in an identifiable manner in the elasticity image.
- a plurality of cursors can be displayed on the histogram image, and the ratio of the elasticity data in the section between the two elasticity data designated by the cursor can be numerically displayed as auxiliary elasticity information. According to this, it is possible to evaluate the elasticity of the tissues at the two sites of interest in a relative ratio.
- a cursor is displayed on the histogram image, the half-value width of the elasticity data of the peak portion of the histogram distribution specified by the cursor is obtained, and the peak portion is supported based on the half-value width!
- the tissue property to be evaluated can be evaluated, and the evaluation result can be displayed as auxiliary elastic information in accordance with the histogram.
- the region of interest set in the elastic image can be automatically changed following the tissue that is displaced according to the pressure, and the histogram in the changed region of interest can be displayed in correspondence.
- the ultrasonic diagnostic apparatus of the present invention is based on the ultrasonic tomographic data measured by applying pressure to the tissue of the subject, and the tissue at a plurality of measurement points of the tomographic site of the subject.
- the displacement vector component is obtained, an elasticity image displaying the displacement vector corresponding to the measurement point of the tomographic site is generated and displayed, and the displacement vector component in a specific direction in the region of interest set in the elasticity image Can be displayed as a histogram.
- the region of interest in addition to an elastic image of a tissue with gradation according to hue and luminance, the region of interest Since the histogram representing the elastic distribution of the tissue in the region is displayed as auxiliary elastic information, the elastic distribution of the tissue in the region of interest including the region of interest can be identified quantitatively and relatively, and tissue differentiation of the lesioned region can be performed. Effective support can be provided.
- ⁇ 1 It is a block configuration diagram of an embodiment of the ultrasonic diagnostic apparatus of the present invention.
- FIG. 2 is a flowchart showing a processing procedure of an embodiment of the ultrasonic diagnostic apparatus of the present invention.
- [3 (A)] is an external view showing an example of an ultrasonic probe.
- [3B] is an external view showing an example of an ultrasonic probe provided with a compression plate.
- [3C] is an external view showing an example of an ultrasonic probe provided with a pressure sensor on the compression plate.
- FIG. 4 is a configuration diagram of a histogram evaluation unit according to an embodiment of the present invention.
- FIG. 5 is a diagram for explaining the concept of elastic frame data.
- FIG. 6 is a diagram illustrating an elastic image according to an embodiment of the present invention.
- FIG. 7 is a diagram for explaining the ROI set in the elastic frame data.
- FIG. 8 is a diagram showing a display example of a histogram according to the first embodiment of the present invention.
- FIG. 9 is a diagram showing an example of an image display mode according to the present invention.
- FIG. 10 is a diagram showing a display example of a histogram according to the fourth embodiment of the present invention.
- FIG. 11 is a diagram showing a display example of a histogram according to the fifth embodiment of the present invention.
- FIG. 12 is a diagram showing a display example of a histogram according to the sixth embodiment of the present invention.
- FIG. 13 is a diagram showing a display example of a histogram according to the seventh embodiment of the present invention.
- FIG. 14 is a diagram showing a display example of a histogram according to the eighth embodiment of the present invention.
- FIG. 15 is a diagram showing a display example of a histogram according to the tenth embodiment of the present invention.
- FIG. 16 is a diagram showing a display example of a histogram for comparison with Example 10.
- FIG. 17 (A)] is a diagram showing a display example of a histogram according to the fourteenth embodiment of the present invention.
- FIG. 17 (B)] is a diagram showing another display example of the histogram of the embodiment 14 of the present invention.
- FIG. 17 (C) is a view showing still another display example of the histogram of the embodiment 14 of the present invention.
- FIG. 18 is a diagram showing a display example of a histogram according to the fifteenth embodiment of the present invention.
- FIG. 19 is a view showing a display example of a histogram according to the embodiment 16 of the present invention.
- FIG. 20 is a view showing another display example of the histogram of the embodiment 16 of the present invention.
- FIG. 21 is a diagram showing a display example of a histogram according to the embodiment 17 of the present invention.
- FIG. 22 is a diagram showing a display example of a histogram according to the example 18 of the present invention.
- FIG. 23 is a diagram showing a display example of a histogram according to the nineteenth embodiment of the present invention.
- FIG. 24 is a view showing a display example of a histogram according to the embodiment 20 of the present invention.
- FIG. 25 is a diagram showing a display example of a histogram according to the embodiment 21 of the present invention.
- FIG. 26 is a diagram showing a display example of a histogram according to the twenty-second embodiment of the present invention.
- FIG. 1 is a block diagram of an ultrasonic diagnostic apparatus according to an embodiment of the present invention.
- FIG. 2 shows a flowchart of an embodiment of the processing procedure of the ultrasonic diagnostic apparatus according to the feature of the present invention.
- the ultrasonic probe 2 used in contact with the subject 1 transmits and receives ultrasonic waves to and from the subject 1 as shown in FIG. 3 (A).
- a plurality of receiving transducers are formed having an ultrasonic transmission / reception surface 21 in which the transducers are aligned.
- the probe 2 is driven by ultrasonic pulses supplied from the transmission circuit 3, and performs beam scanning mechanically or electronically.
- the transmission / reception control circuit 4 controls the transmission timing of the ultrasonic pulses that drive the plurality of transducers of the probe 2 so as to form an ultrasonic beam toward the focal point set in the subject 1. It has become.
- the transmission / reception control circuit 4 electronically scans the ultrasonic beam in the direction in which the transducers of the probe 2 are arranged!
- the probe 2 receives a reflected echo signal generated from within the subject 1 and outputs it to the receiving circuit 5.
- the reception circuit 5 captures the reflected echo signal and performs reception processing such as amplification.
- the reflected echo signal received and processed by the receiving circuit 5 controls the phase of the reflected echo signal received by a plurality of transducers in the phasing and adding circuit 6 and receives ultrasonic waves at one or more convergence points. Form a beam. Reflected echo signal phased and added by phasing and adding circuit 6 This is called ultrasonic tomographic data.
- the radio frequency (RF) signal of the ultrasonic tomographic data generated by the phasing addition circuit 6 may be I and Q signals that have been demodulated in combination.
- the probe 2 scans the ultrasonic beam in the body of the subject 1 in a certain direction, and obtains ultrasonic tomographic data corresponding to one tomographic image!
- the ultrasonic tomographic data processed by the signal processing unit 7 is guided to a black and white scan converter 8, where it is converted into a digital signal and converted into two-dimensional tomographic image data corresponding to the scanning plane of the ultrasonic beam. Converted.
- the black and white scan converter 8 acquires RF signal frame data in the subject 1 including the moving tissue in an ultrasonic cycle, and converts the frame data into an image for display in order to read out in synchronization with the television. It is a means for controlling the tomographic scanning means and the system.
- the A / D converter converts the reflected echo signal from the signal processing unit 7 into a digital signal, and the A / D converter digitizes the A / D converter.
- the signal processing unit 7 and the monochrome scan converter 8 constitute a tomographic image reconstruction means.
- the tomographic image data output from the black-and-white scan converter 8 is supplied to the image display device 10 via the switching adder 9 so that the tomographic image is displayed.
- the image display 10 displays time-series tomographic image data obtained by the black and white scan converter 8, and converts the image data input via the switching adder 9 into an analog signal.
- ⁇ and a color TV monitor that receives the analog video signal from this D / A conversion and displays it as an image.
- the ultrasonic tomographic data output from the phasing addition circuit 6 is guided to the RF signal frame data selection unit 11.
- the RF signal frame data selection unit 11 selects an RF signal group corresponding to the scanning plane (tomographic plane) of the ultrasonic beam as a frame data and stores it in a memory or the like.
- the displacement measurement unit 12 sequentially captures a plurality of pairs of frame data with different acquisition times stored in the RF signal frame data acquisition unit 11, and based on the acquired pair of frame data, displacements of a plurality of measurement points on the tomographic plane.
- the vector is obtained and output to the strain Z modulus calculator 13 as displacement frame data.
- the strain Z elastic modulus calculation unit 13 calculates strains at a plurality of measurement points on the tomographic plane based on the input displacement frame data, and outputs the strain to the elasticity data processing unit 14 as elasticity frame data. Yes. Further, the strain Z elastic modulus calculation unit 13 takes in pressure measurement data applied to the subject from the pressure measurement unit 19, obtains the stress distribution of each part of the subject, and determines the elasticity from the strain frame data and the stress distribution obtained earlier. The rate is obtained and output to the elastic data processing unit 14 as elastic frame data. As shown in FIG. 3 (C), the pressure measurement unit 19 takes in the output signals of a plurality of pressure sensors 23 mounted on the compression plate 22 and measures the pressure applied to the body surface of the subject 1.
- the measured pressure data is sent to the strain Z modulus calculator 13.
- the pressure measurement is based on the displacement frame data of the reference deformation body 24 using the probe 2 with the reference deformation body 24 attached to the front surface of the compression plate 22 as shown in FIG.
- the pressure on the body surface can be measured.
- the elastic data processing unit 14 performs the smoothing process, contrast optimization process in the coordinate plane, and the time axis direction between frames on the elastic frame data of the strain or elastic modulus input from the strain Z elastic modulus calculation unit 13. It performs various image processing such as smoothing processing and sends it to the color scan converter 15! /.
- the color scan converter 15 converts the elastic frame data output from the elastic data processing unit 14 to generate a color elastic image, and displays it on the image display 10 via the switching addition unit 9. ing.
- the color scan converter 15 performs red, green, and gray scales (for example, 256 gradations) on an elastic image based on a range of preset upper and lower limits of elasticity (strain or elastic modulus).
- a hue code such as blue is given.
- a hard region where the elastic modulus of the elastic frame data is measured to be large is converted to a blue code, and conversely, a soft region measured with a small elastic modulus force S is converted to a red code.
- a black and white scan converter can be used in place of the color scan converter 15.
- the distribution of the elastic modulus can be expressed, for example, by increasing the brightness of the hard area measured with a large modulus of elasticity and decreasing the brightness of the soft area measured with a small modulus of elasticity.
- the switching addition unit 9 inputs black and white tomographic image data output from the black and white scan converter 8 and color elastic image data output from the color scan converter 15, and switches both images. It has a function to display one of them, a function to make one of both images translucent, add and combine them and display them on the image display 10, and a function to display both images side by side. Yes. Further, the cine memory unit 18 stores the image data output from the switching addition unit 9 in the memory, and recalls the past image data and displays it on the image display 10 in accordance with a command from the device control interface unit 17. It is like that. In addition, the selected image data can be transferred to a recording medium such as an MO.
- the histogram evaluation unit 20 that is a feature of the present embodiment takes in elastic frame data such as displacement frame data or strain frame data output from the displacement measurement unit 12 and the strain Z elastic modulus calculation unit 13, and will be described later. Thus, a histogram is generated and statistical characteristics are evaluated based on the histogram. Then, histogram image data is generated, and a histogram image is displayed on the image display device 10 via the switching addition unit 9.
- ultrasonic tomographic data is continuously acquired while applying pressure to the tissue of the subject 1 (Sl).
- the ultrasound transmission / reception surface 21 of the probe 2 shown in FIG. 3 (A) is brought into contact with the body surface of the subject 1 and the ultrasound transmission / reception surface 21 is manually moved up and down.
- the compression plate 22 is attached so that the surface is aligned with the ultrasonic wave transmitting / receiving surface 21 of the probe 2.
- the RF signal frame data selection unit 11 sends a frame of the ultrasonic diagnostic apparatus from the phasing addition circuit 6.
- the ultrasonic tomographic data (RF signal frame data) that is output one after another at the frame rate is sequentially secured in the frame memory.
- the latest RF signal frame data is secured as N.
- one RF signal frame data X is selected from the past RF signal frame data N—1, N—2, N—3,.
- the latest RF signal frame data N and the selected RF signal frame data X are output to the displacement measuring unit 12 as a set of RF signal frame data.
- the displacement measurement unit 12 performs one-dimensional or two-dimensional correlation processing on the selected set of RF signal frame data N and X, and displacements or displacement vectors (displacement direction and magnitude) of each measurement point on the tomogram Is measured and displacement frame data is generated.
- a method for detecting the displacement vector there are a block matching method and a gradient method, as described in, for example, Japanese Patent Laid-Open No. 5-317313.
- the block matching method divides the image into multiple blocks with n x n pixel forces, for example, and also searches the previous frame force for the block closest to the block of interest in the current frame to obtain the displacement vector. Is the method.
- the strain Z elastic modulus calculation unit 13 calculates the strain and elastic modulus at each measurement point on the tomographic image based on the displacement frame data and pressure data obtained by the displacement measurement unit 12 and the pressure measurement unit 19. Then, elastic frame data of strain or elastic modulus corresponding to the tomographic image is generated and output to the elastic data processing unit 14.
- the strain is calculated by spatially differentiating the displacement without requiring pressure data.
- the tangential modulus Ym which is one of the elastic moduli, is obtained by calculating the stress (pressure) at each measurement point based on the pressure data as shown in the following formula, and measuring the stress for each measurement. Calculate by dividing by the amount of distortion at the point.
- the indices i and j represent the coordinates of the measurement points of the frame data.
- the elastic data processing unit 14 applies smoothing processing, contrast optimization processing, and inter-frame processing in the coordinate plane to the elastic frame data input from the strain Z elastic modulus calculation unit 13. Performs image processing such as smoothing in the time axis direction and sends it to the color scan converter 15.
- the image processing in the elasticity data processing unit 14 is described in, for example, Japanese Patent Application Laid-Open No. 2004-261198.
- the color scan converter 15 inputs the elastic frame data output from the elastic data processing unit 14 and the upper limit value and the lower limit value that specify the gradation selection range of the elastic frame data, and converts the elastic frame data to the hue. Elastic image data converted into information is generated.
- the upper limit value and the lower limit value are specified by a command from a device control unit (not shown), or a gradation selection range in the elastic frame data as many as 14 elastic data processing units.
- the color scan converter 15 converts an area where the elastic modulus is measured to a blue code, and conversely converts an area where the elastic modulus is measured to a red code.
- a black and white scan converter may be used. In this case, it is possible to brighten the luminance in the region where the distortion is large, and darken the luminance in the region where the distortion is small.
- the elastic image data generated by the color scan converter 15 in this way is output to the switching addition unit 9.
- the switching addition unit 9 adds the black and white tomographic image data output from the black and white scan converter 8 and the color elastic image data output from the color scan converter 15 or outputs them to the image display 10 by switching.
- To display For example, only black and white tomographic image data or color elastic image data is displayed, or both image data are added and combined for display.
- a black and white tomographic image and a color or black and white elastic image can be simultaneously displayed on two screens.
- a black and white tomographic image can be displayed by translucently displaying a color elastic image.
- the display image data output from the switching addition unit 9 is input to the cine memory unit 18 at the same time as being sent to the image display 10.
- the cine memory unit 18 stores input display image data in a memory.
- the cine memory unit 18 calls the past display image data according to the control signal given from the device control interface unit 17.
- the display image data displayed on the image display 10 or selected display data is transferred to a recording medium such as MO.
- histogram evaluation processing according to the feature of the present invention is executed.
- the histogram evaluation unit 20 uses the elastic frame data output from the strain Z elastic modulus calculation unit 13 to determine the strain and elastic modulus at each measurement point in the region of interest at the current time, and obtains a histogram in the region of interest. . Then, image information reflecting the obtained histogram is generated, displayed on the image display 10 via the switching addition unit 9, and provided to the examiner as diagnostic information.
- the image information reflecting this histogram can include the statistical characteristics of the elasticity of each part of the region of interest including the part of interest.
- the histogram evaluation unit 20 includes a memory unit 31, a histogram calculation unit 32, a statistical evaluation unit 33, and an image construction unit 34.
- the memory unit 31 stores the elastic frame data output from the strain Z elastic modulus calculation unit 13 and outputs it to the histogram calculation unit 32.
- the histogram calculation unit 32 receives the elastic frame data output from the memory unit 31 and the region of interest (ROI) coordinate data output from the device control interface unit 17.
- ROI region of interest
- the ROI is input and set from the device control interface unit 17 on the color elastic image displayed on the image display unit 10 in step S5 (S6).
- the histogram calculation unit 32 calculates a histogram using the elastic frame data distributed in the ROI as a population, generates histogram data such as numerical data, and outputs it to the statistical evaluation unit 33. That is, the value of strain or elastic modulus in ROI is divided into a plurality of sections, and the number of measurement points having a strain or elastic modulus value corresponding to each section is calculated. For example, a plurality of sections of strain or elastic modulus are calculated. Histogram data is generated by assigning the number of measurement points on the horizontal axis to the vertical axis.
- the statistical evaluation unit 33 analyzes the histogram data output from the histogram calculation unit 32, and extracts statistical features of the histogram data. Then, statistical evaluation data of statistical features is generated and output to the image construction unit 34 together with the histogram data. Picture The image construction unit 34 constructs an image reflecting the statistical evaluation data and the histogram data output from the statistical evaluation unit 33 as histogram image data, and outputs the histogram image data to the switching addition unit 9.
- Ei, j (t) be the numerical data group of the elastic modulus of the elastic frame data at time t calculated by the strain Z elastic modulus calculator 13 as shown in the conceptual diagram of FIG.
- i l, 2,3, '.', Natural number of N
- j l, 2,3, '.', Natural number of M. That is, in FIG. 5, the index i indicates the coordinate in the horizontal axis direction of the elastic image, and the index j indicates the coordinate in the vertical axis direction of the elastic image.
- the memory unit 31 secures elastic frame data Ei, j (t) and outputs it to the histogram calculation unit 32.
- the region of interest ROI via the input means (for example, a device such as a trackball) of the device control interface unit 17.
- the device control interface unit 17 has coordinates of the set ROI, for example, in the case of a rectangular ROI as shown in FIG. 6, the coordinates of the four corners Pl (xl, yl), P2 (x2, yl), P3 (xl , y2), P4 (x2, y2) are read and output to the histogram calculation unit 32 as ROI coordinate data.
- the histogram calculation unit 32 calculates histogram data using the elastic frame data and the ROI coordinate data. First, referring to the ROI coordinate data, a numerical data group of elastic modulus of measurement points belonging to the inside of the ROI is extracted from the elastic frame data as shown in FIG. 7, for example.
- histogram data is calculated using the numerical data group of the elastic modulus of the measurement points belonging to the inside of the ROI as a population.
- the processing for calculating the histogram data is the same as the procedure for constructing a general histogram. For example, a specific description will be given with reference to FIG. If the histogram display range is elastic modulus, the O (kPa) force is 700 (kPa), and the bin width is 25 (kPa), 28 bins are set. That is, the width of each bin is as follows.
- Bin 1 0—24.99 ⁇ ⁇ ⁇ ⁇ (kPa)
- Bin 2 25-49.99 ⁇ ⁇ ⁇ (kPa)
- Bin 3 50-74.99 ⁇ ⁇ ⁇ (kPa)
- Bin 27 650-674.99
- the medium force of the population also counts the number of measurement points having the elastic modulus in the range of each bin. That is, for example, it is counted as follows
- Bin 1 17 (pieces)
- Bin 3 142 (pieces)
- Bin 27 74 (pieces)
- the relationship between the bin number (bin range) and the number of measurement points is generated as histogram data.
- the statistical evaluation unit 33 calculates statistical characteristics of the generated histogram data, and creates statistical evaluation data.
- This statistical feature includes the number of measurement points that are the samples of the elastic modulus data group as a population, the average value of the elastic modulus data, median value, variance value, standard deviation, half width of the peak, skewness, From the kurtosis etc., it can select and set suitably.
- the statistical characteristics are not limited to these.
- the image construction unit 34 constructs histogram image data based on the histogram data obtained by the histogram calculation unit 32 and the statistical evaluation data obtained by the statistical evaluation unit 33.
- the constructed histogram image data is displayed on the image display device 10 through the switching addition unit 9 and provided to the examiner as auxiliary elasticity information for tissue discrimination.
- the present invention is not limited to this, and the transrectal probe, transesophageal probe are measured.
- any ultrasonic probe such as a probe, an intraoperative probe, and an intravascular probe.
- Example 1 is a specific example of generating a histogram image as shown in FIG.
- ROI includes Organization 1, Organization 2, and Organization 3.
- Information on the hardness inherent to the tissues 1 to 3 is displayed on the image display 10 in real time as an elastic image together with the range display of the tissues 1 to 3.
- a histogram generated based on the elastic modulus data in the ROI is a histogram image as shown in FIG. That is, the vertical axis of the histogram image in the figure correlates with the number of measurement points belonging to the elastic modulus of each bin. Therefore, the larger the number, the longer the bin length.
- the population in the elastic frame data specified by ROI has its own elastic modulus data in organizations 1 to 3, respectively. Therefore, a histogram is formed in which the number of measurement points in the corresponding elastic modulus range is high.
- the statistical evaluation unit 33 for example, the number of measurement points belonging to the population, the average value, the median value, and the variance value of the population's elastic modulus are calculated and adjusted to the histogram as shown in FIG. Displayed on the image display 10 in real time.
- the switching addition unit 9 displays an image on the image display device 10 in various display modes as shown in FIG. 9, for example.
- the tomographic image output from the monochrome scan converter 8, the elastic image output from the color scan converter 15, and the histogram image output from the histogram evaluation unit 20 are combined to construct one display image data. And displayed on the image display 10 in real time.
- the diagnosis depending on the subjectivity of the examiner may give different diagnostic results among the examiners.
- Power it is possible to avoid such an event, establish an objective and universal diagnosis, and provide a clinically useful ultrasonic diagnostic apparatus.
- tissue characteristics such as cancerous tumors that are solid tissues, fibrous gland species that are soft solid tissues, or cysts that are fluid cystic tissues, etc., and provide clinically useful diagnostic information Can do.
- the force showing an example using the elastic modulus is not limited to this, and can be applied to strain data, displacement data, and a displacement vector as inertial data. That is, the strain Z elastic modulus calculation unit 13 obtains the same histogram and statistical characteristic data based on the elastic frame data of the strain calculated in the above manner, and displays the image. I can do it. In addition, based on the displacement frame data calculated by the displacement measuring unit 12, similar histograms and statistical feature data can be obtained to display an image. Also, if the displacement frame data is composed of displacement vectors, similar histograms and statistical feature data can be obtained based on the solid data, and the image can be displayed.
- the histogram and the statistical feature data are obtained based on the elastic frame data composed of the elastic modulus values calculated by the strain Z elastic modulus calculation unit 13 .
- the present invention is not limited to this, and the elastic frame data after the image processing is performed in the elastic data processing unit 14 is input to the histogram evaluation unit 20 to obtain the histogram and statistical feature data. Also, the same effects as in Example 1 can be achieved.
- Example 1 based on the elastic image data converted by the color scan converter 15 Thus, even if the histogram and statistical feature data are obtained, the same effect as in Example 1 can be obtained.
- FIG. 10 shows a display example of the histogram image of the fourth embodiment.
- the histogram is generated and displayed based on the elastic modulus, which is elastic data correlated with the elasticity.
- the present invention is not limited to this, and the statistical evaluation unit 33 approximates (fits) the obtained histogram with an appropriate function such as a double Gaussian function or a plurality of Gaussian functions, as shown in FIG. Images can be displayed. In this case, the average value and standard deviation of each Gaussian function can be obtained and output as statistical evaluation data.
- the fitting function and the elastic modulus ml to m3 at the peak position of the function or the average value M of the elastic modulus and the width information such as the standard deviations ⁇ 1 to ⁇ 3 are obtained.
- Example 2 As compared with Example 1, it is more objective and quantitative evaluation of how hard the region of interest is combined with the surrounding tissue. It is possible to support definitive tissue differentiation.
- FIG. 11 shows a display example of the histogram image of the fifth embodiment.
- the histogram of the arbitrary probability ratio of the histogram shown in the first embodiment is designated and specific histogram information is displayed. That is, a plurality of cursors (Cl, C2,...) Are displayed on the histogram image data via the input means of the device control interface unit 17 so that the examiner can freely move the cursor.
- Information on the histogram of the cursor Cl, C2, ⁇ position specified by the examiner, for example, numerical data such as the number Nl, N2, the average value ml, m2, etc. is displayed in synchronization with the movement of the cursor.
- FIG. 12 shows a display example of the histogram image of the sixth embodiment.
- the range of the histogram is designated with the cursors Cl and C2 ⁇ 'shown in the fifth embodiment, and the statistical characteristics of the range are calculated and displayed as an image.
- FIG. 13 shows a display example of the histogram image of the seventh embodiment.
- region of the elastic image the group of measurement points recorded in the range from the cursor C1 to C2 in the selected histogram V in Example 6 corresponds.
- FIG. 13 it is displayed as a region 41 indicated by a mesh.
- FIG. 14 shows a display example of the histogram image of the eighth embodiment.
- the ratio (m2Zml) of the elastic modulus of the two peak positions CI and C2 of the histogram is calculated and displayed.
- C1 is a normal tissue region, it is possible to immediately grasp how much C2 is harder than the normal tissue.
- the diagnosis depending on the subjectivity of the examiner such as the magnitude of the luminance (color tone) in the elastic image and how it spreads, may result in different diagnosis results among the examiners.
- an objective and universal diagnosis can be established, and a clinically useful ultrasonic diagnostic apparatus can be provided.
- the following functions can be stored in the histogram evaluation unit 20.
- the half-value width of the elastic modulus of the peak portion designated by the cursors Cl, C2, etc. in Fig. 14 is obtained, the tissue properties are evaluated based on the value of the half-value width of the elastic modulus, and the result is output and displayed. That is, for example, when the half width of the elastic modulus is 1 to 20 (kPa), it is evaluated as a cancer tumor which is a very hard solid tissue. When the half width of the elastic modulus is 21 to 70 (kPa), it is evaluated as a fibrous gland type that is a soft solid tissue. If the half-value width of the elastic modulus is 71 (kPa) or more, it is evaluated as a cyst, which is a fluid cystic tissue.
- the ratio of the dispersion value Z average value of the peak portion is obtained to make it dimensionless.
- the ratio of the variance Z average value can be evaluated to be small (for example, 0.1).
- the dispersion is medium and the average is also medium. Therefore, the ratio of the dispersion value Z average value is also evaluated to be medium (for example, 0.3).
- the ratio of the variance Z average is large (for example, 0.5).
- the result of evaluating the tissue property can be output and displayed, so that the reliability of the tissue discrimination can be further improved.
- FIG. 15 shows a display example of the histogram image of the tenth embodiment.
- the ROI set in the elastic image is automatically changed in accordance with the movement of the tissue displaced according to the pressure, and the histogram in the region of interest changed in accordance with the change is displayed. And As a result, ROI is set by tracking the same tissue region, so that appropriate histogram information can be obtained.
- the coordinate data of the ROI is updated in response to the compression operation status of the subject and the biological motion status such as pulsation. That is, in the histogram evaluation unit 20, as a representative point indicating the ROI region set by the device control interface unit 17, for example, as shown in FIG. 6, Pl (tO) at the ROI setting time tO. (xl, yl), P2 (tO) (x2, yl), P3 (t0) (xl, y2), P4 (t0) (x2, y2) force
- the displacement frame data Di, j (tl) at time tl is input from the displacement measuring unit 12 in FIG.
- the position of the tissue displaced at time tl can be obtained. Therefore, the obtained coordinates Pl (tl), P2 (tl), P3 (tl), and P4 (tl) are updated as new ROI coordinate data. Then, using the population distributed within the updated ROI coordinates, the histogram and statistical features are obtained and displayed.
- the force that the tissue in the living body approaches the ultrasonic probe in the compression process by compressing the living body is performed at an arbitrary time.
- the ROI coordinate data is updated according to the compression operation status of the subject and the movement status of the living body such as pulsation.
- the position and size of the ROI changes with time.
- the same tissue region in the living body is extracted at an arbitrary time.
- the population is configured based on the elastic frame data corresponding to the same tissue region, and the histogram and statistical features are calculated and displayed. Therefore, in the compression process, There is no other organization in the ROI. As a result, the number of samples in the population varies. Statistical values such as mean and variance are important information to distinguish tissue characteristics. The features are kept stable at any time, and the values specific to the organization can be evaluated deterministically.
- the ROI coordinate data can be updated according to the compression process, the same tissue region can be extracted at any time, and its statistical characteristics can be evaluated. Accurate diagnosis can be performed with high accuracy.
- the histogram evaluation unit 20 holds the histogram data obtained in the past, and configures the histogram data at the current time and the cumulative histogram data obtained by temporally accumulating the past histogram data for each bin. It is characterized by doing so. In addition, statistical characteristics of accumulated histogram data can be obtained and displayed as images.
- the histogram evaluation unit 20 holds the histogram data obtained in the past, and for each bin, the histogram data at the current time and the smooth histogram data obtained by smoothing the past histogram data with time are processed. Even if you make it up. In this case as well, the statistical characteristics of the smooth histogram data can be obtained and displayed.
- the power described by taking as an example the case of calculating and displaying histograms and statistical features in real time is not limited to this, and the image data stored in the cine memory unit 18 is not limited thereto. Can be used to perform the required processing of FIG.
- the apparatus is frozen through the apparatus control interface unit 17.
- the cine memory unit 18 calls the past display image data according to the control signal of the device control interface unit 17, and designates the display image data to be evaluated by the medium histogram of the display image data of an arbitrary frame and the ROI.
- the ROI coordinate data is output to the histogram evaluation unit 20. It also distorts the information that specifies the frame to be processed.
- the data is output to the rate calculation unit 13 and secured in the strain Z elastic modulus calculation unit 13, and the past elastic frame data is read out and sent to the histogram evaluation unit 20 to calculate the histogram and statistical features.
- this invention can be implemented offline.
- FIG. 17 shows a display example of Example 14 in which a histogram is generated using a displacement vector which is one of elasticity data correlated with elasticity of the present invention.
- tissue displacement vectors at a plurality of measurement points of the tomographic site of the subject are obtained, and the measurement points of the tomographic site are obtained.
- the elasticity image that displays the displacement vector is generated and displayed correspondingly, and the distribution of the displacement vector component in the depth direction, which is the specific direction in the region of interest set in the elasticity image, is displayed as a histogram.
- the horizontal axis of the histogram represents the depth direction component of the displacement.
- FIG. 17 (A) shows the case of a cancer tumor, which is a solid tissue with a very hard region of interest force set in the ROI.
- the displacement vectors at each measurement point point in the same direction as shown by the arrows in the upper diagram of FIG. Therefore, the histogram showing the distribution of the displacement vector in the ROI is a histogram concentrated in a narrow range, as shown in the lower figure of Fig. 17 (A).
- the tissue elements are strongly coupled, so it can be understood that even if subjected to compression, the whole body moves without being deformed.
- FIG. 17 (B) shows a case where the attention area set in the ROI is a fibrous gland type that is a soft solid tissue.
- the displacement vectors at each measurement point in this case are slightly dispersed in direction, as indicated by the arrows in the upper diagram of Fig. 17 (B). Therefore, the histogram showing the distribution of the displacement vector within the ROI has a slightly wider distribution as shown in the lower figure of Fig. 17 (B).
- this example shows some coupling between the tissue elements, It can be understood that they are deformed to some extent while being bound to each other.
- FIG. 17 (C) shows a case where the region of interest set in the ROI is a cyst that is a fluid cystic tissue.
- the displacement vector of each measurement point is largely dispersed in the direction of displacement as indicated by the arrow in the upper diagram of FIG. Therefore, the histogram showing the distribution of the displacement vector within the ROI extends over a wide range as shown in the lower figure of Fig. 17 (C).
- FIG. 18 shows a display example of elasticity information according to Embodiment 15 of the present invention.
- the allocation range of the hue or luminance gradation (for example, 256 gradations) of the elastic image is changed in accordance with the width of the elasticity data of the histogram in the region of interest ROI.
- the gradation allocation range is automatically changed according to the elastic data width of the histogram of the measurement points in the ROI of interest, and the elasticity of the ROI Differences are displayed with a strong contrast.
- FIG. 19 shows a display example of elasticity information according to Embodiment 16 of the present invention.
- a display example of elasticity information in Example 1 shown in FIG. 9 an example in which an elasticity image and an image of auxiliary elasticity information including a histogram (histogram image) are displayed side by side on one screen is shown.
- the histogram image does not need to be overlooked at the same size as the tomographic image or the elastic image because the histogram image is not overlooked even if the resolution is sacrificed.
- the reduced histogram image can be displayed in a small window within the image area of the elastic image. According to this, since the elasticity information of the elasticity image can be observed with sufficient resolution, information loss can be reduced.
- the information of the elastic image of the portion where the histogram image is superimposed is hidden, and there is a possibility that important elastic information may be lost in some cases.
- the histogram image is composed of a translucent image and is displayed superimposed on the elastic image. According to this, it is possible to avoid overlooking the elasticity information due to the superimposed display of the histogram image.
- E (i, j) (Er (i, j), Eg (i, j), Eb (i, j))
- the color information H (i, j) of the pixel at the coordinates (i, j) of the color histogram image is a function of the following equation.
- H (i, j) (Hr (i, j), Hg (i, j), Hb (i, j))
- the color information G (i, j) of the pixel at the coordinates (i, j) of the semi-transparent superposition image of the region where the histogram image is superimposed Is obtained by the following equation.
- the color information of the pixel at the coordinates (i, j) of the semi-transparent superimposed image has both the image information of the elastic image and the histogram image at a ratio of a, and the color information of each image is translucent. Can be recognized at the same time.
- a plurality of images can be equivalently displayed within a limited monitor display screen. If they are arranged so as not to interfere with each other in size, it is possible to avoid the possibility of reducing the accuracy of image diagnosis even when each image becomes smaller.
- FIG. 21 shows a display example of elasticity information according to Embodiment 17 of the present invention.
- 19 and 20 show examples in which a reduced histogram image is displayed in a window in a diagnostic image region such as a tomographic image or an elasticity image.
- a diagnostic image region such as a tomographic image or an elasticity image.
- the size relationship between the display area of the diagnostic image and the histogram image is switched as necessary.
- the size relationship between the display area of the diagnostic image and the histogram image can be configured to switch.
- the image switching key is not limited to the keys and buttons of the device control interface unit 17 and can be displayed on the monitor image and set as an icon.
- FIG. 22 shows a display example of elasticity information when the present invention is applied to a diagnosis for differentiating follicular cancer of the thyroid gland.
- the differentiation of follicular thyroid cancer is very important. It is known that follicular cancer has a distribution of cell density inside the tumor, and the periphery of the tumor has a higher cell density than the center of the tumor.
- the tumor margin tends to be displayed as a hard region compared to the center of the tumor. In other words, it is considered that a region having a high cell density is measured as a region harder than a region having a low cell density.
- FIG. 22 is a display example of an elastic image of a thyroid tumor and a histogram image in the case of follicular cancer. As shown in the figure, to distinguish follicular cancer objectively, tumor center ROI1 and ring tumor By analyzing the histogram of the tumor margin ROI2 as an independent population, information useful for discrimination is acquired and displayed.
- this embodiment is characterized by the following display.
- ROIl and ROI2 There are various methods for setting ROIl and ROI2.
- an input means such as a mouse or a trackball of the device control interface 17 is operated, and the examiner designates the outer periphery and inner periphery of the annular edge portion ROI1 by drawing a line. Thereby, the annular region between the drawing lines can be recognized and set as the edge portion ROI1.
- a plurality of elliptical templates are prepared in advance, and the mouse and the trackball are operated to specify that the two elliptical templates are aligned with the outer periphery and inner periphery of the edge portion ROI1.
- the annular region sandwiched between the outer peripheries of the two elliptical templates can be recognized as the edge portion ROI1, and the inner region of the small elliptical template can be recognized as the RROI2 to be set.
- the annular region sandwiched between the two torus ROIs can be set as the edge portion ROI1 and the inner region of the small torus ROI can be recognized as RROI2.
- a dynamic test is an assessment of how much the tumor of interest is bound to the surrounding tissue. Specifically, when the entire breast is deformed to a large extent, the tumor of interest moves by being dragged by the movement of surrounding tissue, or shows independent movement without being dragged by the movement of surrounding tissue. This is a diagnostic method that evaluates a tomogram while observing it in real time.
- an invasive tumor among breast cancers grows by infiltrating the surrounding tissue, and thus has a high degree of restraint with the surrounding tissue and low mobility.
- benign tumors do not show invasiveness, so they do not show behavior bound to surrounding tissues and are said to be highly mobile.
- assessments of mobility depend on the examiner's experience and objectivity, and are therefore less objective.
- FIGS 23 (a) and 23 (b) show examples of display of elasticity information in this example.
- the part of interest shown in (a) of this figure assumes the case of a highly benign tumor such as a benign tumor.
- a slipping behavior appears at the boundary between the tumor and the mammary gland.
- an ROI is set in the region including the boundary between the tumor and the mammary gland, and the compression operation when a displacement in the depth direction of 0.05 mm, for example, is given within the ROI is evaluated.
- the magnitude of the lateral displacement in the ROI at this time is displayed as a histogram, as shown in the figure, it is the same in the positive direction (for example, rightward in the image) and in the negative direction (for example, leftward in the image).
- Two histogram peaks with a displacement of the order of magnitude appear.
- This histogram can be evaluated as a guide that gives the spread of the histogram (half-width of peak, standard deviation, variance, etc.) and the distance between peaks (mm) 1S.
- the part of interest shown in (b) of the figure is the case of a tumor with low mobility of malignant tumor.
- the tumor is strongly constrained by the surrounding mammary gland, when the whole breast is deformed greatly from the outside, the boundary between the tumor and the mammary gland moves equally.
- the difference between the direction and size of the displacement is small in the region on the tumor side of the boundary and the region on the mammary gland in the boundary, they are completely coupled and displaced.
- the boundary between the tumor and the mammary gland The ROI is set in the area including the boundary, and the compression operation when a depth displacement of 0.05 mm is given in the ROI is used as the evaluation standard, and the magnitude of the lateral displacement in the ROI is a histogram. Is displayed in this figure (this is the case. That is, it has a displacement of the same magnitude in one of the positive directions (for example, rightward in the image) and the negative direction (for example, leftward in the image). Histogram peaks appear, and the spread of the histogram (peak half-width, standard deviation, variance, etc.) can be evaluated as a measure of the degree of mobility.
- Example 19 as one piece of information for benign / malignant discrimination, a histogram based on displacement information is used as auxiliary elasticity information in order to objectively evaluate the mobility related to the degree of invasion to the periphery of the tumor.
- a histogram as auxiliary elastic information for evaluating adhesion between tissues is shown as one information for benign / malignant discrimination.
- Fig. 24 shows an example of an elasticity image when a mammary gland tissue has a tumor.
- this is an example of supplementary elasticity information suitable for performing an operation to remove a tumor near the boundary between the mammary gland and the pectoralis major muscle.
- the tumor may be approached between the mammary gland and the greater pectoral muscle.
- it is very important whether the tumor and the pectoralis major muscle and the mammary gland and the great pectoral muscle are not attached and can be removed.
- an ROI is set in a region where the adhesion between the mammary gland and the great pectoral muscle is to be evaluated.
- the mammary gland and the great pectoral muscle move independently of the deformation of the breast, and a slip occurs at the boundary. Therefore, for example, when a predetermined displacement in the depth direction is given and the lateral displacement is displayed as a histogram image, two histogram peaks appear as shown in FIG.
- the histogram based on the displacement information is used for evaluating the peelable force. It can be used as important information for making a surgical plan in the operation of removing the tumor.
- the evaluation of adhesion is not limited to the case of excision surgery, but also when a bag (such as a saline bag or silicone bag) is inserted between the mammary gland and the greater pectoral muscle during breast augmentation. Can be used effectively to evaluate adhesions between tissues. The effect is similar in external surgery using an endoscope.
- a bag such as a saline bag or silicone bag
- the present embodiment applied the quantitative auxiliary elasticity information of the present invention to the determination of the therapeutic effect of high-intensity focused ultrasound (HIFU).
- HIFU high-intensity focused ultrasound
- the treatment with HIFU is a probe force for ultrasonic therapy.
- Ultrasound is radiated to the treatment site in the subject, high-density ultrasound is focused on a small area of the treatment site, and the focal site is set to 60 to 90, for example.
- Heating to a high temperature of ° C kills the tissue of the lesion (for example, a lesion such as cancer). Then, the focal point is gradually moved within the treatment site, and the entire tissue of the lesion is killed and treated.
- the degree of hardness can be associated with the therapeutic effect.
- the treatment has been applied uniformly throughout the entire tissue of the treatment site. it can.
- the gradation difference can be clearly identified for the minute part of the region where the color gradation approximates. It may not be possible to clearly distinguish whether treatment has been given. In such a case, accurate and objective diagnosis can be supported quantitatively and objectively by the histogram which is the auxiliary elasticity information of the present invention.
- Fig. 25 (a) shows a state in which the cross section of the prostate including the prostate cancer before treatment is depicted by an elastic image. At this time, when the ROI is set in the prostate cancer region, as shown in the figure, for example, it can be seen that it has elasticity centered around 100 (kPa).
- Fig. (B) shows the state of the intermediate process of HIFU treatment in the prostate cancer area.
- Regions of prostate cancer that have been sufficiently irradiated with ultrasound have hardened due to thermal denaturation, for example As shown in the figure (b), it is recognized from the histogram that the hardness has changed to about 300 (kPa). However, as can be seen from the histogram force in (b) of the figure, a tissue region with a hardness of 100 (kPa) representing the hardness of untreated prostate cancer remains inside the ROI, and these regions are It can be seen that this is an incompletely treated region that is not sufficiently irradiated with ultrasound.
- 200 (kPa) or more can also be set as a threshold value as a criterion for determining the completion of treatment.
- a threshold value as a criterion for determining the completion of treatment.
- the method of treating cancer tissue by ablation with heat is not limited to HIFU, and the present invention can also be applied to RFA (percutaneous radiofrequency ablation therapy).
- RFA percutaneous radiofrequency ablation therapy
- cryotherapy in which cancer tissue is rapidly cooled to treat it only involves curing by heat denaturation, and therefore the present invention can be applied in the same manner.
- the treatment treatment is universal for treatment methods whose hardness changes due to tissue degeneration.
- the present invention can be used as a method for evaluating whether or not the treatment has been performed or whether the therapeutic effect has been obtained evenly.
- Example 22 of the present invention suitable for breast cancer screening will be described.
- breast cancer screening is usually performed by scanning an ultrasound probe against the breast and scanning a tomogram of the breast in real time. ing. In such examinations, it is difficult to allocate sufficient examination time to one subject, and an efficient diagnosis is required. However, observing all areas in the measurement cross section of a tomographic image scanned in real time without oversight is a very time-consuming task. Also depends heavily.
- a tomographic image such as a B-mode image is information that is constructed as an image of the structure of a tissue in a living body, and forms a tumor image, the back echo is enhanced, or the marginal area is rough. It uses the form information of the organization such as. In other words, because tumors are detected and differentiated based on many complicated viewpoints, the experience of the examiner is required for efficient diagnosis.
- the tumor is harder than the surrounding fat tissue and mammary gland tissue, for example, when performing an examination based on the elastic modulus as elasticity information, is there a region having a large elastic modulus in the measurement cross section? I will observe it.
- malignant tissue is harder, it is necessary to concentrate attention only on whether the measurement section includes an area with abnormal hardness, so that it can be recognized efficiently even in real time. Is possible.
- the distribution of the elasticity information is histogrammed and displayed as auxiliary elasticity information, thereby realizing an efficient examination.
- elastic information such as strain and elastic modulus is measured in real time during the cross-sectional scan, and the distribution of the elastic information is displayed as a histogram.
- the threshold Eth of suspected malignancy is assumed to be 200 (kPa)
- the distribution of measurement points of elastic modulus higher than Eth existing in the current elasticity image can be immediately judged by just looking at the histogram distribution. Can do. In other words, every corner of an elastic image is confirmed, but according to the present embodiment, such a burden is eliminated.
- the scale of the vertical axis (frequency) of the histogram is changed from the linear scale to the log scale. You may be able to change.
- a log scale even if a hard tumor that is suspected of being malignant is very small or low in frequency, it can be detected with high sensitivity without overlooking it.
- the auxiliary elasticity information including the histogram of the elasticity distribution is added to the elasticity image based on the displacement, strain, elastic modulus, or the like, or instead of the elasticity image, as described in each embodiment.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Human Computer Interaction (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
被検体1の組織に圧力を加えて計測された超音波断層データに基づいて、被検体の断層部位の複数の計測点における組織の歪みや弾性率などの弾性データを求め、その弾性データに基づいて断層部位の弾性画像を生成して表示し、その弾性画像に設定される関心領域ROIにおける歪みや弾性率などの弾性データの分布をヒストグラムで表示することにより、色相や輝度により階調を付けた組織の弾性に関する情報以外に、組織性状を反映した新たな定量的な情報を提供して組織の鑑別を支援する。
Description
明 細 書
超音波診断装置
技術分野
[0001] 本発明は、被検体に圧迫を付与して生体組織の弾性 (硬さ又は軟らかさ)を計測し て生成される弾性画像に加え、組織の弾性に関係する定量的な他の情報をユーザ に提供して、的確な診断を支援することができる超音波診断装置に関する。
背景技術
[0002] 超音波診断装置の技術分野では、特許文献 1、 2等に記載されているように、被検 体に用手法又は機械的な方法により、超音波探触子の超音波送受信面にて被検体 に圧力を加え、時系列的に前後する超音波受信信号の相関演算を利用して生体各 部の変位を求め、求めた変位に基づいて組織の弾性に関する弾性データを求め、そ の弾性データを画像ィ匕することが提案されている。また、組織の弾性データとしては 、生体各部の変位及び変位を空間微分した歪み、さらには、外部からの圧迫により生 じた生体内部の応力分布と歪みとから得られる生体組織のヤング率等に代表される 弾性率などが提案されている。このような弾性データを計測して、それらの大きさに応 じて色相や輝度により階調を付けた弾性画像を表示することにより、生体組織の弾性 を容易に認識することができ、ガン腫瘍等の病変組織を適切に鑑別することが可能 になる。
[0003] つまり、弾性画像は、各計測点における変位、歪み又は弾性率などの弾性の程度 に応じて、色相や輝度により階調を付けた弾性画像であるから、生体組織の硬さや 軟らかさを容易に認識することができる
しかし、弾性画像に基づいて注目部位の弾性の違いを判別して、病変部を鑑別し ようとする場合、異なる色階調の境界領域、あるいは、色階調が近似する領域につい て、弾性の差異を明確に区別できない場合があり、検者の経験や熟練度によって鑑 別に誤認が生ずるおそれがある。
[0004] また、高密度焦点式超音波(以下、 HIFU: high-intensity focused ultrasound と いう。)を病変部に照射して、病変部の組織を破壊して治療する方法が提案されてい
る。これによれば、破壊により死滅された病変部の組織は熱変性により硬くなるから、 その硬さの程度と治療効果とを対応付けることができる。そこで、治療部位内の弾性 を計測して弾性画像を観察すれば、万遍なく均一に治療が施された力どうかを判定 できるが、上述したように、弾性画像では色階調が近似する領域の微小部位につい て階調差を明確に識別できず、治療が施されたかどうかを明確に判断できない場合 がある。
[0005] そこで、弾性画像に基づ!/ヽて病変部を鑑別する場合、あるいは治療の効果を判断 する場合に、色階調等による弾性情報だけでなぐ組織の弾性に関係する定量的な 他の情報を合わせて考慮することができれば、的確な診断の支援をすることが期待 できる。
[0006] 特許文献 1 :特開平 5— 317313号公報
特許文献 2 :特開 2000— 60853号公報
発明の開示
[0007] 本発明は、色相や輝度により階調を付けた組織の弾性画像以外に、的確な診断を 支援する定量的な補助弾性情報を提供することを課題とする。
[0008] 上記の課題を解決する手段を検討した結果、第一に、圧迫などの外力による組織 変形は、組織を構成する組織要素群によって固有の集団的挙動を示すことに着目し 、その固有の集団的挙動を統計的な特徴として把握すれば、的確な診断を支援する 補助弾性情報となり得ることに思 、至った。
[0009] そこで、上記の課題を解決するため、本発明の超音波診断装置は、被検体の組織 に圧力を加えて計測された超音波断層データに基づ 、て、前記被検体の断層部位 の複数の計測点における組織の弾性データを求める弾性データ演算手段と、該弹 性データに基づいて前記断層部位の弾性画像を生成する弾性画像生成手段と、前 記弾性データに基いて弾性分布を表すヒストグラムを生成するヒストグラム評価手段と 、前記弾性画像と前記ヒストグラムの画像の少なくとも一方を表示する画像表示手段 を備えてなることを特徴とする。この場合、弾性画像に関心領域を設定し、その関心 領域における弾性分布をヒストグラム化することができる。
[0010] すなわち、注目部位を含む関心領域等にっ 、て、変位、歪み又は弾性率などの弹
性データの分布をヒストグラム化して表示することにより、組織性状を反映した集団的 挙動に基づいた定量的な補助弾性情報として提供することができる。ここで、ヒストグ ラムは、関心領域に含まれる全ての各計測点における組織の変位、歪み、あるいは 弾性率などの弾性データを母集団とし、同一又は同程度の弾性データを有する計測 点の個数をグラフ化したものである。したがって、検者は、注目部位を含む関心領域 における組織の弾性分布を、ヒストグラムを観察することにより、定量的かつ相対的に 識別できるから、病変部位の組織鑑別を効果的に支援することができる。また、治療 効果の判断を的確に行うことができる。
[0011] 特に、ヒストグラムには、組織の集団的挙動の特徴が現れることから、これにより、病 変部などの注目部位の組織の弾性の均一性などを評価することができ、非常に硬 、 充実性組織であるガン腫瘍、軟らかい充実性組織である繊維腺種、あるいは流動的 な嚢胞性組織である嚢胞、等々の組織性状の識別が可能となり、臨床上有用な診断 情報を提供することができる。
[0012] 上記の場合にぉ ヽて、弾性画像に設定された関心領域を、被検体に加えられる圧 力に応じて変位する関心領域内の組織に追従させて自動的に変更することができる 。これによれば、任意の時刻において、組織の弾性情報の分布がどのように変化した かを、正しく認識することができる。
[0013] また、超音波断層データに基づいて生成された断層部位の断層画像と、弾性画像 と、ヒストグラムとを同一の画面に並べて表示することができる。これによれば、断層画 像に表れた組織分布像と、弾性画像による個々の領域の弾性の情報と、ヒストグラム の補助弾性情報を対比観察できるから、一層、組織鑑別の精度を向上できる。
[0014] ここで、ヒストグラムの画像には、ヒストグラムだけに限られず、ヒストグラム力 派生し て求められる特徴量を補助弾性情報として含めることができる。例えば、ヒストグラム の弾性データの平均、標準偏差、中央値の少なくとも一つを含む。また、ヒストグラム の分布を近似した関数と、ヒストグラムに表れたピークの弾性データと、ピークの標準 偏差と、ピークの半値幅の少なくとも一つをヒストグラムと共に表示することができる。 さら〖こ、ヒストグラムの画像上にカーソルを表示し、そのカーソルにより指定されるヒスト グラムの 1つの区間における計測点の個数と、その区間の弾性データを数値で表示
することができる。また、ヒストグラムの画像上に複数のカーソルを表示し、 2つのカー ソルに挟まれるヒストグラムの複数の区間における計測点の個数と、その複数の区間 の弾性データの平均値と、その複数の区間の弾性データの標準偏差と、ピークの半 値幅の少なくとも一つを数値で表示することができる。
[0015] さらに、超音波断層データに基づいて生成された断層部位の断層画像と、弾性画 像と、ヒストグラムとを同一の画面に並べて表示し、ヒストグラムの画像上に複数の力 一ソルを表示し、 2つのカーソルに挟まれる弾性データの複数の区間における計測 点の個数と、その複数の区間の弾性データの平均値と、その複数の区間の弹性デ一 タの標準偏差の少なくとも一つを、補助弾性情報として数値で表示すると共に、 2つ のカーソルに挟まれた区間に属する複数の計測点の位置を弾性画像に識別可能に 表示することができる。
[0016] また、ヒストグラムの画像上に複数のカーソルを表示し、該カーソルにより指定される 2つの弾性データの区間における弾性データの比を、補助弾性情報として数値で表 示することができる。これによれば、注目する 2つの部位の組織の弾性を相対比で評 価することができる。
[0017] また、ヒストグラムの画像上にカーソルを表示し、該カーソルにより指定されるヒストグ ラム分布のピーク部の弾性データの半値幅を求め、その半値幅に基づ!/、てピーク部 に対応する組織性状を評価し、評価結果を前記ヒストグラムに合わせて補助弾性情 報として表示することができる。
[0018] さらに、圧力に応じて変位する組織に追従させて弾性画像に設定された関心領域 を自動的に変更すると共に、変更された関心領域におけるヒストグラムを対応させて 表示することができる。
[0019] また、本発明の超音波診断装置は、被検体の組織に圧力を加えて計測された超音 波断層データに基づ 、て、前記被検体の断層部位の複数の計測点における組織の 変位ベクトルを求め、前記断層部位の前記計測点に対応させて前記変位ベクトルを 表示する弾性画像を生成して表示し、該弾性画像に設定される関心領域における特 定方向の前記変位ベクトル成分の分布をヒストグラムで表示することができる。
[0020] 本発明によれば、色相や輝度により階調を付けた組織の弾性画像以外に、関心領
域における組織の弾性分布を表わすヒストグラムを補助弾性情報として表示するよう にしているから、注目部位を含む関心領域における組織の弾性分布を、定量的かつ 相対的に識別でき、病変部位の組織鑑別を効果的に支援することができる。
図面の簡単な説明
圆 1]本発明の超音波診断装置の一実施形態のブロック構成図である。
[図 2]本発明の超音波診断装置の一実施形態の処理手順を示すフローチャートであ る。
圆 3(A)]超音波の探触子の一例を示す外観図である。
圆 3(B)]圧迫板を備えた超音波の探触子の一例を示す外観図である。
圆 3(C)]圧迫板に圧力センサーを備えた超音波の探触子の一例を示す外観図であ る。
圆 3(D)]圧迫板に参照変形体を備えた超音波の探触子の一例を示す外観図である 図 4]本発明の一実施形態のヒストグラム評価部の構成図である。
図 5]弾性フレームデータの概念を説明する図である。
図 6]本発明の一実施形態の弾性画像を説明する図である。
図 7]弾性フレームデータに設定された ROIを説明する図である。
図 8]本発明の実施例 1のヒストグラムの表示例を示す図である。
図 9]本発明に係る画像の表示形態の一実施例を示す図である。
図 10]本発明の実施例 4のヒストグラムの表示例を示す図である。
図 11]本発明の実施例 5のヒストグラムの表示例を示す図である。
図 12]本発明の実施例 6のヒストグラムの表示例を示す図である。
図 13]本発明の実施例 7のヒストグラムの表示例を示す図である。
図 14]本発明の実施例 8のヒストグラムの表示例を示す図である。
図 15]本発明の実施例 10のヒストグラムの表示例を示す図である。
図 16]実施例 10と比較のためのヒストグラムの表示例を示す図である。
図 17(A)]本発明の実施例 14のヒストグラムの一表示例を示す図である。
図 17(B)]本発明の実施例 14のヒストグラムの他の表示例を示す図である
[図 17(C)]本発明の実施例 14のヒストグラムのさらに他の表示例を示す図である。
[図 18]本発明の実施例 15のヒストグラムの表示例を示す図である。
[図 19]本発明の実施例 16のヒストグラムの表示例を示す図である。
[図 20]本発明の実施例 16のヒストグラムの他の表示例を示す図である。
[図 21]本発明の実施例 17のヒストグラムの表示例を示す図である。
[図 22]本発明の実施例 18のヒストグラムの表示例を示す図である。
[図 23]本発明の実施例 19のヒストグラムの表示例を示す図である。
[図 24]本発明の実施例 20のヒストグラムの表示例を示す図である。
[図 25]本発明の実施例 21のヒストグラムの表示例を示す図である。
[図 26]本発明の実施例 22のヒストグラムの表示例を示す図である。
発明を実施するための最良の形態
[0022] 以下、本発明を実施形態に基づいて説明する。図 1に、本発明の一実施形態の超 音波診断装置のブロック構成図を示す。図 2に、本発明の特徴部に係る超音波診断 装置の処理手順の一実施形態のフローチャートを示す。
[0023] 図 1に示すように、被検体 1に当接して用いられる超音波の探触子 2は、図 3 (A)に 示すように、被検体 1との間で超音波を送信及び受信する複数の振動子が整列され た超音波送受信面 21を有して形成されている。探触子 2は、送信回路 3から供給さ れる超音波パルスにより駆動され、機械式又は電子的にビーム走査を行うようになつ ている。送受信制御回路 4は、探触子 2の複数の振動子を駆動する超音波パルスの 送信タイミングを制御して、被検体 1内に設定される焦点に向けて超音波ビームを形 成するようになっている。また、送受信制御回路 4は、探触子 2の振動子の配列方向 に電子的に超音波ビームを走査するようになって!/、る。
[0024] 一方、探触子 2は、被検体 1内から発生する反射エコー信号を受信して受信回路 5 に出力する。受信回路 5は、送受信制御回路 4から入力されるタイミング信号に従つ て、反射エコー信号を取り込んで増幅などの受信処理を行う。受信回路 5により受信 処理された反射エコー信号は、整相加算回路 6において複数の振動子により受信さ れた反射エコー信号の位相を制御し、一点又は複数の収束点に対して超音波受波 ビームを形成する。整相加算回路 6において整相加算された反射エコー信号 (以下
、超音波断層データという。 )は、信号処理部 7に入力され、ゲイン補正、ログ圧縮、 検波、輪郭強調、フィルタ処理等の信号処理がなされる。なお、整相加算回路 6にお いて生成される超音波断層データの高周波 (RF)信号は、複合復調した I、 Q信号で あっても良い。このようにして、探触子 2で超音波ビームを被検体 1の体内で一定方 向に走査し、一枚の断層像に対応した超音波断層データを得るようになって!/、る。
[0025] 信号処理部 7により処理された超音波断層データは白黒スキャンコンバータ 8に導 かれ、ここにおいてディジタル信号に変換されるとともに、超音波ビームの走査面に 対応した 2次元の断層像データに変換される。すなわち、白黒スキャンコンバータ 8 は、運動組織を含む被検体 1内の RF信号フレームデータを超音波周期で取得し、そ のフレームデータを画像に変換して表示するために、テレビ同期で読み出すための 断層走査手段及びシステムの制御を行うための手段となるもので、信号処理部 7から の反射エコー信号をディジタル信号に変換する A/D変換器と、この A/D変換器でデ イジタルィ匕された複数枚の断層像データを時系列に記憶するフレームメモリと、これら の動作を制御するコントローラなど力 構成されている。これらの信号処理部 7と白黒 スキャンコンバータ 8によって断層像の画像再構成手段が構成される。 白黒スキャン コンバータ 8から出力される断層像データは、切替加算部 9を介して画像表示器 10 に供給されて断層像が表示されるようになって ヽる。
[0026] 画像表示器 10は、白黒スキャンコンバータ 8によって得た時系列の断層像データを 表示するもので、切替加算部 9を介して入力される画像データをアナログ信号に変換 する D/A変^^と、この D/A変 からのアナログビデオ信号を入力して画像とし て表示するカラーテレビモニタとからなる。
[0027] 一方、整相加算回路 6から出力される超音波断層データは、 RF信号フレームデー タ選択部 11に導かれる。 RF信号フレームデータ選択部 11は、超音波ビームの走査 面(断層面)に対応する RF信号群を、フレームデータとして複数フレーム分を選択し てメモリなどに格納する。変位計測部 12は、 RF信号フレームデータ取得部 11に格 納されている取得時刻が異なる複数対のフレームデータを順次取り込み、取り込ん だ一対のフレームデータに基づいて断層面における複数の計測点の変位ベクトルを 求め、変位フレームデータとして歪み Z弾性率演算部 13に出力するようになって ヽ
る。
[0028] 歪み Z弾性率演算部 13は、入力される変位フレームデータに基づいて断層面に おける複数の計測点の歪みを求め、弾性フレームデータとして弾性データ処理部 14 に出力するようになっている。また、歪み Z弾性率演算部 13は、圧力計測部 19から 被検体に加えられた圧力計測データを取り込み、被検体各部の応力分布を求め、先 に求めた歪みフレームデータと応力分布とから弾性率を求め、弾性フレームデータと して、弾性データ処理部 14に出力するようになっている。圧力計測部 19は、図 3 (C) に示すように、圧迫板 22などに装着された複数の圧力センサー 23の出力信号を取り 込み、被検体 1の体表面に加えられた圧力を計測し、計測した圧力データを歪み Z 弾性率演算部 13に送出する。なお、圧力計測は、図 3 (D)に示すように、圧迫板 22 の前面に装着された参照変形体 24を装着した探触子 2を用い、参照変形体 24の変 位フレームデータに基づいて体表面にカ卩えられた圧力を計測するようにすることがで きる。これらの圧力計測法は、特開 2005— 13283号公報又は特開 2005— 66041 号公報に記載されている。
[0029] 弾性データ処理部 14は、歪み Z弾性率演算部 13から入力される歪み又は弾性率 の弾性フレームデータに座標平面内におけるスムージング処理、コントラスト最適化 処理や、フレーム間における時間軸方向のスムージング処理などの様々な画像処理 を施して、カラースキャンコンバータ 15に送出するようになって!/、る。
[0030] カラースキャンコンバータ 15は、弾性データ処理部 14から出力される弾性フレーム データを変換してカラーの弾性画像を生成し、切替加算部 9を介して画像表示器 10 に表示させるようになつている。つまり、カラースキャンコンバータ 15は、予め設定さ れた弾性 (歪み又は弾性率)の上限値及び下限値の範囲に基づいて、弾性画像に 階調化 (例えば、 256階調)された赤、緑、青などの色相コードを付与する。例えば、 弾性フレームデータの弾性率が大きく計測された硬い領域は青色コードに変換し、 逆に弾性率力 S小さく計測された軟らかい領域は赤色コードに変換する。なお、カラー スキャンコンバータ 15に代えて、白黒スキャンコンバータを用いることができる。この 場合は、弾性率が大きく計測された硬い領域は輝度を明るぐ逆に弾性率が小さく計 測された軟らかい領域は輝度を暗くするなどにより、弾性率の分布を表すことができ
る。
[0031] また、切替加算部 9は、白黒スキャンコンバータ 8から出力される白黒の断層像デー タと、カラースキャンコンバータ 15から出力されるカラーの弾性画像データとを入力し 、両画像を切り替えていずれか一方を表示させる機能と、両画像の一方を半透明に して加算合成して画像表示器 10に重ねて表示させる機能と、両画像を並べて表示さ せる機能を有して形成されている。また、シネメモリ部 18は、切替加算部 9から出力さ れる画像データをメモリに格納し、装置制御インターフェイス部 17からの指令に従つ て、過去の画像データを呼び出して画像表示器 10に表示するようになっている。さら に、選択された画像データを MOなどの記録メディアへ転送することが可能になって いる。
[0032] さらに、本実施形態の特徴であるヒストグラム評価部 20は、変位計測部 12と歪み Z 弾性率演算部 13から出力される変位フレームデータあるいは歪みフレームデータ等 の弾性フレームデータを取り込み、後述するように、ヒストグラムを生成すると共に、ヒ ストグラムに基づいて統計的特徴を評価するようになっている。そして、ヒストグラム画 像データを生成して、切替加算部 9を介して画像表示器 10にヒストグラム画像を表示 させるようになつている。
[0033] このように構成される実施形態の詳細構成を動作とともに、図 2に示したフローチヤ ートを参照して説明する。
(ステップ S1)
まず、被検体 1の組織に圧力を加えながら、連続的に超音波断層データを取得す る(Sl)。このとき、図 3 (A)に示した探触子 2の超音波送受信面 21を被検体 1の体表 に接触させて超音波送受信を行いつつ、超音波送受信面 21を用手的に上下動させ て被検体 1を圧迫する。このとき、被検体 1の診断部位に効果的に応力分布を与える ため、図 3 (B)に示すように、探触子 2の超音波送受信面 21に面を合わせて圧迫板 2 2を装着し、超音波送受信面 21と圧迫板 22により構成される圧迫面で被検体 1を圧 迫することが好ましい。
(ステップ S 2)
RF信号フレームデータ選択部 11は、整相加算回路 6から超音波診断装置のフレ
ームレートで経時的に次々と出力される超音波断層データ (RF信号フレームデータ) をフレームメモリ内に順次確保する。このとき、最新の RF信号フレームデータを Nとし て確保する。そして、制御命令に従って時間的に過去の RF信号フレームデータ N— 1、 N— 2、 N— 3、 · · ·、 N— Mの中から 1つの RF信号フレームデータ Xを選択する。次 Vヽで、最新の RF信号フレームデータ Nと選択した RF信号フレームデータ Xを 1組の R F信号フレームデータとして、変位計測部 12に出力する。
(ステップ S3)
変位計測部 12は、選択された 1組の RF信号フレームデータ N, Xを 1次元もしくは 2 次元相関処理して、断層像上の各計測点の変位もしくは変位ベクトル (変位の方向と 大きさ)を計測し、変位フレームデータを生成する。この変位ベクトルの検出法として は、例えば特開平 5— 317313号公報に記載されているように、ブロックマッチング法 とグラジェント法とがある。ブロックマッチング法は、画像を例えば n X n画素力もなる 複数のブロックに分け、現フレーム中の着目しているブロックに最も近似しているブロ ックを前フレーム力も探索して、変位ベクトルを求める方法である。
(ステップ S4)
歪み Z弾性率演算部 13は、変位計測部 12及び圧力計測部 19で求めた変位フレ ームデータ及び圧力データに基づ 、て、断層像上の各計測点の歪み及び弾性率を 演算する。そして、断層像に対応させた歪み又は弾性率の弾性フレームデータを生 成して、弾性データ処理部 14に出力する。ここで、歪みは、圧力データを必要とせず 、変位を空間微分することによって計算される。また、弾性率の一つである、例えばャ ング率 Ymの演算は、以下の式に示したように、圧力データに基づいて各計測点に おける応力(圧力)を求め、その応力を各計測点における歪み量で除することにより 求める。ここで、 i,jの指標は、フレームデータの計測点の座標を表す。
Ymi,j =圧力(応力) ijZ (歪み量 i,j)
ここで、 (i,j=l, 2, 3, - · -)
(ステップ S 5)
弾性データ処理部 14は、歪み Z弾性率演算部 13から入力される弾性フレームデ ータに、座標平面内におけるスムージング処理、コントラスト最適化処理、フレーム間
における時間軸方向のスムージング処理などの画像処理を施して、カラースキャンコ ンバータ 15に送出する。この弾性データ処理部 14における画像処理については、 例えば、特開 2004— 261198号公報に記載されている。
[0035] カラースキャンコンバータ 15は、弾性データ処理部 14から出力される弾性フレーム データと、弾性フレームデータの階調化選択範囲を指定する上限値及び下限値を入 力し、弾性フレームデータを色相情報に変換した弾性画像データを生成する。この 上限値及び下限値は、図示していない装置制御部からの命令、又は弾性データ処 理部 14力もの弾性フレームデータの中の階調化選択範囲によって指定される。例え ば、カラースキャンコンバータ 15は、入力される弾性フレームデータにおいて、弾性 率が大きく計測された領域は青色コードに変換し、逆に弾性率が小さく計測された領 域は赤色コードに変換する。
[0036] ここで、カラースキャンコンバータ 15に代えて、白黒スキャンコンバータを用いること もできる。この場合は、歪みが大きい領域は輝度を明るぐ逆に歪みが小さい領域は 輝度を暗くさせるようにすることができる。
[0037] このようにしてカラースキャンコンバータ 15で生成された弾性画像データは、切替 加算部 9に出力される。切替加算部 9は、白黒スキャンコンバータ 8から出力される白 黒の断層像データと、カラースキャンコンバータ 15から出力されるカラーの弾性画像 データとを加算して、又は切り替えて画像表示器 10に出力して表示させる。例えば、 白黒の断層像データだけ、又はカラーの弾性画像データだけを表示させ、あるいは 両画像データを加算合成して表示させるようになつている。また、例えば、特開 2000 — 60853号公報に記載されているように、白黒断層像とカラー又は白黒の弾性画像 を同時に 2画面表示するようにすることができる。さらに、例えば、特開 2004— 1359 29号公報に記載されているように、白黒断層像にカラーの弾性画像を半透明的に重 畳して表示するようにすることができる。
[0038] また、切替加算部 9から出力される表示画像データは、画像表示器 10に送られると 同時に、シネメモリ部 18に入力されるようになっている。シネメモリ部 18は、入力され る表示画像データをメモリに格納する。また、シネメモリ部 18は、装置制御インターフ エイス部 17から与えられる制御信号に従って、過去の表示画像データを呼び出して
画像表示器 10に表示し、又は選択された表示画像データを MOなどの記録メディア へ転送するようになって!/、る。
(ステップ S6、 S7)
ステップ S6、 S7において、本発明の特徴に係るヒストグラム評価処理を実行する。 ヒストグラム評価部 20は、歪み Z弾性率演算部 13から出力される弾性フレームデー タを用い、現時刻における関心領域内の各計測点における歪みや弾性率を求めて、 関心領域内のヒストグラムを求める。そして、求めたヒストグラムを反映した画像情報を 生成し、切替加算部 9を介して画像表示器 10に表示して検者に診断情報として提供 する。このヒストグラムを反映した画像情報には、注目部位を含む関心領域の各部位 の弾性の統計的な特徴を含めることができる。
[0039] ここで、ヒストグラム評価部 20における一連の処理の流れを、図 4を参照して説明す る。ヒストグラム評価部 20は、メモリ部 31と、ヒストグラム演算部 32と、統計評価部 33と 、画像構築部 34を備えて構成されている。メモリ部 31は、歪み Z弾性率演算部 13か ら出力される弾性フレームデータをー且格納して、ヒストグラム演算部 32に出力する ようになっている。
[0040] ヒストグラム演算部 32は、メモリ部 31から出力される弾性フレームデータと、装置制 御インターフェイス部 17から出力される関心領域 (ROI)の座標データを入力する。こ こで、 ROIは、ステップ S5で画像表示器 10に表示されたカラーの弾性画像上で、装 置制御インターフェイス部 17から入力設定する(S6)。
[0041] そして、ヒストグラム演算部 32は、 ROI内に分布する弾性フレームデータを母集団と するヒストグラムを演算し、数値データカゝらなるヒストグラムデータを生成して、統計評 価部 33に出力する。すなわち、 ROIにおける歪み又は弾性率の値を複数の区間に 分割し、各区間に対応する歪み又は弾性率の値を有する計測点の個数を算出し、 例えば、歪み又は弾性率の複数の区間を横軸に、計測点の個数を縦軸に割り付け てヒストグラムデータを生成する。
[0042] 次に、統計評価部 33は、ヒストグラム演算部 32から出力されるヒストグラムデータを 解析して、ヒストグラムデータの統計的な特徴を抽出する。そして、統計的な特徴の 統計評価データを生成して、ヒストグラムデータと共に画像構築部 34に出力する。画
像構築部 34は、統計評価部 33から出力される統計評価データとヒストグラムデータ を反映した画像をヒストグラム画像データとして構築し、切替加算部 9に出力する。
[0043] ここで、ヒストグラム評価部 20における処理を、図 5〜図 8を参照して具体的に説明 する。まず、歪み Z弾性率演算部 13において演算された時刻 tにおける弾性フレー ムデータの弾性率の数値データ群を、図 5の概念図に示すように、 Ei,j (t)とする。こ こで、 i= l,2,3, ' . '、 Nの自然数、 j=l,2,3, ' . '、 Mの自然数である。つまり、図 5にお いて、指標 iは弾性画像の横軸方向の座標を、指標 jは弾性画像の縦軸方向の座標 を示している。
[0044] メモリ部 31では弾性フレームデータ Ei,j (t)を確保して、ヒストグラム演算部 32に出 力する。一方、画像表示器 10に表示された弾性画像上において、装置制御インター フェイス部 17の入力手段 (例えばトラックボールなどのデバイス)を介して検者力 例 えば図 6に示すように、関心領域 ROIを設定する。装置制御インターフェイス部 17は 、設定された ROIの座標、例えば図 6に示すように矩形の ROIの場合は、 4つの角の 座標 Pl (xl,yl)、 P2 (x2,yl)、 P3 (xl,y2)、 P4 (x2,y2)を読み取って、 ROI座標 データとしてヒストグラム演算部 32に出力する。
[0045] ヒストグラム演算部 32は、弾性フレームデータと ROI座標データを用いて、ヒストグラ ムデータを演算する。まず、 ROI座標データを参照して、弾性フレームデータの中か ら ROIの内部に属する計測点の弾性率の数値データ群を、例えば、図 7に示すよう に抽出する。
[0046] 次に、 ROIの内部に属する計測点の弾性率の数値データ群を母集団として、ヒスト グラムデータを演算する。このヒストグラムデータを演算する処理は、一般的なヒストグ ラムを構築する際の手順と同様である。例えば、図 8を参照して具体的に説明する。 ヒストグラムの表示範囲を弾性率で、 O (kPa)力も 700 (kPa)とし、ヒストグラムの区間 であるビン幅を 25 (kPa)にすると、 28個のビンが設定される。すなわち、各ビンの幅 は、次のとおりとなる。
[0047] ビン 1 : 0— 24.99 · · · (kPa)
ビン 2 : 25-49.99 · · · (kPa)
ビン 3 : 50- 74.99 · · · (kPa)
ビン 27 : 650-674.99 · · · (kPa)
ビン 28 : 675-699.99 · · · (kPa)
このように設定されたビンに対応させて、母集団の中力もそれぞれのビンの範囲の 弾性率を有する計測点の個数を計上する。すなわち、例えば、次のように計上される
[0048] ビン 1 : 17 (個)
ビン 2 : 87 (個)
ビン 3 : 142 (個) ビン 27 : 74 (個)
ビン 28 : 63 (個)
このような処理によって、ビン番号 (ビンの範囲)と計測点の個数の関係がヒストグラ ムデータとして生成される。
[0049] 統計評価部 33は、生成されたヒストグラムデータの統計的特徴を演算して、統計評 価データを作成する。この統計的特徴としては、母集団としての弾性率データ群の標 本である計測点の個数、弾性率データの平均値、中央値、分散値、標準偏差、ピー クの半値幅、歪度、尖度などから、適宜選択して設定できる。但し、統計的特徴を表 すものはこれらに限らな 、ことは 、うまでもな!/、。
(ステップ S8)
画像構築部 34は、ヒストグラム演算部 32で求められたヒストグラムデータと、統計評 価部 33で求められた統計評価データとに基づ 、て、ヒストグラム画像データを構築 する。構築されたヒストグラム画像データは、切替加算部 9を介して画像表示器 10〖こ 表示され、組織鑑別の補助弾性情報として検者に提供される。
[0050] なお、本実施形態では、被検体 1の体表面に探触子 2を接触させて計測する場合 について説明したが、本発明はこれに限らず、経直腸探触子、経食道探触子、術中 用探触子、血管内探触子など、任意の超音波探触子にて同様に適用できる。
[0051] 以下、ヒストグラム評価部 20により生成されるヒストグラム画像データの具体的な実
施例について説明する。なお、以下の実施例では、 ROIを図 6に示すように設定した 場合を例に説明する。
実施例 1
[0052] 本実施例 1は、図 8に示すようなヒストグラム画像を生成する具体的例である。図 6に 示すように、 ROIの中には、組織 1、組織 2、組織 3が含まれている。これらの組織 1〜 3の組織固有の硬さの情報は、組織 1〜3の範囲表示とともに、弾性画像として画像 表示器 10にリアルタイムに表示される。これにより、
(組織 1の硬さ) < (組織 2の硬さ) < (組織 3の硬さ)
の関係を容易に把握することができる。
[0053] この ROI内の弾性率データに基づいて生成されるヒストグラムは、図 8に示すような ヒストグラム画像になる。すなわち、同図のヒストグラム画像の縦軸は、それぞれのビン の弾性率に属する計測点の個数に相関する。したがって、個数が多いほどビンの長 さが長くなつている。 ROIにより指定される弾性フレームデータの中の母集団は、組 織 1〜3にそれぞれ固有の弾性率データを保有している。したがって、対応する弾性 率の範囲の計測点の個数が高くなつたヒストグラムが形成される。
[0054] また、統計評価部 33において、例えば母集団に属する計測点の個数、母集団の 弾性率の平均値、中央値、分散値が演算され、図 8に示すように、ヒストグラムに合わ せて、画像表示器 10にリアルタイムに表示される。
[0055] また、切替加算部 9は、例えば図 9に示すように、種々の表示形態により、画像表示 器 10に画像を表示する。図示例は、白黒スキャンコンバータ 8から出力された断層像 と、カラースキャンコンバータ 15から出力された弾性画像と、ヒストグラム評価部 20か ら出力されたヒストグラム画像が組み合わされ、 1つの表示画像データが構築されて 画像表示器 10にリアルタイムに表示される。
[0056] これにより、本実施例によれば、注目する領域が、その周囲の組織とどのような硬さ の関係で結合している力を客観的、定量的に評価することが可能となり、確定的な組 織鑑別を支援することができる。
[0057] すなわち、輝度 (色調)の大小やその広がり方などを表示する弾性画像の場合は、 検者の主観に依存した診断により、検者間で異なる診断結果が下されることがあった
力 本実施例によれば、そのような事象を回避して、客観的、普遍的な診断を確立し て、臨床上有用な超音波診断装置を提供することができる。
[0058] 特に、ヒストグラムには、組織の集団的挙動の特徴が現れることから、これにより、病 変部などの注目部位の組織の硬さの均一性などを評価することができ、非常に硬 、 充実性組織であるガン腫瘍、軟らかい充実性組織である繊維腺種、あるいは流動的 な嚢胞性組織である嚢胞、等々の組織性状の識別が可能となり、臨床上有用な診断 情報を提供することができる。
実施例 2
[0059] 実施例 1においては、弾性率を用いた例を示した力 本発明は、これに限らず、弹 性データとして歪みデータ、変位データ、変位ベクトルにも適用できる。つまり、歪み Z弾性率演算部 13にお 、て演算された歪みの弾性フレームデータに基づ!/、て、同 様のヒストグラム及び統計的特徴データを求めて、画像を表示するようにすることがで きる。また、変位計測部 12において演算された変位フレームデータに基づいて、同 様のヒストグラム及び統計的特徴データを求めて、画像を表示するようにすることがで きる。また、変位フレームデータが変位ベクトルにより構成されていれば、そのべタト ルデータに基づいて、同様のヒストグラム及び統計的特徴データを求めて、画像を表 示するよう〖こすることがでさる。
[0060] つまり、組織の硬さに相関する弾性データに基づ!/、て、ヒストグラム及び統計的特 徴データを求めて画像を表示することにより、実施例 1と同様の効果を奏することがで きる。
実施例 3
[0061] 実施例 1では、歪み Z弾性率演算部 13において演算された弾性率の数値により構 成される弾性フレームデータに基づ 、て、ヒストグラム及び統計的特徴データを求め る例を説明した。しかし、本発明はこれに限らず、弾性データ処理部 14において画 像処理が施された後の弾性フレームデータをヒストグラム評価部 20に入力するように して、ヒストグラム及び統計的特徴データを求めても、実施例 1と同様の効果を奏する ことができる。
[0062] 更に、カラースキャンコンバータ 15において変換された弾性画像データに基づい
て、ヒストグラム及び統計的特徴データを求めても、実施例 1と同様の効果を奏するこ とがでさる。
実施例 4
[0063] 図 10に、実施例 4のヒストグラム画像の表示例を示す。実施例 1においては、弾性 に相関する弾性データである弾性率に基づ 、て、ヒストグラムを生成して表示する例 を説明した。しかし、本発明はこれに限らず、統計評価部 33において、求められたヒ ストグラムを、例えばダブルガウシアンや複数のガウシアン関数などの適切な関数で 近似(フィッティング)して、図 10に示すように画像表示することができる。この場合、 それぞれのガウシアン関数の平均値、標準偏差を求め、それらを統計評価データと して出力することができる。つまり、同図に示すように、フィッティング関数と、その関数 のピークの位置の弾性率 ml〜m3と、あるいは弾性率の平均値 Mと、標準偏差 σ 1 〜 σ 3等の幅情報を求めて、ヒストグラムの特徴量データとして表示することができる
[0064] 本実施例によれば、実施例 1に比べて、注目する領域が、その周囲の組織とどのよ うな硬さの関係で結合しているかを、一層、客観的、定量的に評価することが可能と なり、確定的な組織鑑別を支援することができる。
実施例 5
[0065] 図 11に、実施例 5のヒストグラム画像の表示例を示す。本実施例は、実施例 1に示 したヒストグラムの任意の弹'性率のヒストグラムを指定して、具体的なヒストグラム情報 を表示させるようにしたものである。すなわち、装置制御インターフェイス部 17の入力 手段を介して、ヒストグラム画像データ上に複数のカーソル (Cl、 C2、 · · を表示さ せ、検者がカーソルを自由に動かすことができるようにする。そして、検者が指定した カーソル Cl、 C2、 · · '位置のヒストグラムの情報、例えば、個数 Nl, N2、平均値 ml 、 m2などの数値データがカーソルの動きに同期して表示する。
[0066] 本実施例によれば、注目する領域が、その周囲の組織とどのような硬さの関係で結 合しているかを、一層、客観的、定量的に評価することが可能となり、確定的な組織 鑑別を支援することができる。
実施例 6
[0067] 図 12に、実施例 6のヒストグラム画像の表示例を示す。本実施例は、実施例 5に示 したカーソル Cl、 C2- · 'でヒストグラムの範囲を指定し、その範囲の統計的特徴を演 算して画像表示する例である。
[0068] すなわち、例えば図 12に示すように、カーソル C1〜C2までの範囲をヒストグラムに 指定すると、この範囲に対応するヒストグラムデータの計測点の個数 N、平均値 m、標 準偏差 σなどの統計的特徴を演算して、数値データとして表示する。
[0069] また、検者によって自由に指定されたカーソル C1〜C2までの範囲に限定して、実 施例 4で説明した近似関数によるフィッティングを施して、その統計的特徴を抽出して 画像表示することができる。
実施例 7
[0070] 図 13に、実施例 7のヒストグラム画像の表示例を示す。本実施例は、実施例 6にお V、て、選択されたヒストグラムのカーソル C1〜C2までの範囲に計上された計測点の 集団が、弾性画像のどの領域に対応しているのかを、図 13に示すように、網目で示 した領域 41のように表示するようにしたものである。これにより、検者はヒストグラムと 組織の対応関係を容易に認識できる。
実施例 8
[0071] 図 14に、実施例 8のヒストグラム画像の表示例を示す。本実施例は、ヒストグラムの 2 つのピーク位置 CI、 C2の弾性率の比率 (m2Zml)を演算して表示するようにした 例である。本実施例によれば、例えば、 C1が正常組織の領域であるとすると、 C2が 正常組織に対して何倍硬いかを即座に把握できる。
[0072] すなわち、従来は、生体組織を圧迫し、内部の各領域に生じた歪みや弾性率を演 算し、その大きさを階調化した弾性画像を構築して表示している。しかし、注目する 領域力 その周囲の組織とどのような硬さの関係で結合しているかを把握することで きない。この点、本実施例によれば、注目する領域力 その周囲の組織とどのような 硬さの関係で結合している力を客観的、定量的に評価することが可能となり、確定的 な組織鑑別を行うことができる。
[0073] また、同時に、弾性画像における輝度 (色調)の大小やその広がり方など、検者の 主観に依存した診断により、検者間で異なる診断結果が下されることがあつたが、こ
のような事象を回避して、客観的、普遍的な診断を確立して、臨床上有用な超音波 診断装置を提供することができる。
実施例 9
[0074] 本実施例は、実施例 8に加えて、ヒストグラム評価部 20に次の機能をカ卩えることがで きる。
[0075] 図 14のカーソル Cl、 C2等で指定したピーク部の弾性率の半値幅を求め、弾性率 の半値幅の値に基づいて、組織性状を評価し、その結果を出力表示する。すなわち 、例えば、弾性率の半値幅が l〜20 (kPa)の場合は、非常に硬い充実性組織である ガン腫瘍と評価する。弾性率の半値幅が 21〜70 (kPa)の場合は、軟らかい充実性 組織である繊維腺種と評価する。弾性率の半値幅が 71 (kPa)以上の場合は、流動 的な嚢胞性組織である嚢胞と評価する。
[0076] また、ピーク部の分散値 Z平均値の比を求めて無次元化する。そして、例えば、ガ ン腫瘍の場合は、分散が小さぐ平均が大きいことから、分散値 Z平均値の比が小さ く(例えば、 0. 1)なることで評価できる。繊維腺種の場合は、分散が中程度で、平均 も中程度であるから、分散値 Z平均値の比も中程度 (例えば、 0. 3)なることで評価 する。嚢胞の場合は、分散が大きぐ平均は小さいことから、分散値 Z平均値の比が 大きく(例えば、 0. 5)なることで評価する。
[0077] このように、本実施例によれば、組織性状を評価した結果を出力表示することがで きるから、組織鑑別の信頼度を一層向上させることができる。
実施例 10
[0078] 図 15に、実施例 10のヒストグラム画像の表示例を示す。本実施例は、圧力に応じ て変位する組織の移動に追従させて、弾性画像に設定した ROIを自動的に変更し、 合わせて変更された関心領域におけるヒストグラムを対応させて表示することを特徴 とする。これにより、同一組織領域をトラッキングして ROIが設定されることから、適切 なヒストグラム情報を得ることができる。
[0079] つまり、一般に、超音波を用いた弾性の画像ィ匕においては、被検体の診断部位の 体表に接触させた超音波探触子を用手的に上下動させて被検体を圧迫しつつ、そ の圧迫過程の弾性画像を演算して表示する。この場合に、図 16に示すように、設定
した ROIの座標データを固定したままにすると、圧力に応じて変位する組織が ROIか ら外れ、又は ROI外の組織力 ¾OI内に入り、 ROI内の母集団が変化することになる から、そのまま、 ROI内の計測点のヒストグラムを演算しても、客観的に対比すること ができない。したがって、時間の経過とともに圧迫力が相違する弾性画像を対比して も、的確な鑑別情報が得られないことになる。
[0080] そこで、本実施例では、 ROIの座標データを、被検体の圧迫動作状況や、拍動な どの生体の運動状況に応答して更新するようにしている。つまり、ヒストグラム評価部 2 0には、装置制御インターフェイス部 17により設定された ROIの領域を示す代表点と して、例えば、図 6に示したように、 ROI設定時刻 tOにおいて、 Pl (tO) (xl,yl)、 P2 (tO) (x2,yl)、 P3 (t0) (xl,y2)、 P4 (t0) (x2,y2)力 S人力される。次に、時亥 Ijtlにお いて、図 1の変位計測部 12から時刻 tlにおける変位フレームデータ Di,j (tl)が入力 される。時刻 tlに取得された変位フレームデータ Di,j (tl)において、対応する座標 の変位データを参照することにより、時刻 tOにおいて ROI座標 Pl (tO) (xl,yl)、 P2 (tO) (x2,yl)、 P3 (t0) (xl,y2)、 P4 (t0) (x2,y2)に位置した組織が、時刻 tlにお いてどの位置に変位したかを求めることができる。そこで、求められた座標 Pl (tl)、 P2 (tl)、 P3 (tl)、 P4 (tl)を新たな ROI座標データとして更新する。そして、更新さ れた ROI座標内に分布する母集団を用いて、ヒストグラムと統計的特徴を求めて表示 する。
[0081] このように、本実施例によれば、生体を圧迫することにより、生体内の組織は、圧迫 過程において超音波探触子に近づくことになる力 本実施例の処理を任意の時刻に おいて繰り返すことにより、被検体の圧迫動作状況や、拍動などの生体の運動状況 に応じて ROI座標データが更新される。これにより、例えば、図 15に示すように、 ROI の位置と大きさは時刻により変化する力 任意の時刻において生体内の同一組織領 域を抽出することになる。
[0082] このように、本実施例によれば、同一組織領域に対応する弾性フレームデータに基 づいて母集団を構成し、ヒストグラムと統計的特徴を演算して表示するから、圧迫過 程において ROI内に別の組織が入らない。その結果、母集団の標本数は変化する 力 組織の特徴を区別するのに重要な情報としての平均値や分散値などの統計的
特徴は任意の時刻においても安定的な値が保持され、組織固有の値を確定的に評 価することができる。
[0083] つまり、図 16に示したように、 ROIの座標を固定した場合は、演算されたヒストグラム は時刻によって異なったものとなり、またヒストグラムに基づいて演算される統計的特 徴を表す指標も時刻により変動するから、客観的、確定的な診断を妨げることがある 。しかし、本実施例によれば、圧迫過程に応じて ROI座標データを更新し、任意の時 刻において同一組織領域を抽出してその統計的特徴を評価することができるから、 客観的、確定的な診断を高精度に行うことができる。
実施例 11
[0084] ここで、ヒストグラム画像の表示例の他の実施例 11を説明する。本実施例は、ヒスト グラム評価部 20において、過去に求めたヒストグラムデータを保持しておき、各ビン 毎に現時刻におけるヒストグラムデータと、過去のヒストグラムデータを時間的に累積 した累積ヒストグラムデータを構成するようにすることを特徴とする。また、累積ヒストグ ラムデータの統計的特徴を求めて、画像表示することができる。
[0085] 更に、ヒストグラム評価部 20において、過去に求めたヒストグラムデータを保持して おき、各ビン毎に現時刻におけるヒストグラムデータと、過去のヒストグラムデータを時 間的に平滑処理した平滑ヒストグラムデータを構成するようにしてもょ 、。この場合も、 平滑ヒストグラムデータの統計的特徴を求めて、画像表示するようにしてもょ 、。
[0086] 本実施例によれば、安定的に高精度な統計的特微量を抽出できる。
実施例 12
[0087] 上記の各実施例では、ヒストグラムと統計的特徴をリアルタイムで演算して、かつ表 示する場合を例に説明した力 本発明はこれに限らず、シネメモリ部 18に格納した画 像データを用いて図 2の所要の処理を行うようにすることができる。
[0088] この場合は、装置制御インターフェイス部 17を介して装置をフリーズして行う。シネ メモリ部 18は、装置制御インターフェイス部 17の制御信号に従って、過去の表示画 像データを呼び出し、任意のフレームの表示画像データの中力 ヒストグラムによる 評価を行う表示画像データと、 ROIを指定する。そして、 ROIの座標データをヒストグ ラム評価部 20に出力する。また、処理対象のフレームを指定する情報を歪み Z弾性
率演算部 13に出力し、歪み Z弾性率演算部 13に確保されて 、る過去の弾性フレー ムデータを読み出して、ヒストグラム評価部 20に送ってヒストグラムと統計的特徴の演 算を行わせる。これにより、オフラインで本発明を実施することができる。
実施例 13
[0089] 上記実施例で説明したヒストグラムの表示の採否、あるいはヒストグラムの表示設定
(ビン幅、縦軸'横軸の範囲など)は、任意に設定、変更することができる。また、統計 的特徴の平均値、標準偏差などのどの統計情報を演算して表示するかは、装置制 御インターフェイス部 17によって自由に指定できる。
実施例 14
[0090] 図 17に、本発明の弾性に相関する弾性データの一つである変位ベクトルを用いて ヒストグラムを生成する実施例 14の表示例を示す。本実施例は、被検体の組織に圧 力を加えて計測された超音波断層データに基づいて、被検体の断層部位の複数の 計測点における組織の変位ベクトルを求め、断層部位の計測点に対応させて変位べ タトルを表示する弾性画像を生成して表示し、弾性画像に設定される関心領域にお ける特定方向である深度方向の変位ベクトル成分の分布をヒストグラムで表示するも のである。なお、同図におけるヒストグラムの横軸は、変位の深度方向成分を表して いる。
[0091] 図 17 (A)は、 ROI内に設定された注目領域力 非常に硬い充実性組織であるガン 腫瘍の場合を示している。この場合の各計測点の変位ベクトルは、図 17 (A)の上図 に矢印で示したように、同一の方向を向いている。したがって、 ROI内の変位ベクトル の分布を示すヒストグラムは、図 17 (A)の下図に示すように、狭い範囲に集中したヒ ストグラムになる。つまり、この例は、組織要素間で強く結合しているから、圧迫を受け ても、全体として変形せずに、運動することが理解できる。
[0092] 図 17 (B)は、 ROI内に設定された注目領域が、軟らかい充実性組織である繊維腺 種の場合を示している。この場合の各計測点の変位ベクトルは、図 17 (B)の上図に 矢印で示したように、やや方向が分散している。したがって、 ROI内の変位ベクトルの 分布を示すヒストグラムは、図 17 (B)の下図に示すように、分布の幅がやや広くなつ ている。つまり、この例は、組織要素間である程度結合しており、圧迫を受けたときに
、互いに束縛しあいながらも、ある程度変形することが理解できる。
[0093] 図 17 (C)は、 ROI内に設定された注目領域が、流動的な嚢胞性組織である嚢胞の 場合を示している。この場合の各計測点の変位ベクトルは、図 17 (C)の上図に矢印 で示したように、変位の方向が大きく分散している。したがって、 ROI内の変位べタト ルの分布を示すヒストグラムは、図 17 (C)の下図に示すように、分布の幅が大きな範 囲に広がっている。つまり、この例は、組織要素間で全く結合していないため、圧迫 を受けると、各要素群はそれぞれ自由な運動をすることが理解できる。
[0094] これらのことから、変位ベクトルの深度方向成分又は横方向成分でヒストグラムを構 成すれば、組織性状の違いを識別することができる。特に、硬さの違いだけでなぐ 分布幅を観察することにより、組織要素間の結合の程度が識別できるから、組織鑑別 の有効な補助弾性情報が得られる。
実施例 15
[0095] 図 18に、本発明の実施例 15に係る弾性情報の表示例を示す。本実施例は、弾性 画像の色相又は輝度の階調 (例えば、 256階調)の割り付け範囲を、関心領域 ROI 内のヒストグラムの弾性データの幅に対応させて変更することにより、関心領域 ROI 内の弾性の違いのコントラストを際立たせるようにした実施例である。すなわち、図 18 (a)は視野内の弾性画像とヒストグラムを表示したものであり、同図(b)は弾性画像中 に関心領域 ROIをマニュアル又は境界検出によって自動で設定し、関心領域 ROI 内の計測点のヒストグラムを表示したものである。これに対し、同図(c)に示すように、 関心領域 ROI内の計測点のヒストグラムの弾性データの幅に対応させて、階調の割り 付け範囲が自動で変更され、 ROI内の弾性の違いがコントラストを際立たせて表示さ れる。
[0096] 本実施例によれば、 ROI内の組織の硬さ分布のコントラストが強く表示されるから、 組織鑑別の有効な補助弾性情報が得られる。
実施例 16
[0097] 図 19に、本発明の実施例 16に係る弾性情報の表示例を示す。図 9に示した実施 例 1の弾性情報の表示例では、弾性画像とヒストグラムを含む補助弾性情報の画像( ヒストグラム画像)を 1つの画面に並べて表示する例を示した。本実施例では、図 19
に示すように、ヒストグラム画像を小さなウィンドウ表示により、弾性画像の一部に重ね て表示している。
[0098] すなわち、ヒストグラム画像は、解像度を犠牲にしても見落としが生じることはないの で、断層像や弾性画像と同等の大きさで対等に表示する必要はない。そこで、本実 施例のように、弾性画像の画像領域内に、縮小したヒストグラム画像を小さなウィンド ゥで表示することができる。これによれば、弾性画像の弾性情報を十分な解像度で観 察できるから、情報の損失を低減できる。
[0099] なお、図 19の例によれば、ヒストグラム画像が重畳表示された部分の弾性画像の情 報が隠されてしまい、場合によっては重要な弾性情報を失ってしまうおそれがある。こ のような場合は、図 20に示すように、ヒストグラム画像を半透明画像で構成し、弾性画 像に重畳させて表示することが好ましい。これによれば、ヒストグラム画像を重畳して 表示したことによる弾性情報の見落としを回避することができる。
[0100] ここで、カラーの半透明画像を生成するアルファブレンデイングと呼ばれる一般的な 画像構築方法について説明する。例えば、画素のカラー情報を三原色 (R, G, B)で 表現すると、カラーの弾性画像の座標 (i, j)の画素のカラー情報を E(i,j)は、次式の ように、 3原色のカラー情報 Er(i,j), Eg(i,j), Eb(i,j)の関数となる。
E(i,j) = (Er(i,j), Eg(i,j), Eb(i,j))
同様に、カラーのヒストグラム画像の座標 (i, j)の画素のカラー情報を H(i, j)は、次 式の関数となる。
H (i,j) = (Hr (i,j) , Hg (i,j) , Hb (i,j) )
ここで、両画像間の重畳比率を α (0≤ α≤ 1)とすると、ヒストグラム画像を重畳した 領域の半透明重畳画像の座標 (i, j)の画素のカラー情報 G(i, j)は、次式で求めるこ とがでさる。
C(i,j) = a XE(i,j) + (l-a) XH(i,j)
つまり、半透明重畳画像の座標 (i, j)の画素のカラー情報は、弾性画像とヒストグラム 画像の両方の画像情報を aの比率で有するものとなり、それぞれの画像のカラー情 報を半透明的に同時に認識することができる。
[0101] このように、本実施例によれば、限られたモニタ表示画面内に複数の画像を同等の
大きさで干渉しな 、ように配置すると、それぞれの画像が小さくなつてしまう場合でも 、画像診断の精度を落とすおそれを回避できる。
実施例 17
[0102] 図 21に、本発明の実施例 17に係る弾性情報の表示例を示す。図 19及び図 20で は、断層像や弾性画像等の診断画像の領域内に縮小したヒストグラム画像をウィンド ゥで表示する例を示した。本実施例では、図 21に示すように、必要に応じて、診断画 像とヒストグラム画像の表示領域の大小関係を切り替えるようにして 、る。
[0103] 例えば、超音波診断装置の装置制御インターフェイス部 17に備えられた所定のキ 一やボタンに設定された画像切り替えキーを操作することにより、診断画像とヒストグ ラム画像の表示領域の大小関係が切り替わるように構成することができる。なお、画 像切り替えキーは、装置制御インターフェイス部 17のキーやボタンに限らず、モニタ の画像上に表示されアイコンに設定することができる。
[0104] 本実施例によれば、図 11、図 12の実施例において、ヒストグラム画像上における力 一ソル移動などの操作を、大きなヒストグラム画像上で効率的に実行することができる 実施例 18
[0105] 図 22に、本発明を甲状腺の濾胞癌を鑑別する診断に適用した場合の弾性情報の 表示例を示す。甲状腺の画像診断においては、甲状腺濾胞癌の鑑別が非常に重要 になっている。濾胞癌は腫瘍内部の細胞密度に分布があり、腫瘍辺縁部は腫瘍中 心部よりも高い細胞密度を有することが知られている。一方、弾性画像により濾胞癌 を評価した場合、腫瘍辺縁が腫瘍中心部に比較して硬!ヽ領域として画像表示される 傾向がある。つまり、細胞密度が高い部位を細胞密度が低い部位よりも硬い部位とし て計測されるものと考えられる。
[0106] 本実施例は、腫瘍内部の弾性分布の特徴を評価する場合に適切な弾性情報を表 示するものである。例えば、甲状腺腫瘍のように、腫瘍辺縁部と中心部の間に硬さの 違いがある場合でも、通常の弾性画像ではその違いを認識できない場合も多い。図 22は、甲状腺腫瘍の弾性画像と濾胞癌の場合のヒストグラム画像の表示例である。 同図に示すように、濾胞癌を客観的に鑑別するには、腫瘍中心部 ROI1と環状の腫
瘍辺縁部 ROI2をそれぞれ独立した母集団としてヒストグラム解析することにより、鑑 別に有効な情報を取得して表示する。
[0107] すなわち、本実施例では、次の表示を特徴とする。
( 1)濾胞癌と思われる領域の辺縁部と中心部に、それぞれ独立に ROI 1と ROI2を設 定して表示する。
(2) ROIlと ROI2のヒストグラムを、 ROIごとに色分けして合成して表示する。
(3)それぞれのヒストグラムの統計情報を演算して表示する。図示例では、平均、標 準偏差、中央値が表示されている。
[0108] ROIl、 ROI2の設定法は様々考えられる。例えば、装置制御インターフェイス 17の マウスやトラックボールなどの入力手段を操作して、検者により環状の辺縁部 ROI1 の外周と内周を線で描いて指定する。これにより、その描画線の間の環状領域が辺 縁部 ROI1として認識し、設定するようにすることができる。
[0109] また、楕円形の複数のテンプレートを予め準備しておき、マウスやトラックボールを 操作して、辺縁部 ROI1の外周と内周に 2つの楕円形テンプレートを合わせるように 指定する。これにより、 2つの楕円形テンプレートの外周線に挟まれる環状領域が辺 縁部 ROI1とし、小さい楕円形テンプレートの内側領域を RROI2として認識させて、 設定するよう〖こすることができる。
[0110] さらに、辺縁領域指定用に適した複数のトーラス (輪環体)状の ROIを予め準備して おき、マウスやトラックボールなどで操作して、 2つのトーラス ROIを画像上に配置す る。これ〖こより、 2つのトーラス ROIに挟まれる環状領域が辺縁部 ROI1とし、小さいト 一ラス ROIの内側領域を RROI2として認識させて、設定するようにすることができる。
[0111] 本実施例によれば、腫瘍の辺縁部と中心部の硬さの違いを定量的に認識すること ができな!/ヽ弾性画像であっても、定量的で客観的な濾胞癌の診断を行うことができる 実施例 19
[0112] 例 は、文献「Ueno E, et al: Dynamic tests in real-time oreast echography. Ultras ound Med Biol, 14: 53-57, 1988.」に記載されているように、乳腺領域における画像 診断では、腫瘍の可動性 (モピリティーと呼ばれる。)を評価するダイナミックテストと
いわれる検査が知られている。ダイナミックテストとは、関心腫瘍がその周囲の組織に 対してどの程度拘束されているかを評価するものである。具体的には、***全体を大 きく変形させたときに、注目対象の腫瘍が周囲組織の動きに引きずられて動くか、そ れとも、周囲組織の動きに引きずられずに、独立した動きを示すかを、断層像をリア ルタイムで観察しながら評価する診断法である。例えば、乳癌の中でも浸潤性の腫 瘍は、周囲の組織に浸潤して発育するため、周囲組織との拘束度が高ぐ可動性が 低い。一方、良性の腫瘍は浸潤性を示さないため、周囲組織に結合された挙動は示 さず、可動性が高いといわれる。しかし、これらの可動性の評価は、検者の経験、主 観に依存しているため、客観性が低い。
[0113] そこで、本実施例では、ヒストグラムを用いて可動性の評価を行うようにしている。図 23 (a) , (b)に本実施例の弾性情報の表示例を示す。同図 (a)に示す関心部は、良 性腫瘍などで可動性が高い腫瘍の場合を想定している。このとき、腫瘍は周囲の乳 腺部には拘束されていないので、外部力も***全体を大きく変形させたとき、腫瘍と 乳腺の境界部ではスリップする挙動が現れる。つまり、境界の腫瘍側の領域と境界の 乳腺側の領域において、変位の方向と大きさに違いが現れ、場合によっては逆方向 に変位する。そこで、同図(a)に示すように、腫瘍と乳腺の境界を含む領域に ROIを 設定し、その ROI内に例えば 0. 05mmの深度方向の変位が与えられたときの圧迫 操作を、評価の基準として設定する。このときの ROI内の横方向の変位の大きさをヒ ストグラムとして表示すると、図に示すように、正の方向(例えば、画像で右向き)と負 の方向(例えば、画像で左向き)に、同程度の大きさの変位を持ったヒストグラムのピ ークが 2つ現れる。このヒストグラムの広がり方 (ピークの半値幅、標準偏差や分散な ど)や、ピーク間の距離 (mm) 1S 可動性の程度を与える目安として評価することがで きる。
[0114] 一方、同図 (b)に示す関心部は、悪性腫瘍の可動性が低い腫瘍の場合である。こ の場合、腫瘍は周囲の乳腺部には強く拘束されているので、外部から***全体を大 きく変形させたとき、腫瘍と乳腺の境界部は同等の運動をする。つまり、境界の腫瘍 側の領域と境界の乳腺側の領域において、変位の方向と大きさの違いが小さぐ場 合によっては完全に結合して変位する。そこで、同図(a)と同様に、腫瘍と乳腺の境
界を含む領域に ROIを設定し、その ROI内に 0. 05mmの深度方向の変位が与えら れたときの圧迫操作を評価の基準とし、 ROI内の横方向の変位の大きさをヒストグラ ムを表示すると、同図(このようになる。すなわち、正の方向(例えば、画像で右向き) 力 負の方向(例えば、画像で左向き)のいずれか一方に、同程度の大きさの変位を 持ったヒストグラムのピークが現れ、ヒストグラムの広がり方 (ピークの半値幅、標準偏 差や分散など)を可動性の程度を与える目安として評価することができる。
実施例 20
[0115] 実施例 19では、良悪性鑑別の 1つの情報として、腫瘍の周囲への浸潤度に関係す る可動性を客観的に評価するために、変位情報に基づいたヒストグラムを補助弾性 情報として提供する例を説明した。本実施例では、良悪性鑑別の 1つの情報として、 組織間の癒着を評価するための補助弾性情報として、ヒストグラムを適用する方法を 示す。
[0116] 図 24は、乳腺組織に腫瘍がある場合の弾性画像の例を示している。特に、乳腺と 大胸筋の境界部の近くにある腫瘍を摘出する手術を行う場合に適した補助弾性情報 の例である。乳腺と大胸筋の境界部の近くにある腫瘍を摘出する場合、乳腺と大胸 筋の間から腫瘍にアプローチすることがある。このような手術を行う場合、腫瘍と大胸 筋、及び乳腺と大胸筋が癒着しておらず、剥離可能かどうかが非常に重要になる。
[0117] そこで、本実施例では、同図(a)に示すように、乳腺と大胸筋との癒着を評価したい 領域に ROIを設定する。図 23の場合と同様、乳腺と大胸筋が癒着していなければ、 乳腺と大胸筋は***の変形に対してそれぞれ独立の運動を行 、、境界部ではスリツ プする現象が起こる。したがって、例えば所定の大きさの深度方向の変位を与え、横 方向の変位をヒストグラム画像として表示すると、同図(a)のようにヒストグラムのピーク が 2つ現れる。
[0118] 一方、乳腺と大胸筋が強く癒着をしている場合は、同図 (b)に示すように、乳腺と大 胸筋が結合して運動するために、一定の変位のところにピークを持った分布として得 られる。したがって、 ROIにおけるヒストグラムの広がり方 (ピークの半値幅、標準偏差 や分散など)によって、癒着あるいは剥離の可能性を評価する目安になる。
[0119] このように、剥離可能力どうかの評価に変位情報に基づいたヒストグラムを用いるこ
とができ、腫瘍の摘出手術における手術計画を立てるための重要な情報として利用 することができる。
[0120] なお、癒着の評価は摘出手術の場合に限らず、豊胸手術などにおいて乳腺と大胸 筋との間にバッグ (生理食塩水バッグやシリコンバッグなど)を挿入するような場合にも 、組織間の癒着を評価するために効果的に用いることができる。内視鏡を用いた外 科手術においてもその効果は同様である。
実施例 21
[0121] 図 25を参照して、本発明の定量的な補助弾性情報を、高密度焦点式超音波 (以 下、 HIFU: high-intensity focused ultrasound という。 )の治療効果判定に適用した 本実施例 21について説明する。
[0122] HIFUによる治療は、超音波治療用プローブ力 被検体内の治療部位に超音波を 放射し、治療部位の小さな領域に高密度の超音波を集束させ、焦点部位を例えば 6 0〜90°Cの高温に加熱して病変部(例えば、癌などの病巣)の組織を死滅させるもの である。そして、焦点を少しずつ治療部位内で移動させて、病変部の組織全体を死 滅させて治療する。
[0123] この治療により死滅された病変部の組織は熱変性により硬くなるから、その硬さの程 度と治療効果とを対応付けることができる。ところで、治療部位に関心領域 ROIを設 定し、 ROI内の弾性情報を計測して弾性画像を観察すれば、治療部位の組織全域 に渡って万遍なく均一に治療が施されたかどうかを判定できる。しかし、硬さに応じた 階調で表示される弾性画像によれば、背景技術の欄で述べたように、色階調が近似 する領域の微小部位については、階調差を明確に識別できず、治療が施されたかど うかを明確に区別できない場合がある。このような場合に、本発明の補助弾性情報で あるヒストグラムにより、定量的及び客観的に的確な診断を支援することができる。
[0124] 図 25 (a)は、治療前の前立腺癌を含む前立腺の断面を弾性画像で描出している 様子を示す。このとき、前立腺癌の領域に ROIを設定すると、図に示すように、例え ば 100 (kPa)程度を中心とした弾性を有することが分かる。
[0125] 次に、前立腺癌の領域に HIFU治療を施した途中過程の様子を同図 (b)に示す。
十分に超音波が照射された前立腺癌の領域は、熱変性により硬化しており、例えば
同図(b)に示すように 300 (kPa)程度の硬さに変化していることがヒストグラムにより 認識される。しかし、同図(b)のヒストグラム力も把握されるように、未治療の前立腺癌 の硬さを表す 100 (kPa)の硬さをもつ組織領域が ROI内部に残っており、これらの領 域は十分に超音波が照射されな力つた治療不完全な領域であることが分かる。
[0126] このとき、図 13の実施例と同様にして、治療不完全な領域をカーソル Cl、 C2で選 択し、弾性画像にその計測点をフィードバックして表示する処理を行えば、治療不完 全な領域を画像上に示すことができ、検者が認識できる。
[0127] そして、最後に、同図(c)に示すように、熱変性が不十分な治療不完全な領域に対 して、再度、 HIFU治療を施して、ヒストグラムを確認することにより、 ROI内のすべて の領域において均一に万遍なく治療が行われたことを確認することができる。
[0128] なお、治療完了の判断基準として 200 (kPa)以上を閾値として設定することもでき る。これにより、治療効果判定に共通条件を設定することができ、一定の治療基準に 従った確実な治療を実現することができる。
[0129] また、熱で癌組織を焼灼して治療する方法は HIFUに限らず、 RFA (経皮的ラジオ 波焼灼療法)の場合にも、本発明を適用できる。また、逆に癌組織を急速に冷却して 治療するクライオ療法でも、熱変性による硬化を対応づけるだけであるから、本発明 を同様に適用することができる。
[0130] さらに、治療方法に関係なぐ例えば薬剤を利用したホルモン療法なども含めて、 癌組織等の治療において、組織変性により硬さが変化する治療法に対して、その治 療処理が万遍なく行われたかどうか、又は治療効果が万遍なく得られているかどうか を評価する方法として、本発明を利用することができる。
実施例 22
[0131] 図 26を参照して、乳癌検診に好適な本発明の実施例 22を説明する。乳癌の検診 においては、通常、図 26 (a)、(b)に示すように超音波探触子を***に当てて走査し 、***の断層像をリアルタイムでスキャンして、乳癌検診が行われている。このような 検診においては、一人の被験者に十分な検査時間を割くことが困難であり、効率的 な診断が求められる。しかし、リアルタイムにスキャンされる断層像の計測断面内のす ベての領域を、見落としなく観察することは非常に時間を要する作業であり、経験に
も大きく依存する。すなわち、 Bモード像等の断層像は、生体内の組織の構造を画像 として構築した情報であり、腫瘤像を形成しているとか、後方のエコーが増強している とか、辺縁部が粗雑であるとかなど、組織の形態情報を鑑別に利用する。つまり、非 常に多くの複雑な視点に基づいて腫瘍を検出したり、鑑別したりするため、効率的に 診断を行うためには検者の経験が求められる。
[0132] ところで、腫瘍は、周囲の脂肪組織や乳腺組織に比較して硬いため、例えば弾性 情報として弾性率に基づいた検査を行うときには、計測断面内に大きな弾性率を持 つた領域があるかどうかを観察することになる。また、悪性の組織ほど、硬いことが分 力つているので、異常な硬さを持った領域が計測断面に含まれるかどうかだけに注意 を集中すればよいため、リアルタイムでも効率的に認識することが可能である。
[0133] しかし、背景技術の欄で述べたように、弾性画像に基づ!、て乳癌などの腫瘍を弾 性の違いを判別して鑑別しょうとする場合、異なる色階調の境界領域、あるいは、色 階調が近似する領域について、弾性の差異を明確に区別できない場合がある。した がって、弾性画像に基づいて計測断面のすべての領域をくまなぐかつ注意深く観 察する必要があり、リアルタイムに計測断面をスイープしてこれを行うことは検者の負 担が大きい。
[0134] そこで、本実施例では、弾性画像に加えて、あるいは弾性画像に代えて、弾性情報 の分布をヒストグラム化して補助弾性情報として表示することにより、効率的な検診を 実現するようにしている。すなわち、図 26 (c)に示すように、断面のスキャン中にリア ルタイムで歪みや弾性率などの弾性情報を計測し、その弾性情報の分布をヒストグラ ム化して表示する。例えば、悪性を疑う弾性率の閾値 Ethを 200 (kPa)としたとき、ヒ ストグラム分布を見るだけで、現在の弾性画像に存在する Eth以上の弾性率の計測 点の分布を即座に判断することができる。つまり、弾性画像であれば隅々まで確認す ることになるが、本実施例によれば、そのような負担がなくなる。
[0135] さらに、ヒストグラムの縦軸 (度数)〖こも閾値 Cthを設定し、 Eth以上の計測点が Cth 以上計測された場合には、ヒストグラムの表示画像上に「弾性異常領域あり」等の警 告を表示するようにすれば、見落としを防ぐことができる。
[0136] また、ヒストグラムの縦軸(度数)のスケールは、リニアスケールからログスケールに変
更できるようになつていてもよい。ログスケールにすることにより、悪性を疑う硬い腫瘍 が非常に小さい、度数が少ない場合でも、それを見落とすことなく高感度に検出する ことができる。
以上、各実施例で説明したように、変位、歪み、あるいは弾性率等に基づ 、た弾性 画像に加えて、あるいは弾性画像に代えて、弾性分布のヒストグラムを含む補助弾性 情報を画像表示器 10の表示画面に表示することにより、効率的で検者負担の少な V、診断を実現することができる。
Claims
[1] 被検体の組織に圧力を加えて計測された超音波断層データに基づいて、前記被 検体の断層部位の複数の計測点における組織の弾性データを求める弾性データ演 算手段と、該弾性データに基づ!/ヽて前記断層部位の弾性画像を生成する弾性画像 生成手段とを備えた超音波診断装置において、前記弾性データより、前記被検体を 構成する組織の集団的挙動を表す統計的な特徴量を把握するための集団的挙動把 握手段を備えたことを特徴とする超音波診断装置。
[2] 前記集団的挙動把握手段は、前記弾性データに基づいて弾性分布を表して、前 記集団的挙動を把握することを特徴とする請求項 1に記載の超音波診断装置。
[3] 前記集団的挙動は、前記組織の組織性状を反映して ヽることを特徴とする請求項 1 に記載の超音波診断装置。
[4] 前記集団的挙動把握手段は、関心領域に含まれる全ての各計測点の弾性データ を母集団とし、同一又は同程度の弾性データを有する計測点の個数をヒストグラムと してグラフ化する手段を備えたことを特徴とする請求項 1に記載の超音波診断装置。
[5] 前記集団的挙動把握手段は、前記ヒストグラムに限られず、ヒストグラム力 発生し て求められる特徴量を補助弾性情報として表示する手段を備えたことを特徴とする請 求項 4に記載の超音波診断装置。
[6] 前記集団的挙動把握手段は、前記弾性データに画像処理を施した後のデータに つ!ヽて、組織の集団的挙動を表す統計的な特徴量を把握することを特徴とする請求 項 1に記載の超音波診断装置。
[7] 前記集団的挙動把握手段により、注目する前記被検体の領域がその周囲の領域と どのような硬さの関係で結合しているかを評価する手段を備えたことを特徴とする請 求項 1に記載の超音波診断装置。
[8] 前記集団的挙動把握手段は、前記弾性画像の弾性分布を表すヒストグラムを生成 するヒストグラム評価手段と、前記弾性画像と前記ヒストグラムの画像の少なくとも一方 を表示する画像表示手段を備えたことを特徴とする請求項 1に記載の超音波診断装 置。
[9] 前記ヒストグラム評価手段は、前記弾性画像に設定される関心領域における弾性
分布を表すヒストグラムを生成することを特徴とする請求項 8に記載の超音波診断装 置。
[10] 前記弾性データは、前記組織の変位、歪み及び弾性率の 、ずれか一つであること を特徴とする請求項 8又は 9に記載の超音波診断装置。
[11] 前記ヒストグラムは、弾性データを複数の区間に分割し、各区間に対応する前記弾 性データを有する計測点の個数を長さで表示したものであることを特徴とする請求項
8又は 9に記載の超音波診断装置。
[12] 前記弾性画像生成手段は、前記弾性画像に設定された前記関心領域を、前記被 検体に加えられる圧力に応じて変位する前記組織に追従させて自動的に変更するこ とを特徴とする請求項 9に記載の超音波診断装置。
[13] 前記ヒストグラムの画像に、前記ヒストグラムの前記弾性データの平均、標準偏差、 中央値の少なくとも一つを併せて表示する手段を備えたことを特徴とする請求項 8又 は 9に記載の超音波診断装置。
[14] 前記画像表示手段は、前記超音波断層データに基づいて生成された前記断層部 位の断層画像と、前記弾性画像と、前記ヒストグラムの画像とを同一画面に表示する ことを特徴とする請求項 8又は 9に記載の超音波診断装置。
[15] 前記ヒストグラムの画像に、前記ヒストグラムの分布を近似した関数と、前記ヒストグラ ムに表れたピークの前記弾性データと、前記ピークの標準偏差の少なくとも一つを併 せて表示する手段を備えたことを特徴とする請求項 8又は 9に記載の超音波診断装 置。
[16] 前記ヒストグラム評価手段は、前記ヒストグラムの画像上に表示されたカーソルによ り指定される前記弾性データの区間における前記計測点の個数と該区間の弾性デ ータを求め、前記ヒストグラムの画像に数値で表示することを特徴とする請求項 11に 記載の超音波診断装置。
[17] 前記ヒストグラム評価手段は、前記ヒストグラムの画像上に表示された複数のカーソ ルに挟まれる前記弾性データの複数の区間における前記計測点の個数と、該複数 の区間の弾性データの平均値と、該複数の区間の弾性データの標準偏差の少なくと も一つを求め、前記ヒストグラムの画像に数値で表示することを特徴とする請求項 11
に記載の超音波診断装置。
[18] 前記画像表示手段は、前記超音波断層データに基づいて生成された前記断層部 位の断層画像と、前記弾性画像と、前記ヒストグラムとを同一の画面に並べて表示し 前記ヒストグラム評価手段は、前記ヒストグラムの画像上に表示された複数のカーソ ルに挟まれる前記弾性データの複数の区間における前記計測点の個数と、該複数 の区間の弾性データの平均値と、該複数の区間の弾性データの標準偏差の少なくと も一つを求め、前記ヒストグラムの画像に数値で表示すると共に、前記カーソルに挟 まれた区間に属する複数の前記計測点の位置を前記弾性画像に識別可能に表示 することを特徴とする請求項 11に記載の超音波診断装置。
[19] 前記ヒストグラム評価手段は、前記ヒストグラムの画像上に表示された 2つのカーソ ルにより指定される 2つの前記弾性データの区間における前記弾性データの比を求 め、前記ヒストグラムの画像に数値で表示することを特徴とする請求項 11に記載の超 音波診断装置。
[20] 前記ヒストグラム評価手段は、前記ヒストグラムの画像上に表示されたカーソルによ り指定される前記ヒストグラム分布のピーク部の前記弾性データの半値幅を求め、該 半値幅に基づ ヽて前記ピーク部に対応する組織性状を評価し、該評価結果を前記ヒ ストグラムの画像に表示することを特徴とする請求項 11に記載の超音波診断装置。
[21] 前記画像生成手段は、前記圧力に応じて変位する前記組織に追従させて前記弾 性画像に設定された前記関心領域を自動的に変更し、
前記ヒストグラム評価手段は、前記変更された関心領域における前記ヒストグラムを 生成することを特徴とする請求項 9に記載の超音波診断装置。
[22] 前記画像生成手段は、前記関心領域内の前記ヒストグラムが分布する弾性データ の幅に対応させて、前記弾性画像の色相と輝度の少なくとも一方による階調を割り付 け直して、前記弾性画像のコントラストを高くすることを特徴とする請求項 8又は 9に記 載の超音波診断装置。
[23] 前記画像表示手段は、前記ヒストグラムを小さなウィンドウ表示で前記弾性画像に 重ねて又は半透明で表示することを特徴とする請求項 8又は 9に記載の超音波診断
装置。
[24] 前記画像表示手段は、前記断層画像と前記弾性画像と前記ヒストグラムをウィンド ゥ表示するとともに、ウィンドウの大きさを入力指令に応じて切り替え可能に形成され ていることを特徴とする請求項 14に記載の超音波診断装置。
[25] 前記ヒストグラム評価手段は、前記弾性画像の設定される複数の関心領域に対応 する複数のヒストグラムを関心領域ごとに色分けして合成することを特徴とする請求項
8又は 9に記載の超音波診断装置。
[26] 前記弾性データ演算手段は、前記超音波断層データに基づ!、て前記被検体の断 層部位の複数の計測点における組織の変位ベクトルを求め、
前記弾性画像生成手段は、前記断層部位の前記計測点に対応させて前記変位べ タトルを表示する弾性画像を生成し、
前記ヒストグラム評価手段は、前記弾性画像に設定される関心領域における特定 方向の前記変位ベクトル成分の分布のヒストグラムを生成する請求項 8に記載の超音 波診断装置。
[27] 前記関心領域は、異なる組織の境界を含む領域に設定されることを特徴とする請 求項 26に記載の超音波診断装置。
[28] 前記弾性データ演算手段は、被検体の計測断層面を異ならせて計測された超音 波断層データに基づいて、前記被検体の複数の断層面における複数の弾性フレー ムデータを求め、
前記ヒストグラム評価手段は、前記複数の弾性フレームデータの弾性分布を表すヒ ストグラムを生成することを特徴とする請求項 8に記載の超音波診断装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2006800388383A CN101291629B (zh) | 2005-10-19 | 2006-10-10 | 超声波诊断装置 |
US12/083,732 US20090143676A1 (en) | 2005-10-19 | 2006-10-10 | Ultrasonograph for Creating Elastic Image |
EP06811533A EP1938754A4 (en) | 2005-10-19 | 2006-10-10 | ECHOGRAPH FOR CREATING AN ELASTIC IMAGE |
JP2007540932A JP5303147B2 (ja) | 2005-10-19 | 2006-10-10 | 弾性画像を生成する超音波診断装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2005-304347 | 2005-10-19 | ||
JP2005304347 | 2005-10-19 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2007046272A1 true WO2007046272A1 (ja) | 2007-04-26 |
Family
ID=37962368
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2006/320221 WO2007046272A1 (ja) | 2005-10-19 | 2006-10-10 | 弾性画像を生成する超音波診断装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20090143676A1 (ja) |
EP (1) | EP1938754A4 (ja) |
JP (1) | JP5303147B2 (ja) |
CN (1) | CN101291629B (ja) |
WO (1) | WO2007046272A1 (ja) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009063691A1 (ja) * | 2007-11-16 | 2009-05-22 | Hitachi Medical Corporation | 超音波撮像システム |
JP2010022828A (ja) * | 2008-07-16 | 2010-02-04 | Medison Co Ltd | 弾性映像を形成する超音波システム及び方法 |
US20100036250A1 (en) * | 2008-08-05 | 2010-02-11 | Medison Co., Ltd. | Formation Of A Color Map For An Elastic Image |
WO2010044385A1 (ja) * | 2008-10-14 | 2010-04-22 | 株式会社 日立メディコ | 超音波診断装置、及び超音波画像表示方法 |
JP2010124842A (ja) * | 2008-11-25 | 2010-06-10 | Toshiba Corp | 超音波診断装置 |
EP2243431A1 (en) * | 2008-02-07 | 2010-10-27 | Hitachi Medical Corporation | Ultrasonographic device |
WO2011052401A1 (ja) * | 2009-10-30 | 2011-05-05 | 株式会社 日立メディコ | 超音波診断装置、被検体の診断対象部位の疾患の評価用画像生成方法、及び被検体の診断対象部位の疾患の評価用画像生成プログラム |
JP2013146625A (ja) * | 2013-05-08 | 2013-08-01 | Hitachi Medical Corp | 超音波診断装置 |
JP2013188300A (ja) * | 2012-03-13 | 2013-09-26 | Shinji Kume | 超音波診断装置 |
JP2013542046A (ja) * | 2010-11-10 | 2013-11-21 | エコーメトリックス,エルエルシー | 超音波画像処理のシステムおよび方法 |
WO2014002778A1 (ja) * | 2012-06-27 | 2014-01-03 | 日立アロカメディカル株式会社 | 超音波診断装置及び超音波表示方法 |
JP5400773B2 (ja) * | 2008-06-16 | 2014-01-29 | 株式会社日立メディコ | 超音波診断装置、超音波画像表示方法及び超音波診断プログラム |
WO2014103511A1 (ja) | 2012-12-25 | 2014-07-03 | 日立アロカメディカル株式会社 | 超音波診断装置及び弾性解析方法 |
JP2014188149A (ja) * | 2013-03-27 | 2014-10-06 | Seiko Epson Corp | 超音波測定装置、超音波画像装置及び超音波測定方法 |
JP2014193215A (ja) * | 2013-03-28 | 2014-10-09 | Hitachi Aloka Medical Ltd | 超音波撮像装置及び超音波画像表示方法 |
US8948485B2 (en) | 2009-06-10 | 2015-02-03 | Hitachi Medical Corporation | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, ultrasonic image processing program, and ultrasonic image generation method |
JP2015139624A (ja) * | 2014-01-30 | 2015-08-03 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP2015522367A (ja) * | 2012-07-18 | 2015-08-06 | コーニンクレッカ フィリップス エヌ ヴェ | 超音波イメージングデータを処理する方法及びシステム |
WO2015141460A1 (ja) * | 2014-03-19 | 2015-09-24 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP2015186647A (ja) * | 2015-07-16 | 2015-10-29 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP2016515019A (ja) * | 2013-03-15 | 2016-05-26 | セノ メディカル インストルメンツ,インク. | 診断ベクトル分類サポートのためのシステムおよび方法 |
JP2017533773A (ja) * | 2014-11-18 | 2017-11-16 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 組織特性を視覚化する装置 |
US9959622B2 (en) | 2014-07-23 | 2018-05-01 | Samsung Electronics Co. Ltd. | Method and apparatus for supporting diagnosis of region of interest by providing comparison image |
JP2019051192A (ja) * | 2017-09-19 | 2019-04-04 | 株式会社日立製作所 | 超音波医用装置 |
JP2020044044A (ja) * | 2018-09-18 | 2020-03-26 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
Families Citing this family (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008307184A (ja) * | 2007-06-13 | 2008-12-25 | Fujifilm Corp | 画像処理装置および画像処理プログラム |
US9301732B2 (en) * | 2008-03-31 | 2016-04-05 | Hitachi Medical Corporation | Ultrasonic diagnostic arrangements selecting parameters for selecting relevant estimation data for classifying an elasticity image |
WO2011004661A1 (ja) * | 2009-07-07 | 2011-01-13 | 株式会社 日立メディコ | 超音波診断装置及び超音波計測方法 |
US9008382B2 (en) * | 2009-08-03 | 2015-04-14 | Ralph Highnam | Method and system for analysing tissue from images |
JP5665040B2 (ja) * | 2009-09-10 | 2015-02-04 | 学校法人上智学院 | 変位計測方法及び装置、並びに、超音波診断装置 |
WO2011052400A1 (ja) * | 2009-10-28 | 2011-05-05 | 株式会社 日立メディコ | 超音波診断装置及び画像構成方法 |
WO2011087035A1 (ja) | 2010-01-15 | 2011-07-21 | 株式会社 日立メディコ | 超音波診断装置及び超音波画像表示方法 |
WO2012039192A1 (ja) * | 2010-09-21 | 2012-03-29 | 株式会社 日立メディコ | 超音波診断装置および超音波画像の表示方法 |
JP5209026B2 (ja) * | 2010-10-27 | 2013-06-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
JP5209025B2 (ja) * | 2010-10-27 | 2013-06-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
JP2013055567A (ja) * | 2011-09-06 | 2013-03-21 | Olympus Imaging Corp | 撮像装置 |
CN103377475B (zh) * | 2012-04-27 | 2016-08-10 | 通用电气公司 | 血管内斑块的计算机诊断***及方法 |
TWI482956B (zh) * | 2012-07-24 | 2015-05-01 | Univ Nat Cheng Kung | 物體特性感測系統及其控制方法 |
JP5774560B2 (ja) * | 2012-08-20 | 2015-09-09 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置及びその制御プログラム |
JP5710566B2 (ja) * | 2012-09-18 | 2015-04-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置及びその制御プログラム |
KR101512291B1 (ko) * | 2013-05-06 | 2015-04-15 | 삼성메디슨 주식회사 | 의료 영상 장치 및 의료 영상 제공 방법 |
JP2015008733A (ja) * | 2013-06-26 | 2015-01-19 | ソニー株式会社 | 超音波処理装置および方法 |
US20170072227A1 (en) * | 2014-03-28 | 2017-03-16 | Koninklijke Philips N.V., A Corporporation Organized And Existing Under The Laws | Boiling histotripsy methods and systems for uniform volumetric ablation of an object by high-intensity focused ultrasound waves with shocks |
US10646201B2 (en) | 2014-11-18 | 2020-05-12 | C. R. Bard, Inc. | Ultrasound imaging system having automatic image presentation |
EP3220828B1 (en) | 2014-11-18 | 2021-12-22 | C.R. Bard, Inc. | Ultrasound imaging system having automatic image presentation |
KR102325346B1 (ko) * | 2014-12-15 | 2021-11-11 | 삼성전자주식회사 | 의료 영상 진단 장치 및 방법 |
US10580175B2 (en) * | 2015-02-25 | 2020-03-03 | Koninklijke Philips N.V. | Apparatus, method and system for resolution dependent graphical representation of signals |
JP6361528B2 (ja) * | 2015-03-03 | 2018-07-25 | コニカミノルタ株式会社 | 超音波処理装置 |
US20160278741A1 (en) * | 2015-03-24 | 2016-09-29 | Samsung Medison Co., Ltd. | Apparatus and method of measuring elasticity using ultrasound |
US10660604B2 (en) | 2015-07-13 | 2020-05-26 | Otonexus Medical Technologies, Inc. | Apparatus and method for characterization of acute otitis media |
EP3366220A4 (en) | 2015-10-23 | 2019-07-24 | Olympus Corporation | ULTRASONIC OBSERVATION DEVICE, OPERATING METHOD FOR ULTRASONIC OBSERVATION DEVICE AND OPERATING PROGRAM FOR ULTRASONIC OBSERVATION DEVICE |
JP6734079B2 (ja) * | 2016-03-11 | 2020-08-05 | キヤノンメディカルシステムズ株式会社 | 医用診断装置、および医用解析プログラム |
US10675001B2 (en) * | 2016-06-04 | 2020-06-09 | Otonexus Medical Technologies, Inc. | Apparatus and method for characterization of a ductile membrane, surface, and sub-surface properties |
US11138723B2 (en) * | 2016-07-22 | 2021-10-05 | Canon Medical Systems Corporation | Analyzing apparatus and analyzing method |
CN109475343B (zh) * | 2016-08-01 | 2024-04-02 | 深圳迈瑞生物医疗电子股份有限公司 | 剪切波弹性成像测量显示方法及*** |
WO2018116892A1 (ja) * | 2016-12-19 | 2018-06-28 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05317313A (ja) | 1992-05-15 | 1993-12-03 | Ken Ishihara | 超音波診断装置 |
JP2000060853A (ja) | 1998-08-20 | 2000-02-29 | Hitachi Medical Corp | 超音波診断装置 |
JP2004135929A (ja) | 2002-10-18 | 2004-05-13 | Hitachi Medical Corp | 超音波診断装置 |
JP2004261198A (ja) | 2003-01-15 | 2004-09-24 | Hitachi Medical Corp | 超音波診断装置 |
JP2004351062A (ja) * | 2003-05-30 | 2004-12-16 | Hitachi Medical Corp | 超音波診断装置 |
JP2005013283A (ja) | 2003-06-23 | 2005-01-20 | Takeshi Shiina | 超音波探触子及び超音波診断装置 |
WO2005020821A1 (ja) * | 2003-09-01 | 2005-03-10 | Matsushita Electric Industrial Co., Ltd. | 生体信号モニタ装置 |
JP2005066041A (ja) | 2003-08-25 | 2005-03-17 | Hitachi Medical Corp | 超音波探触子及び超音波診断装置 |
JP2005270341A (ja) * | 2004-03-24 | 2005-10-06 | Hitachi Medical Corp | 超音波診断装置 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01254151A (ja) * | 1988-04-04 | 1989-10-11 | Fujitsu Ltd | 超音波診断装置 |
JP3082096B2 (ja) * | 1991-02-08 | 2000-08-28 | 株式会社日立メディコ | 超音波診断装置 |
JP3109747B2 (ja) * | 1991-02-27 | 2000-11-20 | 株式会社日立メディコ | 超音波診断装置 |
JPH10233693A (ja) * | 1996-12-16 | 1998-09-02 | Canon Inc | 画像処理方法及び装置並びに記憶媒体 |
CN1240123A (zh) * | 1998-05-19 | 2000-01-05 | 皇家菲利浦电子有限公司 | 检测弹性变化的方法和用于实施该方法的回波探测装置 |
US6558324B1 (en) * | 2000-11-22 | 2003-05-06 | Siemens Medical Solutions, Inc., Usa | System and method for strain image display |
US7744537B2 (en) * | 2001-08-20 | 2010-06-29 | Japan Science And Technology Agency | Ultrasonic method and system for characterizing arterial tissue using known elasticity data |
DE602004032138D1 (de) * | 2003-06-13 | 2011-05-19 | Panasonic Corp | Ultraschall-diagnosegerät |
WO2005025425A1 (ja) * | 2003-09-12 | 2005-03-24 | Hitachi Medical Corporation | 超音波診断装置 |
CA2457376C (en) * | 2003-10-14 | 2015-09-15 | The University Of British Columbia | Method for imaging the mechanical properties of tissue |
WO2006043529A1 (ja) * | 2004-10-19 | 2006-04-27 | Matsushita Electric Industrial Co., Ltd. | 超音波診断装置 |
-
2006
- 2006-10-10 US US12/083,732 patent/US20090143676A1/en not_active Abandoned
- 2006-10-10 EP EP06811533A patent/EP1938754A4/en not_active Withdrawn
- 2006-10-10 JP JP2007540932A patent/JP5303147B2/ja active Active
- 2006-10-10 CN CN2006800388383A patent/CN101291629B/zh active Active
- 2006-10-10 WO PCT/JP2006/320221 patent/WO2007046272A1/ja active Application Filing
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05317313A (ja) | 1992-05-15 | 1993-12-03 | Ken Ishihara | 超音波診断装置 |
JP2000060853A (ja) | 1998-08-20 | 2000-02-29 | Hitachi Medical Corp | 超音波診断装置 |
JP2004135929A (ja) | 2002-10-18 | 2004-05-13 | Hitachi Medical Corp | 超音波診断装置 |
JP2004261198A (ja) | 2003-01-15 | 2004-09-24 | Hitachi Medical Corp | 超音波診断装置 |
JP2004351062A (ja) * | 2003-05-30 | 2004-12-16 | Hitachi Medical Corp | 超音波診断装置 |
JP2005013283A (ja) | 2003-06-23 | 2005-01-20 | Takeshi Shiina | 超音波探触子及び超音波診断装置 |
JP2005066041A (ja) | 2003-08-25 | 2005-03-17 | Hitachi Medical Corp | 超音波探触子及び超音波診断装置 |
WO2005020821A1 (ja) * | 2003-09-01 | 2005-03-10 | Matsushita Electric Industrial Co., Ltd. | 生体信号モニタ装置 |
JP2005270341A (ja) * | 2004-03-24 | 2005-10-06 | Hitachi Medical Corp | 超音波診断装置 |
Non-Patent Citations (3)
Title |
---|
KANEI H. ET AL.: "Kekkanheki no Danseiritsu Dansozo Imaging", EIZO JOHO MEDICAL, vol. 37, no. 9, 2005, pages 904 - 909, XP003012122 * |
See also references of EP1938754A4 * |
UENO E ET AL.: "Dynamic tests in real-time breast echography", ULTRASOUND MED BIOL., vol. 14, 1988, pages 53 - 57, XP023258532, DOI: doi:10.1016/0301-5629(88)90047-6 |
Cited By (39)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009063691A1 (ja) * | 2007-11-16 | 2009-05-22 | Hitachi Medical Corporation | 超音波撮像システム |
EP2243431A4 (en) * | 2008-02-07 | 2013-10-23 | Hitachi Medical Corp | ULTRASOUND DEVICE |
EP2243431A1 (en) * | 2008-02-07 | 2010-10-27 | Hitachi Medical Corporation | Ultrasonographic device |
JP5400773B2 (ja) * | 2008-06-16 | 2014-01-29 | 株式会社日立メディコ | 超音波診断装置、超音波画像表示方法及び超音波診断プログラム |
JP2010022828A (ja) * | 2008-07-16 | 2010-02-04 | Medison Co Ltd | 弾性映像を形成する超音波システム及び方法 |
US20100036250A1 (en) * | 2008-08-05 | 2010-02-11 | Medison Co., Ltd. | Formation Of A Color Map For An Elastic Image |
JP2010036041A (ja) * | 2008-08-05 | 2010-02-18 | Medison Co Ltd | 弾性情報のカラーマップを形成する超音波システム及びカラーマップ形成方法 |
US9629605B2 (en) * | 2008-08-05 | 2017-04-25 | Samsung Medison Co., Ltd. | Formation of a color map for an elastic image |
US20110194748A1 (en) * | 2008-10-14 | 2011-08-11 | Akiko Tonomura | Ultrasonic diagnostic apparatus and ultrasonic image display method |
WO2010044385A1 (ja) * | 2008-10-14 | 2010-04-22 | 株式会社 日立メディコ | 超音波診断装置、及び超音波画像表示方法 |
JP5479353B2 (ja) * | 2008-10-14 | 2014-04-23 | 株式会社日立メディコ | 超音波診断装置 |
JP2010124842A (ja) * | 2008-11-25 | 2010-06-10 | Toshiba Corp | 超音波診断装置 |
JP5730196B2 (ja) * | 2009-06-10 | 2015-06-03 | 株式会社日立メディコ | 超音波診断装置、超音波画像処理装置、超音波画像生成方法 |
US8948485B2 (en) | 2009-06-10 | 2015-02-03 | Hitachi Medical Corporation | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, ultrasonic image processing program, and ultrasonic image generation method |
JP5560283B2 (ja) * | 2009-10-30 | 2014-07-23 | 株式会社日立メディコ | 超音波診断装置、被検体の診断対象部位の疾患の評価用画像生成方法、及び被検体の診断対象部位の疾患の評価用画像生成プログラム |
WO2011052401A1 (ja) * | 2009-10-30 | 2011-05-05 | 株式会社 日立メディコ | 超音波診断装置、被検体の診断対象部位の疾患の評価用画像生成方法、及び被検体の診断対象部位の疾患の評価用画像生成プログラム |
JP2013542046A (ja) * | 2010-11-10 | 2013-11-21 | エコーメトリックス,エルエルシー | 超音波画像処理のシステムおよび方法 |
JP2013188300A (ja) * | 2012-03-13 | 2013-09-26 | Shinji Kume | 超音波診断装置 |
WO2014002778A1 (ja) * | 2012-06-27 | 2014-01-03 | 日立アロカメディカル株式会社 | 超音波診断装置及び超音波表示方法 |
JP2015522367A (ja) * | 2012-07-18 | 2015-08-06 | コーニンクレッカ フィリップス エヌ ヴェ | 超音波イメージングデータを処理する方法及びシステム |
WO2014103511A1 (ja) | 2012-12-25 | 2014-07-03 | 日立アロカメディカル株式会社 | 超音波診断装置及び弾性解析方法 |
JP2016515019A (ja) * | 2013-03-15 | 2016-05-26 | セノ メディカル インストルメンツ,インク. | 診断ベクトル分類サポートのためのシステムおよび方法 |
US10949967B2 (en) | 2013-03-15 | 2021-03-16 | Seno Medical Instruments, Inc. | System and method for diagnostic vector classification support |
JP2019088838A (ja) * | 2013-03-15 | 2019-06-13 | セノ メディカル インストルメンツ,インク. | 診断ベクトル分類サポートのためのシステムおよび方法 |
JP2014188149A (ja) * | 2013-03-27 | 2014-10-06 | Seiko Epson Corp | 超音波測定装置、超音波画像装置及び超音波測定方法 |
JP2014193215A (ja) * | 2013-03-28 | 2014-10-09 | Hitachi Aloka Medical Ltd | 超音波撮像装置及び超音波画像表示方法 |
JP2013146625A (ja) * | 2013-05-08 | 2013-08-01 | Hitachi Medical Corp | 超音波診断装置 |
WO2015114887A1 (ja) * | 2014-01-30 | 2015-08-06 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP2015139624A (ja) * | 2014-01-30 | 2015-08-03 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP2015177883A (ja) * | 2014-03-19 | 2015-10-08 | 日立アロカメディカル株式会社 | 超音波診断装置 |
WO2015141460A1 (ja) * | 2014-03-19 | 2015-09-24 | 日立アロカメディカル株式会社 | 超音波診断装置 |
US9959622B2 (en) | 2014-07-23 | 2018-05-01 | Samsung Electronics Co. Ltd. | Method and apparatus for supporting diagnosis of region of interest by providing comparison image |
JP2017533773A (ja) * | 2014-11-18 | 2017-11-16 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 組織特性を視覚化する装置 |
JP2022125071A (ja) * | 2014-11-18 | 2022-08-26 | コーニンクレッカ フィリップス エヌ ヴェ | 組織特性を視覚化する装置 |
JP7504158B2 (ja) | 2014-11-18 | 2024-06-21 | コーニンクレッカ フィリップス エヌ ヴェ | 組織特性を視覚化する装置 |
JP2015186647A (ja) * | 2015-07-16 | 2015-10-29 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP2019051192A (ja) * | 2017-09-19 | 2019-04-04 | 株式会社日立製作所 | 超音波医用装置 |
JP2020044044A (ja) * | 2018-09-18 | 2020-03-26 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
JP7162477B2 (ja) | 2018-09-18 | 2022-10-28 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
Also Published As
Publication number | Publication date |
---|---|
US20090143676A1 (en) | 2009-06-04 |
EP1938754A4 (en) | 2011-01-05 |
CN101291629B (zh) | 2010-12-01 |
JP5303147B2 (ja) | 2013-10-02 |
JPWO2007046272A1 (ja) | 2009-04-23 |
EP1938754A1 (en) | 2008-07-02 |
CN101291629A (zh) | 2008-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5303147B2 (ja) | 弾性画像を生成する超音波診断装置 | |
JPWO2007046272A6 (ja) | 弾性画像を生成する超音波診断装置 | |
JP5334413B2 (ja) | 超音波診断装置 | |
JP4919972B2 (ja) | 弾性画像表示方法及び弾性画像表示装置 | |
JP5160227B2 (ja) | 超音波診断装置及び超音波画像表示方法 | |
JP5925438B2 (ja) | 超音波診断装置 | |
JP4966578B2 (ja) | 弾性画像生成方法及び超音波診断装置 | |
JP5028416B2 (ja) | 超音波診断装置 | |
JP5560283B2 (ja) | 超音波診断装置、被検体の診断対象部位の疾患の評価用画像生成方法、及び被検体の診断対象部位の疾患の評価用画像生成プログラム | |
JP6608232B2 (ja) | 医用画像診断装置、医用画像処理装置および医用情報の表示制御方法 | |
JP5285616B2 (ja) | 超音波診断装置とその作動方法及び超音波画像診断プログラム | |
US8118746B2 (en) | Ultrasonic diagnostic apparatus | |
JP5087341B2 (ja) | 超音波診断装置 | |
WO2019131811A1 (ja) | 解析装置 | |
JP2008259605A (ja) | 超音波診断装置 | |
JP6352013B2 (ja) | 超音波診断装置、画像処理装置及び画像処理プログラム | |
JP5113322B2 (ja) | 超音波診断装置 | |
JP5415669B2 (ja) | 超音波診断装置 | |
JP2009261493A (ja) | 超音波診断方法および装置 | |
JP7171228B2 (ja) | 超音波診断装置及び医用情報処理プログラム | |
JP5638641B2 (ja) | 超音波診断装置 | |
JP5663640B2 (ja) | 超音波診断装置 | |
JP2012055742A (ja) | 超音波診断装置 | |
JP7387482B2 (ja) | 医用情報処理装置 | |
JP2013183982A (ja) | 超音波診断装置及び弾性画像生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 200680038838.3 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
ENP | Entry into the national phase |
Ref document number: 2007540932 Country of ref document: JP Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2006811533 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 12083732 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |