WO2001026555A1 - Dispositif d'imagerie ultrasonore - Google Patents

Dispositif d'imagerie ultrasonore Download PDF

Info

Publication number
WO2001026555A1
WO2001026555A1 PCT/JP2000/007114 JP0007114W WO0126555A1 WO 2001026555 A1 WO2001026555 A1 WO 2001026555A1 JP 0007114 W JP0007114 W JP 0007114W WO 0126555 A1 WO0126555 A1 WO 0126555A1
Authority
WO
WIPO (PCT)
Prior art keywords
ultrasonic
transducer
delay time
delay
sound
Prior art date
Application number
PCT/JP2000/007114
Other languages
English (en)
Japanese (ja)
Inventor
Ryuichi Shinomura
Takashi Azuma
Shinichiro Umemura
Yuichi Miwa
Hiroshi Kanda
Hirotaka Baba
Tatsuya Hayashi
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 JP2001529349A priority Critical patent/JP4711583B2/ja
Publication of WO2001026555A1 publication Critical patent/WO2001026555A1/fr

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461—Displaying means of special interest
    • A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13—Tomography
    • A61B8/14—Echo-tomography
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52046—Techniques for image enhancement involving transmitter or receiver
    • G01S7/52049—Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/18—Methods or devices for transmitting, conducting or directing sound
    • G10K11/26—Sound-focusing or directing, e.g. scanning
    • G10K11/34—Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
    • G10K11/341—Circuits therefor
    • G10K11/346—Circuits therefor using phase variation

Definitions

  • the present invention relates to an ultrasonic imaging apparatus for extracting a tissue in a subject as an image using ultrasonic waves, and particularly to an ultrasonic beamformer used for the same.
  • An ultrasonic device for example, an ultrasonic imaging device used for medical image diagnosis obtains a tomographic image of a soft tissue of a living body or an image of a blood flow flowing through the living body in almost real time using an ultrasonic pulse reflection method. It can be displayed on a monitor and observed, and it is said to be highly safe because it does not expose the subject to radiation, unlike a diagnostic imaging device that uses radiation. Applied in the field.
  • an ultrasonic probe is used for transmitting an ultrasonic wave into a subject and receiving an echo signal from inside the subject.
  • One of the scanning methods of the ultrasonic imaging apparatus is an electronic scanning method.
  • the electronic scanning type ultrasonic probe arranges elongated rod-shaped transducers in a one-dimensional array, and gives a predetermined delay time to each transducer for driving.
  • the ultrasonic probe transmits an ultrasonic beam converging at a predetermined depth and a predetermined direction in the subject.
  • the received wave is obtained by synthesizing a received signal from each transducer of the ultrasonic probe with a given delay time for each transducer and combining the received signals from a given depth and direction.
  • the processing part that sets the delay time for each transducer and gives the delay time to each transducer is called an ultrasonic beamformer.
  • the ultrasonic probe uses ultrasonic waves generated by many transducers arranged in a one-dimensional array. A lens layer for converging in a direction orthogonal to the arrangement direction is provided.
  • the phasing circuit of the ultrasonic imaging apparatus since the phasing circuit of the ultrasonic imaging apparatus has recently been digitized, it has become possible to easily and accurately control the delay time and phasing in the transmission and reception of ultrasonic waves.
  • the ultrasonic waves transmitted from each transducer do not converge at the focal point only through a medium with a constant sound velocity.
  • the transmitted pulse is converted from an electric signal into an ultrasonic wave by the piezoelectric vibrator, and reaches a desired focal position through the lens layer on the probe surface and the living tissue.
  • sound refraction occurs at the interface between the two in the arrangement direction of the piezoelectric vibrators.
  • the delay time it is necessary to incorporate the effect of this refraction. Since the refraction angle is obtained by Snell's law, it is possible to calculate the delay time by obtaining the refraction path of the sound.
  • the ultrasonic imaging device captures tomographic images of about 30 frames per second while moving the focus position during the capturing of one frame by 50 in the horizontal direction (azimuth direction) and in the depth direction (distance direction). Needs to be changed by about 20.
  • the number of transducers in a single aperture varies from 32 to 192, depending on the transducer. From these facts, it is necessary to set the delay time for one vibrator at a very high speed even if it is calculated. However, there is a limit to the amount of computation that can be performed in real-time with the MPU (Micro Processor Unit) mounted on a real ultrasonic imaging device, even in light of the arithmetic processing performance of a DSP (Digital Signal Processor).
  • MPU Micro Processor Unit
  • a method of storing data in which the delay time for each focus position is calculated in advance for each probe in a storage medium such as a hard disk, or a method for real-time processing is to perform an approximate calculation using a formula.
  • probes that are specialized for the organ, site, and symptom to be diagnosed are used as probes of the ultrasonic imaging apparatus, and the types thereof are increasing. Therefore, a considerably large-capacity storage device is required to store the delay time data.
  • the number of transducers will increase further. Data amount is increasing.
  • the method of obtaining the delay time by approximation calculation may be insufficient in accuracy with the conventional calculation accuracy depending on conditions. Since the accuracy is determined by the ratio between the period at the center frequency of the probe and the error time, a probe that transmits a pulse with a high center frequency, that is, a pulse with a short period, like a high-definition probe has been developed. The demands on the accuracy of the delay time become stricter. Therefore, a more accurate approximation method is required. In addition, considering the application to a two-dimensional probe, it is necessary to shorten the calculation time for one transducer, so going in a more severe direction is the same as the former case.
  • the image used for ultrasonic diagnosis has been significantly improved in comparison with the conventional image in combination with the adoption of the digital phasing technology, but the image obtained by the ultrasonic imaging apparatus is an X-ray apparatus. There is a demand for further improvement in image quality when compared with images obtained by other modalities such as X-ray CT and MRI.
  • the delay time control data for beam formation is calculated based on the measured or estimated value. Is what you want.
  • This sound velocity correction technique is described in, for example, Japanese Patent Application Laid-Open Nos. 8-317923 and 10-66694 filed by the present applicant.
  • the inventors of the present invention believe that this can be improved by forming a beam by incorporating the refraction in the ultrasonic wave propagation path.
  • the refraction in the ultrasonic wave propagation path the refraction by the lens layer is the first problem, but the refraction of the ultrasonic wave occurs only at the interface between the lens layer and the subject. Les ,.
  • the living body to be diagnosed by the ultrasonic imaging device is composed of various tissues such as fat, muscle fibers, various organs, and blood. Since the speed of sound is different in each of these tissues, subtle refraction of ultrasonic waves occurs at the interface. Of these refractions, it is easiest to incorporate the refraction effect of the fatty layer into beamforming.
  • the ultrasonic pulse penetrates the fat layer when the ultrasound probe is applied percutaneously to see any of the skin. In other words, refraction by the fat layer cannot be avoided unless a transcutaneous probe is used.
  • the fat layer is the outermost part of the living tissue, the effect of refraction on the focal position has the greatest effect.
  • the subcutaneous fat has a certain thickness within the caliber of the probe, so it is easier to incorporate it into the calculation than in the case of other tissues such as blood vessels. is there. Disclosure of the invention
  • the present invention has been made in view of the above, and a first object of the present invention is to further improve the image quality of an ultrasonic image compared to the current state.
  • a second object of the present invention is to provide an ultrasonic imaging apparatus capable of setting a delay time in consideration of the effect of ultrasonic refraction by a lens layer and / or a fat layer of an ultrasonic probe. It is in.
  • a third object of the present invention is to set a delay time in consideration of the influence of refraction by the lens layer or / and fat layer and the ultrasonic wave propagation speed of each subject.
  • An object of the present invention is to provide an ultrasonic imaging apparatus.
  • the present invention employs a method of obtaining a delay time by using an algorithm having a sufficient approximation accuracy and a fast calculation, and taking into account the influence of ultrasonic refraction by the lens layer and the fat layer.
  • an ultrasonic imaging apparatus includes: an ultrasonic probe having an arrayed transducer; and transmitting, focusing, and receiving when transmitting or Z and receiving ultrasonic waves to a subject.
  • Delay control means for controlling a delay time with respect to each transducer for performing focusing; and the transmitting or transmitting by incorporating a refraction effect of ultrasound by an ultrasound propagation medium between the arrayed transducers and a set focal position. It is characterized by comprising a refraction correction delay data generating means for generating a delay time for performing forcing of a received wave and supplying the delay time to the delay control means, and a display unit for displaying an ultrasonic image.
  • the delay control means stores in advance delay time data obtained based on the average sound velocity of the living body, and performs ultrasonic transmission and reception for obtaining refraction correction data using the stored delay time data in advance.
  • the refraction correction delay data generating means uses parameters relating to the ultrasonic probe including a lens layer thickness of the probe, a sound velocity of the lens layer, and an arrangement pitch of the transducer, and calculates a relationship between the transducer and a designated focal point.
  • the delay time to be given to each transducer is calculated by calculation, taking into account the ultrasonic refraction effect in the ultrasonic propagation path between them.
  • the refraction correction delay data generating means includes: a parameter relating to the ultrasonic probe including a lens layer thickness, a lens layer sound velocity, a pitch between transducers; a fat layer thickness of the subject and a sound velocity of the fat layer; The delay time to be given to each transducer is calculated by using the sound velocity data of the weave, taking into account the ultrasonic refraction effect in the ultrasonic propagation path between the transducer and the designated focal point. Further, the refraction correction delay data generation control means may calculate the delay time using a parameter recursively solved from a parameter relating to a sound path from the vibrator next to the vibrator to be calculated to the focal point. Obtained by calculation.
  • the ultrasonic imaging apparatus of the present invention comprises: means for detecting data relating to a layer structure of a subject for refraction correction from a screen of the display unit on which an ultrasonic image is displayed; and an output of the layer structure data detecting means.
  • the layer structure data detecting means includes a caliper for displaying two movable cursors on the screen and measuring a distance between the cursors on the screen.
  • the ultrasonic imaging apparatus of the present invention comprises: an ultrasonic probe having an arrayed transducer; Delay control means for controlling a delay time for each transducer to perform transmission focusing or reception focusing at the time of transmitting or z and receiving ultrasonic waves to the subject, and a structure forming a layer in the subject.
  • Means for measuring the thickness means for measuring the speed of sound of the portion of the layer structure, and the layer thickness measured by the layer thickness measuring means and the sound speed in the layer structure measured by the sound speed measuring means.
  • a delay time to be given to each transducer in consideration of an ultrasonic refraction effect in an ultrasound propagation path between the transducer and a designated focal point, and a refraction correction delay control to be supplied to the delay control means.
  • the sound velocity measuring means calculates the delay time error of each receiving channel by using the output of a delay circuit that performs phasing processing of the echo signals received by the plurality of transducers, and calculates the delay time error in the subject from the delay error. Includes sound speed measurement means for determining sound speed. Further, the sound velocity measuring means includes means for specifying a sound velocity measuring region to a fat layer of the subject and a tissue portion at a part of the fat layer, and measuring a sound velocity for each of the tissue layers.
  • the ultrasonic imaging apparatus of the present invention uses the parameters related to the ultrasonic probe including the lens layer thickness of the ultrasonic probe, the sound velocity of the lens layer, and the transducer array pitch, and the ultrasonic wave exits the transducer.
  • a program to make a computer execute a method to calculate the delay time to be applied to each transducer taking into account the effect of refraction from when the lens reaches the specified focal point, or the lens layer thickness, the sound velocity of the lens layer
  • the ultrasonic The program has a built-in program that allows a computer to execute a method of calculating the delay time given to each transducer taking into account the effect of refraction from exiting the transducer until reaching the designated focal point. are doing.
  • FIG. 1 is a block diagram showing a schematic configuration example of an ultrasonic imaging apparatus according to an embodiment of the present invention
  • FIG. 2 is a block diagram showing an embodiment of an ultrasonic beamformer without the ultrasonic imaging apparatus shown in FIG. 1
  • FIG. 3 is a diagram showing the relationship between the transducer array and the ultrasonic image
  • FIG. 4 is a diagram showing the relationship between the transducer array and the delay time during the transmission and reception of the ultrasound
  • FIG. FIG. 6 is a flowchart showing a flow of tomographic image imaging to which an example of refraction correction is applied.
  • FIG. 6 is a diagram showing an example of a method for measuring a layer thickness based on an image displayed on an image display unit.
  • FIG. 1 is a block diagram showing a schematic configuration example of an ultrasonic imaging apparatus according to an embodiment of the present invention
  • FIG. 2 is a block diagram showing an embodiment of an ultrasonic beamformer without the ultrasonic imaging apparatus shown in FIG. 1
  • FIG. 3
  • FIG. 7 is a diagram showing refraction of an ultrasonic pulse.
  • FIG. 8 is a flowchart illustrating a first example of a procedure for calculating a delay time for refraction correction according to the present invention.
  • FIG. 9 is a flowchart illustrating a procedure for calculating a delay time for refraction correction according to the present invention.
  • FIG. 10 is a flowchart for explaining a second example, FIG. 10 is a diagram for explaining a relationship between a focal position and a transmitting aperture position in a phased array probe, FIG. 11 is an explanatory diagram of a refraction state of an ultrasonic pulse, and FIG. Illustration of the discrete Newton method, FIG.
  • FIG. 13 is a flowchart illustrating a third example of the procedure of delay time calculation for refraction correction according to the present invention
  • FIG. 14 is a procedure of delay time calculation for refraction correction according to the present invention.
  • Flowchart illustrating a fourth example of FIG. 15 is an explanatory diagram of a refraction state of an ultrasonic pulse in a convex type probe
  • FIG. 16 is a block diagram showing a schematic configuration of another embodiment of the ultrasonic imaging apparatus according to the present invention
  • FIG. FIG. 18 is a block diagram showing a schematic configuration of still another embodiment of the ultrasonic imaging apparatus
  • FIG. 18 is a graph showing a time required to reach a focal point from each transducer
  • FIG. 19 is provided with a switch for a refraction correction function.
  • FIG. 20 is a block diagram showing a device configuration for obtaining a sound speed for refraction correction
  • FIG. 21 is a diagram showing a screen display example of sound speed measurement.
  • FIG. 1 is a schematic configuration diagram illustrating an example of an ultrasonic imaging apparatus according to the present invention.
  • the ultrasonic imaging apparatus 100 includes a main body 10, an ultrasonic probe 20, and a cap.
  • the main body 10 of the ultrasonic imaging apparatus includes a control unit 11, a memory 12, an image display unit 13, and an ultrasonic beam former 14.
  • the control unit 11 is connected to an input unit 15 such as a keyboard and a pointing device.
  • the ultrasonic probe 20 is detachable from the main body 10 of the ultrasonic imaging apparatus, and an appropriate probe is selected and attached according to a diagnosis target of the subject 30.
  • the ultrasonic probe 20 includes a transducer array 21 and a memory 22.
  • the memory 12 of the main body 10 stores a table of data unique to the probe corresponding to the ID number of the probe, and data such as a sound velocity in a subcutaneous fat layer and a sound velocity in a living tissue.
  • the transducer-specific data referred to here includes, for example, the thickness of the lens layer, the sound velocity in the lens layer, and the transducer pitch. In the case of a hive or convex type probe, it refers to the radius of the probe, and in the case of an oblique type probe, it refers to parameters such as the oblique angle.
  • the ultrasonic beam former 14 controls the operation of the transducer array 21 in the ultrasonic probe 20 under the control of the control unit 11.
  • Ultrasonic image data is obtained by scanning an ultrasonic beam in a subject using this ultrasonic beamformer.
  • the control unit 11 displays an ultrasonic tomographic image of the subject 30 on the image display unit 13 based on the obtained ultrasonic image data.
  • FIG. 2 is a block diagram illustrating an example of an ultrasonic beamformer.
  • the ultrasonic beamformer includes an operation control circuit 41, a pulser 42, a preamplifier and an AD converter 43, a delay circuit 44, and an addition circuit 45.
  • a delay time is calculated and controlled by a calculation and control circuit 41.
  • the arithmetic and control circuit 41 may be realized by a computer and software. In this case, a program for delay time calculation is recorded on a recording medium such as a ROM, and the calculation and control circuit 41 is realized by a configuration in which the program is read by a computer.
  • the beamforming mechanism for transmitting and receiving is as follows. First, each of the transducers 25 of the transducer array 21 included in the ultrasonic probe 20 from the arithmetic and control circuit 41. , 25 ,, 25 2> ..., 25 n , a pulse signal is sent at a timing shifted by the delay time of each oscillator. Each pulser 42 that receives the pulse signal immediately sends the pulse signal to the transducer 25 connected to itself. , 25 ⁇ 25 2 ,. Send to Transducer 25. , 25 1 ; 25 2 ,..., 25 n generate ultrasonic waves corresponding to the voltage due to their piezoelectricity, and transmit the ultrasonic waves into the subject 30.
  • the ultrasonic waves reflected in the subject 30 return to the transducer array 21 of the ultrasonic probe 20 again, and the individual transducers 25. , 25 "25 2) ⁇ , is converted into an electric signal by a piezoelectric at 25 n, are amplified by preamplifier and AZD converter 43 is converted into a digital signal.
  • the digital signal, the arithmetic and control circuit 41 The delay time is adjusted by the delay circuit 44 adjusted by the signal from the controller and the signal is added to the adder circuit 45 to be added, whereby the phasing addition is performed.
  • the ultrasonic image I is a point F nl , F n2 , F n3 , (where n is an integer from 1 to N) It consists of N received beam data. This point F nl , F n2 , F n3 ,..., Is measured by a single ultrasonic beam.
  • FIG. 3 if I is an ultrasonic image, the ultrasonic image I is a point F nl , F n2 , F n3 , (where n is an integer from 1 to N) It consists of N received beam data. This point F nl , F n2 , F n3 ,..., Is measured by a single ultrasonic beam.
  • FIG. 4 is a diagram schematically showing an ultrasonic probe that transmits an ultrasonic pulse toward the focal point F dockand detects a reflected signal from the focal point F religious.
  • Figure 4 shows how the wavefront moves at regular time intervals Ts.
  • the ultrasonic pulse transmitted from the transducer 25 within the bore of the ultrasonic probe is used.
  • An ultrasonic wave is transmitted with a delay time for each vibrator so that it converges at the focal point Flust.
  • transducers 25 located at both ends of the aperture.
  • the delay time ⁇ given to the j-th oscillator 25” is the oscillator 25 as shown. Is the time until a virtual ultrasound wavefront WA at position 25 n reaches the oscillation surface of the vibrator 25 ".
  • the wavefront WB of the reflected wave propagates as shown in Fig. 4 (b), and the oscillators 25 at both ends. , 25n, and arrives at the transducer surface of the transducer 25 ”with a delay time ⁇ ”.
  • the child adds the phased signals by shifting the minute delay time in each of the ultrasonic probe diameter in the vibrator 25 0 to 25 n, it is possible to pick up the reflected signal from the focal point F " .
  • this transmission / reception operation is a method known as a dynamic focus method, that is, transmission is performed by focusing on a specific depth located in a region of interest (R0I: Region of Interest) within a subject. Transmits a sound wave, and changes the focal position of the received wave from shallow to deep, from region to region as the transmitted ultrasonic pulse travels through the subject during reception. Is used to obtain the received beam signal of When the measurement of one beam is completed, the beam position is sequentially shifted to the adjacent direction (azimuth direction), and the same transmission / reception operation is repeated to obtain N reception beam signals to form an image.
  • R0I Region of Interest
  • the former is a scanning method used for sector-type probes. In this method, all transducers of the probe are used for transmission and reception every time, and the beam direction is changed radially. In this case, when the beam direction is in relation to the transmission aperture, in other words, the focal point of transmission and reception is slightly in front of the aperture. This has great significance in the calculation of the delay time later.
  • the latter is a method used for a linear probe and a convex probe.
  • This method differs from the former in that it does not use all transducers provided in the probe for one transmission.
  • a probe with a total number of transducers of 192 channels and a transmission of 64 channels in diameter Is used.
  • transmission and reception are performed with the aperture positioned in front of the beam, and movement of the beam position is performed by moving the aperture, so that the focus is always in front of the aperture, contrary to the former case. How to transmit waves.
  • step 11 recognition of the ultrasonic probe 20 attached to the main body 10 of the ultrasonic imaging apparatus is performed. I do. This recognition is performed by reading a parameter such as a probe ID number stored in the memory 22 of the ultrasonic probe 20.
  • FIG. 6 is a schematic diagram of a diagnostic screen displayed on the image display unit 13 at this time, and an image of a somewhat blurred target is displayed because no correction is applied.
  • the image shown in Fig. 6 is an image taken by a convex type probe, and a subcutaneous fat layer is shown on the side in contact with the probe.
  • step 13 the diagnostician measures and inputs the thickness of this fat layer. Specifically, the measurement is performed using a caliper conventionally provided in an ultrasonic imaging apparatus. In Fig. 6, the cursor of the caliper is positioned at the point where the fat layer indicated by A starts and the point where the fat layer indicated by B ends. Is determined. (In the example shown, 13 corrupt) This value can be used as the thickness of the fat layer. Fatty layer thickness entered by the diagnostician After that, all processes are executed in the diagnostic device.
  • the imaging mode is entered, and the diagnostician instructs the beamformer 14 from the input unit 15 via the control unit 11 of the ultrasonic imaging apparatus to the beamformer 14 in step 14.
  • the beamformer 14 calculates the delay time after the refractive index correction described later has been performed by the arithmetic and control circuit 41 in step 15.
  • the delay time calculated here is the delay time given to each transducer for transmission corresponding to the focal depth of the transmission, and the delay given to each transducer to form one ultrasonic reception beam. This is the time, that is, the delay time for continuously changing the focal position of the received wave as described above.
  • step 16 when the transducer array of the ultrasonic probe is excited using the determined transmission delay time, the ultrasonic array is converged to the focal point specified on the first ultrasonic beam line. An ultrasonic beam is transmitted.
  • step 17 the received signal from the transducer array of the ultrasonic probe is phased using the receiving delay time calculated in step 15, and added to obtain the first received beam signal. can get.
  • step 18 it is determined whether signal acquisition from all points in the imaging screen has been completed. That is, when an ultrasonic image is formed by N ultrasonic beams, it is determined whether an ultrasonic beam of the Nth address has been obtained. If not completed, the flow returns to step 14 to change the beam position and perform transmission / reception. This is repeated, and when the scanning of the entire imaging range has been completed for the scanning in the depth direction and the scanning in the horizontal direction, the diagnostic images have been captured. After taking the diagnostic image, proceed to step 19, display it on the image display section 13, and end. In the image display, instead of displaying all the data for all points of a single diagnostic image after completing the data collection, each time data for one received beam is obtained in step 17, it is displayed as an image. You may make it display in the part 13 sequentially. This is possible by utilizing the function of the digital 'scan' converter provided in the ultrasonic imaging device.
  • the ultrasonic probe is provided with a transducer row arranged one-dimensionally at a transducer pitch p, and a lens layer having a thickness arranged in front of the transducer row.
  • the lens layer is for converging the ultrasonic waves generated from the transducer row in a direction perpendicular to the transducer row (a direction perpendicular to the plane of FIG. 7).
  • the sound velocity in the living tissue existing inside the fat layer of the living body is uniform on average, and that the fat layer has a certain thickness within the transmission aperture of the ultrasonic probe.
  • the ultrasonic wave transmitted from the transducer propagates through the three layers of the lens layer, the fat layer, and the living tissue to reach the focal point, and the refraction becomes a three-layer problem.
  • the electric pulse signal is converted into ultrasonic waves by a piezoelectric vibrator 25 "having a matching layer on the front surface.
  • the ultrasonic waves generated by the vibrator 25 '' travel through the lens layer of thickness by a distance Xl in the array direction of the vibrator row, enter the subcutaneous fat layer at an incident angle, and at the boundary between the lens layer and the fat layer.
  • the refracted (refraction angle of 0 2 ) ultrasonic wave travels through the thick subcutaneous fat layer by a distance x 2 in the direction of the array of transducer rows, and enters the living tissue below the fat layer at an incident angle ⁇ 2 I do.
  • the ultrasonic wave refracted at the interface between the fat layer and the living tissue travels in the living tissue in a diagonal direction of ⁇ 3 in the array direction of the vibrator row and d 3 in the depth direction. Reach the focus within.
  • Equation 1 is established between the refraction angle 0 ⁇ ⁇ , and the above.
  • the sound speed C is determined by the material of the lens, and can be measured. It is known. In the following description, the sound speeds c 2 and c 3 are also empirically assumed to be known as approximate values, and the following description is advanced. An embodiment will be described later. Incidentally, speed of sound Cl in the memory 22 of the ultrasonic probe 20, also, the sound velocity c 2, c 3 are assumed to be held in memory 12 of the ultrasonic imaging apparatus main body 10.
  • the variables left behind when the variables are dropped are the lateral displacement Xl in the lens layer and the lateral displacement x in the fat layer x 2, there is a choice about the or leave any lateral movement amount x 3 in a living body, a most advantageous and Do Runowa chi 3 under a number of conditions when considering the magnitude of the error This is because the speed of sound gradually increases in the lens layer, fat layer, and living tissue, so that the path of the sound wave passing through each layer gradually lays down, and under many imaging conditions.
  • the thickness is the largest, so in the actual system, the sound path in the living tissue has the largest lateral movement. Kikunaru -.. Since turn seek directly computing a large amount of Ru method der to reduce most the relative error, in the normal imaging conditions become better by obtaining the x 3 and teeth force, regardless of the condition the lens layer thin top, thickness imaging conditions thickness and body tissue of the fat layer in relation to the thickness of the focal length of different is the depth of the imaging site. Yotsute fat layer, the x 2 There may be an algorithm that considers whether to leave or leave x 3 and then solve a polynomial of one variable, but here, the explanation of the method is omitted and only the method of obtaining x 3 will be described. I do.
  • Equation 3 Eliminate X x 2 from the simultaneous equations to obtain Equation 3.
  • the calculation unit of the ultrasonic beamformer in the actual ultrasonic imaging device uses the MPU or DSP as described above.
  • DSP digital signal processor
  • all calculations are converted to sum of products.
  • the square root is equivalent to 32 sums of products. Since the division depends on the value, the number of calculations is not known before the calculation. Therefore, transforming Equation 3 into a form without square root and division is important when using a DSP, and Equation 4 is the result of the operation.
  • I (x) ((c , 2 one c 3 2) x_ one c 3 2 d 3 2)
  • Equation 4 is a 12th-order equation, it is a rational polynomial and is differentiable in the entire range, so it can be approximated as a straight line sufficiently near the solution. It is only necessary that X satisfy the condition that it is always near the solution.
  • Equation 4 has an inflection point near the solution because it also has a solution of gl (x) + g 2 (x), so the position where the approximate value is found needs to be closer to the true solution than this inflection point. It is. Using x + pXd 3 ZF as the starting point of the approximate solution, using the notation in the figure as appropriate for that condition. An X at the vibrator position N and we do x N is determined by the following equation 5. From now on, the left side of Equation 4 is expressed as f (x). d 3 f (x N — t ten pd 3 / F)
  • UPF f (x N , + pd, / F), 5) A method of calculating the delay time of each transducer in real time using this recurrence formula will be described with reference to a flowchart.
  • a probe such as a linear probe or a convex probe, whose transmission plane is orthogonal to the line segment connecting the focal point and the center transducer of the aperture, is used.
  • the case of calculating the delay time will be described. Under this condition, since the focal point is in front of the aperture, the delay time of each transducer is symmetrical within the aperture. In other words, it is enough to calculate only one half.
  • the delay time is calculated one by one toward the outside, starting from the transducer at the center of the aperture.
  • the oscillator number N be 1 (the next oscillator with respect to the center oscillator).
  • x. 0, X!
  • For oscillator number N l! Is calculated by Equation 5 as follows.
  • e 3 arctan (x Zd 3 )
  • ⁇ 2 arcsin (c 2 sin ⁇ 3 , C 3 )
  • Te - x tens phi 2 + 2 2 tens dc 2 + x 3 2 tens d 3 _z c 3 then proceeds to step 24, the transducer number is incremented by one, moved to the vibrator the adjacent, step Step 25
  • x 2 is calculated by the recurrence recurrence formula (Equation 5), and the delay time ⁇ with respect to the oscillator of the oscillator number 2 is calculated based on the x 2 by using Formula 6. This operation is performed for all the transducers having a diameter.
  • step 25 If the determination in step 25 is "YES", it means that the delay time for all transducers has been calculated for this depth of focus, so the process proceeds to step 26, where the calculation result is output from the control circuit to the delay circuit. You. If the next depth of focus to be calculated is given to the ultrasonic beamformer, the same procedure is repeated for the all-diameter transducer. By repeating this for all the depths of focus, imaging is performed in which a shift due to refraction has been corrected.
  • the line connecting the focal point and the center transducer of the aperture such as the oblique type probe, the sector type probe, and the phased array type probe, is not orthogonal to the transmission plane.
  • the delay time for a type of probe, or the delay time for a type of probe in which the angle between the line connecting the focal point and the center transducer of the aperture and the transmitting surface changes The method for calculating the following is described.
  • the relationship between the focal position, the aperture, and the vibrator row depends on the magnitude of the oblique angle 0, as schematically shown in FIG. (1)
  • the perpendicular foot A dropped from the focal point F to the plane where the transducers are aligned enters the transmission aperture
  • step 33 X is obtained by the graduation formula (formula 5), and in step 34, the delay time ⁇ is obtained from X by formula 6.
  • step 34 the process moves to the next vibrator in step 35, and the processes in steps 33 to 35 are repeated until it is determined in step 36 that the vibrator has reached the end of the aperture. With this, the depth of focus and the focus direction are delayed by all transducers. Since the delay time has been calculated, the calculation result is output from the operation / control circuit to the delay circuit in step 38. If the next depth of focus and the direction of focus to be calculated are given to the ultrasonic beamformer, the same calculation is performed again for all the transducers within the aperture. By repeating this for all the depths of focus and the focal directions, imaging is performed in which deviation due to refraction is corrected.
  • the approximate value of the delay time can be obtained with sufficient accuracy (1/10 of the center frequency).
  • This method has a sufficient accuracy compared with the conventional approximation formula, and has no loop in the algorithm as compared with the method of obtaining the accuracy by using the iteration, so that the calculation speed is considerably faster.
  • a real-time high-speed delay time calculation algorithm including refraction by a fat layer as in the present invention has been realized. Furthermore, it is significantly advantageous in that errors do not accumulate as compared with the calculation of a general grading formula.
  • Equation 3 can be regarded as a straight line near the solution, there is no need to calculate the slope from the derivative, and we use x + p X d 3 / F, x + p as the two points that clearly sandwich the solution. , ⁇ ⁇ ⁇
  • the solution is obtained as the following equation as a dividing point. Therefore, the processing speed can be improved by using the recurrence equation (Equation 7) instead of the recurrence equation (Equation 5). fix ⁇ + p d./F)
  • FIG. 11 it is assumed that the sound emitted from the transducer N at the transducer position X N at an angle of 0 N into the lens layer reaches the focal point F.
  • this angle 0 N is obtained, the refraction angle in each layer is obtained, so that the sound path is determined and the delay time can be determined.
  • sound is first emitted at an appropriate angle, and the side of the focal point through which the sound passes is calculated. If the sound is increased 0 N if impassable farther than the focal, to reduce the 0 N if impassable the near focus.
  • the algorithm used in the ultrasonic beamformer according to the present embodiment is an improvement of the two-Juton-Raphson method, which can be applied to a system having no explicit function. That is, think as follows. First, assuming that an ultrasonic wave is emitted at an appropriate angle 0, the magnitude of the shift is ⁇ (0). Since this ⁇ ⁇ ⁇ ⁇ becomes 0, 0 is the solution to be found, and in order to find the derivative of ⁇ ⁇ ⁇ , the deviation ⁇ ( ⁇ + d0) when emitted at an angle of ⁇ + d0 is found. As shown in FIG.
  • Equation 8 As the shift amount, besides using ⁇ , a method using ⁇ or a method using the sum of squares of ⁇ and ⁇ can be considered. However, in this embodiment, a case where ⁇ X is used will be described.
  • the angle ⁇ ⁇ ⁇ ⁇ with respect to the next vibrator (vibrator of number 0) is used.
  • step 53 a deviation ⁇ ⁇ (0 + d0) between the path of the ultrasonic wave and the target focal position when the ultrasonic wave is radiated from the transducer of number 1 at an angle of 0 + d0 is determined.
  • the convex type probe unlike the linear type probe, it is calculated as follows.
  • the probe When observing an object with a convex probe, the probe is pressed against the object and observed, so the fat layer also bends along the probe. Assuming that the fat layer is concentric with the structure of the convex probe, the sound path is as shown in Figure 15.
  • the two methods for calculating the deviation ⁇ when leaving the element and passing near the focal point and the method for finding the deviation ⁇ when leaving the focal point and passing near the element are basically the same. This system is not symmetrical because the reflection conditions and the sound path cannot go to the left half of Fig. 15, and the former method has good solution stability. Since it is not, it is desirable to solve by the latter method.
  • Equation 10 which is a modification of Equation 9, and ⁇ 4 0 It can be obtained by using + arctan (-f '(x 3 )) instead.
  • step 54 a new 0 is obtained by extrapolation from equation 8. Subsequently, the process returns from step 55 to step 52, and uses the new 0 obtained in step 54 to shift ⁇ (0). Ask for.
  • step 53 a deviation ⁇ ( ⁇ + d0) between the path of the ultrasonic wave when the ultrasonic wave is emitted at the new angle ⁇ + d0 and the target focal position is obtained.
  • step 54 a new 0 is obtained by extrapolation from equation (8).
  • step 56 the delay time ⁇ is calculated from the obtained ⁇ by the following equation 11.
  • the oscillator number is increased by one, and the processing from step 42 is repeated. If it is determined in step 58 that the calculation has been completed for all the transducers within the working diameter, the process proceeds to step 59 and outputs the calculation result.
  • the flowchart shown in Fig. 13 is based on the calculation of the recurrence formula in the flowchart shown in Fig. 8 and the calculation of the deviation from the focal position and the extrapolation. This is equivalent to what is replaced by the calculation processing of a new 0 by the above and the judgment processing of the number of loops for performing the loop processing of this part. Also in the case of the second method, when the focal point is not in front of the transmission aperture, the transmission aperture is determined in advance to 0 at the end vibrator before imaging, and that value is stored. A plurality of zeros at the transducer at the end of the transmission aperture are obtained for different subcutaneous fat layer thicknesses.
  • This 0 is read and used to solve the recurrence equation as the first term of the recurrence equation.
  • the flow of this calculation is shown in the flowchart of FIG. The details are easily explained by the fact that the relationship between the calculation process in FIG. 8 and the calculation process in FIG. 9 and the relationship between the calculation process in FIG. 13 and the calculation process in FIG. 14 are exactly the same in the first method.
  • This method of pre-calculating the parameters for the transducer at the end of the aperture is a measure to prevent the DSP or MPU mounted on the current intermediate-level ultrasonic imaging device from having a margin, and a high-speed arithmetic processing unit can be mounted.
  • the first method all calculations for each fat layer thickness are performed in advance within a conceivable range, and all calculation results before imaging are stored in the memory 22 of the ultrasonic probe 20 as data. It is a way to keep it. Then, when the ultrasonic probe 20 is connected to the ultrasonic imaging apparatus main body 10, the contents are transferred to the memory 12 of the ultrasonic imaging apparatus main body 10. During imaging, a delay time is calculated by the ultrasonic beamformer 14 with respect to this parameter and the depth of focus given by the control unit 11, and the delay time is given to the transducer array 21, and transmitted to the subject 30. I do. The reception is also phased by this delay time, and the control unit 11 calculates a diagnostic image and outputs it to the image display unit 13.
  • the second method data calculated in advance is stored in a medium such as a CD-ROM and attached to the ultrasonic probe.
  • the method is to install data from the CD-ROM or the like into the memory 12 of the ultrasonic imaging apparatus body 10 before or when using the probe.
  • the operation after installation is exactly the same as the first method.
  • FIG. 16 is a schematic configuration diagram showing another example of the ultrasonic imaging apparatus according to the present invention.
  • the same functional portions as in FIG. 1 are denoted by the same reference numerals as in FIG. 1, and redundant description will be omitted.
  • the input of the thickness of the subcutaneous fat layer has been performed by the diagnostician. That is, as shown in FIG. 6, when the diagnostician specifies a start portion and an end portion of the fat layer on the diagnostic screen of the diagnostic apparatus, the distance is displayed on the screen. This value was entered by the diagnostician.
  • the ultrasonic imaging apparatus shown in FIG. 16 includes a fat layer thickness calculation unit 16 in the control unit 11 so that screen output related to the subcutaneous fat layer can be directly input to the ultrasonic beamformer 14. It is something that has been done.
  • the fat layer thickness calculation mode as shown in FIG. 6, when the diagnostician designates a portion A where the subcutaneous fat layer starts and a portion B where the subcutaneous fat layer starts by using a pointing device or the like, the fat layer thickness calculating section 16 calculates a distance between them. The calculation result is sent to beamformer 14.
  • FIG. 17 is a schematic configuration diagram showing another example of the ultrasonic imaging apparatus according to the present invention.
  • This ultrasonic imaging apparatus makes it possible to fine-tune the parameters relating to the fat layer, and to cope with a fat layer having a uniform thickness and non-constant thickness.
  • the ultrasonic imaging apparatus having the configuration shown in FIG. 1 was provided with an input device and a parameter calculation unit 17 so that the sound speed of the fat layer held by the ultrasonic beamformer could be varied.
  • a parameter corresponding to an input value is calculated in an input device and a parameter calculation unit 17, and the value is input to the ultrasonic beam former 14.
  • FIG. 19 shows an embodiment in which the examiner can select them.
  • the table having the thickness shown in Fig. 19 (b) and the table of sound speed are prepared in advance and installed in the device.
  • the examiner observes the image 170 displayed on the monitor 120, and obtains the image by referring to the depth scale 160.
  • the thickness of the layer is read from the force, and a similar thickness is selected with switch SW2 (140) or the like, and the sound speed can be similarly selected with switch SW1 (150).
  • switch SW1 and switch SW2 can be operated independently, so that the examiner can select values that are appropriate for each value.
  • a numerical value may be used, or a selection may be made between a fat layer and a muscle layer.
  • the fat layer 170 can be read as approximately 2 cm from the depth scale 160 shown on the screen of FIG. 19, 2 cm is selected by the layer thickness setting switch 140 provided on the console 130. Then, the screen Thickness is displayed at the top as 2cm. Next, the sound speed is set by operating the sound speed setting switch 150 in a predetermined direction. It is advisable to display those input values on the screen to confirm the selected values.
  • the indication may be a numerical value such as 1450m / s, but may be an expression such as hard muscle, muscle, normal, fat, or high fat.
  • the switch can be a rotary type or anything. Of course, a touch panel may be used.
  • FIG. 20 shows another embodiment of the present invention.
  • the sound speed of a living body having a layer structure is obtained, and the refraction correction processing of the above-described embodiment is performed based on the sound speed.
  • FIG. 20 is a block diagram showing an embodiment of a part for obtaining the sound speed in the ultrasonic imaging apparatus.
  • reference numeral 200 denotes a digitally controllable digital delay unit that delay-controls a plurality of ultrasonic signals received by the probe and outputs a reception beam signal, and has a circuit having a number of channels corresponding to the number of transducers used for reception. have.
  • Reference numeral 210 denotes a delay data error estimator that inputs a plurality of signals delayed by the digital delay unit 200 and estimates an error of the digital delay data subjected to delay control with respect to true delay data by calculation
  • 220 denotes a digital delay unit A digital delay control unit that controls the operation of each of the 200 channels.
  • the ultrasonic imaging apparatus of the present embodiment further includes a sound speed corresponding delay time recording unit 230 that stores delay times due to a plurality of medium sound speeds in advance, and a delay obtained by the delay error estimation unit 210.
  • a new delay time is calculated from the error, and the calculated delay time is compared with a value stored in the sound speed corresponding delay time recording unit 230, and a delay time comparing unit 240 which outputs data closest to the stored value, and a sound speed corresponding delay
  • the sound speed data recording unit 250 that records the medium sound speed based on the delay time data stored in the time recording unit 230, and the delay time from the delay time recording location that matches the output of the delay time comparison unit 240
  • a medium sound speed selection unit 260 for selecting the medium sound speed with reference to the recording unit;
  • the output line of the medium sound speed selection unit 260 is connected to the arithmetic and control circuit 41 shown in FIG.
  • the digital delay control unit 220 is connected so as to be controlled by the arithmetic and control circuit
  • the delay circuit is digitally controllable, a delay error estimator 210 is provided at the output thereof, and a digital delay controller 220 is added to the arithmetic and control circuit 41. Further, as described above, by newly providing the sound speed corresponding delay time recording unit 230, the delay time comparing unit 240, the sound speed data recording unit 250, and the medium sound speed selecting unit 260, the ultrasonic imaging apparatus of the present embodiment can realize the true medium sound speed. Refraction correction with a sound speed almost equal to that of the above becomes possible.
  • delay time data obtained by assuming the sound velocity in the living body to be, for example, the equal sound velocity of the average value of the living body is calculated and output from the control circuit 42 to the digital delay unit 2 and the ultrasonic wave is transmitted into the subject. .
  • the focal point of the transmission is set to an appropriate depth.
  • the delay time data D based on the same sound speed as the transmission wave is supplied from the calculation / control circuit 41 to the digital delay unit 200 via the digital delay unit 220, and the received signal is delay-controlled.
  • the correlation processing method disclosed in the above-mentioned Japanese Patent Application Laid-Open No. 8-317923 can be used.
  • D cl output from delay error estimation section 210 is input to delay time comparison section 240.
  • the delay time comparing unit 240 compares the data recorded in the sound speed corresponding delay time recording unit 230 with the input data D cl , selects the delay time data closest to D cl, and sends it to the medium sound speed selecting unit 260. Output.
  • the medium sound speed selection unit 260 selects the sound speed of the input data. This selection can be made by associating the storage address of the data of the sound speed corresponding delay time recording unit 230 with the sound speed data of the sound speed recording unit 250. Therefore, the sound speed corresponding delay time recording unit 230 and the sound speed recording unit 250 can be integrated into one. ,
  • the sound speed data selected by the medium sound speed selection unit 260 is fed back to the digital delay control unit via the arithmetic and control circuit 41.
  • This feedback circuit Is for repetition of the above operation in order to obtain an accurate sound speed. If necessary, the arithmetic and control circuit 41 issues a command to repeatedly execute the above operation. With the above operation, the sound velocity in the living body can be measured as an estimated value.
  • the data D cl output by the delay error estimation unit 210 and the data recorded in the sound speed corresponding delay time recording unit 230 are two-dimensional distribution data. It is useful to be able to use the curve fitting method because the amount of information to be handled can be reduced. Furthermore, it is useful to obtain the difference of the delay time distribution and to fit a first-order straight line to the difference delay time sequence because the amount of information to be handled is further reduced.
  • FIG. 21 is an ultrasonic tomographic image of the subject displayed on the monitor of the ultrasonic imaging apparatus.
  • 170 is a fat layer
  • 190 is a tissue region deeper than the fat layer 170.
  • the sound speed of the fat layer 170 is determined.
  • the focal position is set using the above-mentioned cursor of the caliber at an appropriate depth position inside the fat layer 170, for example, a measurement position near the boundary between the fat layer 170 and the tissue region, and the value is measured.
  • a signal for giving the focal position is supplied to the arithmetic / control circuit 41.
  • the arithmetic control circuit 41 is provided with a gating function.
  • the gating function may be linked to the caliber function.
  • the reflected signal from the fat layer 170 is taken into the digital delay unit 200 by this gating function, and the sound speed is obtained in accordance with the operation description of the configuration shown in FIG.
  • the desired sound velocity may be at a certain point, but it is desirable to determine the sound velocity at multiple points and then determine the average sound velocity from them.
  • the sound velocity in the tissue region deeper than the fat layer 170 is determined.
  • the focal position for measurement is set to the position 190 shown in Fig. 21.
  • the above-mentioned caliber function is also used for this setting.
  • a signal for giving a focal position is supplied to the arithmetic and control circuit 41, and thereafter, an ultrasonic pulse is transmitted from the probe and a reflected signal thereof is received.
  • a reflection signal from the focal position is captured by the above-mentioned gating function.
  • the speed of sound is obtained from the reflected signal.
  • the sound speed obtained here indicates the average sound speed between the probe surface and the measurement point.
  • the sound velocity in the tissue region can be obtained from the above two measurements.
  • the average speed of sound was measured first fat layer, then the obtained fat layer and the average speed of sound including both tissue region put c a and code.
  • the sound velocity c 2 in the tissue region can be obtained by Expression 13.
  • This embodiment considers the correction of the distance measurement function of the carrier used for the sound velocity measurement. That is, the distance measurement function of the caliper incorporated in the ultrasonic imaging apparatus is performed by an operation based on the speed of sound initially set in the apparatus. Therefore, it is desirable that the distance measured by the caliper of l f , 1 ⁇ used in Equation 13 be used after being corrected to the value based on the actual sound speed.
  • the present embodiment corresponds to this.
  • the ultrasonic imaging apparatus is driven to acquire an ultrasonic tomographic image of a section including the region of interest (R0I), and the tomographic image is displayed on a monitor.
  • the refraction correction execution switch 310 arranged on the operation panel of the ultrasonic imaging apparatus shown in FIG. 19 is turned on.
  • This refraction correction execution switch 310 performs refraction correction when turned on, and refraction when turned off. This enables normal imaging without correction.
  • a caliper 300 consisting of two points of force is displayed on the screen of the monitor 120.
  • an input operating device such as a trackball or a joystick
  • one cursor of the caliper 300 is moved to the body surface of the subject, and the other cursor is moved to the end of the fat layer to fix the input information.
  • Operating the key (Enter key) 320 specifies the fat layer to be measured.
  • the force sol positioned at the end of the fat layer is moved into the region of interest at a point deeper than that, and the sound velocity in the tissue is operated by operating the key 320. Identify measurement points for measurement.
  • the data of the measurement points input by the above two-step operation are read into the control unit 11, and the calculation of the above-described sound velocity measurement method is executed.
  • the device setting sound speed is v. , Its sound velocity V. Caliber Depth at x. , Assuming the measured average sound velocity of the fat layer c fn , the true thickness l ft of the fat layer can be obtained from Eq.
  • Various methods other than those described above can be used to determine the sound speed and distance of a living body having a layered structure.
  • the refraction correction method described above is programmed and incorporated into the device as an automatic sequence, and repeatedly executed every few frames during imaging to update the set values, if there are multiple inspection sites Move the probe as shown The refraction correction is performed automatically, so the image is always good.
  • a normal image is displayed on the left of the screen, and a right image is displayed on the left side of the screen using a method in which a displayed image is marked to distinguish between a normal image and a refraction corrected image.
  • a refraction-corrected image is displayed. Even if the image at the point of interest is good, the image may be distorted in other regions in the case of a complex biological structure. Therefore, simultaneous display makes it possible to compare images with and without refraction correction in real time, which is effective for the user. It is also useful to specify the diagnosis site as R0I, acquire an image with refractive correction applied only to that site, fit it into a normal image, and display it.
  • the effect of refraction due to the lens layer and the fat layer (or the muscle layer, etc.) can be reduced, so that the image quality of an ultrasonic image can be improved.
  • the sound velocity of each layer of the subject having a layered structure can be measured, and the influence of the refraction of the ultrasonic wave in each layer can be reduced by taking the value into account, thereby further improving the image quality. can do.

Landscapes

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

Abstract

Cette invention a trait à un dispositif d'imagerie ultrasonore possédant des moyens de production de données de retard de correction de réfraction permettant d'établir un temps de retard aux fins de la focalisation d'une onde émise ou reçue, en fonction de l'effet de réfraction d'une onde ultrasonore par un milieu de propagation ultrasonore entre des dispositifs vibreurs en réseau et le point focal prédéterminé. Ces moyens font parvenir le temps de retard à un dispositif de commande de retard et ce, afin de produire une image ultrasonore d'excellente qualité tout en évitant tout décalage focal de l'émission/réception de l'onde ultrasonore survenant au niveau de l'interface entre la couche d'objectif d'une sonde ultrasonore et la couche adipeuse du sujet ainsi qu'au niveau de l'interface entre la couche adipeuse et un tissu biologique présent sous cette couche adipeuse.
PCT/JP2000/007114 1999-10-15 2000-10-13 Dispositif d'imagerie ultrasonore WO2001026555A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001529349A JP4711583B2 (ja) 1999-10-15 2000-10-13 超音波撮像装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP11/294498 1999-10-15
JP29449899 1999-10-15

Publications (1)

Publication Number Publication Date
WO2001026555A1 true WO2001026555A1 (fr) 2001-04-19

Family

ID=17808556

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2000/007114 WO2001026555A1 (fr) 1999-10-15 2000-10-13 Dispositif d'imagerie ultrasonore

Country Status (2)

Country Link
JP (1) JP4711583B2 (fr)
WO (1) WO2001026555A1 (fr)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009061086A (ja) * 2007-09-06 2009-03-26 Fujifilm Corp 超音波診断装置、並びに、画像処理方法及びプログラム
JP2010167258A (ja) * 2008-12-25 2010-08-05 Canon Inc 生体情報取得装置
JP2010167257A (ja) * 2008-12-25 2010-08-05 Canon Inc 生体情報取得装置
CN101900808A (zh) * 2010-06-30 2010-12-01 深圳市蓝韵实业有限公司 一种多通道扫描延时参数计算方法及装置
EP2347704A1 (fr) * 2010-01-22 2011-07-27 Canon Kabushiki Kaisha Système de mesure d'ondes acoustiques et procédé de génération de données d'image
WO2011105593A1 (fr) 2010-02-23 2011-09-01 Canon Kabushiki Kaisha Appareil d'imagerie échographique et procédé de contrôle de délai
WO2011122044A1 (fr) * 2010-03-31 2011-10-06 富士フイルム株式会社 Procede et dispositif d'acquisition de la vitesse d'un son environnemental
CN103913747A (zh) * 2012-12-31 2014-07-09 深圳市蓝韵实业有限公司 多通道扫描延时计算方法及装置
CN107205729A (zh) * 2015-02-18 2017-09-26 株式会社日立制作所 超声波成像装置以及超声波信号的处理方法
JP2019033822A (ja) * 2017-08-10 2019-03-07 コニカミノルタ株式会社 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
JP2021536276A (ja) * 2018-09-05 2021-12-27 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 超音波画像による脂肪層の識別

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11779312B2 (en) * 2019-06-28 2023-10-10 Siemens Medical Solutions Usa, Inc. Ultrasound medical imaging with optimized speed of sound based on fat fraction
US11864949B2 (en) * 2020-07-30 2024-01-09 Fujifilm Sonosite, Inc. Beamforming in ultrasound imaging systems to correct for refraction of ultrasound by transducer elements

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4835689A (en) * 1987-09-21 1989-05-30 General Electric Company Adaptive coherent energy beam formation using phase conjugation
JPH0595946A (ja) * 1991-10-09 1993-04-20 Yokogawa Medical Syst Ltd 超音波の音速測定法及び音速測定手段を備えた超音波診断装置
JPH08308832A (ja) * 1995-05-17 1996-11-26 Toshiba Corp 超音波診断装置
JPH09187456A (ja) * 1996-01-10 1997-07-22 Fujitsu Ltd 超音波診断装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4989143A (en) * 1987-12-11 1991-01-29 General Electric Company Adaptive coherent energy beam formation using iterative phase conjugation
JP2619446B2 (ja) * 1987-12-21 1997-06-11 株式会社日立製作所 超音波診断装置
JPH01244737A (ja) * 1988-03-26 1989-09-29 Toshiba Corp アニュラアレイ超音波診断装置
JP3474278B2 (ja) * 1994-03-16 2003-12-08 フクダ電子株式会社 超音波診断装置
JPH11104130A (ja) * 1997-10-07 1999-04-20 Hitachi Medical Corp 超音波診断装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4835689A (en) * 1987-09-21 1989-05-30 General Electric Company Adaptive coherent energy beam formation using phase conjugation
JPH0595946A (ja) * 1991-10-09 1993-04-20 Yokogawa Medical Syst Ltd 超音波の音速測定法及び音速測定手段を備えた超音波診断装置
JPH08308832A (ja) * 1995-05-17 1996-11-26 Toshiba Corp 超音波診断装置
JPH09187456A (ja) * 1996-01-10 1997-07-22 Fujitsu Ltd 超音波診断装置

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009061086A (ja) * 2007-09-06 2009-03-26 Fujifilm Corp 超音波診断装置、並びに、画像処理方法及びプログラム
JP2010167258A (ja) * 2008-12-25 2010-08-05 Canon Inc 生体情報取得装置
JP2010167257A (ja) * 2008-12-25 2010-08-05 Canon Inc 生体情報取得装置
US8876717B2 (en) 2008-12-25 2014-11-04 Canon Kabushiki Kaisha Biological information acquisition apparatus
EP2381837A1 (fr) * 2008-12-25 2011-11-02 Canon Kabushiki Kaisha Appareil d'acquisition de données biologiques et méthode d'acquisition de données biologiques
EP2381836A1 (fr) * 2008-12-25 2011-11-02 Canon Kabushiki Kaisha Appareil d'acquisition de données biologiques
EP2347704A1 (fr) * 2010-01-22 2011-07-27 Canon Kabushiki Kaisha Système de mesure d'ondes acoustiques et procédé de génération de données d'image
US8611179B2 (en) 2010-01-22 2013-12-17 Canon Kabushiki Kaisha Acoustic wave measuring system and method of generating image data thereof
US9165552B2 (en) 2010-02-23 2015-10-20 Canon Kabushiki Kaisha Ultrasonic imaging apparatus and method of controlling delay
WO2011105593A1 (fr) 2010-02-23 2011-09-01 Canon Kabushiki Kaisha Appareil d'imagerie échographique et procédé de contrôle de délai
JP2011172611A (ja) * 2010-02-23 2011-09-08 Canon Inc 超音波画像装置および遅延制御方法
US20120281902A1 (en) * 2010-02-23 2012-11-08 Canon Kabushiki Kaisha Ultrasonic imaging apparatus and method of controlling delay
CN102811667A (zh) * 2010-03-31 2012-12-05 富士胶片株式会社 环境音速取得方法和装置
JP2011212066A (ja) * 2010-03-31 2011-10-27 Fujifilm Corp 環境音速取得方法および装置
WO2011122044A1 (fr) * 2010-03-31 2011-10-06 富士フイルム株式会社 Procede et dispositif d'acquisition de la vitesse d'un son environnemental
US9291601B2 (en) 2010-03-31 2016-03-22 Fujifilm Corporation Ambient sound velocity obtaining method and apparatus
CN101900808A (zh) * 2010-06-30 2010-12-01 深圳市蓝韵实业有限公司 一种多通道扫描延时参数计算方法及装置
CN103913747A (zh) * 2012-12-31 2014-07-09 深圳市蓝韵实业有限公司 多通道扫描延时计算方法及装置
CN107205729A (zh) * 2015-02-18 2017-09-26 株式会社日立制作所 超声波成像装置以及超声波信号的处理方法
EP3260047A4 (fr) * 2015-02-18 2018-09-26 Hitachi, Ltd. Dispositif de capture d'image à ultrasons et procédé de traitement de signal ultrasonore
CN107205729B (zh) * 2015-02-18 2019-12-31 株式会社日立制作所 超声波成像装置以及超声波信号的处理方法
US11076827B2 (en) 2015-02-18 2021-08-03 Hitachi, Ltd. Ultrasound image capturing device and method of processing ultrasound signal
JP2019033822A (ja) * 2017-08-10 2019-03-07 コニカミノルタ株式会社 超音波信号処理装置、超音波診断装置、および、超音波信号処理方法
JP2021536276A (ja) * 2018-09-05 2021-12-27 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 超音波画像による脂肪層の識別
JP7358457B2 (ja) 2018-09-05 2023-10-10 コーニンクレッカ フィリップス エヌ ヴェ 超音波画像による脂肪層の識別

Also Published As

Publication number Publication date
JP4711583B2 (ja) 2011-06-29

Similar Documents

Publication Publication Date Title
US6705994B2 (en) Tissue inhomogeneity correction in ultrasound imaging
CN100515346C (zh) 超声波诊断装置、超声波诊断方法及其图像处理程序
EP2816958B1 (fr) Détermination de la raideur d'un matériau à l'aide d'ultrasons à multiples ouvertures
JP5435751B2 (ja) 超音波診断装置、超音波送受信方法、および超音波送受信プログラム
CN105407809A (zh) 超声波诊断装置以及弹性评价方法
EP4018936B1 (fr) Appareil de diagnostic ultrasonore et procédé de commande d'un appareil de diagnostic ultrasonore
CN113679425A (zh) 一种超声弹性检测方法及其***
WO2001026555A1 (fr) Dispositif d'imagerie ultrasonore
JP2012217611A (ja) 超音波診断装置および超音波画像生成方法
US20120310093A1 (en) Ultrasound image producing method and ultrasound image diagnostic apparatus
JP6871478B2 (ja) 超音波診断装置および超音波診断装置の制御方法
JP7476376B2 (ja) 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ
US20120004551A1 (en) Ultrasonic diagnosis apparatus and ultrasonic diagnosis method
JP2023155494A (ja) 超音波診断装置
JP5588924B2 (ja) 超音波診断装置
JP2008073423A (ja) 超音波診断装置、診断パラメータ計測装置及び診断パラメータ計測方法
JP6746895B2 (ja) 超音波診断装置、及び超音波信号処理方法
JP6081744B2 (ja) 超音波診断装置
JP7434529B2 (ja) 超音波診断装置、超音波診断装置の制御方法および超音波診断装置用プロセッサ
JP6876869B2 (ja) 超音波診断装置および超音波診断装置の制御方法
JP4908915B2 (ja) 超音波診断装置
US20180168542A1 (en) Ultrasound diagnosis apparatus and method of controlling the same
CN104684486A (zh) 超声波诊断装置、超声波图像生成方法及程序
JP2012239545A (ja) 超音波診断装置
JP2013102960A (ja) 超音波診断装置及び方法

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): JP US

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
122 Ep: pct application non-entry in european phase