CN110368030A - Ultrasonic elastograph imaging device and method - Google Patents

Abstract

This application discloses a kind of ultrasonic elastograph imaging device and methods, it selectes the shearing wave that tissue location generates several specific frequency in human body using acoustic radiation force impact method, and calculate shear wave velocity, establish velocity-frequency straight line, speed conditions of the different frequency shearing wave by identical tissue when can be observed, to which the slope value qualitative reaction by calculating the straight line goes out the sticky situation of tissue, does not need to know tissue density values as prior art.

Description

Ultrasonic elastograph imaging device and method

Technical field

The present invention relates to ultrasonic imaging technique fields, and in particular to a kind of ultrasonic elastograph imaging device and method, the deviceAnd method is based on acoustic radiation force impact and realizes that biological tissue's adhesive characteristics are extracted.

Background technique

Shearing wave elastography is to impact biological tissue by transmitting acoustic radiation to apply excitation, utilizes Mach cone principleThe shearing wave for generating sufficient intensity in the tissue passes through inspection since shearing wave is variant in the organizational communication speed of different hardnessThe soft or hard situation of the position can be reflected in the spread speed of different location indirectly by surveying shearing wave, finally by shear wave velocity meterTissue elasticity modulus hardness value is calculated, during which, such as 20,000 frames ultrahigh speed imaging system per second captures, tracking is sheared by reachingWave reaches the accuracy of surveyed shear wave velocity such as 1mm/s, then goes out tissue elasticity figure with Color Coding Technology real-time display.Using shearing wave elastography, the real-time hardness quantitative detection of biological tissue may be implemented, be clinical judgment lesion tissue feelingsCondition provides foundation.

Major part elastograph imaging method is carried out in the case where biological tissue is assumed to be ideal linearity elasticity at present,However, the calculating of this elastograph imaging method has ignored the adhesive characteristics feelings of tissue since biological tissue has viscoelasticity characteristicCondition.Moreover, there is also certain relationships, such as alcoholic liver, fatty liver for lesion situation and the adhesive characteristics variation of tissue.So increasingAdd adhesive characteristics to calculate to be of great significance as the supplement of elastogram.Meanwhile in uniform pure elastic fluid, shearingThe spread speed of wave and the frequency of shearing wave are unrelated, i.e. the shearing velocity of wave propagation of different frequency is identical；But special containing toughnessIn the tissue of sign, the propagation of shearing wave will receive the influence of viscosity.Organizational coherence is bigger, and shearing wave frequency rate is higher, spread speedIt is lower；Conversely, shearing wave frequency rate is lower, spread speed is higher.Therefore, increase when carrying out elastogram to tissue and consider that viscosity is specialThe case where sign, can provide more comprehensively and more true organizational information for clinical diagnosis.

The calculating of the adhesive characteristics situation of tissue needs to know the frequency values and Media density value of shearing wave.Existing method mentionsThe shearing wave for stablizing single-frequency is generated by external mechanical oscillation out, is propagated from body surface to in-vivo tissue, while Media densityValue takes 1000 gram/cubic metre of approximation to carry out adhesive characteristics calculating.Although this method can generate more stable single-frequency and cutWave is cut, but since production method is external mechanical oscillation, shearing wave is very big from body surface to internal decaying, and shearing wave is propagated throughSituations such as division can also be generated by a variety of media in journey, so all bigger to position selection and propagation depth limit.Also,The Media density actual value of in-vivo tissue is not all 1000 grams/cubic metre yet, and different tissues are all used with this density value meetingIntroduce error.So existing method is larger to the sticky situation assessment limitation of tissue, accuracy is low.

Summary of the invention

The application provides a kind of ultrasonic elastograph imaging device and method, wherein realizing biological tissue based on acoustic radiation force impactAdhesive characteristics are extracted, without assuming the density value of tissue as prior art.

According to a first aspect of the present application, this application provides a kind of ultrasonic elastograph imaging devices, comprising:

Ultrasonic probe comprising at least one energy converter carries out multiple groups acoustically radiating for the area-of-interest to biological tissueIt penetrates power and focuses impact, to generate the shearing wave of different frequency, be also used to be returned along the detection position line acquisition chosen by biological tissueThe ultrasonic echo for the multiple groups frequency returned；

Transmitting and reception sequence control module, for generating transmitting sequence and the control of control respective transducer transmitting ultrasonic waveRespective transducer processed receives the reception sequence of ultrasonic echo；

Echo Processing module obtains the corresponding multiple groups frequency in each detection position for handling the ultrasonic echoEcho data；

Transmitting and receiving module, respectively with ultrasonic probe, transmitting and receive sequence control module and Echo Processing moduleElectrical connection, for transmitting with the transmitting sequence that sequence control module generates is received to be transferred to ultrasonic probe, and by ultrasonic probeReceived ultrasonic echo is transferred to Echo Processing module；

Data processing module is connect with Echo Processing module, estimates operation, shape for carrying out deformation according to echo dataAt the corresponding deformation-time graph in each detection position, the deformation-time graph is related to the relationship that deformation values change over time,Shear-wave velocity estimation is carried out according to the deformation-time graph, and is fitted and obtains shear wave velocity-frequency straight line, passes through instituteIt states shear wave velocity-frequency straight line slope and determines adhesive characteristics reference value, and then calculate Resilience Results；

Display module, for showing the echo data and/or the Resilience Results.

Further, the ultrasonic probe carries out first group of sound by preset pulse train to the area-of-interestRadiant force focuses impact, generates the shearing wave of first frequency, acquires the returned by biological tissue along the detection position line chosenOne frequency ultrasound echo；It repeats previous process at least twice, obtains the ultrasonic echo of at least three groups different frequencies.

Further, the data processing module estimates echo data progress deformation using normalized crosscorrelation operationMeter, forms deformation-time graph of each detection position, and records peak value and its corresponding time in every deformation time graphValue carries out straight line fitting according to the peak value and its corresponding time value and check bit glove reason distance value, is fittedWhen it is m- apart from straight line, when seeking this m- slope apart from straight line inverse to obtain shear wave velocity angle value, in conjunction with corresponding frequencyThus rate obtains shear wave velocity-frequency straight line.

According to a second aspect of the present application, this application provides a kind of ultrasonic elastograph imaging methods, comprising:

Transmitting and receiving step: the focusing of multiple groups acoustic radiation force is carried out by area-of-interest of the ultrasonic probe to biological tissueImpact surpassing for the multiple groups frequency to generate the shearing wave of different frequency, returned along the detection position line acquisition chosen by biological tissueSound echo；

Echo Processing step: Echo Processing is carried out to the ultrasonic echo, obtains the corresponding multiple groups frequency in each detection positionEcho data；

Data processing step: carrying out deformation according to the echo data and estimate operation, forms the corresponding shape in each detection positionChange-time graph, the deformation-time graph are related to the relationship that deformation values change over time, according to the deformation-time graphShear-wave velocity estimation is carried out, and is fitted and obtains shear wave velocity-frequency straight line, passes through the shear wave velocity-frequency straight lineSlope determine adhesive characteristics reference value, and then calculate Resilience Results；

Export echo data and/or the Resilience Results.

Further, the transmitting and receiving step include: to the area-of-interest by preset pulse trainIt carries out first group of acoustic radiation force and focuses impact, generate the shearing wave of first frequency, acquire along the detection position line chosen by biologyOrganize the first frequency ultrasonic echo returned；Repeat previous process at least twice, the ultrasound for obtaining at least three groups different frequencies is returnedWave.

Further, in the data processing step, the echo data is carried out using normalized crosscorrelation operationDeformation estimation.

Further, the data of the adjacent acquisition on the same detection position line are normalized mutually using following formulaRelated operation:

In formula, s_rFor the echo data of current ultrasonic echo, s_dFor the echo data of ultrasonic echo next time, R is cross-correlationCoefficient, τ are the ultrasonic echo time delay of corresponding shift value, and T is signal segment cycle length；Wherein, the maximum value institute for seeking R is rightThe result τ answered is adjacent echoes time delay, calculates deformation values u according to the following formula_τ: u_ττ/2=c, c is transonic in formulaSpeed.

Further, in the data processing step, deformation-time graph of each detection position is formed, and is recorded everyPeak value and its corresponding time value in deformation time graph.

Further, in the data processing step, according to the peak value and its corresponding time value and check bitGlove manage distance value carry out straight line fitting, be fitted when it is m- apart from straight line, seek the m- slope apart from straight line when thisInverse is to obtain shear wave velocity angle value, in conjunction with corresponding frequency, thus obtains shear wave velocity-frequency straight line.

According to the third aspect of the application, the application provides a kind of computer readable storage medium, including program, the programIt can be executed by processor to realize above-mentioned method.

It is selected in human body using acoustic radiation force impact method according to the ultrasonic elastograph imaging device and method of above-described embodimentTissue location generates the shearing wave of several specific frequency, and calculates shear wave velocity, establishes velocity-frequency straight line, Ke YiguanSpeed conditions of the different frequency shearing wave by identical tissue when are examined, thus, by the slope value qualitative reaction for calculating the straight lineThe sticky situation of tissue out, does not need to know tissue density values as prior art.

Detailed description of the invention

Fig. 1 is the structural schematic diagram of the ultrasonic elastograph imaging device of the embodiment of the present application；

Fig. 2 is flow processing schematic diagram involved in the ultrasonic elastograph imaging method of the embodiment of the present application；

Fig. 3 is the schematic diagram of acoustic radiation force shock pulse sequence in the embodiment of the present application；

Fig. 4 is in the embodiment of the present application along the schematic diagram of detection station acquisition ultrasound echo signal；

Fig. 5 is in the embodiment of the present application along the deformation time curve signal of the collected each detection position of the detection position lineFigure；

Fig. 6 is schematic diagram of the fit time-apart from straight line in the embodiment of the present application；

Fig. 7 is the schematic diagram that velocity-frequency straight line is fitted in the embodiment of the present application.

Specific embodiment

Below by specific embodiment combination attached drawing, invention is further described in detail.Wherein different embodimentsMiddle similar component uses associated similar element numbers.In the following embodiments, many datail descriptions be in order toThe application is better understood.However, those skilled in the art can recognize without lifting an eyebrow, part of featureIt is dispensed, or can be substituted by other elements, material, method in varied situations.In some cases, this ShenPlease it is relevant it is some operation there is no in the description show or describe, this is the core in order to avoid the application by mistakeMore descriptions are flooded, and to those skilled in the art, these relevant operations, which are described in detail, not to be necessary, theyRelevant operation can be completely understood according to the general technology knowledge of description and this field in specification.

It is formed respectively in addition, feature described in this description, operation or feature can combine in any suitable wayKind embodiment.Meanwhile each step in method description or movement can also can be aobvious and easy according to those skilled in the art instituteThe mode carry out sequence exchange or adjustment seen.Therefore, the various sequences in the description and the appended drawings are intended merely to clearly describe a certainA embodiment is not meant to be necessary sequence, and wherein some sequentially must comply with unless otherwise indicated.

It is herein component institute serialization number itself, such as " first ", " second " etc., is only used for distinguishing described object,Without any sequence or art-recognized meanings.And " connection ", " connection " described in the application, unless otherwise instructed, include directly andIt is indirectly connected with (connection).

Shearing wave measurement method based on acoustic emanation power passes through ultrasonic acoustic radiation force and generates shearing wave inside biological tissuePropagation, then by inside biological tissue generate shearing wave propagation and detect its propagation parameter (such as spread speed)Method reflects the nonhomogeneous hardness between biological tissue.

Referring to FIG. 1, ultrasonic elastograph imaging device 100 includes ultrasonic probe 101, transmitting and receives sequence control module102, transmitting and receiving module 103, Echo Processing module 104, data processing module 105 and display module 106.Transmitting and receptionSequence control module 102 is connect by transmitting and receiving module 103 with 101 signal of ultrasonic probe, and ultrasonic probe 101 passes through transmittingIt is connect with receiving module 103 with 104 signal of Echo Processing module, the output end and data processing module of Echo Processing module 104105 connections, the output end of data processing module 105 are connect with display module 106.

Ultrasonic probe 101 includes at least one energy converter, and energy converter is used to emit according to electric signal ultrasonic wave, or will receiveUltrasonic wave be transformed to electric signal.When there are multiple energy converters, multiple energy converters can be arranged in a row, such as linear array probe；Alternatively, multiple energy converters can be arranged into two-dimensional matrix, for example, two-dimensional ultrasound probe or face battle array probe.One of ultrasonic probe 101Energy converter is divided to can be used for emitting to biological tissue 108 ultrasonic wave, another part energy converter can be used for receiving to be returned by biological tissueUltrasonic echo.In the embodiment of the present application, ultrasonic probe 101 is used to carry out multiple groups acoustically radiating to the area-of-interest of biological tissueIt penetrates power and focuses impact, to generate the shearing wave of different frequency, be also used to be returned along the detection position line acquisition chosen by biological tissueThe ultrasonic echo for the multiple groups frequency returned.Wherein, area-of-interest can be selected by user, such as when display screen display ultrasound imageWhen, user can select area-of-interest on the ultrasound image, and ultrasonic elastograph imaging device 100 can be according to selected area-of-interestIt calculates transmitting and receives sequence.In the embodiment having, ultrasonic elastograph imaging device defaults the tissue table that ultrasonic probe is contactedSubcutaneous setting range is area-of-interest.

Transmitting is used to generate transmitting sequence and reception sequence with sequence control module 102 is received, and transmitting sequence is for providingThe energy converter number of transmitting in ultrasonic probe 101 and to biological tissue transmitting ultrasonic wave parameter (such as amplitude, frequency, hair waveNumber, hair wave angle degree, wave mode etc.), it receives sequence and is used to provide energy converter number and its reception of reception in ultrasonic probe 101The parameter (such as received angle, depth etc.) of echo.When ultrasonic probe 101 is fixed on some position in biological tissue 108Afterwards, corresponding energy converter emits certain amplitude, frequency and angle into biological tissue 108 according to transmitting sequence in ultrasonic probe 101The ultrasonic wave of degree, ultrasonic wave in biological tissue 108 can Formation cross-section, while also can receive the ultrasonic echo in the section.

Transmitting and receiving module 103 are connected to ultrasonic probe and transmitting and receive sequence control module 102, Echo Processing mouldBetween block 104, for transmitting to be transferred to ultrasonic probe 101 with the transmitting sequence for receiving sequence control module 102, and will be ultrasonic101 received ultrasonic echos of popping one's head in are transferred to Echo Processing module 104.

Echo Processing module 104 is used to handle ultrasonic echo, such as is filtered, amplifies to ultrasonic echo, waveThe processing such as Shu Hecheng.

Data processing module 105 receives 104 output end of echo processing module treated echo-signal, and is calculated using relatedParameter required for method obtains or image.In the embodiment of the present application, data processing module 105 can carry out shape according to echo dataBecome estimation operation, formed the corresponding deformation-time graph in each detection position (that reflects the relationships that deformation values change over time),Shear-wave velocity estimation is carried out according to deformation-time graph, and is fitted and obtains shear wave velocity-frequency straight line, passes through shearing waveThe slope of velocity-frequency straight line determines adhesive characteristics reference value, and then calculates Resilience Results.This point will be apparent from below.

Display module 106 is used as output module, and display module 106 receives the various visual of the output of data processing module 105Change data, various images, figure, chart, text or data is shown in the display interface, including various elastic images.?In some embodiments, data processing module 105 can directly export calculated Resilience Results after calculating Resilience ResultsIt is shown to display module.In some embodiments, display module can be touch display screen or multiple displays etc..

Based on the application previous embodiment, another embodiment of the present invention additionally provides a kind of ultrasonic elastograph imaging method,The adhesive characteristics of shock analysis biological tissue are focused based on acoustic radiation force.This method includes transmitting and receiving step, Echo ProcessingStep, data processing step and output step.In transmitting and receiving step, by ultrasonic probe to the interested of biological tissueRegion carries out multiple groups acoustic radiation force and focuses impact to generate the shearing wave of different frequency, acquires along the detection position line chosen by giving birth toThe ultrasonic echo for the multiple groups frequency that object tissue returns.In Echo Processing step, Echo Processing is carried out to the ultrasonic echo of acquisition,Obtain the echo data of the corresponding multiple groups frequency in each detection position.In data processing step, deformation is carried out according to echo dataEstimate operation, form the corresponding deformation-time graph in each detection position, shear-wave velocity is carried out according to deformation-time graph and is estimatedMeter, and be fitted and obtain shear wave velocity-frequency straight line, adhesive characteristics are determined by shear wave velocity-frequency straight line slopeReference value, and then calculate Resilience Results.In output step, echo data and/or Resilience Results are exported.

Ultrasonic elastograph imaging method chooses elasticity in biological tissue elasticity measurement process in one of the embodiments,Viewing area then proceedes to process flow as shown in Figure 2, i.e. following step 201 to step 206.

In step 201, according to the viewing area position of selection, acoustic radiation force is carried out according to preset pulse trainImpact is focused, pulse train can be for as shown in Figure 3.Acoustic radiation force shock pulse length is chosen as such as 50us, acoustic radiationPower shock pulse repetition rate Push_Prf is chosen as such as 100Hz, and acoustic radiation force impact focus quantity is chosen as such as 3 cokesPoint, acquisition pulse repetition rate are chosen as such as 2kHz, and acquisition pulse quantity is chosen as such as 200 times.Fig. 4 is shown along detectionThe position line acquires ultrasound echo signal.By this step, the shearing wave that frequency is Push_Prf can produce.

In step 202, repeat step 201 process at least twice, acoustic radiation force shock pulse repetition rate Push_Prf is chosen as such as 200Hz, 300Hz, and acquisition pulse repetition rate is chosen as such as 2kHz, 2.5kHz, acquisition pulse quantityIt is chosen as such as 200 times.It can produce shearing wave the returning by same tissue location of at least three groups different frequencies by step 202Wave number evidence.

In step 203, echo data is set according to each check bit, carries out deformation and estimates operation.Deformation algorithm for estimating is optionalSelect such as normalized crosscorrelation operation.For example, as shown in following formula (1): to the same detection position line, the data of adjacent acquisition intoRow normalized crosscorrelation operation, wherein s_rFor current ultrasonic rf echo signal, s_dFor ultrasonic radio frequency echo-signal next time, R isCross-correlation coefficient, τ are the ultrasonic echo time delay of corresponding shift value, and T is signal segment cycle length.

Seeking result τ corresponding to the maximum value of R is adjacent echoes time delay, and deformation values u is calculated according to formula (2)_τ,C is ultrasonic propagation velocity in formula.

u_ττ/2=c (2)

In step 204, the deformation-time graph for making each detection position, as shown in figure 5, it illustrates formed whereinOne depth respectively detects deformation-time graph of position.Every deformation-time graph peak value in search graph, and record peak valueCorresponding lateral coordinates value, i.e. time value.

In step 205, time value is corresponded to according to peak value in step 204 and check bit glove manages distance value, make fittingWhen it is m- apart from straight line, as shown in Figure 6；The inverse of the m- slope apart from straight line, can be obtained shear wave velocity angle value when seeking.

In step 206, step 203-205 is repeated, the shear wave velocity estimation of at least 3 groups different frequencies is finally obtainedValue.In conjunction with corresponding frequency, velocity-frequency straight line can be made, as shown in fig. 7, straight slope can be found out in turn.The slope valueReference value i.e. as organizational coherence feature.

By abovementioned steps as can be seen that first passing through during analyzing the adhesive characteristics of biological tissue to tissueSelected location carries out first group of acoustic radiation force and focuses impact, generates the shearing wave of the first frequency；Then it is seen by ultrasound acquisitionExamine the echo-signal of region designated position；It repeats above-mentioned acoustic radiation force impact and echo acquirement process at least twice, finally obtainsThe shearing wave of at least three groups different frequencies passes through the echo data in specified region；The echo-signal of each detection position continues deformationEstimate operation, deformation-time graph is made according to deformation estimated result；The velocity of wave of shearing wave is carried out according to deformation-time graphEstimation；Velocity amplitude of the different shearing wave frequency rates by same tissue when is finally obtained, and linearly fits shear wave velocity-frequencyStraight line acquires straight slope.This straight-line frequency reflects variation of the shearing wave by speed when same tissue location and frequency and closesSystem, qualitative can reflect the sticky situation of tissue.That is, slope is bigger, organizational coherence is bigger；Conversely, slope is closer to zero, groupIt is smaller to knit viscosity.To can analyze out organizational coherence based on acoustic radiation force impact when carrying out elastogram to biological tissueFeature enables elastogram more fully and is truly reflected out organizational information, and then can be clinical judgment lesion tissue feelingsCondition provides more accurate foundation.

In device and method provided by the embodiments of the present application, by selected location carry out multiple groups acoustic radiation force impact andEcho acquirement generates the shearing wave of different frequency in target area, and collects returning for tissue regions of the shearing wave by observationThen wave signal obtains deformation-time graph of each detection position by carrying out cross-correlation calculation to echo data, according to shapeChange-time graph calculates shear wave velocity angle value；It is straight that velocity-frequency is finally fitted according to the shear wave velocity of different frequencyLine calculates the slope of straight line to obtain the adhesive characteristics value of tissue.This method is selected using acoustic radiation force impact method in human bodyDetermine tissue location and generate the shearing wave of several specific frequencies, and calculate shear wave velocity, establishes velocity-frequency straight line, Ke YiguanSpeed conditions of the different frequency shearing wave by identical tissue when are examined, thus, by the slope value qualitative reaction for calculating the straight lineThe sticky situation of tissue out, does not need to know tissue density values as prior art.

It is illustrated herein by reference to various example embodiments.However, those skilled in the art will recognize that, it is not taking offIn the case where from broad scope hereof, exemplary embodiment can be made a change and be corrected.For example, various operating procedures and being used forThe component of operating procedure is executed, according to the specific application or can consider associated with the operation of system any amount of costFunction is realized in different ways.

In addition, as understood by those skilled in the art, the principle of this paper can reflect in computer readable storage mediumOn computer program product in, which is preinstalled with computer readable program code.Any tangible, nonvolatileThe computer readable storage medium of property can all be used, including magnetic storage apparatus (hard disk, floppy disk etc.), optical storage apparatus (CD-ROM, DVD, BluRay disk etc.), flash memory and/or the like.These computer program instructions can be loaded into general purpose computer,To form machine on special purpose computer or other programmable data processing devices, so that these are on computers or other are programmableThe device for realizing specified function can be generated in the instruction executed on data processing equipment.These computer program instructions can also be withIt is stored in computer-readable memory, which can indicate computer or the processing of other programmable datasEquipment is run in a particular manner, and the instruction being stored in computer-readable memory in this way can form a manufacture,Realization device including realizing specified function.Computer program instructions can also be loaded at computer or other programmable datasIt manages in equipment, to execute series of operation steps on a computer or other programmable device to generate a computer realizationProcess, instruction executed on a computer or other programmable device is provided for realizing the step of specified functionSuddenly.

Although having been illustrated with the principle of this paper in various embodiments, many is especially suitable for specific environment and behaviourThe modification for structure, arrangement, ratio, element, material and the component being required can be in the principle and range for not departing from present disclosureIt uses.The above modification and other change or corrections will be comprised within this context.

Aforementioned illustrate is described with reference to various embodiments.However, it would be recognized by those skilled in the art that canTo carry out various modifications and variations in the case where not departing from the range of present disclosure.Therefore, for present disclosure the considerations of, will be saidBright property and not restrictive in the sense that, and all such modifications will be all included within its scope.Equally, about eachThe advantages of kind embodiment, other advantages and solution to the problem are as already identified above.However, the solution party of benefit, advantage, problemCase and it is any can generate these element, or the solution for becoming more apparent it be all not necessarily to be construed as it is crucial, mustIt is needing or necessary.Term " includes " used herein and its any other variant, coming under nonexcludability includes, and wraps in this wayProcess, method, article or the equipment for including element list not only include these elements, further include not expressly listed or are not belonging to thisProcess, method, system, the other element of article or equipment.In addition, term " coupling " used herein and its any otherVariant all refers to physical connection, electrical connection, magnetic connects, light connects, communication connection, function connects and/or any other connection.

People with art technology, can be to upper it will be recognized that without departing from the basic principles of the present inventionThe details for stating embodiment carries out many changes.Therefore, the scope of the present invention should be determined according to following following claims.

Claims (10)

1. a kind of ultrasonic elastograph imaging device characterized by comprising

Ultrasonic probe comprising at least one energy converter carries out multiple groups acoustic radiation force for the area-of-interest to biological tissueImpact is focused, to generate the shearing wave of different frequency, is also used to be returned along the detection position line acquisition chosen by biological tissueThe ultrasonic echo of multiple groups frequency；

Transmitting and reception sequence control module, for generating the transmitting sequence and control phase of control respective transducer transmitting ultrasonic waveEnergy converter is answered to receive the reception sequence of ultrasonic echo；

Echo Processing module obtains returning for each corresponding multiple groups frequency in detection position for handling the ultrasonic echoWave number evidence；

Transmitting and receiving module are electrically connected with ultrasonic probe, transmitting and reception sequence control module and Echo Processing module respectivelyIt connects, for transmitting to be transferred to ultrasonic probe with the transmitting sequence that sequence control module generates is received, and ultrasonic probe is receivedUltrasonic echo be transferred to Echo Processing module；

Data processing module is connect with Echo Processing module, is estimated operation for carrying out deformation according to echo data, is formed eachCorresponding deformation-the time graph in position is detected, the deformation-time graph is related to the relationship that deformation values change over time, according toThe deformation-time graph carries out shear-wave velocity estimation, and is fitted and obtains shear wave velocity-frequency straight line, is cut by describedIt cuts wave velocity-frequency straight line slope and determines adhesive characteristics reference value, and then calculate Resilience Results；

Display module, for showing the echo data and/or the Resilience Results.

2. device as described in claim 1, which is characterized in that the ultrasonic probe is to the area-of-interest by presettingPulse train carry out first group of acoustic radiation force and focus impact, the shearing wave of first frequency is generated, along the detection position line chosenAcquire the first frequency ultrasonic echo returned by biological tissue；It repeats previous process at least twice, obtains at least three groups of different frequenciesThe ultrasonic echo of rate.

3. device as claimed in claim 2, which is characterized in that the data processing module uses normalized crosscorrelation operation pairThe echo data carries out deformation estimation, forms deformation-time graph of each detection position, and records every deformation time graphMiddle peak value and its corresponding time value are carried out according to the peak value and its corresponding time value and check bit glove reason distance valueStraight line fitting, be fitted when it is m- apart from straight line, when seeking this m- slope apart from straight line inverse to obtain shearing waveThus velocity amplitude obtains shear wave velocity-frequency straight line in conjunction with corresponding frequency.

4. a kind of ultrasonic elastograph imaging method characterized by comprising

Transmitting and receiving step: multiple groups acoustic radiation force is carried out by area-of-interest of the ultrasonic probe to biological tissue and focuses impactTo generate the shearing wave of different frequency, the ultrasound of the multiple groups frequency returned along the detection position line acquisition chosen by biological tissue is returnedWave；

Echo Processing step: carrying out Echo Processing to the ultrasonic echo, obtains returning for each corresponding multiple groups frequency in detection positionWave number evidence；

Data processing step: carrying out deformation according to the echo data and estimate operation, formed each detection position it is corresponding deformation-whenHalf interval contour, the deformation-time graph are related to the relationship that deformation values change over time, and are carried out according to the deformation-time graphShear-wave velocity estimation, and be fitted and obtain shear wave velocity-frequency straight line, it is oblique by the shear wave velocity-frequency straight lineRate determines adhesive characteristics reference value, and then calculates Resilience Results；

Export step: output echo data and/or the Resilience Results.

5. method as claimed in claim 4, which is characterized in that described emit with receiving step includes: to the region of interestDomain carries out first group of acoustic radiation force by preset pulse train and focuses impact, generates the shearing wave of first frequency, and edge is chosenThe detection position line acquire the first frequency ultrasonic echo that is returned by biological tissue；Repeat previous process at least twice, obtain toThe ultrasonic echo of few three groups of different frequencies.

6. method as claimed in claim 5, which is characterized in that in the data processing step, using normalized crosscorrelationOperation carries out deformation estimation to the echo data.

7. method as claimed in claim 6, which is characterized in that adopted to the data of the adjacent acquisition on the same detection position lineComputing cross-correlation is normalized with following formula:

In formula, s_rFor the echo data of current ultrasonic echo, s_dFor the echo data of ultrasonic echo next time, R is cross correlationNumber, τ are the ultrasonic echo time delay of corresponding shift value, and T is signal segment cycle length；

Wherein, seeking result τ corresponding to the maximum value of R is adjacent echoes time delay, calculates deformation values according to the following formulau_τ:

u_ττ/2=c

C is ultrasonic propagation velocity in formula.

8. the method for claim 7, which is characterized in that in the data processing step, form each detection positionDeformation-time graph, and record peak value and its corresponding time value in every deformation time graph.

9. method according to claim 8, which is characterized in that in the data processing step, according to the peak value and itsCorresponding time value and check bit glove reason distance value carry out straight line fitting, be fitted when it is m- apart from straight line, seek thisWhen the m- slope apart from straight line inverse to obtain shear wave velocity angle value, in conjunction with corresponding frequency, thus obtain shear wave velocityDegree-frequency straight line.

10. a kind of computer readable storage medium, which is characterized in that including program, described program can be executed by processor withRealize the method as described in any one of claim 4-9.