CN101645167B - Method and device for frame correlation coefficient control and frame correlation processing, and ultrasonic imaging system - Google Patents
Method and device for frame correlation coefficient control and frame correlation processing, and ultrasonic imaging system Download PDFInfo
- Publication number
- CN101645167B CN101645167B CN200810142787XA CN200810142787A CN101645167B CN 101645167 B CN101645167 B CN 101645167B CN 200810142787X A CN200810142787X A CN 200810142787XA CN 200810142787 A CN200810142787 A CN 200810142787A CN 101645167 B CN101645167 B CN 101645167B
- Authority
- CN
- China
- Prior art keywords
- frame
- frequency response
- filter amplitude
- difference
- frame rate
- 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.)
- Active
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a method and a device for frame correlation coefficient control and frame correlation processing, and an ultrasonic imaging system. The method for the frame correlation coefficient control comprises the following steps: comparing the amplitude-frequency response of a current filter at a current imaging frame rate with the amplitude-frequency response of a reference filter at a reference imaging frame rate; acquiring the difference between the amplitude-frequency response of the current filter and the amplitude-frequency response of the reference filter in an overlapped frequency band; and judging the frame correlation coefficient corresponding to the amplitude-frequency response of the current filter when the difference is minimum to be the best frame correlation coefficient corresponding to the current imaging frame rate. The method and the device for the frame correlation coefficient control and the frame correlation processing can obtain the frame correlation coefficient accurately, make the frame correlation processing results of images very close to each other when the imaging frame rate changes dramatically, ensure the noise filtering effect of an ultrasonic image, and improve the signal-to-noise ratio of the image effectively.
Description
Technical field
The present invention relates to the technical field of optimizing of ultrasonoscopy time domain smothing filtering coefficient in the ultrasonic system, in particular, a kind of control method of frame correlation coefficient and device thereof and the frame correlation processing method and the ultrasonic image-forming system that utilize this method.
Background technology
Ultrasonoscopy is owing to exist bigger circuit noise and spot (Speckle) noise (wherein; Speckle noise is owing to the scattered signal coherence stack of taking measurements greatly less than the particle of ultrasound wave wavelength forms; This is intrinsic problem in the ultrasonoscopy system); So need smoothing processing, to improve the signal to noise ratio (S/N ratio) of image through some time domains and spatial domain.Wherein, the smoothing processing of time domain is also claimed the frame relevant treatment usually.The mode that frame relevant treatment commonly used adopts previous frame output image and present frame image to do weighted sum realizes that be a kind of simple recursive filtering method, it can be expressed as:
y(t)＝ay(t-1)+(1-a)x(t) (1)
Wherein, the image of y (t) for exporting after the present frame frame relevant treatment, x (t) is the present frame image, and y (t-1) is the output image of previous frame, and a is a frame correlation coefficient.
In fact, the frame relevant treatment belongs to a kind of filtering method, and it can be described as a first order IIR (Infinite Impulse Response) digital filter in the signal Processing field, and the transfer function H (z) in this digital filter Z territory can be expressed as:
With z=exp{j2 π f ＇ }, the above-mentioned H of substitution (z) promptly obtains the frequency response H (f ＇) of wave filter.Generally, only consider the amplitude-frequency response of this wave filter during the frame relevant treatment, i.e. the amplitude of said frequencies response | H (f) |, wherein, f ＇ is the normalization numerical frequency.
Frequency response according to above-mentioned wave filter can know that this wave filter is a low-pass filter, because the circuit noise of different frame is normally independent; And the speckle noise of different frame moving owing to tissue or probe; Cause the speckle noise of different frame also to have more weak correlativity, and the image of the tissue of different frame or blood flow has stronger correlativity, therefore can effectively suppress the noise in the image through this wave filter; And the enhancing useful signal, thereby improve the contrast resolution of image.Can realize squelch in various degree through different frame correlation coefficient a is set, along with the increase of frame correlation coefficient a, noise suppression effect is remarkable more, but also causes the fuzzy of motor tissue image easily.Same frame correlation coefficient a shows different squelch and motor tissue blur effect under the pictures different imaging frame rate; As when imaging frame rate is higher; Bigger frame correlation coefficient a is set then can reduces picture noise effectively; And can not cause motion blur, but when imaging frame rate was lower, same frame correlation coefficient a then can cause tangible motion blur.
In the existing compuscan, normal frame correlation coefficient a is set to the exponential function of imaging frame rate, and it specifically is expressed as
a＝exp{-1/k×Fr} (3)
Wherein, k is the relevant control coefrficient of frame, and Fr is an imaging frame rate.It is thus clear that along with the reduction of imaging frame rate Fr, frame correlation coefficient a also decreases, improved the imaging frame rate Fr influence level and smooth to a certain extent to the image time domain.
The frame correlation coefficient a that formula (3) is represented is updated in the transport function in the represented digital filter Z territory of aforementioned formula (2),
Wherein,
Can obtain corresponding wave filter amplitude-frequency response | H (f) |.As an example; Under the different imaging frame rates, the relevant pairing wave filter amplitude-frequency response of control coefrficient k of same number of frames is as shown in Figure 1, and it has provided the relevant control coefrficient k of frame respectively is 0.02; Frame per second is 15Hz; 30Hz is under the situation of 45Hz and 60Hz, through four amplitude-frequency responses of said method acquisition.The frequency range that the amplitude-frequency response of real coefficient digital filter is corresponding is zero to 1/2nd of current imaging frame rate, i.e. f=0～Fr
t/ 2; Visible by Fig. 1; Amplitude-frequency response under the different frame per second is along with the reduction institute acquisition wave filter of the frame per second inhibition ability to high frequency noise also decreases, when frame per second changes more greatly the noise suppression effect difference also more greatly, therefore when frame per second changes greatly; The high frequency attenuation of frame relevant treatment obviously descends, thereby has reduced the relevant inhibiting effect to speckle noise and circuit noise of frame.
According to above-mentioned analysis, imaging frame rate Fr change little in, adopt above-mentioned formula (3) simply to find the solution frame correlation coefficient a when carrying out the frame relevant treatment; Its frame relevant effect is close, still, and when imaging frame rate Fr changes greatly; The situation that the noise suppression effect variation just possibly occur; Switch to the color flow angiography pattern such as compuscan from organizing gray scale (B-mode) imaging pattern, be generally image-to-frame ratio Fr and drop to several frames, organize this moment the noise of gray level image to increase suddenly from tens frames; Carry out the frame relevant treatment if at this moment also adopt above-mentioned formula (3) simply to find the solution frame correlation coefficient a; Just can't guarantee the frame relevant effect of image, image fault can occur, reduce the filtration result of picture noise.
It is thus clear that also there is certain drawback in the method that obtains frame correlation coefficient a in the prior art, need to improve further.
Summary of the invention
The frame correlation processing method and the ultrasonic image-forming system that the object of the present invention is to provide a kind of control method and device thereof of frame correlation coefficient and utilize this method; It can obtain frame correlation coefficient; And when imaging frame rate alters a great deal; The frame correlation coefficient that uses the present invention to obtain can make the frame relevant treatment result of image very approaching, guarantees the ultrasonic image noise filtration result, and the signal to noise ratio (S/N ratio) of effectively improving image.
To achieve these goals, the present invention adopts following technical scheme:
The invention provides a kind of control method of frame correlation coefficient; Said method comprises: the reference filter amplitude-frequency response under a pre-filter amplitude-frequency response and a benchmark imaging frame rate under A, the more current imaging frame rate obtains said when pre-filter amplitude-frequency response and the difference of reference filter amplitude-frequency response in overlapping bands; B, with said difference be hour when the pairing frame correlation coefficient of pre-filter amplitude-frequency response, be judged to be the pairing optimum frame related coefficient of current imaging frame rate.
On the basis of said method, the present invention also provides a kind of frame correlation coefficient controller, and said controller comprises: first module is used for calculating the frame correlation coefficient under the benchmark imaging frame rate according to preset benchmark imaging frame rate; Unit second, the frame correlation coefficient that is used for obtaining according to first module calculates the reference filter amplitude-frequency response; And Unit the 3rd, be used for more current imaging frame rate under when pre-filter amplitude-frequency response and the difference of said reference filter amplitude-frequency response in overlapping bands, that exports said difference and be hour works as the pairing frame correlation coefficient of pre-filter amplitude-frequency response.
On the basis of said method, the present invention also provides a kind of ultrasonoscopy frame correlation processing method, said method comprising the steps of:
The reference filter amplitude-frequency response under a pre-filter amplitude-frequency response and a benchmark imaging frame rate under S1, the more current imaging frame rate obtains said when pre-filter amplitude-frequency response and the difference of reference filter amplitude-frequency response in overlapping bands;
S2, with said difference be hour when the pairing frame correlation coefficient of pre-filter amplitude-frequency response, be judged to be the pairing optimum frame related coefficient of current imaging frame rate;
The optimum frame related coefficient that S3, utilization obtain is carried out the frame related operation, the current frame image after obtaining to handle.
To said method, the present invention also provides a kind of ultrasonic image-forming system, and said system comprises: ultrasonoscopy processing module and the probe that is used to obtain ultrasonoscopy; Be provided with a frame correlation coefficient controller and a frame relevant treatment unit in the said ultrasonoscopy processing module, this controller comprises: first module is used for calculating the frame correlation coefficient under the benchmark imaging frame rate according to preset benchmark imaging frame rate; Unit second, the frame correlation coefficient that is used for obtaining according to first module calculates the reference filter amplitude-frequency response; And Unit the 3rd, be used for more current imaging frame rate under when pre-filter amplitude-frequency response and the difference of said reference filter amplitude-frequency response in overlapping bands, that exports said difference and be hour works as the pairing frame correlation coefficient of pre-filter amplitude-frequency response; Said frame relevant treatment unit is used to receive the frame correlation coefficient of said Unit the 3rd output, and utilizes this frame correlation coefficient to carry out the frame related operation, the current frame image after obtaining to handle.
In addition, the present invention also provides a kind of ultrasonic image-forming system, and said system comprises: ultrasonoscopy processing module and the probe that is used to obtain ultrasonoscopy; Be provided with a storage medium in the said ultrasonoscopy processing module, search a unit and a frame relevant treatment unit; Said storage medium is preserved a relation table; This table record has optimum frame related coefficient under current imaging frame rate, the current imaging frame rate, and characterizes the corresponding relation between the relevant control coefrficient three of frame of frame relevant effect, the difference minimum of reference filter amplitude-frequency response in overlapping bands under pre-filter amplitude-frequency response and a benchmark imaging frame rate that said optimum frame related coefficient is corresponding; Said frame relevant treatment unit links to each other with said storage medium through the said unit of searching; The said unit of searching is used for searching said relation table and current imaging frame rate and the corresponding optimum frame related coefficient of the relevant control coefrficient of present frame, and this optimum frame related coefficient is sent in the said frame relevant treatment unit; Said frame relevant treatment unit is used to utilize said optimum frame related coefficient to carry out the frame related operation, the current frame image after obtaining to handle.
Can know from technique scheme; The present invention is through setting up the minimum difference principle of reference filter amplitude-frequency response in overlapping bands under a pre-filter amplitude-frequency response and a benchmark imaging frame rate; Select best frame correlation coefficient, in order to the frame relevant treatment of ultrasonoscopy, in terms of existing technologies; The present invention adopt reverse thinking from the frame correlation coefficient of the effect calculating optimum that will obtain; Thereby can obtain the frame correlation coefficient more suitable, optimize ultrasonoscopy time domain smothing filtering coefficient, feasible frame correlation processing method and the ultrasonic image-forming system thereof that utilizes this method than prior art; Have better filter effect and when imaging frame rate alters a great deal, also can guarantee the display effect of image and make the frame relevant treatment result of image very approaching.
Description of drawings
Fig. 1 is the wave filter amplitude-frequency response characteristic figure that prior art obtains;
Fig. 2 is the process flow diagram of frame correlation coefficient control method of the present invention;
Fig. 3 is the method flow diagram of the embodiment of the invention 1 to 3;
Fig. 4 be the present invention to utilize the relevant control coefrficient k of the frame that method obtained of embodiment 1 be 0.02, frame per second is followed successively by the pairing amplitude-frequency response characteristic figure of 15Hz, 30Hz, 45Hz and 60Hz;
Fig. 5 be the present invention to utilize the relevant control coefrficient k of the frame that method obtained of embodiment 4 be 0.02, frame per second is followed successively by the pairing amplitude-frequency response characteristic figure of 15Hz, 30Hz, 45Hz and 60Hz;
Fig. 6 is the structural representation of frame correlation coefficient controller of the present invention;
Fig. 7 is the inner structure synoptic diagram of Unit the 3rd among Fig. 6;
Fig. 8 is the structural representation that relates to image processing section in a kind of ultrasonic image-forming system provided by the invention;
Fig. 9 is the structural representation that relates to image processing section in the another kind of ultrasonic image-forming system provided by the invention.
Embodiment
For when imaging frame rate alters a great deal, can guarantee that also the frame relevant treatment result of image is very approaching, the present invention has provided a kind of control method of frame correlation coefficient, its accurate Calculation the frame correlation coefficient under the current frame rate.As shown in Figure 2, the main process of the frame correlation coefficient control method that the present invention is given is following:
At first, more current imaging frame rate Fr
tUnder when the pre-filter amplitude-frequency response | H (f) | with a benchmark imaging frame rate Fr
0Under the reference filter amplitude-frequency response | H
0(f) |, obtain when the pre-filter amplitude-frequency response | H (f) | with the reference filter amplitude-frequency response | H
0(f) | the difference E in overlapping bands, promptly E=F (| H (f) |, | H
0(f) |), F representes difference function;
Then, with said difference E be hour when the pre-filter amplitude-frequency response | H (f) | pairing frame correlation coefficient is judged to be current imaging frame rate Fr
tPairing optimum frame related coefficient.
It is thus clear that; The present invention on the basis of existing technology; Provided a kind of new frame correlation coefficient computing method, this method has defined a minimum difference principle, and it calculates the frame correlation coefficient that obtains and makes when the pre-filter amplitude-frequency response | H (f) | with the reference filter amplitude-frequency response | H
0(f) | difference minimum, that is to say this moment in overlapping bands, wave filter is at current imaging frame rate Fr
tWith benchmark imaging frame rate Fr
0Has close amplitude-frequency response characteristic down.Only guarantee when Flame Image Process that the wave filter amplitude-frequency response under the different imaging frame rates all approaches the benchmark amplitude-frequency response; Could guarantee in overlapping bands, to have close amplitude-frequency response characteristic between the wave filter amplitude-frequency response under the different imaging frame rates, thereby effectively guarantee the smoothing processing effect and the Noise Suppression effect of Flame Image Process.
As shown in Figure 2, frame correlation coefficient control method of the present invention need be set the frequency Fr of a benchmark
0Frame correlation coefficient a (Fr under (like real-time imaging frame rate, i.e. 50HZ or 60H)
0).This frame correlation coefficient a (Fr
0) can be according to aforementioned formula (3), confirm to calculate the relevant gear of frame of the corresponding ultrasonic image-forming system of the relevant control coefrficient k of frame, the frame relevant effect of its token image through selecting the relevant control coefrficient k of frame.With frame correlation coefficient a (Fr
0) corresponding amplitude-frequency response | H
0(f) | function as a reference, as current imaging frame rate Fr
tGreater than Fr
0The time, smoothing processing requires current imaging frame rate Fr
tFollowing frame correlation coefficient a (Fr
t) corresponding wave filter amplitude-frequency response | H (f) | at 0～Fr
0Has close response in/2 the scope; And, as current imaging frame rate Fr
tLess than Fr
0The time, require current imaging frame rate Fr
tFollowing frame correlation coefficient a (Fr
t) corresponding wave filter amplitude-frequency response | H (f) | at 0～Fr
tHas close response in/2 scopes.It is thus clear that, the overlapping bands of two wave filter amplitude-frequency responses can be expressed as [0, min (Fr
0/ 2, Fr
t/ 2)], through setting current imaging frame rate Fr
tUnder frame correlation coefficient a (Fr
t) make corresponding wave filter amplitude-frequency response frequency band [0, min (Fr
0/ 2, Fr
t/ 2) close], realize the smoothing processing of image, thereby guarantee to change under the situation greatly at imaging frame rate, also can obtain the frame relevant treatment effect of basically identical.
The control method of frame correlation coefficient shown in Figure 2, according to difference function E=F (| H (f) |, | H
0(f) |) the mode of setting up, can have multiple implementation, below be elaborated to several kinds of commonly used modes.
Embodiment 1: as shown in Figure 3, the E of difference described in the present embodiment is: the reference filter amplitude-frequency response | H
0(f) | with when the pre-filter amplitude-frequency response | H (f) | difference overlapping bands [0, min (Fr
0/ 2, Fr
t/ 2) quadratic sum], promptly
As shown in Figure 3, the detailed process of present embodiment is following:
Step 200, predefined benchmark imaging frame rate Fr
0, the benchmark imaging frame rate Fr here
0Be the real time imagery frame per second, like 50HZ or 60HZ.
Step 210 is with predefined benchmark imaging frame rate Fr
0, calculate its frame correlation coefficient a (Fr in the substitution aforementioned formula (3)
0), the result is as follows:
a(Fr
0)＝exp{-1/k×Fr
0}
Wherein, k is the relevant control coefrficient of frame, and it is confirmed that by experience this k value has determined the frame relevant effect of image, and corresponding with the relevant gear of frame of supersonic sounding instrument, and its span is 0～1.
Step 220 obtains frame correlation coefficient a (Fr with aforementioned calculation
0) be input in the transfer function model of digital filter the benchmark amplitude-frequency response of calculating filter | H
0(f) |, its corresponding frequency band range is 0～Fr
0/ 2.
Here, the transfer function model that belongs to single order IIR (InfiniteImpulse Response) digital filter Z territory with the transfer function model of digital filter is an example, specifies the benchmark amplitude-frequency response | H
0(f) | computation process.
According to aforementioned formula (2)
With z=exp{j2 π f ＇ }, f ＇=f/Fr, f ∈ [0, Fr/2] promptly obtains the frequency response H (f) of wave filter among the above-mentioned H of substitution (z), the amplitude-frequency response that its amplitude is corresponding | H (f) | expression formula is:
Wherein, f ＇ is the normalization numerical frequency, and f is an actual frequency, and Fr is an imaging frame rate.
According to above-mentioned formula, with frame correlation coefficient a (Fr
0), benchmark imaging frame rate Fr
0Be updated in the above-mentioned formula (5), can obtain about benchmark imaging frame rate Fr
0Under the benchmark amplitude-frequency response | H
0(f) | function expression.
Step 230 is set current imaging frame rate Fr
tUnder N candidate frame related coefficient a` (Fr
t), wherein N is the natural number greater than 2; And then respectively with this N candidate frame related coefficient a` (Fr
t) be input in the transfer function model (transfer function model in the single order infinite impulse response digital filter Z territory shown in above-mentioned formula (5)) of digital filter, calculate each candidate frame related coefficient a` (Fr
t) the corresponding current amplitude-frequency response of wave filter | H (f) |, its corresponding frequency band range is 0～Fr
t/ 2.
The N here candidate frame related coefficient a` (Fr
t) can rule of thumb come to confirm, also can come in such a way to confirm:
The first step is calculated current imaging frame rate Fr according to conventional method
tUnder the frame correlation coefficient estimated value
For example, with current imaging frame rate Fr
tIn the exponential function of frame correlation coefficient shown in the substitution aforementioned formula (3) and imaging frame rate, calculate its frame correlation coefficient estimated value
The result is as follows:
Can find out the candidate frame related coefficient a` (Fr of the corresponding different frame relevant effect of the relevant control coefrficient k of different frames from above-mentioned formula
t), calculating this estimated value
The time, its k value can be with reference to selected value in the step 210.
Second step is in said estimated value
Near N of selection value as candidate frame related coefficient a` (Fr
t), this N is the natural number more than or equal to 2.For example, with
Value is got 100 numerical value in the interval at center, as candidate frame related coefficient a` (Fr
t).
Step 240, at the benchmark amplitude-frequency response | H
0(f) | with current amplitude-frequency response | H (f) | overlapping bands in, more said respectively benchmark amplitude-frequency response | H
0(f) | with each candidate frame related coefficient a` (Fr
t) the current amplitude-frequency response of pairing wave filter | H (f) | between difference E, promptly according to above-mentioned formula (4) calculate [0, min (Fr
0/ 2, Fr
t/ 2) each a` (Fr]
t) pairing E value, according to the minimum principle of difference, from each candidate frame related coefficient a` (Fr
t) the middle current imaging frame rate Fr of conduct that selects an optimum
tUnder optimum frame related coefficient a (Fr
t), be about to the E value hour | H (f) | pairing candidate frame related coefficient a` (Fr
t) as current imaging frame rate Fr
tUnder optimum frame related coefficient a (Fr
t).
The present embodiment filter coefficient is according to when the pre-filter amplitude-frequency response | H (f) | frequency band [0, min (Fr
0/ 2, Fr
t/ 2)] interior and a (Fr
0) corresponding wave filter amplitude-frequency response | H
0(f) | have that minimum square error principle confirms.It is thus clear that the amplitude-frequency response function is the nonlinear function of frame correlation coefficient a and frequency f, current imaging frame rate Fr
tFollowing corresponding frame correlation coefficient a (Fr
t) can pass through numerical methods of solving.Utilize the amplitude-frequency response of the different imaging frame rates frame correlation coefficient correspondence down that the implementation of present embodiment obtains as shown in Figure 4.Visible by Fig. 4, the amplitude-frequency response in overlapping bands, under the different frame per second is also very close, has guaranteed that the fog-level of different frame per second situation undertissue motion is close, and noise also had suppresses effect preferably.
The account form of certain above-mentioned difference has more than a kind of, asks for an interview following embodiment.
Embodiment 2: the difference of present embodiment and embodiment 1 is that in step 240, the E of difference described in the present embodiment is: the reference filter amplitude-frequency response | H
0(f) | with when the pre-filter amplitude-frequency response | H (f) | difference overlapping bands [0, min (Fr
0/ 2, Fr
t/ 2) the absolute value sum], promptly
Through the size of this difference value E relatively, with this difference value E be hour when the pre-filter amplitude-frequency response | H (f) | pairing candidate frame related coefficient a` (Fr
t) confirm as current imaging frame rate Fr
tUnder optimum frame related coefficient a (Fr
t).All the other steps of present embodiment are identical with embodiment 1.The present embodiment filter coefficient is according to when the pre-filter amplitude-frequency response | H (f) | frequency band [0, min (Fr
0/ 2, Fr
t/ 2)] interior and a (Fr
0) corresponding wave filter amplitude-frequency response | H
0(f) | have that the principle of minimum Error Absolute Value sum confirms.With embodiment 1, current frame rate Fr
tFollowing corresponding frame correlation coefficient a (Fr
t) can find the solution through numerical method.Utilize amplitude-frequency response and Fig. 4 of the different imaging frame rates frame correlation coefficient correspondence down that the implementation of present embodiment obtains similar.
Embodiment 3: the difference of present embodiment and embodiment 1 is that in step 240, the E of difference described in the present embodiment is: overlapping bands [0, min (Fr
0/ 2, Fr
t/ 2)], at a setpoint frequency f
1Following reference filter amplitude-frequency response | H
0(f) | with when the pre-filter amplitude-frequency response | H (f) | the absolute value of difference, promptly
E＝‖H(f
1)|-|H
0(f
1)‖，(7)
Wherein, setpoint frequency f
1For above-mentioned overlapping bands (0, min (Fr
0/ 2, Fr
t/ 2) the arbitrary value] except that zero.Through setting a frequency f
1, calculate according to above-mentioned formula (7) then and the size of comparing difference value E, with this difference value E be hour work as the pre-filter amplitude-frequency response | H (f) | pairing candidate frame related coefficient a` (Fr
t) confirm as current imaging frame rate Fr
tUnder optimum frame related coefficient a (Fr
t).All the other steps of present embodiment are identical with embodiment 1.The present embodiment filter coefficient be according to frequency band [0, min (Fr
0/ 2, Fr
t/ 2) work as the pre-filter amplitude-frequency response] | H (f) | with a (Fr
0) corresponding wave filter amplitude-frequency response | H
0(f) | the minimum principle of absolute value of difference confirm.With embodiment 1, current frame rate Fr
tFollowing corresponding frame correlation coefficient a (Fr
t) can find the solution through numerical method.
Embodiment 4: the difference of present embodiment and embodiment 3 is that in step 240, difference described in the present embodiment is meant for minimum situation: at a setpoint frequency f
1Following reference filter amplitude-frequency response | H
0(f) | with when the pre-filter amplitude-frequency response | H (f) | difference be zero, promptly
E＝|H(f
1)|-|H
0(f
1)|＝0， (8)
Just, at setpoint frequency f
1On make | H
0(f) |=| H (f) |, this moment is when the pre-filter amplitude-frequency response | H (f) | pairing candidate frame related coefficient a` (Fr
t) be judged as current imaging frame rate Fr
tUnder optimum frame related coefficient a (Fr
t), f wherein
1=(0, min (Fr
0/ 2, Fr
t/ 2)].All the other steps of present embodiment are identical with embodiment 3, and these all the other steps are meant except that step 230 setting candidate frame related coefficient, reach the step relevant with the candidate frame related coefficient.
This instance is through setting a (Fr
t) make when the pre-filter amplitude-frequency response | H (f) | in certain predefined frequency f
1(f
1Be not equal to 0) on have and the reference filter amplitude-frequency response | H
0(f) | identical numerical value, even | H
0(f) |=| H (f) |, wherein, frequency f
1Optimum may be selected to be min (Fr
0/ 2, Fr
t/ 2), i.e. Fr
0/ 2 and Fr
t/ 2 smallers among both.And, through with a (Fr
0), f
1Substitution amplitude-frequency response function | H
0(f) |, and with f
1Substitution amplitude-frequency response function | H (f) |, find the solution above-mentioned equality (8) through parsing and promptly obtain current imaging frame rate Fr
tUnder optimum frame related coefficient a (Fr
t).
Present embodiment belongs to a kind of special circumstances among the embodiment 3; After the parsing of employing present embodiment is found the solution; Can simplify above-mentioned computation rule as shown in Figure 3; As can reduce step of setting the candidate frame related coefficient or the like, because above-mentioned formula (8) is a definite equality, it can directly calculate optimum frame related coefficient a (Fr
t).Present embodiment is for embodiment 1 to 3, and the implementation of present embodiment can be by imaging system software according to current imaging frame rate Fr because calculated amount is less
tThe relevant control coefrficient k of the frame that relevant gear with frame is corresponding calculates in real time and obtains current imaging frame rate Fr
tFrame correlation coefficient a (Fr
t).Utilize the amplitude-frequency response of the different imaging frame rates frame correlation coefficient correspondence down that the implementation of present embodiment obtains as shown in Figure 5.Visible by Fig. 5, amplitude-frequency response is also very close in overlapping bands, has guaranteed that the fog-level of different frame per second situation undertissue motion is close, and noise is also had the better inhibited effect.
Embodiment 5: in order to simplify calculation procedure, present embodiment is built the table step below also increasing on the basis of the foregoing description 1,2,3 or 4:
Confirmed current imaging frame rate Fr through embodiment 1,2,3 or 4
tUnder optimum frame related coefficient a (Fr
t) after, with this optimum frame related coefficient a (Fr
t) be stored in the relation table the current imaging frame rate Fr of this table record
t, current imaging frame rate Fr
tUnder optimum frame related coefficient a (Fr
t), and characterize the corresponding relation between the relevant control coefrficient k three of frame of frame relevant effect;
Then, when realtime graphic is handled, can be according to current imaging frame rate Fr
tAnd the relevant control coefrficient k of frame, search said relation table, confirm current imaging frame rate Fr in real time
tUnder optimum frame related coefficient a (Fr
t).The relation table about the optimum frame related coefficient that utilizes present embodiment to set up can be stored in the computer-readable recording medium.
Because the numerical methods of solving operand of the foregoing description 1 to 3 is bigger, so present embodiment is stored as a coefficient table through calculating the optimum frame related coefficient under the good different imaging frame rates in advance with these frame correlation coefficients.When needs are provided with different frame relevant effect, can be through change benchmark frame per second Fr
0Under the relevant control coefrficient k of frame, obtain different | H
0(f) |, numerical evaluation goes out corresponding different imaging frame rate Fr then
tUnder frame correlation coefficient a, thereby formed a two-dimentional coefficient table of different frame relevant effect and different imaging frame rates, will utilize the optimum frame related coefficient a (Fr under the current imaging frame rate that embodiment 1,2 or 3 solves then
t) add wherein, can obtain the above-mentioned relation table.Ultrasonic image-forming system is tabled look-up to the above-mentioned relation table according to the gear and the imaging frame rate of different frame relevant effect, can obtain the frame correlation coefficient under the current imaging frame rate.
The account form of common several species diversity has just been enumerated in above-mentioned analysis, and how the present invention calculates difference and unrestrictedly.Such as, say that according to what the difference mathematical model of above-mentioned each embodiment can be summarized the mathematical model of difference E can be expressed as the reference filter amplitude-frequency response | H
0(f) | with when the pre-filter amplitude-frequency response | H (f) | difference absolute value overlapping bands [0, min (Fr
0/ 2, Fr
t/ 2) the n power sum], promptly
Wherein, n is the natural number greater than zero.In the above-mentioned formula (9) during n=1, i.e. situation shown in the foregoing description 2; When n=2, i.e. situation shown in the foregoing description 1; When n=3, also can adopt, just calculate relative complex some.In addition, according to mathematical model and the formula (9) of above-mentioned each difference E, also have a kind of situation of equivalence in fact, promptly the mathematical model of difference E can be expressed as the reference filter amplitude-frequency response | H
0(f) | with when the pre-filter amplitude-frequency response | H (f) | difference overlapping bands [0, min (Fr
0/ 2, Fr
t/ 2) the 2n power sum], promptly
Wherein, n is the natural number greater than zero.In the above-mentioned formula (10) during n=1, i.e. situation shown in the foregoing description 1; When n=2, equally also be desirable account form.
Therefore; The control method of the frame correlation coefficient that the present invention proposes is not limited to organizes the gray scale imaging processing; Can also be applicable to all kinds of imaging patterns such as color flow angiography, tissue doppler imaging; It can confirm the optimum frame related coefficient accurately in order in the frame relevant treatment, guarantees the effect of picture smooth treatment.
Usually effectively image information changes between consecutive frame not quite; Invalid noise information then changes bigger usually; Utilize the relevant low-pass filter characteristic that is adopted of frame; Can better keep interframe to change little low frequency AP information, and suppress the interframe variation invalid noise of high frequency greatly.Adopt in the frame correlation filter amplitude-frequency response shown in Figure 1 of prior art acquisition, can find out significantly between each curve to differ greatly, have more noise information; And adopt in Fig. 4 that method of the present invention obtains and the frame correlation filter amplitude-frequency response shown in Figure 5, have consistent frame relevant effect between each adjacent curve basically, and each curve has similar amplitude versus frequency characte in overlapped frequency band.It is thus clear that the present invention passes through different imaging frame rates, thereby the corresponding frame correlation coefficient of optimized choice obtains the frame relevant effect of the basically identical under the different imaging frame rates, has guaranteed the smoothing processing effect and the Noise Suppression effect of Flame Image Process effectively.
Control method according to above-mentioned frame correlation coefficient; The present invention is also corresponding to provide a kind of frame correlation coefficient controller; The optimum frame related coefficient that it can directly obtain under the current imaging frame rate guarantees that the wave filter amplitude-frequency response under the different imaging frame rates has close frequency response characteristic in overlapping bands.As shown in Figure 6, said controller comprises:
The 3rd unit 320 is used for more current imaging frame rate Fr
tUnder when the pre-filter amplitude-frequency response | H (f) | with said reference filter amplitude-frequency response | H
0(f) | the difference E in overlapping bands, export said difference E for hour when the pre-filter amplitude-frequency response | H (f) | pairing frame correlation coefficient a (Fr
t).
The parameter a (Fr here
0), Fr
0, | H
0(f) |, Fr
t, E, | H (f) |, a (Fr
t), the calculating or the setting principle of overlapping bands, can be referring to the implementation process of the control method of above-mentioned frame correlation coefficient.
As shown in Figure 7, for calculated difference E, can in the 3rd unit 320 of controller, set up with lower unit:
Comparing unit 322 is used for more said difference E, export said difference E for hour when the pre-filter amplitude-frequency response | H (f) | pairing frame correlation coefficient.Adopt this structure, and utilize the foregoing description 4, promptly can realize calculating in real time the optimum frame related coefficient under the current imaging frame rate; If but adopt the foregoing description 1,2 or at 3 o'clock; Relate in the time of from a plurality of candidate frame related coefficients, to select; Said the 3rd unit 320 just need increase following unit again on the basis of differential pattern unit 321 and comparing unit 322; As shown in Figure 7, said the 3rd unit 320 also comprises:
Candidate set computing module 323 is used to generate current imaging frame rate Fr
tUnder candidate frame related coefficient a` (Fr
t);
Wave filter amplitude-frequency response computing module 324 is used for said candidate frame related coefficient a` (Fr
t) be input in the transfer function model 325 of digital filter calculated candidate frame correlation coefficient a` (Fr
t) corresponding to the pre-filter amplitude-frequency response | H (f) |, with this when the pre-filter amplitude-frequency response | H (f) | be input to and carry out difference in the said differential pattern unit 321 and calculate, and relatively obtain corresponding frame correlation coefficient in the back through comparing unit 322.
Need to preserve and set up quick demand with relation table based on data result, frame correlation coefficient controller of the present invention can increase by a storage unit 330 on the basis of structure shown in Figure 6, and it is used to write down current imaging frame rate Fr
tUnder optimum frame related coefficient a (Fr
t), this optimum frame related coefficient a (Fr
t) be said difference E for hour when the pre-filter amplitude-frequency response | H (f) | pairing frame correlation coefficient, promptly from selecting frame correlation coefficient a` (Fr
t) in the frame correlation coefficient selected according to the minimum principle of difference.The storage unit 330 here can be provided with a relation table, and with the optimum frame related coefficient a (Fr that writes down
t) add wherein, utilize this relation table to write down current imaging frame rate Fr
t, current imaging frame rate Fr
tUnder optimum frame related coefficient a (Fr
t), and characterize the corresponding relation between the relevant control coefrficient k three of frame of frame relevant effect.
Can find out from the structrual description of above-mentioned frame correlation coefficient controller; It mainly is to utilize the control method of above-mentioned frame correlation coefficient and set up a kind of product; This product adopts the control method of above-mentioned frame correlation coefficient can obtain the optimum frame related coefficient under the current frame rate; Utilizing this frame correlation coefficient to carry out ultrasonoscopy handles; Can guarantee that the wave filter amplitude-frequency response under the different imaging frame rates has close frequency response characteristic in overlapping bands, and effectively guarantee the smooth effect and the noise suppression effect of Flame Image Process.
Based on the application of control method in the frame relevant treatment of above-mentioned frame correlation coefficient, the present invention also provides a kind of ultrasonoscopy frame correlation processing method, and this method may further comprise the steps:
Step S1 adopts the control method of above-mentioned frame correlation coefficient to obtain current imaging frame rate Fr
tPairing optimum frame related coefficient a (Fr
t);
Step S2 is with the optimum frame related coefficient a (Fr that obtains
t) be input in the mathematical model of frame relevant treatment, carry out the frame related operation, the current frame image after obtaining to handle.The frame relevant treatment mathematical model here can adopt aforementioned formula (1), and its specific embodiment can be referring to Fig. 3.Visible by Fig. 3, calculate optimum frame related coefficient a (Fr at the foregoing description 1,2 or 3
t) after, as long as utilize this frame correlation coefficient to carry out the frame relevant treatment of image.
In the ultrasonoscopy frame correlation processing method proposed by the invention, about optimum frame related coefficient a (Fr
t) computation process with and the technical scheme expanded referring to above-mentioned explanation, in this no longer repeat specification about the frame correlation coefficient control method.
On the basis of above-mentioned frame correlation coefficient control method and frame correlation processing method thereof; The present invention also provides a kind of new ultrasonic image-forming system; As shown in Figure 8, said ultrasonic image-forming system comprises: ultrasonoscopy processing module and the probe that is used to obtain ultrasonoscopy; Be provided with a frame correlation coefficient controller 400 and a frame relevant treatment unit 410 in the said ultrasonoscopy processing module, the frame correlation coefficient controller 400 here adopts structure as shown in Figure 6, specifically can be referring to above stated specification; Said frame relevant treatment unit 410 is used for the frame correlation coefficient of Unit the 3rd 320 outputs of received frame related coefficient controller 400; And this frame correlation coefficient is input in the mathematical model of frame relevant treatment and carries out the frame related operation, obtain the current frame image after the smoothing processing.In this ultrasonic image-forming system, about the content of frame correlation coefficient controller 400 with and the technical scheme expanded referring to above-mentioned explanation, in this no longer repeat specification about frame correlation coefficient controller shown in Figure 6.
Based on the implementation that above-mentioned optimum frame related coefficient can adopt lookup table mode to carry out, the present invention also provides a kind of new ultrasonic image-forming system, and is as shown in Figure 9, comprises in this system: ultrasonoscopy processing module and the probe that is used to obtain ultrasonoscopy; And; In the ultrasonoscopy processing module, be provided with a storage medium 430, search a unit 420 and a frame relevant treatment unit 410 shown in Figure 6; Said storage medium 430 is preserved a relation table; This table record has optimum frame related coefficient under current imaging frame rate, the current imaging frame rate, and characterizes the corresponding relation between the relevant control coefrficient three of frame of frame relevant effect, the difference minimum of reference filter amplitude-frequency response in overlapping bands under pre-filter amplitude-frequency response and a benchmark imaging frame rate that said optimum frame related coefficient is corresponding; Said frame relevant treatment unit 410 links to each other with said storage medium 430 through the said unit 420 of searching; The said unit 420 of searching is used for searching said relation table and current imaging frame rate and the corresponding optimum frame related coefficient of the relevant control coefrficient of present frame, and this optimum frame related coefficient is sent in the said frame relevant treatment unit 410; Said frame relevant treatment unit 410 is used for the mathematical model that said optimum frame related coefficient is input to the frame relevant treatment is carried out the frame related operation, the current frame image after obtaining to handle.
Ultrasonic image-forming system shown in Figure 9 is: the optimum frame related coefficient that adopts above-mentioned frame correlation coefficient control method to be obtained is set up a relation table; The corresponding frame correlation coefficient of acquisition is used for when Flame Image Process, tabling look-up in real time; Do like this and can reduce calculated amount, improve the speed of Flame Image Process.How the key issue of building table is the calculating optimum frame correlation coefficient, and it can be referring to above-mentioned specifying about the frame correlation coefficient control method.Such as: adopt the mode of the foregoing description 1,2 or 3 like figure, the optimum frame related coefficient can adopt following method to obtain:
At first, set a plurality of candidate frame related coefficients under the current imaging frame rate, the system of selection of this candidate frame related coefficient is following:
According to the exponential function of frame correlation coefficient and imaging frame rate, calculate the frame correlation coefficient estimated value under the current imaging frame rate; Near said estimated value, select N value as the candidate frame related coefficient, this N is the natural number more than or equal to 2;
Then, respectively said N candidate frame related coefficient is input in the transfer function model of digital filter, calculate each candidate frame related coefficient corresponding when the pre-filter amplitude-frequency response;
At last; Respectively with corresponding the comparing when pre-filter amplitude-frequency response and said reference filter amplitude-frequency response (its computation process can referring to the related description of the foregoing description 1) of said N candidate frame related coefficient; According to working as pre-filter amplitude-frequency response and the difference of said reference filter amplitude-frequency response in overlapping bands is minimum principle; From said N candidate frame related coefficient, select one as the pairing optimum frame related coefficient of current imaging frame rate, promptly said difference be hour when the pairing candidate frame related coefficient of pre-filter amplitude-frequency response.If implement 4 or 5 but adopt, then can specify referring to above-mentioned.
In sum; The present invention is through setting the frame correlation coefficient under the different frame per second; Feasible corresponding wave filter amplitude-frequency response has close frequency response characteristic in overlapping frequency band requirement; A kind of computing method of frame correlation coefficient are provided; Based on these computing method, the present invention has derived a kind of frame CCU that can process stand-alone product and be used to obtain the optimum frame related coefficient again, has a kind ofly utilized the frame correlation processing method of above-mentioned frame correlation coefficient computing method and utilize the ultrasonic image-forming system of this method and the ultrasonic image-forming system that a kind of employing contains the storage medium of optimum frame related coefficient relation table.For last ultrasonic image-forming system, the optimum frame related coefficient relation table that contains in its storage medium adopts above-mentioned frame correlation coefficient computing method to obtain.Disclosed the whole bag of tricks of the present invention and system; All can calculate the optimum frame related coefficient that obtains under the current imaging frame rate according to the minimum difference principle; And be used for subsequent image processing; Make the present invention when imaging frame rate takes place than about-face, the treatment effect that ultrasonic image frame is relevant remains unchanged basically.
Should be understood that, concerning those of ordinary skills, can improve or conversion, and all these improvement and conversion all should belong to the protection domain of accompanying claims of the present invention according to above-mentioned explanation.
Claims (24)
1. the control method of a frame correlation coefficient is characterized in that, said method comprises:
The reference filter amplitude-frequency response under a pre-filter amplitude-frequency response and a benchmark imaging frame rate under A, the more current imaging frame rate obtains said when pre-filter amplitude-frequency response and the difference of reference filter amplitude-frequency response in overlapping bands;
B, with said difference be hour when the pairing frame correlation coefficient of pre-filter amplitude-frequency response, be judged to be the pairing optimum frame related coefficient of current imaging frame rate.
2. method according to claim 1 is characterized in that, the scope of said overlapping bands is zero closed interval to the T value, and this T value does
Fr 0/ 2 draws
Fr t/ 2 smallers among both, wherein
Fr 0Be the benchmark imaging frame rate,
Fr tBe current imaging frame rate.
3. method according to claim 1 is characterized in that, said steps A may further comprise the steps:
A plurality of candidate frame related coefficients under A1, the current imaging frame rate of setting;
A2, according to the transport function of digital filter, calculate each candidate frame related coefficient corresponding when the pre-filter amplitude-frequency response;
A3, respectively with corresponding the comparing when pre-filter amplitude-frequency response and said reference filter amplitude-frequency response of said each candidate frame related coefficient, obtain each and work as pre-filter amplitude-frequency response and the difference of reference filter amplitude-frequency response in overlapping bands;
And, among the said step B,, obtain minimum said difference through comparing the pairing said difference of each candidate frame related coefficient.
4. according to claim 1,2 or 3 described methods, it is characterized in that said difference is: the absolute value of said reference filter amplitude-frequency response and said difference when the pre-filter amplitude-frequency response is in overlapping bands
nThe power sum, wherein,
nFor greater than zero natural number.
5. method according to claim 1 is characterized in that, said difference is: said reference filter amplitude-frequency response under a setpoint frequency and said poor when the pre-filter amplitude-frequency response.
6. method according to claim 2 is characterized in that, said difference is: said reference filter amplitude-frequency response under a setpoint frequency and said poor when the pre-filter amplitude-frequency response.
7. method according to claim 3 is characterized in that, said difference is: said reference filter amplitude-frequency response under a setpoint frequency and said poor when the pre-filter amplitude-frequency response.
8. according to claim 5,6 or 7 described methods, it is characterized in that said setpoint frequency is the arbitrary value except that zero in the said overlapping bands.
9. method according to claim 6 is characterized in that, said difference is meant for minimum situation: said reference filter amplitude-frequency response under a setpoint frequency and the said difference of working as the pre-filter amplitude-frequency response are zero.
10. method according to claim 9 is characterized in that said setpoint frequency does
Fr 0/ 2 draws
Fr t/ 2 smallers among both.
11. method according to claim 3 is characterized in that, the system of selection of said candidate frame related coefficient may further comprise the steps:
According to the exponential function of frame correlation coefficient and imaging frame rate, calculate the frame correlation coefficient estimated value under the current imaging frame rate;
Near said estimated value, select N value as the candidate frame related coefficient, this N is the natural number more than or equal to 2.
12., it is characterized in that said method is further comprising the steps of according to claim 1,2 or 3 described methods:
C, said optimum frame related coefficient is stored in the relation table, the optimum frame related coefficient under the current imaging frame rate of this table record, the current imaging frame rate, and characterizes the corresponding relation between the relevant control coefrficient three of frame of frame relevant effect;
D, according to the relevant control coefrficient of current imaging frame rate and frame, search said relation table, confirm the optimum frame related coefficient under the current imaging frame rate.
13. a frame correlation coefficient controller is characterized in that, said controller comprises:
First module is used for calculating the frame correlation coefficient under the benchmark imaging frame rate according to preset benchmark imaging frame rate;
Unit second, the frame correlation coefficient that is used for obtaining based on first module calculates the reference filter amplitude-frequency response; And
Unit the 3rd, be used under the more current imaging frame rate when pre-filter amplitude-frequency response and the difference of said reference filter amplitude-frequency response in overlapping bands, export said difference for hour work as the pairing frame correlation coefficient of pre-filter amplitude-frequency response.
14. controller according to claim 13 is characterized in that, said Unit the 3rd comprises:
The differential pattern unit is used to calculate said difference when pre-filter amplitude-frequency response and said reference filter amplitude-frequency response; And
Comparing unit is used for more said difference, export said difference for hour when the pairing frame correlation coefficient of pre-filter amplitude-frequency response.
15. controller according to claim 14 is characterized in that, said Unit the 3rd also comprises:
The candidate set computing module is used to generate the candidate frame related coefficient under the current imaging frame rate; And
Wave filter amplitude-frequency response computing module; Be used for said candidate frame related coefficient is input to the transport function of digital filter; The calculated candidate frame correlation coefficient corresponding when the pre-filter amplitude-frequency response, this is input in the said differential pattern unit when the pre-filter amplitude-frequency response.
16., it is characterized in that said controller also comprises according to claim 13 or 14 described controllers:
Be provided with the storage unit of a relation table, said relation table is used to write down optimum frame related coefficient under current imaging frame rate, the current imaging frame rate, and characterizes the corresponding relation between the relevant control coefrficient three of frame of frame relevant effect; This optimum frame related coefficient be said difference for hour when the pairing frame correlation coefficient of pre-filter amplitude-frequency response.
17. a ultrasonoscopy frame correlation processing method is characterized in that, said method comprising the steps of:
The reference filter amplitude-frequency response under a pre-filter amplitude-frequency response and a benchmark imaging frame rate under S1, the more current imaging frame rate obtains said when pre-filter amplitude-frequency response and the difference of reference filter amplitude-frequency response in overlapping bands;
S2 with said difference be hour when the pairing frame correlation coefficient of pre-filter amplitude-frequency response, be judged to be the pairing optimum frame related coefficient of current imaging frame rate;
The optimum frame related coefficient that S3, utilization obtain is carried out the frame related operation, the current frame image after obtaining to handle.
18. method according to claim 17 is characterized in that, may further comprise the steps between said step S2 and the step S3:
Said optimum frame related coefficient is stored in the relation table, the optimum frame related coefficient under the current imaging frame rate of this table record, the current imaging frame rate, and characterizes the corresponding relation between the relevant control coefrficient three of frame of frame relevant effect;
According to current imaging frame rate and the relevant control coefrficient of frame, search said relation table, confirm the optimum frame related coefficient under the current imaging frame rate.
19. a ultrasonic image-forming system, said system comprises: ultrasonoscopy processing module and the probe that is used to obtain ultrasonoscopy; It is characterized in that be provided with a frame correlation coefficient controller and a frame relevant treatment unit in the said ultrasonoscopy processing module, this controller comprises:
First module is used for calculating the frame correlation coefficient under the benchmark imaging frame rate according to preset benchmark imaging frame rate;
Unit second, the frame correlation coefficient that is used for obtaining based on first module calculates the reference filter amplitude-frequency response; And
Unit the 3rd, be used under the more current imaging frame rate when pre-filter amplitude-frequency response and the difference of said reference filter amplitude-frequency response in overlapping bands, export said difference for hour work as the pairing frame correlation coefficient of pre-filter amplitude-frequency response;
Said frame relevant treatment unit is used to receive the frame correlation coefficient of said Unit the 3rd output, and utilizes this frame correlation coefficient to carry out the frame related operation, the current frame image after obtaining to handle.
20. ultrasonic image-forming system according to claim 19 is characterized in that, said Unit the 3rd comprises:
The differential pattern unit is used to calculate said difference when pre-filter amplitude-frequency response and said reference filter amplitude-frequency response; And
Comparing unit is used for more said difference, export said difference for hour when the pairing candidate frame related coefficient of pre-filter amplitude-frequency response.
21. ultrasonic image-forming system according to claim 20 is characterized in that, said Unit the 3rd also comprises:
The candidate set computing module is used to generate the candidate frame related coefficient under the current imaging frame rate; And
Wave filter amplitude-frequency response computing module; Be used for said candidate frame related coefficient is input to the transport function of digital filter; The calculated candidate frame correlation coefficient corresponding when the pre-filter amplitude-frequency response, this is input in the said differential pattern unit when the pre-filter amplitude-frequency response.
22., it is characterized in that said ultrasonic image-forming system also comprises according to claim 19 or 20 described ultrasonic image-forming systems:
Be provided with the storage unit of a relation table, said relation table is used to write down optimum frame related coefficient under current imaging frame rate, the current imaging frame rate, and characterizes the corresponding relation between the relevant control coefrficient three of frame of frame relevant effect; This optimum frame related coefficient be said difference for hour when the pairing frame correlation coefficient of pre-filter amplitude-frequency response.
23. a ultrasonic image-forming system, said system comprises: ultrasonoscopy processing module and the probe that is used to obtain ultrasonoscopy; It is characterized in that, be provided with a storage medium in the said ultrasonoscopy processing module, search a unit and a frame relevant treatment unit,
Said storage medium is preserved a relation table; This table record has optimum frame related coefficient under current imaging frame rate, the current imaging frame rate, and characterizes the corresponding relation between the relevant control coefrficient three of frame of frame relevant effect, the difference minimum of reference filter amplitude-frequency response in overlapping bands under pre-filter amplitude-frequency response and a benchmark imaging frame rate that said optimum frame related coefficient is corresponding;
Said frame relevant treatment unit links to each other with said storage medium through the said unit of searching; The said unit of searching is used for searching said relation table and current imaging frame rate and the corresponding optimum frame related coefficient of the relevant control coefrficient of present frame, and this optimum frame related coefficient is sent in the said frame relevant treatment unit;
Said frame relevant treatment unit is used to utilize said optimum frame related coefficient to carry out the frame related operation, the current frame image after obtaining to handle.
24. ultrasonic image-forming system according to claim 23 is characterized in that, the computation process of said optimum frame related coefficient may further comprise the steps:
A plurality of candidate frame related coefficients under C1, the current imaging frame rate of setting, the system of selection of this candidate frame related coefficient is following:
According to the exponential function of frame correlation coefficient and imaging frame rate, calculate the frame correlation coefficient estimated value under the current imaging frame rate; Near said estimated value, select N value as the candidate frame related coefficient, this N is the natural number more than or equal to 2;
C2, according to the transport function of digital filter, calculate each candidate frame related coefficient corresponding when the pre-filter amplitude-frequency response;
C3, respectively with corresponding the comparing when pre-filter amplitude-frequency response and said reference filter amplitude-frequency response of said N candidate frame related coefficient; With said difference be hour when the pairing candidate frame related coefficient of pre-filter amplitude-frequency response, be judged to be said optimum frame related coefficient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200810142787XA CN101645167B (en) | 2008-08-06 | 2008-08-06 | Method and device for frame correlation coefficient control and frame correlation processing, and ultrasonic imaging system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN200810142787XA CN101645167B (en) | 2008-08-06 | 2008-08-06 | Method and device for frame correlation coefficient control and frame correlation processing, and ultrasonic imaging system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101645167A CN101645167A (en) | 2010-02-10 |
| CN101645167B true CN101645167B (en) | 2012-08-15 |
Family
ID=41657045
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN200810142787XA Active CN101645167B (en) | 2008-08-06 | 2008-08-06 | Method and device for frame correlation coefficient control and frame correlation processing, and ultrasonic imaging system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101645167B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102908163B (en) * | 2011-08-02 | 2016-01-20 | 深圳迈瑞生物医疗电子股份有限公司 | A kind of frame correlation technique and device thereof and ultrasonic image-forming system |
| CN103929143A (en) * | 2014-04-11 | 2014-07-16 | 深圳市东微智能科技有限公司 | Method for acquiring amplitude-frequency response characteristic of filter |
| CN106603885B (en) * | 2015-10-20 | 2019-07-19 | 腾讯科技（深圳）有限公司 | Method of video image processing and device |
| WO2019075621A1 (en) * | 2017-10-16 | 2019-04-25 | 北京深迈瑞医疗电子技术研究院有限公司 | Ultrasonic imaging device and system, and image enhancement method for ultrasonic radiography and imaging |
| KR102527296B1 (en) * | 2018-05-04 | 2023-05-02 | 삼성디스플레이 주식회사 | Display system and method of synchronizing a frame driving timing for the same |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5060515A (en) * | 1989-03-01 | 1991-10-29 | Kabushiki Kaisha Toshiba | Image signal processing circuit for ultrasonic imaging apparatus |
| US6056691A (en) * | 1998-06-24 | 2000-05-02 | Ecton, Inc. | System for collecting ultrasound imaging data at an adjustable collection image frame rate |
| CN1541622A (en) * | 2003-03-27 | 2004-11-03 | Ge医药系统环球科技公司 | Ultrasonic imaging method and ultrasonic diagnostic apparatus |
| CN1909836A (en) * | 2004-01-08 | 2007-02-07 | 松下电器产业株式会社 | Ultrasonic diagnosis apparatus |
-
2008
- 2008-08-06 CN CN200810142787XA patent/CN101645167B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5060515A (en) * | 1989-03-01 | 1991-10-29 | Kabushiki Kaisha Toshiba | Image signal processing circuit for ultrasonic imaging apparatus |
| US6056691A (en) * | 1998-06-24 | 2000-05-02 | Ecton, Inc. | System for collecting ultrasound imaging data at an adjustable collection image frame rate |
| CN1541622A (en) * | 2003-03-27 | 2004-11-03 | Ge医药系统环球科技公司 | Ultrasonic imaging method and ultrasonic diagnostic apparatus |
| CN1909836A (en) * | 2004-01-08 | 2007-02-07 | 松下电器产业株式会社 | Ultrasonic diagnosis apparatus |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101645167A (en) | 2010-02-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101527034B (en) | Method and device for correlating adaptive frames | |
| CN101645167B (en) | Method and device for frame correlation coefficient control and frame correlation processing, and ultrasonic imaging system | |
| CN100544677C (en) | Handle the method for Doppler signal gap | |
| JP4722283B2 (en) | Method and apparatus for motion visualization in ultrasonic flow imaging using continuous data acquisition | |
| US20080021945A1 (en) | Method of processing spatial-temporal data processing | |
| US5642732A (en) | Apparatus and method for estimating missing doppler signals and spectra | |
| US5097836A (en) | Untrasound diagnostic equipment for calculating and displaying integrated backscatter or scattering coefficients by using scattering power or scattering power spectrum of blood | |
| JP4627366B2 (en) | Method and apparatus for motion visualization in ultrasonic flow imaging using packet data acquisition | |
| US20100185093A1 (en) | System and method for processing a real-time ultrasound signal within a time window | |
| JPS6034434A (en) | Vingmed as | |
| Flynn et al. | Estimation and display for vector Doppler imaging using planewave transmissions | |
| Mozumi et al. | Utilization of singular value decomposition in high-frame-rate cardiac blood flow imaging | |
| US6306093B1 (en) | Method and system for ultrasound enhanced-resolution spectral Doppler | |
| US8162838B2 (en) | Doppler signal spectrum calculating method and apparatus | |
| CN105708496A (en) | Blood flow information multi-dimensional imaging system based on ultrasound | |
| CN107122724A (en) | A kind of method of the online denoising of sensing data based on adaptive-filtering | |
| JPH11155860A (en) | Ultrasonic color doppler video system and method for filtering ultrasonic doppler signal | |
| JPH0616782B2 (en) | Method and apparatus for combining continuous estimated signals in real time | |
| US4691707A (en) | AM/FM hybrid ultrasonic imaging | |
| Ukte et al. | Two empirical methods for improving the performance of statistical multirate high-resolution signal reconstruction | |
| CN101807298A (en) | Method for determining intensity of speckle noise in images | |
| JP2801734B2 (en) | Ultrasound diagnostic equipment | |
| Yamakoshi et al. | Novel estimation method for shear wave displacement amplitude excited by vibrator | |
| Valeau et al. | Instantaneous frequency tracking of a sine wave phase modulation signal | |
| CN105212963A (en) | Realize the wall filtering method and system of ultrasonic color doppler blood flow imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20100210 Assignee: Shenzhen Mindray Animal Medical Technology Co.,Ltd. Assignor: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd. Contract record no.: X2022440020009 Denomination of invention: Frame correlation coefficient control and frame correlation processing method, device and ultrasonic imaging system Granted publication date: 20120815 License type: Common License Record date: 20220804 |