CN107505233A - Medium viscoplasticity quantitative approach and device - Google Patents

Abstract

The invention discloses a method and a device for quantifying the viscoelasticity of a medium. After the medium is excited by the vibration, a position-time map of the propagation of the vibration in the medium is acquired. The slope of the slope line of the position-time graph is obtained by image segmentation. A viscoelastic parameter of the medium is obtained according to the slope. The slope of the position-time diagram is the propagation velocity of the vibration in the medium. Since the propagation speed of vibration in the medium is related to the viscoelasticity of the medium, the viscoelastic parameters of the medium can be quantitatively calculated after the slope of the position-time diagram is obtained. The present invention does not need to select characteristic points from the position-time diagram to calculate the slope of the position-time diagram, is free from noise interference and has a small amount of calculation, and can efficiently and accurately quantify the viscoelasticity of the medium.

Description

Medium viscoelasticity quantitative method and device

technical field

The invention relates to the field of measurement technology, in particular to a method and device for quantifying medium viscoelasticity.

Background technique

When the medium is excited by vibration, the propagation characteristics of the vibration in the medium are related to the viscoelasticity of the medium. By measuring the propagation characteristics of the vibration, the viscoelasticity of the medium can be quantified.

The above principles have been applied to many technical fields. Taking medical detection as an example, when detecting organs or tissues such as liver, thyroid and muscle, the lesion can be located by quantifying the viscoelasticity of the medium.

Therefore, how to carry out efficient and accurate quantification of medium viscoelasticity is a problem that needs to be solved.

Contents of the invention

The embodiment of the present invention provides a method and a device for quantifying the viscoelasticity of a medium. In order to provide a basic understanding of some aspects of the disclosed embodiments, a brief summary is presented below. This summary is not an overview, nor is it intended to identify key/critical elements or delineate the scope of these embodiments. Its sole purpose is to present some concepts in a simplified form as a prelude to the more detailed description that is presented later.

According to the first aspect of the embodiments of the present invention, a method for quantifying the viscoelasticity of a medium is provided, the method comprising:

After the medium is excited by the vibration, obtaining a position-time diagram of the vibration propagating in the medium;

Obtaining the slope of the slope line of the position-time graph by image segmentation;

According to the slope, the viscoelastic parameters of the medium are obtained.

On the basis of the method, as an optional embodiment 1, the acquisition of the slope of the slope line of the position-time graph by using image segmentation includes:

performing image segmentation on the position-time map;

Extract image features;

Using the image features to perform linear fitting to obtain the slope of the slope line of the position-time graph.

On the basis of the first embodiment, as an optional second embodiment, the image feature is at least one of the following:

Central axis, peak value, valley value, zero crossing point, maximum value of the slope of the rising or falling part of the curve, minimum value of the slope of the rising or falling part of the curve, maximum value of the second derivative, and minimum value of the second derivative.

On the basis of the first embodiment, as an optional third embodiment, the method further includes:

Before the linear fitting, the extracted image features are screened to remove image features that deviate from the linear region by a set amount.

On the basis of the method, as an optional fourth embodiment, the acquiring a position-time diagram of vibration propagating in the medium includes:

performing frequency-domain transformation on the vibration detection signal to obtain a frequency-domain signal;

removing signals outside the set vibration speed range from the frequency domain signals to obtain processed signals;

A position-time map of the vibration is obtained using the processed signal.

On the basis of the fourth embodiment, as an optional embodiment five, the removal of signals outside the set vibration speed range in the frequency domain signal to obtain the processed signal includes:

performing filtering or eigenvalue selection on the frequency domain signal to obtain a processed signal;

The parameters of the filtering are related to the set vibration speed range;

The selection of the characteristic value is related to the set vibration speed range.

On the basis of any one of the method, the first embodiment to the fifth embodiment, as an optional embodiment six, the method further includes:

Before the image segmentation, the reflected waves in the position-time diagram are filtered out.

According to the second aspect of the embodiments of the present invention, there is provided a medium viscoelastic quantitative device, the device comprising:

A first acquisition module, configured to acquire a position-time diagram of vibration propagating in the medium after the medium is excited by vibration;

The second acquisition module is used to acquire the slope of the slope line of the position-time graph by using image segmentation;

The quantitative module is used to obtain the viscoelastic parameters of the medium according to the slope.

On the basis of the device, as an optional embodiment 1, the second acquisition module includes:

A segmentation sub-module, configured to perform image segmentation on the position-time map;

Extraction sub-module for extracting image features;

The fitting submodule is used to perform linear fitting by using the image features, and obtain the slope of the slope line of the position-time diagram.

On the basis of the first embodiment, as an optional second embodiment, the image feature is at least one of the following:

Central axis, peak value, valley value, zero crossing point, maximum value of the slope of the rising or falling part of the curve, minimum value of the slope of the rising or falling part of the curve, maximum value of the second derivative, and minimum value of the second derivative.

On the basis of the first embodiment, as an optional third embodiment, the fitting submodule screens the extracted image features before the linear fitting, and removes the deviation from the linear region to reach the set amount of image features.

On the basis of the device, as an optional fourth embodiment, the first acquisition module includes:

The first processing submodule is used to perform frequency domain transformation on the vibration detection signal to obtain a frequency domain signal;

The second processing submodule is used to remove signals outside the set vibration speed range in the frequency domain signal to obtain processed signals;

The obtaining submodule is used to obtain the position-time diagram of the vibration by using the processed signal.

On the basis of the fourth embodiment, as an optional fifth embodiment, the second processing submodule performs filtering or feature value selection on the frequency domain signal to obtain a processed signal;

The parameters of the filtering are related to the set vibration speed range;

The selection of the characteristic value is related to the set vibration speed range.

On the basis of the device, any one of the first embodiment to the fifth embodiment, as an optional fifth embodiment, the device further includes: a filtering module, configured to perform Before image segmentation, the reflected waves in the position-time diagram are filtered out.

According to a third aspect of an embodiment of the present invention, a medium viscoelastic quantitative device is provided, the device comprising:

memory, storing execution instructions;

The processor is configured to read the execution instruction and complete the following operations:

After the medium is excited by the vibration, obtaining a position-time diagram of the vibration propagating in the medium;

Obtaining the slope of the slope line of the position-time graph by image segmentation;

According to the slope, the viscoelastic parameters of the medium are obtained.

The technical solutions provided by the embodiments of the present invention may include the following beneficial effects:

The slope of the slope line of the position-time diagram is obtained by image segmentation, and the slope of the position-time diagram is the propagation speed of the vibration in the medium. Since the propagation speed of vibration in the medium is related to the viscoelasticity of the medium, the viscoelastic parameters of the medium can be quantitatively calculated after the slope of the position-time diagram is obtained. The embodiment of the present invention does not need to select feature points from the position-time diagram to calculate the slope of the position-time diagram, will not be disturbed by noise and the calculation amount is small, and can efficiently and accurately quantify the viscoelasticity of the medium.

It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the invention.

Description of drawings

The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the invention and together with the description serve to explain the principles of the invention.

Fig. 1 is a flowchart of a method for quantifying the viscoelasticity of a medium according to an exemplary embodiment;

Fig. 2 is a flowchart of a method for quantifying viscoelasticity of a medium according to an exemplary embodiment;

Fig. 3 is a flowchart of a method for quantifying viscoelasticity of a medium according to an exemplary embodiment;

Fig. 4 is a block diagram of a medium viscoelastic quantitative device according to an exemplary embodiment;

Fig. 5 is a block diagram of the second acquisition module shown in Fig. 4;

Fig. 6 is a block diagram of the first acquisition module shown in Fig. 4;

Fig. 7 is a block diagram of a medium viscoelastic quantitative device according to an exemplary embodiment;

Fig. 8 is a block diagram of a medium viscoelastic quantitative device according to an exemplary embodiment.

detailed description

The following description and drawings illustrate specific embodiments of the invention sufficiently to enable those skilled in the art to practice them. The examples merely represent possible variations. Individual components and functions are optional unless explicitly required, and the order of operations may vary. Portions and features of some embodiments may be included in or substituted for those of other embodiments. The scope of embodiments of the present invention includes the full scope of the claims, and all available equivalents of the claims. Herein, various embodiments may be referred to individually or collectively by the term "invention", which is for convenience only and is not intended to automatically limit the scope of this application if in fact more than one invention is disclosed. A single invention or inventive concept. Herein, relational terms such as first and second etc. are used only to distinguish one entity or operation from another without requiring or implying any actual relationship or relationship between these entities or operations. order. Furthermore, the term "comprises", "comprises" or any other variation thereof is intended to cover a non-exclusive inclusion such that a process, method or apparatus comprising a set of elements includes not only those elements but also other elements not expressly listed elements. Various embodiments herein are described in a progressive manner, each embodiment focuses on the differences from other embodiments, and the same and similar parts of the various embodiments may be referred to each other. As for the structures, products, etc. disclosed in the embodiments, since they correspond to the parts disclosed in the embodiments, the description is relatively simple, and for relevant parts, please refer to the description of the method part.

Fig. 1 is a flow chart of a method for quantifying viscoelasticity of a medium according to an exemplary embodiment. As shown in Figure 1, the method includes the following steps.

In step 11, after the medium is excited by the vibration, a position-time map of the vibration propagating in the medium is acquired.

In step 12, the slope of the slope line of the position-time diagram is obtained by image segmentation.

In step 13, according to the slope, the viscoelastic parameters of the medium are obtained.

The viscoelastic parameter may include at least one of a viscous parameter and an elastic parameter.

The slope of the position-time diagram is determined by the distance traveled by the vibration per unit time, that is, the speed at which the vibration travels in the medium. In a homogeneous medium, the speed of vibration propagation is related to the viscoelasticity of the medium. After obtaining the slope of the position-time diagram, the viscoelastic parameters of the medium can be calculated quantitatively. Therefore, how to obtain the above-mentioned slope efficiently and accurately becomes a quantitative measure of the viscoelasticity of the medium. key. In this exemplary embodiment, the region where vibration signals appear intensively is determined through image segmentation, and then the slope of the slope line of the position-time diagram is obtained. Compared with the traditional method of selecting feature points from the position-time diagram to calculate the above slope, the method in this exemplary embodiment is free from noise interference and has a small amount of calculation. It is an efficient and accurate method for quantifying the viscoelasticity of a medium.

In an exemplary embodiment, after the medium is excited to vibrate through mechanical vibration, acoustic radiation force or other means that can generate vibration, the medium vibrates, and the vibration propagates in the medium. Since the propagation speed of the above-mentioned vibrations in the medium is limited, the detection signal can be used to perform dynamic imaging on the medium. The above-mentioned detection signal may be light wave, ultrasonic wave or the like. The aforementioned dynamic imaging may be one-dimensional imaging, two-dimensional imaging, or three-dimensional imaging.

Since the above-mentioned vibration propagates in the medium, at different moments, the wave front will reach different positions along the propagation direction, and a position-time diagram can be formed according to this characteristic. The position-time plot is usually a straight, diagonally-striped plot. There are many ways to form a position-time map, and two examples are given below.

As an optional first implementation mode, the echo signal generated by the detection signal imaging the medium will undergo phase decorrelation. Using the characteristics of this phase decorrelation, the medium can be obtained through algorithms such as cross-correlation, autocorrelation, and optical flow. Motion information, along the set vibration propagation direction, to obtain a position-time diagram. The above-mentioned algorithm can be various block-matching-based methods or non-block-matching methods. In this embodiment, it is equivalent to first performing motion estimation on vibration propagation through information such as displacement and strain of the medium, and then obtaining the position-time diagram of the vibration.

As an optional second implementation manner, the presence or absence of vibration may also be determined by processing in the frequency domain, so as to obtain a position-time diagram. First, perform frequency domain transformation on the above detection signal to obtain the frequency domain signal, then remove the signal outside the set vibration speed range in the frequency domain signal to obtain the processed signal, and finally use the processed signal to obtain the vibration along the set vibration propagation direction position time map. In this implementation manner, no complicated calculation is required, and a position-time diagram that is not characterized by displacement or strain can be obtained by processing in the frequency domain. This method is equivalent to no need to estimate the motion of the vibration propagation first, and only needs to judge the presence or absence of the vibration to obtain the position-time map, which is a relatively efficient method for obtaining motion information.

Further, in the second implementation manner above, the frequency domain transformation can be implemented by selecting Fourier transform or singular value decomposition.

Further, in the above second implementation manner, filtering or feature value selection may be used to remove signals outside the set vibration speed range in the frequency domain signal. If the filtering method is used, the parameters of the filtering are related to the set vibration speed range, and if the eigenvalue selection method is used, the eigenvalue selection is related to the set vibration speed range.

Further, in the above-mentioned first and second embodiments, the above-mentioned set vibration propagation direction is the actual propagation direction of the vibration when the vibration propagates in only one propagation direction, and is the actual propagation direction of the vibration when the vibration propagates in multiple propagation directions. , it is a selected propagation direction. For example, when the medium is a uniform sheet, after the medium is vibrated, the vibration will propagate along the extension direction of the sheet, and the set vibration propagation direction at this time is the actual propagation direction of the vibration. For another example, when the medium has a three-dimensional irregular shape, the wavefront of vibration propagation is a three-dimensional shape, for example, the wavefront of vibration propagation is an ellipsoid, and the position-time diagrams obtained along different vibration propagation directions are different. A given vibration propagation direction is a selected propagation direction of interest. The propagation direction of interest is determined according to the actual direction to be measured, for example, it may be at least one of the fastest vibration propagation direction, the slowest vibration propagation direction, and the direction with a certain range of vibration propagation speeds.

In an exemplary embodiment, when the vibration propagates in the medium, reflected waves will be generated when encountering the edge of the medium or a foreign object, in order to improve the accuracy of subsequent processing, as shown in Figure 2, before image segmentation, you can also Step 11' is included, that is, filtering the reflected waves in the position-time diagram. There are many ways to filter out, and directional filtering is one of the implementation ways.

In an exemplary embodiment, when acquiring the slope of the slope line of the position-time diagram, image segmentation is performed first, and then image features are extracted. The image features here can be at least one of the central axis, peak value, valley value, and zero-crossing point, or other available feature points, such as the maximum and minimum values of the slope of the rising or falling part of the curve, and the maximum value of the second derivative. value and minimum value etc. The above-mentioned central axis refers to the skeleton of the pattern on the position-time diagram, and the above-mentioned zero-crossing point refers to the maximum point of the slope or the maximum point of the second-order derivative. The extracted image features all contain vibration information. Obtain the slope of the slope line of the position-time map by performing a linear fit on the extracted image features. This method of obtaining the slope of the slope line of the position-time graph avoids a large number of calculations, and it is a relatively efficient method to fit the slope of the slope line of the position-time graph through image feature points.

As an optional implementation manner, before performing linear fitting, the extracted image features may also be screened to remove image features that deviate greatly from the linear region, so as to further improve the accuracy of linear fitting. A set amount that deviates from the linear region can be set in advance, and when the removal operation is performed, image features reaching the set amount are removed.

As an optional implementation manner, image segmentation may be implemented using various methods, such as threshold-based segmentation methods, region-based segmentation methods, and edge-based segmentation methods.

As an optional implementation manner, the linear fitting can be implemented using the least square method.

According to the principle of mechanics, the viscoelasticity of the medium determines the propagation speed of the vibration in it. Therefore, by obtaining the slope of the position-time diagram, the velocity of the vibration propagating in the medium can be known, and then according to the principle of mechanics, the viscoelastic parameters of the medium can be quantitatively obtained . The viscoelastic parameters here may include shear modulus, viscoelastic modulus, shear viscoelasticity, viscous modulus, shear viscosity, mechanical impedance, mechanical relaxation time, anisotropy, and the like.

The application of the medium viscoelasticity quantification method in the embodiment of the present invention is given below in a specific application scenario.

When performing non-destructive viscoelasticity testing on viscoelastic media such as human liver, it is necessary to quantify the viscoelasticity of the media. The detection equipment includes an excitation device and an imaging device, wherein the excitation device performs vibration excitation on the medium to be detected, and the imaging device uses ultrasonic waves to image the vibration-excited medium. When the vibration propagates in the medium, at different moments, the wave front will reach different positions along the propagation direction, forming a position-time diagram. The aforementioned wave fronts may be one of crests, troughs, or the same phase of vibration.

As shown in FIG. 3 , the method for quantifying the viscoelasticity of media in this specific application scenario may include the following steps.

In step 31, the medium is vibrated.

In step 32, frequency-domain transform is performed on the ultrasonic signal of the detected vibration to obtain a frequency-domain signal.

In step 33, signals outside the set vibration speed range in the frequency domain signal are removed to obtain a processed signal.

In step 34, a position-time map of the vibration is obtained using the processed signal.

In step 35, direction filtering is performed on the position-time map.

In step 36, image segmentation is performed on the position-time map.

In step 37, the central axis is extracted.

In step 38, linear fitting is performed on the central axis to obtain the slope of the slope line of the position-time graph.

In step 39, according to the determined slope and the mechanical principle, the viscoelastic parameters of the medium are calculated.

In each exemplary embodiment of the above medium viscoelasticity quantitative method, when there are at least two set vibration propagation directions, each set vibration propagation direction will correspond to a position-time diagram, and then the position-time diagram will be obtained The graph corresponds to the viscoelastic parameters of the medium. At least two sets of viscoelasticity parameters obtained through synthesis can evaluate the viscoelasticity of the medium more comprehensively.

The various exemplary embodiments of the method for quantifying the viscoelasticity of the medium given above can be combined according to the situation, and the combination relationship between the various exemplary embodiments is not limited here.

Fig. 4 is a block diagram of a medium viscoelasticity quantitative device shown in an exemplary embodiment. The device can be located in the control host of the medium viscoelasticity testing equipment, for example, in the field of medical testing, it can be located in the control host of the liver non-destructive testing equipment. The device can also be located in the cloud, and the detection data of the medium viscoelasticity detection equipment needs to be processed in the cloud. The same concept in the previous method involved in the device described below will not be repeated here.

The device shown in FIG. 4 includes: a first acquisition module 41 , a second acquisition module 42 and a quantitative module 43 .

The first acquisition module 41 is configured to acquire a position-time diagram of vibration propagating in the medium after the medium is excited by vibration.

The second obtaining module 42 is configured to use image segmentation to obtain the slope of the slope line of the position-time graph.

The quantitative module 43 is used to obtain the viscoelastic parameters of the medium according to the slope.

In an exemplary embodiment, as shown in FIG. 5 , the second obtaining module 42 uses linear fitting to obtain the slope of the slope line of the position-time graph. At this point, the second acquisition module 42 includes: a segmentation submodule 421 , an extraction submodule 422 and a fitting submodule 423 .

The segmentation sub-module 421 is configured to perform image segmentation on the position-time map.

The extraction sub-module 422 is used to extract image features.

The fitting sub-module 423 is configured to perform linear fitting using image features to obtain the slope of the slope line of the position-time graph.

As an optional implementation, the fitting sub-module 423 may also screen the extracted image features before performing linear fitting, and remove image features that deviate from the linear region by a set amount.

In an exemplary embodiment, as shown in FIG. 6 , the first obtaining module 41 obtains a position-time map by processing the detection signal in the frequency domain. At this point, the first acquisition module 41 includes: a first processing submodule 411 , a second processing submodule 412 and an acquisition submodule 413 .

The first processing sub-module 411 is configured to perform frequency domain transformation on the vibration detection signal to obtain a frequency domain signal. The frequency domain transformation can adopt various time-frequency transformation methods such as Fourier transform and singular value decomposition.

The second processing sub-module 412 is configured to remove signals outside the set vibration speed range in the frequency domain signal to obtain processed signals.

The obtaining sub-module 413 is used to obtain the position-time map of the vibration by using the processed signal.

Further, the second processing sub-module 412 performs filtering or feature value selection on the frequency domain signal to obtain a processed signal. At this time, the parameters of the filter are related to the set vibration speed range, and the feature value selection is related to the set vibration speed range.

Of course, in other optional implementation manners, the first acquisition module 41 can also obtain the motion information of the medium through algorithms such as cross-correlation, auto-correlation, and optical flow, and obtain a position-time map along the set vibration propagation direction.

In an exemplary embodiment, as shown in FIG. 7 , the medium viscoelasticity quantification device further includes: a filter module 44 , configured to filter reflected waves in the position-time diagram before the second acquisition module 42 performs image segmentation.

Fig. 8 is a block diagram of a medium viscoelasticity quantitative device according to an exemplary embodiment. The device can be located in the control host of the medium viscoelasticity testing equipment, for example, in the field of medical testing, it can be located in the control host of the liver non-destructive testing equipment. The device can also be located in the cloud, and the detection data of the medium viscoelasticity detection equipment needs to be processed in the cloud.

The device shown in FIG. 8 includes: a memory 81 and a processor 82 .

The execution instruction is stored in the memory 81 .

The processor 82 is configured to read the execution instructions in the memory 81, and execute some or all of the steps in the various exemplary embodiments of the above-mentioned method for quantifying the viscoelasticity of the medium. Processor 82 may be implemented by a chip.

If the medium viscoelasticity quantitative device shown in Figure 8 is located in the control host of the medium viscoelasticity detection equipment, it can be connected with the excitation device and the imaging device in the medium viscoelasticity quantitative equipment through bus, wireless, etc. It has the interface and corresponding communication mechanism to realize the above connection.

If the medium viscoelasticity quantitative device shown in FIG. 8 is located in the cloud, it can communicate with the medium viscoelasticity detection equipment through the network.

It should be understood that the present invention is not limited to the processes and structures that have been described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from the scope thereof. The scope of the invention is limited only by the appended claims.

Claims (15)

1. a medium viscoelasticity quantitative method, is characterized in that, described method comprises: After the medium is excited by the vibration, obtaining a position-time diagram of the vibration propagating in the medium; Obtaining the slope of the slope line of the position-time graph by image segmentation; According to the slope, the viscoelastic parameters of the medium are obtained. 2. The method according to claim 1, wherein said utilizing image segmentation to obtain the slope of the slope line of said position-time diagram comprises: performing image segmentation on the position-time map; Extract image features; Using the image features to perform linear fitting to obtain the slope of the slope line of the position-time graph. 3. The method according to claim 2, wherein the image feature is at least one of the following: a central axis, a peak value, a valley value, a zero-crossing point, a maximum slope of a rising or falling part of a curve, a rising or falling curve The minimum value of the slope, the maximum value of the second derivative, and the minimum value of the second derivative of the section. 4. The method of claim 2, further comprising: Before the linear fitting, the extracted image features are screened to remove image features that deviate from the linear region by a set amount. 5. The method according to claim 1, wherein said obtaining a position-time diagram of vibration propagating in said medium comprises: performing frequency domain transformation on said vibration detection signal to obtain a frequency domain signal; removing signals outside the set vibration speed range from the frequency domain signals to obtain processed signals; A position-time map of the vibration is obtained using the processed signal. 6. The method according to claim 5, wherein said removing signals outside the set vibration velocity range in said frequency domain signal to obtain a processed signal comprises: performing filtering or eigenvalue selection on the frequency domain signal to obtain a processed signal; The parameters of the filtering are related to the set vibration speed range; The selection of the characteristic value is related to the set vibration speed range. 7. The method according to any one of claims 1 to 6, further comprising: before performing image segmentation, filtering reflected waves in the position-time diagram. 8. A medium viscoelastic quantitative device, characterized in that the device comprises: A first acquisition module, configured to acquire a position-time diagram of vibration propagating in the medium after the medium is excited by vibration; The second acquisition module is used to acquire the slope of the slope line of the position-time graph by using image segmentation; The quantitative module is used to obtain the viscoelastic parameters of the medium according to the slope. 9. The device according to claim 7, wherein the second acquiring module comprises: A segmentation sub-module, configured to perform image segmentation on the position-time map; Extraction sub-module for extracting image features; The fitting submodule is used to perform linear fitting by using the image features, and obtain the slope of the slope line of the position-time graph. 10. The device according to claim 9, wherein the image feature is at least one of the following: a central axis, a peak value, a valley value, a zero-crossing point, a maximum slope of a rising or falling part of a curve, a rising or falling curve The minimum value of the slope, the maximum value of the second derivative, and the minimum value of the second derivative of the section. 11. The device according to claim 9, wherein the fitting submodule filters the extracted image features before the linear fitting, and removes images that deviate from the linear region to a set amount feature. 12. The device according to claim 7, wherein the first acquiring module comprises: The first processing submodule is used to perform frequency domain transformation on the vibration detection signal to obtain a frequency domain signal; The second processing submodule is used to remove signals outside the set vibration speed range in the frequency domain signal to obtain processed signals; The obtaining submodule is used to obtain the position-time diagram of the vibration by using the processed signal. 13. The device according to claim 12, wherein the second processing submodule performs filtering or feature value selection on the frequency domain signal to obtain a processed signal; The parameters of the filtering are related to the set vibration speed range; The selection of the characteristic value is related to the set vibration speed range. 14. The device according to any one of claims 7 to 13, further comprising: a filtering module, configured to filter out the position-time map before the second acquisition module performs image segmentation reflected waves in . 15. A medium viscoelastic quantitative device, characterized in that the device comprises: memory, storing execution instructions; The processor is configured to read the execution instruction and complete the following operations: After the medium is excited by the vibration, obtaining a position-time diagram of the vibration propagating in the medium; Obtaining the slope of the slope line of the position-time graph by image segmentation; According to the slope, the viscoelastic parameters of the medium are obtained.