CN110101409A - Beam synthesizing method, ultrasonic imaging method, device and equipment - Google Patents
Beam synthesizing method, ultrasonic imaging method, device and equipment Download PDFInfo
- Publication number
- CN110101409A CN110101409A CN201910204036.4A CN201910204036A CN110101409A CN 110101409 A CN110101409 A CN 110101409A CN 201910204036 A CN201910204036 A CN 201910204036A CN 110101409 A CN110101409 A CN 110101409A
- Authority
- CN
- China
- Prior art keywords
- signal
- depth
- wave beam
- collection
- mentioned
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/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
-
- 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
Abstract
This application provides a kind of beam synthesizing method, ultrasonic imaging method, device and equipment, include the following steps: to calibrate the corresponding depth of focus value Z of every wave beam by given step, and the depth of focus value Z2 after calibration is replaced into the corresponding original depth of focus value Z of wave beam；Signal alignment is carried out according to the depth of focus value Z2 beam collection after calibration, and obtains the BF1 signal collection in the beam collection after signal alignment, then each BF1 signal that BF1 signal is concentrated carries out specified superposition and obtains BF2 signal.By realizing the signal alignment of multi-beam, can preferably suppressed sidelobes, improve image resolution ratio and signal-to-noise ratio；It is superimposed by multi-beam and greatly improves line density, improve the spatial resolution of image；Device keeps data after Beam synthesis more accurate by increasing calibration module and synthesis module in Beam synthesis module；And is judged by sound field, to effective sound field data investigation, realize preferable focusing effect compared to conventional method, improve picture quality.
Description
Technical field
This application involves field of computer technology, more particularly to a kind of beam synthesizing method, ultrasonic imaging method, device
And equipment.
Background technique
Ultrasonic imaging is widely used in clinical doctor because it is with the advantages such as safe, real-time, portable, noninvasive and at low cost
Learn diagnosis.And Beam synthesis is to play conclusive effect to image quality in core position in ultrasonic image-forming system.Wave beam
Main lobe width and sidelobe magnitudes are used to judge the height of formed beam quality in synthesis, and the width of usual main lobe is narrower, then
The lateral resolution of imaging is higher；The amplitude of secondary lobe is smaller, then the contrast being imaged is bigger, and artifact noise is fewer.
In addition Wave beam forming reception is divided into simple beam and multi-beam receives, in locomotor real-time diagnosis, multi-beam
Reception can greatly improve frame per second, however phase is misaligned between the multi-beam that transmitting receives every time, and superposition in this way obtains
RF signal artifact noise equally with higher.
Conventional method is by carrying out delay superposition wave beam, and poor image quality, secondary lobe is higher ranked, and spatial resolution is low, though
So beam main lobe width and sidelobe magnitudes can be controlled by three kinds of dynamic focusing, amplitude apodization and dynamic aperture methods, from
And suppressed sidelobes improves the quality of ultrasonic imaging, but also inhibit main lobe simultaneously, have an impact to image resolution ratio；And it is existing its
His method is also to improve resolution ratio to reduce sidelobe magnitudes, and improving signal-to-noise ratio etc. improves focusing effect, but improves frame per second,
The method for improving spatial resolution is still existing technology blind area.
Summary of the invention
In view of the above problems, it proposes the embodiment of the present invention and overcomes the above problem or at least partly in order to provide one kind
The method and device of the Beam synthesis to solve the above problems.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of beam synthesizing method, include the following steps:
The corresponding depth of focus value Z of every wave beam is calibrated by given step, and the depth of focus value Z2 after calibration is replaced
Commutation corresponds to the original depth of focus value Z of wave beam；
Signal alignment is carried out according to the depth of focus value Z2 beam collection after above-mentioned calibration, and obtains the wave beam after signal alignment
The BF1 signal collection of concentration, then each BF1 signal that above-mentioned BF1 signal is concentrated is subjected to specified superposition and obtains BF2 signal.
Further, the above-mentioned depth of focus value Z2 beam collection according to after above-mentioned calibration carries out signal alignment, and obtains letter
The BF1 signal collection in beam collection after number alignment, then each BF1 signal that above-mentioned BF1 signal is concentrated is subjected to specified superposition and is obtained
Further include following steps before the step of BF2 signal:
Useful signal filtering is carried out to above-mentioned BF1 signal collection, and separator goes out effective BF1 that above-mentioned BF1 signal is concentrated
Signal and invalid BF1 signal.
Further, the corresponding depth of focus value Z of every wave beam is calibrated above by given step, and will be poly- after calibration
Depth of focus angle value Z2 replaces the step of corresponding wave beam original depth of focus value Z, includes the following steps:
It is changed by array number N_elements, adjacent array element spacing pitch, transmitting line number nlines and moving step length step
Calculate the spacing dx of each wave beam line in linear array；
Pass through sample on the spacing dx of wave beam line each in linear array, every time transmitting or the number of beams beam received, emission lines
The distance z_sample1 and focal length F of point to launch point are conversed respectively on the sample point in the emission lines of center to each wave beam line
Virtual focus point lateral distance Dxibeam and fore-and-aft distance Dyibeam；
According to the lateral distance Dxibeam of the virtual focus point on the sample point in the emission lines of center to each wave beam line and indulge
To distance Dyibeam converse the sample point in the emission lines of center to each wave beam line virtual focus point distance dr1ibeam；
It is conversed according to the distance dr1ibeam and focal length F of the sample point in the emission lines of center to each wave beam line virtual focus point
Depth of focus value Z2 after the corresponding calibration of each wave beam line.
Further, above-mentioned that useful signal filtering is carried out to above-mentioned BF1 signal collection, and separator goes out above-mentioned BF1 signal
The step of effective BF1 signal and invalid BF1 signal for concentrating, include the following steps:
It is converted by the aperture A of ultrasonic transmitter, port number N_channel, adjacent array element spacing pitch and focal length F
Emit the angle, θ 1 of sound field out；
Calculate the angle angle1ibeam that each sample point on each wave beam line deviates emission lines；
Judge whether above-mentioned angle angle1ibeam is greater than θ 1/2；
If it is not, then determining that sample point corresponding with above-mentioned angle angle1ibeam is within the scope of effective sound field.
Further, the above-mentioned depth of focus value Z2 beam collection according to after above-mentioned calibration carries out signal alignment, and obtains letter
The BF1 signal collection in beam collection after number alignment, then each BF1 signal that above-mentioned BF1 signal is concentrated is subjected to specified superposition and is obtained
The step of BF2 signal, include the following steps:
Pass through the depth of focus value Z2 conversion after the calibration of coordinate ref1, the array element coordinate xe1 and each wave beam line that receive line
Distance dn1ibeam of each array element to focusing sample point out；
Depth of focus value Z2 after calibration and each battle array by receiving the coordinate ref1 of line, array element coordinate xe1, each wave beam line
Member converses the delay τ 1 of each array element to the distance dn1ibeam for focusing sample point；
Corresponding wave beam line is subjected to phase alignment by the delay τ 1 of each array element, and the wave beam after signal is aligned
The BF1 signal collection of line is overlapped to obtain above-mentioned BF2 signal.
Further, corresponding wave beam line is subjected to phase alignment above by the delay τ 1 of each array element, and will letter
The BF1 signal collection of wave beam line after number alignment is overlapped the step of obtaining above-mentioned BF2 signal, includes the following steps:
Correspondence is conversed according to the delay τ 1 of each array element and the above-mentioned angle angle1ibeam within the scope of effective sound field
Sample point echo-signal；
The echo-signal of above-mentioned corresponding sample point is collected to the BF1 signal collection after forming above-mentioned signal alignment.
Further, the corresponding depth of focus value Z of every wave beam is calibrated above by given step, and will be poly- after calibration
Depth of focus angle value Z2 replaces the step of corresponding wave beam original depth of focus value Z, includes the following steps:
Convex battle array central angle is conversed according to array number N_elements, adjacent array element spacing pitch and convex battle array center of circle radius R
beta；
Angle between receiving line is conversed according to transmitting line number nlines, moving step length step and above-mentioned convex battle array central angle beta
dx_beta；
Reception wave beam is conversed according to angle dx_beta between each number of beams beam and above-mentioned reception line for emitting or receiving
To the included angle A ngle_beam of launching beam；
It is changed respectively according to focal length F, above-mentioned convex battle array center of circle radius R and the included angle A ngle_beam for receiving wave beam to launching beam
Sample point is calculated to the lateral distance Hxibeam and fore-and-aft distance Hyibeam for receiving wave beam virtual focus point, and according to above-mentioned
Convex battle array center of circle radius R and fore-and-aft distance Hyibeam converses correction focal length FA；
According to the distance z_sample2 of sample point in emission lines to the center of circle, above-mentioned convex battle array center of circle radius R, lateral distance
Hxibeam and fore-and-aft distance Hyibeam converses the sample point in the emission lines of center to the distance of each wave beam line virtual focus point
dr2ibeam；
According to the sample point in above-mentioned center emission lines to the distance dr2ibeam of each wave beam line virtual focus point, focal length F and
Correction focal length FA converses the depth of focus value Z2 after the corresponding calibration of each wave beam line.
Further, above-mentioned that useful signal filtering is carried out to above-mentioned BF1 signal collection, and separator goes out above-mentioned BF1 signal
The step of effective BF1 signal and invalid BF1 signal for concentrating, include the following steps:
Transmitting is conversed according to convex battle array center of circle radius R, port number N_channel, focal length F and adjacent array element spacing pitch
The angle, θ 2 of sound field；
Obtain the angle angle2ibeam that sample point on each wave beam line deviates emission lines；
Judge whether above-mentioned angle angle2ibeam is greater than θ 2/2
If it is not, then determining that sample point corresponding with above-mentioned angle angle2ibeam is within the scope of effective sound field.
Further, the above-mentioned depth of focus value Z2 beam collection according to after above-mentioned calibration carries out signal alignment, and obtains letter
The BF1 signal collection in beam collection after number alignment, then each BF1 signal that above-mentioned BF1 signal is concentrated is subjected to specified superposition and is obtained
The step of BF2 signal, include the following steps:
It is changed according to the depth of focus value Z2 after the reception starting position coordinates M of line, the position coordinates N of array element, above-mentioned calibration
Distance dn2ibeam of the calculating array element to focus point；
Each array element is conversed according to the distance dn2ibeam of depth of focus value Z2 and array element to focus point after above-mentioned calibration
Delay τ 2；
Corresponding wave beam line is subjected to phase alignment by the delay τ 2 of each array element, and the wave beam after signal is aligned
The BF1 signal collection of line is overlapped to obtain above-mentioned BF2 signal.
Further, corresponding wave beam line is subjected to phase alignment above by the delay τ 2 of each array element, and will letter
The BF1 signal collection of wave beam line after number alignment is overlapped the step of obtaining above-mentioned BF2 signal, includes the following steps:
Correspondence is conversed according to the delay τ 2 of each array element and the above-mentioned angle angle2ibeam within the scope of effective sound field
Sample point echo-signal；
The echo-signal of above-mentioned corresponding sample point is collected to the BF1 signal collection after forming above-mentioned signal alignment.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of ultrasonic imaging method, include the following steps,
It obtains ultrasound data and is converted into corresponding digital signal；
Above-mentioned digital signal is synthesized into radiofrequency signal by specified correction；
Above-mentioned radiofrequency signal is isolated into carrier signal by specified signal processing；
Above-mentioned carrier signal is obtained into ultrasound image by specified image procossing；
Above-mentioned the step of above-mentioned digital signal is synthesized into radiofrequency signal by specified correction, including any of the above-described implementation
Beam synthesizing method described in example.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of beam synthesizers, including following specific module:
Calibration module, for calibrating the corresponding depth of focus value Z of every wave beam by given step, and will be poly- after calibration
Depth of focus angle value Z2 replaces the corresponding original depth of focus value Z of wave beam；
Synthesis module for carrying out signal alignment according to the depth of focus value Z2 beam collection after above-mentioned calibration, and obtains letter
The BF1 signal collection in beam collection after number alignment, then each BF1 signal that above-mentioned BF1 signal is concentrated is subjected to specified superposition and is obtained
BF2 signal.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of supersonic imaging devices, including following specific module:
Module is obtained, for obtaining ultrasound data and being converted into corresponding digital signal；
Beam synthesis module, for above-mentioned digital signal to be synthesized radiofrequency signal by specified correction；
Separation module, for above-mentioned radiofrequency signal to be isolated carrier signal by specified signal processing；
Image-forming module, for above-mentioned carrier signal to be obtained ultrasound image by specified image procossing,
Above-mentioned Beam synthesis module, including beam synthesizer described in any of the above-described embodiment.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of supersonic imaging apparatus, including generate ultrasonic wave letter
It number is radiated tested tissue and neutralizes the signal transceiver for absorbing reflected sonic signals, received sound wave is converted into electric signal
Vibrating sensor will receive signal sampling and digitized analog-digital converter (A/D), compensate the ultrasound wave amplitude as caused by depth
Be worth decaying time gain compensation device (TGC), convert digital signal to the reception beam synthesizer of RF signal, by RF signal into
Row envelope extraction and demodulation process isolate the signal processor of carrier signal, and carrier signal is scanned conversion and rear end figure
As the digital scan converter DSC handled and for the display of the image shown in finally.
Wherein, above-mentioned reception beam synthesizer, including phase correction module, the first time wave beam for synthesizing BF1 signal
Beam synthesis module, effective sound field judgment module and second of beam synthesizer for synthesizing BF2 signal, and above-mentioned phase school
Positive module, the first time multi-beam beam synthesis module for synthesizing BF1 signal, effective sound field judgment module and for synthesizing BF2
Second of beam synthesizer of signal is sequentially connected.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of computer equipment, including memory, processor with
And the computer program that can be run on a memory and on a processor is stored, the processor realizes this when executing described program
Beam synthesizing method described in any one of inventive embodiments.
Compared with prior art, the application includes following advantages:
In the embodiment of the present invention, by realizing the spacial alignment and phase alignment of multi-beam, can preferably suppressed sidelobes,
Improve image resolution ratio and signal-to-noise ratio；It is superimposed by multi-beam and greatly improves line density, improve the spatial resolution of image；
This method realizes that algorithm is simple, and multirate, high s/n ratio, high-resolution can be realized in the case where consuming a small amount of process resource
Effect；More data are more by increasing calibration module and synthesis module in Beam synthesis module, after making Beam synthesis for device
Accurately；And judged by sound field, to effective sound field data investigation, realize preferable focusing effect compared to conventional method,
Improve picture quality.
Detailed description of the invention
Fig. 1 is the step flow diagram of the beam synthesizing method of one embodiment of the invention；
Fig. 2 is the step flow diagram of the beam synthesizing method of one embodiment of the invention；
Fig. 3 is the step flow diagram of the beam synthesizing method of one embodiment of the invention；
Fig. 4 is the step flow diagram of the beam synthesizing method of one embodiment of the invention；
Fig. 5 is that the linear array depth of focus correction of one embodiment of the invention calculates schematic diagram；
Fig. 6 is that the effective sound field of linear array of one embodiment of the invention calculates schematic diagram；
Fig. 7 is that the linear array delay distance of one embodiment of the invention calculates schematic diagram；
Fig. 8 is that the convex battle array depth of focus correction of one embodiment of the invention calculates schematic diagram；
Fig. 9 is that the effective sound field of convex battle array of one embodiment of the invention calculates schematic diagram；
Figure 10 is that the convex battle array delay distance of one embodiment of the invention calculates schematic diagram；
Figure 11 is the modular structure schematic diagram of the beam synthesizer of one embodiment of the invention；
Figure 12 is a kind of structural schematic diagram of computer equipment of one embodiment of the invention.
1, calibration module；2, synthesis module；12, computer equipment；14, external equipment；16, processing unit；18, bus；
20, network adapter；22, (I/O) interface；24, display；28, system storage；30, random access memory (RAM)；32,
Cache memory；34, storage system；40, program/utility；42, program module.
The object of the invention is realized, the embodiments will be further described with reference to the accompanying drawings for functional characteristics and advantage.
Specific embodiment
In order to make the above objects, features, and advantages of the present application more apparent, with reference to the accompanying drawing and it is specific real
Applying mode, the present application will be further described in detail.
Obviously, described embodiment is only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
It in addition, the technical solution between each embodiment can be combined with each other, but must be with ordinary skill
Based on personnel can be realized, this technical side will be understood that when the combination of technical solution appearance is conflicting or cannot achieve
The combination of case is not present, also not the present invention claims protection scope within.
Finally, all the embodiments in this specification are described in a progressive manner, what each embodiment stressed
It is the difference from other embodiments, the same or similar parts between the embodiments can be referred to each other.
It should be noted that signal alignment includes spacial alignment and phase alignment in any embodiment of the present invention；Space
Alignment: refer to that reception line position is identical；Phase alignment: refer to that different reception lines (the reception line after spacial alignment) are prolonged using different
When the time, it is more identical (focusing on a bit) to enable signals to reach spatial position simultaneously, enables signal by correction
Enough gather same position.
It should be noted that for convenience of the expression of formula, in the calculating parameter mentioned in following any embodiment
“Xibeam" corresponding " X ", the X that is expressed as any wave beam line i of parameteribeam(j) any point j being expressed as in any wave beam line i is corresponding
" X ", such as: by the depth of focus value Z2 after the corresponding calibration of any wave beam line with Z2ibeamIt is indicated in formula.
Referring to Fig.1, a kind of beam synthesizing method of the application is shown, is included the following steps:
S1, the corresponding depth of focus value Z of every wave beam is calibrated by given step, and by the depth of focus value Z2 after calibration
Replace the original depth of focus value Z of corresponding wave beam；
S2, according to after above-mentioned calibration depth of focus value Z2 beam collection carry out signal alignment, and obtain signal alignment after
BF1 signal collection in beam collection, then each BF1 signal that above-mentioned BF1 signal is concentrated is subjected to specified superposition and obtains BF2 signal.
As described in above-mentioned steps S1, the corresponding depth of focus value Z of every wave beam is calibrated by given step, and will be after calibration
Depth of focus value Z2 replace the corresponding original depth of focus value Z of wave beam, it should be noted that the depth of focus value into
The mode of row calibration includes two kinds in embodiments of the present invention, respectively according to linear array correction method and convex battle array correction method, wherein line
Battle array correction method passes through array number N_elements, adjacent array element spacing pitch, transmitting line number nlines, moving step length step, every
Sample point is to the distance z_sample1 and focal length F of launch point to described on the number of beams beam of secondary transmitting or receiving, emission lines
Depth of focus value Z is corrected；Convex battle array correction method passes through according to array number N_elements, adjacent array element spacing pitch, convex battle array
Center of circle radius R, transmitting line number nlines, moving step length step, the number of beams beam for emitting every time or receiving, focal length F, transmitting
The distance z_sample2 in sample point to the center of circle is corrected the depth of focus value Z on line.
As described in above-mentioned steps S2, signal alignment is carried out according to the depth of focus value Z2 beam collection after above-mentioned calibration, and obtain
The BF1 signal collection in beam collection after the number of winning the confidence alignment, then each BF1 signal that above-mentioned BF1 signal is concentrated is subjected to specified superposition
Obtain BF2 signal, it should be noted that the implementation method of step S2 includes two kinds, respectively according to linear array calibration method and convex battle array
Calibration method, wherein the depth of focus value Z2 after coordinate ref1, array element coordinate xe1, calibration of the linear array calibration method by receiving line will
Corresponding wave beam line carries out phase alignment, and the BF1 signal collection after signal is aligned is overlapped to obtain the BF2 signal；It is convex
Battle array calibration hair will by the depth of focus value Z2 after receiving the starting position coordinates M of line, the position coordinates N of array element, the calibration
Corresponding wave beam line carries out phase alignment, and the BF1 signal collection after signal is aligned is overlapped to obtain the BF2 signal.?
Before corresponding wave beam line is carried out phase alignment by the depth of focus value Z2 after calibration, spacial alignment finishes wave beam line,
It is delayed by the depth of focus value Z2 after above-mentioned calibration to wave beam.BF1 is obtained, judges the validity of data, then is superimposed
To BF2 since the phase of BF1 signal has been aligned, so only needing the same space data superimposed.
Referring to Fig. 2-4, a kind of step flow chart of beam synthesizing method embodiment 2 of the application is shown, specifically can wrap
Include following steps:
In embodiments of the present invention, the above-mentioned depth of focus value Z2 beam collection according to after above-mentioned calibration carries out signal alignment,
And the BF1 signal collection in the beam collection after signal alignment is obtained, then each BF1 signal that above-mentioned BF1 signal is concentrated is specified
Further include following steps before the step of superposition obtains BF2 signal:
S3, useful signal filtering is carried out to above-mentioned BF1 signal collection, and separator goes out the effective of above-mentioned BF1 signal concentration
BF1 signal and invalid BF1 signal.
As described in above-mentioned steps S3, useful signal filtering is carried out to above-mentioned BF1 signal collection, and separator goes out above-mentioned BF1
The effective BF1 signal and invalid BF1 signal that signal is concentrated are it should be noted that the implementation method of step S3 includes two kinds, respectively
For according to linear array calibration method and convex battle array calibration method, wherein linear array calibration method passes through the aperture A by ultrasonic transmitter, channel
Each sample point deviates the angle of emission lines on number N_channel, adjacent array element spacing pitch, focal length F and each wave beam line
angle1ibeamIt calculates and filters out the BF1 signal within the scope of effective sound field；Convex battle array calibration method is according to convex battle array center of circle radius R, channel
Sample point deviates the angle of emission lines on number N_channel, focal length F, adjacent array element spacing pitch and each wave beam line
angle2ibeamIt calculates and filters out the BF1 signal within the scope of effective sound field.
In embodiments of the present invention, the corresponding depth of focus value Z of every wave beam, and high-ranking officers are calibrated above by given step
Depth of focus value Z2 after standard replaces the step of corresponding wave beam original depth of focus value Z, includes the following steps:
S111, pass through array number N_elements, adjacent array element spacing pitch, transmitting line number nlines and moving step length
Step converses the spacing dx of each wave beam line in linear array；
S112, the number of beams beam for emitting or receiving by the spacing dx of wave beam line each in linear array, every time, in emission lines
The distance z_sample1 and focal length F of sample point to launch point converse the sample point in the emission lines of center to each wave beam respectively
The lateral distance Dx of virtual focus point on lineibeamWith fore-and-aft distance Dyibeam；
S113, according to the lateral distance Dx of the virtual focus point on the sample point in the emission lines of center to each wave beam lineibeam
With fore-and-aft distance DyibeamConverse the sample point in the emission lines of center to each wave beam line virtual focus point distance dr1ibeam；
S114, according to the sample point in the emission lines of center to the distance dr1 of each wave beam line virtual focus pointibeamIt is changed with focal length F
Depth of focus value Z2 after calculating the corresponding calibration of each wave beam line.
As described in above-mentioned steps S111, pass through array number N_elements, adjacent array element spacing pitch, transmitting line number
Nlines and moving step length step converses the spacing dx of each wave beam line in linear array, it should be noted that above-mentioned calculating process institute
The formula of use is preferred in embodiments of the present invention are as follows:
Wherein, it when moving step length step is 1, receives line spacing and is equal to emission lines spacing.Step and final line density
Correlation, when the bigger line density of step is higher, but the quantity of multi-beam superposition is just reduced, and physical relationship is as follows:
Beam superposition number × step=beam
In formula, beam is the number of beams for emitting or receiving every time.
As described in above-mentioned steps S112, pass through the spacing dx of wave beam line each in linear array, every time transmitting or the numbers of beams received
The distance z_sample1 and focal length F of sample point to launch point converse the sample in the emission lines of center respectively on amount beam, emission lines
The lateral distance Dx of virtual focus point on this to each wave beam lineibeamWith fore-and-aft distance Dyibeam, it should be noted that it is above-mentioned
Formula used by calculating process is preferred in embodiments of the present invention are as follows:
Dx=abs (- (beam-1) × dx × 1/2+ (ibeam-1) * dx) ibeam=1,2.....beam formula (2)
Dy=abs (z_sample-F)
In formula, abs is represented as the calculated result in formula and takes absolute value.
As described in above-mentioned steps S113, according to the virtual focus point on the sample point in the emission lines of center to each wave beam line
Lateral distance DxibeamWith fore-and-aft distance DyibeamConverse sample point in the emission lines of center to each wave beam line virtual focus point away from
From dr1ibeam, it should be noted that formula used by above-mentioned calculating process is preferred in embodiments of the present invention are as follows:
By taking Fig. 5 as an example: P1 in figure, Pn are that two sample points near field and far field connect by taking sample point P1 as an example in emission lines
It receives wave beam line to be illustrated by taking bm2 as an example, S is the virtual focus point on bm2 with transmitting focus with position.It is each on S to emission lines
The distance of sample point is drbm2, Dx1For sample point P any in emission lines to the lateral distance for receiving virtual focus point S；Dy1For longitudinal direction
Depth.By formula (3), the distance of sample point P1 to wave beam line bm2 virtual focus point S in emission lines can be calculated, i.e.,
dr1。
As described in above-mentioned steps S114, according to the sample point in the emission lines of center to the distance of each wave beam line virtual focus point
dr1ibeamDepth of focus value Z2 after conversing the corresponding calibration of each wave beam line with focal length F, it should be noted that above-mentioned to calculate
Formula used by journey is preferred in embodiments of the present invention are as follows:
Z2ibeam(j)=F ± dribeam(j) formula (4)
It should be noted that in formula, when z_sample1 is less than F use '-'；When z_sample1 is greater than F use '+'；When
When z_sample1 is equal to F, i.e. focal position.At this time focal point sound field concentrate, by other transmitting sound fields influenced it is smaller, because
This, when z_sample1 is equal to F, only with the data at current line focus.
By above-mentioned formula (4) can be obtained any bar will on received wave beam line ibeam any one point j correction it is poly-
Depth of focus angle value Z2, it is ready to receive to arrive the consistent multi-beam BF1 of phase using corresponding Z2 value when delay calculates.
In embodiments of the present invention, above-mentioned that useful signal filtering is carried out to above-mentioned BF1 signal collection, and separator is above-mentioned out
The step of effective BF1 signal and invalid BF1 signal that BF1 signal is concentrated, include the following steps:
S311, pass through the aperture A of ultrasonic transmitter, port number N_channel, adjacent array element spacing pitch and focal length F
Converse the angle, θ 1 of transmitting sound field；
S312, the angle angle1 that each sample point on each wave beam line deviates emission lines is calculatedibeam；
S313, judge above-mentioned angle angle1ibeamWhether θ 1/2 is greater than；
S314, if it is not, then determine with above-mentioned angle angle1ibeamCorresponding sample point is within the scope of effective sound field.
As described in above-mentioned steps S311, by between the aperture A, port number N_channel, adjacent array element of ultrasonic transmitter
The angle, θ 1 of transmitting sound field is conversed away from pitch and focal length F, it should be noted that the sound field shape of linear array is similar to hourglass-shaped, such as
Shown in Fig. 6.It is that A emits ultrasonic wave with pore size, carries out multi-beam reception, wherein P2, P3 is a specified reception wave beam
Two sample points on line, do vertical line to center emission lines from P2, P3, and the point that hangs down is N1, N2.
Therefore the angle of transmitting sound field is preferably calculated by following equation:
As described in above-mentioned steps S312, the angle angle1 that each sample point on each wave beam line deviates emission lines is calculatedibeam,
It should be noted that formula used by above-mentioned calculating process is preferred in embodiments of the present invention are as follows:
As described in above-mentioned steps S313, above-mentioned angle angle1 is judgedibeamWhether θ 1/2 is greater than, as shown in Figure 6, it is known that P2
Angle angle1 not within transmitting sound field, with center emission lines1Greater than θ 1/2；Sample point P3 is emitting within sound field,
The angle angle1 of itself and center emission lines2Less than θ 1/2.It therefore can be by judging each point and center emission lines on wave beam line
Whether the angle of focus judges each point within effective sound field.
As described in above-mentioned steps S314, if it is not, then determining and above-mentioned angle angle1ibeamCorresponding sample point is in effective
Within the scope of sound field, specifically, corresponding sample point is marked after judgement, it is specific to mark are as follows:
In embodiments of the present invention, the above-mentioned depth of focus value Z2 beam collection according to after above-mentioned calibration carries out signal alignment,
And the BF1 signal collection in the beam collection after signal alignment is obtained, then each BF1 signal that above-mentioned BF1 signal is concentrated is specified
The step of superposition obtains BF2 signal, includes the following steps:
S211, pass through the depth of focus value Z2 after the calibration of coordinate ref1, the array element coordinate xe1 and each wave beam line of reception line
Each array element is conversed to the distance dn1 for focusing sample pointibeam；
S212, pass through the depth of focus value Z2 after the calibration that receives the coordinate ref1 of line, array element coordinate xe1, each wave beam line
With each array element to the distance dn1 for focusing sample pointibeamConverse the delay τ 1 of each array element；
S213, corresponding wave beam line is carried out by phase alignment by the delay τ 1 of each array element, and after signal is aligned
The BF1 signal collection of wave beam line be overlapped to obtain above-mentioned BF2 signal.
As described in above-mentioned steps S211, after the calibration by coordinate ref1, the array element coordinate xe1 and each wave beam line that receive line
Depth of focus value Z2 converse each array element to focus sample point distance dn1ibeam, it should be noted that the meter for the distance that is delayed
For calculation method as shown in fig. 7, sample point P0 is moved on wave beam line, dn is each array element to focusing sample point distance, therefore, above-mentioned to change
Formula used by calculation process is preferred in embodiments of the present invention are as follows:
As described in above-mentioned steps S212, after the calibration by receiving the coordinate ref1 of line, array element coordinate xe1, each wave beam line
Depth of focus value Z2 and each array element to focus sample point distance dn1ibeamThe delay τ 1 for conversing each array element, needs to illustrate
It is that formula used by above-mentioned conversion process is preferred in embodiments of the present invention are as follows:
In formula, c is the bulk sound velocity in tissue, c=1540m/s.
As described in above-mentioned steps S213, corresponding wave beam line is carried out by phase alignment by the delay τ 1 of each array element,
And the BF1 signal collection of the wave beam line after being directed at signal is overlapped to obtain above-mentioned BF2 signal, it should be noted that believes in BF1
The specially following process of the step of number carrying out phase alignment:
S213-1, the delay τ 1 and above-mentioned angle angle1 within the scope of effective sound field according to each array elementibeamConversion
Sample point echo-signal is corresponded to out, it should be noted that formula used by above-mentioned conversion process is excellent in embodiments of the present invention
It is selected as:
In formula, S1ibeamIt (j) is the receives echo-signal of j-th of sample point on i-th beam received wave bunch, s1j(τ1)
For the echo that array element j is received in sub-aperture, WibeamIt (j) is signal amplitude trace-changing coefficient.
The multi-beam echo-signal of phase alignment can be obtained by above-mentioned formula (10), that is, BF1 signal.Wherein, BF1
The line number of signal is related with beam, emits 1 time and receives beam time, therefore transmitting nlines times, receive to obtain altogether nlines ×
Beam wave beam.
It should be noted that before the BF1 signal collection after being directed at signal is overlapped to obtain above-mentioned BF2 signal, also
Include the following steps:
S213-2, the echo-signal of above-mentioned corresponding sample point is collected to the BF1 signal collection after forming above-mentioned signal alignment, it will
The BF1 signal of nlines × beam wave beam collects for BF1 signal collection, during which is this Beam synthesis
Total collection of BF1 signal in all effective sound fields.
The BF1 signal concentrated after having collected BF1 signal collection to signal, which is overlapped, eventually forms BF2 signal, due to BF1
The phase of signal has been aligned, so only needing the same space data are superimposed to can be obtained BF2 signal.
In embodiments of the present invention, the corresponding depth of focus value Z of every wave beam, and high-ranking officers are calibrated above by given step
Depth of focus value Z2 after standard replaces the step of corresponding wave beam original depth of focus value Z, includes the following steps:
S121, convex battle array is conversed according to array number N_elements, adjacent array element spacing pitch and convex battle array center of circle radius R
Central angle beta；
S122, reception line is conversed according to transmitting line number nlines, moving step length step and above-mentioned convex battle array central angle beta
Between angle dx_beta；
Angle dx_beta, which is conversed, between number of beams beam and above-mentioned reception line that S123, basis emit or receive every time connects
Included angle A ngle_beam of the receipts wave beam to launching beam；
S124, according to focal length F, above-mentioned convex battle array center of circle radius R and receive wave beam to launching beam included angle A ngle_beam
Sample point is conversed respectively to the lateral distance Hx for receiving wave beam virtual focus pointibeamWith fore-and-aft distance Hyibeam, and according to above-mentioned
Above-mentioned convex battle array center of circle radius R and fore-and-aft distance HyibeamConverse correction focal length FA；
S125, according to the distance z_sample2 of sample point in emission lines to the center of circle, above-mentioned convex battle array center of circle radius R, laterally away from
From HxibeamWith fore-and-aft distance HyibeamConverse the sample point in the emission lines of center to each wave beam line virtual focus point distance
dr2ibeam；
S126, according to the sample point in above-mentioned center emission lines to the distance dr2 of each wave beam line virtual focus pointibeam, focal length F
The depth of focus value Z2 after the corresponding calibration of each wave beam line is conversed with correction focal length FA.
As described in above-mentioned steps S121, according to array number N_elements, adjacent array element spacing pitch and the convex battle array center of circle half
Diameter R converses convex battle array central angle beta, it should be noted that formula used by above-mentioned conversion process is in embodiments of the present invention
It is preferred that are as follows:
As described in above-mentioned steps S122, according to transmitting line number nlines, moving step length step and above-mentioned convex battle array central angle
Beta converses angle dx_beta between reception line, it should be noted that formula used by above-mentioned conversion process is implemented in the present invention
In example preferably are as follows:
It should be noted that same linear array, as step=1, angle is equal to angle between reception line between emission lines, and step is bigger, connects
Angle is smaller between take-up, then line density is higher, and in mono- timing of beam, then the numbers of beams being superimposed is fewer.General line density is higher, superposition
Number more multi-focus effect is better, but step and stacking fold are in inverse relation, therefore generally can be according to reality in practical application
Border situation weighs step, beam, superposition number, to obtain preferable focusing effect.
As described in above-mentioned steps S123, according to angle between each number of beams beam and above-mentioned reception line for emitting or receiving
Dx_beta converses the included angle A ngle_beam for receiving wave beam to launching beam, it should be noted that above-mentioned conversion process is adopted
Formula is preferred in embodiments of the present invention are as follows:
As described in above-mentioned steps S124, according to focal length F, above-mentioned convex battle array center of circle radius R and wave beam is received to launching beam
Included angle A ngle_beam converses sample point to the lateral distance Hx for receiving wave beam virtual focus point respectivelyibeamAnd fore-and-aft distance
Hyibeam, and according to above-mentioned convex battle array center of circle radius R and fore-and-aft distance HyibeamCorrection focal length FA is conversed, needs to illustrate
It is to converse sample point to the lateral Hx for receiving wave beam virtual focus pointibeam, fore-and-aft distance HyibeamUsed formula is in this hair
In bright embodiment preferably are as follows:
Hxibeam=| (F+R) × cos (Angle_beamibeam)|
Hyibeam=| (F+R) × sin (Angle_beamibeam) | formula (14).
As shown in figure 8, P4, P5 are two sample points near field and far field on transmitted wave bunch in figure by taking bm4 as an example, connect
Receipts wave beam virtual focus point is S, and F_curve is focal length circular arc, is done from the virtual focus point S of F_curve to center transmitted wave bunch
Vertical line, the point that hangs down are F_Al, and the transverse and longitudinal distance of sample point to S point is as shown in Figure 8；It is since Hy is invariable for linear array
Focal length, convex battle array virtual focus point vertical centre emission lines, the point that hangs down is the focus of correction, then correcting focal length is FA, therefore
FAibeam=Hyibeam- R formula (15).
As described in above-mentioned steps S125, according to the distance z_sample2 of sample point in emission lines to the center of circle, above-mentioned convex battle array circle
Heart radius R, lateral distance HxibeamWith fore-and-aft distance HyibeamThe sample point conversed in the emission lines of center is virtual to each wave beam line
The distance dr2 of focusibeam, it should be noted that formula used by above-mentioned conversion process is preferred in embodiments of the present invention are as follows:
As described in above-mentioned steps S126, according to the sample point in above-mentioned center emission lines to each wave beam line virtual focus point away from
From dr2ibeam, focal length F with correction focal length FA converse the depth of focus value Z2 after the corresponding calibration of each wave beam line, need to illustrate
It is that formula used by above-mentioned conversion process is preferred in embodiments of the present invention are as follows:
In embodiments of the present invention, above-mentioned that useful signal filtering is carried out to above-mentioned BF1 signal collection, and separator is above-mentioned out
The step of effective BF1 signal and invalid BF1 signal that BF1 signal is concentrated, include the following steps:
S321, it is conversed according to convex battle array center of circle radius R, port number N_channel, focal length F and adjacent array element spacing pitch
Emit the angle, θ 2 of sound field；
S322, the angle angle2 that sample point on each wave beam line deviates emission lines is obtainedibeam；
S323, judge above-mentioned angle angle2ibeamWhether θ 2/2 is greater than；
S324, if it is not, then determine with above-mentioned angle angle2ibeamCorresponding sample point is within the scope of effective sound field.
As described in above-mentioned steps S321, according to convex battle array center of circle radius R, port number N_channel, focal length F and adjacent array element
Spacing pitch converses the angle, θ 2 of transmitting sound field, it should be noted that passes through the aperture A of ultrasonic transmitter, port number N_
Channel, adjacent array element spacing pitch and focal length F converse the angle, θ 2 of transmitting sound field, it should be noted that the sound of convex battle array
Field such as linear array is hourglass-shaped, as shown in figure 9, emitting ultrasonic wave with aperture A, multi-beam is received.Wherein, it is sent out from focus F to center
Ray does vertical line, and the point that hangs down is N3.
Therefore the angle of transmitting sound field is preferably calculated by following equation:
As described in above-mentioned steps S322, the angle angle2 that sample point on each wave beam line deviates emission lines is obtainedibeam, need
Illustrate, the angle angle2 that sample point deviates emission lines on wave beam line is preferably calculated by the following method to be obtained, and example is passed through
Son explains step in detail below, specifically:
1. before sample point is located at and hangs down point N3
It calculates in triangle OFP1, ∠ OFP1 is angle21, wherein OF=a1 is enabled, FP1=b1, OP1=c1
A1=F+R
B1=dr (1:markibeam, ibeam) and ibeam=1,2,3.......beam
C1=z_sample+R
Mark is the point coordinate of each wave beam line colonel's positive focal length.
2. sample point is located at after vertical point N3
It calculates in triangle OFP2, ∠ OFP2 is angle22, wherein OF=a2 is enabled, FP2=b2, OP2=c2
A2=F+R
B2=dr (markibeam: end, ibeam) ibeam=1,2,3.......beam
C2=z_sample+R
As described in above-mentioned steps S323, above-mentioned angle angle2 is judgedibeamWhether θ 2/2 is greater than, as shown in figure 9, to receive
Wave beam line bm4 for, judge that each sample point is whether within effective sound field on wave beam line.P6, P7 are near field on bm4 and far
Two sample points of field, P6 is not emitting within sound field as seen from the figure, the angle angle2 with center emission lines1Greater than θ 2/
2；Angle angle2 of the sample point P7 within transmitting sound field, with center emission lines2Less than θ 2/2.Therefore it can pass through judgement
Whether the angle of each point and center transmitting line focus judges each point within effective sound field on wave beam line.
As described in above-mentioned steps S324, if it is not, then determining and above-mentioned angle angle2ibeamCorresponding sample point is in effective
Within the scope of sound field, specifically, corresponding sample point is marked after judgement, it is specific to mark are as follows:
In embodiments of the present invention, the above-mentioned depth of focus value Z2 beam collection according to after above-mentioned calibration carries out signal alignment,
And the BF1 signal collection in the beam collection after signal alignment is obtained, then each BF1 signal that above-mentioned BF1 signal is concentrated is specified
The step of superposition obtains BF2 signal, includes the following steps:
S221, according to receiving the starting position coordinates M of line, the position coordinates N of array element, the depth of focus value after above-mentioned calibration
Z2 converses array element to the distance dn2 of focus pointibeam；
S222, according to after above-mentioned calibration depth of focus value Z2 and array element to focus point distance dn2ibeamConverse each battle array
The delay τ 2 of member；
S223, corresponding wave beam line is carried out by phase alignment by the delay τ 2 of each array element, and after signal is aligned
The BF1 signal collection of wave beam line be overlapped to obtain above-mentioned BF2 signal.
As described in above-mentioned steps S221, according to reception the starting position coordinates M of line, the position coordinates N of array element, above-mentioned calibration
Depth of focus value Z2 afterwards converses array element to the distance dn2 of focus pointibeam, it should be noted that array element to focus point away from
From dn2ibeamCalculation formula pass through following procedure obtain:
As shown in Figure 10, M is the initial position for receiving line, and N point is the position of array element, focus point
Q is moved on wave beam line, and dn is array element to focus point distance, as seen from the figure:
In formula, (ref-xe) is distance of the array element to reception line initial position.
Since the distance of array element to focus point is dn
In formula, | MP |=Z2ibeam(j), ∠ NMP=δ
Therefore it can obtain:
As described in above-mentioned steps S222, according to after above-mentioned calibration depth of focus value Z2 and array element to focus point distance
dn2ibeamConverse the delay τ 2 of each array element, it should be noted that formula used by above-mentioned conversion process is implemented in the present invention
In example preferably are as follows:
τ2j=(Z2ibeam(j)+dnibeam(j))/c formula (24)
As described in above-mentioned steps S223, corresponding wave beam line is carried out by phase alignment by the delay τ 2 of each array element,
And the BF1 signal collection of the wave beam line after being directed at signal is overlapped to obtain above-mentioned BF2 signal, it should be noted that believes in BF1
The specially following process of the step of number carrying out phase alignment:
S223-1, the delay τ 2 and above-mentioned angle angle2 within the scope of effective sound field according to each array elementibeamConversion
Sample point echo-signal is corresponded to out, it should be noted that formula used by above-mentioned conversion process is excellent in embodiments of the present invention
It is selected as:
In formula, S2ibeamIt (j) is the receives echo-signal of j-th of sample point on i-th beam received wave bunch, s2j(τ2)
For the echo that array element j is received in sub-aperture.
The multi-beam echo-signal of phase alignment can be obtained by above-mentioned formula (25), that is, BF1 signal.Wherein, BF1
The line number of signal is related with beam, emits 1 time and receives beam time, therefore transmitting nlines times, receive to obtain altogether nlines ×
Beam wave beam.
It should be noted that before the BF1 signal collection after being directed at signal is overlapped to obtain above-mentioned BF2 signal, also
Include the following steps:
S223-2, the echo-signal of above-mentioned corresponding sample point is collected to the BF1 signal collection after forming above-mentioned signal alignment, it will
The BF1 signal of nlines × beam wave beam collects for BF1 signal collection, during which is this Beam synthesis
Total collection of BF1 signal in all effective sound fields.
The BF1 signal concentrated after having collected BF1 signal collection to signal, which is overlapped, eventually forms BF2 signal, due to BF1
The phase of signal has been aligned, so only needing the same space data are superimposed to can be obtained BF2 signal.
In practical applications, as shown in fig. 7, to receive 16 wave beams in figure, the multi-beam superposition in the case where step=4 is shown
It is intended to, other situations, such as wave beam 32, step=4 etc. can be obtained similarly.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of ultrasonic imaging method, include the following steps,
It obtains ultrasound data and is converted into corresponding digital signal；
Above-mentioned digital signal is synthesized into radiofrequency signal by specified correction；
Above-mentioned radiofrequency signal is isolated into carrier signal by specified signal processing；
Above-mentioned carrier signal is obtained into ultrasound image by specified image procossing,
Above-mentioned the step of above-mentioned digital signal is synthesized into radiofrequency signal by specified correction, including any of the above-described implementation
Beam synthesizing method described in example.
For device embodiment, since it is basically similar to the method embodiment, related so being described relatively simple
Place illustrates referring to the part of embodiment of the method.
Referring to Fig.1 1, show a kind of beam synthesizer of the invention, including following specific module:
Calibration module 1, for calibrating the corresponding depth of focus value Z of every wave beam by given step, and will be after calibration
Depth of focus value Z2 replaces the corresponding original depth of focus value Z of wave beam；
Synthesis module 2 for carrying out signal alignment according to the depth of focus value Z2 beam collection after above-mentioned calibration, and obtains letter
The BF1 signal collection in beam collection after number alignment, then each BF1 signal that above-mentioned BF1 signal is concentrated is subjected to specified superposition and is obtained
BF2 signal.
Above-mentioned calibration module 1 is generally used for calibrating the corresponding depth of focus value Z of every wave beam by given step, and will
Depth of focus value Z2 after calibration replaces the corresponding original depth of focus value Z of wave beam, it should be noted that deep to the focusing
The mode that angle value is calibrated includes two kinds in embodiments of the present invention, respectively according to linear array correction method and convex battle array correction method,
Wherein, linear array correction method passes through array number N_elements, adjacent array element spacing pitch, transmitting line number nlines, moving step length
Distance z_sample1 and focal length of the sample point to launch point on step, every time the number of beams beam of transmitting or receiving, emission lines
F is corrected the depth of focus value Z；Convex battle array correction method passes through according to array number N_elements, adjacent array element spacing
Pitch, convex battle array center of circle radius R, transmitting line number nlines, moving step length step, each transmitting or the number of beams beam received,
The distance z_sample2 in sample point to the center of circle is corrected the depth of focus value Z on focal length F, emission lines.
Above-mentioned synthesis module 2 is generally used for carrying out signal alignment according to the depth of focus value Z2 beam collection after above-mentioned calibration,
And the BF1 signal collection in the beam collection after signal alignment is obtained, then each BF1 signal that above-mentioned BF1 signal is concentrated is specified
Superposition obtains BF2 signal, it should be noted that the implementation method of step S2 includes two kinds, respectively according to linear array calibration method and
Convex battle array calibration method, wherein the depth of focus value after coordinate ref1, array element coordinate xe1, calibration of the linear array calibration method by receiving line
Corresponding wave beam line is carried out phase alignment by Z2, and the BF1 signal collection after signal is aligned is overlapped to obtain the BF2 letter
Number；Convex battle array calibration hair passes through the depth of focus value after receiving the starting position coordinates M of line, the position coordinates N of array element, the calibration
Corresponding wave beam line is carried out phase alignment by Z2, and the BF1 signal collection after signal is aligned is overlapped to obtain the BF2 letter
Number.Before corresponding wave beam line is carried out phase alignment by the depth of focus value Z2 after calibration, wave beam line spacial alignment
It finishes, is delayed by the depth of focus value Z2 after above-mentioned calibration to wave beam.BF1 is obtained, judges the validity of data, then is folded
Add to obtain BF2 since the phase of BF1 signal has been aligned, so only needing the same space data superimposed.
Show a kind of supersonic imaging device of the invention, including following specific module:
Module is obtained, for obtaining ultrasound data and being converted into corresponding digital signal；
Beam synthesis module, for above-mentioned digital signal to be synthesized radiofrequency signal by specified correction；
Separation module, for above-mentioned radiofrequency signal to be isolated carrier signal by specified signal processing；
Image-forming module, for above-mentioned carrier signal to be obtained ultrasound image by specified image procossing；
Above-mentioned Beam synthesis module, including beam synthesizer described in any of the above-described embodiment.
A kind of supersonic imaging apparatus of the invention is shown, including generates ultrasonic signal and is radiated tested tissue neutralization
Received sound wave, is converted to the vibrating sensor of electric signal, will receive signal by the signal transceiver for absorbing reflected sonic signals
Sampling and digitized analog-digital converter (A/D), the time gain compensation of the compensation decaying of the ultrasonic amplitude as caused by depth
Device (TGC) converts digital signal to the reception beam synthesizer of RF signal, and RF signal is carried out envelope extraction and demodulation process
The signal processor for isolating carrier signal converts the digital scan that carrier signal is scanned conversion and back-end image processing
The display of device DSC and the image for finally showing.
Wherein, above-mentioned reception beam synthesizer, including phase correction module, the first time wave beam for synthesizing BF1 signal
Beam synthesis module, effective sound field judgment module and second of beam synthesizer for synthesizing BF2 signal, and above-mentioned phase school
Positive module, the first time multi-beam beam synthesis module for synthesizing BF1 signal, effective sound field judgment module and for synthesizing BF2
Second of beam synthesizer of signal is sequentially connected.
Referring to Fig.1 2, a kind of computer equipment for realizing the beam synthesizing method of the invention is shown, it specifically can be with
Include the following:
Above-mentioned computer equipment 12 is showed in the form of universal computing device, the component of computer equipment 12 may include but
Be not limited to: one or more processor or processing unit 16, system storage 28, connecting different system components (including is
Unite memory 28 and processing unit 16) bus 18.
Bus 18 indicates one of a few 18 structures of class bus or a variety of, including memory bus 18 or memory control
Device, peripheral bus 18, graphics acceleration port, processor or the office using 18 structure of any bus in a variety of 18 structures of bus
Domain bus 18.For example, these architectures include but is not limited to industry standard architecture (ISA) bus 18, microchannel
Architecture (MAC) bus 18, enhanced isa bus 18, audio-video frequency electronic standard association (VESA) local bus 18 and outer
Enclose component interconnection (PCI) bus 18.
Computer equipment 12 typically comprises a variety of computer system readable media.These media can be it is any can be by
The usable medium that computer equipment 12 accesses, including volatile and non-volatile media, moveable and immovable medium.
System storage 28 may include the computer system readable media of form of volatile memory, such as arbitrary access
Memory (RAM) 30 and/or cache memory 32.Computer equipment 12 may further include other movement/it is not removable
Dynamic, volatile/non-volatile computer decorum storage medium.Only as an example, storage system 34 can be used for read and write can not
Mobile, non-volatile magnetic media (commonly referred to as " hard disk drive ").Although being not shown in Figure 12, can provide for can
The disc driver of mobile non-volatile magnetic disk (such as " floppy disk ") read-write, and to removable anonvolatile optical disk (such as CD~
ROM, DVD~ROM or other optical mediums) read-write CD drive.In these cases, each driver can pass through one
A or multiple data mediums interface is connected with bus 18.Memory may include at least one program product, the program product
With one group of (for example, at least one) program module 42, these program modules 42 are configured to perform the function of various embodiments of the present invention
Energy.
Program/utility 40 with one group of (at least one) program module 42, can store in memory, for example,
Such program module 42 includes --- but being not limited to --- operating system, one or more application program, other program moulds
It may include the realization of network environment in block 42 and program data, each of these examples or certain combination.Program mould
Block 42 usually executes function and/or method in embodiment described in the invention.
Computer equipment 12 can also with one or more external equipments 14 (such as keyboard, sensing equipment, display 24,
Camera etc.) communication, the equipment interacted with the computer equipment 12 can be also enabled a user to one or more to be communicated, and/
Or with enable the computer equipment 12 and one or more other calculate any equipment that equipment are communicated (such as network interface card,
Modem etc.) communication.This communication can be carried out by interface input/output (I/O) 22.Also, computer equipment
12 can also by network adapter 20 and one or more network (such as local area network (LAN)), wide area network (WAN) and/or
Public network (such as internet) communication.As shown, network adapter 20 passes through other of bus 18 and computer equipment 12
Module communication.It should be understood that although being not shown in Figure 12 other hardware and/or software can be used in conjunction with computer equipment 12
Module, including but not limited to: microcode, device driver, redundant processing unit 16, external disk drive array, RAID system,
Tape drive and data backup storage system 34 etc..
Processing unit 16 by the program that is stored in system storage 28 of operation, thereby executing various function application and
Data processing, such as realize beam synthesizing method provided by the embodiment of the present invention.
That is, above-mentioned processing unit 16 is realized when executing above procedure: by carrying out conspicuousness inspection to original galactophore image
It surveys, adjusts the contrast between the target area and nontarget area of above-mentioned original galactophore image, obtain Saliency maps picture；To upper
It states Saliency maps picture and carries out the processing of mammary gland contours segmentation, determine effective mammary region in above-mentioned Saliency maps picture；Determination is had
Above-mentioned Saliency maps picture after imitating mammary region carries out region coarse segmentation, retains the coarse segmentation region for meeting specified gray threshold,
Obtain coarse segmentation image；The specified coarse segmentation region in above-mentioned coarse segmentation image is finely divided using K-means clustering method
It cuts, obtains thin segmented image；False positive region is carried out to above-mentioned thin segmented image according to designated modality feature to filter out, and is obtained true
Positive region image；True positives zone marker in true positives area image is subjected to mammary gland profile to Saliency maps picture above-mentioned
In image after dividing processing, the galactophore image for having marked lesion region is obtained.
It, can be preferably by realizing the spacial alignment and phase alignment of multi-beam in any of the above-described inventive embodiments
Suppressed sidelobes improves image resolution ratio and signal-to-noise ratio；It is superimposed by multi-beam and greatly improves line density, improve the sky of image
Between resolution ratio；This method realizes that algorithm is simple, and multirate, high noise can be realized in the case where consuming a small amount of process resource
Than high-resolution effect；Device makes Beam synthesis by increasing calibration module and synthesis module in Beam synthesis module
More data are more accurate afterwards；And is judged by sound field, to effective sound field data investigation, realized preferably compared to conventional method
Focusing effect improves picture quality.
Above to a kind of beam synthesizing method provided herein and device, it is described in detail, it is used herein
The principle and implementation of this application are described for specific case, and the above embodiments are only used to help understand
The present processes and its core concept；At the same time, for those skilled in the art is having according to the thought of the application
There will be changes in body embodiment and application range, in conclusion the content of the present specification should not be construed as to the application
Limitation.
Claims (10)
1. a kind of beam synthesizing method, which comprises the steps of:
The corresponding depth of focus value Z of every wave beam is calibrated by given step, and the depth of focus value Z2 after calibration is replaced into phase
The corresponding original depth of focus value Z of wave beam；
Signal alignment is carried out according to the depth of focus value Z2 beam collection after the calibration, and is obtained in the beam collection after signal alignment
BF1 signal collection, then each BF1 signal that the BF1 signal is concentrated is subjected to specified superposition and obtains BF2 signal.
2. the method according to claim 1, wherein the depth of focus value Z2 wave beam according to after the calibration
Collection carry out signal alignment, and obtain signal alignment after beam collection in BF1 signal collection, then by the BF1 signal concentrate it is each
Further include following steps before BF1 signal carries out the step of specified superposition obtains BF2 signal:
Useful signal filtering is carried out to the BF1 signal collection, and separator goes out effective BF1 signal that the BF1 signal is concentrated
With invalid BF1 signal.
3. according to the method described in claim 2, it is characterized in that, described by given step to calibrate every wave beam corresponding poly-
Depth of focus angle value Z, and the step of depth of focus value Z2 after calibration is replaced into corresponding wave beam original depth of focus value Z, including
Following steps:
It is conversed by array number N_elements, adjacent array element spacing pitch, transmitting line number nlines and moving step length step
The spacing dx of each wave beam line in linear array；
It is arrived by sample point on the spacing dx of wave beam line each in linear array, every time transmitting or the number of beams beam received, emission lines
The distance z_sample1 and focal length F of launch point converse the void on the sample point in the emission lines of center to each wave beam line respectively
The lateral distance Dx of quasi- focusibeamWith fore-and-aft distance Dyibeam；
According to the lateral distance Dx of the virtual focus point on the sample point in the emission lines of center to each wave beam lineibeamAnd fore-and-aft distance
DyibeamConverse the sample point in the emission lines of center to each wave beam line virtual focus point distance dr1ibeam；
According to the sample point in the emission lines of center to the distance dr1 of each wave beam line virtual focus pointibeamEach wave beam is conversed with focal length F
Depth of focus value Z2 after the corresponding calibration of line.
4. according to the method described in claim 2, it is characterized in that, it is described to the BF1 signal collection carry out useful signal filtering,
And separator goes out the step of effective BF1 signal and invalid BF1 signal that the BF1 signal is concentrated, includes the following steps:
It is set out by the aperture A of ultrasonic transmitter, port number N_channel, adjacent array element spacing pitch and focal length F conversion
Penetrate the angle, θ 1 of sound field；
Calculate the angle angle1 that each sample point on each wave beam line deviates emission linesibeam；
Judge the angle angle1ibeamWhether θ 1/2 is greater than；
If it is not, then determining and the angle angle1ibeamCorresponding sample point is within the scope of effective sound field.
5. according to the method described in claim 4, it is characterized in that, the depth of focus value Z2 wave beam according to after the calibration
Collection carry out signal alignment, and obtain signal alignment after beam collection in BF1 signal collection, then by the BF1 signal concentrate it is each
BF1 signal carries out the step of specified superposition obtains BF2 signal, includes the following steps:
It is conversed respectively by the depth of focus value Z2 after the calibration of coordinate ref1, the array element coordinate xe1 and each wave beam line that receive line
Distance dn1 of the array element to focusing sample pointibeam；
It is arrived by depth of focus value Z2 after receiving the calibration of the coordinate ref1 of line, array element coordinate xe1, each wave beam line and each array element
Focus the distance dn1 of sample pointibeamConverse the delay τ 1 of each array element；
Corresponding BF1 signal is subjected to phase alignment by the delay τ 1 of each array element, and the BF1 signal after signal is aligned
Collection is overlapped to obtain the BF2 signal.
6. a kind of ultrasonic imaging method, includes the following steps,
It obtains ultrasound data and is converted into corresponding digital signal；
The digital signal is synthesized into radiofrequency signal by specified correction；
The radiofrequency signal is isolated into carrier signal by specified signal processing；
The carrier signal is obtained into ultrasound image by specified image procossing,
It is characterized in that, the described the step of digital signal is synthesized into radiofrequency signal by specified correction, including right
It is required that beam synthesizing method described in 1-5 any one.
7. a kind of beam synthesizer, which is characterized in that including following specific module:
Calibration module, for calibrating the corresponding depth of focus value Z of every wave beam by given step, and the focusing after calibration is deep
Angle value Z2 replaces the corresponding original depth of focus value Z of wave beam；
Synthesis module for carrying out signal alignment according to the depth of focus value Z2 beam collection after the calibration, and obtains signal pair
BF1 signal collection in beam collection after standard, then each BF1 signal that the BF1 signal is concentrated is subjected to specified superposition and obtains BF2 letter
Number.
8. a kind of supersonic imaging device, including following specific module:
Module is obtained, for obtaining ultrasound data and being converted into corresponding digital signal；
Beam synthesis module, for the digital signal to be synthesized radiofrequency signal by specified correction；
Separation module, for the radiofrequency signal to be isolated carrier signal by specified signal processing；
Image-forming module, for the carrier signal to be obtained ultrasound image by specified image procossing；
It is characterized in that, the Beam synthesis module, including beam synthesizer as claimed in claim 7.
9. a kind of supersonic imaging apparatus, which is characterized in that absorbed including generating ultrasonic signal and being radiated tested tissue and neutralize
Received sound wave is converted to the vibrating sensor of electric signal by the signal transceiver of reflected sonic signals, will receive signal sampling
With digitized analog-digital converter (A/D), the time gain compensation device of the decaying of the ultrasonic amplitude as caused by depth is compensated
(TGC), digital signal is converted to the reception beam synthesizer of RF signal, RF signal is subjected to envelope extraction and demodulation process point
Carrier signal is scanned the digital scan converter of conversion and back-end image processing by the signal processor for separating out carrier signal
The display of DSC and the image for finally showing.
Wherein, the reception beam synthesizer, including phase correction module, the first time multi-beam beam for synthesizing BF1 signal
Synthesis module, effective sound field judgment module and second of beam synthesizer for synthesizing BF2 signal, and the phasing mould
Block, the first time multi-beam beam synthesis module for synthesizing BF1 signal, effective sound field judgment module and for synthesizing BF2 signal
Second of beam synthesizer be sequentially connected.
10. a kind of computer equipment, can run on a memory and on a processor including memory, processor and storage
Computer program, which is characterized in that the processor is realized when executing described program such as any one of claim 1~6 institute
The method stated.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910204036.4A CN110101409B (en) | 2019-03-18 | 2019-03-18 | Beam forming method, ultrasonic imaging method, device and equipment |
| CN202111470011.2A CN114129185B (en) | 2019-03-18 | 2019-03-18 | Beam forming method, ultrasonic imaging method, device and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910204036.4A CN110101409B (en) | 2019-03-18 | 2019-03-18 | Beam forming method, ultrasonic imaging method, device and equipment |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111470011.2A Division CN114129185B (en) | 2019-03-18 | 2019-03-18 | Beam forming method, ultrasonic imaging method, device and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110101409A true CN110101409A (en) | 2019-08-09 |
| CN110101409B CN110101409B (en) | 2022-03-22 |
Family
ID=67484350
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910204036.4A Active CN110101409B (en) | 2019-03-18 | 2019-03-18 | Beam forming method, ultrasonic imaging method, device and equipment |
| CN202111470011.2A Active CN114129185B (en) | 2019-03-18 | 2019-03-18 | Beam forming method, ultrasonic imaging method, device and equipment |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111470011.2A Active CN114129185B (en) | 2019-03-18 | 2019-03-18 | Beam forming method, ultrasonic imaging method, device and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN110101409B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110974304A (en) * | 2019-12-13 | 2020-04-10 | 山东大学齐鲁医院 | Ultrasonic beam synthesis system and method based on wearable flexible ultrasonic transducer |
| CN111772679A (en) * | 2020-08-04 | 2020-10-16 | 深圳蓝韵医学影像有限公司 | Ultrasonic imaging method and device, computer equipment and computer readable medium |
| CN112842382A (en) * | 2019-11-12 | 2021-05-28 | 通用电气精准医疗有限责任公司 | Method and system for streaming channel data to apply non-linear beamforming |
| WO2021114071A1 (en) * | 2019-12-10 | 2021-06-17 | 深圳大学 | Transcranial ultrasonic imaging method and apparatus, and computer-readable storage medium |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH01207042A (en) * | 1987-12-11 | 1989-08-21 | General Electric Co <Ge> | Method and apparatus for reducing repeated adaptation of phase differential effect |
| CN1788685A (en) * | 2004-12-15 | 2006-06-21 | 深圳迈瑞生物医疗电子股份有限公司 | Receiving method and its device based on double beam and synthetic aperture |
| CN101097256A (en) * | 2006-06-28 | 2008-01-02 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic beam synthetic focusing parameter real-time algorithm and device thereof |
| CN101664321A (en) * | 2009-09-07 | 2010-03-10 | 无锡祥生科技有限公司 | Group sound velocity real-time adjustable ultrasound diagnostic equipment and wave beam synthesis method thereof |
| CN101865998A (en) * | 2010-05-19 | 2010-10-20 | 深圳市蓝韵实业有限公司 | Delay parameter real-time computation and production line type multi-beam synthesis method |
| CN101915916A (en) * | 2010-08-04 | 2010-12-15 | 深圳市蓝韵实业有限公司 | Method for processing ultrasonic emission delay |
| CN102670250A (en) * | 2012-05-24 | 2012-09-19 | 深圳市开立科技有限公司 | Array calibration system and method of ultrasonic transducer |
| CN102727256A (en) * | 2012-07-23 | 2012-10-17 | 重庆博恩富克医疗设备有限公司 | Dual focusing beam forming method and device based on virtual array elements |
| CN102846336A (en) * | 2011-06-27 | 2013-01-02 | 深圳市蓝韵实业有限公司 | Ultrasonic imaging delay control method and device |
| CN103654732A (en) * | 2013-12-31 | 2014-03-26 | 南京大学 | Photoacoustic image optimization method based on linear delay compensation |
| US20150359521A1 (en) * | 2014-06-12 | 2015-12-17 | General Electric Company | Method and systems for adjusting a pulse generated for ultrasound multi-line transmit |
| CN105832366A (en) * | 2016-03-18 | 2016-08-10 | 深圳安盛生物医疗技术有限公司 | Real-time delay computing method applied to beam-forming process |
| CN106971055A (en) * | 2016-01-12 | 2017-07-21 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic contrast imaging method and system |
| CN107569254A (en) * | 2016-07-05 | 2018-01-12 | 柯尼卡美能达株式会社 | Ultrasonic signal processing unit, ultrasonic signal processing method and diagnostic ultrasound equipment |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6432056B1 (en) * | 1999-10-08 | 2002-08-13 | Atl Ultrasound | Ultrasonic diagnostic imaging system with high frame rate synthetic transmit focus |
| CN102895000B (en) * | 2012-11-06 | 2014-06-18 | 重庆大学 | Double-focusing beamforming method based on self-adaptive weighting |
| WO2016125509A1 (en) * | 2015-02-04 | 2016-08-11 | 日立アロカメディカル株式会社 | Ultrasound imaging device and ultrasound signal processing method |
| JP6380194B2 (en) * | 2015-03-27 | 2018-08-29 | コニカミノルタ株式会社 | Ultrasonic signal processing apparatus and ultrasonic diagnostic apparatus |
| CN109363714B (en) * | 2018-09-17 | 2020-12-15 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic imaging equipment and ultrasonic imaging method thereof |
-
2019
- 2019-03-18 CN CN201910204036.4A patent/CN110101409B/en active Active
- 2019-03-18 CN CN202111470011.2A patent/CN114129185B/en active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH01207042A (en) * | 1987-12-11 | 1989-08-21 | General Electric Co <Ge> | Method and apparatus for reducing repeated adaptation of phase differential effect |
| CN1788685A (en) * | 2004-12-15 | 2006-06-21 | 深圳迈瑞生物医疗电子股份有限公司 | Receiving method and its device based on double beam and synthetic aperture |
| CN101097256A (en) * | 2006-06-28 | 2008-01-02 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic beam synthetic focusing parameter real-time algorithm and device thereof |
| CN101664321A (en) * | 2009-09-07 | 2010-03-10 | 无锡祥生科技有限公司 | Group sound velocity real-time adjustable ultrasound diagnostic equipment and wave beam synthesis method thereof |
| CN101865998A (en) * | 2010-05-19 | 2010-10-20 | 深圳市蓝韵实业有限公司 | Delay parameter real-time computation and production line type multi-beam synthesis method |
| CN101915916A (en) * | 2010-08-04 | 2010-12-15 | 深圳市蓝韵实业有限公司 | Method for processing ultrasonic emission delay |
| CN102846336A (en) * | 2011-06-27 | 2013-01-02 | 深圳市蓝韵实业有限公司 | Ultrasonic imaging delay control method and device |
| CN102670250A (en) * | 2012-05-24 | 2012-09-19 | 深圳市开立科技有限公司 | Array calibration system and method of ultrasonic transducer |
| CN102727256A (en) * | 2012-07-23 | 2012-10-17 | 重庆博恩富克医疗设备有限公司 | Dual focusing beam forming method and device based on virtual array elements |
| CN103654732A (en) * | 2013-12-31 | 2014-03-26 | 南京大学 | Photoacoustic image optimization method based on linear delay compensation |
| US20150359521A1 (en) * | 2014-06-12 | 2015-12-17 | General Electric Company | Method and systems for adjusting a pulse generated for ultrasound multi-line transmit |
| CN106971055A (en) * | 2016-01-12 | 2017-07-21 | 深圳迈瑞生物医疗电子股份有限公司 | Ultrasonic contrast imaging method and system |
| CN105832366A (en) * | 2016-03-18 | 2016-08-10 | 深圳安盛生物医疗技术有限公司 | Real-time delay computing method applied to beam-forming process |
| CN107569254A (en) * | 2016-07-05 | 2018-01-12 | 柯尼卡美能达株式会社 | Ultrasonic signal processing unit, ultrasonic signal processing method and diagnostic ultrasound equipment |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112842382A (en) * | 2019-11-12 | 2021-05-28 | 通用电气精准医疗有限责任公司 | Method and system for streaming channel data to apply non-linear beamforming |
| WO2021114071A1 (en) * | 2019-12-10 | 2021-06-17 | 深圳大学 | Transcranial ultrasonic imaging method and apparatus, and computer-readable storage medium |
| CN110974304A (en) * | 2019-12-13 | 2020-04-10 | 山东大学齐鲁医院 | Ultrasonic beam synthesis system and method based on wearable flexible ultrasonic transducer |
| CN111772679A (en) * | 2020-08-04 | 2020-10-16 | 深圳蓝韵医学影像有限公司 | Ultrasonic imaging method and device, computer equipment and computer readable medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114129185B (en) | 2023-11-07 |
| CN110101409B (en) | 2022-03-22 |
| CN114129185A (en) | 2022-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110101409A (en) | Beam synthesizing method, ultrasonic imaging method, device and equipment | |
| Perrot et al. | So you think you can DAS? A viewpoint on delay-and-sum beamforming | |
| US8038620B2 (en) | Fresnel zone imaging system and method | |
| US7740583B2 (en) | Time delay estimation method and system for use in ultrasound imaging | |
| US8045777B2 (en) | Clutter suppression in ultrasonic imaging systems | |
| US8968199B2 (en) | Spatial compound imaging in an ultrasound system | |
| JPH10506800A (en) | Adjustable frequency scanning method and apparatus in ultrasound images | |
| US20130172745A1 (en) | Providing vector doppler image based on decision data in ultrasound system | |
| US11406362B2 (en) | Providing user interface in ultrasound system | |
| JP2009505771A (en) | Ultrasound imaging system and method for flow imaging with real-time space synthesis | |
| US20070167802A1 (en) | Accurate time delay estimation method and system for use in ultrasound imaging | |
| JP2000157548A (en) | Method and system for imaging ultrasonic wave scattered body | |
| US20130165784A1 (en) | Providing motion profile information of target object in ultrasound system | |
| CN104272134A (en) | Clutter suppression in ultrasonic imaging systems | |
| US9474510B2 (en) | Ultrasound and system for forming an ultrasound image | |
| US9232932B2 (en) | Providing motion mode image in ultrasound system | |
| US20130172749A1 (en) | Providing doppler spectrum images corresponding to at least two sample volumes in ultrasound system | |
| US20130172755A1 (en) | Providing turbulent flow information based on vector doppler in ultrasound system | |
| JP5069022B2 (en) | Method and system for accurate time delay estimation for use in ultrasound imaging | |
| US20130172747A1 (en) | Estimating motion of particle based on vector doppler in ultrasound system | |
| EP2610638B1 (en) | Forming vector information based on vector doppler in ultrasound system | |
| US20040054285A1 (en) | Elevation compounding for ultrasound imaging | |
| US20130172744A1 (en) | Providing particle flow image in ultrasound system | |
| KR20070089101A (en) | Accurate time delay estimation method and system for use in ultrasound imaging | |
| US20130165793A1 (en) | Providing doppler information of target object based on vector doppler in ultrasound system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address |
Address after: No.103, baguang District Service Center, No.2 BaiShaWan Road, baguang community, Kuiyong street, Dapeng New District, Shenzhen, Guangdong 518000 Patentee after: Shenzhen Lanying Medical Technology Co.,Ltd. Address before: 518000 Shenzhen, Baoan District, Guangdong Shiyan street, Lang Xin neighborhood committee, brick factory village, Liyuan Industrial Zone, hung lung hi tech Industrial Park (blue Yun Industrial Park), 1 1 buildings. Patentee before: SHENZHEN LANYUN MEDICAL IMAGE CO.,LTD. |
|
| CP03 | Change of name, title or address |