CN101371789A - Real-time computing method and device for focalization delay parameters of beam-forming - Google Patents
Real-time computing method and device for focalization delay parameters of beam-forming Download PDFInfo
- Publication number
- CN101371789A CN101371789A CNA2007100766203A CN200710076620A CN101371789A CN 101371789 A CN101371789 A CN 101371789A CN A2007100766203 A CNA2007100766203 A CN A2007100766203A CN 200710076620 A CN200710076620 A CN 200710076620A CN 101371789 A CN101371789 A CN 101371789A
- Authority
- CN
- China
- Prior art keywords
- delay
- parameter
- interpolation
- real
- focusing
- 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.)
- Pending
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention relates to a focus delay parameter real-time computing method for beam forming and a device, used for the beam forming of an ultrasonic imaging device; according to the requirements of a near-field region for compensating delay coefficient and aiming at each receiving point of each receiving channel, real-time computing is carried out to provide focus delay parameters so as not to cause periodic delay error. The real-time computing is carried out by adopting linear interpolation processing and can greatly simplify the hardware design of the device when the point number N of interpolation is the exponent of 2. Furthermore, the method and the device have the advantages of low cost and easy implementation.
Description
Technical field the present invention relates to ultrasonic technique, relates in particular to the signal processing in the medical ultrasound imaging system, particularly relates to ultrasonic signal and receives synthetic parameter determination method and device.
The theory diagram of background technology typical case ultrasonic image-forming system as shown in Figure 1.Front end comprises a probe of being made up of a plurality of array elements, and each array element all has the function of the signal of telecommunication and the mutual conversion of acoustical signal.Described probe is connected with emission/reception conversion portion.When system is in emission state, pulse generator produces corresponding pumping signal according to predetermined transmitted waveform and emission delay parameter, be sent to probe by described emission/reception conversion portion, each array element in this probe transmitting aperture is that acoustical signal forms the emission sound field with electrical signal conversion.When system is in accepting state, it is the signal of telecommunication from the echo signal of detected body and conversion that the array element of probe in the receiving aperture receives, be transferred to the system handles part by described emission/reception conversion portion again, further be converted to digital signal by A/D (mould/number) transducer; Because each array element connects a passage, the electrical signal data of a plurality of like this passages enters beam synthesizer, be synthesized into the rf data that receives line by this beam synthesizer according to different time-delays, utilize detector to finish the signal demodulation again, ask works of treatment such as envelope and double sampling to obtain scan-line data, being processed into through digital scan converter (DSC) can be for the rectangular coordinate ultrasonic image data of display demonstration.
In above-mentioned reception processing procedure, focusing on retardation coefficient is the synthetic indispensable important parameter of wave beam.Related method thereof commonly used is delay to be divided into thick time-delay finish with thin time-delay, wherein the precision of thick time-delay is decided by A/D converter, the precision of thin time-delay then realizes by interpolation, after depositing the echo data of the higher delay precision that obtains after the interpolation in memorizer, read the address by control and read this memorizer, thereby realize that the echo data that different passages are received postpones the purpose that focuses on.Therefore, produce the described method of reading the address and directly determined the synthetic focusing accuracy of wave beam.
What existing real-time implementation focused on that the method that postpones adopts is, does once and revises reading the address every N point (N〉8), that is to say, whenever does once accurate the delay at a distance of N point and focuses on, and the N point is read the address all equally with interior consecutive points.Its weak point is: for far-field signal, the variation that focuses on retardation coefficient is less, and the delay error that brings thus almost can be ignored; But at near-field region, the variation of retardation coefficient is bigger, adopts this method will inevitably introduce tangible delay error.Especially when the port number of system, when the array number that just participates in receiving increased, situation can be more obvious.
The summary of the invention the technical problem to be solved in the present invention is at above-mentioned the deficiencies in the prior art part, and a kind of method and device are proposed, be used in the wave beam of ultrasonic image-forming system synthetic in, the synthetic focusing delay parameter of compute beam in real time is to reduce the delay error of near-field region.
For solving the problems of the technologies described above, of the present inventionly be contemplated that substantially: consider the demand of near-field region, take pointwise to calculate and obviously to reduce delay error the retardation coefficient precision; Therefore, the delay focusing parameter at each receive path will be calculated and provide in system in real time, realize the vernier focusing of pointwise.
As the technical scheme that realizes the present invention's design be, a kind of method of real-time compute beam synthetic focusing delay parameter is provided, the wave beam that is used for supersonic imaging apparatus is synthetic, comprises step:
A. will be sent to beam synthesizer from the ultrasonic echo electrical signal data of a plurality of receive paths of same received beam;
B. the data to these passages postpone to handle with the corresponding focusing parameter that postpones respectively;
C. the synthetic ultrasonic echo data that obtains a scanning line;
Especially, the described delay focusing parameter among the step B is that each acceptance point at each receive path calculates in real time and provides, thereby reduces the delay error of near-field region.
In the such scheme, described real-time calculating adopts interpolation processing to calculate, and comprises step:
A. to the real-time respectively delay focusing parameter that calculates on putting two times of reception of each receive path; If between these two acceptance points N-1 acceptance point arranged, N is a natural number;
B. postponing focusing parameter with these two simultaneously is that starting point and terminal point carry out interpolation calculation, obtains the real time delay focusing parameter of N-1 interpolation point, makes a corresponding described N-1 acceptance point respectively.
In the such scheme, described interpolation processing adopts linear interpolation method, thereby is that to postpone focusing parameters with described two be that starting point and terminal point carry out equally spaced linear interpolation and calculate among the step b.
In the such scheme, the described interpolation N-1 that counts is decided by that device hardware realizes the precision of thin time-delay, and the interpolation error of being brought by described interpolation calculation; The former is greater than the latter.
In the such scheme, described interval N is taken as 2 power.
As the technical scheme that realizes the present invention's design still, a kind of device of real-time compute beam synthetic focusing delay parameter is provided, the wave beam synthesis unit that is used for supersonic imaging apparatus, this device comprises a parameter calculation unit, be used for calculating the focusing delay parameter of being scheduled to acceptance point according to the probe element position, especially, this device also comprises a plurality of interpolation circuits, and each described interpolation circuit is corresponding to a receive path.
In the such scheme, described interpolation circuit comprises three depositors, and first depositor is used for receiving and keeping in the focusing delay parameter corresponding to this passage of described parameter calculation unit output under the control of first data enable signal; The outfan of this depositor connects the input of second depositor, the outfan of described second depositor connects the input of the 3rd depositor, thereby this second and the 3rd depositor receives and exports the data from first and second depositors respectively under the control of second data enable signal; Described first and second data enable signals produce and export by a timing control module.
In the such scheme, described interpolation circuit also comprises one the 4th subtractor, be used for to from described second and the data of the 3rd depositor ask poor, and the output that this difference DELTA is sent to accumulator and this accumulator is added up; Described interpolation circuit also comprises a multiplier, receive and calculate the long-pending of described the 3rd depositor dateout and constant N, output is toward second adder again, and is sent to a divider after the data summation from described accumulator, parameter divided by constant N and after exporting the interpolation of each point; Wherein, N is a natural number; When N is 2 power, described multiplier can be reduced to one first shift register, and log will move to left from the data of described the 3rd depositor
2Be sent to second adder after the operation of N position, described divider can be reduced to one second shift register, with the data of the described second adder output log that moves to right
2Parameter after the interpolation of the operation of N position and output each point.
In the such scheme, described interpolation circuit also comprises one the 4th subtractor, be used for to from described second and the data of the 3rd depositor ask poor, and the output that this difference DELTA is sent to accumulator and this accumulator is added up; Described interpolation circuit also comprises a divider, will give second adder after divided by N from the data of described accumulator output, and from parameter after the interpolation of the dateout summation back output each point of described the 3rd depositor; Wherein, N is a natural number; As N 2 power; Described divider can be reduced to a shift register, and log will move to right from the data of described accumulator output
2Give second adder after the operation of N position.
In the such scheme, described parameter calculation unit comprises: timing control module control produces collectiong focusing point coordinates value, and (Xi is Yi) with collectiong focusing degree of depth Fi; Channel counter is by counting, and order is read the array element coordinate Ex of each receive path correspondence from the array element coordinate storage, and Ey, first subtractor calculate the difference of described coordinate Xi and Ex, and second subtractor calculates the difference of coordinate Yi and Ey; First squaring circuit and second squaring circuit are squared to above-mentioned two differences respectively; First adder and sqrt circuit with these two square values summations after sqrt again, obtain described focus point and the distance that receives array element; The focusing delay distance after the recovery that the 3rd subtractor is obtained this distance and the described collectiong focusing degree of depth is sent to first divider, be divided by with the velocity of sound constant c that is sent to this divider, thereby this divider output corresponding focus on delay parameter Delay (i, j); Described timing control module control also produces and exports the accumulator reset signal toward each described interpolation circuit.
Adopt technique scheme, realize that the pointwise that focuses on delay parameter in the wave beam building-up process calculates in real time, can not produce interim delay time error, and it is simple to have a circuit, is easy to realize, realizes the low advantage of cost.
Description of drawings Fig. 1 is the theory diagram of existing typical ultrasonic image-forming system
Fig. 2 is that the present invention receives array element delay calculating principle sketch map
Fig. 3 is the time-delay coefficient curve bunch sketch map that collectiong focusing point of the present invention is tried to achieve when receiving aperture
Fig. 4 is the time-delay coefficient curve bunch sketch map that Fig. 3 data of the present invention obtain by interpolation for the basis
Fig. 5 is a first passage retardation coefficient difference change curve sketch map before and after the interpolation of the present invention
Fig. 6 is the structural representation of the real-time accountant parameter calculation unit of the present invention
Fig. 7 is the structural representation of the real-time accountant interpolation circuit of the present invention
Fig. 8 is the structural representation of the real-time accountant interpolation circuit of the present invention most preferred embodiment
Below the specific embodiment, the most preferred embodiment shown in is further set forth the present invention in conjunction with the accompanying drawings.
The present invention is that example has illustrated the computational methods that focus on delay with the protruding battle array (situation of linear array and phased array is suitable equally) of Fig. 2.Can be because receive focus not at the center of receiving aperture, for being without loss of generality, the reception focus shown in the figure is positioned at the left side of receiving aperture.
Then focus point with the focusing delay distance difference that receives array element is
Wherein, j is probe array element sequence number, (E_X
j, E_Y
j) be the array element coordinate, i is for receiving burnt period, (x
i, y
i) for receiving focal coordinates, F
iFor receiving focal depth.Because every initial focal coordinates (x that receives line
0, y
0), and the adjustment amount between two focuses (dx, dy) all known, therefore receive focal coordinates and be:
x
i＝x
0+i*dx y
i＝y
0+i*dy (2)
In order hard-wired the time, to avoid multiplication, can replace following formula to finish the calculating that receives focus with the method that adds up one by one, shown in (3) formula:
x
i+1＝x
i+dx y
i+1＝y
i+dy (3)
In like manner define F and represent the focus spacing, then
F
i+1＝F
i+F (4)
In (4) formula, when i=0, F
0=0.
After the probe type is given, array element coordinate (E_X
j, E_Y
j) just determined that also focus spacing F also can pre-set, therefore can calculate each and receive the delay inequality of focus for each array element, the time-delay that obtains each array element then is
Delay(i，j)＝distance_diff(i，j)/c (5)
Wherein c is the velocity of sound.
In said method, the retardation coefficient of the array element of a passage correspondence of calculating need be done 4 sub-additions, and (subtraction is treated as and is addition, for each passage, receives focal coordinates Xi, and the calculating of Yi and focal length Fi is the same; Once just passable in each focus point calculating, its amount of calculation does not count), twice square, and once extract square root; And the retardation coefficient of the array element of each passage correspondence all needs to calculate, if the retardation coefficient of each focus point all calculates like this in the receiving course, then each passage all needs the such circuit of a cover, for multi-beam, more need increase circuit at double, spending of hardware is very big.In order to reduce the consumption of hardware resource, need take easier method.
Fig. 3 has illustrated collectiong focusing point for (0,20), (0,44.4141), (0,68.8281), (0,93.2422), (0,117.6562), (0,142.0703), the time-delay coefficient curve of calculating with formula (1) when (0,166.4844) bunch.As seen along with the increase that receives focal length, has the linear trend that reduces the time delay of the corresponding array element of each passage, thereby the method that can adopt interpolation is asked for the retardation coefficient between the collectiong focusing point that N is ordered, the method of interpolation comprises spline interpolation, interpolation methods such as linear interpolation, this method embodiment is that example is asked for the retardation coefficient between the focus point that N is ordered with the linear interpolation, but the present invention can also adopt other interpolation method.Described linear interpolation relies on following formula and carries out:
Delay(k，j)＝(Delay(i，j)*N+Δ*k)/N
(6)
Wherein Δ=Delay (i+N, j)-Delay (i, j)
Wherein, j represents is to receive array element number, and k is the collectiong focusing point between from i to i+N, when the N value is 2 power, multiplies each other with N and can be reduced to shifting function divided by N on hardware is realized, Δ * k can realize by the mode that adds up.Through type (1) calculate Delay (i, j) and Delay (i+N, j) after, for each passage, when the retardation coefficient of real-time each acceptance point of calculating, only need an accumulator, an adder, twice shifting function is just passable, compares with formula (1) with direct, and hardware spending greatly reduces.
Fig. 4 has provided when interpolation and has counted when being 6, with seven delay curves among Fig. 3 respectively as starting point and terminal point, six groups of (every group of 6 curves) retardation coefficient set of curves that obtain after the interpolation.As seen from Figure 4, the retardation coefficient that each group interpolation obtains all is evenly distributed between the virgin curve, hence one can see that, when the distance that reduces between the described virgin curve that is used for interpolation, increase N simultaneously to calculating in real time till each acceptance point retardation coefficient, with the retardation coefficient set of curves that obtains distributing more intensive.
Whether the retardation coefficient that adopts the method for interpolation to ask between consecutive points is fit to, and the difference between the time delay that obtains after time delay that need be by assessment actual reception point and the interpolation is judged.Fig. 5 adopts interpolation algorithm (to be example at a distance of points N=64) time delay that obtains and the difference of actual delay time with the change in depth curve for each acceptance point in the first passage.As seen, except the starting point that interpolation is used, the time delay that all the other each points calculate is all variant, but described difference is all less than 1ns.When the echo-signal sample rate was 40M, thick delay precision was 25ns, and the delay precision of setting after carefully delaying time is 6.25ns, so the delay precision that adopts interpolation to obtain satisfies hard-wired needs.The delay variance of first passage is maximum, so the error of other passages is littler as can be known, thereby the retardation coefficient that adopts linear interpolation to ask for intermediate point is feasible, and the principle that the interpolation points N is chosen is to guarantee the precision of interpolation error less than the thin time-delay of hardware.
Thereby the method for the real-time compute beam synthetic focusing of the present invention delay parameter is can brief summary as follows, comprises step:
A. will be sent to beam synthesizer from the ultrasonic echo electrical signal data of a plurality of receive paths of same received beam;
B. the data to these passages postpone to handle with the corresponding focusing parameter that postpones respectively;
C. the synthetic ultrasonic echo data that obtains a scanning line;
Especially, the described delay focusing parameter among the step B is that each acceptance point at each receive path calculates in real time and provides, thereby reduces the delay error of near-field region.This method is applicable to various port numbers, includes but not limited to the supersonic imaging apparatus of 32,64 or 128 passages; Not only, can also be applicable to the multi-beam receiving system applicable to the simple beam receiving system.
Ultrasonic image-forming system with one 64 passage simple beam is example (carrying out need carrying out the delay parameter calculating of maximum 64 passages when wave beam is synthetic) below, and the device of realizing this method is described.
Apparatus of the present invention comprise a parameter calculation unit as shown in Figure 6 as the part of beam synthesizer, according to element position calculate a predetermined acceptance point focusing delay parameter Delay (i, j).This unitary timing control module control generation collectiong focusing point coordinates value (Xi, Yi), collectiong focusing degree of depth Fi; Channel counter is by counting, and order is read the array element coordinate Ex of each receive path correspondence from the array element coordinate storage, and Ey, first subtractor calculate the difference of described coordinate Xi and Ex, and second subtractor calculates the difference of coordinate Yi and Ey; First squaring circuit and second squaring circuit are squared to above-mentioned two differences respectively; First adder and sqrt circuit with these two square values summations after sqrt again, obtain described focus point and the distance that receives array element; The focusing delay distance after the recovery that the 3rd subtractor is obtained this distance and the described collectiong focusing degree of depth is sent to first divider, be divided by with the velocity of sound constant c that is sent to this divider, thereby this divider output corresponding focus on delay parameter Delay (i, j).
Apparatus of the present invention also comprise a plurality of interpolation circuits, and each described interpolation circuit is corresponding to a receive path, the calculating of perfect in this passage (6).Therefore for 64 passages, need 64 interpolation circuits.
Among the present invention embodiment as shown in Figure 7, described interpolation circuit comprises three depositors, and first depositor is used for receiving and keeping in the focusing delay parameter corresponding to this passage of described parameter calculation unit output under the control of first data enable signal; The outfan of this depositor connects the input of second depositor, the outfan of described second depositor connects the input of the 3rd depositor, thereby this second and the 3rd depositor receives and exports the data from first and second depositors respectively under the control of second data enable signal; Described first and second data enable signals are all exported from described timing control circuit.
In the initial work process before beam synthesizer begins to receive, described parameter calculation unit is calculated earlier the delay parameter of each passage on predetermined collectiong focusing point in 64 passages one by one, for example Delay (0, j), wherein (0≤j＜64), and the effect by described first data enable signal, these parameters are deposited in first depositor of respective channel interpolation circuit; System is by the control of described second data enable signal, and these parameters are restored in each described second depositor; Described parameter calculation unit calculate one by one again each passage delay parameter Delay (N, j), and same as above depositing in each described first depositor; System is the control by described second data enable signal once more, make each passage parameter Delay (0, j) enter the 3rd depositor, parameter Delay (N by second depositor, j) enter second depositor by first depositor after, the initial work of apparatus of the present invention is finished.
Described interpolation circuit also comprises one the 4th subtractor, be used for to from described second and the data of the 3rd depositor ask poor, this difference DELTA is sent to accumulator, adding up with the output of this accumulator replaces carrying out the mlultiplying circuit that Δ * k calculates; In order to save resource, interval N between two interpolation points of present embodiment is taken as 2 power, thereby multiplication in (6) and division are replaced with shifting function, and therefore described interpolation circuit also comprises first shift register, and log will move to left from the data of the 3rd depositor
2Be sent to second adder after the operation of N position, and from being sent to second shift register after the summation of the data of described accumulator, log moves to right
2The operation of N position.After thereby beam synthesizer begins to receive, i=0, in the receiving course of 0 to N-1 acceptance point, point of every reception, described accumulator adds up once, parameter after the interpolation of a point of described second shift register output.
Meanwhile, described parameter calculation unit is calculated the delay parameter Delay of each passage one by one (2*N j), and is deposited in each first depositor by the control to described first data enable signal; After acceptance point changes to the N point, system is by the control to described second data enable signal, with described second and the data of the 3rd depositor be updated to the current data Delay (2*N of described first and second depositors respectively, j) and Delay (N, j), simultaneously, described timing control circuit produces an accumulator reset signal and is sent to described accumulator this accumulator that resets, so enter reception and the real-time interpolation calculation process of delay parameter that N is ordered to 2*N-1.
The rest may be inferred, in receiving course, described parameter calculation unit always constantly calculate in advance 64 each passages of passage following N point place delay parameter and be temporarily stored in each first depositor, upgrade described second and the value of the 3rd depositor before coming to begin for the next round interpolation process.So circulate, just can constantly obtain the retardation coefficient value of each acceptance point of each passage in real time.
Be noted that described interpolation circuit can replace first shift register with multiplier, receive and calculate the long-pending of described the 3rd depositor dateout and constant N, output is toward described second adder again; In the time of can replacing second shift register with divider, the dateout of described second adder divided by constant N, is obtained the retardation coefficient value of each acceptance point.Two shift registers in Fig. 7 interpolation circuit only are that described multiplier and divider are substituting under 2 the power condition at N, can simplify circuit.In addition, replacing accumulator to realize Δ * k with multiplier, also is a kind of replacement scheme that is equal to.
In order further to reduce the complexity of circuit, the present invention can further be reduced to formula (6)
Delay(k，j)＝Delay(i，j)+Δ*k/N
Wherein Δ=Delay (i+N, j)-Delay (i, j)
Thereby, shown in the most preferred embodiment of Fig. 8, described interpolation circuit can reduce by a shift register, only need a subtractor, an accumulator (or multiplier), adder and a shift register (simplification of divider when N is 2 power) to get final product, be specially: described accumulator dateout is through the shift register log that moves to right
2Give second adder after the operation of N position, after suing for peace, export the retardation coefficient value of each acceptance point from the dateout of described the 3rd depositor.
Claims (11)
1. the method for a real-time compute beam synthetic focusing delay parameter, the wave beam that is used for supersonic imaging apparatus is synthetic, comprises step:
A. will be sent to beam synthesizer from the ultrasonic echo electrical signal data of a plurality of receive paths of same received beam;
B. the data to these passages postpone to handle with the corresponding focusing parameter that postpones respectively;
C. the synthetic ultrasonic echo data that obtains a scanning line;
It is characterized in that the described delay focusing parameter among the step B is that each acceptance point at each receive path calculates in real time and provides, thereby reduce the delay error of near-field region.
2. according to the method for the described real-time compute beam synthetic focusing delay parameter of claim 1, it is characterized in that described real-time calculating adopts interpolation processing to calculate, and comprises step:
A. to the real-time respectively delay focusing parameter that calculates on putting two times of reception of each receive path; If between these two acceptance points N-1 acceptance point arranged, N is a natural number;
B. postponing focusing parameter with these two simultaneously is that starting point and terminal point carry out interpolation calculation, obtains the real time delay focusing parameter of N-1 interpolation point, makes a corresponding described N-1 acceptance point respectively.
3. according to the method for the described real-time compute beam synthetic focusing delay parameter of claim 2, it is characterized in that:
Described interpolation processing adopts linear interpolation method, thereby is that to postpone focusing parameters with described two be that starting point and terminal point carry out equally spaced linear interpolation and calculate among the step b.
4. according to the method for the described real-time compute beam synthetic focusing delay parameter of claim 3, it is characterized in that the real-time results of two delay focusing parameters of putting the time of reception are obtained by following formula respectively described in the step a
Delay(i，j)＝distance_diff(i，j)/c；
Wherein, c is the velocity of sound,
Wherein, j is the probe array element sequence number of described passage correspondence, (E_X
j, E_Y
j) be this array element coordinate, i is for receiving burnt period, (x
i, y
i) for receiving focal coordinates, F
iFor receiving focal depth.
5. according to the method for the described real-time compute beam synthetic focusing delay parameter of claim 3, it is characterized in that:
The described interpolation N-1 that counts is decided by that device hardware realizes the precision of thin time-delay, and the interpolation error of being brought by described interpolation calculation; The former is greater than the latter.
6. according to the method for the described real-time compute beam synthetic focusing delay parameter of claim 5, it is characterized in that:
Described interval N is 2 power.
7. the device of a real-time compute beam synthetic focusing delay parameter, the wave beam synthesis unit that is used for supersonic imaging apparatus, this device comprises a parameter calculation unit, is used for calculating the focusing delay parameter of being scheduled to acceptance point according to the probe element position, it is characterized in that:
This device also comprises a plurality of interpolation circuits, and each described interpolation circuit is corresponding to a receive path.
8. according to the device of the described real-time compute beam synthetic focusing delay parameter of claim 7, it is characterized in that:
Described interpolation circuit comprises three depositors, and first depositor is used for receiving and keeping in the focusing delay parameter corresponding to this passage of described parameter calculation unit output under the control of first data enable signal; The outfan of this depositor connects the input of second depositor, the outfan of described second depositor connects the input of the 3rd depositor, thereby this second and the 3rd depositor receives and exports the data from first and second depositors respectively under the control of second data enable signal; Described first and second data enable signals produce and export by a timing control module.
9. the device of described according to Claim 8 real-time compute beam synthetic focusing delay parameter is characterized in that:
Described interpolation circuit also comprises one the 4th subtractor, be used for to from described second and the data of the 3rd depositor ask poor, and the output that this difference DELTA is sent to accumulator and this accumulator is added up; Described interpolation circuit also comprises a multiplier, receive and calculate the long-pending of described the 3rd depositor dateout and constant N, output is toward second adder again, and is sent to a divider after the data summation from described accumulator, parameter divided by constant N and after exporting the interpolation of each point; Wherein, N is a natural number;
When N is 2 power, described multiplier can be reduced to one first shift register, and log will move to left from the data of described the 3rd depositor
2Be sent to second adder after the operation of N position, described divider can be reduced to one second shift register, with the data of the described second adder output log that moves to right
2Parameter after the interpolation of the operation of N position and output each point.
10. the device of described according to Claim 8 real-time compute beam synthetic focusing delay parameter is characterized in that:
Described interpolation circuit also comprises one the 4th subtractor, be used for to from described second and the data of the 3rd depositor ask poor, and the output that this difference DELTA is sent to accumulator and this accumulator is added up; Described interpolation circuit also comprises a divider, will give second adder after divided by N from the data of described accumulator output, and from parameter after the interpolation of the dateout summation back output each point of described the 3rd depositor; Wherein, N is a natural number;
As N 2 power; Described divider can be reduced to a shift register, and log will move to right from the data of described accumulator output
2Give second adder after the operation of N position.
11. the device according to the described real-time compute beam synthetic focusing delay parameter of claim 7 is characterized in that,
Described parameter calculation unit comprises: timing control module control produces collectiong focusing point coordinates value, and (Xi is Yi) with collectiong focusing degree of depth Fi; Channel counter is by counting, and order is read the array element coordinate Ex of each receive path correspondence from the array element coordinate storage, and Ey, first subtractor calculate the difference of described coordinate Xi and Ex, and second subtractor calculates the difference of coordinate Yi and Ey; First squaring circuit and second squaring circuit are squared to above-mentioned two differences respectively; First adder and sqrt circuit with these two square values summations after sqrt again, obtain described focus point and the distance that receives array element; The focusing delay distance after the recovery that the 3rd subtractor is obtained this distance and the described collectiong focusing degree of depth is sent to first divider, be divided by with the velocity of sound constant c that is sent to this divider, thereby this divider output corresponding focus on delay parameter Delay (i, j);
Described timing control module control also produces and exports the accumulator reset signal toward each described interpolation circuit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNA2007100766203A CN101371789A (en) | 2007-08-22 | 2007-08-22 | Real-time computing method and device for focalization delay parameters of beam-forming |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNA2007100766203A CN101371789A (en) | 2007-08-22 | 2007-08-22 | Real-time computing method and device for focalization delay parameters of beam-forming |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN101371789A true CN101371789A (en) | 2009-02-25 |
Family
ID=40446306
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNA2007100766203A Pending CN101371789A (en) | 2007-08-22 | 2007-08-22 | Real-time computing method and device for focalization delay parameters of beam-forming |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101371789A (en) |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101858972A (en) * | 2010-03-23 | 2010-10-13 | 深圳市蓝韵实业有限公司 | Multi-beam synthesis method and device based on delay parameter real-time calculation and assembly line |
| CN102770079A (en) * | 2010-02-23 | 2012-11-07 | 佳能株式会社 | Ultrasonic imaging apparatus and method of controlling delay |
| CN103654853A (en) * | 2013-11-19 | 2014-03-26 | 深圳先进技术研究院 | Base band beam forming ultrasonic imaging method and system |
| CN104349719A (en) * | 2012-05-25 | 2015-02-11 | 富士胶片株式会社 | Ultrasound diagnostic device, and data processing method |
| CN104434218A (en) * | 2014-12-15 | 2015-03-25 | 飞依诺科技(苏州)有限公司 | Real-time calculation method and device of focusing delays for ultrasonic beam-forming |
| CN104655728A (en) * | 2015-01-29 | 2015-05-27 | 中国科学院声学研究所 | Acoustic phased array imaging method |
| CN106529121A (en) * | 2016-09-30 | 2017-03-22 | 飞依诺科技(苏州)有限公司 | Ultrasonic dynamic focus delay data recovery method and recovery system |
| CN107948573A (en) * | 2017-11-22 | 2018-04-20 | 深圳市华星光电技术有限公司 | A kind of linear interpolated value method and device of digital signal |
| CN109087626A (en) * | 2018-08-08 | 2018-12-25 | 京东方科技集团股份有限公司 | Phased array transmit circuit, launching technique, phased array device and preparation method thereof |
| US10217217B2 (en) | 2015-01-20 | 2019-02-26 | Indian Institute Of Technology Bombay | Systems and methods for obtaining 3-D images from X-ray information |
| CN109975814A (en) * | 2017-12-28 | 2019-07-05 | 深圳先进技术研究院 | Ultrasonic imaging method, system and equipment |
| CN110613476A (en) * | 2018-06-19 | 2019-12-27 | 青岛海信医疗设备股份有限公司 | Ultrasonic signal processing method and device |
| CN110720947A (en) * | 2018-07-16 | 2020-01-24 | 青岛海信医疗设备股份有限公司 | Control method and control device for delayed signal transmission of ultrasonic probe |
| CN113180730A (en) * | 2021-03-31 | 2021-07-30 | 上海深至信息科技有限公司 | Distributed beam imaging method |
-
2007
- 2007-08-22 CN CNA2007100766203A patent/CN101371789A/en active Pending
Cited By (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102770079A (en) * | 2010-02-23 | 2012-11-07 | 佳能株式会社 | Ultrasonic imaging apparatus and method of controlling delay |
| CN102770079B (en) * | 2010-02-23 | 2015-03-11 | 佳能株式会社 | Ultrasonic imaging apparatus and method of controlling delay |
| US9165552B2 (en) | 2010-02-23 | 2015-10-20 | Canon Kabushiki Kaisha | Ultrasonic imaging apparatus and method of controlling delay |
| CN101858972A (en) * | 2010-03-23 | 2010-10-13 | 深圳市蓝韵实业有限公司 | Multi-beam synthesis method and device based on delay parameter real-time calculation and assembly line |
| CN101858972B (en) * | 2010-03-23 | 2013-01-30 | 深圳市蓝韵实业有限公司 | Multi-beam synthesis method and device based on delay parameter real-time calculation and assembly line |
| CN104349719A (en) * | 2012-05-25 | 2015-02-11 | 富士胶片株式会社 | Ultrasound diagnostic device, and data processing method |
| CN103654853A (en) * | 2013-11-19 | 2014-03-26 | 深圳先进技术研究院 | Base band beam forming ultrasonic imaging method and system |
| CN104434218A (en) * | 2014-12-15 | 2015-03-25 | 飞依诺科技(苏州)有限公司 | Real-time calculation method and device of focusing delays for ultrasonic beam-forming |
| US10217217B2 (en) | 2015-01-20 | 2019-02-26 | Indian Institute Of Technology Bombay | Systems and methods for obtaining 3-D images from X-ray information |
| CN104655728B (en) * | 2015-01-29 | 2017-09-19 | 中国科学院声学研究所 | A kind of acoustics phased array imaging method |
| CN104655728A (en) * | 2015-01-29 | 2015-05-27 | 中国科学院声学研究所 | Acoustic phased array imaging method |
| CN106529121B (en) * | 2016-09-30 | 2019-01-29 | 飞依诺科技(苏州)有限公司 | The restoration methods and recovery system of ultrasonic dynamic focusing delay data |
| CN106529121A (en) * | 2016-09-30 | 2017-03-22 | 飞依诺科技(苏州)有限公司 | Ultrasonic dynamic focus delay data recovery method and recovery system |
| CN107948573B (en) * | 2017-11-22 | 2020-11-03 | 深圳市华星光电技术有限公司 | Digital signal linear interpolation method and device |
| CN107948573A (en) * | 2017-11-22 | 2018-04-20 | 深圳市华星光电技术有限公司 | A kind of linear interpolated value method and device of digital signal |
| CN109975814A (en) * | 2017-12-28 | 2019-07-05 | 深圳先进技术研究院 | Ultrasonic imaging method, system and equipment |
| CN109975814B (en) * | 2017-12-28 | 2020-09-22 | 深圳先进技术研究院 | Ultrasound imaging method, system and device |
| CN110613476A (en) * | 2018-06-19 | 2019-12-27 | 青岛海信医疗设备股份有限公司 | Ultrasonic signal processing method and device |
| CN110720947A (en) * | 2018-07-16 | 2020-01-24 | 青岛海信医疗设备股份有限公司 | Control method and control device for delayed signal transmission of ultrasonic probe |
| CN109087626A (en) * | 2018-08-08 | 2018-12-25 | 京东方科技集团股份有限公司 | Phased array transmit circuit, launching technique, phased array device and preparation method thereof |
| CN109087626B (en) * | 2018-08-08 | 2023-08-25 | 京东方科技集团股份有限公司 | Phased array transmitting circuit, transmitting method, phased array device and preparation method thereof |
| CN113180730A (en) * | 2021-03-31 | 2021-07-30 | 上海深至信息科技有限公司 | Distributed beam imaging method |
| CN113180730B (en) * | 2021-03-31 | 2022-10-18 | 上海深至信息科技有限公司 | Distributed beam imaging method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101371789A (en) | Real-time computing method and device for focalization delay parameters of beam-forming | |
| CN101410061B (en) | Delay controller for ultrasound receive beamformer | |
| CN103796594B (en) | Ultrasonic imaging apparatus | |
| US7800979B2 (en) | Beamforming method and apparatus used in ultrasonic imaging system | |
| US6231511B1 (en) | Ultrasonic signal focusing method and apparatus for ultrasonic imaging system | |
| US8545406B2 (en) | Dynamic aperture control and normalization for apodization in beamforming | |
| JPH01207042A (en) | Method and apparatus for reducing repeated adaptation of phase differential effect | |
| DE69628810T2 (en) | Real-time distributed calculation arrangement for beam delays in an ultrasound imaging system | |
| CN103547220A (en) | Beamforming method and diagnostic ultrasound apparatus | |
| CN103033816B (en) | Synthetic aperture focused ultrasonic imaging implementation method based on arc scanning transition | |
| CN101900808B (en) | Delay parameter computation method for multi-channel scanning and device thereof | |
| EP1060470B1 (en) | Method and apparatus for distributed, agile calculation of beamforming time delays and apodization values | |
| CN105832366B (en) | A kind of delay real-time computing technique for during Beam synthesis | |
| JPH0870404A (en) | Delay causing device in ultrasonic beam forming apparatus | |
| JPH024290B2 (en) | ||
| EP0256282B1 (en) | Ultrasonic wave diagnostic apparatus employing interpolated values of weighting data | |
| JP2575603B2 (en) | Scanning device | |
| US8926514B2 (en) | Iterative time delay values for ultrasound beamforming | |
| US5501219A (en) | Real-time dynamic time-of-flight calculator | |
| US20080092660A1 (en) | Multi-line beamforming extention using sub-arrays | |
| CN101961251B (en) | Method and device for computing apodization curve in real time in medical ultrasonic diagnosis system | |
| US20070239013A1 (en) | Delay controller for ultrasound receive beamformer | |
| EP0022966A2 (en) | Ultrasonic diagnosis apparatus | |
| JPS60127457A (en) | Ultrasonic medium characteristic value measuring apparatus | |
| CN104142608B (en) | Nonlinear dynamic focusing control method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C12 | Rejection of a patent application after its publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20090225 |