WO2013121842A1 - 超音波撮像装置 - Google Patents
超音波撮像装置 Download PDFInfo
- Publication number
- WO2013121842A1 WO2013121842A1 PCT/JP2013/051305 JP2013051305W WO2013121842A1 WO 2013121842 A1 WO2013121842 A1 WO 2013121842A1 JP 2013051305 W JP2013051305 W JP 2013051305W WO 2013121842 A1 WO2013121842 A1 WO 2013121842A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- unit
- phasing
- reception
- output
- imaging apparatus
- Prior art date
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 98
- 230000003044 adaptive effect Effects 0.000 claims abstract description 129
- 238000012545 processing Methods 0.000 claims abstract description 81
- 230000003111 delayed effect Effects 0.000 claims abstract description 60
- 239000002131 composite material Substances 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims description 211
- 238000004364 calculation method Methods 0.000 claims description 148
- 230000002093 peripheral effect Effects 0.000 claims description 144
- 238000009826 distribution Methods 0.000 claims description 73
- 230000015572 biosynthetic process Effects 0.000 claims description 65
- 238000003786 synthesis reaction Methods 0.000 claims description 65
- 230000002194 synthesizing effect Effects 0.000 claims description 20
- 239000000203 mixture Substances 0.000 claims description 11
- 238000012285 ultrasound imaging Methods 0.000 claims description 6
- 230000001934 delay Effects 0.000 claims description 4
- 230000006978 adaptation Effects 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 72
- 230000015556 catabolic process Effects 0.000 abstract description 4
- 238000006731 degradation reaction Methods 0.000 abstract description 4
- 239000000523 sample Substances 0.000 description 56
- 230000005540 biological transmission Effects 0.000 description 42
- 239000011159 matrix material Substances 0.000 description 41
- 230000006870 function Effects 0.000 description 38
- 238000010586 diagram Methods 0.000 description 21
- 230000000694 effects Effects 0.000 description 18
- 230000000875 corresponding effect Effects 0.000 description 14
- 230000008569 process Effects 0.000 description 12
- 230000033001 locomotion Effects 0.000 description 11
- 238000002604 ultrasonography Methods 0.000 description 8
- 238000012935 Averaging Methods 0.000 description 7
- 238000000926 separation method Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 230000002596 correlated effect Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 239000011295 pitch Substances 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005094 computer simulation Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000001727 in vivo Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241001272567 Hominoidea Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
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/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
-
- 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/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
-
- 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
-
- 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
-
- 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/465—Displaying means of special interest adapted to display user selection data, e.g. icons or menus
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8993—Three dimensional imaging systems
Definitions
- the present invention relates to an ultrasound imaging technique for capturing an image in a subject using ultrasound.
- the ultrasound imaging technique is a technique for non-invasively imaging the inside of a subject such as a human body using ultrasound (a sound wave not intended to be heard, generally a sound wave having a high frequency of 20 kHz or higher). It is.
- a medical ultrasonic imaging apparatus will be briefly described.
- the ultrasonic probe transmits an ultrasonic wave toward a patient's body and receives an echo signal reflected from the patient's body.
- the received signal is subjected to signal processing in one or both of the ultrasonic probe and the ultrasonic imaging apparatus main body, and then delivered to the image display unit to display an ultrasonic image.
- a transmission beam former is generated in a transmission beam former in the main body of the ultrasonic imaging apparatus, and is transmitted to an ultrasonic probe after passing through a transmission / reception separation circuit (T / R).
- the ultrasonic probe transmits ultrasonic waves.
- the ultrasonic probe receives an echo signal from the body and then transmits the signal to the imaging apparatus body.
- the received signal passes through the transmission / reception separation circuit again, is subjected to phasing processing in the reception beamformer, and then transmitted to the image processing unit.
- various image processes such as various filters and a scan converter are executed.
- an ultrasonic image is displayed on the image display unit.
- a general ultrasonic diagnostic apparatus is configured by three techniques of transmission beam forming, reception beam forming, and back-end image processing.
- the beamformer at the time of transmission and reception performs signal processing at the RF (high frequency) level
- the beamformer algorithm and the mounting architecture determine the basic image quality of the ultrasonic image. For this reason, the beamformer is the backbone of the apparatus.
- the receive beamformer has a delay time in which the amount of delay is distributed in a concave shape according to the relationship between the focal position and the position of the element for each received signal (received data) of a plurality of elements constituting the ultrasonic probe. Given, after focusing on one point in the virtual space, the received signal data is added. This method is called phasing by the delay addition method. In this delay addition method, the received data received by a plurality of elements of the ultrasonic diagnostic apparatus is multiplied by a fixed weight vector stored in the diagnostic apparatus, and weighted before addition. The same applies to the transmission beamformer as well as the reception beamformer.
- Non-Patent Documents 1 to 6 These adaptive methods are realized by adaptively changing the complex component of the weight vector for delay addition based on the correlation matrix of the received data. That is, conventionally, the weight vector has been a fixed value, but in the adaptive method, a weight vector is obtained by calculation using the received signal for each sampling point in the time direction of the received signal, and this is multiplied by the received signal.
- Patent Documents 1 to 3 relate to delay time adjustment using a plurality of beams
- Patent Document 4 relates to sum of delay outputs
- Documents 5 are known respectively.
- the directivity of a beam is sharp in an ultrasonic diagnostic apparatus using an adaptive beamformer.
- the transmitted ultrasonic signal is reflected at the focal point and scattered around the focal point, so that information on the object at the focal point is also present around the focal point.
- an ultrasonic signal existing around the focal point cannot be acquired, and image quality deteriorates.
- wavefront distortion due to non-uniform sound speed or the like occurs. Due to the wavefront distortion, reflected waves from the medium around the focus are mixed in the received sound wave. The reflected wave from the surrounding medium is a noise signal correlated with information from the focal point, and there is also a problem that the image quality is deteriorated due to the correlated noise.
- An object of the present invention is to provide an ultrasonic imaging apparatus that can compensate for image quality degradation caused by inhomogeneity of a subject medium.
- the ultrasonic imaging apparatus of the present invention delays a plurality of elements that receive an ultrasonic signal from a subject and a reception signal received by the plurality of elements according to a position of a predetermined reception focus, and outputs a delayed reception signal.
- a delay unit to be generated, a reception focal point and a peripheral region of the reception focal point, a peripheral information calculation unit that acquires information about a plurality of points from the received signal after delay, and information acquired for each of the plurality of points are combined
- a peripheral information synthesis unit that generates a phasing output using the synthesized information
- an image processing unit that generates image data using the phasing output.
- ultrasonic signals scattered around the focal point can be acquired, it is possible to reduce signal loss, high robustness against correlation noise caused by wavefront distortion, and a large processing load. It is possible to compensate for image quality degradation without increasing it.
- (A) The perspective view which shows schematic structure of the ultrasonic imaging device of this embodiment
- (A)-(f) Explanatory drawing which shows the effect of embodiment of this invention.
- (A)-(c) Explanatory drawing which shows the steering vector of the receiving beam former in 4th Embodiment.
- movement and effect of the signal processing in 8th Embodiment The block diagram which shows the structure of the temporary storage part in the receiving beam former of 9th Embodiment.
- movement and effect of signal processing in 9th Embodiment (A)-(g) Explanatory drawing which shows the steering vector setting operation
- (d) Explanatory drawing which shows the luminance distribution of the ultrasonic image (B mode image) obtained by the steering angle set in 11th Embodiment.
- FIG. 27 is a graph showing one point scatterer and its surrounding luminance distribution for the ultrasonic image of FIG. 26.
- the ultrasonic imaging apparatus delays a plurality of elements that receive an ultrasonic signal from a subject and a reception signal received by the plurality of elements according to the position of a predetermined reception focus.
- a delay unit that generates a reception signal, a peripheral information calculation unit, a peripheral information synthesis unit, and an image processing unit.
- the peripheral information calculation unit obtains information on a plurality of points from the reception signal after delay from the reception focal point and the peripheral region of the reception focal point.
- the peripheral information combining unit combines the information acquired for each of the plurality of points, and generates a phasing output using the combined information.
- the image processing unit generates image data using the phasing output. This will be specifically described below.
- FIG. 1A is a perspective view of the apparatus
- FIG. 1B is a block diagram showing a schematic configuration inside.
- the ultrasonic imaging apparatus includes an ultrasonic probe 101, an apparatus main body 102, an image display unit 103, and a console 110.
- a transmission beam former 104 As shown in FIG. 1B, a transmission beam former 104, a transmission / reception separation circuit (T / R) 107, a reception beam former 108, an image processing unit 109, and operations thereof are controlled.
- a control unit 111 is arranged.
- the ultrasonic probe 101 includes a plurality of elements (ultrasonic transducers) 106 arranged in an array.
- a transmission beam signal is generated in the transmission beam former 104 and passed to the ultrasonic probe 101 through the transmission / reception separation circuit 107.
- the ultrasonic probe 101 transmits ultrasonic waves from a plurality of elements toward the inside of the subject 100.
- the echo signal reflected in the body is received by the ultrasonic probe 101.
- the received signal is subjected to phasing calculation processing and the like in the reception beam former 108 again through the transmission / reception separation circuit 107.
- the received signal after the phasing calculation process is transferred to the image processing unit 109, and various image processes such as various filters and a scan converter are executed to generate an ultrasonic image.
- the ultrasonic image is transferred to the image display unit 103 and displayed.
- FIG. 2 is a block diagram showing a configuration of the reception beamformer 108.
- the reception beamformer 108 includes a delay circuit 204, a peripheral information calculation unit 205, a peripheral information combining unit 206, and a bypass line 207 from the output of the delay circuit 204 to the peripheral information combining unit.
- Each unit of the reception beamformer 108 can be configured by an independent circuit, and the operation of each unit is realized by a memory that stores a program in advance and a CPU that reads and executes the program. It is also possible.
- an active channel setting unit sets an active channel in a part of the finite aperture of the ultrasound probe 101 under the control of the control unit 111, and receives beamforming processing. I do.
- a part of the elements 106 are defined as the active channel 201, and the reception signal of the active channel 201 is represented by Using this, one piece of image data (raster: phasing output y (n)) in the ultrasonic propagation direction is generated. As shown in FIG.
- the active channel 202, the active channel 201, and the active channel 203 are sequentially configured while shifting the position of the element 106 little by little, and rasters are generated for the respective active channels 202, 201, and 203, and the results are arranged.
- the image becomes an ultrasonic image.
- the delay circuit 204 is arranged as a preceding stage of the peripheral information calculation unit 205, and gives a delay time to each received signal (received data) of the plurality of elements 106 constituting the active channel 201 in accordance with the position of the element. This block performs processing for focusing on a certain point in space.
- the delay time provided for each element 106 is prepared in advance as a delay time set corresponding to each target position in the imaging target space.
- the delay circuit 204 selects a delay time set according to the set focal position, and gives different delay times to the reception signals of the plurality of elements. As a result, it is possible to perform focus processing that covers the entire space to be imaged.
- the peripheral information calculation unit 205 performs an operation of collecting information scattered around the reception focus due to reflection from the reception focus. For example, the peripheral information calculation unit 205 obtains adaptive weights for two or more points of the reception focus and its peripheral part (peripheral region) by performing adaptive signal processing. Thereby, the peripheral information of the reception focus can be collected.
- the peripheral portion is a region where signals are scattered due to reflection of signals at the reception focal point.
- the receiving element, transmission frequency, reception frequency, ultrasonic pulse waveform, element array size / pitch, and dots The size varies depending on conditions such as the spread function, but roughly, it has a long axis that extends within ⁇ 5 ° in the active channel direction with the axis connecting the reception focus and the center of the active channel as the central axis, and is active
- the ultrasonic wave propagation direction orthogonal to the channel direction is a region of an approximately elliptical shape (elliptical sphere in the case of three dimensions) having a minor axis that is as long as or slightly shorter than the active channel direction.
- the focal point is located at the intersection of the short axis and the long axis.
- peripheral portion is a noise in which reflected waves from the medium around the focus are mixed in the received sound wave due to wavefront distortion, and the reflected waves from the surrounding medium are correlated with information from the focus. It is also an area that appears as a signal in the received signal.
- the peripheral information synthesis unit 206 synthesizes (compounds) the peripheral information collected by the peripheral information calculation unit 205 to obtain a final phasing output y (n). For example, a composite weight obtained by combining two or more adaptive weights obtained by the peripheral information calculation unit 205 is obtained, and the signal of each element 106 subjected to delay processing by the delay circuit 204 is weighted and phased by the composite weight. Peripheral information can be synthesized.
- the final phasing output y (n) is sequentially transmitted to the image processing unit 109 for each raster corresponding to the active channels 202, 201, and 203, and an ultrasonic image is generated based on the synthesized phasing output. Generated.
- the peripheral information calculation unit 205 collects information from the periphery of the reception focal point, so that it is reflected by the focal point and scattered while using an adaptive beamformer with high directivity. Information can also be collected. Furthermore, by using the collected peripheral information, the asymmetry of the correlation signal can be used, and the correlation noise caused by wavefront distortion can be canceled (decorrelated). Therefore, even if the subject medium is inhomogeneous, image quality deterioration can be prevented, and an ultrasonic imaging apparatus with high robustness can be obtained that is hardly affected by correlation noise due to wavefront distortion.
- the peripheral information calculation unit 205 obtains an adaptive weight as information by performing adaptive beamforming.
- the peripheral information calculation unit 205 uses a steering vector, which is a direction vector connecting a predetermined position on the surface of the array composed of the plurality of elements 106 shown in the first embodiment and a plurality of points, for a plurality of points.
- the adaptive weight is obtained.
- the peripheral information calculation unit 205 includes a matrix calculation unit that generates a covariance matrix using the delayed received signal, and a weight vector calculation unit that obtains adaptive weight vectors for the plurality of points from the covariance matrix and the steering vector. It is set as the structure containing these.
- the peripheral information combining unit 206 adds a weighting unit that generates adaptive weights by adding adaptive weights for a plurality of points obtained by the peripheral information calculating unit 205, and multiplies the combined weight and the delayed received signal to output a phasing output.
- a fixed apodization multiplication unit that multiplies the adaptive weights for a plurality of points obtained by the peripheral information calculation unit by predetermined fixed weights, respectively. Is also possible.
- the inner product calculation unit multiplies the delayed received signal by the combined weight and then adds the delayed received signal to generate a phasing output.
- the ultrasonic imaging apparatus according to the second embodiment of the present invention will be specifically described with reference to FIG. Note that description of the same configuration as in the first embodiment is omitted.
- the peripheral information calculation unit 205 collects the peripheral information of the reception focus by obtaining adaptive weights for two or more points of the reception focus and the peripheral portion. Specifically, when calculating the adaptive weight, a plurality of steering vectors a are used, and the point where the direction from the active channel element of the ultrasonic probe deviates from the reception focus (that is, the direction of the steering vector is the focus direction). The adaptive weights are also obtained for points different from the above.
- the peripheral information combining unit 206 combines two or more adaptive weights obtained by the peripheral information calculation unit 205, and performs phasing processing of the received signal using the obtained combined weights.
- the steering vector is a direction vector that connects the center position of the active channel and the point. The steering vector will be described in detail with reference to equation (4) described later. Furthermore, it demonstrates concretely.
- FIG. 3 is a block diagram showing the configuration of the receive beamformer 108 of the second embodiment.
- an embodiment using an adaptive beamformer in which an adaptive signal processing technique is applied to the beamformer. Indicates. It should be noted that images other than the reception beamformer 108, that is, the ultrasonic probe 101, the transmission beamformer 104, the transmission / reception circuit 107, and the reception beamformer 108, which are located in the previous stage of the reception beamformer 108, are located in the subsequent stage. Since the processing unit 109 and the image display unit 103 are the same as those in the first embodiment, description and illustration are omitted.
- the peripheral information calculation unit 205 in the reception beamformer 108 includes a matrix calculation unit 300 and an adaptive weight calculation unit 301.
- the adaptive weight calculation unit 301 includes a plurality of weight vector calculation units 3021, 3022, and 3023, and calculates weight vectors for the reception focus and its peripheral points, respectively.
- weight vectors can be obtained for three points.
- the number of weight vector calculation units is not limited to three, and may be any number as long as it is plural.
- the matrix calculation unit 300 receives an input (delayed reception data x (n)) from the delay circuit 204 and generates a spatial covariance matrix R (n).
- the spatial covariance matrix is input to the adaptive weight calculation unit 301.
- three predetermined points with different steering vector directions in the reception focus and the peripheral portions thereof. The weight vectors w 1 (n), w 2 (n), and w 3 (n) are calculated and output to the peripheral information combining unit 206.
- the steering vector will be described in detail later.
- the peripheral information synthesis unit 206 includes a weight synthesis unit 306 and an inner product calculation unit 307.
- the weight combining unit 306 combines a plurality of weight vectors w 1 (n), w 2 (n), and w 3 (n) to generate a combined weight w sum (n).
- the inner product calculation unit 307 weights and adds the delayed received data x (n) input from the bypass line 207 with the combined weight w sum (n) to obtain a phasing output y (n). Thereby, information from a plurality of points with different steering vector directions in the same raster can be obtained as one phasing output.
- the inner product calculating unit 307 multiplies the combined weight w sum (n) by the delayed received data x (n), and an adder 3072 that calculates the sum of the delayed delayed data. And is configured.
- a fixed apodization multiplication unit 305 can be added to the peripheral information synthesis unit 206 as necessary.
- the fixed apodization multiplication unit 305 multiplies the weight vectors w 1 (n), w 2 (n), and w 3 (n) by predetermined fixed weights to weight them.
- the phasing output y (n) obtained by the peripheral information synthesis unit 206 is transmitted to the image processing unit 109.
- the image processing unit 109 generates an ultrasonic image based on the phasing result y (n).
- each unit when processing the reception data of each element 106 of one active channel 201 corresponding to one transmission / reception and generating phasing data of one raster will be described.
- FIG. 4 is a diagram showing the operation of the delay circuit 204 and the like.
- K pieces of received data are respectively received by the K elements 106 constituting the active channel 201.
- the K pieces of received data are also referred to as K channel received data u (n).
- the reception data u (n) is input to the delay circuit 204 through the transmission / reception separation circuit 107.
- the delay circuit 204 receives the received signals (received data 404) u 1 (n), u 2 (n),..., U K (n) of the K elements 106 constituting the active channel 201 in the space.
- n is a certain time (snapshot time) in the time direction of the ultrasonic reception signal (the depth direction in the ultrasonic wave propagation direction). For a total number N of samples of one raster, 0 ⁇ n ⁇ (N + 1).
- the K delayed received data at a snapshot time n can be expressed as a vector x (n) on the left side of the following equation (1). it can.
- the delayed received data x (n) is input to the adaptive processing unit, and is composed of the weight value w K (n) for each of K channels at the snapshot time n.
- Vector w (n) [w 1 (n), w 2 (n), ..., w K (n)] T is generated.
- an inner product operation is performed on the adaptive weight vector w (n) and the delayed received data x (n) by the inner product operation unit, thereby obtaining an adaptive beamformer output y (n) at a certain snapshot time n.
- an adaptive weight vector is also obtained for the peripheral portion of the focal position using a plurality of steering vectors a.
- the delayed received data vector x (n) is input to the peripheral information calculation unit 205 as shown in FIG.
- the matrix calculation unit 300 obtains a spatial covariance matrix R (n) by the following equation (2).
- x (n) may be used as a real signal, or may be used after being converted to complex data by performing Hilbert transform or baseband modulation. .
- x (n) is converted into complex data ⁇ (n) and a spatial covariance matrix R (n) is obtained as in equation (2).
- R (n) in equation (2) is an ensemble average of the product of the complex vector ⁇ (n) represented by equation (3) and its (complex) conjugate transposed vector ⁇ H (n). .
- a method of taking an addition average after multiplying each sample in the time direction by an arbitrary weight such as a trapezoidal weight may be used.
- the spatial covariance matrix R (n) output from the matrix calculation unit 300 is then input to the adaptive weight calculation unit 301.
- the adaptive weight calculator 301 calculates the weight vector w (n) using the MVDR method.
- focus peripheral information is obtained by using a plurality of steering vectors a.
- the steering vector a p is expressed by Expression (4).
- the steering vector a p is a direction vector having K vector elements (0 to (K ⁇ 1)) equal to the number of active channels, and the received frequency f p and the element 106 It is represented by a function of angles (hereinafter referred to as steering angles) ( ⁇ p , ⁇ p ) formed by the normal vector direction of the surface and the steering vector.
- steering angles ⁇ p , ⁇ p
- ⁇ p is an opening angle from the normal vector
- ⁇ p is a turning angle from the arrangement direction of the elements 106.
- the steering angle is expressed by the final expression of Expression (4).
- ⁇ p is the wavelength of the sound wave corresponding to the frequency f p
- d is the distance between the element centers of the element 106 (element pitch).
- the weight vector w (n) by the MVDR method for the direction of the steering vector a p as described above is calculated by the equation (5) here. Therefore, by calculating the weight vector w (n) for the different steering vectors a p by the weight vector calculation units 3021 to 3023, the number P of steering vectors a p (that is, the number of weight vector calculation units 3021 to 3023) is calculated. As many adaptive weight vectors w 1 (n) to w P (n) as possible can be obtained.
- R (n) is the spatial covariance matrix in snapshot n in the time direction generated by equation (2), and the superscript ⁇ 1 represents the inverse matrix.
- the number P of steering vectors may be any value as long as it is an integer of 2 or more.
- the method of selecting the direction of the steering vector is also arbitrary, in this embodiment, the focal point is on the central axis of the active channel 201, that is, when the number of active channels is an even number, the K elements 106 constituting the active channel 201 are focused. It is assumed that it is located on a normal line (normal line orthogonal to the active channel plane) passing through the midpoint between the K / 2nd and (K + 2) / 2nd elements.
- a normal line passing through the center of the (K + 1) / 2th element 106 among the K elements 106 constituting the active channel 201 (a method orthogonal to the active channel plane). Line).
- One of the steering vectors is preferably oriented in the focal direction.
- the other (P-1) steering vectors are set so as to point in an arbitrary direction within a region around the focal point in the imaging space.
- the total number of steering vectors is arbitrary, and the steering angles thereof are irregular regardless of whether they are equiangular pitches or follow non-linear rules.
- a random selection method may be used.
- the adaptive weight vector w (n) for the focal direction matches the adaptive weight vector w (n) according to the conventional MVDR method.
- the configuration matches the configuration of an adaptive beamformer using a normal MVDR method.
- a different plurality of steering vectors a p are set, and adaptive weight vectors w 1 (n) to w P (n) are obtained for each of them.
- the P adaptive weight vectors w 1 (n) to w P (n) output from the adaptive weight calculator 301 are input to the peripheral information synthesizer 206.
- the peripheral information combining unit 206 adds the P adaptive weight vectors w 1 (n) to w P (n) in the weight combining unit 306, takes an arithmetic average, and combines the combined weight vectors w sum as shown in equation (6). (n) is calculated.
- a fixed apodization multiplication unit 305 may be arranged in the previous stage for obtaining the combined weight vector w sum (n), and each of the weight vectors w p (n) may be multiplied by a fixed apodization.
- the fixed apodization multiplication unit 305 provides fixed apodization weights b 1 , b 2 , and b 3 .
- the inner product calculation unit 307 includes a multiplication unit 3071 and an addition unit 3072, and is a combination of the combined weight vector w sum (n) and the delayed received data vector x (n) transmitted from the delay circuit 204 through the bypass line 207.
- the inner product calculation is performed as shown in Equation (8) to obtain a phasing output y (n).
- multiplication section 3071 calculates the product (1 to K) between the respective channel elements of the weight vector and the delayed received data vector.
- the adder 3072 calculates the sum of the K products obtained by the multiplier 3071 and sets it as the final output (phased output y (n), scalar value) of the peripheral information synthesizer 206.
- the phasing output y (n) of the inner product calculation unit 307 is transmitted to the image processing unit 109.
- N adjustments of y (1), y (2),... Y (N) are performed.
- the phase output y (n) is also sequentially transmitted to the image processing unit 109.
- the reception beamformer 108 completes the calculation of the phasing output y (n) of one raster corresponding to the active channel 201.
- the phasing output y (n) of one raster obtained by the equation (8) is acquired while shifting the active channel from the active channel 201 to the active channels 202 and 203 on the receiving array, and is obtained for each raster.
- the image is transferred to the image processing unit 109.
- all rasters are arranged and a two-dimensional image is generated by a scan converter corresponding to the scanning method of the ultrasonic probe 101.
- various back-end image processing such as various filter processing and calculation of measurement application are performed.
- the ultrasonic image and the calculation result of the measurement application are displayed on the image display unit 103.
- a plurality of steering vectors a p are set, adaptive weight vectors w 1 (n) to w P (n) are obtained for each, and a combined weight vector w sum obtained by combining these.
- the received data vector x (n) is phased to obtain one phased output y (n).
- FIG. 5A shows an example in which reception data is received once from the active array 201 while focusing on the reception focal point 2201, and reception beam forming is performed.
- reception beam forming is performed by the conventional delay addition method
- a signal in the area indicated by the point spread function 2202 around the reception focal point 2201 can also be acquired by one reception beam forming.
- adaptive beamforming is performed, the beam directivity is sharper than that of the delay addition method, and therefore the signal range that can be obtained by one reception beamforming is expressed by a point spread function 2204, It becomes smaller than the point spread function 2202 of the addition method.
- the signal of the object at the focal point 2201 is considered to be randomly distributed around the focal point 2201 like the elliptical region 2205 by being reflected at the focal point 2201.
- This elliptical area 2205 corresponds to the peripheral part (peripheral area) of the first embodiment.
- adaptive weight vectors obtained by the adaptive beamformer are synthesized for a plurality of points (here, three points) having different steering angles.
- This processing is equivalent to overlapping the range of the point spread function 2204 around a plurality of points as shown in FIG. Therefore, since the region 2206 obtained by superimposing the three point spread functions 2204 whose positions are shifted is a range in which a signal can be acquired by one reception in the second embodiment, the region 2206 is scattered in the elliptical region 2205 around the focal point 2201. The majority of the signal of the target object can be acquired, and the missed signal can be reduced.
- the asymmetry of the correlation noise on the left and right of the central axis can be used.
- the correlated noise signal in the received signal can be canceled (decorrelated) by combining the adaptive weights of different steering angles, the information on the focus signal is not missed, while the correlation is high
- a noise signal can be selectively removed.
- the second embodiment improves the demerits of being susceptible to signal loss and wavefront distortion due to the sensitive beam directivity of the adaptive beamformer, and provides robustness that is less susceptible to these effects. be able to.
- this reception beamforming is performed for each reception (for each raster) as shown in FIG.
- noise included in the phasing outputs 1504, 1505, 1506 of the rasters 1501, 1502, 1503 can be reduced, and the influence of wavefront distortion can be reduced. Therefore, the ultrasonic image obtained by arranging the phasing outputs of these rasters is a clear image with less noise, in-vivo sound speed non-uniformity, and reduced wavefront distortion due to scatterer distribution and body motion effects. Become.
- Subarray matrix R ⁇ SUBL in space average calculation is expressed by the product of the partial space vector xi] ⁇ l (n) (equation (10)).
- Subspace vector ⁇ ⁇ l (n) is the delayed received data for K active channels (in this case, represented by generalized complex signal vector ⁇ (n), but instead, real signal vector x (n) May be a vector obtained by extracting some components (for L elements). Therefore, the total number of subspace vectors is KL + 1 (0 ⁇ l (l is lowercase L) ⁇ KL + 1).
- a forward-backward spatial averaging method may be used.
- backward subarray matrix R ⁇ SUBL (n) as in equation (13) can be obtained by the product of the backward subspace vectors xi] ⁇ l (n).
- the backward subspace vector is expressed by Expression (14).
- the backward subarray space matrix R 1 to (n) can be calculated as the equation (15) by the same operation as the forward space average.
- the subarray space matrix R ⁇ (n) in the case of the above-mentioned forward spatial average is arithmetically averaged, so that the forward / backward subarray spatial covariance matrix R FB (n ).
- this forward / backward subarray spatial covariance matrix R FB (n) is used by replacing R (n) in Expression (5) in the calculation in the adaptive weight calculation unit 301,
- the weight vector w p (n) can be calculated as in Expression (17).
- the output of the matrix calculation unit 300 is L ⁇ L in size, and the number of elements constituting the weight vector w p (n) is L.
- a plurality (P) of weight vectors w p calculated by Expression (12) or Expression (17) using the spatial averaging method are transmitted to the peripheral information combining unit 206 as in the case where the spatial averaging method is not used.
- the weight synthesizer 306 adds weight vectors to calculate w sum (equation (7), equation (8)). Further, the inner product calculation unit 307 performs inner product calculation with the delayed received data transmitted from the bypass input 207, and is finally delivered to the image processing unit 109 as phasing data y (n).
- a dimension compression unit 308 is arranged in the middle of the bypass line 207 and before the inner product calculation unit 307, and the dimension compression unit 308 starts from ⁇ (n) having K elements. It can be configured to generate a vector g (n) having L elements. Note that the dimension compression unit 308 can also be arranged inside the peripheral information synthesis unit 206.
- Equation (19) The phasing output y (n) when the spatial averaging method is used is expressed by Equation (19).
- an effect of suppressing the correlation noise included in the ultrasonic reception signal can be obtained. Therefore, an ultrasonic image with less influence of noise can be obtained by combining the spatial average processing of the subarray matrix with a configuration for combining a plurality of adaptive weight vectors around the focus.
- the MVDR has been described as an example of the adaptive beamforming technique.
- the algorithm applied in the adaptive weight calculation unit 301 uses the spatial covariance matrix calculated by the matrix calculation unit 300. Any algorithm may be used, and any other method such as the MMSE method, the APES method, the Eigenspace-MV method (ESMV, EIBMV) using the spatial covariance matrix and its eigenvalue / eigenvector, the ESPRIT method, and the MUSIC method may be used.
- the peripheral information synthesis unit 206 multiplies the delayed received signal by adaptive weights for the plurality of points obtained by the peripheral information calculation unit 205, and A plurality of inner product calculation units for generating a pre-synthesis phased output for each adaptive weight, and an output synthesis unit.
- the output combining unit adds and combines the pre-combination phasing outputs for each of the plurality of adaptive weights to generate a phasing output used for generating image data.
- a fixed apodization multiplication unit that multiplies adaptive weights for a plurality of points obtained by the peripheral information calculation unit 205 by predetermined fixed weights, respectively.
- a fixed apodization multiplication unit that multiplies each of the delayed received signals after phasing by a plurality of adaptive weights by a predetermined fixed weight between the plurality of inner product calculation units and the output synthesis unit. It is also possible to do.
- a fixed apodization multiplication unit that multiplies a plurality of signals after multiplication by a predetermined fixed weight between a plurality of multiplication units and a plurality of addition units constituting a plurality of inner product calculation units. It is also possible to arrange.
- the peripheral information synthesizing unit 206 is configured to generate the final phasing output used by the image processing unit for image generation by synthesizing information obtained from the reception signal of one active channel by the peripheral information calculating unit 205. For example, the peripheral information synthesizing unit 206 adds the pre-composition phasing output for each of the plurality of adaptive weights obtained by the peripheral information calculating unit 205 from the received signal of one active channel, and the image processing unit 109 uses it for image generation. A final phasing output is generated for each active channel.
- the ultrasonic imaging apparatus will be specifically described with reference to FIG.
- the configuration of the peripheral information combining unit 206 is different from that of the second embodiment, different parts will be described, and description of other parts will be omitted.
- the peripheral information combining unit 206 combines a plurality of adaptive weight vectors w 1 (n), w 2 (n), and w 3 (n) by the weight combining unit 306.
- the delay-and-sum processing of the delayed received data x (n) is performed by the obtained composite weight w sum (n). Since this calculation is a linear calculation, the order of the weight synthesis and the phasing addition process may be reversed. Therefore, in the plurality of inner product calculation units 307 of the third embodiment, as shown in FIG.
- a plurality of multiplication units (3071-1, 3071-2, and 3071-3) and elements after multiplication (delayed data) )
- One set of each of the multiplication unit 3071 and the addition unit 3072 constitutes one inner product calculation unit for each adaptive weight vector, and a plurality of these sets are prepared to form a plurality of inner product calculation units 307. is doing.
- the delayed received data x (n) is input by the bypass inputs 2071, 2072, and 2073, respectively.
- the output composition unit 500 is arranged after the plurality of inner product calculation units 307.
- the plurality of inner product arithmetic units 307 perform phasing addition of the bypass input 2071 using the weight vector w 1 (n) and phasing addition of the bypass input 2072 using the weight vector w 2 (n). Then, by phasing and adding the bypass input 2073 using the weight vector w 3 (n), a plurality of phasing outputs y 1 (n), y 2 (n), y before synthesis corresponding to each steering vector 3 (n) is calculated (the following formula (20)).
- the output synthesis unit 500 takes the arithmetic average of the phasing outputs y 1 (n), y 2 (n), y 3 (n) corresponding to each steering vector, and synthesizes them as shown in equation (21). A phasing output y sum (n) is obtained.
- a fixed apodization multiplication unit 305 can be added also in this embodiment.
- the fixed apodization multiplication unit 305 is arranged before the inner product calculation unit 307.
- the position of the fixed apodization multiplication unit 305 is not limited to the position of FIG.
- the multiplication units 3071-1 to 37-1 in the inner product calculation unit A fixed apodization multiplication unit 305 may be added between the addition units 3072-1 to 3072-1. In either case, the final combined phasing output y sum (n) can be expressed by Expression (22).
- the difference between the second embodiment and the third embodiment is whether the weight is synthesized and then the inner product operation is performed or the phasing output after the inner product operation is synthesized.
- the embodiment is different.
- the final output is the same because of the linear operation, and the final phasing output of the second embodiment is y (n) in the equation (19) and the final output of the third embodiment.
- the combined phasing output y sum (n) is equal.
- the delay unit (delay circuit 204) generates a delayed reception signal for each of a plurality of different reception focal points
- the peripheral information calculation unit 205 and the peripheral information synthesis unit 206 are configured to receive a plurality of receptions. For each focal point, information is acquired about a plurality of points in the reception focal point and a peripheral region of the reception focal point, and a phasing output is generated.
- the peripheral information calculation unit 205 and the peripheral information combining unit 206 are arranged for each delay unit, and are configured to acquire information from the delayed received signal for each reception focus generated by the delay unit and generate a phasing output.
- a part of the plurality of reception focal points for the reception signal of the active channel at a certain time point can be configured to overlap with a part of the plurality of reception focal points for the active channel at different time points.
- the ultrasonic imaging apparatus according to the fourth embodiment of the present invention will be specifically described with reference to FIGS. A description of the same configuration as that of the second and third embodiments will be omitted, and only different portions will be described.
- the focal point is located on the central axis of the active channel (the normal to the active channel surface passing through the active channel center position) for each received signal.
- a plurality of focal points are set as focal points at a certain depth for each received signal.
- the focal point is set for points 601 and 603 on different axes in the vicinity of the central axis 600 on both sides of the point 602 on the central axis 1600, and information on the subject 100 is collected.
- Focusing on a plurality of focal points 601, 602, and 603 for each received signal means that the distribution of delay times of received signals of the K elements 106 constituting the active channel is changed to the focal points 601, 602. This is realized by setting the concave molds 604, 605, and 606 with 603 as the center.
- peripheral information is collected, and calculation for combining the peripheral information is performed.
- FIG. 9 shows a block diagram of the receive beamformer 108 of the fourth embodiment.
- the reception beamformer 108 has the same number of delay circuits as the number of focal points to be set.
- 2041, 2042, and 2043, peripheral information calculation units 2051, 2052, and 2053, and peripheral information combining units 2061, 2062, and 2063 are provided.
- the configurations of the peripheral information calculation units 2051, 2052, 2053, and the peripheral information combining units 2061, 2062, 2063 are the same as the configurations of the peripheral information calculating unit 205 and the peripheral information combining unit 206 of the second embodiment, respectively. Therefore, the number of phasing outputs y (n) is also the focal point and the number, and in the example of FIG. 9, there are three y1 (n), y2 (n), and y3 (n).
- the delay circuit 2042 focuses on the focal point 602 on the central axis 1600 of the active channel 201 at the time of a certain snapshot n, and the delay circuits 2041 and 2043 respectively focus on focal points 601 and 603 on the axis deviated from the central axis 1600.
- a delayed concave surface is formed, and delayed received data is generated.
- the peripheral information combining unit 2051 performs the calculation described in the second embodiment on the delayed received data from the delay circuit 2041 as illustrated in FIG. 10A, so that the focal point 601 and the direction of the focal point 601 are obtained.
- Weight vectors w1 1 (n), w1 2 (n), and w1 3 (n) are calculated for two points forming the steering angle 1602, respectively.
- the peripheral information combining unit 2061 combines the weight vectors w1 1 (n), w1 2 (n), and w1 3 (n), and outputs the phasing output y1 (n) using the combined weight vector w1 sum (n). obtain.
- the peripheral information combining unit 2052 performs calculation on the delayed received data from the delay circuit 2042 as illustrated in FIG. 10B, and the two points that form the steering angle 1604 with respect to the direction of the focal point 602 and the focal point 602. , W2 1 (n), w2 2 (n), and w2 3 (n) are calculated respectively.
- the peripheral information combining unit 2062 combines the weight vectors w2 1 (n), w2 2 (n), and w2 3 (n), and outputs the phasing output y2 (n) using the combined weight vector w2 sum (n). obtain.
- the peripheral information combining unit 2053 performs calculation on the delayed received data from the delay circuit 2043 as illustrated in FIG. 10C, and forms two steering angles 1606 with respect to the focus 603 and the direction of the focus 603.
- Weight vectors w3 1 (n), w3 2 (n), and w3 3 (n) are calculated for the two points, respectively.
- the peripheral information combining unit 2062 combines the weight vectors w3 1 (n), w3 2 (n), and w3 3 (n), and outputs the phasing output y3 (n) using the combined weight vector w3 sum (n). obtain.
- ultrasonic image data of the three focal points 601, 602, and 603 can be obtained using only one transmission / reception signal for the active channel 201.
- Such processing means that if the direction of the focal point 602 is the main beam direction, an image can also be created in the sub-beam direction (the directions of the focal points 601 and 603). That is, by using the configuration of the fourth embodiment, it can be understood that the configuration of the present invention in which signals of an object scattered around the focal point can be applied to techniques such as sub-beam processing and parallel beam forming. Therefore, it is possible to obtain an effect of reducing signal loss and noise reduction.
- this processing can be applied to a plurality of active channels 201, 202, 203.
- data can be acquired again for the focal point 602 by setting the focal point at the same point as the focal point 602 in the active channel 201 so that the other active channels 202 and 203 look obliquely. Therefore, it is possible to superimpose a plurality of active channels 201, 202, and 203 by acquiring data for viewing the same point 602 from different directions and combining the acquired data.
- Such processing is called aperture synthesis (synthetic aperture) processing. That is, by using the configuration of the fourth embodiment, it is possible to apply the present invention for acquiring the signal of the object scattered around the focal point to the aperture synthesis process.
- the number of focal points may be any number equal to or greater than two.
- the same number of delay circuits, peripheral information calculation units, and peripheral information combining units as the focal points are arranged.
- the present embodiment it is possible to obtain the effect of reducing the wavefront distortion according to the present invention even for off-axis phasing technology for focusing off-axis, parallel beam forming, and aperture synthetic beam forming.
- peripheral information calculation units 2051 to 2053 and the peripheral information combining units 2061 to 2063 in FIG. 9 in this embodiment are the same as those shown in FIG. 3, but are within the scope of the embodiment of the present invention. Any configuration may be used as long as it is present.
- peripheral information synthesis unit 206 of the third embodiment shown in FIG. The configurations of the peripheral information calculation unit 205 and the peripheral information combining unit 206 in each of the fifth to eighth embodiments to be described later can also be used.
- the pre-synthesis phased output for each of the plurality of adaptive weights generated by each inner product calculation unit is stored between the plurality of inner product calculation units and the output synthesis unit of the third embodiment.
- a storage unit temporary storage unit 800
- An ultrasonic imaging apparatus according to the fifth embodiment of the present invention will be specifically described with reference to FIG.
- the fifth embodiment has the same configuration as that of the third embodiment, except that a temporary storage unit 800 is arranged between the plurality of inner product calculation units 307 and the output synthesis unit 500 of the peripheral information synthesis unit 206. This is different from the third embodiment.
- a temporary storage unit 800 is arranged between the plurality of inner product calculation units 307 and the output synthesis unit 500 of the peripheral information synthesis unit 206. This is different from the third embodiment.
- the signal to be phased sequentially flows from the reception on the active channel to the image processing unit 109 with time. Therefore, the synthesizing process in the output synthesizing unit 500 is performed in each of the phasing outputs y obtained by the plurality of inner product calculation units 307 for a plurality of steering vectors for the data of the current snapshot n in the currently phased raster.
- the process is limited to the synthesis of 1 (n), y 2 (n), and y 3 (n).
- the temporary storage unit 800 is arranged as a memory unit for storing the phasing outputs y 1 (n), y 2 (n), and y 3 (n) before the output synthesis unit 500. .
- the temporary storage unit 800 stores the phasing outputs y 1 (n), y 2 (n), and y 3 (n) obtained by the plurality of inner product calculation units 307 every time the snapshot n is updated or the raster is updated. .
- the output synthesis unit 500 can read out and synthesize the phasing outputs y 1 (n), y 2 (n), y 3 (n) stored in the temporary storage unit 800. It is possible to synthesize phasing outputs of different time samples or phasing outputs between different rasters.
- these phasing outputs y 1 (n), y 2 (n), y 3 (n) are weight vectors w 1 (n), w 2 obtained for different steering vectors at each snapshot time n.
- the phasing outputs are phasing with (n) and w 3 (n), respectively. Therefore, the phasing outputs y 1 (n), y 2 (n), and y 3 (n) are phasing outputs in the directions of the respective steering vectors. Therefore, the output combining unit 500 can select and combine the phasing outputs in different steering vector directions in different combinations between different rasters.
- the phasing outputs in different steering vector directions to be selected can be determined in advance or can be set by the operator.
- the angle adjustment unit 502 and the memory output adjustment unit 503 can be arranged in the reception beamformer 108.
- the angle adjustment unit 502 controls the weight vector calculation units 3021, 3022, and 3023 of the adaptive weight calculation unit 301 to set the steering vector direction, number, and density.
- the memory output adjustment unit 503 reads out the phasing output from the temporary storage unit 800 and delivers it to the output synthesis unit 500 (clock), and specifies the address in memory of the phasing output to be read out.
- an adjustment unit 501 that can adjust the parameters of both the angle adjustment unit 502 and the memory output adjustment unit 503 in conjunction with each other may be further provided.
- peripheral information combining unit 206 provided with the temporary storage unit 800 of this embodiment can be arranged in place of the peripheral information combining units 2061, 2062, and 2063 of the fourth embodiment. This makes it possible to combine the phasing outputs of different time samples of the same raster with respect to a plurality of focal points, synthesize the phasing outputs between different rasters, and the like.
- the output synthesizing unit 500 acquires and adds phasing outputs before synthesis that are different in the time direction from the storage unit (temporary storage unit 800) and the plurality of inner product calculation units, and generates a phasing output. .
- the ultrasonic imaging apparatus according to the sixth embodiment of the present invention will be specifically described with reference to FIGS.
- the present embodiment includes a temporary storage unit 800 as in the ultrasonic imaging apparatus of the fifth embodiment, but the internal configuration of the temporary storage unit 800 is different from that of the fifth embodiment.
- the temporary storage unit 800 includes a memory unit 900 and a bypass line 901 for phasing output (hereinafter also referred to as post-phasing data) y 1 (n) to y P (n).
- the memory unit 900 adjusts each of the post-phasing data y 1 (n) to y P (n) corresponding to the P steering vectors for N snapshots in the time direction, that is, as many as P ⁇ N. It has a storage area (address) for storing post-phase data.
- P memory areas 9001, 9002, and 9003 having N addresses are prepared in the memory unit 900. Since other configurations are the same as those of the fifth embodiment, the description thereof is omitted.
- Temporary storage unit 800 receives post-phasing data y 1 (n) to y P (n) from a plurality of inner product calculation units 307 at a snapshot time n.
- the post-phasing data y 1 (n) to y P (n) are sent to both the memory unit 900 and the bypass line 901.
- the memory unit 900 stores post-phasing data y 1 (n) to y P (n) at predetermined write addresses (WA) in the memory areas 9001, 9002, and 9003, respectively.
- y 1 (ni) to y P at sample times n-1 to ni stored in the read address (RA) (ni) is output.
- the output that combines the read output at sample times n-1 to ni and the phasing data y 1 (n) to y P (n) of the current snapshot n via the bypass line is output from the output bus 902. It is transmitted toward the unit 500.
- the phasing outputs of the plurality of steering vectors 1 to P are obtained as n-1...
- Ni information that is, y 1 (n ⁇ 1) to y P (n ⁇ 1) even at the sampling time n. ⁇ ⁇ ⁇ Y 1 (ni) to y P (ni) can be traced back and output along with the current post-phasing data y 1 (n) to y P (n). Can be combined together.
- the information stored in the memory unit 900 is overwritten every time the raster is updated, and is sequentially overwritten by the post-phasing data of the new raster.
- the post-phasing data y 1 (n) to y P (n)... Y 1 (ni) to y P (ni) at a plurality of sample times n to ni are synthesized. Therefore, phasing outputs 1703, 1704, and 1705 obtained by synthesizing the phasing results of a plurality of time samples in the same raster 1501, 1502, 1503 can be obtained. As a result, it is possible to collect inhomogeneous scattering data 1701 and 1702 scattered in the time direction, and it is possible to compensate for image quality degradation of the ultrasonic image in the time direction (ultrasonic propagation direction).
- the read address (RA) can be determined in advance, or can be configured to include the output adjustment unit 503 and be specified by the output adjustment unit 503.
- the output adjustment unit 503 for example, information on which data is output from the current post-phasing data y 1 (n) to y P (n) and the previous data held in the memory unit
- the signal for controlling the output timing and the output data address can be sent to the temporary storage unit 800.
- the plurality of inner product calculation units generate a pre-synthesis phased output for each of the plurality of adaptive weights for each of the active channels
- the output synthesis unit 500 includes a storage unit (temporary storage unit 800) and Pre-combination phasing outputs generated for different active channels are acquired from a plurality of inner product calculation units and added to generate a phasing output.
- the ultrasonic imaging apparatus according to the seventh embodiment of the present invention will be specifically described with reference to FIGS. The description of the same configuration as that of the sixth embodiment is omitted.
- the ultrasonic imaging apparatus according to the seventh embodiment has a phasing output (data after phasing) y 1 (n) to y P (n) corresponding to P steering vectors of three rasters.
- Memory units 900-1, 900-2, and 900-3 are stored in the temporary storage unit 800.
- the structures of the memory units 900-1, 900-2, and 900-3 are the same as those of the memory unit 900 of the sixth embodiment. Thereby, post-phasing data for three rasters can be stored.
- the current post-phasing data y 1 (r, n), y 2 (r, n), y 3 (r, n) at the snapshot time n is written to the write address WA in the figure and at the same time the bypass line 901. To the output composition unit 500.
- the post-phasing data at one snapshot time n has three data y 1 (n) to y 3 for each raster. (n).
- Other configurations are the same as those of the sixth embodiment, and thus description thereof is omitted.
- 1) to y 1 (r, N), y 2 (r, N)... Y p (r, N) can be held in the memory section 900-1, and the r ⁇ 1, r ⁇ th, Data after phasing of two rasters can also be stored in the memory units 900-2 and 900-3, respectively.
- post-phasing data having different rasters, snapshot times, and steering vectors are stored in the memory units 900-1 to 900-3. Therefore, a desired phasing is selected from these post-phasing data.
- One or more pieces of post-phase data can be selected and output to the output combining unit 500.
- a certain steering angle ⁇ from the rth raster 1 (P 1)
- the data can be output to the bus 1007.
- the combining process performed by the output combining unit 500 can be expressed by Expression (23).
- post-phasing data obtained from a plurality of rasters can be synthesized for the same focal point 1500, compared with the case where post-phasing data obtained from one steering vector of a single raster is used.
- the post-phasing data stored in the memory units 900-1, 900-2, 900-3 is overwritten every time the raster is updated, and the post-phasing data for the latest three rasters is always stored in the memory unit 900-. 1, 900-2, 900-3.
- the output adjustment unit 503 can be arranged.
- the output adjustment unit 503 sends a data address to be output to the output synthesis unit 500 and a signal for controlling the timing to the temporary storage unit 800.
- FIG. 16 In the ultrasonic diagnostic apparatus of the eighth embodiment, (m + 1) memory units 900-1, 900-2... 900- (m + 1) are arranged in the temporary storage unit 800. ing. Each of the memory units 900-1 to 900- (m + 1) has the same configuration as the memory unit 900-1 and the like of the seventh embodiment. As a result, one or more selected from all post-phasing data in any combination for the three dimensions of multiple steering angle directions 1 to P, snapshot times 1 to N, and raster movement directions 1 to m + 1 in the same raster, The data can be output to the output combining unit 500 and combined.
- the outputs 1106, 1107, and 1108 to be combined can be selected from the P ⁇ N ⁇ (m + 1) post-phasing data and output to the output combining unit 500 to be combined. Note that a description of the same configuration as that of the seventh embodiment is omitted.
- FIG. 17 it is possible to select and synthesize a phasing output that allows an arbitrary angle from not only the same raster but also from adjacent rasters. That is, similarly to the seventh embodiment, by obtaining and synthesizing the phasing data of the same focal point 1500 from different rasters 1501, 1502, and 1503, the data of the focal point 1500 can be repeatedly obtained, and the accuracy of the data Can be improved. Moreover, in the eighth embodiment, unlike the seventh embodiment, post-phasing data for the points 1701 to 1704 around the focal point 1500 is also selected and output, and synthesized by the output synthesis unit 500. it can. Thereby, the information (data after phasing) of the point spread function 2204 around the focal point 1500 can be synthesized as shown in FIG.
- the (r-1) th composite data is sent to the image processing unit 109 at the time of the rth raster currently being scanned. Therefore, there is a slight delay time in the image processing timing.
- the peripheral information calculation unit 205 generates a pre-synthesis phased output for a plurality of points set in the two-dimensional direction.
- the storage unit (temporary storage unit) 800 stores these pre-synthesis phased outputs.
- the ultrasonic imaging apparatus of the ninth embodiment has a configuration in which Q temporary storage units 800 of FIG. 16 of the eighth embodiment are further arranged in parallel.
- the probe 101 one in which elements 106 are arranged in a two-dimensional direction is used.
- a two-dimensional array probe, a two-dimensional array convex probe, a short axis bundle probe, or the like is used.
- a probe (mechanical scan 3D probe) having a form in which the dimension can be increased by moving the one-dimensional array mechanically and the same information as the two-dimensional array can be acquired may be used.
- the description of the same configuration as in the eighth embodiment is omitted.
- the adaptive weight calculation unit 301 changes the calculated steering angle in a two-dimensional direction in which the elements 106 of the probe 101 are arranged, and obtains an adaptive weight vector for P ⁇ Q steering angles.
- the Q temporary storage units 800 post-phasing data obtained for the adaptive weight vector of each steering vector is stored for each snapshot time. Thereby, peripheral information (data after phasing) can be collected three-dimensionally.
- the steering vector a p represented by the above-described formula (4) is expanded in the two-dimensional direction, and particularly in the case of two-dimensional angles 1 ⁇ p ⁇ P and 1 ⁇ q ⁇ Q.
- the direction of the steering vector can be selected with an arbitrary set of ( ⁇ p , ⁇ q ).
- the steering vector a (p, q) in this case can be expressed as in Expression (24).
- FIG. 19 shows a signal processing operation according to the ninth embodiment.
- the subject 100 is viewed from the center element 2001 on the two-dimensional array of probes 101, information on the target point on the spherical surface with a line segment connecting the reception focus 2002 and the element center 2001 as a radius is collected.
- P ⁇ Q steering vectors a (p, q) corresponding to the two-dimensional angle ( ⁇ p , ⁇ q ) can be set on the spherical surface.
- FIG. 20 shows a setting example of the steering vector by the adjusting unit 501 in FIG.
- the direction of the steering vector may use a fixed angle in advance, or may be configured to be changed under the control of the adjustment unit 501.
- 20 (a) to 20 (f) show embodiments of two steering angles.
- the steering angle 2104 is fixed.
- the depth of the imaging target point (data collection point) increases from FIG. 20A to FIG. 20C, the steering angle 2104 is fixed.
- the range in which signals can be collected around the imaging target point increases.
- FIGS. 20D to 20F are examples in which the steering angle is decreased as the imaging target point becomes deeper, and the horizontal distance 2108 (spreading of the imaging target point) from the reception focal point at which the steering vector is expected is constant. . That is, in one raster, by changing the steering angle of the steering vector for each depth of the imaging target point (snapshot time n), the spread of the imaging target point (point spread function 2204) is changed to the depth of the imaging target point. Regardless, it can be kept constant.
- the adjustment unit 501 can set the number of steering vectors set between the steering angles 2109 and the angles of adjacent steering vectors.
- the adjustment unit 501 can select a suitable combination of steering angle and number and set it to the angle adjustment unit 502 and the output adjustment 503. Further, a configuration in which the operator selects a combination of the steering angle and the number via the console 110, or a configuration in which the operator inputs an arbitrary steering angle value and the number to the adjustment unit 501 via the console 110 may be adopted. Is possible.
- the eleventh embodiment further includes a memory unit that stores in advance the distribution of the spread angle of the steering vector to a plurality of points according to the position of the reception focus in the image generated by the image processing unit.
- the adjustment unit obtains the adaptive weight using the steering vector of the divergence angle read from the memory unit according to the position of the reception focus.
- the divergence angle of the steering vector stored in the memory unit is set to be small at the reception focus near the transmission focus of the ultrasonic signal to the subject and large at the reception focus away from the transmission focus.
- the divergence angle of the steering vector is set smaller at the edge portion of the image than at the central portion of the image.
- the adjustment unit 501 generates an image by applying different steering angles for each imaging target point in one image according to the type, imaging condition, and imaging sequence type of the ultrasound probe 101.
- a set of different steering angles is applied corresponding to a probe such as a linear probe, a convex probe, a sector probe, a two-dimensional array probe, and a mechanical 3D imaging probe.
- Generate an image is generated by applying a set of different steering angles for each imaging sequence and imaging conditions such as transmission / reception focus, transmission / reception frequency, frame rate, number of parallel beam forming, tissue harmonics, contrast contrast, etc. .
- the focal point at the time of transmission is set to one or a plurality of predetermined positions when imaging one ultrasonic image, and the reception focal point is changed by the reception beamformer 104, so that the inside of the imaging region is set. An image of each imaging target point is acquired. Therefore, the focal point 2210 at the time of transmission is a fixed position in the image.
- the signal of the object spreads to the elliptical region 2205 around the reception focal point due to reflection at the reception focal point, and the size is the depth of the reception focal point. It depends on (time). As shown in FIG.
- the size of the elliptical area 2205 is small when it is close to the focal point 2210 at the time of transmission, and is large when it is away from the focal point 2210 at the time of transmission. Therefore, when receiving beam forming is performed with a steering vector having a fixed steering angle (expansion angle) 2211, an elliptical area 2205 indicating the signal spread of the object and a point spread function 2204 capable of acquiring a signal determined by the spread angle of the steering vector. There is a mismatch with the size of.
- the signal intensity acquired at the image edge portion is lower than that at the image center portion.
- the mismatch and the point spread function 2204 protruding from the image appear as uneven brightness in the B-mode image.
- FIG. 22B the structure around the transmission focal point 2210 is depicted with uniform brightness, but in the shallow part, the deep part, and the image edge part, an elliptical area 2205 indicating the signal spread of the target object Due to the mismatch with the point spread function 2204 and the protrusion of the point spread function 2204 from the image, the image luminance is lowered or uneven.
- the steering spread angle is changed in accordance with the signal spread region of the object, and the point spread function 2204 at the time of beam forming is changed to an elliptic region. 2205. Further, in the image edge portion, a steering angle is used for suppressing the point spread function 2204 determined by the spread angle of the steering vector from being unnecessarily spread out of the imaging target.
- the adjustment unit 501 includes a memory 2215 as shown in FIG.
- the memory 2215 stores in advance data such as a table format that determines the distribution of the steering vector spread angle (steering angle) in the raster scan direction and depth (time) direction.
- the steering angle distribution depends on the type of the probe 101 and the imaging sequence (focus on transmission, etc.) It is determined to improve image brightness and reduce unevenness.
- the steering angle distribution is prepared for each combination of the type of the probe 101 and the imaging sequence.
- the steering angle distribution that improves image brightness and reduces unevenness can be obtained in advance by calculation (including computer simulation), experiment, or the like.
- the adjustment unit 501 receives the type of the probe 101 and the information of the imaging sequence from the control unit 111, and reads data defining the distribution of the spread angle of the steering vector corresponding to these from the memory 2215. According to the read data, the adjustment unit 501 sets the steering vector spread angle (steering angle) in the angle adjustment unit 502 and the output adjustment unit 503 of the reception beamformer 108 for each scan direction and depth (time) direction of the raster.
- the point spread function 2204 at the time of beam forming can be matched with the elliptical region 2205 indicating the spread of the signal of the object.
- the image edge portion it is possible to suppress the point spread function 2204 from being unnecessarily spread out of the imaging target. Therefore, as shown in FIG. 22D, an image having no luminance unevenness is formed in the image, which is effective in improving the image quality.
- the configuration of the reception beamformer 108 in the eleventh embodiment is not limited to the configuration in FIG. 21, and the reception beamformer 108 that performs reception beamforming in the first to ninth embodiments is used. it can.
- the transmission focus is fixed to one point or a plurality of points.
- the transmission dynamic focus that changes the transmission focus in the depth direction is used, the unevenness of the elliptical region 2205 that indicates the spread of the signal of the object. Therefore, the tenth embodiment can be applied.
- 22 (a) to 22 (d) show images of convex probes, but any other probe type such as a linear probe, sector probe, 2D array probe, mechanical scan 3D probe, etc.
- the tenth embodiment can be applied to any raster scanning method.
- the steering angle distribution is determined in advance.
- the steering angle distribution is obtained from the signal intensity or luminance distribution of the B-mode image. Receive beam forming is performed using the obtained steering angle distribution. That is, the adjustment unit uses a predetermined function to determine the distribution of luminance or signal intensity of the image obtained by setting the steering vector using the distribution of the divergence angle of the steering vector stored in the memory. The steering vector is converted to a spread angle distribution. The difference between the distribution of the steering vector divergence angle in the memory and the distribution of the steering vector divergence angle obtained using the function corrects the distribution of the steering vector divergence angle in the memory or the brightness or signal strength of the image. To do.
- beam forming of the first to ninth embodiments is performed using a fixed steering angle as shown in FIG. 22A, and a B-mode image as shown in FIG. 22B is captured.
- the adjustment unit 501 obtains the luminance distribution within the imaging target or the signal intensity distribution 230 after phasing as shown in FIG.
- the obtained luminance distribution or signal intensity distribution 230 is applied to a predetermined function for converting the luminance / signal intensity distribution into the steering angle distribution, and the steering angle distribution 231 is obtained as shown in FIG.
- the steering angle distribution 231 and performing reception beam forming as shown in FIG. 22C of the eleventh embodiment it is possible to actively improve the luminance unevenness of the image.
- the function for converting the luminance / signal intensity distribution into the steering angle distribution is obtained in advance by experiment or calculation (including computer simulation) and stored in the memory 2215 or the like as a database such as a table or a mathematical expression.
- the received beam forming is performed as shown in FIG. 22C using the distribution of the steering vector spread angle (steering angle) stored in the memory 2215, and the obtained figure.
- the luminance distribution within the imaging target or the signal intensity distribution 230 after phasing may be obtained as shown in FIG.
- the obtained luminance or signal intensity distribution is converted into the steering angle distribution 231 using the above-described function.
- the obtained steering angle distribution 231 is compared with the steering angle distribution 232 used for generating the B-mode image in FIG. 22D, a difference between the two is obtained, and this difference is reflected in the phasing process.
- a method of correcting the steering angle distribution stored in the memory 2215 according to the difference and performing reception beam forming again using the corrected steering angle distribution, or the brightness of the B-mode image according to the obtained difference Or a method of correcting the distribution of signal intensity.
- the steering angle distribution 232 or the image luminance and signal intensity distribution 230 can be fed back and corrected.
- the steering angle distribution can be obtained actively or corrected by a feedback operation.
- This process is not only for luminance unevenness when using an adaptive beamformer, but also for luminance unevenness due to energy nonuniformity of the acoustic signal that occurs in phasing by conventional DAS (Delay And Sum: delay addition processing). Is also applicable. Therefore, it can contribute to the improvement of the essential ultrasonic B-mode image.
- FIG. 24 shows an example of the console 110 of the ultrasonic imaging apparatus of the present invention.
- an operation for changing the number and density of steering vectors, the angle direction of the phasing result to be synthesized, the time direction, the number of raster directions, etc. on the console 110 of the ultrasonic diagnostic apparatus As the portion, knob portions 1301, 1302 and the like with scales can be arranged.
- the switch unit 1303 can be arranged on the console 110 as an operation unit for performing switching of the steering angle mode, switching of whether to synthesize the phasing result, and the like described in the above embodiments.
- the operator can change various parameters of the phasing / synthesizing process while observing the actual ultrasonic image, and can perform diagnosis under optimum conditions for each subject 100.
- the set value may be displayed on a part 1304 of the display area of the image display unit 103.
- the console mode switching unit may be switched according to switching with different probes, imaging conditions, and imaging sequence. For example, a set of different steering angles is applied corresponding to a probe such as a linear probe, a convex probe, a sector probe, a two-dimensional array probe, and a mechanical 3D imaging probe. Even a switching unit capable of switching a mode capable of generating an image, for example, transmission / reception focus, transmission / reception frequency, frame rate, number of parallel beam forming, tissue harmonics, contrast contrast A switching unit that can generate an image by applying a different set of steering angles for each imaging condition or imaging sequence.
- FIG. 25 shows a perspective view of the console 110 and the image display unit 103 of another specific example of the ultrasonic diagnostic apparatus of the present invention.
- the operator performs a manual operation for setting a specific ROI (region of interest) 1401 while referring to a normal ultrasound image (an image not subjected to the synthesis processing of the present invention) 103.
- a unit 1403 (for example, a mouse) is prepared.
- the surgeon can generate an image 1402 that has been subjected to the adaptive processing of the present invention only for a specific ROI 1401.
- the image 1402 that has been subjected to the adaptive processing of the present invention can also be displayed in another area of the image display unit 103.
- the adjustment unit 501 shown in FIG. 11 receives the calculation parameter set by the operation unit 1403, outputs a control signal to the angle adjustment unit 502 and the combined output adjustment unit 503, and the calculation parameter in the adaptive weight calculation unit 301 The output data of the temporary storage unit 800 is adjusted.
- FIG. 26 shows an ultrasonic image (B mode image) obtained by simulation in order to show the effect of the embodiment of the present invention.
- These images are obtained by assuming the six point scatterers as the object 100 and calculating the ultrasonic images by a computer.
- An image 2300 and an image 2303 in FIG. 26A are cases where no wavefront distortion has occurred in the ultrasonic beam.
- An image 2300 is an image obtained by conventional delay-added reception beamforming
- an image 2303 is an image obtained by conventional adaptive beamforming.
- images 2301, 2302, 2304, and 2305 in FIG. 26B are cases where wavefront distortion has occurred in the ultrasonic beam.
- An image 2301 is an image obtained by conventional delay-added reception beamforming
- an image 2304 is an image obtained by conventional adaptive beamforming.
- An image 2302 is a delayed addition method, but is an image obtained by synthesizing a steering vector by a method similar to the method of the second or third embodiment.
- An image 2305 is an image obtained by performing adaptive beamforming and steering vector synthesis by the method of the second or third embodiment.
- the image 2301 obtained by the conventional delay addition method is compared with the image 2302 obtained by synthesizing the steering vector by the conventional delay addition method, even if the steering vector is synthesized by the delay addition method, it is scattered. It can be seen that the band-like high luminance region around the body is hardly improved, and the influence of wavefront distortion cannot be removed.
- FIG. 27 is a graph showing one point scatterer of the images 2301, 2304, and 2305 in FIG.
- the horizontal axis represents the azimuth direction of the receiving array (raster scanning direction), and the vertical axis represents the luminance level of the image.
- the position at which the horizontal axis is exactly zero is the position of the dot-scattering body.
- the dotted line 2401 in the graph of FIG. 27 is the luminance distribution of the conventional delayed addition image 2301, the alternate long and short dash line 2402 is the luminance distribution of the conventional adaptive beamformer image 2304, and the solid line 2403 is the second or third.
- the luminance distribution of the image 2305 of the adaptive beamformer subjected to the steering composition of the embodiment is shown.
- the luminance peak position is deviated from the true position of the object as indicated by a dotted line 2401.
- the peak position approaches the true position of the object as indicated by the alternate long and short dash line 2402, but the luminance level is reduced.
- the luminance distribution of the image 2305 according to the second or third embodiment of the present invention maintains the signal intensity equivalent to the luminance distribution (dotted line 2401) of the delay addition method, while maintaining the true position of the object ( The peak intensity is located at the zero position on the horizontal axis, and the effect of wavefront distortion compensation can be confirmed.
- the adaptive beamformer is provided with the ability to avoid information loss due to its sharp beam directivity. Moreover, the ability to cancel and reduce unnecessary correlation noise from the medium around the focal point is provided.
- the present invention can provide an ultrasonic imaging apparatus that is robust against wavefront distortion due to in-vivo sound velocity non-uniformity, scatterer distribution, body motion effects, and the like. Further, the present invention can increase the resolution of a point image with a plurality of prospective angles without performing focus calculation (delay calculation) for each target point in the subject 100. Therefore, according to the present invention, it is possible to provide an ultrasonic imaging apparatus including an adaptive beamformer that solves the above problems with a relatively small processing load.
- the device configuration including the peripheral information calculation unit and the peripheral information combining unit is shown. That is, a plurality of elements that receive ultrasonic signals from the subject, a delay unit that delays the reception signals received by the plurality of elements according to the position of a predetermined reception focus, and generates a delayed reception signal, Of the reception focal point and the peripheral region of the reception focal point, a peripheral information calculation unit that acquires information about a plurality of points from the delayed received signal, and the information acquired for each of the plurality of points are combined and combined.
- the ultrasonic imaging apparatus includes a peripheral information synthesis unit that generates a phasing output using the information and an image processing unit that generates image data using the phasing output.
- an apparatus configuration for performing aperture synthesis in the apparatus of (1) is shown. That is, the delay unit generates the delayed received signal for each of a plurality of different reception focal points, and the peripheral information calculation unit and the peripheral information synthesis unit each of the reception focal points and In the ultrasonic imaging apparatus, the information is acquired for the plurality of points in the peripheral region of the reception focal point, and the phasing output is generated.
- a configuration in which a plurality of delay units are provided in the device of (2) is also shown. That is, there are a plurality of delay units, and the delayed reception signal is generated for a reception focus that is different for each of the plurality of delay units.
- the peripheral information calculation unit and the peripheral information synthesis unit are arranged for each delay unit, acquire the information from the delayed received signal for each reception focus generated by the delay unit, and generate the phasing output To do.
- the configuration for performing aperture synthesis in the apparatus (2) or (3) is also shown. That is, the active channel setting unit is further configured to sequentially set active channels for the plurality of elements at different positions in time series, and to pass the received signals of the elements included in the active channel to the delay unit. did. At this time, a part of the plurality of reception focal points for the reception signal of the active channel at a certain time point overlaps with a part of the plurality of reception focal points for the active channel at different time points.
- the peripheral information calculation unit is configured to obtain an adaptive weight as the information by performing adaptive beamforming.
- the surrounding information calculation unit uses a steering vector that is a direction vector connecting a predetermined element of the plurality of elements and the plurality of points.
- a configuration is also shown in which the adaptive weights for the plurality of points are used.
- the peripheral information calculation unit includes a matrix calculation unit that generates a covariance matrix using the delayed received signal, and a weight vector for obtaining an adaptive weight vector for the plurality of points from the covariance matrix and the steering vector. And a calculation unit.
- the peripheral information combining unit includes a weight combining unit that generates a combined weight by adding adaptive weights for the plurality of points obtained by the peripheral information calculating unit, and an inner product of the combined weight and the delayed received signal. It is the structure containing the inner product calculating part which calculates and produces
- the configuration for weighting and adding the delayed received signal in the apparatus of (8) is also shown. That is, the inner product calculation unit is configured to generate the phasing output by multiplying the delayed received signal by the combined weight and then adding the delayed received signal.
- the surrounding information combining unit calculates an inner product of the adaptive weights for the plurality of points obtained by the surrounding information calculating unit with the delayed received signal, and outputs a pre-composition phased output for each of the plurality of adaptive weights.
- the configuration of multiplying a plurality of adaptive weights by fixed weights in the apparatus of (11) is also shown. That is, between the peripheral information calculation unit and the plurality of inner product calculation units, the adaptive weights for the plurality of points obtained by the peripheral information calculation unit are respectively multiplied by predetermined fixed weights.
- This is a configuration in which a fixed apodization multiplication unit is arranged.
- the configuration of weighting and adding the delayed received signal in the apparatus of (11) is also shown. That is, the plurality of inner product calculation units are configured by a set of a plurality of multiplication units and a plurality of addition units, and each of the plurality of multiplication units multiplies the delayed received signal by an adaptive weight for the plurality of points. Then, the plurality of adders add the delayed received signals to generate the pre-synthesis phased output.
- the active channel setting unit further includes an active channel setting unit that sequentially sets active channels at different positions in time series with respect to the plurality of elements, and delivers the reception signals of the elements included in the active channel to the delay unit.
- the peripheral information combining unit generates the final phasing output used by the image processing unit for image generation by combining the information obtained from the reception signal of the active channel by the peripheral information calculating unit. Like that.
- the active channel setting unit further includes an active channel setting unit that sequentially sets active channels at different positions in time series with respect to the plurality of elements, and delivers the reception signals of the elements included in the active channel to the delay unit.
- the peripheral information synthesizing unit adds the pre-synthesis phasing output for each of the plurality of adaptive weights obtained from the reception signal of the active channel by the peripheral information calculating unit, and the image processing unit performs image generation. The final phased output to be used is generated for each active channel.
- the adder unit is configured to generate the phasing output by acquiring and adding the pre-synthesis phasing outputs that differ in the time direction from the storage unit and the phasing unit.
- the active channel setting unit further includes an active channel setting unit that sequentially sets active channels at different positions in time series with respect to the plurality of elements, and delivers the reception signals of the elements included in the active channel to the delay unit.
- the phasing unit generates a pre-composition phasing output for each of the plurality of adaptive weights for each of the active channels, and the addition unit generates the different active channels from the storage unit and the phasing unit. The obtained pre-synthesis phased output is acquired and added to generate the phased output.
- the plurality of elements are two-dimensional arrays arranged in a two-dimensional manner, and the peripheral information calculation unit is set in the two-dimensional direction.
- the pre-synthesis phased output for the plurality of points is generated, and the storage unit stores the pre-synthesis phased output.
- the configuration of the apparatus of (6) above has a configuration including an adjustment unit that adjusts at least one of the number and direction of the steering vectors.
- the image processing unit further includes a memory unit that preliminarily stores steering angle spread angle distributions to the plurality of points according to the position of the reception focus in the image generated by the image processing unit, and the adjustment unit includes: The adaptive weight is obtained by using a steering vector of the divergence angle read from the memory unit according to the position of the reception focus.
- the relationship between the steering angle distribution and the transmission focal point is shown in the apparatus (21). That is, the divergence angle of the steering vector stored in the memory unit is set to be small at the reception focus near the transmission focus of the ultrasonic signal to the subject and large at the reception focus away from the transmission focus. .
- the relationship between the steering angle distribution and the image edge is shown in the apparatus (21) or (22). That is, the divergence angle of the steering vector stored in the memory unit is set smaller at the edge portion of the image than at the central portion of the image.
- the adjustment unit converts the brightness or signal intensity distribution of the B-mode image of the subject imaged in advance into a distribution of the steering vector divergence angle using a predetermined function. Then, the distribution of the spread angle of the obtained steering vector is used.
- the adjustment unit has previously determined a distribution of luminance or signal intensity of the image obtained by setting the steering vector using the distribution of the divergence angle of the steering vector stored in the memory.
- a function is used to convert the steering vector into a spread angle distribution. Then, the difference between the distribution of the steering vector divergence angle in the memory and the distribution of the steering vector divergence angle obtained using the function, or the distribution of the steering vector divergence angle in the memory or the Correct image brightness or signal strength.
- the apparatus described above further includes an image display unit and an operation unit for displaying an image, and accepts setting of a steering vector for the adjustment unit from an operator in at least one of the image display unit and the operation unit. It is good also as a structure that the reception part is arrange
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Acoustics & Sound (AREA)
- Molecular Biology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Multimedia (AREA)
- Gynecology & Obstetrics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Human Computer Interaction (AREA)
Abstract
Description
第1の実施形態の超音波撮像装置は、被検体からの超音波信号を受信する複数の素子と、所定の受信焦点の位置に応じて複数の素子が受信した受信信号をそれぞれ遅延させ遅延後受信信号を生成する遅延部と、周辺情報演算部と、周辺情報合成部と、画像処理部とを有する。周辺情報演算部は、受信焦点および受信焦点の周辺領域のうち、複数の点についての情報を遅延後受信信号からそれぞれ取得する。周辺情報合成部は、複数の点についてそれぞれ取得された情報を合成し、合成した情報を用いて整相出力を生成する。画像処理部は、整相出力を用いて画像データを生成する。以下、具体的に説明する。
第2の実施形態では、第1の実施形態の超音波撮像装置において、周辺情報演算部205が、適応ビームフォーミングを行うことにより、情報として適応重みを求める。
第3の実施形態では、第2の実施形態とは異なり、周辺情報合成部206は、周辺情報演算部205が求めた複数の点についての適応重みをそれぞれ遅延後受信信号に乗算し、複数の適応重みごとの合成前整相出力を生成する複数の内積演算部と、出力合成部とを有する。出力合成部は、複数の適応重みごとの合成前整相出力を加算して合成し、画像データの生成に用いる整相出力を生成する。
第4の実施形態では、遅延部(遅延回路204)は、異なる複数の受信焦点について、それぞれ遅延後受信信号を生成する構成とする周辺情報演算部205および周辺情報合成部206は、複数の受信焦点ごとに、受信焦点およびその受信焦点の周辺領域のうちの複数の点について情報を取得し、整相出力を生成する。例えば、遅延部(遅延回路204)は、複数であり、複数の遅延部ごとに異なる受信焦点について遅延後受信信号を生成する構成とする。周辺情報演算部205および周辺情報合成部206は、遅延部ごとに配置され、遅延部が生成した受信焦点ごとの遅延後受信信号から情報を取得し、整相出力を生成する構成とする。
第5の実施形態では、第3の実施形態の複数の内積演算部と出力合成部との間に、それぞれの内積演算部で生成された、複数の適応重みごとの合成前整相出力を格納する格納部(一時格納部800)を配置する。
図11を用いて、本発明の第5の実施形態の超音波撮像装置について具体的に説明する。第5の実施形態は、第3の実施形態と同様の構成であるが、周辺情報合成部206の複数の内積演算部307と出力合成部500との間に、一時格納部800が配置されている点が第3の実施形態とは異なる。第3の実施形態においては、図7のように、整相すべき信号は、アクティブチャンネルでの受信から画像処理部109まで逐次、時間とともに流れている。そのため、出力合成部500における合成処理は、現在整相を行っているラスタ内で、かつ、現スナップショットnのデータについて複数ステアリングベクトルについて複数の内積演算部307が求めたそれぞれの整相出力y1(n)、y2(n)、y3(n)を合成する処理に限られている。
第6の実施形態では出力合成部500は、格納部(一時格納部800)および複数の内積演算部から、時間方向に異なる合成前整相出力を取得して加算し、整相出力を生成する。
第7の実施形態では、複数の内積演算部は、アクティブチャンネルのそれぞれについて、複数の適応重みごとの合成前整相出力を生成し、出力合成部500は、格納部(一時格納部800)および複数の内積演算部から、異なるアクティブチャンネルについて生成された合成前整相出力を取得して加算し、整相出力を生成する。
図16、図17を用いて、本発明の第8の実施形態の超音波撮像装置について説明する。図16のように、第8の実施形態の超音波診断装置は、一時格納部800内に、(m+1)個のメモリ部900-1、900-2・・・900-(m+1)を配置している。各メモリ部900-1~900-(m+1)は、第7の実施形態のメモリ部900-1等と同様の構成である。これにより、同一ラスタ内の複数ステアリング角方向1~P、スナップショット時刻1~N、ラスタ移動方向1~m+1の3次元について、すべての整相後データから任意の組み合わせで1以上選択して、出力合成部500に出力し、合成することができる。すなわち、P×N×(m+1)個の整相後データの中から合成すべき出力1106,1107,1108を選んで出力合成部500に出力し、合成することができる。なお、第7の実施形態と同様の構成については説明を省略する。
第9の実施形態では、複数の素子106は、2次元に配列された2次元アレイである場合について説明する。周辺情報演算部205は、2次元方向に設定された複数の点についての合成前整相出力を生成する。格納部(一時格納部)800は、これら合成前整相出力を格納する。
第10の実施形態では、ステアリングベクトルの本数および方向のうちの少なくとも一方を調整する調整部501について説明する。図20は、図11の調整部501によるステアリングベクトルの設定例を示している。ステアリングベクトルの方向は、あらかじめ固定の角度を用いてもよいし、調整部501の制御により、変更される構成としてもよい。
第11の実施形態では、画像処理部で生成される画像内の受信焦点の位置に応じて、複数の点へのステアリングベクトルの拡がり角の分布を予め記憶するメモリ部をさらに有する構成とする。調整部は、受信焦点の位置に応じてメモリ部から読み出した拡がり角のステアリングベクトルを用いて、適応重みを求める。メモリ部に格納されたステアリングベクトルの拡がり角は、被検体への超音波信号の送信時焦点に近い受信焦点では小さく、送信時焦点から離れた受信焦点では大きく設定されている。しかも、ステアリングベクトルの拡がり角は、画像のエッジ部分では、画像の中央部よりも小さく設定されている。
第12の実施形態の超音波診断装置を図23(a),(b)を用いて説明する。第11の実施形態ではステアリング角度の分布を予め定めておく構成であったが、第12の実施形態ではBモード画像の信号強度または輝度分布から求める。求めたステアリング角度分布を用いて受信ビームフォーミングを行う。すなわち、調整部は、メモリに格納された前記ステアリングベクトルの拡がり角の分布を用いてステアリングベクトルを設定して得られた画像の輝度または信号強度の分布を、予め定めておいた関数を用いて、ステアリングベクトルの拡がり角の分布に変換する。メモリ内のステアリングベクトルの拡がり角の分布と、関数を用いて得たステアリングベクトルの拡がり角の分布との差分により、メモリ内のステアリングベクトルの拡がり角の分布または前記画像の輝度または信号強度を補正する。
図24に本発明の超音波撮像装置のコンソール110の一例を示す。上述してきた実施形態を実現するため、超音波診断装置のコンソール110に、ステアリングベクトルの本数や密度、合成すべき整相結果の角度方向、時間方向、ラスタ方向の個数などを変更するための操作部として、目盛りをともなったつまみ部1301,1302等を配置することができる。また、上述の各実施形態で説明したステアリング角のモードの切り替え、整相結果の合成を行うか否かの切り替え部などを行うための操作部として、スイッチ部1303をコンソール110に配置できる。これによって、操作者は実際の超音波画像を見ながら、整相合成処理の各種パラメータを変更し、被検体100ごとに最適な条件で診断を行うことができる。また、設定した値を画像表示部103の表示領域の一部1304に表示する構成にすることもできる。
図26に、本発明の実施形態の効果を示すために、シミュレーションにより得た超音波画像(Bモード像)を示す。これらの画像は、被検体100として6つの点散乱体を想定し、その超音波画像を計算機によって求めたものである。図26(a)の画像2300と画像2303は、超音波ビームに波面歪が生じていない場合である。画像2300は、従来の遅延加算方式の受信ビームフォーミングにより得た画像であり、画像2303は、従来の適応ビームフォーミングにより得た画像である。
一方、図26(b)の画像2301、2302、2304、2305は、超音波ビームに波面歪が生じている場合である。画像2301は、従来の遅延加算方式の受信ビームフォーミングにより得た画像であり、画像2304は、従来の適応ビームフォーミングにより得た画像である。画像2302は、遅延加算方式であるが、第2もしくは第3の実施形態の手法と同様の手法にステアリングベクトルの合成を行った得た画像である。画像2305は、第2もしくは第3の実施形態の手法により、適応ビームフォーミングとステアリングベクトルの合成を行って得た画像である。
101 超音波探触子
102 超音波撮像装置
103 画像表示部
104 送信ビームフォーマ
106 アレイ状の素子(超音波振動子)
107 送受信分離回路(T/R)
108 受信ビームフォーマ
109 画像処理部
110 コンソール
111 制御部
201、202、203 アクティブチャンネル
204 遅延回路
205 周辺情報演算部
206 周辺情報合成部
207 遅延後受信データのバイパス線
300 行列演算部
301 適応重み演算部
305 固定アポダイゼーション乗算部
306 重み合成部
307 内積演算部
3071 乗算部
3072 加算部
308 次元圧縮部
401、402、403 空間内の一点
404 遅延前の受信信号
405、406、407 凹面フォーカス
408、409、410 波面のそろった受信信号(遅延後受信データ)
500 出力合成部
501 調整部
502 角度調整部
503 出力調整部
601、603 中心軸近傍の異なる軸上の点
602 中心軸上の点
604、605,606 遅延凹面フォーカス
800 一時格納部
900 メモリ部
901 整相後データのバイパス線
902 出力バス
9001、9002、9003 N個のアドレスを持ったメモリ領域
1301、1302 目盛りを伴ったつまみ部
1303 スイッチ部
1304 表示領域の一部
1401 関心領域(ROI)
1402 適応的処理を施した画像
1403 用手的操作部
1501、1502、1503 ラスタ
1504、1505、1506 整相出力
1600 アクティブチャンネルの中心軸
1602、1604、1606 ステアリング角度
2071、2072、2073 遅延後受信データのバイパス線
2201 受信フォーカス
2202 遅延加算方式の点拡がり関数
2204 適応ビームフォーマの点拡がり関数
2205 焦点のデータが存在する楕円領域
3021、3022、3023 重みベクトル演算部
Claims (15)
- 被検体からの超音波信号を受信する複数の素子と、所定の受信焦点の位置に応じて前記複数の素子が受信した受信信号をそれぞれ遅延させ遅延後受信信号を生成する遅延部と、前記受信焦点および受信焦点の周辺領域のうち、複数の点についての情報を前記遅延後受信信号からそれぞれ取得する周辺情報演算部と、前記複数の点についてそれぞれ取得された前記情報を合成し、合成した前記情報を用いて整相出力を生成する周辺情報合成部と、前記整相出力を用いて画像データを生成する画像処理部とを有することを特徴とする超音波撮像装置。
- 請求項1に記載の超音波撮像装置において、前記遅延部は、異なる複数の受信焦点について、それぞれ前記遅延後受信信号を生成し、
前記周辺情報演算部および前記周辺情報合成部は、複数の前記受信焦点ごとに、前記受信焦点およびその受信焦点の周辺領域のうちの前記複数の点について前記情報を取得し、前記整相出力を生成することを特徴とする超音波撮像装置。 - 請求項2に記載の超音波撮像装置において、前記複数の素子に対してアクティブチャンネルを時系列に異なる位置に順次設定し、前記アクティブチャンネルに含まれる前記素子の受信信号を前記遅延部に受け渡すアクティブチャンネル設定部をさらに有し、
ある時点の前記アクティブチャンネルの前記受信信号についての複数の前記受信焦点の一部は、異なる時点の前記アクティブチャンネルについて複数の前記受信焦点の一部と、位置が重なっていることを特徴とする超音波撮像装置。 - 請求項1ないし3のいずれか1項に記載の超音波撮像装置において、前記周辺情報演算部は、適応ビームフォーミングを行うことにより、前記情報として適応重みを求めることを特徴とする超音波撮像装置。
- 請求項4に記載の超音波撮像装置において、
前記周辺情報演算部は、前記複数の素子のうちの所定の素子と前記複数の点とを結ぶ方向ベクトルであるステアリングベクトルを用いて前記複数の点についての前記適応重みを求めることを特徴とする超音波撮像装置。 - 請求項4または5に記載の超音波撮像装置において、前記周辺情報合成部は、前記周辺情報演算部が求めた前記複数の点についての適応重みを加算して合成重みを生成する重み合成部と、前記合成重みと前記遅延後受信信号を内積演算して前記整相出力を生成する内積演算部を含むことを特徴とする超音波撮像装置。
- 請求項4または5に記載の超音波撮像装置において、前記周辺情報合成部は、前記周辺情報演算部が求めた前記複数の点についての前記適応重みをそれぞれ前記遅延後受信信号と内積演算し、複数の前記適応重みごとの合成前整相出力を生成する複数の内積演算部と、前記複数の適応重みごとの合成前整相出力を加算して合成し、前記画像データの生成に用いる前記整相出力を生成する出力合成部とを含むことを特徴とする超音波撮像装置。
- 請求項1ないし7のいずれか1項に記載の超音波撮像装置において、前記複数の素子に対してアクティブチャンネルを時系列に異なる位置に順次設定し、前記アクティブチャンネルに含まれる前記素子の受信信号を前記遅延部に受け渡すアクティブチャンネル設定部をさらに有し、
前記周辺情報合成部は、前記周辺情報演算部が一つの前記アクティブチャンネルの前記受信信号から得た前記情報を合成して前記画像処理部が画像生成に用いる最終的な前記整相出力を生成することを特徴とする超音波撮像装置。 - 請求項6に記載の超音波撮像装置において、前記複数の素子に対してアクティブチャンネルを時系列に異なる位置に順次設定し、前記アクティブチャンネルに含まれる前記素子の受信信号を前記遅延部に受け渡すアクティブチャンネル設定部をさらに有し、
前記周辺情報合成部は、前記周辺情報演算部が一つの前記アクティブチャンネルの前記受信信号から得た、前記複数の適応重みごとの合成前整相出力を加算し、前記画像処理部が画像生成に用いる最終的な前記整相出力を前記アクティブチャンネルごとに生成することを特徴とする超音波撮像装置。 - 請求項7に記載の超音波撮像装置において、前記整相部と前記加算部との間には、前記整相部で生成された、前記複数の適応重みごとの前記合成前整相出力を格納する格納部が配置されていることを特徴とする超音波撮像装置。
- 請求項10に記載の超音波撮像装置において、前記複数の素子に対してアクティブチャンネルを時系列に異なる位置に順次設定し、前記アクティブチャンネルに含まれる前記素子の受信信号を前記遅延部に受け渡すアクティブチャンネル設定部をさらに有し、
前記整相部は、それぞれの前記アクティブチャンネルについて、複数の前記適応重みごとの合成前整相出力を生成し、
前記加算部は、前記格納部および前記整相部から、異なる前記アクティブチャンネルについて生成された前記合成前整相出力を取得して加算し、前記整相出力を生成することを特徴とする超音波撮像装置。 - 請求項10に記載の超音波撮像装置において、前記複数の素子は、2次元に配列された2次元アレイであり、前記周辺情報演算部は、前記2次元方向に設定された前記複数の点についての前記合成前整相出力を生成し、前記格納部は、前記合成前整相出力を格納することを特徴とする超音波撮像装置。
- 請求項5に記載の超音波撮像装置において、前記ステアリングベクトルの本数および方向のうちの少なくとも一方を調整する調整部を有することを特徴とする超音波撮像装置。
- 請求項13に記載の超音波撮像装置において、前記画像処理部で生成される画像内の前記受信焦点の位置に応じて、前記複数の点へのステアリングベクトルの拡がり角の分布を予め記憶するメモリ部をさらに有し、
前記調整部は、前記受信焦点の位置に応じて前記メモリ部から読み出した拡がり角のステアリングベクトルを用いて、前記適応重みを求めることを特徴とする超音波撮像装置。 - 請求項13に記載の超音波撮像装置において、前記画像を表示する画像表示部および操作部をさらに有し、
前記画像表示部および操作部の少なくとも一方に、前記調整部へのステアリングベクトルの設定を操作者から受け付ける受け付け部が配置されていることを特徴とする超音波撮像装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201380009474.6A CN104114099B (zh) | 2012-02-15 | 2013-01-23 | 超声波摄像装置 |
JP2014500136A JP5913557B2 (ja) | 2012-02-15 | 2013-01-23 | 超音波撮像装置 |
US14/378,507 US20150025385A1 (en) | 2012-02-15 | 2013-01-23 | Ultrasonic imaging device |
EP13749172.6A EP2815701A4 (en) | 2012-02-15 | 2013-01-23 | ULTRASONIC IMAGING DEVICE |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012-030991 | 2012-02-15 | ||
JP2012030991 | 2012-02-15 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013121842A1 true WO2013121842A1 (ja) | 2013-08-22 |
Family
ID=48983967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2013/051305 WO2013121842A1 (ja) | 2012-02-15 | 2013-01-23 | 超音波撮像装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20150025385A1 (ja) |
EP (1) | EP2815701A4 (ja) |
JP (1) | JP5913557B2 (ja) |
CN (1) | CN104114099B (ja) |
WO (1) | WO2013121842A1 (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104414687A (zh) * | 2013-09-10 | 2015-03-18 | 精工爱普生株式会社 | 超声波测定装置、超声波图像装置及超声波测定方法 |
JP2015071028A (ja) * | 2013-09-05 | 2015-04-16 | セイコーエプソン株式会社 | 超音波測定装置、超音波画像装置及び超音波測定方法 |
JP2015217172A (ja) * | 2014-05-19 | 2015-12-07 | 株式会社東芝 | 超音波診断装置 |
WO2016072541A1 (ko) * | 2014-11-07 | 2016-05-12 | 삼성전자주식회사 | 초음파 영상 장치 및 초음파 영상 장치의 제어 방법 |
JP2016086875A (ja) * | 2014-10-30 | 2016-05-23 | セイコーエプソン株式会社 | 超音波測定装置、超音波診断装置及び超音波測定方法 |
WO2016132924A1 (ja) * | 2015-02-18 | 2016-08-25 | 日立アロカメディカル株式会社 | 超音波撮像装置および超音波信号の処理方法 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3154435B1 (en) | 2014-06-12 | 2018-08-15 | Koninklijke Philips N.V. | Ultrasound transducer assembly |
CN107205730B (zh) * | 2015-02-12 | 2020-02-07 | 株式会社日立制作所 | 超声波成像装置、发送间权重调整方法和超声波成像方法 |
JP2017164408A (ja) * | 2016-03-18 | 2017-09-21 | セイコーエプソン株式会社 | 画像生成装置および画像生成方法 |
JP6747108B2 (ja) * | 2016-07-05 | 2020-08-26 | コニカミノルタ株式会社 | 超音波信号処理装置、超音波信号処理方法、及び、超音波診断装置 |
JP6819184B2 (ja) * | 2016-09-26 | 2021-01-27 | セイコーエプソン株式会社 | 超音波測定装置および超音波測定装置の制御方法 |
CN109843180A (zh) * | 2016-10-09 | 2019-06-04 | 柯惠有限合伙公司 | 用于驱动超声成像换能器的系统和方法 |
US11719797B2 (en) * | 2017-05-11 | 2023-08-08 | Koninklijke Philips N.V. | Methods and systems for controlling the generation of a compound ultrasound image |
CN108324324A (zh) * | 2018-03-12 | 2018-07-27 | 西安交通大学 | 一种超声低频经颅容积超分辨率三维造影成像方法及系统 |
JP7008549B2 (ja) * | 2018-03-16 | 2022-01-25 | 富士フイルムヘルスケア株式会社 | 超音波診断装置 |
CN109223035B (zh) * | 2018-08-21 | 2021-09-28 | 青岛海信医疗设备股份有限公司 | 超声信号处理方法、装置、设备及存储介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5249088A (en) | 1975-09-08 | 1977-04-19 | Stahli Christian | Pipetter and pippetting method |
JPH07303640A (ja) | 1994-03-16 | 1995-11-21 | Fujitsu Ltd | 超音波診断装置 |
JP2002336249A (ja) | 2000-06-10 | 2002-11-26 | Medison Co Ltd | 多段構造の遅延素子を用いる超音波受信ビーム成形装置 |
JP2003079623A (ja) * | 2001-09-05 | 2003-03-18 | Medison Co Ltd | 超音波撮像システム及び超音波撮像方法 |
JP2003225237A (ja) * | 2002-02-05 | 2003-08-12 | Hitachi Medical Corp | 超音波撮像断装置及び撮像方法 |
JP2004049926A (ja) * | 2002-07-23 | 2004-02-19 | Medison Co Ltd | デジタル受信集束装置 |
JP2008212492A (ja) * | 2007-03-06 | 2008-09-18 | Hitachi Medical Corp | 超音波診断装置 |
JP2010063875A (ja) | 2008-08-11 | 2010-03-25 | Canon Inc | 遅延調整モジュールおよび超音波受信ビーム成形装置 |
JP2010082371A (ja) | 2008-10-02 | 2010-04-15 | Canon Inc | 超音波受信ビーム成形装置 |
JP2011045708A (ja) * | 2009-07-28 | 2011-03-10 | Toshiba Corp | 超音波診断装置、超音波画像処理装置、超音波診断装置制御プログラム及び超音波画像処理プログラム |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2290240C (en) * | 1999-11-24 | 2008-03-11 | Stergios Stergiopoulos | High resolution 3d ultrasound imaging system deploying a multi-dimensional array of sensors and method for multi-dimensional beamforming sensor signals |
US7207942B2 (en) * | 2003-07-25 | 2007-04-24 | Siemens Medical Solutions Usa, Inc. | Adaptive grating lobe suppression in ultrasound imaging |
US7744532B2 (en) * | 2004-03-31 | 2010-06-29 | Siemens Medical Solutions Usa, Inc. | Coherence factor adaptive ultrasound imaging methods and systems |
JP2009089940A (ja) * | 2007-10-10 | 2009-04-30 | Toshiba Corp | 超音波診断装置 |
KR101798081B1 (ko) * | 2011-04-20 | 2017-11-15 | 삼성전자주식회사 | 진단영상을 생성하는 방법, 이를 수행하는 장치, 진단시스템 및 의료영상시스템 |
US9754185B2 (en) * | 2011-09-15 | 2017-09-05 | Hitachi, Ltd. | Ultrasound imaging apparatus |
JP6039305B2 (ja) * | 2012-08-23 | 2016-12-07 | キヤノン株式会社 | 被検体情報取得装置、情報処理装置および被検体情報取得方法 |
JP6103857B2 (ja) * | 2012-08-28 | 2017-03-29 | キヤノン株式会社 | 被検体情報取得装置、表示方法、及びプログラム |
JP6025456B2 (ja) * | 2012-08-28 | 2016-11-16 | キヤノン株式会社 | 被検体情報取得装置、表示方法、及びプログラム |
JP6253296B2 (ja) * | 2012-08-28 | 2017-12-27 | キヤノン株式会社 | 被検体情報取得装置、表示方法、及びプログラム |
JP6021520B2 (ja) * | 2012-08-28 | 2016-11-09 | キヤノン株式会社 | 被検体情報取得装置、表示方法、及びプログラム |
KR20140056597A (ko) * | 2012-10-30 | 2014-05-12 | 삼성전자주식회사 | 빔 포밍 계수 결정 방법, 빔 포밍 방법 및 초음파 이미징 장치 |
US20150351720A1 (en) * | 2013-01-11 | 2015-12-10 | Hitachi Aloka Medical, Ltd. | Ultrasonic imaging device |
KR102185415B1 (ko) * | 2013-01-11 | 2020-12-02 | 삼성전자주식회사 | 빔 포밍 모듈, 상기 빔 포밍 모듈을 이용하는 초음파 이미징 장치, 상기 빔 포밍 모듈을 이용한 빔 포밍 방법 및 상기 빔 포밍 모듈을 이용한 초음파 이미징 장치의 제어 방법 |
KR101832835B1 (ko) * | 2013-07-11 | 2018-02-28 | 삼성전자주식회사 | 영상 처리 모듈, 초음파 영상 장치, 영상 처리 방법 및 초음파 영상 장치의 제어 방법 |
KR20150008261A (ko) * | 2013-07-11 | 2015-01-22 | 삼성전자주식회사 | 빔 포머, 빔 포밍 방법, 상기 초음파 영상 장치 및 초음파 영상 장치의 제어 방법 |
-
2013
- 2013-01-23 US US14/378,507 patent/US20150025385A1/en not_active Abandoned
- 2013-01-23 EP EP13749172.6A patent/EP2815701A4/en not_active Withdrawn
- 2013-01-23 JP JP2014500136A patent/JP5913557B2/ja active Active
- 2013-01-23 CN CN201380009474.6A patent/CN104114099B/zh not_active Expired - Fee Related
- 2013-01-23 WO PCT/JP2013/051305 patent/WO2013121842A1/ja active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5249088A (en) | 1975-09-08 | 1977-04-19 | Stahli Christian | Pipetter and pippetting method |
JPH07303640A (ja) | 1994-03-16 | 1995-11-21 | Fujitsu Ltd | 超音波診断装置 |
JP2002336249A (ja) | 2000-06-10 | 2002-11-26 | Medison Co Ltd | 多段構造の遅延素子を用いる超音波受信ビーム成形装置 |
JP2003079623A (ja) * | 2001-09-05 | 2003-03-18 | Medison Co Ltd | 超音波撮像システム及び超音波撮像方法 |
JP2003225237A (ja) * | 2002-02-05 | 2003-08-12 | Hitachi Medical Corp | 超音波撮像断装置及び撮像方法 |
JP2004049926A (ja) * | 2002-07-23 | 2004-02-19 | Medison Co Ltd | デジタル受信集束装置 |
JP2008212492A (ja) * | 2007-03-06 | 2008-09-18 | Hitachi Medical Corp | 超音波診断装置 |
JP2010063875A (ja) | 2008-08-11 | 2010-03-25 | Canon Inc | 遅延調整モジュールおよび超音波受信ビーム成形装置 |
JP2010082371A (ja) | 2008-10-02 | 2010-04-15 | Canon Inc | 超音波受信ビーム成形装置 |
JP2011045708A (ja) * | 2009-07-28 | 2011-03-10 | Toshiba Corp | 超音波診断装置、超音波画像処理装置、超音波診断装置制御プログラム及び超音波画像処理プログラム |
Non-Patent Citations (6)
Title |
---|
F. VIGNON ET AL.: "Capon beamforming in medical ultrasound imaging with focused beams", IEEE TRANS. ULTRASON. FERROELECTR. FREQ. CONTROL, vol. 55, no. 3, 2008, pages 619 - 628, XP055091879, DOI: doi:10.1109/TUFFC.2008.686 |
I. K. HOLFORT ET AL.: "Adaptive receive and transmit apodization for synthetic aperture ultrasound imaging", PROC. IEEE ULTRASON. SYMP., 2009, pages 1 - 4, XP031654941 |
I. K. HOLFORT ET AL.: "Broadband minimum variance beamforming for ultrasound imaging", IEEE TRANS. ULTRASON. FERROELECTR. FREQ. CONTROL., vol. 56, no. 2, 2009, pages 314 - 325, XP011252292 |
J. CAPON: "High-resolution frequency wavenumber spectrum analysis", PROC. IEEE, vol. 57, August 1969 (1969-08-01), pages 1408 - 1418, XP008025691 |
J. F. SYNNEVAG ET AL.: "Benefits of Minimum-Variance Beamforming in Medical Ultrasound Imaging", IEEE TRANS. ULTRASON, FERROELECTR. FREQ. CONTROL, vol. 56, no. 9, 2009, pages 1868 - 1879, XP011283394, DOI: doi:10.1109/TUFFC.2009.1263 |
Z. WANG ET AL.: "Time-delay- and time-reversal-based robust capon beamformers for ultrasound imaging", IEEE TRANS. MED. IMAG., vol. 24, no. 10, 2005, pages 1308 - 1322, XP055098581, DOI: doi:10.1109/TMI.2005.857222 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015071028A (ja) * | 2013-09-05 | 2015-04-16 | セイコーエプソン株式会社 | 超音波測定装置、超音波画像装置及び超音波測定方法 |
CN104414687A (zh) * | 2013-09-10 | 2015-03-18 | 精工爱普生株式会社 | 超声波测定装置、超声波图像装置及超声波测定方法 |
JP2015077393A (ja) * | 2013-09-10 | 2015-04-23 | セイコーエプソン株式会社 | 超音波測定装置、超音波画像装置、及び超音波測定方法 |
JP2015217172A (ja) * | 2014-05-19 | 2015-12-07 | 株式会社東芝 | 超音波診断装置 |
JP2016086875A (ja) * | 2014-10-30 | 2016-05-23 | セイコーエプソン株式会社 | 超音波測定装置、超音波診断装置及び超音波測定方法 |
WO2016072541A1 (ko) * | 2014-11-07 | 2016-05-12 | 삼성전자주식회사 | 초음파 영상 장치 및 초음파 영상 장치의 제어 방법 |
WO2016132924A1 (ja) * | 2015-02-18 | 2016-08-25 | 日立アロカメディカル株式会社 | 超音波撮像装置および超音波信号の処理方法 |
CN107205729A (zh) * | 2015-02-18 | 2017-09-26 | 株式会社日立制作所 | 超声波成像装置以及超声波信号的处理方法 |
JPWO2016132924A1 (ja) * | 2015-02-18 | 2017-11-16 | 株式会社日立製作所 | 超音波撮像装置および超音波信号の処理方法 |
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 |
Also Published As
Publication number | Publication date |
---|---|
EP2815701A4 (en) | 2015-04-22 |
EP2815701A1 (en) | 2014-12-24 |
JPWO2013121842A1 (ja) | 2015-05-11 |
US20150025385A1 (en) | 2015-01-22 |
CN104114099B (zh) | 2016-02-17 |
CN104114099A (zh) | 2014-10-22 |
JP5913557B2 (ja) | 2016-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5913557B2 (ja) | 超音波撮像装置 | |
US6685641B2 (en) | Plane wave scanning reception and receiver | |
JP5813776B2 (ja) | 超音波撮像装置 | |
US7789831B2 (en) | Synthetic elevation aperture for ultrasound systems and methods | |
JP6189867B2 (ja) | 超音波撮像装置 | |
EP2534502B1 (en) | Ultrasound imaging system using beamforming techniques for phase coherence grating lobe suppression | |
KR101942595B1 (ko) | 영상 획득속도 최적화를 구비한 영상장치 | |
JP4130114B2 (ja) | 超音波イメージング装置及び超音波信号処理方法 | |
CN110101411B (zh) | 超声成像空间复合方法及系统 | |
WO2012049124A2 (en) | Methods and systems for producing compounded ultrasound images | |
CN109363714B (zh) | 一种超声成像设备及其超声成像方法 | |
JP6352050B2 (ja) | 超音波診断装置 | |
US20070083109A1 (en) | Adaptive line synthesis for ultrasound | |
CN107569254B (zh) | 超声波信号处理装置、超声波信号处理方法以及超声波诊断装置 | |
WO2015025654A1 (ja) | 超音波撮像装置 | |
WO2019202593A1 (en) | Sparse convolutional beamforming for ultrasound imaging | |
JP6863817B2 (ja) | 超音波撮像装置 | |
Bottenus et al. | Implementation of swept synthetic aperture imaging | |
Wang et al. | Generalized sidelobe canceller beamforming method for ultrasound imaging | |
JP7387249B2 (ja) | 超音波診断装置、医用画像処理装置及び医用画像処理プログラム | |
Matrone et al. | Spatial Coherence Based Beamforming in Multi-Line Transmit Echocardiography | |
JP6469746B2 (ja) | 超音波診断装置及び制御プログラム | |
Soozande et al. | Virtually Extended Array Imaging Improves Lateral Resolution in High Frame Rate Volumetric Imaging | |
JP6143541B2 (ja) | 超音波診断装置及び制御プログラム | |
Xu et al. | Transmit aperture function for large depth focusing combined with phase coherence imaging for interference suppressing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 201380009474.6 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 13749172 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2014500136 Country of ref document: JP Kind code of ref document: A |
|
REEP | Request for entry into the european phase |
Ref document number: 2013749172 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 14378507 Country of ref document: US Ref document number: 2013749172 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |