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

超音波診断装置 Download PDF

Info

Publication number
WO2009131028A1
WO2009131028A1 PCT/JP2009/057470 JP2009057470W WO2009131028A1 WO 2009131028 A1 WO2009131028 A1 WO 2009131028A1 JP 2009057470 W JP2009057470 W JP 2009057470W WO 2009131028 A1 WO2009131028 A1 WO 2009131028A1
Authority
WO
WIPO (PCT)
Prior art keywords
elastic
ultrasonic
coupler
elastic coupler
diagnostic apparatus
Prior art date
Application number
PCT/JP2009/057470
Other languages
English (en)
French (fr)
Inventor
松村　剛
Original Assignee
株式会社 日立メディコ
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 株式会社 日立メディコ filed Critical 株式会社 日立メディコ
Priority to JP2010509144A priority Critical patent/JP5329532B2/ja
Priority to US12/989,527 priority patent/US8747320B2/en
Publication of WO2009131028A1 publication Critical patent/WO2009131028A1/ja

Links

Images

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/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6843—Monitoring or controlling sensor contact pressure
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06—Measuring blood flow
    • 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/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
    • A61B8/4455—Features of the external shape of the probe, e.g. ergonomic aspects
    • 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/488—Diagnostic techniques involving Doppler signals

Definitions

  • the present invention relates to an ultrasonic diagnostic apparatus, in particular, at the time of examination such as a tomographic image (B mode image) or elasticity information of a biological tissue of a subject, blood flow diagnosis by Doppler measurement or color flow mode (CFM), etc.
  • the present invention relates to a technique capable of detecting an absolute pressure applied to an object by an ultrasonic probe.
  • an ultrasonic diagnostic apparatus transmits an ultrasonic wave from an ultrasonic probe (hereinafter simply referred to as a probe) into a living body, which is a subject, and is a reflected echo signal of an ultrasonic wave reflected from the living body. Is received by a probe, and based on the received reflected echo signal (RF signal), an image suitable for examination of a tissue or function in a living body is generated and diagnosed.
  • a probe an ultrasonic wave from an ultrasonic probe
  • RF signal received reflected echo signal
  • the inspection using the B mode image presses the living tissue by pressing the probe against the subject with a relatively strong force.
  • blood flow tests such as Doppler measurement and CFM, if the probe is pressed too hard and the living tissue is compressed, the blood vessel cross-section is collapsed and correct blood flow information cannot be obtained. It is desirable to carry out the inspection under a tight compression condition.
  • the absolute pressure applied to the living tissue may be hindered from promptly making an appropriate diagnosis if the examination proceeds in an inappropriate compression state when switching between examination methods or examination methods. It is desirable to measure and display in real time.
  • Patent Document 1 discloses a pair of RF signal frame data having different acquisition times in order to measure an actual pressure applied to a living tissue of a subject, that is, an absolute pressure applied to the living tissue (hereinafter simply referred to as an absolute pressure). Based on the above, it has been proposed to obtain the strain of the elastic coupler using a known displacement / strain calculation, and to convert the obtained strain into an absolute pressure using the elastic modulus of a preset elastic coupler.
  • the absolute pressure applied to the subject can be measured, and the elasticity information of the hardness or softness of the living tissue can be acquired.
  • the boundary between the body surface of the subject is appropriate In some cases, the thickness of the elastic coupler cannot be detected properly. In addition, if there is contact unevenness, the deformation of the contact part around the non-contact part will increase, so if you measure the deformation or distortion of the boundary part and convert it to absolute pressure, you can not measure the appropriate absolute pressure distribution Sometimes.
  • the problem to be solved by the present invention is to enable detection of absolute pressure with high accuracy.
  • another object is to improve the usability by simplifying the operation of detecting the absolute pressure using an appropriate elastic coupler according to the inspection method or the like.
  • a first aspect of the present invention is an ultrasonic probe that transmits and receives an ultrasonic wave in contact with a subject, a transmission unit that drives the ultrasonic probe, and the ultrasonic probe that is received by the ultrasonic probe.
  • an ultrasonic diagnostic apparatus comprising: a reception unit that receives and processes an RF signal that is a reflected echo signal; and an image generation unit that generates an ultrasonic image based on the RF signal output from the reception unit.
  • a pressure calculation unit that obtains a pressure applied to the subject based on deformation of an elastic coupler mounted on an ultrasonic transmission / reception surface of the ultrasonic probe, and the elastic coupler has at least two layers;
  • the pressure calculation unit detects a position of a boundary surface between the two layers, and changes a position of the boundary surface based on a detection position of the boundary surface and an initial position of the boundary surface obtained in advance.
  • the absolute pressure is obtained based on the position change and the elastic characteristics of the elastic coupler.
  • the elastic coupler is formed with at least two layers having different ultrasonic reflection characteristics due to an elastic material having flexibility, and a boundary surface of the layers having different ultrasonic reflection characteristics is attached to the ultrasonic transmission / reception surface.
  • the pressure calculation unit is based on the RF signal output from the reception unit, The position of the boundary surface in the thickness direction of the elastic coupler is detected, the position change of the boundary surface is obtained based on the detection position of the boundary surface and the initial position of the boundary surface obtained in advance, and the position change and the The absolute pressure applied to the subject is obtained based on the set elastic characteristic of the elastic coupler.
  • the position of the boundary surface formed inside the elastic coupler is detected, and the absolute pressure is obtained based on the change in position, so the boundary between the elastic coupler and the body surface of the subject is detected.
  • the absolute pressure can be detected with high accuracy because the inner boundary surface can be stably detected regardless of the contact between the elastic coupler and the body surface of the subject.
  • the initial position of the boundary surface in the initial state where the compression force is not applied to the elastic coupler is the RF output from the detection reception unit in the same manner as the position detection of the boundary surface in the pressurized state where the compression force is applied. Detection can be based on the signal.
  • the characteristic configuration of the second aspect of the present invention is the same as that of the first aspect in the elastic coupler, but the configuration of the pressure calculation unit is different. That is, the pressure calculation unit detects the displacement of the position of the boundary surface in the thickness direction of the elastic coupler based on a pair of the RF signal frame data having different acquisition times output from the reception unit, and the displacement The strain in the thickness direction of the boundary surface is obtained based on the above, and the absolute pressure applied to the subject is obtained based on the strain in the thickness direction and a predetermined elastic property of the elastic coupler.
  • the internal boundary surface can be stably detected, so that the distortion of the boundary surface can be detected with high accuracy, and the distortion and elastic characteristics can be detected. Based on the relationship, the absolute pressure can be accurately detected.
  • the pressure calculation unit accumulates strain in the thickness direction of the boundary surface over time from the initial state where the elastic coupler is not pressurized, and sets the strain integrated value in advance.
  • the absolute pressure applied to the subject can be obtained based on the elastic characteristics of the elastic coupler. According to this, since the strain in the thickness direction of the boundary surface is integrated, the strain on the boundary surface can be detected with higher accuracy, and the absolute pressure can be detected with higher accuracy based on the relationship between the integrated strain value and the elastic characteristic. It can be detected.
  • the elastic coupler is formed with a thin intermediate layer sandwiched between the boundary surfaces of the two layers, and the ultrasonic reflection characteristics of the intermediate layer are different from those of the other two layers.
  • the intermediate layer includes a plurality of linear ultrasonic reflectors extending in a direction perpendicular to the scanning direction of the ultrasonic beam of the ultrasonic probe and spaced apart from each other. can do.
  • the pressure calculation unit detects the position of the intermediate layer in the thickness direction, and sets the detection position of the intermediate layer and the initial position of the intermediate layer determined in advance. On the basis of this, it is possible to obtain a position change of the intermediate layer, and to obtain an absolute pressure applied to the subject based on the position change and a preset elastic characteristic of the elastic coupler.
  • the pressure calculation unit is configured so that the intermediate of the elastic coupler is based on a pair of the RF signal frame data having different acquisition times output from the reception unit. A displacement in a thickness direction of the layer is detected, a strain in the thickness direction of the intermediate layer is obtained based on the displacement, and the subject based on the strain in the thickness direction and a predetermined elastic characteristic of the elastic coupler It can be set as the structure which calculates
  • the elastic coupler is formed to have a two-layer structure, the layer on the contact surface side in contact with the body surface of the subject is thin, and the ultrasonic reflection characteristics are strongly formed. can do.
  • the elastic coupler can be configured to have a large ultrasonic attenuation characteristic by mixing an ultrasonic scatterer in the elastic material.
  • the elastic coupler of the first or second aspect In order to simplify the operation of detecting the absolute pressure using an appropriate elastic coupler according to the inspection method, etc., and to solve the problem of improving usability, the elastic coupler of the first or second aspect,
  • the position in the thickness direction of the boundary surface of the two layers is made different according to the type of the elastic coupler, and an identification code that can identify the type of the elastic coupler is formed by the RF signal.
  • the pressure calculation unit detects the identification code based on the RF signal or the RF signal frame data, identifies the type of the elastic coupler, and is set according to the type of the elastic coupler. It can be set as the structure which calculates
  • the ultrasonic probe In the mounted state the type of elastic coupler can be automatically identified by the pressure calculation unit. Therefore, if an elastic characteristic corresponding to the type of the elastic coupler is stored in the pressure calculation unit in advance, the absolute pressure can be calculated in accordance with the elastic characteristic of the type of the elastic coupler.
  • the identification code may be formed by changing the position in the thickness direction of the boundary surface of the two layers.
  • the pressure calculation unit may be configured to identify the type of the elastic coupler based on the RF signal or the RF signal frame data based on the depth distribution pattern of the RF signal in the coupler echo area.
  • the identification code may be formed by encoding and dispersing ultrasonic scatterers in at least one of the scanning direction and the thickness direction in the region portions at both ends of the elastic coupler in the scanning direction.
  • the pressure calculation unit can be configured to identify the type of the elastic coupler based on the RF signal or the RF signal frame data based on the RF signal pattern in the region at both ends of the elastic coupler.
  • the elastic properties set in the coupler database are the elastic modulus, the relationship curve between the deformation in the thickness direction and the elastic modulus, the relationship curve between the strain in the thickness direction and the elastic modulus, the integrated value of the deformation or strain in the thickness direction and the elastic modulus. It may be at least one of a modulus of elasticity correction coefficient for a relationship curve, deformation or strain in the thickness direction.
  • absolute pressure can be detected with high accuracy.
  • the operation is simple and the usability is improved.
  • FIG. 1 is a block configuration diagram of an ultrasonic diagnostic apparatus according to an embodiment of the present invention.
  • the block block diagram which shows the detailed structure of the pressure calculating part of embodiment.
  • FIG. 3 is a diagram showing an ultrasonic image of the elastic coupler of FIG.
  • FIG. 1 shows a functional block configuration diagram of an embodiment in which the present invention is applied to an ultrasonic diagnostic apparatus that performs an examination based on elasticity information.
  • the ultrasonic diagnostic apparatus of this embodiment obtains an elastic image representing the hardness or softness of a living tissue while obtaining a tomographic image of a diagnostic region of a subject using ultrasonic waves. is there.
  • an ultrasonic probe hereinafter simply referred to as a probe
  • the probe 1 is formed by arranging a large number of transducers which are ultrasonic wave generation sources and receive reflected echoes in a strip shape.
  • Each transducer generally has a function of converting an input pulse wave or continuous wave ultrasonic signal into an ultrasonic wave and emitting it, and a reflected echo emitted from the inside of the subject 100 into an electric echo signal. It has the function of converting and outputting.
  • the transmission circuit 2 and the reception circuit 3 are cooperatively controlled by the ultrasonic transmission / reception control circuit 4.
  • the ultrasonic transmission / reception control circuit 4 controls the timing of transmitting and receiving ultrasonic waves so that an ultrasonic wave transmission beam suitable for the intended ultrasonic inspection is driven from the probe 1 to the subject 100.
  • the reception circuit 3 is controlled to receive a desired ultrasonic reception beam from the probe 1.
  • the ultrasonic transmission / reception control circuit 4 of the present embodiment scans the ultrasonic transmission beam along the tomographic plane so as to form an ultrasonic transmission beam and an ultrasonic reception beam suitable for measurement of a B-mode tomographic image.
  • the transmission circuit 2 and the reception circuit 3 are controlled.
  • the receiving circuit 3 amplifies the reflected echo signal received by the probe 1 with a predetermined gain.
  • the number of reflected echo signals corresponding to the number of amplified transducers is input to the adder circuit 5.
  • the adder circuit 5 adds the phases of the plurality of reflected echo signals amplified by the receiving circuit 3 together and generates RF signal frame data corresponding to the tomographic plane.
  • the signal processing unit 6 receives the RF signal frame data output from the adder circuit 5 and performs various signal processing such as gain correction, log correction, detection, contour enhancement, and filter processing to generate image data.
  • the black and white scan converter 7 acquires image data output from the signal processing unit 6 at a sound wave cycle, and controls tomographic scanning means and a system for reading out the ultrasonic image at a television system cycle.
  • an A / D converter that converts image data from the signal processing unit 6 into a digital signal
  • a plurality of frame memories that store the image data digitized by the A / D converter in time series
  • a controller for controlling these operations.
  • the time-series image data of, for example, the B mode image generated by the monochrome scan converter 7 is output to the image display 9 via the switching adder 8.
  • the image display 9 is a D / A converter that converts the image data output from the black and white scan converter 7 into an analog signal, and a color television monitor that displays the analog video signal output from the D / A converter as an image. It consists of.
  • the RF signal frame data generated by the adder circuit 5 is input to the RF signal frame data selection unit 10.
  • the RF signal frame data selection unit 10 sequentially stores the RF signal frame data sequentially output from the adder circuit 5 in the frame memory. Then, according to a command from a control unit (not shown), for example, one RF signal frame data is selected from the latest RF signal frame data and temporally past RF signal frame data, and a displacement / distortion calculation unit In FIG. 11, a pair of RF signal frame data having different acquisition times is output.
  • the RF signal may be a signal in the form of a code demodulated I or Q signal.
  • the displacement / distortion calculation unit 11 performs one-dimensional or two-dimensional correlation processing on the pair of RF signal frame data output from the RF signal frame data selection unit 10, and scans the ultrasonic beam of the RF signal frame data and Displacement (or displacement vector) is calculated for a plurality of measurement points i and j set in the depth direction.
  • the calculated displacement data of the plurality of measurement points is generated as displacement frame data.
  • a well-known block matching method or gradient method can be used as the displacement calculation method as described in Patent Document 1.
  • the block matching method sets a target block consisting of multiple pixels centered on the pixel of the target measurement point, and moves the block whose image information of the target block is approximate to the pre-displacement frame This is a method of searching and assuming that the position of the closest block is displaced to the current position.
  • the displacement / strain calculator 11 uses the displacement frame data to calculate the strain at each measurement point by spatially differentiating the displacement at each measurement point, as is well known. Strain frame data is generated from the strain at each measurement point obtained by the calculation and output to the elastic modulus calculation unit 12.
  • the elastic modulus calculation unit 12 is based on the absolute pressure (stress) acting on each measurement point i, j given from the pressure calculation unit 30 described later based on the strain at each measurement point of the input strain frame data. Then, the elastic modulus (for example, Young's modulus) at each measurement point is calculated by a known method.
  • Elastic modulus frame data is generated from the elastic modulus of each measurement point obtained by the calculation, and is output to the elastic data processing unit 13.
  • the elastic data processing unit 13 performs, for example, smoothing processing within the frame, contrast optimization processing, or smoothing processing in the time axis direction between frames on the elastic frame data of the strain or elastic modulus generated by the elastic modulus calculation unit 12. Various image processing is performed and output to the color scan converter 14.
  • the color scan converter 14 generates elastic image data by adding hue information such as red, green, and blue to pixels corresponding to each measurement point based on the elastic frame data output from the elastic data processing unit 13. For example, in the strain frame data output from the elasticity data processing unit 13, for a pixel with a large distortion measured, the pixel is converted into a red code in the elasticity image data, and conversely for a pixel with a small distortion measured. Converts the pixel into a blue code in the elastic image data. Similarly, in the case of elastic modulus frame data, hue information is assigned to each pixel to generate elastic image data. Note that, instead of the color scan converter 14, a black and white scan converter may be used to provide a luminance hierarchy corresponding to the magnitude of distortion or the like.
  • the elastic image data generated by the color scan converter 14 is output to the switching addition unit 8.
  • the switching addition unit 8 inputs black and white tomographic image data output from the black and white scan converter 7 and color elastic image data output from the color scan converter 14, and adds or adds both images according to an input command. Switch to output to the image display 9. Specifically, only the black and white tomographic image data, only the color elastic image data, or both the images can be arranged and displayed on the image display 9. In addition, the output image can be switched in response to an input command, for example, by adding and synthesizing both image data, that is, generating a translucent superimposed image and displaying it on the image display 9. .
  • FIG. 2 shows an embodiment of the elastic coupler 20.
  • the elastic coupler 20 is formed in a plate shape from a gel-like substance that is a flexible elastic material, and a collar portion 21 is formed around one surface.
  • the material of the gel-like substance that forms the elastic coupler 20 has a low ultrasonic attenuation, a sound velocity and an acoustic impedance, such as an acoustic coupling material and an acoustic lens material. And those excellent in acoustic coupling characteristics with a living body are preferable.
  • the elastic coupler 20 of the present embodiment is formed to have two layers 20A and 20B having different ultrasonic reflection characteristics as shown in the figure, and the boundary surface 22 of these layers 20A and B is the probe 1.
  • the ultrasonic wave transmitting / receiving surface is provided between the mounting surface mounted on the ultrasonic transmission / reception surface and the contact surface in contact with the body surface of the subject 100.
  • the elastic coupler 20 when the elastic coupler 20 is attached to the ultrasonic transmission / reception surface of the probe 1, even if the elastic coupler 20 is in contact with an uneven portion such as the contact surface with the ultrasonic transmission / reception surface and the body surface of the subject 100, the elastic coupler 20 The material which does not produce is preferable. That is, when a gap is generated between the ultrasonic transmission / reception surface and the body surface, the ultrasonic wave irradiated from the probe 1 is reflected at the boundary of the air in the gap and becomes noise of the ultrasonic image.
  • gelatin for example, gelatin, agar, oil gel, aqueous gel (hydrogel) composed of an aqueous polymer such as acrylamide and polyvinyl alcohol and water, polyurethane, oil gel obtained by crosslinking a composition containing rubber and an oil component, and plasticizer for raw rubber
  • a rubber obtained by molding and crosslinking a composition containing a low molecular weight rubber can be used.
  • the elastic coupler 20 formed in this way is used by being mounted on the ultrasonic transmission / reception surface 24 of the probe 1 using the mounting tool 23. That is, as shown in FIG. 3 (A), the mounting tool 23 has a groove 25 formed in accordance with the flange portion 21 of the elastic coupler 20, and is formed in a frame shape with resin or the like. A locking claw 26 that engages with the outer peripheral portion of the ultrasonic wave transmitting / receiving surface 24 is formed inside the frame of the mounting tool 23. Further, a pair of gripping pieces 27 are formed on the opposing long sides of the frame, and an opening 28 into which the elastic coupler 20 is inserted is formed on the top surface. As shown in FIG.
  • the elastic coupler 20 is inserted into the opening 28 from the bottom side of the mounting tool 23 formed in this way, and the flange portion 21 of the elastic coupler 20 is pushed into the groove 25 of the frame body.
  • the probe 1 is mounted on the ultrasonic transmission / reception surface 24 and used. At this time, so that there is no gap between the contact surface between the lower surface of the elastic coupler 20 and the ultrasonic transmission / reception surface 24, it is applied to the ultrasonic transmission / reception surface 24 of the probe 1 with jelly or the like interposed therebetween, and the locking claw 26 is Hook it on the outer periphery of the sound wave transmitting / receiving surface 24 and fix it firmly.
  • the grasping portion 1A of the probe 1 to which the elastic coupler 20 is mounted is grasped, the exposed surface of the elastic coupler 20 is pressed against the body surface of the subject 100, and an ultrasonic examination is performed.
  • the subject pressing mechanism 18 holds the grip 1A of the probe 1 so that the subject pressing mechanism 18 applies a pressing force to the body surface of the subject 100. it can.
  • the pressure calculation unit 30 includes the blocks shown in FIG. That is, a coupler ID identification unit 31, an initialization processing unit 32, a coupler distortion calculation unit 33, a pressure conversion unit 34, a compression state image construction unit 35, and a coupler database 36 are configured.
  • the elastic coupler 20 applied to the present embodiment has the same shape as that shown in FIG. 2 (A), but for each type of the elastic coupler 20 depending on the application so that the type of the elastic coupler 20 can be automatically identified.
  • the same type of elastic coupler 20 with an identification code (ID) is used by making the depth of the boundary surface 22 or the reflection characteristics of the layers 20A and B different.
  • the coupler ID identification unit 31 captures the RF signal frame data from the RF signal frame data selection unit 10, detects the presence or absence of the boundary surface 22, and detects that the elastic coupler 20 is attached to the probe. Furthermore, the depth distribution of the RF signal or the like in the coupler echo area is determined to detect the depth of the boundary surface 22, and the ID code of the elastic coupler 20 is automatically identified with reference to the coupler database 36. The identified ID code is output to the pressure converter 34. In the coupler database 36, an ID code is set corresponding to the depth pattern of the boundary surface 22 and the position pattern in the depth direction by an input means (not shown), and further elastic characteristics are set corresponding to the ID code. Has been.
  • the initialization processing unit 32 detects that the elastic coupler 20 is in a non-pressurized initial state based on the intensity change of the RF signal of the RF signal frame data, and the elastic coupler 20 is based on the RF signal in the initial state. The initial thickness of the elastic coupler 20 is obtained.
  • the coupler distortion calculation unit 33 detects the position in the thickness direction of the boundary surface 22 of the elastic coupler 20 based on the RF signal frame data, and detects the change in the thickness direction position of the boundary surface 22, The displacement and distortion are obtained. The obtained displacement or distortion is output to the pressure conversion unit 33 together with the ID code.
  • the pressure converter 34 reads the elastic characteristic corresponding to the input ID code from the coupler database 36, and converts the displacement or strain input from the coupler distortion calculator 33 into an absolute pressure based on the read elastic characteristic. It has become.
  • the compression state image construction unit 35 constructs a compression state image in order to display the absolute pressure output from the pressure conversion unit 34 on the image display 9.
  • the inspection mode is, for example, an inspection using a B-mode tomographic image.
  • the elastic coupler 20 is attached to the ultrasonic wave transmitting / receiving surface 24 of the probe 1 by manual operation as shown in FIG. Detection of whether or not the elastic coupler 20 is attached can be automatically recognized based on whether or not the presence / absence of the boundary surface 22 can be detected by capturing the RF signal frame data in the coupler ID identification unit 31 and the intensity distribution of the RF signal. In addition, the operator may input that the elastic coupler 20 has been attached to the coupler ID identification unit 31 via an input unit (not shown).
  • the coupler ID identifying unit 31 takes in the RF signal frame data from the RF signal frame data selecting unit 10, detects the presence of the boundary surface 22, and then determines the depth distribution of the RF signal or the like in the coupler echo region to determine the boundary surface 22 And the ID code of the elastic coupler 20 is automatically identified with reference to the coupler database 36. The identified ID code is output to the pressure converter 34.
  • FIG. 6 shows an example of a B-mode image 101 obtained by attaching the elastic coupler 20. As shown in the figure, echo images of the two layers 20A and 20B and the boundary surface 22 of the elastic coupler 20 appear in the upper part of the image. Therefore, by detecting the position where the RF signal in the coupler echo region changes rapidly, the position can be detected on the boundary surface 22 and the ID code can be automatically identified.
  • the pressure calculation unit 30 is activated when the coupler ID identification unit 31 detects the ID code of the elastic coupler 20.
  • the initial position D (0) in the depth direction of the boundary surface 22 of the elastic coupler 20 can be measured in advance, but the liquid component of the gel-like substance that is the material of the elastic coupler 20 may evaporate, etc. It is conceivable that the initial position D (0) of 22 becomes small. Therefore, in order to accurately detect the absolute pressure, it is desirable to automatically measure the initial position D (0) of the boundary surface in the initial state every time the ultrasonic inspection is performed.
  • the initialization processing unit 32 searches for an operator.
  • the initial state in which the toucher 1 is held and held in the air can be automatically recognized.
  • This initial state without pressure can be detected by a multiple echo signal included in the RF signal. That is, multiple echoes are generated due to a sudden change in acoustic impedance at the boundary between the exposed surface of the acoustic lens and air or the boundary between the exposed surface of the elastic coupler 20 and air.
  • an acoustic lens (thickness of about 1 mm) is placed in the depth region close to the ultrasonic transmission / reception surface. ) Is received, which includes periodic multiple echoes originating from.
  • the elastic coupler 20 for example, about 5 to 10 mm thick
  • the multiple echoes derived from the acoustic lens disappear and the multiple echoes of a relatively long period derived from the elastic coupler 20 are received. Therefore, it is possible to automatically recognize that the elastic coupler 20 is maintained in the non-pressurized initial state by detecting the presence / absence and period of multiple echoes.
  • the operator can hold the probe 1 and hold it in the air, and input the interface from the input means (not shown) to the coupler thickness calculator 32 manually. Automatic measurement of 22 initial positions D (0) can be commanded.
  • the initialization processing unit 32 automatically measures the initial position D (0) of the boundary surface when the initial state of the elastic coupler 20 is recognized or when a command is input.
  • the initialization processing unit 32 takes in the luminance data or signal intensity of the B-mode image output from the signal processing unit 6 instead of the RF signal output from the RF frame signal selection unit 10, and performs similar processing.
  • the initial position D (0) of the boundary surface of the elastic coupler 20 can be detected.
  • the time ti (0) when the intensity of the RF signal changes greatly is obtained, and the initial position D (0) of the boundary surface is detected.
  • ti (0) is the round-trip time of the ultrasonic wave reflected from the boundary surface of the elastic coupler 20
  • multiplying ti (0) by the speed of sound C and dividing by 2 results in the boundary surface being a one-dimensional distribution in the scanning direction.
  • the initial position Di (0) can be obtained.
  • Di (0) can be obtained by obtaining a period of multiple echoes (period above the threshold) T in which the intensity of the RF signal greatly changes and multiplying 1/2 of T by the sound speed C.
  • the initial position Di (0) of the boundary surface 22 of the elastic coupler 20 can be obtained as a two-dimensional distribution taking into account the distribution in the direction orthogonal to the scanning direction.
  • S7 Coupler deformation (distortion) measurement
  • the process in S7 is a process in the coupler distortion calculation unit 33.
  • the coupler distortion calculation unit 33 can detect that the elastic coupler 20 is in a pressurized state by the disappearance of the multiple echo caused by the elastic coupler 20 described above.
  • Coupler deformation (distortion) measurement method 1 The coupler distortion calculator 33 detects the depth position of the boundary surface 22 of the elastic coupler 20 based on the RF signal at an arbitrary time t in the pressurized state of the elastic coupler 20, and scans the elastic coupler 20 in the pressurized state.
  • the boundary position distribution Di (t) in the direction is obtained.
  • the intensity of the RF signal Qi (t) changes greatly.
  • the boundary surface position distribution Di (t) in the scanning direction orthogonal to the ultrasonic beam is obtained.
  • the boundary position change distribution ⁇ Di (t) is obtained by the following equation (1)
  • the total strain distribution Si (t) in the scanning direction of the boundary surface 22 of the elastic coupler 20 is obtained by the following equation (2).
  • Coupler deformation (distortion) measurement method 2 the coupler strain calculation unit 33 can measure the strain of the elastic coupler 20 and obtain the absolute pressure Pi (t) based on the measured strain. That is, the coupler distortion calculation unit 33 includes the RF signal frame data Qij (0) in the coupler echo area in the initial state output from the RF signal frame data selection unit 10, and the RF signal frame data at an arbitrary time t in the pressurized state.
  • i is a coordinate in the scanning direction of the elastic coupler 20
  • j is a coordinate in the thickness direction (depth direction) of the elastic coupler 20.
  • the coupler distortion calculation unit 33 takes in the RF signal Qij (t) of the coupler echo region that changes in real time in the pressurized state, and based on Qij (0) in the initial state and Qij (t) at an arbitrary time t Then, the displacement frame data is generated by obtaining the displacement of each measurement point i, j by a known displacement calculation method. Then, the displacement frame data is spatially differentiated to obtain distortion frame data composed of the total distortion amounts Sij (t) at the respective measurement points i and j.
  • the average value Si * () of the total distortion amount Sij (t) over the entire range of the depth j of the coupler echo area. t) is obtained, and the average value Si * (t) is evaluated as the total distortion amount Sij (t) at the coordinate position i in the ultrasonic scan direction.
  • the absolute pressure Pi (t) is obtained by the pressure conversion unit 38 using the evaluated total strain amount Sij (t).
  • the processing for obtaining the absolute pressure Pi (t) can be executed simultaneously with the calculation of the strain of the living tissue in the elasticity calculation unit 12 of FIG.
  • Coupler deformation (distortion) measurement method 3 In the coupler compression evaluation method 3, a pair of RF signal frame data having different acquisition times is acquired from the RF signal frame data selection unit 10 continuously from the initial state for each measurement time point in the pressurized state. And every time a pair of RF signal frame data is acquired, the strain change ⁇ Sij (tk),..., ⁇ Sij (t) of each measurement point i, j is calculated for the entire region up to the boundary surface 22 of the elastic coupler 20. Ask. Then, for example, distortion changes ⁇ Si, j (t) at time (t ⁇ 1) and time (t) that are adjacent in time are obtained.
  • the distortion change ⁇ Si, j (tk),..., ⁇ Sij (t) is sequentially accumulated for a pair of time-sequential RF signal frame data, and the distortion change integrated value ⁇ Sij (t) at the current time is sequentially accumulated.
  • ⁇ Sij (t) at the measurement points i and j is averaged in the direction of the coordinate j in the coupler echo area, and a distortion change integrated value ⁇ Si * (t) at the measurement point i is obtained.
  • the pressure conversion unit 34 reads the elastic characteristic (for example, Young's modulus E) corresponding to the ID of the elastic coupler 20 identified by the coupler ID identification unit 31 from the coupler database 36, and obtained by the coupler distortion calculation unit 33 by the method 1.
  • the strain Si (t) is converted into the absolute pressure distribution Pi (t) applied to the living tissue of the subject 100 by the following equation (3).
  • Pi (t) Si (t) ⁇ E (3)
  • the current absolute pressure distribution Pi (t) obtained by the conversion is output to the elastic modulus calculator 12 in FIG.
  • the elastic modulus calculation unit 12 performs the measurement of each measurement point i, j based on the strain ⁇ ij (t) obtained for each measurement point i, j of the living tissue by a known calculation process.
  • the elastic modulus (for example, Young's modulus) Eij (t) is obtained by the following equation (4) and output to the elastic data processing unit 13.
  • the total distortion amount Sij (t) obtained based on the average value Si * (t) of the method 2 by the coupler distortion calculation unit 33 or the distortion change integrated value ⁇ Si * (t) obtained by the method 3 is expressed by the following equation: Convert to absolute pressure distribution Pi (t) by (5) and (6).
  • the pressure converter 34 may be formed integrally with the coupler distortion calculator 33.
  • the compression state image constructing unit 35 images the absolute pressure Pi (t) obtained in the absolute pressure conversion process of S8 and displays the image on the image display 9, so that the operator can perform the same while performing the ultrasonic examination. In this screen, it is possible to immediately determine whether or not the compression state is appropriate according to the inspection item.
  • the compression state image constructing unit 38 constructs at least one image such as a numerical display, a graph display of change over time, a bar chart display, and the like by using the color scan converter 14 to generate a color image, for example, the absolute pressure Pi (t)
  • the data is converted into data and displayed side by side or superposed on a part of the ultrasonic image displayed on the image display 9.
  • FIG. 7 shows an example in which the compression state image is displayed on the elastic image side by side.
  • an elastic image 110 is displayed at the center of the screen, and an elastic modulus E (kPa) and an absolute pressure ⁇ (kPa) are displayed numerically in the display window 111 near a rectangular region of interest (ROI).
  • a bar chart 112 is displayed in which the absolute pressure ⁇ at the current time is an average value of the absolute pressure Pi.
  • a graph 113 of the absolute pressure distribution Pi in the scanning direction is displayed at the top of the screen, and a graph 114 of the average value of the absolute pressure Pi over time is displayed at the bottom.
  • a color bar 115 with an elastic modulus E is displayed at the right end of the screen.
  • the operator can make a diagnosis by observing the image of FIG. 7 to evaluate the elasticity of the ROI living tissue under an appropriate absolute pressure.
  • FIG. 8 shows another example in which the compression state image is displayed on the elastic image so as to overlap or be displayed side by side.
  • an elastic image 110 is displayed at the center of the screen, and elastic modulus E (kPa) and absolute pressure ⁇ (kPa) are displayed as bar charts 116 and 117 near a rectangular region of interest (ROI).
  • a bar chart 118 representing the average value of the thickness change ⁇ Di (t) of the elastic coupler from the initial state is displayed.
  • the boundary between the elastic coupler and the subject is displayed by a plurality of points 119 at the top of the screen.
  • the absolute pressure which is the calculation result of the pressure calculation unit 30
  • the reference pressure range is compared, and the reference pressure range is exceeded, for example, the display color of the boundary line indicating the ROI is displayed. By changing or blinking, it is possible to easily identify whether or not the compression state is appropriate.
  • the absolute pressure applied to the living tissue of the subject 100 by the probe 1 can be accurately measured in real time in the normal elasticity image measurement processing.
  • the ID code of the elastic coupler can be read during the elastic image measurement process, even if various different elastic couplers are arbitrarily replaced, the elastic coupler is automatically identified and the elastic characteristics of the elastic coupler are detected. Therefore, the absolute pressure can be measured with high accuracy, and the labor of the operator can be saved and the usability can be improved.
  • an ultrasonic diagnostic apparatus that generates and displays an elastic image has been described as an example.
  • the present invention is not limited to this, and diagnosis based on a tomographic image (B-mode image) of a living tissue of a subject is performed.
  • the present invention can be applied to an ultrasonic diagnostic apparatus for performing tests such as blood flow diagnosis by Doppler measurement or color flow mode (CFM). Thereby, the evaluation of the compression state suitable for various inspections can be realized with high accuracy.
  • FIG. 9A is a perspective view of the elastic coupler 20 according to the second embodiment.
  • the same attachment tool as that shown in FIG. 3 can be used.
  • the second embodiment is different from the elastic coupler 20 of the first embodiment shown in FIG. 2 (A) in that the ultrasonic reflection characteristics are thin on the boundary surface between the two layers 20A and 20B. This is because the intermediate layer 20C is formed therebetween.
  • the ultrasonic reflection characteristics of the intermediate layer 20C are different from those of the other two layers 20A and 20B.
  • the reflection intensity of the intermediate layer 20C may be higher or lower than the other two layers 20A and 20B.
  • An example of the B mode image 101 in this case is shown in FIG.
  • an echo image of the intermediate layer 20C appears at the boundary between the two layers 20A and B of the elastic coupler 20 at the top of the image.
  • the coupler ID identifying unit 31 can identify the ID code of the elastic coupler based on the distribution in the thickness direction of the RF signal, as in the first embodiment.
  • the coupler strain calculation unit 33 adjusts the position of the intermediate layer 20C in the thickness direction, the strain in the thickness direction of the intermediate layer 20C, or the elastic coupler according to the passage of time from the initial state in which no pressure is applied.
  • the integrated value of the strain obtained by integrating the strain in the thickness direction of the intermediate layer 20C can be obtained to obtain the absolute pressure applied to the subject.
  • FIG. 10A shows a perspective view of the elastic coupler 20 according to the third embodiment.
  • the same attachment tool as that shown in FIG. 3 can be used.
  • the difference of the third embodiment from the elastic coupler 20 of the second embodiment shown in FIG. 9A is that, instead of the intermediate layer 20C, the probe 1 is extended in a direction perpendicular to the scanning direction of the ultrasonic beam.
  • a plurality of linear ultrasonic reflectors 20D are provided at intervals.
  • the ID code of the elastic coupler can be identified based on the position in the depth direction of the ultrasonic reflector 20D. Based on the change in position, the absolute pressure can be detected.
  • the layer on the contact surface side that is in contact with the body surface of the subject 100 (20A in the example of FIG. 2 (A)) is thin, and Ultrasonic reflection characteristics can be strongly formed.
  • the ID encoding may be performed by dispersing and mixing the ultrasonic scatterer throughout the elastic coupler 20 and changing the dispersion concentration of the scatterer.
  • the ultrasonic scatterer a material such as graphite powder or polyethylene powder having an acoustic impedance different from that of the elastic coupler 20 can be used.
  • At least one of the end portions in the scanning direction of the elastic coupler 20 that is out of the region of interest (ROI) is formed with a code region in which the scatterer is encoded and dispersed in the scanning direction, and the type of the elastic coupler 20 A corresponding ID code can also be assigned.
  • an ultrasonic scatterer can be mixed in the entire elastic coupler 20 to increase the attenuation of the ultrasonic waves.
  • the detection principle of the absolute pressure of the present invention is based on the fact that the thickness of the elastic coupler 20 changes in relation to the applied compression force (pressure), and the correlation depends on the elastic characteristics of the elastic coupler 20. . Therefore, in order to obtain the absolute pressure, it is necessary to measure the elastic characteristics for each type of the elastic coupler 20 in advance and set it in the coupler database in association with the ID code of the elastic coupler 20.
  • the elastic characteristic set in the coupler database is that the deformation (strain) in the thickness direction of the boundary surface or intermediate layer is linear with respect to the absolute pressure, or the elastic modulus (for example, Young's modulus) of the elastic coupler 20 is constant.
  • the elastic modulus E may be set corresponding to the ID code.
  • the relationship curve between deformation (strain) in the thickness direction and absolute pressure when the change in position in the thickness direction of the boundary surface or intermediate layer is non-linear with respect to absolute pressure, when the elastic modulus E is non-linear with respect to absolute pressure, the relationship curve between deformation (strain) in the thickness direction and absolute pressure.
  • a relationship curve between deformation (strain) in the thickness direction and elastic modulus E is set.
  • it can be a relationship curve between the integrated value of deformation or strain in the thickness direction and the elastic modulus, or an elastic modulus correction coefficient for deformation (strain) in the thickness direction.
  • the absolute pressure applied to the living tissue of the subject 100 by the probe 1 can be accurately measured in real time in the normal elasticity image measurement processing.
  • the ID code of the elastic coupler can be read during the elastic image measurement process, even if various different elastic couplers are arbitrarily replaced, the elastic coupler is automatically identified and the elastic characteristics of the elastic coupler are detected. Therefore, the absolute pressure can be measured with high accuracy, and it is possible to improve the usability without the need for the operator.
  • an ultrasonic diagnostic apparatus that generates and displays an elastic image has been described as an example.
  • the present invention is not limited to this, and diagnosis based on a tomographic image (B-mode image) of a living tissue of a subject is performed.
  • the present invention can be applied to an ultrasonic diagnostic apparatus for performing tests such as blood flow diagnosis by Doppler measurement or color flow mode (CFM). Thereby, the evaluation of the compression state suitable for various inspections can be realized with high accuracy.

Landscapes

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

Abstract

　弾性カプラと被検体の体表との接触性が悪くても絶対圧力の検出を精度よく行えるようにする。弾性カプラ20を柔軟性を有する弾性材料により超音波反射特性が少なくとも異なる2つの層を有して形成して超音波送受信面に装着し、圧力演算部30は、RF信号フレームデータ選択部10から出力されるRF信号フレームデータに基づいて、それらの層の境界面22の位置を検出し、境界面の検出位置と予め求められた境界面の初期位置とに基づいて境界面の位置変化を求め、位置変化と予め設定された弾性カプラの弾性特性に基づいて被検体に加えられる絶対圧力を求める。この際に、例えば弾性カプラの境界面の初期位置などを弾性カプラの種類に応じて異ならせることにより弾性カプラにID符号を付与しておき、圧力演算部30はID符号を識別し、カプラデータベースを参照して弾性カプラの種類を識別してID符号に対応した弾性特性を読み出す。

Description

超音波診断装置
　本発明は、超音波診断装置に係り、特に、被検体の生体組織の断層画像(Bモード像)又は弾性情報による診断、ドプラ計測又はカラーフローモード(CFM)による血流診断等の検査時に、超音波探触子によって被検体に加えられる絶対圧力を検出可能にする技術に関する。
　一般に、超音波診断装置は、超音波探触子(以下、単に探触子という。)から超音波を被検体である生体内に送波し、生体内から反射される超音波の反射エコー信号を探触子により受波し、受波した反射エコー信号(RF信号)に基づいて、生体内の組織や機能等の検査に適した画像等を生成して診断を行なうものである。
　このような超音波診断装置において、超音波は生体組織内部の伝播過程において減衰することから、Bモード像による検査は、比較的強い力で探触子を被検体に押し当てて生体組織を圧迫して変形させ、深部の組織を探触子側に近づけて撮像することにより高画質な画像を得ることが望ましい。一方、ドプラ計測やCFMなどの血流検査では、あまり強く探触子を押し当てて生体組織を圧迫すると、血管の断面が潰れて正しい血流情報を取得できないため、Bモード診断時よりも軽微な圧迫状態で検査を行うことが望ましい。また、生体組織の硬さ又は軟らかさの弾性情報を画像化するエラストグラフィによる検査でも、生体組織は圧迫の強さに応じて組織の硬さが変化する非線形性を有するため、一定の絶対圧力による圧迫状態下で得られた弾性画像に基づいて診断を行うことが重要とされる。
　したがって、各検査法あるいは検査法を切り替えたときなどに、不適切な圧迫状態のまま検査を進めてしまうと、適切な診断を速やかに行なうことを妨げるおそれがあるから、生体組織に加わる絶対圧力をリアルタイムで計測して表示することが望まれる。
　特許文献1には、被検体の生体組織に加わる実際の圧力、つまり生体組織に加わる絶対的な圧力(以下、単に絶対圧力という。)を計測するため、取得時間が異なる一対のRF信号フレームデータに基づいて周知の変位・歪み演算を用いて弾性カプラの歪みを求め、予め設定されている弾性カプラの弾性率を用いて、求めた歪みを絶対圧力に換算することが提案されている。
　したがって、特許文献1に記載の圧力演算方法によれば、例えば、被検体に加えられた絶対圧力を計測して、生体組織の硬さ又は軟らかさの弾性情報を取得することができる。
特開2005-66041号公報
　しかしながら、特許文献1の弾性カプラにより絶対圧力を検出する方法によれば、基本的に、絶対圧力を高感度かつ高精度で検出することができるが、次に述べるような改善すべき余地がある。
　例えば、弾性カプラの被検体側の接触面と被検体の体表との接触性が悪く、微小な凹凸や空気層などが介在する接触ムラがあると、被検体の体表との境界を適正に検出することができず、弾性カプラの厚みの検出を適正にできないことがある。また、接触ムラがあると、非接触部の周囲の接触部の変形が大きくなることから、その境界部分の変形や歪みを計測して絶対圧力に換算すると、適正な絶対圧力の分布を計測できないことがある。
　本発明が解決しようとする課題は、絶対圧力の検出を精度よく行えるようにすることにある。　
　また、上記の課題に加えて、検査法等に応じて適切な弾性カプラを用いて絶対圧力を検出する操作を簡単化して、使い勝手を向上することにある。
　本発明の第1の態様は、被検体に当接させて超音波を送受する超音波探触子と、該超音波探触子を駆動する送信部と、前記超音波探触子により受信された反射エコー信号であるRF信号を受信処理する受信部と、該受信部から出力されるRF信号に基づいて超音波像を生成する画像生成部とを備えてなる超音波診断装置において、
　前記超音波探触子の超音波送受信面に装着された弾性カプラの変形に基づいて前記被検体に加えられた圧力を求める圧力演算部を備え、前記弾性カプラは、少なくとも2つの層を有して形成され、　前記圧力演算部は、前記2つの層間の境界面の位置を検出し、該境界面の検出位置と予め求められた境界面の初期位置とに基づいて前記境界面の位置変化を求め、該位置変化と前記弾性カプラの弾性特性に基づいて絶対圧力を求める。
　すなわち、前記弾性カプラは、柔軟性を有する弾性材料により超音波反射特性が異なる少なくとも2つの層を有して形成され、前記超音波反射特性が異なる層の境界面が前記超音波送受信面に装着される装着面と前記被検体の体表に当接される接触面との間に位置させて設けられてなり、前記圧力演算部は、前記受信部から出力される前記RF信号に基づいて、前記弾性カプラの厚み方向の前記境界面の位置を検出し、該境界面の検出位置と予め求められた境界面の初期位置とに基づいて前記境界面の位置変化を求め、該位置変化と予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求めることを特徴とする。
　このように、弾性カプラの内部に形成された境界面の位置を検出し、その位置変化に基づいて絶対圧力を求めるようにしているから、弾性カプラと被検体の体表との境界を検出する場合に比べて、弾性カプラと被検体の体表との接触性にかかわらず、内部の境界面を安定して検出できるから、絶対圧力の検出を精度よく行える。なお、弾性カプラに圧迫力が加えられていない初期状態における境界面の初期位置は、圧迫力が加えられている加圧状態における境界面の位置検出と同様に、検出受信部から出力されるRF信号に基づいて検出することができる。
　また、絶対圧力の検出を精度よく行えるようにするため、本発明の第2の態様の特徴構成は、弾性カプラは、第1の態様と同一であるが、圧力演算部の構成が異なる。つまり、前記圧力演算部は、前記受信部から出力される取得時刻が異なる一対の前記RF信号フレームデータに基づいて、前記弾性カプラの厚み方向の前記境界面の位置の変位を検出し、該変位に基づいて前記境界面の厚み方向の歪みを求め、該厚み方向の歪みと予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求めることを特徴とする。
　これによれば、弾性カプラと被検体の体表との接触性にかかわらず、内部の境界面を安定して検出できるから、その境界面の歪みを精度よく検出でき、歪みと弾性特性との関係に基づいて、絶対圧力を精度よく検出できる。
　さらに、第2の態様において、前記圧力演算部は、前記弾性カプラが無加圧の初期状態から時間経過に合わせて前記境界面の厚み方向の歪みを積算し、該歪みの積算値と予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求めるようにすることができる。これによれば、境界面の厚み方向の歪みを積算していることから、一層精度よく境界面の歪みを検出でき、歪み積算値と弾性特性との関係に基づいて、絶対圧力を一層精度よく検出できる。
　第1又は第2の態様において、前記弾性カプラは、前記2層の境界面に薄い中間層を挟んで形成され、該中間層の超音波反射特性は他の2層と異ならせて形成された構成とすることができる。この場合において、中間層は、前記超音波探触子の超音波ビームのスキャン方向に直交する方向に延在させて、かつ間隔をあけて設けられた線状の複数の超音波反射体から構成することができる。
　弾性カプラに中間層を設けた場合、第1の態様の圧力演算部は、前記中間層の厚み方向の位置を検出し、該中間層の検出位置と予め求められた中間層の初期位置とに基づいて前記中間層の位置変化を求め、該位置変化と予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求める構成とすることができる。
　また、弾性カプラに中間層を設けた場合、第2の態様の圧力演算部は、前記受信部から出力される取得時刻が異なる一対の前記RF信号フレームデータに基づいて、前記弾性カプラの前記中間層の厚み方向の位置の変位を検出し、該変位に基づいて前記中間層の厚み方向の歪みを求め、該厚み方向の歪みと予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求める構成とすることができる。さらに、第2の態様の圧力演算部は、前記弾性カプラが無加圧の初期状態から時間経過に合わせて前記中間層の厚み方向の歪みを積算し、該歪みの積算値と予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求める構成とすることができる。
　さらに、前記弾性カプラは、2層構造を有して形成され、前記被検体の体表に当接される接触面側の層が薄く、かつ前記超音波反射特性が強く形成されてなる構成とすることができる。また、前記弾性カプラは、前記弾性材料の中に超音波散乱体を混入して超音波減衰特性を大きく形成してなる構成とすることができる。
　一方、検査法等に応じて適切な弾性カプラを用いて絶対圧力を検出する操作を簡単化して、使い勝手を向上するという課題を解決するため、第1又は第2の態様の前記弾性カプラを、前記2層の境界面の厚み方向の位置を当該弾性カプラの種類に応じて異ならせて、前記RF信号により当該弾性カプラの種類を識別可能な識別符合が形成された構成とする。この場合の圧力演算部は、前記RF信号又は前記RF信号フレームデータに基づいて前記識別符合を検出して前記弾性カプラの種類を識別し、前記弾性カプラの種類に対応して設定されている前記弾性特性に基づいて前記絶対圧力を求める構成とすることができる。
　これにより、弾性特性が異なる複数の弾性カプラを用意しておき、検査法あるいは関心部位の深度などに応じた適切な絶対圧力の計測に適した弾性カプラに交換した場合、超音波探触子に装着された状態で弾性カプラの種類を圧力演算部により自動的に識別できる。したがって、予め弾性カプラの種類に応じた弾性特性を圧力演算部に記憶させておけば、弾性カプラの種類の弾性特性に合わせて絶対圧力を演算することができる。
　また、前記識別符号は、前記2層の境界面の厚み方向の位置を異ならせて形成された構成とすることができる。この場合、前記圧力演算部は、前記RF信号又は前記RF信号フレームデータに基づいてカプラエコー領域内のRF信号の深度分布パターンにより、弾性カプラの種類を識別する構成とすることができる。
　また、前記識別符号は、前記弾性カプラのスキャン方向の両端の領域部に、スキャン方向と厚み方向の少なくとも一方に超音波散乱体を符号化して分散させて形成された構成とすることができる。この場合、前記圧力演算部は、前記RF信号又は前記RF信号フレームデータに基づいて前記弾性カプラの両端の領域部のRF信号のパターンにより、弾性カプラの種類を識別する構成とすることができる。
　ここで、カプラデータベースに設定する弾性特性は、弾性率、厚み方向の変形と弾性率の関係曲線、厚み方向の歪みと弾性率の関係曲線、厚み方向の変形又は歪みの積算値と弾性率の関係曲線、厚み方向の変形又は歪みに対する弾性率補正係数の少なくとも一つとすることができる。
　本発明によれば、絶対的圧力の検出を精度よく行うことができる。また、操作が簡単で、使い勝手が向上する。
本発明の実施形態の超音波診断装置のブロック構成図。 弾性カプラの一実施例の構成図。 弾性カプラの装着具の一実施例の構成図。 実施形態の圧力演算部の詳細構成を示すブロック構成図。 実施形態の特徴部の処理手順を説明するフローチャート。 図2の弾性カプラの超音波像を示す図。 圧迫状態画像の一例を示す図。 圧迫状態画像の他の例を示す図。 弾性カプラの他の実施例の構成図及び弾性カプラの超音波像を示す図。 弾性カプラのさらに他の実施例の構成図及び弾性カプラの超音波像を示す図。
　1　探触子、2　送信回路、3　受信回路、8　切替加算器、9　画像表示器、10　RF信号フレームデータ選択部、11　変位・歪み演算部、12　弾性率演算部、13　弾性データ処理部、20　弾性カプラ、30　圧力演算部、31　カプラID識別部、32　初期化処理部、33　カプラ歪み演算部、34　圧力変換部、35　圧迫状態画像構築部、36　カプラデータベース
　以下、本発明の超音波診断装置の実施形態について、図面を参照して説明する。　
　図1に、本発明を弾性情報による検査を行う超音波診断装置に適用した一実施形態の機能ブロック構成図を示す。
　図1に示すように、本実施形態の超音波診断装置は、超音波を利用して被検体の診断部位について断層画像を得るとともに生体組織の硬さ又は軟らかさを表す弾性画像を得るものである。同図に示すように、超音波探触子(以下、単に、探触子という。)1は送信回路2及び受信回路3と電気的に接続されている。探触子1は、周知のように、超音波の発生源となり、かつ反射エコーを受波する振動子を短冊状に多数配列して形成されている。各振動子は、一般に、入力されるパルス波又は連続波の超音波信号を超音波に変換して発射する機能と、被検体100の内部から発射する反射エコーを電気信号である反射エコー信号に変換して出力する機能を有している。
　送信回路2は、送信回路2と受信回路3は、超音波送受信制御回路4によって協調制御されるようになっている。超音波送受信制御回路4は、超音波を送信及び受信するタイミングを制御するものであり、目的とする超音波検査に適した超音波送波ビームを探触子1から被検体100に打ち込むように送信回路2を制御するとともに、探触子1から所望の超音波受波ビームを受波するように受信回路3を制御するようになっている。本実施形態の超音波送受信制御回路4は、Bモード断層画像の計測に適した超音波送波ビーム及び超音波受波ビームを形成するように、超音波送波ビームを断層面に沿って走査するように送信回路2と受信回路3を制御する。
　受信回路3は、探触子1で受信した反射エコー信号を所定のゲインで増幅する。増幅された各振動子の数に対応した数の反射エコー信号が加算回路5に入力される。加算回路5は、受信回路3で増幅された複数の反射エコー信号の位相を合わせて加算して、断層面に対応したRF信号フレームデータを生成する。
　信号処理部6は、加算回路5から出力されるRF信号フレームデータを入力して、ゲイン補正、ログ補正、検波、輪郭強調、フィルタ処理等の各種信号処理を行なって画像データを生成する。白黒スキャンコンバータ7は、信号処理部6から出力される画像データを音波周期で取得し、この超音波像を表示するためテレビジョン方式の周期で読み出すための断層走査手段及びシステムの制御を行うための手段、例えば、信号処理部6からの画像データをディジタル信号に変換するA/D変換器と、このA/D変換器でディジタル化された画像データを時系列に記憶する複数枚のフレームメモリと、これらの動作を制御するコントローラなどを含んで構成される。白黒スキャンコンバータ7によって生成された時系列の例えばBモード像の画像データは、切替加算器8を介して画像表示器9に出力される。画像表示器9は、白黒スキャンコンバータ7から出力される画像データをアナログ信号に変換するD/A変換器と、このD/A変換器から出力されるアナログビデオ信号を画像として表示するカラーテレビモニタとからなる。
　次に、本実施形態の弾性画像を生成して画像表示器9に表示させる処理について説明する。一般に、弾性画像を生成するために、探触子1により被検体100を圧迫し、診断部位の生体組織に圧力(応力)を加え、その応力によって生体組織が変形することを利用している。つまり、生体組織は同一の応力が作用しても、生体組織の硬さ又は軟らかさである弾性に応じて変形の度合が異なる。そこで、探触子1で超音波送受信を行なうとともに、探触子1により被検体100を異なる圧力で圧迫したときの一対のBモード像に対応するRFフレームデータを用いて、組織の変形の度合を求めて、正常部位や病変部位を識別できる弾性画像を生成するようにしている。
　本実施形態では、加算回路5で生成されたRF信号フレームデータは、RF信号フレームデータ選択部10に入力されている。RF信号フレームデータ選択部10は、加算回路5から順次出力されるRF信号フレームデータをフレームメモリ内に順次蓄積する。そして、図示していない制御部からの指令に従って、例えば最新のRF信号フレームデータと、時間的に過去のRF信号フレームデータの中から1つのRF信号フレームデータを選択して、変位・歪み演算部11に取得時刻が異なる一対のRF信号フレームデータを出力するようになっている。なお、RF信号は符号復調したI，Q信号の形式になった信号であってもよい。
　変位・歪み演算部11は、RF信号フレームデータ選択部10から出力される一対のRF信号フレームデータについて1次元又は2次元の相関処理を実行し、RF信号フレームデータの超音波ビームのスキャン方向及び深度方向にそれぞれ設定された複数の計測点i，jについて、変位(又は変位ベクトル)を演算する。演算された複数の計測点の変位データは、変位フレームデータとして生成される。なお、変位演算の方法としては、相関処理の他に、例えば、特許文献1に記載されたように、周知のブロック・マッチング法やグラジェント法を用いることができる。ブロック・マッチング法は、注目している計測点の画素を中心として複数画素からなる注目ブロックを設定し、その注目ブロックの画像情報が近似しているブロックを変位前のフレームに対して移動しながら探索し、最も近似しているブロックの位置から現在位置に変位したものとする方法である。
　次いで、変位・歪み演算部11は、変位フレームデータを用い、周知のように、各計測点の変位を空間微分して各計測点の歪みを演算する。演算により得られた各計測点の歪みにより歪みフレームデータを生成して弾性率演算部12に出力する。
　弾性率演算部12は、入力される歪みフレームデータの各計測点の歪みに基づいて、後述する圧力演算部30から与えられる各計測点i，jに作用している絶対圧力(応力)に基づいて、周知の方法により、各計測点の弾性率(例えば、ヤング率)を演算する。演算により得られた各計測点の弾性率により弾性率フレームデータを生成して、弾性データ処理部13に出力する。弾性データ処理部13は、弾性率演算部12で生成された歪み又は弾性率の弾性フレームデータに対し、フレーム内におけるスムージング処理、コントラスト最適化処理、あるいはフレーム間における時間軸方向のスムージング処理などの様々な画像処理を行ない、カラースキャンコンバータ14に出力する。
　カラースキャンコンバータ14は、弾性データ処理部13から出力された弾性フレームデータに基づいて、各計測点に対応する画素に赤、緑、青などの色相情報を付与して弾性画像データを生成する。例えば、弾性データ処理部13から出力される歪みフレームデータにおいて、歪みが大きく計測された画素については、弾性画像データ内でその画素を赤色コードに変換し、逆に歪みが小さく計測された画素については、弾性画像データ内でその画素を青色コードに変換する。弾性率フレームデータの場合も同様に各画素に色相情報を付与して弾性画像データを生成する。なお、カラースキャンコンバータ14に代えて、白黒スキャンコンバータを使用して、歪みの大きさ等に対応して輝度階層を付与するようにすることもできる。
　カラースキャンコンバータ14で生成された弾性画像データは、切替加算部8に出力される。切替加算部8は、白黒スキャンコンバータ7から出力される白黒の断層画像データと、カラースキャンコンバータ14から出力されるカラーの弾性画像データとを入力し、入力指令に応じて、両画像を加算又は切り替えて画像表示器9に出力する。具体的には、白黒の断層画像データだけ、カラーの弾性画像データだけ、又は両画像を並べて、画像表示器9に切替表示させることができる。また、両画像データを加算合成して、つまり半透明の重ね合わせ画像を生成して、画像表示器9に出力表示させるなど、入力指令に応じて出力画像を切り替えることができるようになっている。
　次に、本実施形態の特徴部である弾性カプラ20と、圧力演算部30の詳細構成を説明する。　
　図2に、弾性カプラ20の一実施例を示す。弾性カプラ20は、図2(A)に示す斜視図のように、柔軟性を有する弾性材料であるゲル状物質により板状に形成され、一方の面の周囲に鍔部21が形成されている。弾性カプラ20を形成するゲル状物質の素材は、特許文献1に記載されているように、音響結合材料や音響レンズ素材などのように、超音波減衰が小さく、かつ、音速及び音響インピーダンスが生体に近く、生体との音響結合特性に優れたものが好ましい。また、同時に、保形性、柔軟性、適度な弾性、形状復元性に優れた素材を用いることが好ましい。特に、本実施例の弾性カプラ20は、図示のように、超音波反射特性が異なる2つの層20A、Bを有して形成され、それらの層20A、Bの境界面22が探触子1の超音波送受信面に装着される装着面と被検体100の体表に当接される接触面との間に位置させて設けられている。
　また、弾性カプラ20は、探触子1の超音波送受信面に装着した場合、超音波送受信面との接触面及び被検体100の体表のように、凹凸のある部分に接触させても隙間を生じない材料が好ましい。すなわち、超音波送受信面及び体表との間に隙間が生じると、探触子1から照射された超音波がその隙間の空気の境界で反射されて、超音波像のノイズとなる。例えば、ゼラチン、寒天、オイルゲル、アクリルアミド、ポリビニルアルコール等の水性高分子と水とからなる水性ゲル(ハイドロゲル)、ポリウレタン、ゴムと油状成分とを含む組成物を架橋成形したオイルゲル、生ゴムに可塑剤として低分子量のゴムを配合したものを含む組成物を成形架橋したゴムを用いることができる。
　このように形成された弾性カプラ20を、装着具23を用いて探触子1の超音波送受信面24に装着して用いる。すなわち、装着具23は、図3(A)に示すように、弾性カプラ20の鍔部21に合わせて形成された溝25を有して樹脂等により枠状に形成されている。また、装着具23の枠体の内側に超音波送受信面24の外周部に係合する係止爪26が形成されている。また、枠体の対向する長辺に一対の把持片27が形成され、頂面には弾性カプラ20が挿入される開口28が形成されている。このように形成された装着具23に、図3(B)に示すように、底部側から弾性カプラ20を開口28に挿入し、弾性カプラ20の鍔部21を枠体の溝25に押し込んで組み付け、図2(B)に示すように、探触子1の超音波送受信面24に装着して用いる。このとき、弾性カプラ20の下面と超音波送受信面24との接触面に隙間ができないように、ゼリーなどを介在させて探触子1の超音波送受信面24に当て、係止爪26を超音波送受信面24の外周部に引っ掛けて、しっかりと固定する。
　このようにして弾性カプラ20が装着された探触子1の把持部1Aを把持して、弾性カプラ20の露出面を被検体100の体表に押し当てて、超音波検査を行なう。なお、図1に示すように、探触子1の把持部1Aを被検体圧迫機構18に把持させて、被検体圧迫機構18により被検体100の体表に圧迫力を加えるようにすることができる。
　圧力演算部30は、図4に示すブロックを備えて構成されている。すなわち、カプラID識別部31と、初期化処理部32と、カプラ歪み演算部33と、圧力変換部34と、圧迫状態画像構築部35と、カプラデータベース36を備えて構成される。
　本実施形態に適用する弾性カプラ20は、図2(A)に示したものと同一の形状であるが、弾性カプラ20の種類を自動識別できるように、用途に応じた弾性カプラ20の種類ごとに、境界面22の深さ、あるいは層20A，Bの反射特性を異ならせることにより、同一種類の弾性カプラ20に識別符合(ID)を付与したものを用いる。
　カプラID識別部31は、RF信号フレームデータをRF信号フレームデータ選択部10から取り込み、境界面22の有無を検出して、弾性カプラ20が探触子に装着されていることを検知する。さらに、カプラエコー領域内のRF信号等の深度分布を判別して境界面22の深さを検出し、カプラデータベース36を参照して、弾性カプラ20のID符号を自動で識別する。識別したID符号は、圧力変換部34に出力される。なお、カプラデータベース36には、図示していない入力手段によって、境界面22の深さ及び深さ方向の位置パターンに対応させてID符号が設定され、さらにID符号に対応させて弾性特性が設定されている。
　初期化処理部32は、RF信号フレームデータのRF信号の強度変化に基づいて、弾性カプラ20が無加圧の初期状態にあることを検出し、弾性カプラ20が初期状態におけるRF信号に基づいて弾性カプラ20の初期厚みを求めるようになっている。
　カプラ歪み演算部33は、RF信号フレームデータに基づいて弾性カプラ20の境界面22の厚み方向の位置を検出するとともに、境界面22の厚み方向の位置の変化を検出して、境界面22の変位及び歪みを求めるようになっている。求めた変位又は歪みはID符号とともに圧力変換部33に出力される。
　圧力変換部34は、入力されるID符号に対応した弾性特性をカプラデータベース36から読み出し、読み出した弾性特性に基づいて、カプラ歪み演算部33から入力される変位又は歪みを絶対圧力に換算するようになっている。
　圧迫状態画像構築部35は、圧力変換部34から出力される絶対圧力を画像表示器9に表示させるために、圧迫状態画像を構築するようになっている。
　以下、図5に示したフローチャートに沿って、本実施形態の圧力演算部50における圧力演算処理を詳細に説明する。
　［S1：超音波診断装置の起動］
　超音波診断装置を手動により起動する。検査モードは、例えば、Bモード断層画像による検査とする。
　［S2：カプラ装着及び検知］
　手操作で図2(B)のように探触子1の超音波送受信面24に弾性カプラ20を装着する。弾性カプラ20が装着されたか否かの検出は、カプラID識別部31においてRF信号フレームデータを取り込み、RF信号の強度分布等により境界面22の有無を検出できたか否かにより自動で認識できる。また、図示していない入力手段を介して、操作者によりカプラID識別部31に弾性カプラ20を装着したことを入力するようにしてもよい。
　［S3：カプラIDの識別］
　カプラID識別部31は、RF信号フレームデータをRF信号フレームデータ選択部10から取り込み、境界面22の存在を検出した後、カプラエコー領域内のRF信号等の深度分布を判別して境界面22の深さ位置を検出し、カプラデータベース36を参照して、弾性カプラ20のID符号を自動で識別する。識別したID符号は、圧力変換部34に出力する。図6に、弾性カプラ20を装着して得られるBモード像101の一例を示す。図のように、画像の上部に弾性カプラ20の2層20A，Bと境界面22のエコー画像が現れる。したがって、カプラエコー領域内のRF信号が急激に変化する位置を検出することにより、境界面22に位置を検出して、ID符号を自動で識別できる。
　［S4：圧力演算部の起動］
　圧力演算部30は、カプラID識別部31において弾性カプラ20のID符号を検知して起動される。
　［S5：初期状態の生成］
　弾性カプラ20の境界面22の深さ方向の初期位置D(0)は予め計測できるが、弾性カプラ20の素材であるゲル状物質の液体成分が蒸発等することがあり、経時変化により境界面22の初期位置D(0)が小さくなることが考えられる。そこで、精度よく絶対圧力を検出するためには、超音波検査の都度、初期状態の境界面の初期位置D(0)を自動計測することが望ましい。
　本実施形態では、弾性カプラ20の露出面を空中に保持した無加圧の初期状態で境界面の初期位置D(0)を自動計測するために、初期化処理部32において、操作者が探触子1を把持して空中に保持した初期状態を自動認識可能にしている。この無加圧の初期状態は、RF信号に含まれる多重エコー信号により検知できる。つまり、音響レンズの露出面と空気との境界、又は弾性カプラ20の露出面と空気との境界において音響インピーダンスが急激に変化することに起因して多重エコーが発生する。例えば、探触子1に弾性カプラ20を装着されていない状態で空中に保持してBモード検査用の超音波ビームを送波すると、超音波送受信面に近い深度領域に音響レンズ(1mm程度厚)に由来する周期的な多重エコーを含むRF信号が受信される。他方、弾性カプラ20(例えば、5～10mm程度厚)が装着されると、音響レンズ由来の多重エコーが消失して、弾性カプラ20に由来する比較的長い周期の多重エコーが受信される。そこで、多重エコーの有無及び周期を検出することにより、弾性カプラ20が無加圧の初期状態の保持されていることを自動で認識できる。
　なお、初期状態を自動認識することに代えて、操作者が探触子1を把持して空中に保持した状態で、図示していない入力手段から、手動操作によりカプラ厚み演算部32に境界面22の初期位置D(0)の自動計測を指令することができる。
　［S6：初期化処理の実行］
　初期化処理部32は、弾性カプラ20の初期状態を認識したとき、又は指令が入力されたとき、境界面の初期位置D(0)の自動計測を実行する。
　初期化処理部32は、弾性カプラ20の探触子1の振動子の配列方向(長軸方向)である超音波ビームのスキャン方向の座標i(i＝0～n)、及び深度方向の座標j(j＝0～m)における計測点ijついて、RF信号(又は輝度)の強度Qijが予め設定された閾値Qthresより大きくなる又は小さくなる座標jの平均値Jを求め、J＝境界面22の初期位置D(0)とする。なお、座標Jは、超音波送受信面から超音波が送波されから座標Jに到達して戻って来るまでの時間の概念であるから、弾性カプラ20内の音速をCとすると、弾性カプラ20の境界面の初期位置D(0)は、D(0)＝J×C/2で算出できる。
　また、初期化処理部32は、RFフレーム信号選択部10から出力されるRF信号に代えて、信号処理部6から出力されるBモード像の輝度データ又は信号強度を取り込み、同様の処理にて弾性カプラ20の境界面の初期位置D(0)を検出することができる。
　言い換えれば、送信回路2から探触子1に超音波を送波してから、RF信号の強度が大きく変化する時間ti(0)を求めて、境界面の初期位置D(0)を検出することができる。ti(0)は弾性カプラ20の境界面から反射する超音波の往復の時間であるから、ti(0)に音速Cを乗じて2で割ると、スキャン方向の1次元分布である境界面の初期位置Di(0)を求めることができる。また、RF信号の強度が大きく変化する多重エコーの周期(閾値以上の周期)Tを求め、Tの1/2に音速Cを乗じてもDi(0)を求めることができる。
　なお、弾性カプラ20の境界面22の初期位置Di(0)は、スキャン方向に直交する方向の分布を加味した2次元分布を求めることができる。
　［S7：カプラ変形(歪み)計測］
　S7における処理は、カプラ歪み演算部33における処理である。まず、カプラ歪み演算部33は、前述した弾性カプラ20に起因する多重エコーが消失していることにより、弾性カプラ20が加圧状態にあることを検出することができる。
　カプラ変形(歪み)計測の方式1
　カプラ歪み演算部33は、弾性カプラ20の加圧状態における任意の時刻tのRF信号に基づいて弾性カプラ20の境界面22の深さ位置を検出して、弾性カプラ20の加圧状態におけるスキャン方向の境界面位置分布Di(t)を求める。つまり、操作者が弾性カプラ20を介して被検体100に探触子1を押し付けて圧迫を加え、その加圧状態で超音波を送信してから、RF信号Qi(t)の強度が大きく変化するまでの時間ti(t)の1/2と音速Cに基づいて、超音波ビームに直交するスキャン方向の境界面位置分布Di(t)を求める。
　次いで、次式(1)により境界面位置変化分布ΔDi(t)を求め、次式(2)により弾性カプラ20の境界面22のスキャン方向の全歪み量分布Si(t)を求める。
　　ΔDi(t)＝Di(0)-Di(t)　　　　　　(1)
　　Si(t)＝ΔDi(t)/Di(0)　　　　　　(2)
　カプラ変形(歪み)計測の方式2
　上記方式1に代えて、カプラ歪み演算部33は、弾性カプラ20の歪みを計測し、これに基づいて絶対圧力Pi(t)を求めることができる。すなわち、カプラ歪み演算部33は、RF信号フレームデータ選択部10から出力される初期状態におけるカプラエコー領域のRF信号フレームデータQij(0)と、加圧状態における任意の時刻tにおいてRF信号フレームデータ選択部10から出力されるRF信号フレームデータQij(t)とに基づいて、各計測時点における弾性カプラ内部の計測点の変位を求めて全歪み量Sij(t)を求めることができる。ここで、前述したように、iは弾性カプラ20のスキャン方向の座標、jは弾性カプラ20の厚み方向(深度方向)の座標である。
　すなわち、カプラ歪み演算部33は、加圧状態においてリアルタイムで変化するカプラエコー領域のRF信号Qij(t)を取り込み、初期状態におけるQij(0)と、任意の時刻tにおけるQij(t)に基づいて周知の変位演算方法によって、各計測点i，jの変位を求めて変位フレームデータを生成する。そして、変位フレームデータを空間微分して各計測点i，jの全歪み量Sij(t)からなる歪みフレームデータを求める。
　次いで、歪みフレームデータの全歪み量Sij(t)に基づいて、同一のスキャン方向の座標iごとに、カプラエコー領域の深度jの全範囲にわたって全歪み量Sij(t)の平均値Si＊(t)を求め、その平均値Si＊(t)を超音波スキャン方向の座標位置iにおける全歪み量Sij(t)として評価する。そして、評価した全歪み量Sij(t)を用いて、圧力変換部38において絶対圧力Pi(t)を求める。
　本方式によれば、図1の弾性演算部12における生体組織の歪み計算と同時に、絶対圧力Pi(t)を求める処理を実行できる。
　本方式の場合、弾性カプラ20に超音波の散乱体を分散して混入することにより、RF信号の強度を強くできるので、厚みあるいは歪みの演算精度が向上する。
　カプラ変形(歪み)計測の方式3
　カプラ圧迫評価の方式3は、初期状態から加圧状態における各計測時点について継続して、RF信号フレームデータ選択部10から取得時間が異なる一対のRF信号フレームデータを取得する。そして、一対のRF信号フレームデータを取得するたびに、弾性カプラ20の境界面22までの全領域について、各計測点i，jの歪み変化ΔSij(t-k)、・・・、ΔSij(t)を求める。そして、時間的に隣り合う、例えば時刻(t-1)と時刻(t)の歪み変化ΔSi、j(t)を求める。さらに、時間的に連続する一対のRF信号フレームデータについて、その歪み変化ΔSi、j(t-k)、・・・、ΔSij(t)を順次積算して、現時刻の歪み変化積算値ΣΔSij(t)を求める。次いで、計測点i、jのΣΔSij(t)について、カプラエコー領域の座標jの方向について平均し、計測点iについての歪み変化積算値ΣΔSi＊(t)を求める。
　［S8：絶対圧力の換算処理］
　圧力変換部34は、カプラID識別部31で識別された弾性カプラ20のIDに対応する弾性特性(例えば、ヤング率E)をカプラデータベース36から読み出し、カプラ歪み演算部33で方式1により求めた歪みSi(t)を、次式(3)により、被検体100の生体組織に加えられている絶対圧力分布Pi(t)に換算する。
　　　Pi(t)＝Si(t)×E　　　　　　　　　　(3)
　換算により得られた現在の絶対圧力分布Pi(t)は、図1の弾性率演算部12に出力される。これにより、弾性率演算部12は、前述したように、周知の演算処理により、生体組織の各計測点i，jについて得られた歪みεij(t)に基づいて、各計測点i，jの弾性率(例えば、ヤング率)Eij(t)を、次式(4)により求めて、弾性データ処理部13に出力する。
　　　Eij(t)＝Pi(t)/εij(t)　　　　(4)
　また、カプラ歪み演算部33で方式2の平均値Si＊(t)に基づいて求めた全歪み量Sij(t)、又は方式3により求めた歪み変化積算値ΣΔSi＊(t)を、次式(5)、(6)により、絶対圧力分布Pi(t)に換算する。なお、圧力変換部34はカプラ歪み演算部33と一体に形成してもよい。
　　　Pi(t)＝Si＊(t)×E　　　　　　　　　(5)
　　　Pi(t)＝ΣΔSi＊(t)×E　　　　　　　　(6)
　［S9：圧力分布の表示処理］
　圧迫状態画像構築部35は、S8の絶対圧力の換算処理で求めた絶対圧力Pi(t)を、画像化して画像表示器9に表示させることにより、操作者は超音波検査を実行しながら同一の画面で検査項目に応じた適切な圧迫状態か否かを即座に判断できる。
　すなわち、圧迫状態画像構築部38は、例えば、絶対圧力Pi(t)を、数値表示、経時変化のグラフ表示、バーチャート表示などの少なくとも一つの画像を構築して、カラースキャンコンバータ14によってカラー画像データに変換して、画像表示器9に表示される超音波像に並べて又は一部に重ねて表示させる。
　図7に、圧迫状態画像を弾性画像に重ねて並びに並べて表示した一例を示す。図において、画面中央に弾性画像110が表示され、矩形の関心領域(ROI)の近くに弾性率E(kPa)と絶対圧力σ(kPa)が数値で表示窓111に表示されている。また、現時刻の絶対圧力σが絶対圧力Piの平均値などのバーチャート112が表示されている。さらに、画面の上部にスキャン方向の絶対圧力分布Piのグラフ113が、下部に絶対圧力Piの平均値の経時変化のグラフ114が表示されている。さらに、画面の右端に、弾性率Eのカラーバー115が表示されている。
　したがって、操作者は、図7の画像を観察することにより、ROIの生体組織について適切な絶対圧力の下における弾性を評価して、診断を行うことができる。特に、絶対圧力Piの平均値などのバーチャート112を見て、圧迫力が適正か否かを即座に判断することができる。
　図8に、圧迫状態画像を弾性画像に重ねて又は並べて表示した他の例を示す。図において、画面中央に弾性画像110が表示され、矩形の関心領域(ROI)の近くに弾性率E(kPa)と絶対圧力σ(kPa)がバーチャート116，117で表示されている。また、初期状態からの弾性カプラの厚み変化ΔDi(t)の平均値などを表すバーチャート118が表示されている。さらに、画面の上部に弾性カプラと被検体との境界が複数の点119により表示されている。
　また、圧力演算部30の演算結果である絶対圧力を検査法に対応させて設定され基準圧力範囲を比較して、基準圧力範囲を超えた場合は、例えば、ROIを示す境界線の表示色を変化させたり、又は点滅させることにより、圧迫状態の適否を容易に識別することができる。
　以上説明したように、本実施形態によれば、通常の弾性画像計測の検査の処理において、探触子1により被検体100の生体組織に加えられる絶対圧力をリアルタイムで精度よく計測できる。
　また、弾性カプラのID符号を弾性画像計測の処理の中で読み取ることができるから、種々の異なる弾性カプラを任意に交換しても、自動的に弾性カプラを識別して、弾性カプラの弾性特性に応じて絶対圧力を精度よく計測でき、操作者の手間が省けて使い勝手をよくすることができる。
　また、本実施形態では、弾性画像を生成表示する超音波診断装置を例に説明したが、本発明はこれに限られるものではなく、被検体の生体組織の断層画像(Bモード像)による診断、ドプラ計測又はカラーフローモード(CFM)による血流診断等の検査を行う超音波診断装置に適用することができる。これにより、種々の検査に適した圧迫状態の評価を精度よく実現できる。
　ここで、弾性カプラの他の実施例について説明する。　
　図9(A)に、弾性カプラ20の実施例2の斜視図を示す。装着具は、図3と同一のものを用いることができる。図9(A)に示すように、本実施例2が図2(A)の実施例1の弾性カプラ20と異なる点は、超音波反射特性が同一の2層20A，Bの境界面に薄い中間層20Cを挟んで形成したことにある。中間層20Cの超音波反射特性は、他の2層20A，Bと異ならせている。この場合、中間層20Cの反射強度が、他の2層20A，Bに比べて高くても、低くてもよい。この場合のBモード像101の一例を図9(B)に示す。図示のように、画像の上部に弾性カプラ20の2層20A，Bの境界に中間層20Cのエコー画像が現れる。
　この場合、カプラID識別部31は、実施例1の場合と同様に、RF信号の厚み方向の分布に基づいて、弾性カプラのID符号を識別するようにできる。また、カプラ歪み演算部33は、実施例1と同様に、中間層20Cの厚み方向の位置変化、中間層20Cの厚み方向の歪み、あるいは弾性カプラが無加圧の初期状態から時間経過に合わせて中間層20Cの厚み方向の歪みを積算した歪みの積算値を求めて、被検体に加えられる絶対圧力を求めるようにすることができる。
　また、図10(A)に、弾性カプラ20の実施例3の斜視図を示す。装着具は、図3と同一のものを用いることができる。本実施例3が図9(A)の実施例2の弾性カプラ20と異なる点は、中間層20Cに代えて、探触子1の超音波ビームのスキャン方向に直交する方向に延在させて、かつ間隔をあけて線状の複数の超音波反射体20Dを設けたことにある。
　本実施例3の場合も、実施例2と同様に、超音波反射体20Dの深さ方向の位置に基づいて、弾性カプラのID符号を識別するようにでき、また、超音波反射体20Dの位置変化に基づいて、絶対圧力を検出することができる。
　さらに、図示していないが、実施例1の層20A，Bのうち、被検体100の体表に当接される接触面側の層(図2(A)の例では20A)を薄く、かつ超音波反射特性を強く形成して構成することができる。
　また、超音波散乱体を弾性カプラ20の全体に分散混入し、散乱体の分散濃度を異ならせることによりID符号化してもよい。超音波散乱体としては、弾性カプラ20の素材と音響インピーダンスが異なるグラファイト粉末やポリエチレン粉末等の材料を用いることができる。
　さらに、関心領域(ROI)から外れた弾性カプラ20のスキャン方向の両端の領域部の少なくとも一方に、スキャン方向に散乱体をコード化して分散させた符合領域を形成し、弾性カプラ20の種類に応じたID符号を付与することもできる。
　また、弾性カプラ20の実施例1～3のいずれの場合も、弾性カプラ20の全体に、超音波散乱体を混入して超音波の減衰を大きくすることができる。
　ところで、本発明の絶対圧力の検出原理は、加えられる圧迫力(圧力)に相関して弾性カプラ20の厚みが変化すること、その相関は弾性カプラ20の弾性特性に依存することに基づいている。そこで、絶対圧力を求めるには、弾性カプラ20の種類ごとに弾性特性を予め計測し、弾性カプラ20のID符号に対応付けてカプラデータベースに設定しておく必要がある。
　ここで、カプラデータベースに設定する弾性特性は、境界面又は中間層の厚み方向の変形(歪み)が絶対圧力に対して線形の場合、又は弾性カプラ20の弾性率(例えば、ヤング率)が一定(線形)の場合は、ID符号に対応させて弾性率Eを設定すればよい。
　また、境界面又は中間層の厚み方向の位置変化が絶対圧力に対して非線形の場合、弾性率Eが絶対圧力に対して非線形の場合は、厚み方向の変形(歪み)と絶対圧力の関係曲線、又厚み方向の変形(歪み)と弾性率Eの関係曲線を設定する。あるいは、厚み方向の変形又は歪みの積算値と弾性率の関係曲線、又は厚み方向の変形(歪み)に対する弾性率補正係数とすることができる。
　以上説明したように、本実施形態によれば、通常の弾性画像計測の検査の処理において、探触子1により被検体100の生体組織に加えられる絶対圧力をリアルタイムで精度よく計測できる。
　また、弾性カプラのID符号を弾性画像計測の処理の中で読み取ることができるから、種々の異なる弾性カプラを任意に交換しても、自動的に弾性カプラを識別して、弾性カプラの弾性特性に応じて絶対圧力を精度よく計測でき、操作者の手間を省いて使い勝手をよくすることができる。
　また、本実施形態では、弾性画像を生成表示する超音波診断装置を例に説明したが、本発明はこれに限られるものではなく、被検体の生体組織の断層画像(Bモード像)による診断、ドプラ計測又はカラーフローモード(CFM)による血流診断等の検査を行う超音波診断装置に適用することができる。これにより、種々の検査に適した圧迫状態の評価を精度よく実現できる。

Claims (13)

  1. 　被検体に当接させて超音波を送受する超音波探触子と、該超音波探触子を駆動する送信部と、前記超音波探触子により受信された反射エコー信号であるRF信号を受信処理する受信部と、該受信部から出力されるRF信号に基づいて超音波像を生成する画像生成部とを備えてなる超音波診断装置において、
    　前記超音波探触子の超音波送受信面に装着された弾性カプラの変形に基づいて前記被検体に加えられた圧力を求める圧力演算部を備え、
    　前記弾性カプラは、少なくとも2つの層を有して形成され、
    　前記圧力演算部は、前記2つの層間の境界面の位置を検出し、該境界面の検出位置と予め求められた境界面の初期位置とに基づいて前記境界面の位置変化を求め、該位置変化と前記弾性カプラの弾性特性に基づいて絶対圧力を求めることを特徴とする超音波診断装置。
  2. 　請求項1に記載の超音波診断装置において、
    　前記弾性カプラは、超音波反射特性が異なる少なくとも2つの層を有して形成され、前記超音波反射特性が異なる層の境界面が前記超音波送受信面に装着される装着面と前記被検体の体表に当接される接触面との間に位置させて設けられてなり、
    　前記圧力演算部は、前記受信部から出力される取得時刻が異なる一対の前記RF信号のフレームデータに基づいて、前記弾性カプラの厚み方向の前記境界面の位置の変位を検出し、該変位に基づいて前記境界面の厚み方向の歪みを求め、該厚み方向の歪みと予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求めることを特徴とする超音波診断装置。
  3. 　請求項2に記載の超音波診断装置において、
    　前記圧力演算部は、前記弾性カプラが無加圧の初期状態から時間経過に合わせて前記境界面の歪みを積算し、該歪みの積算値と予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求めることを特徴とする超音波診断装置。
  4. 　請求項1に記載の超音波診断装置において、
    　前記弾性カプラは、前記境界面に中間層を挟んで形成され、該中間層の超音波反射特性は他の層と異ならせて形成されてなり、
    　前記圧力演算部は、前記中間層の厚み方向の位置を検出し、該中間層の検出位置と予め求められた中間層の初期位置とに基づいて前記中間層の位置変化を求め、該位置変化と予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求めることを特徴とする超音波診断装置。
  5. 　請求項2に記載の超音波診断装置において、
    　前記弾性カプラは、前記境界面に中間層を挟んで形成され、該中間層の超音波反射特性は他の層と異ならせて形成されてなり、
    　前記圧力演算部は、前記受信部から出力される取得時刻が異なる一対の前記RF信号フレームデータに基づいて、前記弾性カプラの前記中間層の厚み方向の位置の変位を検出し、該変位に基づいて前記中間層の厚み方向の歪みを求め、該厚み方向の歪みと予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求めることを特徴とする超音波診断装置。
  6. 　請求項3に記載の超音波診断装置において、
    　前記弾性カプラは、前記境界面に中間層を挟んで形成され、該中間層の超音波反射特性は他の層と異ならせて形成されてなり、
    　前記圧力演算部は、前記弾性カプラが無加圧の初期状態から時間経過に合わせて前記中間層の厚み方向の歪みを積算し、該歪みの積算値と予め設定された前記弾性カプラの弾性特性に基づいて前記被検体に加えられる絶対圧力を求めることを特徴とする超音波診断装置。
  7. 　請求項4乃至6のいずれか一項に記載の超音波診断装置において、
    　前記中間層は、前記超音波探触子の超音波ビームのスキャン方向に直交する方向に延在させて、かつ間隔をあけて設けられた線状の複数の超音波反射体からなることを特徴とする超音波診断装置。
  8. 　請求項1又は2のいずれかに記載の超音波診断装置において、
    　前記被検体の体表に当接される接触面側の層は、前記超音波送受信面に当接される接触面側の層より薄い、又は前記超音波反射特性が強く形成されてなることを特徴とする超音波診断装置。
  9. 　請求項1に記載の超音波診断装置において、
    　前記弾性カプラは、前記弾性材料の中に超音波散乱体が混入され、超音波減衰特性を大きく形成してなることを特徴とする超音波診断装置。
  10. 　請求項1又は2のいずれかに記載の超音波診断装置において、
    　前記弾性カプラは、前記境界面の厚み方向の位置を当該弾性カプラの種類に応じて異ならせて、前記RF信号により当該弾性カプラの種類を識別可能な識別符合が形成され、　前記圧力演算部は、前記RF信号又は前記RF信号フレームデータに基づいて前記識別符合を検出して前記弾性カプラの種類を識別し、前記弾性カプラの種類に対応して設定されている前記弾性特性に基づいて前記絶対圧力を求めることを特徴とする超音波診断装置。
  11. 　請求項10に記載の超音波診断装置において、
    　前記識別符号は、前記境界面の厚み方向の位置を異ならせて形成され、
    　前記圧力演算部は、前記RF信号又は前記RF信号フレームデータに基づいてカプラエコー領域内のRF信号の深度分布パターンにより、弾性カプラの種類を識別することを特徴とする超音波診断装置。
  12. 　請求項10に記載の超音波診断装置において、
    　前記識別符号は、前記弾性カプラのスキャン方向の両端の領域部に、スキャン方向と厚み方向の少なくとも一方に超音波散乱体を符号化して分散させて形成され、
    　前記圧力演算部は、前記RF信号又は前記RF信号フレームデータに基づいて前記弾性カプラの両端の領域部のRF信号のパターンにより、弾性カプラの種類を識別することを特徴とする超音波診断装置。
  13. 　請求項1に記載の超音波診断装置において、
    　前記弾性カプラの弾性特性は、弾性率、厚み方向の変形と弾性率の関係曲線、厚み方向の歪みと弾性率の関係曲線、厚み方向の変形又は歪みの積算値と弾性率の関係曲線、厚み方向の変形又は歪みに対する弾性率補正係数の少なくとも一つであることを特徴とする超音波診断装置。
PCT/JP2009/057470 2008-04-25 2009-04-14 超音波診断装置 WO2009131028A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2010509144A JP5329532B2 (ja) 2008-04-25 2009-04-14 超音波診断装置
US12/989,527 US8747320B2 (en) 2008-04-25 2009-04-14 Ultrasonic diagnostic apparatus

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2008116388 2008-04-25
JP2008-116388 2008-04-25

Publications (1)

Publication Number Publication Date
WO2009131028A1 true WO2009131028A1 (ja) 2009-10-29

Family

ID=41216765

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2009/057470 WO2009131028A1 (ja) 2008-04-25 2009-04-14 超音波診断装置

Country Status (3)

Country Link
US (1) US8747320B2 (ja)
JP (1) JP5329532B2 (ja)
WO (1) WO2009131028A1 (ja)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012042794A1 (ja) * 2010-09-27 2012-04-05 パナソニック株式会社 超音波診断用アダプタ、超音波診断装置、及び、超音波診断方法
JP2012176197A (ja) * 2011-02-28 2012-09-13 Nagasaki Univ フィルム状エコーゲル及び超音波センサユニット
JP2014135962A (ja) * 2013-01-15 2014-07-28 Seiko Epson Corp 超音波測定装置、超音波診断装置及び超音波測定用シート
JP5920746B1 (ja) * 2015-01-08 2016-05-18 学校法人早稲田大学 穿刺支援システム
WO2019167592A1 (ja) * 2018-03-01 2019-09-06 有限会社ユーマンネットワーク 平面形素子を用いた集束式音波治療装置
KR20200042228A (ko) * 2018-10-15 2020-04-23 한국과학기술연구원 연조직의 물성을 정량적으로 측정하기 위한 시스템
JP2020114289A (ja) * 2019-01-17 2020-07-30 伊藤超短波株式会社 超音波デバイス用固定具

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011044421A1 (en) 2009-10-08 2011-04-14 C. R. Bard, Inc. Spacers for use with an ultrasound probe
CA2806353A1 (en) 2010-08-09 2012-02-16 C.R. Bard Inc. Support and cover structures for an ultrasound probe head
JP2013116138A (ja) * 2011-12-01 2013-06-13 Sony Corp 画像処理装置及び方法
WO2013188833A2 (en) * 2012-06-15 2013-12-19 C.R. Bard, Inc. Apparatus and methods for detection of a removable cap on an ultrasound probe
WO2015149805A1 (en) * 2014-04-04 2015-10-08 Aarhus Universitet Reference material element for ultrasound scanning probe
RU2707831C2 (ru) * 2014-11-25 2019-11-29 Конинклейке Филипс Н.В. Система для контроля использования медицинского устройства
KR102669618B1 (ko) 2016-04-12 2024-05-27 삼성메디슨 주식회사 초음파 진단 장치 및 그 제어 방법
WO2023186550A1 (en) * 2022-03-31 2023-10-05 Koninklijke Philips N.V. Generating an ultrasound device interaction indicator
EP4252667A1 (en) * 2022-03-31 2023-10-04 Koninklijke Philips N.V. System and method for generating an ultrasound device interaction indicator

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005066041A (ja) * 2003-08-25 2005-03-17 Hitachi Medical Corp 超音波探触子及び超音波診断装置
WO2005120358A1 (ja) * 2004-06-09 2005-12-22 Hitachi Medical Corporation 弾性画像表示方法及び超音波診断装置
JP2006247203A (ja) * 2005-03-11 2006-09-21 Matsushita Electric Ind Co Ltd ３次元超音波探触子及び３次元超音波診断装置
WO2006121031A1 (ja) * 2005-05-09 2006-11-16 Hitachi Medical Corporation 超音波診断装置及び超音波画像表示方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19637392A1 (de) * 1996-09-13 1998-03-26 Siemens Ag Drucksensor unter Verwendung von mit akustischen Oberflächenwellen arbeitenden Elementen - OFW-Elemente -
WO2001082225A2 (en) * 2000-04-24 2001-11-01 Washington University Method and apparatus for probabilistic model of ultrasonic images
US6500126B1 (en) * 2001-12-20 2002-12-31 Koninklijke Philips Electronics N.V. Ultrasound system transducer adapter
US7369458B2 (en) * 2004-05-10 2008-05-06 Airmar Technology Corporation Transducer identification
CN101489490B (zh) * 2006-07-31 2011-03-23 株式会社日立医药 压迫装置及使用该压迫装置的超声波探头和超声波诊断装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005066041A (ja) * 2003-08-25 2005-03-17 Hitachi Medical Corp 超音波探触子及び超音波診断装置
WO2005120358A1 (ja) * 2004-06-09 2005-12-22 Hitachi Medical Corporation 弾性画像表示方法及び超音波診断装置
JP2006247203A (ja) * 2005-03-11 2006-09-21 Matsushita Electric Ind Co Ltd ３次元超音波探触子及び３次元超音波診断装置
WO2006121031A1 (ja) * 2005-05-09 2006-11-16 Hitachi Medical Corporation 超音波診断装置及び超音波画像表示方法

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012042794A1 (ja) * 2010-09-27 2012-04-05 パナソニック株式会社 超音波診断用アダプタ、超音波診断装置、及び、超音波診断方法
CN102665568A (zh) * 2010-09-27 2012-09-12 松下电器产业株式会社 超声波诊断用适配器、超声波诊断装置以及超声波诊断方法
JP5888229B2 (ja) * 2010-09-27 2016-03-16 コニカミノルタ株式会社 超音波診断用アダプタ、超音波診断装置、及び、超音波測定方法
US9445783B2 (en) 2010-09-27 2016-09-20 Konica Minolta, Inc. Ultrasound diagnostic adapter, ultrasound diagnostic apparatus, and ultrasound diagnostic method
JP2012176197A (ja) * 2011-02-28 2012-09-13 Nagasaki Univ フィルム状エコーゲル及び超音波センサユニット
JP2014135962A (ja) * 2013-01-15 2014-07-28 Seiko Epson Corp 超音波測定装置、超音波診断装置及び超音波測定用シート
JP2016123794A (ja) * 2015-01-08 2016-07-11 学校法人早稲田大学 穿刺支援システム
JP5920746B1 (ja) * 2015-01-08 2016-05-18 学校法人早稲田大学 穿刺支援システム
US10004532B2 (en) 2015-01-08 2018-06-26 Waseda University Puncture assistance system
WO2019167592A1 (ja) * 2018-03-01 2019-09-06 有限会社ユーマンネットワーク 平面形素子を用いた集束式音波治療装置
KR20200042228A (ko) * 2018-10-15 2020-04-23 한국과학기술연구원 연조직의 물성을 정량적으로 측정하기 위한 시스템
KR102174348B1 (ko) 2018-10-15 2020-11-04 한국과학기술연구원 연조직의 물성을 정량적으로 측정하기 위한 시스템
JP2020114289A (ja) * 2019-01-17 2020-07-30 伊藤超短波株式会社 超音波デバイス用固定具

Also Published As

Publication number Publication date
JPWO2009131028A1 (ja) 2011-08-18
JP5329532B2 (ja) 2013-10-30
US8747320B2 (en) 2014-06-10
US20110040186A1 (en) 2011-02-17

Similar Documents

Publication Publication Date Title
JP5329532B2 (ja) 超音波診断装置
JP5329533B2 (ja) 超音波診断装置
US8105240B2 (en) Ultrasonic imaging apparatus and low attenuation medium with a prescribed pattern for apparatus localization
JP4769715B2 (ja) 超音波診断装置および弾性画像表示方法
JP5028416B2 (ja) 超音波診断装置
JP4455003B2 (ja) 超音波診断装置
JP4966578B2 (ja) 弾性画像生成方法及び超音波診断装置
US20110194748A1 (en) Ultrasonic diagnostic apparatus and ultrasonic image display method
WO2009098961A1 (ja) 超音波診断装置
WO2009104657A1 (ja) 超音波診断装置、超音波画像処理方法及び超音波画像処理プログラム
JP2005066041A5 (ja)
JP2004135929A (ja) 超音波診断装置
JP5156421B2 (ja) 超音波診断装置
JP2007105400A (ja) 超音波診断装置及び画像処理装置
JP5113322B2 (ja) 超音波診断装置
WO2009104525A1 (ja) 超音波診断装置、超音波弾性情報処理方法及び超音波弾性情報処理プログラム
JP5473527B2 (ja) 超音波診断装置
JP5523019B2 (ja) 弾性映像を形成する超音波システム及び弾性映像形成方法
JP5623609B2 (ja) 超音波診断装置
CN109069117B (zh) 超声波诊断装置
JP4368185B2 (ja) 超音波診断装置
JP5638641B2 (ja) 超音波診断装置
JP6169396B2 (ja) 参照変形体、超音波探触子、及び超音波撮像装置
JP2013183982A (ja) 超音波診断装置及び弾性画像生成方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 09735346

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2010509144

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 12989527

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 09735346

Country of ref document: EP

Kind code of ref document: A1