WO2005025425A1 - 超音波診断装置 - Google Patents

超音波診断装置 Download PDF

Info

Publication number
WO2005025425A1
WO2005025425A1 PCT/JP2004/013140 JP2004013140W WO2005025425A1 WO 2005025425 A1 WO2005025425 A1 WO 2005025425A1 JP 2004013140 W JP2004013140 W JP 2004013140W WO 2005025425 A1 WO2005025425 A1 WO 2005025425A1
Authority
WO
WIPO (PCT)
Prior art keywords
area
image
elastic
ultrasonic diagnostic
scoring
Prior art date
Application number
PCT/JP2004/013140
Other languages
English (en)
French (fr)
Inventor
Takeshi Matsumura
Ryuichi Shinomura
Tsuyoshi Mitake
Mitsuhiro Oshiki
Original Assignee
Hitachi Medical Corporation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hitachi Medical Corporation filed Critical Hitachi Medical Corporation
Priority to US10/571,578 priority Critical patent/US8118746B2/en
Priority to EP04787803.8A priority patent/EP1665987B1/en
Priority to JP2005513893A priority patent/JP4685633B2/ja
Publication of WO2005025425A1 publication Critical patent/WO2005025425A1/ja

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S15/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8977Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/5206Two-dimensional coordinated display of distance and direction; B-scan display
    • G01S7/52063Sector scan display
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52071Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging

Definitions

  • the present invention relates to an ultrasonic diagnostic apparatus capable of displaying a diagnostic site in a subject using an ultrasonic wave as an elastic image representing the hardness or softness of a living tissue.
  • a conventional general ultrasonic diagnostic apparatus includes an ultrasonic transmission / reception control unit for controlling ultrasonic transmission / reception, an ultrasonic transmission / reception unit for transmitting / receiving ultrasonic waves to / from a subject, and an ultrasonic transmission / reception unit.
  • Tomographic scanning means for repeatedly obtaining tomographic image data in a subject including a moving tissue using a reflected echo signal at predetermined intervals, and image display means for displaying time-series tomographic image data obtained by the tomographic scanning means. It was constituted having. Then, the structure of the living tissue inside the subject was displayed, for example, as a B-mode tomographic image.
  • Patent Document 1 JP-A-5-317313
  • Patent Document 2 JP-A-2000-060853
  • the present invention has been made in view of the above-mentioned points, and is useful for diagnosis of a diseased part. It is an object of the present invention to provide an ultrasonic diagnostic apparatus capable of giving a certain index.
  • an ultrasonic diagnostic apparatus of the present invention causes an ultrasonic wave to be transmitted from a probe to a subject, receives a reflected echo signal corresponding to the transmission of the ultrasonic wave, and receives the echo signal.
  • a displacement measurement unit configured to measure a displacement of the living tissue of the subject based on a projection echo signal; an elasticity image forming unit configured to determine an amount of distortion or an elastic modulus of the displacement power to form an elasticity image; and display the elasticity image.
  • a scoring means for specifying a distortion state or an elastic state based on information of the elasticity image output from the elasticity image forming unit.
  • the scoring means specifies a distortion state or an elastic state based on the area or shape of the elastic image.
  • a tomographic image forming unit configured to form a tomographic image of the reflected echo signal; and the scoring unit determines a distortion state based on a relationship between a specific region on the tomographic image and a region of a hard part of the elastic image.
  • the elastic state is specified.
  • FIG. 1 is a block diagram showing an embodiment of an ultrasonic diagnostic apparatus according to the present invention.
  • FIG. 2 is a diagram showing an example of a method of attaching a pressure measuring unit (pressure sensor) to an ultrasonic probe and measuring a pressure between a head of the probe and a test object.
  • a pressure measuring unit pressure sensor
  • FIG. 3 is a diagram showing an example of a scoring process executed by the score input means of FIG. 1.
  • FIG. 4 is a diagram showing a display example of an image display when a scoring process is performed using a score input unit.
  • FIG. 5 is a diagram showing scoring by a stress signal.
  • FIG. 6 is a diagram showing scoring based on a change in a region of an elastic image due to compression.
  • FIG. 7 is a diagram showing scoring of an elasticity image based on a hue pattern.
  • FIG. 8 is a view showing another embodiment of the scoring process executed by the score input means of FIG. 1.
  • FIG. 9 is a diagram showing still another embodiment of the scoring process executed by the score input means of FIG. 1. It is.
  • FIG. 10 is a flowchart illustrating an example of an automatic scoring process using an area extraction application.
  • FIG. 11 is a diagram schematically showing an example of an operation of the automatic scoring process.
  • FIG. 12 is a diagram showing a timing of acquiring an elasticity image for scoring.
  • FIG. 13 is a diagram showing one embodiment of an RF signal frame data selection unit in FIG. 1.
  • FIG. 14 is a diagram showing another example of the RF signal frame data selection unit in FIG. 1.
  • FIG. 1 is a block diagram showing an embodiment of an ultrasonic diagnostic apparatus according to the present invention.
  • This ultrasonic diagnostic apparatus obtains a tomographic image of a diagnostic site of the subject 100 using ultrasonic waves, and displays an elastic image representing the hardness or softness of a living tissue.
  • this ultrasonic diagnostic apparatus uses an ultrasonic probe 101, a transmitting circuit 102, a receiving circuit 103, a phasing addition circuit 104, and a signal processing unit 105.
  • An adder 114, a score input means 115 for inputting a score, and a state of distortion from an elasticity image or an elasticity image and a tomographic image! ⁇ is provided with a scoring unit 116 for specifying the elastic state.
  • the ultrasonic probe 101, the transmission circuit 102, the reception circuit 103, the phasing addition circuit 104, and the signal processing unit 105 constitute an ultrasonic transmission / reception unit.
  • the ultrasonic transmitting and receiving means obtains one tomographic image by scanning the ultrasonic beam in a fixed direction in the body of the subject 100 using the ultrasonic probe 101.
  • the ultrasonic probe 101 is formed by arranging a large number of transducers in a strip shape, and transmits and receives ultrasonic waves to the subject 100 by performing beam scanning mechanically or electronically. Although not shown in the figure, a vibrator that generates a supersonic wave and receives a reflected echo is incorporated therein.
  • Each transducer In general, a function of converting an input pulse wave or continuous wave transmission signal into an ultrasonic wave and emitting the ultrasonic wave, and receiving an ultrasonic wave reflected from the inside of the subject 100 into a reception signal of an electric signal. It is formed with the function of converting and outputting.
  • the transmission circuit 102 generates a transmission pulse for driving the ultrasonic probe 101 to generate an ultrasonic wave, and transmits an ultrasonic wave transmitted by a built-in transmission phasing addition circuit.
  • the convergence point of is set to a certain depth.
  • the receiving circuit 103 amplifies the reflected echo signal received by the ultrasonic probe 101 with a predetermined gain.
  • the number of received signals corresponding to the number of the amplified vibrators are input to the phasing addition circuit 104 as independent received signals.
  • the phasing addition circuit 104 inputs the received signal amplified by the receiving circuit 103, controls the phase thereof, and forms an ultrasonic beam at one or a plurality of convergence points.
  • the signal processing unit 105 receives the received signal from the phasing addition circuit 104 and performs various signal processing such as gain correction, log compression, detection, contour enhancement, and filter processing.
  • the black-and-white scan converter 106 uses the reflected echo signal output from the signal processing unit 105 of the ultrasonic transmission / reception means to generate RF signal frame data in the subject 100 including the moving tissue using the reflected echo signal.
  • the RF signal frame data is acquired at a sound wave cycle, and displayed on the image display 107 via the switching adder 114. Therefore, the black-and-white scan converter 106 is a means for controlling the tomographic scanning means and the system for sequentially reading out the RF signal frame data at the cycle of the television system, for example, the reflected echo signal from the signal processing unit 105 is digitally converted. It includes an A / D converter that converts it into a signal, a plurality of frame memories that store the tomographic image data digitized by the A / D converter in time series, and a controller that controls these operations. Be composed.
  • the image display 107 displays the time-series tomographic image data obtained by the black-and-white scan converter 106, that is, the B-mode tomographic image, and is output from the black-and-white scan converter 106 via the switching adder 114. It consists of a D / A converter that converts image data into analog signals, and a color TV monitor that inputs analog video signals from the D / A converter and displays them as images.
  • the output-side force of the phasing addition circuit 104 is branched, and an RF signal frame data selection unit 108 and a displacement measurement unit 109 are provided.
  • a pressure measuring unit 110 is provided in parallel with this.
  • a strain and elasticity calculating unit 111 At the subsequent stage of the displacement measuring unit 109 and the pressure measuring unit 110, there are provided a strain and elasticity calculating unit 111, an elasticity data processing unit 112, and a color scan converter 113, and the output of the color scan converter 113 and the monochrome scan converter 106 is provided.
  • a switching adder 114 is provided on the side. Then, the outputs of the black-and-white scan converter 106 and the color scan converter 113 are introduced into a scoring unit 116, where a scoring process is performed. The details of this scoring process will be described later.
  • the displacement measuring unit 109 performs one-dimensional or two-dimensional correlation processing based on the set of RF signal frame data selected by the RF signal frame data selecting unit 108, and performs displacement or displacement of each point on the tomographic image. It measures the movement vector (direction and magnitude of displacement).
  • a method of detecting the movement vector for example, there are a block matching method and a gradient method as described in Patent Document 1.
  • the block matching method divides an image into blocks composed of, for example, NXN pixels, searches for the block closest to the block of interest in the current frame in the previous frame, and performs predictive coding based on this. It is.
  • the pressure measurement unit 110 measures or estimates the pressure in the body cavity of the diagnostic site of the subject 100.
  • This ultrasonic diagnostic apparatus is installed in the probe head 1011 while transmitting and receiving ultrasonic waves under the control of the control unit 200 using the ultrasonic probe 101 provided in the probe head 1011.
  • a pressure is applied or reduced by a pressurizer (not shown) to apply a stress distribution to the body cavity of the diagnostic site of the subject 100.
  • the pressure applied to the rod-like member is detected.
  • the pressure sensor 1012 is attached to the side of the probe head 1011 and the pressure between the probe head 1011 and the subject 100 is measured at an arbitrary time phase. To be sent to In FIG. 2, a pressurizer that pressurizes and depressurizes the probe head 1011 is omitted.
  • the strain and elastic modulus calculating unit 111 calculates the strain or the elastic modulus of each point on the tomographic image based on the displacement (displacement) and the pressure output from the displacement measuring unit 109 and the pressure measuring unit 110, respectively. To generate numerical data of strain or elastic modulus (elastic frame data). This is output to the polarity data processing unit 112.
  • the calculation of strain performed by the strain and elastic modulus calculation unit 111 does not require pressure data, for example, and is obtained by calculation by spatially differentiating the displacement.
  • the calculation of the Young's modulus, one of the elastic moduli is obtained by calculation by dividing the change in pressure by the change in the amount of movement.
  • the color scan converter 113 converts the elasticity frame data output from the elasticity data processing unit 112 with the command output from the control unit 200 of the ultrasonic diagnostic apparatus or the elasticity frame data output from the elasticity data processing unit 112.
  • a hue information conversion unit for inputting an upper limit value and a lower limit value as a middle gradation selection range, and providing hue information such as red, green, and blue as elastic image data based on the elastic frame data; You.
  • This hue information conversion means converts, for example, a region where large distortion is measured in the elasticity frame data output from the elasticity data processing unit 112 into a red code in the elasticity image data, and conversely, measures a small distortion. It operates to convert the set area into a blue code in the elasticity image data.
  • the color scan converter 113 may be constituted by the above-described monochrome scan converter 106.
  • the brightness may be increased in the elasticity image data in an area where the distortion is measured to be large, and may be decreased in the elasticity image data in the area where the distortion is measured to be small.
  • the switch adder 114 is a means for inputting the black-and-white tomographic image data from the black-and-white scan converter 106 and the color elastic image data from the color scan converter 113 and adding or switching both images. The output is switched so that only black-and-white tomographic image data or only color elastic image data is output, or both image data are added and synthesized. Further, for example, as described in Patent Document 2, a monochrome tomographic image and a monochrome elastic image by a color or monochrome scan converter may be simultaneously displayed on a two-screen display. The image data output from the switching adder 114 is output to the image display 107.
  • FIG. 3 is a diagram showing a first embodiment of the scoring process performed by the scoring unit 116 of FIG. 1.
  • FIG. 3 (A) shows a monochrome tomographic image output from the monochrome scan converter 106.
  • FIG. 4 is a diagram illustrating an example of data (B-mode tomographic image).
  • Figure 3 (B)-(F) shows the color scan
  • FIG. 9 is a diagram showing an example of color elasticity image data output from the data 113, and schematically shows a case where each symptom is different. These elasticity image data actually perform color display according to the hardness.
  • a soft part (a part with a large distortion) is displayed as red
  • a hard part (a part with a small distortion) is displayed as blue
  • the middle between the two is displayed as a continuous color change.
  • Fig. 3 it is not possible to display in such colors, so that the state is indicated by shading, and the soft areas are displayed by rough shading, the hard areas are displayed by fineness, and the areas are displayed by shading.
  • the central elliptical portion 61 is a disease site such as a mammary tumor, and is called a low echo area.
  • the scoring process is performed by operating the elasticity image data score input means 115 in accordance with how the elasticity image data is displayed in the low echo area 61 according to the state. Each sex image is classified into a score of 1 to 5.
  • the scores 115 are based on the following criteria. The observer determines whether or not the power satisfies the criteria, and performs a scoring process. Hereinafter, the criteria for each of the scores 115 will be described.
  • Score 2 As shown in Fig. 3 (C), cases where hard regions 63 and 64 based on the elasticity image are partially recognized, avoiding the center of the low echo region 61.
  • Score 4 As shown in Fig. 3 (E), a case in which the image is recognized evenly distributed up to the boundary (contour) of the hard region 66 force low echo region 61 based on the elasticity image
  • the hard region is a certain threshold value of distortion in a distortion image as an elastic image. For example, if a strain of 0% is set as the threshold value Ts, the hard region is a region in which no compression is achieved even by compression from the body surface. (Measurement point cloud). In the case where the pressing force is captured during scanning of the diagnostic region and the elasticity image is used as the elasticity image, the region where the elasticity is measured is larger than a certain elasticity threshold Ty. do it.
  • FIG. 4 is a diagram showing a display example of the image display 107 when performing a scoring process using the score input means 115.
  • a translucent color image of elastic image data showing the distribution of tissue hardness is superimposed on a B-mode tomographic image as shown in Fig. 3 (B)-(F).
  • a column 1072 for inputting a scoring index is provided on the right side of the column, and a display area 1073 for displaying the index of each score is provided.
  • the contents of the respective scores as described above are simply displayed, or the figures such as FIGS. 3 (B) to (F) are also displayed.
  • the examiner compares the scoring index 1073 with the scoring index human power column 1072 and inputs a score that is considered appropriate by looking at the display screen 1071. These operations are executed by the score input unit 115.
  • a B-mode tomographic image may be displayed in parallel in another window as shown in FIG. This makes it possible to improve the visibility during the scoring process.
  • the score entered in the scoring index input field 1072 is also linked to the individual patient reporting function provided in the ultrasound system, so that it can be entered simultaneously in the scoring index input field provided in the reporting function. It may be. Further, a model diagram serving as an index of scoring as shown in FIGS. 3 (B)-(F) may be added to the display area 1073.
  • the scoring unit 116 automatically scores using image processing by the force scoring unit 116 described in the case where the examiner performs scoring processing while viewing the image displayed on the display screen 1071 ( Classify.
  • the scoring unit 116 performs edge detection using well-known edge detection software on the color elastic image data (distorted image) output from the color scan converter 113, and determines the area B of the hard part area 62-67 as a pixel. Calculate by number.
  • scoring is performed using only elasticity image data. For example, an elasticity image serving as a reference image is obtained in advance as reference information, and the number of pixels of the elasticity image area is set as an area B, and a score based on the area of the elasticity image is set as follows.
  • scoring is performed by applying the area of the newly acquired elasticity image data to the area B set as described above. If these elastic images are connected in series, set the minimum score to 2 or more. If the images are connected by a smooth circular system, set the minimum score to 3 or more. The score may be determined based on the shape of the image.
  • FIG. 5 shows a scoring process using the characteristic that a stress signal is output on an elastic image by pressing a subject.
  • FIGS. 5A and 5B are views showing a state where the subject 100 is pressed by the probe 101.
  • FIG. 5 (b) shows a malignant tumor.
  • FIG. 5 (a) shows a benign tumor.
  • scoring section 116 a score corresponding to width 123 of stress signal 122 is set in advance. For example, the score is set as shown in FIG. 4 such that the elastic image is more malignant as the width is larger. Then, by comparing the width set corresponding to the score with the width 123 of the obtained stress signal 122 of the elastic image 121, the elastic image 121 is scored.
  • scoring may be performed using the area of the force stress signal 122 in which scoring is performed using the width of the stress signal 122. In this way, scoring can be performed using the stress signal of the elasticity image.
  • FIG. 6 shows an example in which a scoring process is performed using an elasticity image that changes when a subject is pressed.
  • a malignant tumor there is a characteristic that, when the subject is gradually pressed, for example, the area displayed in blue gradually becomes smaller.
  • a benign tumor there is a characteristic that when the subject is gradually pressed, the area displayed in blue gradually increases. So this Using the characteristic, the score for the characteristic corresponding to the color change of the region caused by the compression is set in the front scoring unit 116 as shown in FIG. Then, scoring is performed according to the color of the elastic image that changes with the compression. Thus, scoring can be performed based on a change in the area of the elastic image due to the compression.
  • FIG. 7 shows an example of performing scoring processing using a cyst pattern.
  • the cyst pattern is a pattern in which a portion 141 of a specific region of the displayed elasticity image with a small depth is displayed in blue, and a region 142 with a large depth is displayed in red in a pair. It should be noted that the area displayed by this cyst pattern is empirically known to be a benign tumor. Therefore, the cyst pattern is separately set as a benign score in the scoring unit 116, and if the cyst pattern is recognized on the elasticity image, it is scored as benign. In this manner, scoring can be performed based on the hue pattern of the elastic image.
  • the scoring unit 116 detects the low echo area 61 of the diseased part in the black-and-white tomographic image data (B-mode tomographic image) output from the black-and-white scan converter 106 using known edge detection software.
  • the area A of the detected area 61 is obtained by the number of pixels or the like.
  • the scoring unit 116 performs similar edge detection on the color elasticity image data (distortion image) output from the color scan converter 113, and performs the same edge detection in the low echo area 61 of the diseased part.
  • the area B of the hard part area 62-67 to be measured is obtained by the number of pixels or the like.
  • the sum of the areas is defined as the area B.
  • the scoring unit 116 performs the following scoring processing based on the obtained relationship between the areas A and B.
  • the ratio of the areas A and B is calculated as B / A force 3 ⁇ 4. Also, in the supersonic ultrasound diagnostic apparatus, the threshold of the area ratio is set as
  • Score 1 threshold Thl (eg 0.1)
  • Score 2 threshold Th2 (for example, 0.3)
  • Score 3 threshold Th3 (eg 0.7)
  • Score 4 threshold Th4 (for example, 1.0)
  • a central portion in the low echo area 61 has a certain ratio smaller than the area of the low echo area (for example, Area AX 0.6)
  • a region of interest 61a is provided, and for the region of interest 61a, in the same manner as described above, the hardness of the elasticity image, the extent to which the region occupies, and the evaluation are compared to the scoring method described above. This can be used to evaluate the scoring of the spatial distribution of hard regions inside and outside the low echo region, not just the ratio of the occupied area of the hard region to the low echo region. It can be expected that the result will increase the accuracy of the evaluation.
  • a plurality of regions of interest 62b are provided inside or outside the boundary with reference to the boundary of the low echo region as well as the central region of interest.
  • the degree of occupancy of the hard region in the elasticity image may be evaluated for the plurality of regions of interest, and the area occupied by each of the plurality of regions of interest 62b may be evaluated.
  • Edge detection can be achieved by using a known snake or region growing method. If the edge is hidden by shadows, it can be handled by a method such as setting the ROI or interpolating surrounding edge information.
  • scoring is performed using the occupancy of the hard region distributed in the low echo region.
  • Statistical processing is performed using a group of gender measurement points as a population, and scoring may be performed based on the statistical features. For example, a method using the average value will be described as follows.
  • N the total number of elements in the measurement point group in the low echo area
  • the scoring unit 116 calculates the area A of the low-coordinate area 61 of the diseased site in the black-and-white tomographic image data (B-mode tomographic image) output from the black-and-white scan converter 106. It is detected by known region extraction software, for example, the region growing method. The area A of the detected area is obtained by the number of pixels or the like.
  • the scoring unit 116 displays the color elastic image data (distorted image) output from the color scan converter 113 as a hard region 62-67 within the area of the area A. Count the number of pixels that can be arbitrarily determined and have a predetermined luminance value or more. The number of pixels is defined as area B. The following scoring process is performed based on the area A and the area B (the number of pixels having a predetermined luminance value or more).
  • the ratio of the area (number of pixels) A, B is calculated as B / A force 3 ⁇ 4.
  • the area ratio threshold is set in advance.
  • Score 1 threshold Thl (eg 0.1)
  • Score 2 threshold Th2 (for example, 0.3)
  • Score 3 threshold Th3 (eg 0.7)
  • Score 4 threshold Th4 (for example, 1.0)
  • a hard region 65 ⁇ 3 is obtained by a known region extraction application (eg, a region growing method). An area indicated as 67 and having a predetermined luminance or more is extracted. this At this time, the extraction region is also a region to be extracted that includes the edge that does not remain in the low echo region 61 described above.
  • the obtained extraction area is defined as C.
  • C / A is calculated as Z 'as the ratio of area A to area C. This value of Z '
  • scoring can be performed up to a score of 3 without using a region extraction process for elasticity image data.
  • FIG. 10 is a flowchart illustrating an example of an automatic scoring process using the region extraction application.
  • FIG. U is a diagram schematically showing an example of the operation of the automatic scoring process.
  • the region extraction algorithm (Region Growing) will be described as an example.
  • step S91 in order to obtain a strain elasticity image, a low echo portion whose region is to be extracted is selected from a B-mode tomographic image which is drawn after being pressed by an arbitrary means.
  • This selection processing is performed by designating an arbitrary position in the low echo area 61 as a source point (denoted as SP in the figure) as shown in FIG.
  • a low echo area 61 as shown in FIG. 11B is extracted based on the source point SP specified in step S91 and a threshold given by any means.
  • the region growing method is a method of extracting an area adjacent to the designated source point SP, where the luminance value difference falls within a preset value (threshold). Therefore, the area of the extracted low echo area 61 is counted by means such as counting the number of pixels. In FIG. 11B, this area is represented as SBW.
  • step S93 in the low echo area, an area having a hardness value and a luminance value representing the region that is harder than the luminance value representing the region arbitrarily determined by the user in the drawn elasticity image.
  • this hard region is denoted as Sbl.
  • step S94 the area Sbl of the hard region is smaller than the area SBWX 0.7 of the low echo region. It is determined whether or not it is. If the result of the determination is that the area Sbl of the hard area is smaller than the area SBWX 0.7 of the low echo area! / (SBWX 0.7> Sbl), the process proceeds to step S95. Then, in step S95, values up to a score of 13 corresponding to each value are returned.
  • step S94 Conversely, as a result of the determination in step S94, the area Sbl of the hard region is changed to the area SBW of the low echo region.
  • step S96 If X is 0.7 or more (SBWX 0.7 ⁇ Sbl), the process proceeds to step S96 because there is a possibility that the area Sbl of the hard part is not within the low echo area.
  • a region (area Sbl ′) is extracted again from the elasticity image data by the region growing method. In this case, extraction is performed up to the area including the area around the low echo area of the B-mode tomographic image.
  • the source point SP may be the one specified in the B-mode tomographic image as it is, or may be newly set separately. Let the area of the region extracted here be Sbl '.
  • values of scores 4 and 5 corresponding to the ratio (SBWZSbl ') between the area SBW of the low echo area and the area Sbl' are returned. Specifically, if Sbl 'is SBW, a score of 4 will be returned; otherwise, a score of 5 will be returned.
  • the ultrasonic Using the interface the examiner inputs the outline of the diseased part on the displayed B-mode tomogram, obtains the area A of the area 61 based on the outline information, and uses the method described above.
  • the scoring may be automatically determined in comparison with the area B according to. Note that the above scoring is an example, and it is important to score using an elastic image including a strain image.
  • the numerical values used in the above scoring are merely examples, and actual clinical It is needless to say that it is desirable to use various optimum values in accordance with the above.
  • the scoring method particularly for the mammary gland region is used as an example. Needless to say, it defines the coring method.
  • the scoring method according to the above-described embodiments may be used not only by itself for scoring but also by combining a plurality of methods and performing one scoring.
  • the scoring unit 116 outputs black-and-white tomographic image data (B-mode tomographic image) output from the black-and-white scan converter 106 and color elastic image data (distorted image data) output from the color scan converter 113. ), The same scoring process is performed using the output from the signal processing unit 105, the output from the elasticity data processing unit 112, or the output from the distortion and elastic modulus calculation unit 111. You may do it. Further, in the above-described embodiment, when the force distortion in which the pressure measuring unit is shown in FIG. 1 is obtained and a distorted image is displayed, the pressure measuring unit may not be provided.
  • FIG. 12 is a diagram showing the timing at which the distortion or elastic modulus is calculated by the distortion and elastic modulus calculator 111 of the above-described embodiment.
  • the state is indicated by 0, and the limit of the pressing is indicated by MAX.
  • the timing of acquiring the elasticity image is the timing of the light compression state, and the elasticity image in the compressed state of about 3 to 20% of the compression limit is acquired, and the acquired elasticity image is scored.
  • the pressure sensor 1012 detects a pressure value at the limit of compression in advance, and stores the detected pressure value in the scoring unit 116 as a pressure maximum value.
  • the pressure value of the elasticity image obtained by pressing is set using the score input means 115.
  • the pressure value to be set may be arbitrarily determined as, for example, a pressure value of 10% of the maximum pressure value. Then, the subject 100 is pressed, an elasticity image at the time when the set pressure value is reached is obtained, and scoring is performed on the obtained elasticity image.
  • An elasticity image in a state where the subject 100 is strongly pressed may be hard or hard to determine the position and softness of the object 100.
  • the scoring accuracy can be improved by ringing.
  • FIG. 13 is a diagram showing an embodiment of the RF signal frame data selection unit of FIG. It is.
  • the RF signal frame data selection unit 108 arbitrarily selects the number of frames going back in the past (the number of frame intervals with the current frame data) as one RF signal frame data serving as a reference for displacement measurement. That is, the RF signal frame data selection unit 108 includes the RF signal frame data output one after another from the phasing addition circuit 104 one after another at the frame rate of the ultrasonic diagnostic apparatus. In the frame memory 1081.
  • the RF signal frame data selection unit 108 sets the data currently secured in the frame memory 1081 as the RF signal frame data N.
  • the RF signal frame data selection unit 108 temporally stores the past RF signal frame data N-1, N-2, N-3, N, M in accordance with a control command from the control unit 200 of the ultrasonic diagnostic apparatus. , One RF signal frame data is selected, and it is temporarily stored in the RF signal frame data selection circuit 1082 as RF signal frame data X.
  • the RF signal frame data selection unit 108 stores the latest RF signal frame data N stored in the frame memory 1081 and the RF signal frame data X stored in the RF signal frame data selection circuit 1082 in parallel. Then, it outputs to the displacement measuring unit 109.
  • FIG. 14 is a diagram illustrating another example of the RF signal frame data selection unit in FIG.
  • the RF signal frame data selection unit 108 in FIG. 14 secures the RF signal frame data P acquired in a certain past phase P in the frame memory 1081 according to a control instruction from the control unit 200 of the ultrasonic diagnostic apparatus. .
  • the RF signal frame data selection circuit 1082 always refers to the RF signal frame data P secured in the frame memory 1081 as past RF signal frame data at an arbitrary time phase without updating. Therefore, the displacement measurement unit 109 indicates whether the currently secured RF signal frame data N and RF signal frame data P A set of RF signal frame data composed of The settings such as whether to adopt the function shown in Fig. 14 and when to acquire the RF signal frame data P when it is adopted can be arbitrarily switched by the user interface of the ultrasonic diagnostic apparatus. You can change it.
  • the interval between past and present RF signal frame data N and P that constitute one set of RF signal frame data is limited to adjacent frames, multiple RF signals acquired during a series of pressurizing or depressurizing operation processes
  • the amount of pressurization or depressurization given in the time interval between the RF signal frame data that constitutes the frame data set is sufficient for the amount of pressurization or decompression (generally about 1%) suitable for the rendering of elasticity image data. May not be reached.
  • the RF signal frame data selection unit as shown in FIGS. 13 and 14
  • the frame interval between the past and present RF signal frame data can be made sufficiently large, and the elasticity can be improved.
  • An elastic image based on frame data can be appropriately drawn. This is especially true in ultrasonography, where physical limitations on the physique of the subject do not allow the rate of pressurization or decompression during the series of pressurization or decompression operations to be sufficiently high. Useful.

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Pathology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Veterinary Medicine (AREA)
  • Acoustics & Sound (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

 診断の指標であるスコアを与えて表示させることにより、診断の一助として認識させるよう、探触子101から超音波を被検体100に送信させ、該超音波の送信に対応する反射エコー信号を受信し、反射エコー信号に基づいて被検体100の生体組織の変位を計測する変位計測部109と、変位から歪み量或いは弾性率を求め弾性画像を構成する弾性画像構成部と、弾性画像を表示させる表示部107とを備えた超音波診断装置において、弾性画像構成部から出力される弾性画像の情報に基づいて、歪み状態或いは弾性状態を特定するスコアリング手段を備えた。

Description

明 細 書
超音波診断装置
技術分野
[0001] 本発明は、超音波を利用して被検体内の診断部位について、生体組織の硬さまた は柔らかさを表す弾性画像として表示することができる超音波診断装置に関する。 背景技術
[0002] 従来の一般的な超音波診断装置は、超音波送受信を制御する超音波送受信制御 手段と、被検体に超音波を送信及び受信する超音波送受信手段と、この超音波送 受信手段からの反射エコー信号を用いて運動組織を含む被検体内の断層像データ を所定周期で繰り返して得る断層走査手段と、この断層走査手段によって得た時系 列断層像データを表示する画像表示手段とを有して構成されていた。そして、被検 体内部の生体組織の構造を例えば Bモード断層像として表示していた。
[0003] これに対して、最近、被検体の体表面から加圧装置もしくは探触子で人為的に外 力を与え、生体内部組織を圧迫し、時系列的に隣接する 2フレーム (連続 2フレーム) の超音波受信信号の相関演算を利用して、各点における変位を求め、さらにその変 位を空間微分することによって歪みを計測し、この歪みデータを画像化する手法、更 には、外力による応力分布と歪みデータから、生体組織のヤング率等に代表される 弾性率データを画像化する手法が提案されて ヽる。このような歪み及び弾性率デー タ(以下、弾性フレームデータ)を基にした弾性画像によれば、生体組織の硬さゃ柔 らかさを計測して弾性画像として表示することができる。
超音波診断装置として、特許文献 1又は特許文献 2に記載されたものなどがある。 特許文献 1:特開平 5-317313号公報
特許文献 2:特開 2000-060853号公報
[0004] しかし、このような従来の超音波診断装置による生体組織の弾性率データの画像 ィ匕においては、実際の疾患に対し、その硬さの程度を認識させる程度にすぎず、疾 患に対応した診断の指標を与えるものではない。
[0005] この発明は、上述の点に鑑みなされたものであり、疾患部位に対して診断の一助と なる指標を与えることのできる超音波診断装置を提供することを目的とする。
発明の開示
[0006] 上記目的を達成するために、本願発明の超音波診断装置は、探触子から超音波を 被検体に送信させ、該超音波の送信に対応する反射エコー信号を受信し、前記反 射エコー信号に基づいて前記被検体の生体組織の変位を計測する変位計測部と、 前記変位カゝら歪み量或いは弾性率を求め弾性画像を構成する弾性画像構成部と、 前記弾性画像を表示させる表示部とを備えた超音波診断装置にお!ヽて、前記弾性 画像構成部から出力される前記弾性画像の情報に基づ ヽて、歪み状態或いは弾性 状態を特定するスコアリング手段を備えた。この前記スコアリング手段は、前記弾性 画像の面積或いは形状に基づいて、歪み状態或いは弾性状態を特定する。また、 前記反射エコー信号カゝら断層像を構成する断層像構成部を備え、前記スコアリング 手段は前記断層像上の特定領域に対する前記弾性画像の硬い部位の領域の関係 に基づいて、歪み状態或いは弾性状態を特定する。
[0007] よって、弾性画像、断層像及び弾性画像を用いて診断の指標であるスコアを与え て表示させることにより、診断の一助として認識させることができる。
図面の簡単な説明
[0008] [図 1]本発明による超音波診断装置の実施例を示すブロック図である。
[図 2]超音波探触子に圧力計測部 (圧力センサ)を取りつけ、探触子のヘッドと被検 体の間の圧力を計測する方法の一例を示す図である。
[図 3]図 1のスコア入力手段の実行するスコアリング処理の一例を示す図である。
[図 4]スコア入力手段を用いてスコアリング処理を行う場合の画像表示器の表示例を 示す図である。
[図 5]ストレス信号によるスコアリングを示す図である。
[図 6]圧迫による弾性画像の領域の変化によるスコアリングを示す図である。
[図 7]弾性画像の色相パターンによるスコアリングを示す図である。
[図 8]図 1のスコア入力手段の実行するスコアリング処理の別の実施例を示す図であ る。
[図 9]図 1のスコア入力手段の実行するスコアリング処理のさらに別の実施例を示す図 である。
[図 10]領域抽出アプリケーションを用いた自動スコアリング処理の一例を示すフロー チャート図である。
[図 11]この自動スコアリング処理の動作の一例を模式的に示す図である。
[図 12]スコアリングするための弾性画像を取得するタイミングを示す図である。
[図 13]図 1の RF信号フレームデータ選択部の一実施例を示す図である。
[図 14]図 1の RF信号フレームデータ選択部の別の実施例を示す図である。
発明を実施するための最良の形態
[0009] 以下、本発明の実施例を添付図面に基づいて詳細に説明する。図 1は、本発明に よる超音波診断装置の実施例を示すブロック図である。この超音波診断装置は、超 音波を利用して被検体 100の診断部位について断層像を得ると共に生体組織の硬さ 又は柔ら力さを表わす弾性画像を表示するものである。
[0010] この超音波診断装置は、図 1に示すように、超音波探触子 101と、送波回路 102と、 受信回路 103と、整相加算回路 104と、信号処理部 105とを用いて断層像を構成する 断層像構成部 118として備え、変位計測部 109と、圧力計測部 110と、歪み及び弾性 率演算部 111と、弾性データ処理部 112と、カラースキャンコンバータ 113とを用いて弹 性画像を構成する弾性画像構成部 119として備え、さらに断層像構成部の出力信号 を画像変換する白黒スキャンコンバータ 106と、画像表示器 107と、 RF信号フレームデ ータ選択部 108と、切替加算器 114と、スコアを入力するスコア入力手段 115と、弾性 画像若しくは弾性画像及び断層像から歪み状態或!ヽは弾性状態を特定するスコアリ ング部 116とを備えて構成されている。
[0011] 超音波探触子 101、送波回路 102、受信回路 103、整相加算回路 104及び信号処理 部 105によって、超音波送受信手段が構成される。この超音波送受信手段は、超音 波探触子 101を用いて超音波ビームを被検体 100の体内で一定方向に走査させるこ とにより、一枚の断層像を得るものである。超音波探触子 101は、多数の振動子を短 冊状に配列して形成されたものであり、機械式または電子的にビーム走査を行って 被検体 100に超音波を送信及び受信するもので、図示は省略したがその中には超音 波の発生源であると共に反射エコーを受信する振動子が内蔵されている。各振動子 は、一般に、入力されるパルス波、または連続波の送波信号を超音波に変換して発 射する機能と、被検体 100の内部から反射する超音波を受けて電気信号の受波信号 に変換して出力する機能を有して形成される。
[0012] 送波回路 102は、超音波探触子 101を駆動して超音波を発生させるための送波パ ルスを生成すると共に、内蔵された送波整相加算回路によって送信される超音波の 収束点をある深さに設定するものである。受信回路 103は、超音波探触子 101で受信 した反射エコー信号を所定のゲインで増幅するものである。増幅された各振動子の 数に対応した数の受波信号がそれぞれ独立した受波信号として整相加算回路 104に 入力される。整相加算回路 104は、受信回路 103で増幅された受波信号を入力し、そ れらの位相を制御し、一点又は複数の収束点に対して超音波ビームを形成するもの である。信号処理部 105は、整相加算回路 104からの受波信号を入力してゲイン補正 、ログ圧縮、検波、輪郭強調、フィルタ処理等の各種信号処理を行うものである。
[0013] 白黒スキャンコンバータ 106は、前述の超音波送受信手段の信号処理部 105から出 力される反射ェコ一信号を用 Vヽて運動組織を含む被検体 100内の RF信号フレーム データを超音波周期で取得し、この RF信号フレームデータを切替加算器 114を介し て画像表示器 107に表示するものである。従って、白黒スキャンコンバータ 106は、テ レビジョン方式の周期で RF信号フレームデータを順次読み出すための断層走査手 段及びシステムの制御を行うための手段、例えば、信号処理部 105からの反射エコー 信号をディジタル信号に変換する A/D変^^と、この A/D変 でディジタルィ匕され た断層像データを時系列に記憶する複数枚のフレームメモリと、これらの動作を制御 するコントローラなどを含んで構成される。
[0014] 画像表示器 107は、白黒スキャンコンバータ 106によって得られた時系列の断層像 データすなわち Bモード断層像を表示するものであり、切替加算器 114を介して白黒 スキャンコンバータ 106から出力される画像データをアナログ信号に変換する D/A変 と、この D/A変 力ゝらのアナログビデオ信号を入力して画像として表示する力 ラーテレビモニタとから構成される。
[0015] この実施形態においては、整相加算回路 104の出力側力 分岐して RF信号フレー ムデータ選択部 108と変位計測部 109とが設けられると共に、組織弾性率を求める場 合は、これと並列に圧力計測部 110が設けられている。変位計測部 109及び圧力計測 部 110の後段には、歪み及び弾性率演算部 111と弾性データ処理部 112とカラースキ ヤンコンバータ 113とが設けられ、このカラースキャンコンバータ 113と白黒スキャンコン バータ 106の出力側には切替加算器 114が設けられている。そして、白黒スキャンコン バータ 106とカラースキャンコンバータ 113の出力は、スコアリング部 116に導入され、 そこでスコアリング処理が実行される。このスコアリング処理の詳細については後述す る。
[0016] 変位計測部 109は、 RF信号フレームデータ選択部 108によって選択された 1組の RF 信号フレームデータに基づいて 1次元もしくは 2次元相関処理を実行し、断層像上の 各点の変位もしくは移動ベクトル (変位の方向と大きさ)を計測するものである。この移 動ベクトルの検出法としては、例えば、特許文献 1に記載されたようなブロックマツチン グ法とグラジェント法とがある。ブロックマッチング法は、画像を例えば N X N画素から なるブロックに分け、現フレーム中の着目しているブロックに最も近似しているブロック を前フレーム力 探し、これを参照して予測符号ィ匕を行うものである。
[0017] 圧力計測部 110は、被検体 100の診断部位の体腔内圧力を計測又は推定するもの である。この超音波診断装置は、探触子ヘッド 1011に設けられた超音波探触子 101を 用いて制御部 200の制御の下で超音波送受信を行 、つつ、探触子ヘッド 1011に設 けられた加圧器(図示せず)によって加圧もしくは減圧し、被検体 100の診断部位の 体腔内に応力分布を与える方法を採用している。この方法では、探触子ヘッド 1011と 被検体 100との間にどの程度の圧力が印加されているかを計測するために、例えば、 図 2に示すように、棒状部材に架力る圧力を検出する圧力センサ 1012を探触子ヘッド 1011の側面に取り付け、探触子ヘッド 1011と被検体 100の間の圧力を任意の時相で 測定し、測定された圧力値を歪み及び弾性率演算部 111に送出するようになって 、る 。なお、図 2では、探触子ヘッド 1011を加圧 '減圧する加圧器については省略してあ る。
[0018] 歪み及び弾性率演算部 111は、変位計測部 109及び圧力計測部 110からそれぞれ 出力される移動量 (変位)及び圧力に基づ!、て断層像上の各点の歪み或いは弾性 率を演算して、歪みもしくは弾性率の数値データ(弾性フレームデータ)を生成し、そ れを弹性データ処理部 112に出力するものである。歪み及び弾性率演算部 111が行う 歪みの演算についは、例えば、圧力のデータを必要とせず、その変位を空間微分す ることによって計算上求める。また、弾性率の内の一つであるヤング率の演算につい ては、圧力の変化を移動量の変化で除することによって計算上求める。
[0019] カラースキャンコンバータ 113は、弾性データ処理部 112から出力される弾性フレー ムデータと、超音波診断装置の制御部 200から出力される命令又は弾性データ処理 部 112から出力される弾性フレームデータの中の階調化選択範囲とする上限値及び 下限値とを入力し、その弾性フレームデータカゝら弾性画像データとして赤、緑、青な どの色相情報を付与する色相情報変換手段を備えて 、る。この色相情報変換手段 は、例えば、弾性データ処理部 112から出力される弾性フレームデータにおいて、歪 みが大きく計測された領域を、弾性画像データ内で赤色コードに変換し、逆に歪み が小さく計測された領域を、弾性画像データ内で青色コードに変換するように動作す る。また、カラースキャンコンバータ 113は、上述の白黒スキャンコンバータ 106で構成 してもよい。この場合には、歪みが大きく計測された領域は、弾性画像データ内で輝 度を明るくし、逆に歪みが小さく計測された領域は、弾性画像データ内で輝度を暗く すればよい。
[0020] 切替加算器 114は、白黒スキャンコンバータ 106からの白黒の断層像データとカラー スキャンコンバータ 113からのカラーの弾性画像データとを入力し、両画像を加算又 は切替る手段となるもので、白黒の断層像データだけ又はカラーの弾性画像データ だけを出力したり、あるいは両画像データを加算合成して出力したりするように切替る ようになつている。また、例えば、特許文献 2に記載してあるように、 2画面表示におい て白黒断層像とカラーもしくは白黒スキャンコンバータによる白黒弾性画像を同時に 表示しても良い。切替加算器 114から出力された画像データは画像表示器 107に出 力されるようになっている。
[0021] 次に、本実施形態に係るスコアリングについて説明する。
[0022] 図 3は図 1のスコアリング部 116の実行するスコアリング処理の第 1の実施形態を示す 図であり、図 3 (A)は、白黒スキャンコンバータ 106から出力される白黒の断層像デー タ(Bモード断層像)の一例を示す図である。図 3 (B)一 (F)は、カラースキャンコンパ ータ 113から出力されるカラーの弾性画像データの一例を示す図であり、それぞれの 症状が異なる場合を模式的に表したものである。これらの弾性画像データは、実際は 、その硬さに応じたカラー表示を行うものである。例えば、柔らかい部分 (歪みの大き な部分)は赤色、硬い部分 (歪みの小さな部分)は青色、両者の中間は緑色として、 連続的な色の変化として表示される。なお、図 3では、このような色で表示することが できないので、網掛けにてその様子を示すこととし、軟らかい領域を荒い網掛けで、 硬 、領域を細力 、網掛けで表示することにする。図 3 (A)の Bモード断層像の場合、 中央の楕円形の部分 61は、乳腺腫瘍等の疾患部位であり、低エコー領域と呼ばれる ものである。また、弾性画像において低エコー領域 61の同一座標領域に対応する領 域を図 3 (B)—(F)において、点線の楕円形で表示することにする。この低エコー領 域 61に対して、弾性画像データがどのように表示されて!、るのかその状態に応じて、 弾性画像データのスコア入力手段 115を操作して、スコアリング処理を行い、それぞ れの弹性画像をスコア 1一 5に分類する。各スコア 1一 5は次のような基準になっており 、この基準を満たす力どうかの判断を観察者が行いスコアリング処理を行う。以下、各 スコア 1一 5の基準について説明する。
スコア 1 :図 3 (B)に示すように、弾性画像による硬い領域 62力 低エコー領域 61の 内部に、明確には認識されないような症例
スコア 2 :図 3 (C)に示すように、弾性画像による硬い領域 63, 64が、低エコー領域 61 の中央部を避けて、 1部分的に認識されるような症例
スコア 3 :図 3 (D)に示すように、弾性画像による硬い領域 65力 低エコー領域 61の 境界部(輪郭)にまでは到達しないが、低エコー領域 61の中央部を含んで認識される ような症例
スコア 4:図 3 (E)に示すように、弾性画像による硬い領域 66力 低エコー領域 61の 境界部 (輪郭)まで一様に分布して認識されるような症例
スコア 5 :図 3 (F)に示すように、弾性画像による硬い領域 67が、低エコー領域 61の 全面を覆うと同時に、低エコー領域 61の境界部(輪郭)を超えて外部の領域にまで拡 力 て認識されるような症例
ここで、前記硬い領域は、弾性画像としての歪み画像においては、ある歪みの閾値 Tsより小さい歪みを計測した領域として識別すればよぐ例えば、歪み 0パーセントを 閾値 Tsとして設定すれば、前記硬い領域としては、体表からの圧迫によっても全く圧 縮が成されな力つた領域 (計測点群)として選別される。なお、診断部位の走査時に その押圧力を取り込み、弾性画像として弾性率画像を用いる場合においては、ある 弾性率の閾値 Tyより大き ヽ弾性率を計測した領域を硬!ヽ領域として識別するように すればよい。
[0024] このようにスコアリング (分類)することで、疾患の判断の指標とすることができ、診断 支援を行うことができる。図 4は、スコア入力手段 115を用いてスコアリング処理を行う 場合の画像表示器 107の表示例を示す図である。画像表示器 107の表示画面 1071に は、図 3 (B)一 (F)のような Bモード断層像上に組織の硬さの分布を示した弾性画像 データの半透明カラー像を重ねたものを表示し、その右横にスコアリング指標を入力 する欄 1072と、各スコアの指標を表示しておく表示領域 1073を設けてある。この表示 領域 1073には、前述のような各スコアの内容を簡単に明記したものを表示したり、図 3 (B)—(F)のような図を併せて表示したりする。検査者は、スコアリング指標 1073と対 比しながら、スコアリング指標人力欄 1072に、表示画面 1071を見て適当と思われるス コアを入力する。これらの動作は、スコア入力手段 115によって実行される。なお、表 示画面 1071に表示される画像は、図 3 (A)のように Bモード断層像を別の窓に並列に 表示してもよい。これによつてスコアリング処理時の視認性を向上することが可能とな る。スコアリング指標入力欄 1072に入力されたスコアは、超音波診断装置に備えられ た患者個々のレポート機能にもリンクし、レポート機能に備えられたスコアリング指標 入力欄にも同時に入力されるようになっていても良い。また、図 3 (B)—(F)のようなス コアリングの指標となるモデル図を表示領域 1073に併記するようにしてもよい。
[0025] 上述の実施形態では、検査者が表示画面 1071に表示された画像を見てスコアリン グ処理を行う場合について説明した力 スコアリング部 116によって画像処理を使って 自動的にスコアリング (分類)する。以下、このスコアリング部 116による自動スコアリン グ処理について説明する。スコアリング部 116は、カラースキャンコンバータ 113から出 力されるカラーの弾性画像データ (歪み像)において公知のエッジ検出ソフトウェアで 検出するエッジ検出を行い、硬い部分の領域 62— 67の面積 Bを画素数等で求める。 [0026] 第 2の実施形態では、弾性画像データのみでスコアリングを行うことを示す。例えば 、基準情報としてリファレンス画像となる弾性画像を予め取得しておき、弾性画像の 領域のピクセル数を面積 Bとして、弾性画像の面積によるスコアを下記のとおり設定し ておく。
[0027] XK B ならば スコア 1
X2< B≤X1 ならば スコア 2
X3< B≤X2 ならば スコア 3
X4< B≤X3 ならば スコア 4
B<X4 ならば スコア 5
そして、新たに取得した弾性画像データの面積を上記の通り設定した面積 Bに当て はめてスコアリングを行う。また、この弾性画像が一連で繋がっている画像であれば、 最低スコア 2以上と設定したり、滑らかな円系で繋がっている画像であれば、最低スコ ァ 3以上と設定したりし、弾性画像の形状でスコアを定めてもよい。
[0028] また、被検体を圧迫することで、弾性画像上にストレス信号が出るという特性を利用 したスコアリング処理を図 5に示す。図 5 (a) (b)は、被検体 100を探触子 101で圧迫し た状態を示した図である。悪性の腫瘍を圧迫すると、図 5 (b)に示されるようにストレス 信号 122が表示される。なお図 5 (a)は良性の腫瘍を示す。
[0029] スコアリング部 116において、予めストレス信号 122の幅 123に対応するスコアを設定 する。例えば、この幅が大きいほど弾性画像は悪性であるというようにスコアを図 4に 示されるように設定する。そして、スコアに対応して設定された幅と、得られた弾性画 像 121のストレス信号 122の幅 123と比較を行うことにより、弾性画像 121スコアリングを 行う。ここではストレス信号 122の幅を用いてスコアリングを行う例を示した力 ストレス 信号 122の面積を用いてスコアリングを行ってもよい。このように弾性画像のストレス信 号によってスコアリングを行うことができる。
[0030] また、被検体を圧迫することで変化する弾性画像を用いてスコアリング処理を行う例 を図 6に示す。悪性の腫瘍の場合、被検体を徐々に圧迫すると例えば青色表示され る領域が徐々に小さくなるという特性がある。逆に良性の腫瘍の場合、被検体を徐々 に圧迫すると青色表示される領域が徐々に大きくなるという特性がある。そこで、この 特性を用 、てスコアリング部 116にお 、て圧迫により生じる領域の色変化に対応する 特性に対するスコアを図 4に示されるように設定する。そして、弾性画像が圧迫に従つ て変化する色にしたがってスコアリングを行う。このように圧迫による前記弾性画像の 前記領域の変化に基づいてスコアリングを行うことができる。
[0031] シストパターンを利用してスコアリング処理を行う例を図 7に示す。シストパターンと は、表示される弾性画像のある特定領域の深度が浅い部分 141が青色、深度が深い 領域 142が赤色で対に表示されるパターンである。なお、このシストパターンで表示さ れる領域は経験的に良性の腫瘍であることが分力つている。そこで、スコアリング部 116においてシストパターンは良性のスコアと別途設定しておき、弾性画像上でシスト パターンを認識したら良性であると 、うスコアリングを行う。このように弾性画像の色相 パターンに基づいてスコアリングを行うこともできる。
[0032] 上述の実施形態では、弾性画像データのみでスコアリング処理を行う場合につい て説明したが、第 3の実施形態では、基準情報として断層像データを用いた弾性画 像データによる自動スコアリング処理について、図 3等を用いて説明する。まず、スコ ァリング部 116は、白黒スキャンコンバータ 106から出力される白黒の断層像データ(B モード断層像)の疾患部位の低エコー領域 61を公知のエッジ検出ソフトウェアで検出 する。検出された領域 61の面積 Aを画素数等で求める。次に、スコアリング部 116は、 カラースキャンコンバータ 113から出力されるカラーの弾性画像データ(歪み像)にお V、ても同様のエッジ検出を行 、、疾患部位の低エコー領域 61の内部に計測される硬 い部分の領域 62— 67の面積 Bを画素数等で求める。ここで、硬い領域が、領域 63, 64のように、複数に散在する場合は、それらの面積和を面積 Bとする。
[0033] スコアリング部 116は、求められたそれぞれの面積 A, Bの関係に基づいて、次のよう なスコアリング処理を行う。
[0034] まず、前記面積 A, Bの比として、 B/A力 ¾として演算される。また、あら力じめ超音波 診断装置において、面積比の閾値として
スコア 1の閾値 Thl (例えば 0. 1)
スコア 2の閾値 Th2 (例えば 0. 3)
スコア 3の閾値 Th3 (例えば 0. 7) スコア 4の閾値 Th4 (例えば 1. 0)
が設定されている。そこで、
Z<Thl ならば スコア 1
Thl <Z≤Th2 ならば スコア 2
Th2<Z≤Th3 ならば スコア 3
Th3<Z≤Th4 ならば スコア 4
Th4<Z ならば スコア 5
の判定を行うようになって!/、る。
[0035] また、上述の実施例にカ卩えて、例えば、図 8に示されるように、低エコー領域 61内の 中央部に、低エコー領域の面積より小さな、一定の比率を有した (例えば面積 A X 0.6 )関心領域 61aを設け、その関心領域 61aに対して、上述と同様に、弾性画像におけ る硬 、領域がどの程度占有するかと 、つた評価を上述のスコアリング方法にカ卩えて 行ってもよぐこれにより、低エコー領域に対する硬い領域の占有面積比率だけでな ぐ低エコー領域内部及び外部において、硬い領域がどのような空間的な分布として 配置されているかをスコアリングの評価に加えることになり、評価の精度を高める結果 を期待することができる。
[0036] 更にまた、例えば、図 9に示されるように、前記中央関心領域のみならず、低エコー 領域の境界を基準として、複数の関心領域 62bを境界の内部、もしくは、内外に設け 、これらの複数の関心領域に対して、上述と同様に、弾性画像における硬い領域が どの程度占有するかといった評価を行ってもよぐまた、複数の関心領域 62bのそれ ぞれについて、それぞれが占める面積の内、弾性画像における硬い領域が占める割 合がある一定の比率を超えて 、るかどうかをまず評価し、その比率を超えて 、た関心 領域が、いくつ計上されたかをスコアリングの評価に利用してもよい。また、エッジ検 出は、公知のスネークやリージヨングローイング法などを用いることで達成できる。シャ ドーなどでエッジが隠れている場合は、 ROIを設定したり、あるいは、周囲のエッジ情 報力も補間するなどの方法で対応することができる。
[0037] また、上述の実施形態では、低エコー領域内に分布する硬い領域の占有率を用い て、スコアリングする例を示した力 別の実施例として、低エコー領域内に含まれる弹 性の計測点群を母集団とした統計処理を施し、その統計的特徴量に基づ ヽてスコア リングしてもよぐ例えば、その平均値を用いる方法を以下のように説明する。
低エコー領域に内在する計測点群の総要素数を Nとし、各計測点の歪みもしくは弾 性率を、 Ei (i=l, 2, 3 · ·Ν)とする。そこで、低エコー領域に内在する計測点の歪みも しくは弾性率の平均値 Emを以下のように演算する。
(平均値Em) = Σ Ei (i= l, 2, 3 · ·Ν)
[0038] また、あら力じめ超音波診断装置において、平均値の閾値として
スコア 1の閾値 Tml
スコア 2の閾値 Tm2
スコア 3の閾値 Tm3
スコア 4の閾値 Tm4
が設定されており、弾性を反映した値として弾性率を演算している場合は、
Tml <Tm2<Tm3<Tm4
の大小関係を保持して設定されており、
Zく Tml ならば スコア 1
Thl <Z≤Tm2 ならば スコア 2
Th2<Z≤Tm3 ならば スコア 3
Th3<Z≤Tm4 ならば スコア 4
Th4<Z ならば スコア 5
の判定を行うようになって!/、る。
[0039] また、弾性を反映した値として歪みを用いている場合は、
Tml >Tm2>Tm3 >Tm4
の大小関係を保持して設定されており、
TmKZ ならば スコア 1
Tm2<Z≤Tml ならば スコア 2
Tm3<Z≤Tm2 ならば スコア 3
Tm4<Z≤Tm3 ならば スコア 4
Z<Tm4 ならば スコア 5 の判定を行うようになって!/、る。
[0040] 以上の説明では、例として平均値を用いたが、本発明はこれに限らず、低エコー領 域内に含まれる弾性の計測点群を母集団とした統計処理を施し、その統計的特徴量 に基づ!/、てスコアリングすることが重要である。
[0041] スコアリング部 116が画像処理を用いて自動的にスコアリング (分類)する場合の具 体例について説明する。まず、スコアリング部 116は、図 3に示されるように、白黒スキ ヤンコンバータ 106から出力される白黒の断層像データ (Bモード断層像)の疾患部位 の低ェコ一領域 61の面積 Aを公知の領域抽出ソフトである、例えばリージョングロー イング法によって検出する。検出された領域の面積 Aを画素数等で求める。次に、ス コアリング部 116は、カラースキャンコンバータ 113から出力されるカラーの弾性画像デ ータ(歪み像)において、面積 Aの領域内における硬い領域 62— 67として示される、 装置使用者が任意に定めることが可能な、所定の輝度値以上をもつピクセル数を数 える。このピクセル数を面積 Bとする。この面積 Aと面積 B (所定の輝度値以上をもつピ クセル数)とに基づき、次のようなスコアリング処理を行う。
[0042] まず、面積 (ピクセル数) A, Bの比として、 B/A力 ¾として演算される。また、予め超音 波診断装置において、面積比の閾値として
スコア 1の閾値 Thl (例えば 0. 1)
スコア 2の閾値 Th2 (例えば 0. 3)
スコア 3の閾値 Th3 (例えば 0. 7)
スコア 4の閾値 Th4 (例えば 1. 0)
が設定されている。
[0043] そこで、今、 Z< 0. 7の時は、次に示す規則に則りスコアを返す。
[0044] Z<Thl ならば スコア 1
Thl<Z≤Th2 ならば スコア 2
Th2≤Z≤Th3 ならば スコア 3
さらに、 Th3≤Z (0.7≤Z)ならば、得られた弾性画像データ(歪弾性画像)において 、前述と同様、公知の領域抽出アプリケーション (例えば、リージヨングローイング法) にて、硬い領域 65— 67として示される、所定の輝度以上をもつ領域を抽出する。この 際、抽出領域は前述の低エコー領域 61内に留まることなぐその辺縁を含む領域をも 抽出対象となる。
[0045] ここで得られた抽出面積を Cとする。面積 Aと面積 Cの比として C/Aが Z'として演算さ れる。この Z'の値が、
Th3<Z' ≤Th4 ならば スコア 4
Th4く Z' ならば スコア 5
を返す。
[0046] 上述の実施形態では、スコア 3までは、弾性画像データにお!、て、領域抽出処理を 要せず、スコアリングが可能となる。
[0047] 図 10は、領域抽出アプリケーションを用いた自動スコアリング処理の一例を示すフ ローチャート図である。図 uは、この自動スコアリング処理の動作の一例を模式的に 示す図である。ここでは、領域抽出アルゴリズムとしてリージヨングローイング法( Region Growing)を例に説明する。
[0048] ステップ S91では、歪弾性像を得るために、任意の手段にて加圧された後、描出さ れる Bモード断層像にて、領域抽出したい低エコー部を選択する。この選択処理は、 図 11 (A)に示すように低エコー領域 61内の任意の位置をソースポイント(図中 SPと記 載)として指定することによって行う。
[0049] ステップ S92では、ステップ S91で指示されたソースポイント SPと、任意の手段にて与 えられた閾値に基づいて、図 11 (B)に示すような低エコー領域 61の抽出を行う。リー ジョングロ一イング法は、指示されたソースポイント SPに隣接する領域で、輝度値差が 予め設定された値(閾値)内に収まる領域を抽出する手法である。従って、抽出され た低エコー領域 61の面積を、ピクセル数を数えるなどの手段により計数する。図 11 (B )では、この面積を SBWとして表記する。
[0050] ステップ S93では、この低エコー領域内において、描出された歪弾性像で、使用者 が任意に定める硬 、領域を示す輝度値よりも、より硬 、領域を表す輝度値をもつ面 積をピク
セル数で数える。図 11 (C)では、この硬い領域を Sblと表記する。
[0051] ステップ S94では、硬い領域の面積 Sblが低エコー領域の面積 SBWX 0.7よりも小さ いか否かの判定を行う。判定の結果、硬い領域の面積 Sblが低エコー領域の面積 SBWX 0.7よりも小さ!/、(SBWX 0.7>Sblである)場合にはステップ S95に進む。そして 、ステップ S95ではそれぞれの値に応じたスコア 1一 3までの値が返されるようになる。
[0052] 逆に、ステップ S94の判定の結果、硬い領域の面積 Sblが低エコー領域の面積 SBW
X 0.7以上(SBWX 0.7≤Sblである)の場合には、低エコー領域内に硬 、部位の面積 Sblが収まっていない可能性があるため、ステップ S96に進む。
[0053] ステップ S96では、弾性画像データにおいて再度、リージヨングローイング法による 領域 (面積 Sbl')の抽出が行われる。この場合は、 Bモード断層像の低エコー領域の 周辺領域を含む領域までの抽出が行われる。ここでの領域抽出処理ではソースボイ ント SPは、前述の Bモード断層像にて指示されたものをそのまま用いても良いし、別途 新たに設定しても良いものとする。ここで抽出された領域の面積を Sbl'とする。そして 、ステップ S97では、低エコー領域の面積 SBWと面積 Sbl'との比(SBWZSbl')に応じ たスコア 4, 5の値が返されるようになる。具体的には、 Sbl 'く SBWであれば、スコア 4 が返され、それ以外ではスコア 5が返されることになる。
[0054] このようにして、公知の領域抽出アプリケーション(例えば、リージヨングローイング法 )にて領域抽出を行うことで、組織の硬さ分布に応じたスコアリング処理を自動で行う ことが可能となる。ここで、上述の閾値は全て一例に過ぎず、使用者の判断によって 、 自由に変更可能であることは言うまでもない。
[0055] また、上述の実施形態では、 Bモード断層像における疾患部位の低エコー領域の 輪郭情報を公知のエッジ検出ソフトウェアで検出する例を示した力 別の実施例とし て、超音波装置のインターフェイス(マウス、トラックボールなど)を用いて、検査者が 表示された Bモード断層像上に疾患部位の輪郭を入力し、その輪郭情報を元にして 領域 61の面積 Aを求め、上述の方法による面積 Bと比較して自動的にスコアリングを 判断するようにしてもよい。なお、上述のスコアリングは、一例であり、ひずみ像を含 む弾性画像を用いてスコアリングすることが重要であり、上述のスコアリングで用いた 数値などは一例に過ぎず、実際の臨床例に従って種々最適な値を用いることが望ま しいことは言うまでもない。また、上述の説明においては、特に乳腺領域を対象とした スコアリング方法を例とした力 それ以外の領域では、その領域の症例にあわせたス コアリング方法を定義するものであることは言うまでもない。また、以上の実施例による スコアリング方法はそれぞれを単独としてスコアリングに利用するのみならず、複数の 方法を組み合わせて利用し、ひとつのスコアリングを行ってもよいのは言うまでもない
[0056] 上述の実施形態では、スコアリング部 116が白黒スキャンコンバータ 106から出力さ れる白黒の断層像データ (Bモード断層像)とカラースキャンコンバータ 113から出力さ れるカラーの弾性画像データ (歪み像)を用いてスコアリング処理を行う場合につい て説明したが、信号処理部 105からの出力、弾性データ処理部 112又は歪み及び弹 性率演算部 111からの出力を用いて同様のスコアリング処理を行うようにしてもょ 、。 また、上述の実施例において、図 1に圧力計測部を記した力 歪を求めて、歪像を表 示する場合には、圧力計測部は、なくてもよい。
[0057] 図 12は、上述した実施形態の歪み及び弾性率演算部 111で歪み又は弾性率を演 算するタイミングを示す図である。探触子 101で被検体 100を圧迫して 、な 、状態を 0 、圧迫の限界を MAXで示している。弾性画像を取得するタイミングは軽い圧迫状態 のタイミングであり、圧迫の限界の 3— 20%程度の圧迫した状態の弾性画像を取得し 、取得した弾性画像に対してスコアリングを行う。
[0058] 具体的には、圧力センサ 1012によって、予め圧迫の限界の圧力値を検出して、そ の圧力値を圧力 MAX値としてスコアリング部 116に記憶しておく。記憶された圧力 MAX値を参考にして、スコア入力手段 115を用いて、圧迫して取得する弾性画像の 圧力値を設定する。なお、設定する圧力値は、例えば圧力 MAX値の 10%の圧力値と 任意に定めておいてもしてもよい。そして、被検体 100を圧迫して、設定した圧力値に 到達した時点の弾性画像を取得し、その取得した弾性画像にっ 、てのスコアリングを 行う。
[0059] 被検体 100を強く圧迫した状態の弾性画像は、硬 、部位と軟らカ 、部位の見極め が難しくなる場合があるため、上述したように軽く圧迫したタイミングで取得した弾性 画像においてスコアリングすることにより、スコアリング精度を向上させることができる。
[0060] 上述の実施形態に係る RF信号フレームデータ選択部 108の動作について、図 13を 用いて説明する。図 13は、図 1の RF信号フレームデータ選択部の一実施例を示す図 である。 RF信号フレームデータ選択部 108は、変位計測の基準となる一つの RF信号 フレームデータとして過去に遡るフレーム数(現フレームデータとのフレーム間隔数) を任意に選択するものである。すなわち、 RF信号フレームデータ選択部 108は、超音 波診断装置のフレームレートで整相加算回路 104から経時的に次々と出力される RF 信号フレームデータを、 RF信号フレームデータ選択部 108に備えられたフレームメモ リ 1081内に順次確保する。 RF信号フレームデータ選択部 108は、フレームメモリ 1081 内に現時点で確保されているものを RF信号フレームデータ Nとする。 RF信号フレーム データ選択部 108は、超音波診断装置の制御部 200からの制御命令に従って時間的 に過去の RF信号フレームデータ N- 1, N-2, N-3, · · · , N- Mの中から 1つの RF信号フ レームデータを選択し、それを RF信号フレームデータ Xとして、 RF信号フレームデー タ選択回路 1082に一時的に格納する。 RF信号フレームデータ選択部 108は、フレー ムメモリ 1081に格納されて 、る最新の RF信号フレームデータ Nと RF信号フレームデ ータ選択回路 1082に格納されている RF信号フレームデータ Xとを、並列的に変位計 測部 109に出力する。
[0061] すなわち、 RF信号フレームデータ選択部 108は、まず、変位計測部 109に送出する 1組の RF信号フレームデータを構成する過去の RF信号フレームデータ Xとして、現在 の RF信号フレームデータ Nに時間的に隣接する RF信号フレームデータ N-1のみなら ず、過去の RF信号フレームデータ Xとして Mフレーム(M=l, 2, 3, · · · )を間引いた RF 信号フレームデータ N-Mを任意に選択することができるものである。なお、間引いた フレーム間隔数 M (M=1, 2, 3, · · · )は、超音波診断装置のユーザーインターフェイス によって任意に設定 ·変更できるようになって!/、る。
[0062] 図 14は、図 1の RF信号フレームデータ選択部の別の実施例を示す図である。図 14 の RF信号フレームデータ選択部 108は、超音波診断装置の制御部 200からの制御命 令に従って、過去のある時相 Pにおいて取得された RF信号フレームデータ Pをフレー ムメモリ 1081内に確保する。 RF信号フレームデータ選択回路 1082は、フレームメモリ 1081内に確保されている RF信号フレームデータ Pを、アップデートすることなぐ任意 の時相における過去の RF信号フレームデータとして常に参照する。従って、変位計 測部 109には、現在確保された RF信号フレームデータ Nと RF信号フレームデータ Pか ら構成される 1組の RF信号フレームデータが取り込まれる。図 14のような機能を採用 するかどうか、また採用した場合に RF信号フレームデータ Pを取得するタイミングはど うするかなどの設定は、超音波診断装置のユーザーインターフェイスによって任意に 切替 '設定'変更できるようになって 、る。
1組の RF信号フレームデータを構成する過去と現在の RF信号フレームデータ N, P の間隔を隣接フレームに限定した場合、一連の加圧もしくは減圧操作過程の間に取 得された複数の RF信号フレームデータの組を構成する RF信号フレームデータの間 の時間間隔において与えられた加圧量もしくは減圧量は、弾性画像データの描出に 適した加圧量もしくは減圧量 (一般に 1%程度)に十分に達することができない場合が ある。これに対して、図 13及び図 14に示すように RF信号フレームデータ選択部を構 成することによって、過去と現在の RF信号フレームデータの間のフレーム間隔を十分 に大きくすることができ、弾性フレームデータによる弾性画像を適切に描出することが 可能となる。これは特に、超音波検査において、被検者の体格の物理的な制約によ り、一連の加圧もしくは減圧操作過程における加圧もしくは減圧の速度を十分に大き くできないような状況において非常に有用である。

Claims

請求の範囲
[1] 探触子から超音波を被検体に送信させ、該超音波の送信に対応する反射エコー 信号を受信し、前記反射エコー信号に基づ 、て前記被検体の生体組織の変位を計 測する変位計測部と、前記変位から歪み量或いは弾性率を求め弾性画像を構成す る弾性画像構成部と、前記弾性画像を表示させる表示部とを備えた超音波診断装置 において、
前記弾性画像構成部から出力される前記弾性画像の情報に基づ!/ヽて、歪み状態 或いは弾性状態を特定するスコアリング手段を備えることを特徴とする超音波診断装 置。
[2] 前記スコアリング手段は、前記弾性画像の面積或!、は形状に基づ!、て、歪み状態 或いは弾性状態を特定することを特徴とする請求項 1記載の超音波診断装置。
[3] 前記面積は、前記弾性画像の画素数で決定されることを特徴とする請求項 2記載 の超音波診断装置。
[4] 前記弾性画像が複数に存在する場合、前記面積は該複数の弾性画像の面積和で あることを特徴とする請求項 2記載の超音波診断装置。
[5] 前記反射エコー信号から断層像を構成する断層像構成部を備え、前記スコアリン グ手段は前記断層像上の特定領域に対する前記弾性画像の硬い領域の関係に基 づ!ヽて、歪み状態或いは弾性状態を特定することを特徴とする請求項 1記載の超音 波診断装置。
[6] 前記断層像或!ヽは前記弾性画像に前記領域を設定する領域設定手段を備えるこ とを特徴とする請求項 5記載の超音波診断装置。
[7] 前記領域設定手段は前記断層像上の低エコー領域を第 1の領域とし、前記弾性画 像上の硬 、部位の領域を第 2の領域として設定して設定し、前記スコアリング手段は 前記第 1の領域と第 2の領域の比に基づいて歪み状態或いは弾性状態を特定するこ とを特徴とする請求項 6記載の超音波診断装置。
[8] 前記領域設定手段は、前記断層像、或いは前記弾性画像の輝度によって領域を 抽出するリージヨングローイング手段又はエッジ検出手段によって領域を求めることを 特徴とする請求項 6記載の超音波診断装置。
[9] 前記リージヨングローイング手段は、前記領域内の位置をソースポイントとして設定 し、該設定されたソースポイントに隣接する領域で輝度差が予め設定された閾値内 に収まる領域を抽出することを特徴とする請求項 8記載の超音波診断装置。
[10] 前記領域設定手段は、前記断層像、或いは前記弾性画像上に輪郭を入力して、 前記輪郭情報を元にして前記領域を求めることを特徴とする請求項 6記載の超音波 診断装置。
[11] 前記スコアリング手段は、前記低エコー領域に複数の関心領域を設定し、前記複 数の関心領域における硬 、領域の占有度合いに基づ 、て、歪み状態或いは弹性状 態を特定することを特徴とする請求項 7記載の超音波診断装置。
[12] 前記スコアリング手段は、前記低エコー領域の中央部に関心領域を設け、前記弹 性画像における硬 、領域の占有度合いに基づ!/、て、歪み状態或いは弾性状態を特 定することを特徴とする請求項 7記載の超音波診断装置。
[13] 前記被検体を圧迫する圧力値を設定する圧力設定手段を備え、該設定した圧力 値に圧迫した際の弾性画像に基づ!/ヽて歪み状態或!ヽは弾性状態を特定することを 特徴とする請求項 1記載の超音波診断装置。
[14] 前記スコアリング手段は、前記領域に基づ!、て、歪み状態或!、は弾性状態を複数 段階に分類することを特徴とする請求項 1記載の超音波診断装置。
[15] 前記表示部にスコアリング入力欄を設け、前記入力手段を用いて前記スコアリング 入力欄にスコアを入力することを特徴とする請求項 1記載の超音波診断装置。
[16] 前記変位計測部に送出する 1組の前記反射エコー信号のフレームデータを選択す るフレームデータ選択部を備え、該選択された 1組のフレームデータに基づ 、て前記 弾性画像を構成することを特徴とする請求項 1記載の超音波診断装置。
[17] 前記フレームデータ選択部は、前記 1組のフレーム間隔を任意に選択し、該選択さ れた 1組のフレームデータに基づいて前記弾性画像を構成することを特徴とする請求 項 16記載の超音波診断装置。
[18] 前記フレームデータ選択部は、 1組の隣接するフレームのフレームデータを選択す ることを特徴とする請求項 16記載の超音波診断装置。
[19] 前記スコアリング手段は、前記弾性画像のストレス信号に基づ!ヽて、歪み状態或 ヽ は弾性状態を特定することを特徴とする請求項 1記載の超音波診断装置。
[20] 前記スコアリング手段は、徐々に圧迫することによる前記弾性画像の領域の変化に 基づ 、て、歪み状態或いは弾性状態を特定することを特徴とする請求項 1記載の超 音波診断装置。
[21] 前記スコアリング手段は、前記弾性画像の色相パターンに基づいて、歪み状態或 いは弾性状態を特定することを特徴とする請求項 1記載の超音波診断装置。
[22] 探触子から超音波を被検体に送信させ、該超音波の送信に対応する反射エコー 信号を受信して断層像を構成するステップと、前記反射エコー信号に基づ 、て前記 被検体の生体組織の変位を計測するステップと、前記変位から歪み量或 1ヽは弾性率 を求め弾性画像を構成するステップと、前記断層像と前記弾性画像を表示させるス テツプと、前記弾性画像の情報に基づ!ヽて歪み状態或!ヽは弾性状態を特定するスコ ァリングするステップとを含むことを特徴とする超音波診断方法。
[23] 前記スコアリングするステップは、前記弾性画像の面積或いは形状に基づ 、て、歪 み状態或いは弾性状態を特定することを特徴とする請求項 19記載の超音波診断方 法。
[24] 前記断層像上の低エコー領域を第 1の領域とし、前記弾性画像上の硬!、部位の領 域を第 2の領域として設定して設定するステップを有し、前記スコアリングするステップ は前記第 1の領域と第 2の領域の比に基づいて歪み状態或いは弾性状態を特定する ことを特徴とする請求項 19記載の超音波診断方法。
PCT/JP2004/013140 2003-09-12 2004-09-09 超音波診断装置 WO2005025425A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US10/571,578 US8118746B2 (en) 2003-09-12 2004-09-09 Ultrasonic diagnostic apparatus
EP04787803.8A EP1665987B1 (en) 2003-09-12 2004-09-09 Ultrasonograph
JP2005513893A JP4685633B2 (ja) 2003-09-12 2004-09-09 超音波診断装置

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2003-322070 2003-09-12
JP2003322070 2003-09-12
JP2004-067983 2004-03-10
JP2004067983 2004-03-10

Publications (1)

Publication Number Publication Date
WO2005025425A1 true WO2005025425A1 (ja) 2005-03-24

Family

ID=34315671

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2004/013140 WO2005025425A1 (ja) 2003-09-12 2004-09-09 超音波診断装置

Country Status (4)

Country Link
US (1) US8118746B2 (ja)
EP (1) EP1665987B1 (ja)
JP (1) JP4685633B2 (ja)
WO (1) WO2005025425A1 (ja)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007312958A (ja) * 2006-05-25 2007-12-06 Hitachi Medical Corp 超音波診断装置
WO2007138881A1 (ja) 2006-05-25 2007-12-06 Hitachi Medical Corporation 超音波診断装置
EP1938754A1 (en) * 2005-10-19 2008-07-02 Hitachi Medical Corporation Ultrasonograph for creating elastic image
JP2008212522A (ja) * 2007-03-07 2008-09-18 Hitachi Medical Corp 超音波診断装置
WO2009063691A1 (ja) * 2007-11-16 2009-05-22 Hitachi Medical Corporation 超音波撮像システム
WO2009123208A1 (ja) * 2008-03-31 2009-10-08 株式会社 日立メディコ 超音波診断装置
JP2010193944A (ja) * 2009-02-23 2010-09-09 Konica Minolta Medical & Graphic Inc 超音波診断装置
JP2011200282A (ja) * 2010-03-24 2011-10-13 Ge Medical Systems Global Technology Co Llc 超音波診断装置及びその制御プログラム
CN1975370B (zh) * 2005-11-22 2012-07-04 通用电气公司 超声弹性模量计算和成像的方法
JP2013146625A (ja) * 2013-05-08 2013-08-01 Hitachi Medical Corp 超音波診断装置
JP2014030774A (ja) * 2013-11-22 2014-02-20 Super Sonic Imagine 関心領域の粘弾性の平均値を測定するための方法および装置
JP5438002B2 (ja) * 2008-06-03 2014-03-12 株式会社日立メディコ 医用画像処理装置及び医用画像処理方法
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
JP2016214868A (ja) * 2015-05-21 2016-12-22 東芝メディカルシステムズ株式会社 超音波診断装置及び画像診断装置
US10660614B2 (en) 2013-06-26 2020-05-26 Sony Corporation Ultrasonic processing apparatus and method

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1864612A4 (en) * 2005-03-30 2009-10-28 Hitachi Medical Corp ULTRASOUND DEVICE
JP4602825B2 (ja) * 2005-04-18 2010-12-22 オリンパスメディカルシステムズ株式会社 画像表示装置
EP1747757B8 (en) * 2005-07-27 2012-02-29 Samsung Medison Co., Ltd. Ultrasound system for displaying an elastic image
FR2899336B1 (fr) 2006-03-29 2008-07-04 Super Sonic Imagine Procede et dispositif pour l'imagerie d'un milieu viscoelastique
JP5479353B2 (ja) * 2008-10-14 2014-04-23 株式会社日立メディコ 超音波診断装置
KR101286331B1 (ko) * 2010-06-09 2013-07-15 삼성메디슨 주식회사 탄성 영상을 제공하는 초음파 시스템 및 방법
KR101512291B1 (ko) * 2013-05-06 2015-04-15 삼성메디슨 주식회사 의료 영상 장치 및 의료 영상 제공 방법
JP6594245B2 (ja) * 2016-03-30 2019-10-23 富士フイルム株式会社 硬さ導出装置、医用撮影システム、硬さ導出方法、及び硬さ導出プログラム
JP6726744B2 (ja) * 2016-07-19 2020-07-22 オリンパス株式会社 超音波観測装置、超音波観測装置の作動方法、及び超音波観測装置の作動プログラム
US20220148172A1 (en) * 2020-11-10 2022-05-12 Henry Ford Health System System and method for detecting medical conditions
CN114926389B (zh) * 2022-03-30 2024-07-09 什维新智医疗科技(上海)有限公司 一种超声图像的弹性信号提取装置
CN114881923A (zh) * 2022-03-30 2022-08-09 什维新智医疗科技(上海)有限公司 一种基于超声图像弹性信号的结节回声分析装置

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4201396B2 (ja) 1998-08-20 2008-12-24 株式会社日立メディコ 超音波診断装置
US6120445A (en) * 1998-10-02 2000-09-19 Scimed Life Systems, Inc. Method and apparatus for adaptive cross-sectional area computation of IVUS objects using their statistical signatures
WO2004017815A2 (en) * 2002-08-23 2004-03-04 The Board Of Regents Of University Of Texas System Novel risk assessment method based upon coronary calcification distribution pattern imaged by computed tomography
JP3932482B2 (ja) * 2002-10-18 2007-06-20 株式会社日立メディコ 超音波診断装置
JP2004246868A (ja) * 2003-01-22 2004-09-02 Canon Inc 画像抽出方法及び該方法をコンピュータに実行させるプログラム
US7257244B2 (en) * 2003-02-24 2007-08-14 Vanderbilt University Elastography imaging modalities for characterizing properties of tissue
CN1805713B (zh) * 2003-06-13 2010-04-28 松下电器产业株式会社 超声波诊断装置
US8123691B2 (en) * 2003-08-19 2012-02-28 Kabushiki Kaisha Toshiba Ultrasonic diagnostic apparatus for fixedly displaying a puncture probe during 2D imaging
CA2457376C (en) * 2003-10-14 2015-09-15 The University Of British Columbia Method for imaging the mechanical properties of tissue

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
HIGASHINO E. ET AL.: "Nyusen shikkan ni taisuru choonpa class bunrui", JOURNAL OF MEDICAL ULTRASONICS, vol. 30, no. SPECIAL EXTRA ISSUE, 15 April 2003 (2003-04-15), pages 125, XP002985734 *
ITO A. ET AL.: "Nyusen shikkan ni okeru elastography no risho oyo", DAI 11 KAI NIPPON NYUSEN KOJOSEN CHOONPA SHINDAN KAIGI SHOROKUSHU, 11 October 2003 (2003-10-11), pages 50 - 51, XP002985736 *
KANAI H. ET AL.: "Kekkanheki katasa bunpu no keihiteki hyokaho to denshi senshoku", RINSHOKU KENSA, vol. 51, no. 8, 25 August 2003 (2003-08-25), pages 805 - 812, XP002985733 *
See also references of EP1665987A4 *
T. SHIINA ET AL.: "Real time tissue elasticity imaging using the combined autocorrelation method", J. MED. ULTRASONICS, JAPAN, SOCIETY OF ULTRASONICS IN MEDICINE, TOKYO, JP, vol. 29, 21 September 2002 (2002-09-21), pages 119 - 128
W.-H. YEH ET AL.: "Elastic modulus measurements of human liver and correlation with pathology", ULTRASOUND IN MEDICINE AND BIOLOGY, NEW YORK, USA, vol. 28, no. 4, 1 April 2002 (2002-04-01), pages 467 - 474, XP004361097, DOI: doi:10.1016/S0301-5629(02)00489-1
YAMAKAWA M. ET AL.: "Nyusen shikkan no shoshiki dansei eizoho ni yoru chushutsu", JOURNAL OF MEDICAL ULTRASONICS, vol. 30, no. SPECIAL EXTRA ISSUE, 15 April 2003 (2003-04-15), pages 468, XP002985735 *
YAMAKAWA M. ET AL.: "Tissue elasticity reconstruction based on 3-dimensional finite-element model", JPN. J. APPL. PHYS., vol. 38, no. 5B, 30 May 1999 (1999-05-30), pages 3393 - 3398, XP002973311 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1938754A1 (en) * 2005-10-19 2008-07-02 Hitachi Medical Corporation Ultrasonograph for creating elastic image
JP5303147B2 (ja) * 2005-10-19 2013-10-02 株式会社日立メディコ 弾性画像を生成する超音波診断装置
EP1938754A4 (en) * 2005-10-19 2011-01-05 Hitachi Medical Corp ECHOGRAPH FOR CREATING AN ELASTIC IMAGE
CN1975370B (zh) * 2005-11-22 2012-07-04 通用电气公司 超声弹性模量计算和成像的方法
WO2007138881A1 (ja) 2006-05-25 2007-12-06 Hitachi Medical Corporation 超音波診断装置
JP2007312958A (ja) * 2006-05-25 2007-12-06 Hitachi Medical Corp 超音波診断装置
JP2008212522A (ja) * 2007-03-07 2008-09-18 Hitachi Medical Corp 超音波診断装置
JP4903271B2 (ja) * 2007-11-16 2012-03-28 株式会社日立メディコ 超音波撮像システム
WO2009063691A1 (ja) * 2007-11-16 2009-05-22 Hitachi Medical Corporation 超音波撮像システム
CN103349557A (zh) * 2008-03-31 2013-10-16 株式会社日立医疗器械 超声波诊断装置
US20110026800A1 (en) * 2008-03-31 2011-02-03 Akiko Tonomura Ultrasonic diagnostic apparatus.
JP2009240464A (ja) * 2008-03-31 2009-10-22 Hitachi Medical Corp 超音波診断装置
WO2009123208A1 (ja) * 2008-03-31 2009-10-08 株式会社 日立メディコ 超音波診断装置
US8718339B2 (en) 2008-03-31 2014-05-06 Hitachi Medical Corporation Ultrasonic diagnostic arrangements selecting parameters for selecting relevant estimation data for classifying an elasticity image
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
JP5438002B2 (ja) * 2008-06-03 2014-03-12 株式会社日立メディコ 医用画像処理装置及び医用画像処理方法
JP2010193944A (ja) * 2009-02-23 2010-09-09 Konica Minolta Medical & Graphic Inc 超音波診断装置
JP2011200282A (ja) * 2010-03-24 2011-10-13 Ge Medical Systems Global Technology Co Llc 超音波診断装置及びその制御プログラム
JP2013146625A (ja) * 2013-05-08 2013-08-01 Hitachi Medical Corp 超音波診断装置
US10660614B2 (en) 2013-06-26 2020-05-26 Sony Corporation Ultrasonic processing apparatus and method
JP2014030774A (ja) * 2013-11-22 2014-02-20 Super Sonic Imagine 関心領域の粘弾性の平均値を測定するための方法および装置
JP2016214868A (ja) * 2015-05-21 2016-12-22 東芝メディカルシステムズ株式会社 超音波診断装置及び画像診断装置

Also Published As

Publication number Publication date
EP1665987A1 (en) 2006-06-07
US8118746B2 (en) 2012-02-21
US20070112267A1 (en) 2007-05-17
EP1665987B1 (en) 2016-11-09
JPWO2005025425A1 (ja) 2007-11-08
EP1665987A4 (en) 2009-12-09
JP4685633B2 (ja) 2011-05-18

Similar Documents

Publication Publication Date Title
JP4685633B2 (ja) 超音波診断装置
JP4455003B2 (ja) 超音波診断装置
JP4314035B2 (ja) 超音波診断装置
JP5264097B2 (ja) 超音波診断装置
JP4657106B2 (ja) 超音波診断装置
JP5479353B2 (ja) 超音波診断装置
EP2263545B1 (en) Ultrasonograph
EP1554982A1 (en) Ultrasonographic device
US8948485B2 (en) Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, ultrasonic image processing program, and ultrasonic image generation method
US20110098563A1 (en) Ultrasonic diagnostic apparatus, ultrasonic image display method and ultrasonic diagnostic program
JP2007105400A (ja) 超音波診断装置及び画像処理装置
JP5016911B2 (ja) 超音波診断装置
KR101629541B1 (ko) 초음파 진단 장치 및 그 제어 프로그램
JP5473527B2 (ja) 超音波診断装置
JP5623609B2 (ja) 超音波診断装置
JP5802790B2 (ja) 超音波診断装置
JP4368185B2 (ja) 超音波診断装置
JP5680703B2 (ja) 超音波診断装置
JP5128149B2 (ja) 超音波診断装置
JP4601413B2 (ja) 超音波診断装置
JP6144990B2 (ja) 超音波画像撮像装置及び超音波画像撮像方法
JP2012228538A (ja) 超音波診断装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200480026263.4

Country of ref document: CN

AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BW BY BZ CA CH CN CO CR CU CZ DK DM DZ EC EE EG ES FI GB GD GE GM HR HU ID IL IN IS JP KE KG KP KZ LC LK LR LS LT LU LV MA MD MK MN MW MX MZ NA NI NO NZ PG PH PL PT RO RU SC SD SE SG SK SY TJ TM TN TR TT TZ UA UG US UZ VN YU ZA ZM

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): BW GH GM KE LS MW MZ NA SD SZ TZ UG ZM ZW AM AZ BY KG MD RU TJ TM AT BE BG CH CY DE DK EE ES FI FR GB GR HU IE IT MC NL PL PT RO SE SI SK TR BF CF CG CI CM GA GN GQ GW ML MR SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2005513893

Country of ref document: JP

REEP Request for entry into the european phase

Ref document number: 2004787803

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2004787803

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 2004787803

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2007112267

Country of ref document: US

Ref document number: 10571578

Country of ref document: US

WWP Wipo information: published in national office

Ref document number: 10571578

Country of ref document: US