CN103829975B

CN103829975B - Method for calculating ultrasonic elastography displacement value

Info

Publication number: CN103829975B
Application number: CN201410096941.XA
Authority: CN
Inventors: 周兰; 周一兰; 凌涛; 陈惠人
Original assignee: Vinno Technology Suzhou Co Ltd
Current assignee: Feiyinuo Technology Co ltd
Priority date: 2014-03-17
Filing date: 2014-03-17
Publication date: 2017-01-18
Anticipated expiration: 2034-03-17
Also published as: CN103829975A

Abstract

The invention provides a method for calculating an ultrasonic elastography displacement value. The method comprises the first step of acquiring echo signals of a front frame and echo signals of a back frame of a compressed organization, wherein the echo signals in each frame comprise a plurality of data lines; the second step of selecting a plurality of displacement estimation points on each data line and dividing the echo signals into a plurality of parts in the axial direction of the data lines; the third step of selecting at least one part, calculating the absolute displacement value of each displacement estimation point in the part, working out the priori values of the data lines according to the absolute displacement values of the data lines, and serving the priori value of each data line as the calculation basis for calculating the absolute displacement values of one or more parts of the same data line or the adjacent data lines. According to the method for calculating the ultrasonic elastography displacement value, the hunting zone can be adjusted in a self-adaptation mode according to the different stresses, the calculated amount is greatly reduced, very high robustness is achieved so that error extension can be avoided, and the real-time performance and accuracy of elastography are guaranteed.

Description

The computational methods of ultrasonic elastograph imaging shift value

Technical field

The present invention relates to a kind of computational methods of ultrasonic elastograph imaging shift value, it is mainly used in medical ultrasonic diagnostic technology Field.

Background technology

Ultrasonic elastograph imaging is the new model of medical ultrasound image, and it mainly utilizes the ultrasonic mechanics parameters to tissue Checked and be imaged.

When tissue occurs pathology, its mechanics parameters generally also changes.For example, when breast cancer, prostatitis When the diseases such as gland cancer occur, its corresponding coefficient of elasticity (Young's modulus shows as hardness) also can change.Elastogram is The inspection of clinically numerous disease, especially tumour and diagnosis provide valuable help.

Elastogram generally first passes through and tissue is applied with excitation (entering row energization using probe compressing or using acoustic pressure), so Ultrasonic signal before and after afterwards tissue being deformed upon is analyzed, and calculates the Displacements Distribution of tissue, then recycles gradient fortune Calculate and calculate Strain Distribution, finally Strain Distribution is shown with pseudo- color form.

Wherein, the Displacements Distribution of computation organization is the core link during elastogram, the quality directly shadow of this algorithm The quality of imaging and the speed of imaging are rung.Usually, as shown in figure 1, Fig. 1 is ultrasonic elastograph imaging shift value in prior art Computational methods flow chart, accordingly, the usual step of elastogram is:

Acoustic pressure power of s1, use probe compressing or ultrasonic generation etc., along a direction urging tissue, obtains group with this Knit two frame of digital signals in pressurized process.

S2, on two frame of digital signals segmentation choose Displacement Estimation point.

On s3, a bit of signal that the Displacement Estimation point on wherein one frame is represented and another frame, the signal of correspondence position does Related operation simultaneously searches for the maximum of coefficient correlation, can obtain the exhausted of this Displacement Estimation point according to coefficient correlation maximum value position To shift value.

S4, the absolute displacement values to acquisition in axial direction carry out gradient algorithm, obtain strain signal.

S5, Pseudo-color technology is carried out to strain signal, as the output of elastogram result.

Wherein, obtain Displacement Estimation point shift value be most important in elastogram be also step the most time-consuming, not only Because its calculating itself is time-consuming more, and due to searching for a segment signal of the first frame Displacement Estimation point representative in the second frame During position, neither one hunting zone, need using all signals in a very wide range in the second frame signal and this displacement The signal of estimation point does matching primitives, considerably increases the calculating time, becomes an important obstruction of real-time elastogram.

Content of the invention

In order to solve the above problems, it is an object of the invention to provide a kind of calculating side of ultrasonic elastograph imaging shift value Method, the priori value searching algorithm that the method is limited by zone boundary calculates the absolute displacement values of Displacement Estimation point, decreases super The calculating time of acoustic elasticity imaging.

Correspondingly, the computational methods of the ultrasonic elastograph imaging shift value of an embodiment of the present invention, methods described include with Lower step:

Two frame echo-signals before and after s1, acquisition compression organization, every frame echo-signal includes a plurality of data lines；

S2, multiple Displacement Estimation points are chosen on every data line, and axially described echo-signal being drawn along data wire It is divided into some；

S3, choose at least a portion, calculate the absolute displacement values in this partly each Displacement Estimation point interior, according to every number Calculate the priori value of every data line according to the absolute displacement values on line, and the priori value of every data line is same as calculating The calculating basis of the absolute displacement values of certain a part of or a few part on data line；

Wherein, some that described echo-signal divides includes Part I and Part II, described s3 step tool Body includes:

Choose Part I, calculate the absolute displacement values of each Displacement Estimation point in described Part I, and according to described The absolute displacement values of each the Displacement Estimation point in Part I calculate the Part I priori value of every data line；

According to Part I and Part II position relationship in the axial direction, by the Part I priori value of every data line As an end value of the region of search of Part II intrinsic displacement estimation point on same data line, to be calculated by searching algorithm The absolute displacement values of each Displacement Estimation point in Part II on same data line.

As a further improvement on the present invention, it is separated with least one portion between described Part I and described Part II Point.

As a further improvement on the present invention, some that described echo-signal divides also includes Part III, described Part III is presented axially between described Part I and Part II；Described s3 step includes:

Calculate the second of every data line according to the absolute displacement values of each the Displacement Estimation point in described Part II Part priori value；

According to Part I, Part II, Part III position relationship in the axial direction, by first of every data line Priori value and Part II priori value is divided to form the region of search of Part III intrinsic displacement estimation point on same data line, with logical Cross the absolute displacement values that searching algorithm calculates each Displacement Estimation point in Part III on same data line.

As a further improvement on the present invention, " calculate every data line according to the absolute displacement values on every data line Priori value " step specifically includes:

Filter the noise in the absolute displacement values having calculated that；

Take the intermediate value of the absolute displacement values after filtering noise or round the most absolute displacement values of rear occurrence number as elder generation Test value.

Correspondingly, the computational methods of the ultrasonic elastograph imaging shift value of another embodiment of the present invention, methods described includes Following steps:

S2, multiple Displacement Estimation points are chosen on every data line；

S3, by searching algorithm calculate at least one data wire on each Displacement Estimation point absolute displacement values；

S4, calculate priori value according to the absolute displacement values of each Displacement Estimation point, and phase is calculated according to described priori value The absolute displacement values of each the Displacement Estimation point on adjacent data wire；

Wherein, " according to the absolute displacement values of each the Displacement Estimation point on described priori value calculating adjacent data line " step Specifically include:

Described echo-signal is divided into some by the axial direction along data wire, by corresponding to adjacent data line juxtaposed position Part and its priori value formed the region of search, so that each Displacement Estimation point on adjacent data line is calculated by searching algorithm Absolute displacement values.

As a further improvement on the present invention, " calculate priori value according to the absolute displacement values of each Displacement Estimation point " Step specifically includes:

Filter the noise in the absolute displacement values having calculated that；

Correspondingly, the computational methods of the ultrasonic elastograph imaging shift value of a further embodiment of the present invention, methods described includes Following steps:

The method comprising the steps of:

S3, the absolute displacement of each the Displacement Estimation point being calculated by searching algorithm in certain part in a certain data line Value；

S4, calculate part priori value according to the absolute displacement values of each Displacement Estimation point, and according to described part priori Value calculates the absolute displacement values of each the Displacement Estimation point in another part of this data line, until obtaining this data line institute There are the absolute displacement values of Displacement Estimation point；

S5, calculate overall priori value according to the absolute displacement values of each Displacement Estimation point on this data line, and according to Described overall priori value calculates the absolute displacement values of each the Displacement Estimation point on adjacent data line.

As a further improvement on the present invention, some that described echo-signal divides includes Part I and second Part, described s3 step includes:

Choose the Part I of a certain data line, calculate the absolute displacement of each Displacement Estimation point in described Part I Value, and calculate Part I priori value according to the absolute displacement values of each the Displacement Estimation point in described Part I；

According to Part I and Part II position relationship in the axial direction, using Part I priori value as Part II One end value of the region of search of intrinsic displacement estimation point, to calculate each Displacement Estimation point in Part II by searching algorithm Absolute displacement values.

As a further improvement on the present invention, some that described echo-signal divides also includes Part III, described Part III is presented axially between described Part I and Part II, and described s3 step includes:

Calculate Part II priori value according to the absolute displacement values of each the Displacement Estimation point in described Part II；

According to Part I, Part II, Part III position relationship in the axial direction, by Part I priori value and Two part priori values form the region of search of Part III intrinsic displacement estimation point, to be calculated in Part III by searching algorithm The absolute displacement values of each Displacement Estimation point,.

As a further improvement on the present invention, " calculate part priori according to the absolute displacement values of each Displacement Estimation point Value " step specifically includes:

Filter the noise in the absolute displacement values having calculated that；

Compared with prior art, the computational methods of the ultrasonic elastograph imaging shift value of the present invention, can answer for different Become size adaptation ground adjustment hunting zone, greatly reduce amount of calculation, there is very high robustness to avoid mistake expansion simultaneously Dissipate it is ensured that the real-time of elastogram and accuracy.

Brief description

Fig. 1 is the flow chart of the computational methods of ultrasonic elastograph imaging shift value in prior art；

Fig. 2 is the flow chart of the computational methods of ultrasonic elastograph imaging shift value in first embodiment of the invention；

Fig. 3 is the flow chart of the computational methods of ultrasonic elastograph imaging shift value in second embodiment of the invention；

Fig. 4 is the flow chart of the computational methods of ultrasonic elastograph imaging shift value in third embodiment of the invention；

Fig. 5 is single frames echo-signal region division schematic diagram；

Fig. 6 is the structural representation of the computing device of ultrasonic elastograph imaging shift value of the present invention.

Specific embodiment

Describe the present invention below with reference to each embodiment shown in the drawings.But these embodiments are not Limit the present invention, structure, method or change functionally that those of ordinary skill in the art is made according to these embodiments Change and be all contained in protection scope of the present invention.

As shown in Fig. 2 the computational methods of the ultrasonic elastograph imaging shift value in first embodiment of the invention, described side Method comprises the following steps:

Specifically, during quasistatic compression tissue, using the edge such as acoustic pressure power of probe compressing or ultrasonic generation A direction urging tissue, obtain two frame of digital signals in tissue pressurized process with this, accordingly, described data signal is Choose two frame rf echo-signals before and after compressing or iq signal；In specific examples below, will using rf echo-signal as Specific embodiment is introducing.

Accordingly, in conjunction with shown in Fig. 5, Fig. 5 is that single frames echo-signal subregion divides schematic diagram.

Specifically, along the axial direction of single frames echo-signal, i.e. transonic direction, press on every data line respectively Choose multiple Displacement Estimation points according to certain interval.Accordingly, due to the signal that obtains, sampled point is closeer in the axial direction, and displacement General all smaller, therefore, generally choose sampled point at certain intervals as the Displacement Estimation point of this regional area.

After multiple Displacement Estimation points are chosen on every data line, by the Displacement Estimation point on each data wire axially Direction is divided into sizable some.It is understood that in Practical Calculation, can be according to Practical Calculation demand, will The numerical value of the plurality of part specifically sets.Accordingly, for convenience of explanation, as a example being divided into 5 parts in the present embodiment.

Accordingly, 5 parts single frames echo-signal being divided, are followed successively by part 1, part 2, third portion, the 4th Point, the 5th part.

S3, choose at least a portion, calculate the absolute displacement values in this partly each Displacement Estimation point interior, according to every number Calculate the priori value of every data line according to the absolute displacement values on line, and the priori value of every data line is same as calculating The calculating basis of the absolute displacement values of certain a part of or a few part on data line.

Accordingly, if only including Part I and Part II in some of described echo-signal division.

Choose Part I, calculate the absolute displacement values of each Displacement Estimation point in described Part I, and according to described The absolute displacement values of each the Displacement Estimation point in Part I calculate the Part I priori value of every data line；According to A part and Part II position relationship in the axial direction, using the Part I priori value of every data line as same data An end value of the region of search of Part II intrinsic displacement estimation point on line, to be calculated on same data line by searching algorithm The absolute displacement values of each Displacement Estimation point in Part II.

If being separated with least one Part III between described Part I and described Part II.

The of every data line is then calculated according to the absolute displacement values of each the Displacement Estimation point in described Part II Two part priori values；According to Part I, Part II, Part III position relationship in the axial direction, by every data line Part I priori value and Part II priori value form the field of search of Part III intrinsic displacement estimation point on same data line Between, to calculate the absolute displacement values of each Displacement Estimation point in Part III on same data line by searching algorithm.

The like, calculate the absolute displacement values of the Displacement Estimation point of each several part on every data line.

Specifically, in conjunction with shown in Fig. 5, in the process that each data line is calculated with all absolute displacement values that it comprises In, according to first calculating part 1, then the 5th part, third portion, part 2, the order of the 4th part is calculated.Whenever one After the completion of the absolute displacement values of individual part calculate, the intermediate value choosing this part absolute displacement values as priori value, for adjacent part Use；Can also choose and absolute displacement values be rounded with the most absolute displacement values of rear occurrence number as priori value.

Preferably, filter the noise in the absolute displacement values having calculated that；Take in the absolute displacement values after filtering noise Value is as priori value.

Or filter the noise in the absolute displacement values having calculated that；Take and filter noise and round most exhausted of rear occurrence number To shift value as priori value.

Accordingly, choose part 1, each Displacement Estimation point in described part 1 is calculated using conventional search algorithm Absolute displacement values, and calculate every data line according to the absolute displacement values of each the Displacement Estimation point in described part 1 Part I priori value.

In order to reduce hunting zone as far as possible, after obtaining Part I priori value, with the Part I priori obtaining Based on value, the last part chosen on current data line is calculated, and in the present embodiment, chooses the 5th part first and enters Row calculates.Certainly, in Practical Calculation, it is also possible to choose the 2nd on current data line after obtaining Part I priori value Point, third portion or the 4th is partly first calculated, and after obtaining Part I priori value, selection other parts carry out calculating and compare Calculated with choosing last part, the hunting zone in its calculating process can be expanded, increased workload.

Accordingly, based on Part I priori value, calculate the 5th partly in comprise all Displacement Estimation points exhausted To shift value, according to the computational methods of Part I, calculate the Part II priori value on every data line.

Accordingly, the position relationship in the axial direction according to current data line each several part, by the Part I of current data line The region of search of third portion intrinsic displacement estimation point on priori value and Part II priori value formation current data line, with by searching Rope algorithm calculates the absolute displacement values of each Displacement Estimation point in third portion on every data line.

In the specific embodiment of the present invention, describe for convenience, do taking a data line as a example and specifically introduce:

Assume with first Displacement Estimation point as reference zero, choose part 1, obtain the exhausted of each Displacement Estimation point To shift value it is: 1.4,13,1.6,2,31 ... 5.3,9,5.6；As well known to those skilled in the art: described absolute displacement values should It is dull, therefore judges that 13,31,9 is the mistake that the reasons such as noise cause, filter making an uproar in the absolute displacement values having calculated that It is assumed that the intermediate value calculating is 4.2 after sound, then the absolute displacement values 4.2 being calculated using part 1 are as the of current data line A part of priori value.

According to part 1, part 2, third portion, the 4th part, and the 5th part position relationship in the axial direction, will One end value of the region of search as the 5th part intrinsic displacement estimation point for the Part I priori value of current data line, with by searching Rope algorithm calculates on current data line the absolute displacement values of the 5th partly each Displacement Estimation point interior.

Calculate the 5th partly during the absolute displacement values of each displacement estimation point, will open after Part I priori value 4.2 Begin to carry out related operation；Accordingly, the absolute displacement values of the 5th partly each displacement estimation point are respectively as follows: 10.2,10.5 ... .13.5, It is assumed that the intermediate value calculating is 12.2 after noise in filtering the absolute displacement values having calculated that；Then calculated with the 5th part The absolute displacement values 12.2 going out are as the Part II priori value of current data line.

Then calculate third portion, in the absolute displacement values of interval [4.2,12.2] each displacement estimation point calculated above, obtain A series of values are it is assumed that the intermediate value that it calculates is 8.5.The absolute displacement values 8.5 then being calculated using third portion are as current data The Part III priori value of line.

Finally calculate part 2 and the 4th part, respectively in [4.2,8.5], the enterprising row operation in interval of [8.5,12.2], Until obtaining all absolute displacement values in all of data lines.Whole absolute displacement values are obtained using this kind of algorithm, greatly reduces Hunting zone.And due to take be the intermediate value that adjacent part is used as priori value, as long as the majority of this conductors inside Displacement Estimation point is correct, can ensure the validity of this priori value.

The described absolute displacement values that s4, foundation obtain in axial direction carry out gradient algorithm, and obtain corresponding dependent variable According to showing strain data image in the form of pseudo- coloured silk.

In the first embodiment of the present invention, it is only to be concrete Jie as a example the data line in a frame echo-signal Continue, during actual operation, the every data line included in a frame echo-signal is calculated simultaneously, finally by every data Whole described absolute displacement values in line in axial direction carry out gradient algorithm, and obtain corresponding strain data, with pseudo- coloured silk Form shows strain data image.

The computational methods of the ultrasonic elastograph imaging shift value of the present invention, can adjust for different strain size adaptation ground Whole hunting zone, greatly reduces amount of calculation, and due to take be the intermediate value that adjacent part is used as priori value, if this Most displacement estimation points of conductors inside are correct, can ensure the validity of this priori value, have very high robustness Avoid error diffusion it is ensured that the real-time of elastogram and accuracy.

As shown in figure 3, the computational methods of the ultrasonic elastograph imaging shift value in the present invention two embodiment.

Described second embodiment is with the difference of first embodiment, calculates at least initially with conventional search algorithm The absolute displacement values of each the Displacement Estimation point on one data line；Absolute displacement values meter according to each Displacement Estimation point afterwards Calculate priori value, and calculate the absolute displacement values of each the Displacement Estimation point on adjacent data line according to described priori value.

Specifically, the computational methods of the ultrasonic elastograph imaging shift value in second embodiment of the invention include following step Rapid:

S2, multiple Displacement Estimation points are chosen on every data line；

Accordingly, ask for the calculation of the absolute displacement values of each Displacement Estimation point on a data line described in step s2, s3 Method, is asked for using traditional algorithm, here is not described in detail.

S4, calculate priori value according to the absolute displacement values of each Displacement Estimation point, and phase is calculated according to described priori value The absolute displacement values of each the Displacement Estimation point on adjacent data wire.

Accordingly, in conventional search algorithm, filter the noise in the absolute displacement values having calculated that；Take after filtering noise Absolute displacement values intermediate value as priori value.

Accordingly, along the axial direction of data wire, described echo-signal is divided into some, institute is right according to adjacent data line The part answered position relationship in the axial direction and its priori value form the region of search, to calculate adjacent data by searching algorithm The absolute displacement values of each Displacement Estimation point on line.

Specifically, in conjunction with shown in Fig. 5, for example: calculated by conventional search algorithm, the displacement on the first data wire is estimated All absolute displacement values of enumeration are followed successively by: 1.4,13,1.6,2,31 ... 5.3,9,5.6；Afterwards will along the axial direction of data wire Described echo-signal is divided into 5 parts, after division, as well known to those skilled in the art: described absolute displacement values should be dull , therefore judge that 13,31,9 is the mistake that the reasons such as noise cause, after filtering the noise in the absolute displacement values having calculated that, The each several part intermediate value asking for the first data wire is followed successively by: the intermediate value of part 1 is 4.2, and the intermediate value of part 2 is the 6.3, the 3rd The intermediate value divided is the intermediate value of the 8.5, the 4th part is 12.2 for the intermediate value of the 10.3, the 5th part.

By each several part intermediate value of the first data wire, that is, the priori value of each several part of described first data wire it may be determined that The hunting zone of adjacent data line each several part Displacement Estimation point region.Accordingly, calculating is adjacent with described first data wire Current data line in all Displacement Estimation points absolute displacement values when, its use described first data wire priori value algorithm Have many kinds, in the present embodiment, using described first data wire n+1 part priori value as with described first data This algorithm of displacement maximum of the adjacent current data line of line, to ask for the current data line adjacent with described first data wire On all Displacement Estimation points absolute displacement values.And in this manner, continuation order is asked for adjacent with current data line The absolute displacement values of all Displacement Estimation points on other data wires, until obtaining all Displacement Estimation points on whole Frame Absolute displacement values.

Specifically, using this kind of algorithm, ask for each several part displacement on the current data line adjacent with described first data wire The interval that the absolute displacement values of estimation point are distributed is: the interval of part 1 is: [0,6.3], part 2 [4.2,8.5], the 3rd Partly [6.3,10.2], the interval of the 4th part is: [8.5,10.3], the 5th is partially larger than 10.3.

With reference to the relevant range at each several part, directly ask on the current data line adjacent with described first data wire Each displacement estimation point absolute displacement values.And in this manner, continuation order asks for other adjacent with current data line The absolute displacement values of all Displacement Estimation points on data wire, on the whole Frame of acquisition, all Displacement Estimation points is absolute Shift value.

The described absolute displacement values that s5, foundation obtain in axial direction carry out gradient algorithm, and obtain corresponding dependent variable According to showing strain data image in the form of pseudo- coloured silk.

Step s5 technology known to those skilled in the art, here is not described in detail.

The computational methods of the ultrasonic elastograph imaging shift value described in second embodiment of the invention, compared to existing skill Art, equally can greatly reduce amount of calculation for different strain size adaptation ground adjustment hunting zones, and due to taking Be the intermediate value that adjacent part is used as priori value, as long as most displacement estimation points of this conductors inside are correct, can Ensure the validity of this priori value, have very high robustness come to avoid error diffusion it is ensured that elastogram real-time and Accuracy.

Accordingly, as shown in figure 4, Fig. 4 is the calculating side of the ultrasonic elastograph imaging shift value of third embodiment of the invention Method.Described 3rd embodiment is combined with first embodiment and second embodiment.Its principle is: by the first embodiment party Formula asks for the priori value of each several part on wherein 1 data line, calculates each on adjacent data line according to described priori value afterwards The absolute displacement values of Displacement Estimation point.

Accordingly, the computational methods of the ultrasonic elastograph imaging shift value of described Third Way comprise the following steps:

S2, multiple Displacement Estimation points are chosen on every data line, and axially described echo-signal being drawn along data wire It is divided into some.

Accordingly, in conjunction with shown in Fig. 5, Fig. 5 is single frames echo-signal region division schematic diagram.

S3, the absolute displacement of each the Displacement Estimation point being calculated by searching algorithm in certain part in a certain data line Value.

Specifically, each Displacement Estimation point of the part 1 of current data line is calculated by traditional searching algorithm Absolute displacement values.

S4, calculate part priori value according to the absolute displacement values of each Displacement Estimation point, and according to described part priori Value calculates the absolute displacement values of each the Displacement Estimation point in the another part on this data line, until obtaining this data line The absolute displacement values of upper all Displacement Estimation points；

It should be noted that the part priori value now obtaining, simply according to wherein one data line institute having calculated that Ask for.

Accordingly, choose at least a portion of wherein one data line, calculate in this partly interior each Displacement Estimation point Absolute displacement values, calculate the priori value of current data line according to the absolute displacement values on current data line, and by described priori Value is as the calculating basis calculating the absolute displacement values of certain a part of or a few part on current data line.

Choose Part I, calculate the absolute displacement of each Displacement Estimation point in the described Part I in current data line Value, and the Part I of current data line is calculated according to the absolute displacement values of each the Displacement Estimation point in described Part I Priori value；

According to Part I and Part II position relationship in the axial direction, by the Part I priori value of current data line As an end value of the region of search of Part II intrinsic displacement estimation point on current data line, worked as with being calculated by searching algorithm The absolute displacement values of each Displacement Estimation point in Part II on front data wire.

The of current data line is then calculated according to the absolute displacement values of each the Displacement Estimation point in described Part II Two part priori values；Afterwards, the position relationship in the axial direction according to Part I, Part II, Part III, by current data The search of Part III intrinsic displacement estimation point on the Part I priori value of line and Part II priori value formation current data line Interval, to calculate the absolute displacement values of each Displacement Estimation point in Part III on current data line by searching algorithm.

Specifically, in conjunction with shown in Fig. 5, during all absolute displacement values that it comprises to current data line computation, According to first calculating part 1, then the 5th part, third portion, part 2, the order of the 4th part is calculated.Whenever one After the completion of partial absolute displacement values calculate, the intermediate value choosing this part absolute displacement values, as priori value, makes for adjacent part With；Can also choose and absolute displacement values be rounded with the most absolute displacement values of rear occurrence number as priori value.

Accordingly, choose part 1, each Displacement Estimation point in described part 1 is calculated using conventional search algorithm Absolute displacement values, and calculate current data line according to the absolute displacement values of each the Displacement Estimation point in described part 1 Part I priori value.

Accordingly, based on Part I priori value, calculate the 5th partly in comprise all Displacement Estimation points exhausted To shift value, according to the computational methods of Part I, calculate the Part II priori value on current data line.

Accordingly, the position relationship in the axial direction according to current data line each several part, by the Part I of current data line The region of search of third portion intrinsic displacement estimation point on priori value and Part II priori value formation current data line, with by searching Rope algorithm calculates the absolute displacement values of each Displacement Estimation point in third portion on current data line.

Specifically, it is assumed that with first Displacement Estimation point as reference zero, choosing part 1, obtaining in current data line The absolute displacement values obtaining each Displacement Estimation point are: 1.4,13,1.6,2,31 ... 5.3,9,5.6；Those skilled in the art Known: described absolute displacement values should be dull, therefore judge that 13,31,9 is the mistake that the reasons such as noise cause, filter It is assumed that the intermediate value calculating is 4.2 after noise in the absolute displacement values calculating, then the absolute position being calculated with part 1 Shifting value 4.2 is as the Part I priori value of current data line.

Finally calculate part 2 and the 4th part, respectively in [4.2,8.5], the enterprising row operation in interval of [8.5,12.2], Until obtaining all absolute displacement values on current data line.Whole absolute displacement values are obtained using this kind of algorithm, greatly reduces Hunting zone.And due to take be the intermediate value that adjacent part is used as priori value, as long as the majority of this conductors inside Displacement Estimation point is correct, can ensure the validity of this priori value.

Specifically, in conjunction with shown in Fig. 5, by each several part intermediate value of the first data wire, i.e. each portion of described first data wire Point priori value it may be determined that search model with described first data wire adjacent data line each several part Displacement Estimation point region Enclose.Accordingly, when calculating the absolute displacement values of all Displacement Estimation points in the current data line adjacent with described first data wire, The algorithm of its priori value using described first data wire has many kinds, in the present embodiment, using described first data wire N+1 part priori value as the current data line adjacent with described first data wire this algorithm of displacement maximum, with Ask for the absolute displacement values of all Displacement Estimation points on the current data line adjacent with described first data wire.And according to this kind of Method, continuation order asks for the absolute displacement values of all Displacement Estimation points on other data wires adjacent with current data line, Until obtaining the absolute displacement values of all Displacement Estimation points on whole Frame.

The described absolute displacement values that s6, foundation obtain in axial direction carry out gradient algorithm, and obtain corresponding dependent variable According to showing strain data image in the form of pseudo- coloured silk.

Step s6 technology known to those skilled in the art, here is not described in detail.

The computational methods of the ultrasonic elastograph imaging shift value of the present invention, compared to prior art, can answer for different Become size adaptation ground adjustment hunting zone, greatly reduce amount of calculation, and be intermediate value using adjacent part due to take As priori value, as long as most displacement estimation points of this conductors inside are correct, the validity of this priori value can be ensured.

Accordingly, as shown in fig. 6, Fig. 6 is the computing device of ultrasonic elastograph imaging shift value in embodiment of the present invention Structural representation.

In the first embodiment of the present invention, described device includes: sensing module 100, for obtaining before compression organization Two frame echo-signals afterwards, every frame echo-signal includes a plurality of data lines；Processing module 200, for selecting on every data line Take multiple Displacement Estimation points, and described echo-signal is divided into some by the axial direction along data wire；Choose at least a portion, Calculate the absolute displacement values in this partly each Displacement Estimation point interior, calculated often according to the absolute displacement values on every data line The priori value of data line, and using the priori value of every data line as certain a part of or a few portion on the same data line of calculating The calculating basis of the absolute displacement values dividing.

Specifically, described sensing module 100, during quasistatic compression tissue, is oppressed or ultrasonic product using probe Raw acoustic pressure power etc., along a direction urging tissue, obtains two frame of digital signals in tissue pressurized process with this, accordingly, Described data signal is two frame rf echo-signals or the iq signal chosen before and after compressing；In specific examples below, will be with Rf echo-signal to be introduced as specific embodiment.

Specifically, described processing module 200 is along the axial direction of single frames echo-signal, i.e. transonic direction, difference Choose multiple Displacement Estimation points according to certain interval on every data line.Accordingly, due to obtain signal in the axial direction Sampled point is closeer, and displacement is typically all smaller, therefore, generally choose sampled point at certain intervals as the position of this regional area Move estimation point.

After described processing module 200 chooses multiple Displacement Estimation points on every data line, by the position on each data wire Move estimation point and be axially divided into sizable at least 3 parts.It is understood that in Practical Calculation, permissible According to Practical Calculation demand, the numerical value of the plurality of part is specifically set.Accordingly, for convenience of explanation, in the present embodiment As a example being divided into 5 parts.

Described processing module 200 chooses Part I, calculates the absolute position of each Displacement Estimation point in described Part I Shifting value, and first of every data line is calculated according to the absolute displacement values of each the Displacement Estimation point in described Part I Divide priori value；According to Part I and Part II position relationship in the axial direction, by the Part I priori of every data line It is worth the end value as the region of search of Part II intrinsic displacement estimation point on same data line, to calculate by searching algorithm Go out the absolute displacement values of each Displacement Estimation point in Part II on same data line.

Described processing module 200 then calculates according to the absolute displacement values of each the Displacement Estimation point in described Part II The Part II priori value of every data line；According to Part I, Part II, Part III position relationship in the axial direction, The Part I priori value of every data line and Part II priori value are formed Part III intrinsic displacement on same data line The region of search of estimation point, to calculate each Displacement Estimation point in Part III on same data line by searching algorithm Absolute displacement values.

Specifically, in conjunction with shown in Fig. 5, described processing module 200 each data line is calculated with it comprises all absolutely During shift value, according to first calculating part 1, then the 5th part, third portion, part 2, the order of the 4th part Calculated.After the completion of the absolute displacement values of a part calculate, the intermediate value choosing this part absolute displacement values is as elder generation Test value, use for adjacent part；Can also choose and absolute displacement values are rounded with the most absolute displacement values conduct of rear occurrence number Priori value.

Accordingly, described processing module 200 chooses part 1, is calculated every in described part 1 using conventional search algorithm The absolute displacement values of individual Displacement Estimation point, and calculated according to the absolute displacement values of each the Displacement Estimation point in described part 1 Go out the Part I priori value of every data line.

In order to reduce hunting zone as far as possible, described processing module 200 after obtaining Part I priori value, with obtain Part I priori value based on, choose current data line on last part calculated, in the present embodiment, first Choose the 5th part to be calculated.Certainly, in Practical Calculation, it is also possible to choose current data after obtaining Part I priori value Part 2 on line, third portion or the 4th is partly first calculated, and after obtaining Part I priori value, chooses other parts Carry out calculating to compare and choose last part and calculated, the hunting zone in its calculating process can be expanded, increase workload.

Accordingly, described processing module 200 is based on Part I priori value, calculate the 5th partly in the institute that comprises There are the absolute displacement values of Displacement Estimation point, according to the computational methods of Part I, calculate the Part II on every data line Priori value.

Accordingly, described processing module 200 according to current data line each several part position relationship in the axial direction, by current number Form the search of third portion intrinsic displacement estimation point on current data line according to the Part I priori value of line and Part II priori value Interval, to calculate the absolute displacement values of each Displacement Estimation point in third portion on every data line by searching algorithm.

Assume with first Displacement Estimation point as reference zero, described processing module 200 chooses part 1, obtain each The absolute displacement values of individual Displacement Estimation point are: 1.4,13,1.6,2,31 ... 5.3,9,5.6；As well known to those skilled in the art: absolutely Should be dull to shift value, therefore judge that 13,31,9 is the mistakes that the reason such as noise causes, filter have calculated that absolute It is assumed that the intermediate value calculating is 4.2 after noise in shift value, then the absolute displacement values 4.2 being calculated using part 1 are as working as The Part I priori value of front data wire.

Described processing module 200 is according to part 1, part 2, third portion, the 4th part, and the 5th part is in axial direction On position relationship, using the Part I priori value of current data line as the region of search of the 5th part intrinsic displacement estimation point One end value, to calculate on current data line the absolute displacement values of the 5th partly each Displacement Estimation point interior by searching algorithm.

Described processing module 200 calculating the 5th, will be from Part I partly during the absolute displacement values of each displacement estimation point Related operation is proceeded by after priori value 4.2；Accordingly, the absolute displacement values of the 5th partly each displacement estimation point are respectively as follows: 10.2,10.5 ... .13.5, it is assumed that the intermediate value calculating is after the noise in filtering the absolute displacement values having calculated that 12.2；The absolute displacement values 12.2 then being calculated using the 5th part are as the Part II priori value of current data line.

Then calculate third portion, described processing module 200 is in interval [4.2,12.2] each displacement estimation point calculated above Absolute displacement values, obtain a series of values it is assumed that the intermediate value that it calculates is 8.5.The absolute displacement values then being calculated with third portion 8.5 as current data line Part III priori value.

Described processing module 200 finally calculates part 2 and the 4th part, described processing module 200 respectively [4.2, 8.5], the enterprising row operation in the interval of [8.5,12.2], until obtaining all absolute displacement values in all of data lines.Using this kind Algorithm obtains whole absolute displacement values, greatly reduces hunting zone.And be intermediate value using adjacent part due to take As priori value, as long as most displacement estimation points of this conductors inside are correct, the validity of this priori value can be ensured.

Accordingly, the described absolute displacement values that described processing module 200 is additionally operable to according to obtaining in axial direction carry out ladder Degree computing, and obtain corresponding strain data, show strain data image in the form of pseudo- coloured silk.

In the first embodiment of the present invention, it is only to be concrete Jie as a example the data line in a frame echo-signal Continue, during actual operation, the every data line included in a frame echo-signal is calculated simultaneously, finally described process mould Whole described absolute displacement values in every data line are in axial direction carried out gradient algorithm by block 200, and obtain corresponding answering Become data, in the form of pseudo- coloured silk, show strain data image.

In the corresponding present invention two embodiment, the computing device of ultrasonic elastograph imaging shift value includes: sensing module 100, For obtaining two frame echo-signals before and after compression organization, every frame echo-signal includes a plurality of data lines；Processing module 200, uses In multiple Displacement Estimation points are chosen on every data line, and axially described echo-signal is divided into multiple portions along data wire Point；Calculate the absolute displacement values of each the Displacement Estimation point in certain part in a certain data line by searching algorithm；According to The absolute displacement values of each Displacement Estimation point calculate part priori value, calculate each in the another part on this data line The absolute displacement values of Displacement Estimation point, until obtain the shift value of all Displacement Estimation points on this data line；According to this number Absolute displacement values according to each Displacement Estimation point on line calculate overall priori value, and are calculated adjacent according to described overall priori value The absolute displacement values of each the Displacement Estimation point on data wire.

Specifically, described sensing module 100, during quasistatic compression tissue, is oppressed or ultrasonic using probe Acoustic pressure power producing etc., along a direction urging tissue, obtains two frame of digital signals in tissue pressurized process, accordingly with this , described data signal is two frame rf echo-signals or the iq signal chosen before and after compressing；In specific examples below, will To be introduced using rf echo-signal as specific embodiment.

Accordingly, described processing module 200 asks for the calculation of the absolute displacement values of each Displacement Estimation point on a data line Method, is asked for using traditional algorithm, here is not described in detail.

Described processing module 200 is additionally operable in conventional search algorithm, filters making an uproar in the absolute displacement values having calculated that Sound；The intermediate value taking the absolute displacement values after filtering noise is as priori value.

Or described processing module 200 is used for filtering the noise in the absolute displacement values having calculated that；Take and filter noise and take The most absolute displacement values of whole rear occurrence number are as priori value.

Accordingly, described echo-signal is divided into some along the axial direction of data wire by described processing module 200, according to Part corresponding to adjacent data line position relationship in the axial direction and its priori value form the region of search, with by searching algorithm Calculate the absolute displacement values of each Displacement Estimation point on adjacent data line.

Specifically, in conjunction with shown in Fig. 5, for example: described processing module 200 is calculated by conventional search algorithm, first All absolute displacement values of the Displacement Estimation point on data wire are followed successively by: 1.4,13,1.6,2,31 ... 5.3,9,5.6；Afterwards Described echo-signal is divided into 5 parts by the axial direction along data wire, after division, as well known to those skilled in the art: described absolute Shift value should be dull, therefore judges that 13,31,9 is the mistake that the reasons such as noise cause, and filters the absolute position having calculated that After noise in shifting value, each several part intermediate value asking for the first data wire is followed successively by: the intermediate value of part 1 is 4.2, part 2 Intermediate value is 6.3, and the intermediate value of third portion is 12.2 for the intermediate value of the 8.5, the 4th part for the intermediate value of the 10.3, the 5th part.

Described processing module 200 by each several part intermediate value of the first data wire, i.e. each several part of described first data wire Priori value is it may be determined that the hunting zone of adjacent data line each several part Displacement Estimation point region.Accordingly, calculate with described In the adjacent current data line of first data wire during the absolute displacement values of all Displacement Estimation points, it uses described first data wire The algorithm of priori value have many kinds, in the present embodiment, using the priori value conduct of the n+1 part of described first data wire This algorithm of displacement maximum of the current data line adjacent with described first data wire, to ask for and described first data wire phase The absolute displacement values of all Displacement Estimation points on adjacent current data line.And in this manner, described processing module 200 continues Continuous order asks for the absolute displacement values of all Displacement Estimation points on other data wires adjacent with current data line, until obtaining The absolute displacement values of all Displacement Estimation points on whole Frame.

Specifically, described processing module 200 adopts this kind of algorithm, asks for the current data adjacent with described first data wire The interval that on line, the absolute displacement values of each several part Displacement Estimation point are distributed is: the interval of part 1 is: [0,6.3], the 2nd Divide [4.2,8.5], third portion [6.3,10.2], the interval of the 4th part is: [8.5,10.3], the 5th is partially larger than 10.3.

With reference to the relevant range at each several part, described processing module 200 is directly asked for and described first data wire phase The absolute displacement values of each displacement estimation point on adjacent current data line.And in this manner, described processing module 200 continues Order asks for the absolute displacement values of all Displacement Estimation points on other data wires adjacent with current data line, whole until obtaining The absolute displacement values of all Displacement Estimation points on individual Frame.

Described processing module 200 is additionally operable to: in axial direction carry out gradient algorithm according to the described absolute displacement values obtaining, And obtain corresponding strain data, show strain data image in the form of pseudo- coloured silk.

In the corresponding present invention three embodiment, the computing device of ultrasonic elastograph imaging shift value includes: sensing module 100, For obtaining two frame echo-signals before and after compression organization, every frame echo-signal includes a plurality of data lines；Processing module 200 is used In: multiple Displacement Estimation points are chosen on every data line, and multiple along being axially divided into described echo-signal of data wire Part；Calculate the absolute displacement values of each the Displacement Estimation point in certain part in a certain data line by searching algorithm；Root Absolute displacement values according to each Displacement Estimation point calculate part priori value, calculate every in the another part on this data line The absolute displacement values of individual Displacement Estimation point, until obtain the shift value of all Displacement Estimation points on this data line；According to this On data wire, the absolute displacement values of each Displacement Estimation point calculate overall priori value, and calculate phase according to described overall priori value The absolute displacement values of each the Displacement Estimation point on adjacent data wire.

Specifically, described sensing module 100 is used for, during quasistatic compression tissue, oppressing using probe or super Acoustic pressure power that sound produces etc., along a direction urging tissue, obtains two frame of digital signals in tissue pressurized process, phase with this Answer, described data signal is two frame rf echo-signals or the iq signal chosen before and after compressing；In specific examples below, To be introduced using rf echo-signal as specific embodiment.

Specifically, in conjunction with shown in Fig. 5, Fig. 5 is single frames echo-signal region division schematic diagram.

After described processing module 200 chooses multiple Displacement Estimation points on every data line, by the position on each data wire Move estimation point and be axially divided into sizable some.It is understood that in Practical Calculation, can be according to Practical Calculation demand, the numerical value of the plurality of part is specifically set by described processing module 200.Accordingly, for convenience of explanation, As a example being divided into 5 parts in the present embodiment.

Accordingly, 5 parts that single frames echo-signal is divided by described processing module 200, are followed successively by part 1, the 2nd Point, third portion, the 4th part, the 5th part.

Specifically, described processing module 200 calculates the every of the part 1 of current data line by traditional searching algorithm The absolute displacement values of individual Displacement Estimation point.

Described processing module 200 calculates part priori value according to the absolute displacement values of each Displacement Estimation point, and according to Described part priori value calculates the absolute displacement values of each the Displacement Estimation point in the another part on this data line, until obtaining Obtain the absolute displacement values of all Displacement Estimation points on this data line；It should be noted that the part priori value now obtaining, only It is according to striked by wherein one data line having calculated that.

Accordingly, described processing module 200 chooses at least a portion of wherein one data line, calculates partly interior every at this The absolute displacement values of individual Displacement Estimation point, calculate the priori of current data line according to the absolute displacement values on current data line Value, and using described priori value as the calculating base calculating the absolute displacement values of certain a part of or a few part on current data line Plinth.

Described processing module 200 chooses Part I, calculates each displacement in the described Part I in current data line The absolute displacement values of estimation point, and calculated currently according to the absolute displacement values of each the Displacement Estimation point in described Part I The Part I priori value of data wire；

Described processing module 200 position relationship in the axial direction according to Part I and Part II, by current data line Part I priori value as the region of search of Part II intrinsic displacement estimation point on current data line an end value, to pass through Searching algorithm calculates the absolute displacement values of each Displacement Estimation point in Part II on current data line.

Described processing module 200 then calculates according to the absolute displacement values of each the Displacement Estimation point in described Part II The Part II priori value of current data line；

Afterwards, described processing module 200 according to Part I, Part II, Part III position relationship in the axial direction, The Part I priori value of current data line and Part II priori value are formed Part III intrinsic displacement on current data line estimate The region of search of enumeration, to calculate the absolute of each Displacement Estimation point in Part III on current data line by searching algorithm Shift value.

Specifically, in conjunction with shown in Fig. 5, described processing module 200 to current data line computation its comprise all definitely During shift value, according to first calculating part 1, then the 5th part, third portion, part 2, the order of the 4th part is entered Row calculates.After the completion of the absolute displacement values of a part calculate, the intermediate value choosing this part absolute displacement values is as priori Value, uses for adjacent part；Intermediate value can be rounded it is also possible to choose to absolute displacement values as the algorithm of priori value using intermediate value The most absolute displacement values of occurrence number are as priori value afterwards.

Preferably, described processing module 200 filters the noise in the absolute displacement values having calculated that；Take after filtering noise The intermediate value of absolute displacement values is as priori value.

Or described processing module 200 filters the noise in the absolute displacement values having calculated that；Take after filtering noise and rounding The most absolute displacement values of occurrence number are as priori value.

Accordingly, described processing module 200 chooses part 1, is calculated every in described part 1 using conventional search algorithm The absolute displacement values of individual Displacement Estimation point, and calculated according to the absolute displacement values of each the Displacement Estimation point in described part 1 Go out the Part I priori value of current data line.

Accordingly, described processing module 200 is based on Part I priori value, calculate the 5th partly in the institute that comprises There are the absolute displacement values of Displacement Estimation point, according to the computational methods of Part I, calculate the Part II on current data line Priori value.

Accordingly, described processing module 200 according to current data line each several part position relationship in the axial direction, by current number Form the search of third portion intrinsic displacement estimation point on current data line according to the Part I priori value of line and Part II priori value Interval, to calculate the absolute displacement values of each Displacement Estimation point in third portion on current data line by searching algorithm.

Specifically, it is assumed that with first Displacement Estimation point as reference zero in current data line, described processing module 200 selection part 1s, the absolute displacement values obtaining each Displacement Estimation point are: 1.4,13,1.6,2,31 ... 5.3,9,5.6 ；As well known to those skilled in the art: absolute displacement values should be dull, therefore judge that 13,31,9 is that the reasons such as noise cause Mistake, it is assumed that the intermediate value calculating is 4.2 after filtering the noise in the absolute displacement values having calculated that, is then calculated with part 1 The absolute displacement values 4.2 going out are as the Part I priori value of current data line.

Finally calculate part 2 and the 4th part, described processing module 200 respectively in [4.2,8.5], [8.5,12.2] Interval enterprising row operation, until obtaining all absolute displacement values on current data line.Described processing module 200 adopts this kind to calculate Method obtains whole absolute displacement values, greatly reduces hunting zone.And be to be made using the intermediate value of adjacent part due to take For priori value, as long as most displacement estimation points of this conductors inside are correct, the validity of this priori value can be ensured.

Described processing module 200 calculates overall elder generation according to the absolute displacement values of each Displacement Estimation point on this data line Test value, and calculate the absolute displacement values of each the Displacement Estimation point on adjacent data line according to described overall priori value.

Specifically, in conjunction with shown in Fig. 5, described processing module 200 passes through each several part intermediate value of the first data wire, and that is, described the The priori value of each several part of one data wire is it may be determined that the search model of adjacent data line each several part Displacement Estimation point region Enclose.Accordingly, described processing module 200 calculates all Displacement Estimation points in the current data line adjacent with described first data wire Absolute displacement values when, the algorithm of its priori value using described first data wire has many kinds, in the present embodiment, described Processing module 200 adopts the priori value of the n+1 part of described first data wire as adjacent with described first data wire current This algorithm of displacement maximum of data wire, to ask for all displacements on the current data line adjacent with described first data wire The absolute displacement values of estimation point.And in this manner, described processing module 200 continuation order is asked for adjacent with current data line Other data wires on all Displacement Estimation points absolute displacement values, until obtain whole Frame on all Displacement Estimation points Absolute displacement values.

Specifically, described processing module 200 adopts this kind of algorithm, asks for the current data adjacent with described first data wire The interval that on line, the absolute displacement values of each several part Displacement Estimation point are distributed is: the interval of part 1 is: [0,6.3], the 2nd Divide [4.2,8.5], third portion [6.3,10.2], the interval of the 4th part is: [8.5,10.3], the 5th is partially larger than 10.3.Reference Relevant range at each several part, described processing module 200 directly asks for the current data adjacent with described first data wire The absolute displacement values of each displacement estimation point on line.And in this manner, described processing module 200 continuation order is asked for and is worked as The absolute displacement values of all Displacement Estimation points on other adjacent data wires of front data wire, until obtaining institute on whole Frame There are the absolute displacement values of Displacement Estimation point.

Compared with prior art, the computational methods of ultrasonic elastograph imaging shift value of the present invention and device, can be for not Same strain size adaptation ground adjustment hunting zone, greatly reduces amount of calculation, has very high robustness to avoid simultaneously Error diffusion is it is ensured that the real-time of elastogram and accuracy.

For convenience of description, it is divided into various modules to be respectively described with function when describing apparatus above.Certainly, implementing this The function of each module can be realized in same or multiple softwares and/or hardware during application.

As seen through the above description of the embodiments, those skilled in the art can be understood that the application can Mode by software plus necessary general hardware platform to be realized.Based on such understanding, the technical scheme essence of the application On in other words prior art is contributed partly can be embodied in the form of software product, this computer software product Can be saved in preservation medium, such as rom/ram, magnetic disc, CD etc., include some instructions use so that a computer equipment (can be personal computer, Information Push Server, or network equipment etc.) execution each embodiment of the application or reality Apply some partly described methods of mode.

Device embodiments described above are only the schematically wherein said module illustrating as separating component Can be or may not be physically separate, as the part that module shows can be or may not be physics mould Block, you can with positioned at a place, or can also be distributed on multiple mixed-media network modules mixed-medias.Can be selected it according to the actual needs In the purpose to realize present embodiment scheme for some or all of module.Those of ordinary skill in the art are not paying creation Property work in the case of, you can to understand and to implement.

The application can be used in numerous general or special purpose computing system environment or configuration.For example: personal computer, information Push server computer, handheld device or portable set, laptop device, multi-processing module system, be based on microprocessor mould The system of block, set top box, programmable consumer-elcetronics devices, network pc, minicom, mainframe computer, include to take up an official post DCE of what system or equipment etc..

The application can be described in the general context of computer executable instructions, such as program Module.Usually, program module includes execution particular task or the routine realizing particular abstract data type, program, object, group Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these DCEs, by The remote processing devices connected by communication network are executing task.In a distributed computing environment, program module is permissible It is located at and preserve in medium in the local and remote computer including preservation equipment.

It should be understood that although this specification is been described by according to embodiment, but not each embodiment only comprises one Individual independent technical scheme, only for clarity, those skilled in the art should will say this narrating mode of specification As an entirety, the technical scheme in each embodiment can also be through appropriately combined, and forming those skilled in the art can for bright book With the other embodiment understanding.

The a series of detailed description of those listed above is only for the feasibility embodiment of the present invention specifically Bright, they simultaneously are not used to limit the scope of the invention, all equivalent implementations made without departing from skill spirit of the present invention Or change should be included within the scope of the present invention.

Claims

1. a kind of computational methods of ultrasonic elastograph imaging shift value are it is characterised in that the method comprising the steps of:

S2, multiple Displacement Estimation points are chosen on every data line, and axially described echo-signal being divided into along data wire Some；

S3, choose at least a portion, calculate the absolute displacement values in this partly each Displacement Estimation point interior, according to every data line On absolute displacement values calculate the priori value of every data line, and using the priori value of every data line as calculating same number Calculating basis according to the absolute displacement values of certain a part of or a few part on line;

Wherein, some that described echo-signal divides includes Part I and Part II, and described s3 step is specifically wrapped Include:

Choose Part I, calculate the absolute displacement values of each Displacement Estimation point in described Part I, and according to described first Partly the absolute displacement values of each interior Displacement Estimation point calculate the Part I priori value of every data line；

According to Part I and Part II position relationship in the axial direction, using the Part I priori value of every data line as An end value of the region of search of Part II intrinsic displacement estimation point on same data line, same to be calculated by searching algorithm The absolute displacement values of each Displacement Estimation point in Part II on data line.

2. the computational methods of ultrasonic elastograph imaging shift value according to claim 1 are it is characterised in that described Part I It is separated with least one portion between and described Part II between.

3. the computational methods of ultrasonic elastograph imaging shift value according to claim 2 are it is characterised in that described echo-signal The some dividing also includes Part III, described Part III be presented axially in described Part I and Part II it Between；

Described s3 step includes:

Calculate the Part II of every data line according to the absolute displacement values of each the Displacement Estimation point in described Part II Priori value；

According to Part I, Part II, Part III position relationship in the axial direction, the Part I of every data line is first Test value and Part II priori value forms the region of search of Part III intrinsic displacement estimation point on same data line, with by searching Rope algorithm calculates the absolute displacement values of each Displacement Estimation point in Part III on same data line.

4. the computational methods of ultrasonic elastograph imaging shift value according to claim 1 are it is characterised in that " according to every number Calculating the priori value of every data line according to the absolute displacement values on line " step specifically includes:

Filter the noise in the absolute displacement values having calculated that；

Take the intermediate value of the absolute displacement values after filtering noise or round the most absolute displacement values of rear occurrence number as priori value.

5. a kind of computational methods of ultrasonic elastograph imaging shift value are it is characterised in that the method comprising the steps of:

S2, multiple Displacement Estimation points are chosen on every data line；

S4, calculate priori value according to the absolute displacement values of each Displacement Estimation point, and consecutive number is calculated according to described priori value Absolute displacement values according to each the Displacement Estimation point on line；

Wherein, step is concrete " to calculate the absolute displacement values of each the Displacement Estimation point on adjacent data line according to described priori value " Including:

Described echo-signal is divided into some by the axial direction along data wire, by the portion corresponding to adjacent data line juxtaposed position Point and its priori value form the region of search, so that the absolute of each Displacement Estimation point on adjacent data line is calculated by searching algorithm Shift value.

6. the computational methods of ultrasonic elastograph imaging shift value according to claim 5 are it is characterised in that " according to each position The absolute displacement values moving estimation point calculate priori value " step specifically includes:

Filter the noise in the absolute displacement values having calculated that；

7. a kind of computational methods of ultrasonic elastograph imaging shift value are it is characterised in that the method comprising the steps of:

S3, the absolute displacement values of each the Displacement Estimation point being calculated by searching algorithm in certain part in a certain data line；

S4, calculate part priori value according to the absolute displacement values of each Displacement Estimation point, and according to described part priori value meter Calculate the absolute displacement values of each the Displacement Estimation point in another part of this data line, until obtaining all positions of this data line Move the absolute displacement values of estimation point；

8. the computational methods of ultrasonic elastograph imaging shift value according to claim 7 are it is characterised in that described echo-signal The some dividing includes Part I and Part II, and described s3 step includes:

Choose the Part I of a certain data line, calculate the absolute displacement values of each Displacement Estimation point in described Part I, And calculate Part I priori value according to the absolute displacement values of each the Displacement Estimation point in described Part I；

According to Part I and Part II position relationship in the axial direction, using Part I priori value as position in Part II Move an end value of the region of search of estimation point, so that the absolute of each Displacement Estimation point in Part II is calculated by searching algorithm Shift value.

9. the computational methods of ultrasonic elastograph imaging shift value according to claim 8 are it is characterised in that described Part I It is separated with least one portion between and described Part II between.

10. the computational methods of ultrasonic elastograph imaging shift value according to claim 9 are it is characterised in that described echo is believed Number divide some also include Part III, described Part III is presented axially in described Part I and Part II Between, described s3 step includes:

According to Part I, Part II, Part III position relationship in the axial direction, by Part I priori value and second Priori value is divided to form the region of search of Part III intrinsic displacement estimation point, to calculate each in Part III by searching algorithm The absolute displacement values of Displacement Estimation point.

The computational methods of 11. ultrasonic elastograph imaging shift values according to claim 7 are it is characterised in that " according to each position The absolute displacement values moving estimation point calculate part priori value " step specifically includes:

Filter the noise in the absolute displacement values having calculated that；