CN111401108B - Slope parameter extraction method and device and computer readable storage medium - Google Patents
Slope parameter extraction method and device and computer readable storage medium Download PDFInfo
- Publication number
- CN111401108B CN111401108B CN201910002361.2A CN201910002361A CN111401108B CN 111401108 B CN111401108 B CN 111401108B CN 201910002361 A CN201910002361 A CN 201910002361A CN 111401108 B CN111401108 B CN 111401108B
- Authority
- CN
- China
- Prior art keywords
- slope parameter
- target
- filtering
- coordinate
- propagation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
Abstract
The application provides a slope parameter extraction method, a device and a storage medium, which comprises the steps of obtaining a variable of a test signal which is transmitted along a preset direction along with time at a position to be tested; obtaining target coordinates corresponding to each propagation distance in the preset direction; and fitting the target coordinates corresponding to the propagation distances to obtain the slope parameters of the position to be measured. According to the scheme, the target coordinate is adopted to extract the slope parameter, so that the error caused by noise interference can be effectively weakened, and the accuracy and reliability of slope parameter extraction are improved.
Description
Technical Field
The present application relates to the field of data processing, and in particular, to a method and an apparatus for extracting a slope parameter, and a computer-readable storage medium.
Background
In applications of tracking wave propagation speed such as ultrasonic elastography, geological exploration, seismic wave detection and the like, wave speed detection is involved, namely, the propagation process of wave features along the interested propagation direction in a certain time is detected.
At present, a slope parameter is extracted by adopting an image processing method to obtain a propagation speed, and because the process of wave propagation under actual conditions is complex, the imaging quality of strain data is generally poor, and thus the slope parameter obtained by the existing method is often inaccurate.
Disclosure of Invention
The application provides a slope parameter extraction method, a slope parameter extraction device and a computer-readable storage medium, which are used for solving the problem that the slope parameter obtained by the existing method is inaccurate.
A first aspect of the present application provides a slope parameter extraction method, including: acquiring a variable of a test signal propagating along a preset direction at a position to be tested along with time; obtaining target coordinates corresponding to each propagation distance in the preset direction; and fitting the target coordinates corresponding to the propagation distances to obtain the slope parameters of the position to be measured.
Another aspect of the present application provides a slope parameter extraction apparatus, including: the acquisition module is used for acquiring a variable of the propagation of the test signal along a preset direction at the position to be tested along with time; the determining module is used for obtaining a target coordinate corresponding to each propagation distance in the preset direction; and the processing module is used for fitting and processing the target coordinates corresponding to the propagation distances to obtain the slope parameters of the position to be detected.
Another aspect of the present application provides a slope parameter extracting apparatus, including: at least one processor and memory; the memory stores computer-executable instructions; the at least one processor executes computer-executable instructions stored by the memory to perform the method as previously described.
It is a further aspect of the present application to provide a computer readable storage medium having stored therein a computer program which, when executed, implements a method as previously described.
According to the slope parameter extraction method, the slope parameter extraction device and the computer-readable storage medium, the target coordinates corresponding to each propagation distance in the preset direction are obtained by obtaining the variable of the propagation of the test signal in the position to be tested along the preset direction along with the time, and the target coordinates corresponding to each propagation distance are subjected to fitting processing, so that the slope parameter of the position to be tested is obtained. According to the scheme, the target coordinate is adopted to extract the slope parameter, so that the error caused by noise interference can be effectively weakened, and the accuracy and reliability of slope parameter extraction are improved.
Drawings
In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a schematic flow chart illustrating a slope parameter extraction method according to the present application in an exemplary embodiment;
FIG. 2a is a schematic flow chart illustrating a slope parameter extraction method according to the present application in accordance with another exemplary embodiment;
FIG. 2b is a schematic diagram of the signal propagation diagram in the embodiment of FIG. 2 a;
fig. 3 is a schematic structural diagram of a slope parameter extraction apparatus according to an exemplary embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are some embodiments of the present application, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
Fig. 1 is a schematic flow chart of a slope parameter extraction method according to the present application, shown in fig. 1, according to an exemplary embodiment, the slope parameter extraction method according to the present application includes:
102, obtaining a target coordinate corresponding to each propagation distance in the preset direction;
and 103, fitting and processing the target coordinates corresponding to each propagation distance to obtain a slope parameter of the position to be measured.
The test signal may have a variety of forms, for example, a sound wave, a shear wave, and the like, and the preset direction may be set in combination with an actual scene, for example, when the application scene is elasticity imaging, the test signal may be a shear wave, and the preset direction may be a depth direction. The propagation distance refers to a distance from a current arriving position to a propagation starting point when the test signal propagates along a preset direction. For example, when the test signal propagates in the depth direction, each propagation distance is each depth. The variable may be any variable reflecting the characteristics of the fluctuation, such as deformation data and the like. The destination coordinates are a series of coordinate points determined according to/based on variables and block matching methods. The slope parameter extracting method is applied to a slope parameter extracting device in various implementation forms, for example, the implementation form of the slope parameter extracting device may be a computer program stored in a storage medium, or software installed in an entity device, or an entity device, where the entity device includes but is not limited to an ultrasound imaging device, a computer, a terminal, a chip, and the like.
For example, the slope parameter extraction method of this embodiment may adopt an ultrasound scanning system to perform elastography on a tissue, and excite a test signal in the tissue, for example, excite a shear wave in the tissue by pressing, then acquire ultrasound echo data including propagation information of the shear wave propagating in the tissue in a depth direction, and obtain deformation estimation data using the ultrasound echo data, where the data includes deformation data of the shear wave at different depths at different times at a position to be measured. The deformation data includes variables reflecting propagation of the shear wave in the depth direction at the to-be-measured position of the to-be-measured tissue along with time, and the expression form of the deformation data may be a signal propagation diagram, for example, a two-dimensional data diagram including deformation information at different times and at different depths. The deformation data is in a stripe shape in the two-dimensional data map, the slope of the stripe is the slope parameter of the position to be measured and is also the propagation speed of the shear wave in the tissue, and then the selection is carried out by utilizing the formula of Young modulus, wherein E is 3 rho VS 2Wherein E is a hardness value, rho is the density of the tissue to be measured, and VSAnd (4) solving the elastic hardness value of the tissue for the speed value of the shear wave at the position to be measured, and further judging the pathology of the tissue. Wherein the variable canStrain values, displacement values, etc.
As an example, step 103 may specifically include: fitting the target coordinates corresponding to the propagation distances to obtain a fitted straight line; and taking the slope of the fitting straight line as a slope parameter of the position to be measured. The method for fitting the target coordinate may be various, such as fitting by using least square method, linear fitting, polynomial fitting, and the like.
Specifically, as described above, the variable corresponding to each propagation distance may be represented by a stripe in a signal propagation map, and in order to obtain a more accurate slope parameter, the selection of the stripe is particularly important, in this scheme, in order to select a stripe with better quality, a target coordinate corresponding to each propagation distance needs to be determined first, where the target coordinate is a coordinate point determined by a preset algorithm, for example, a coordinate point on a stripe corresponding to a different propagation distance (for example, tissue depth) obtained by a normalized cross-correlation algorithm based on a block matching method, and then, a slope parameter of the position to be measured is obtained by fitting a series of target coordinates corresponding to different propagation distances. The scheme of solving the tissue slope by the target coordinate can effectively weaken the error caused by noise interference, so that the result is more accurate.
In this embodiment, the propagation distances may be selected from the signal propagation maps mentioned below, and the propagation distances corresponding to the middle segments with less interference are still selected in combination with the ultrasound imaging scene, for example, the depth of the fat layer near the skin (e.g. the depth is 0cm to 1cm), the depth with the deepest part affected by noise is removed, and the like, that is, the depths at both ends are usually removed, and the depths corresponding to the middle segments are selected. The size of the selected interval between the propagation distances may be determined according to the size of the ultrasonic frequency, for example, if the ultrasonic frequency is larger, the interval between the propagation distances is smaller; if the ultrasonic frequency is small, the interval between the propagation distances is large.
In the slope parameter extraction method of this embodiment, the target coordinates corresponding to each propagation distance in the preset direction are obtained by obtaining the variable of the propagation of the test signal in the position to be tested along the preset direction with time, and the target coordinates corresponding to each propagation distance are subjected to fitting processing, so as to obtain the slope parameter of the position to be tested. According to the scheme, the target coordinate is adopted to extract the slope parameter, so that the error caused by noise interference can be effectively weakened, and the accuracy and reliability of slope parameter extraction are improved.
Fig. 2a is a schematic flow chart of a slope parameter extraction method according to the present application, shown in another exemplary embodiment, and as shown in fig. 2a, on the basis of the previous embodiment, step 102 may specifically include:
and 203, determining a destination coordinate corresponding to each propagation distance in the signal propagation diagram based on a block matching method and the current destination coordinate.
Taking an elastography scene as an example, as described in the previous embodiment, the signal propagation map (as shown in fig. 2 b) is a two-dimensional data map in which the horizontal axis is a time axis (e.g., ms) and the vertical axis is a depth axis (e.g., mm) representing a propagation distance, the two-dimensional data map reflects propagation information of shear waves at different depths at different times at the position to be measured, the horizontal coordinate of each coordinate point in the two-dimensional data map is a time coordinate, and the vertical coordinate is a depth coordinate representing the propagation distance. The two-dimensional data map may be presented as a bar graph, and the color of each bar in the bar graph may represent the value of a certain depth variable at a certain time. Of course, the signal propagation diagram can also be represented by other forms of graphs, such as a three-dimensional graph. This embodiment is not particularly limited thereto. Wherein, the dotted line in fig. 2b is the tissue slope of the position to be measured obtained after the fitting process.
Still exemplified in connection with elastography scenarios: the range and position of the initial region to be measured can be specified in the signal propagation diagram according to requirements. The initial destination coordinate may be a maximum coordinate, a minimum coordinate, or a median coordinate determined based on the variable. When the initial target coordinate is the maximum coordinate, the initial target coordinate is the coordinate of the maximum variable value in the initial region to be measured; when the initial target coordinate is the minimum coordinate, the initial target coordinate is the coordinate of the minimum variable value in the initial region to be measured; and when the initial target coordinate is the median coordinate, the initial target coordinate is a coordinate of a middle value determined according to the maximum value coordinate and the minimum value coordinate corresponding to the same propagation distance in the initial region to be measured. Specifically, the intermediate value coordinate is determined according to the intermediate value of the abscissa of the maximum value coordinate and the minimum value coordinate.
Alternatively, when the initial target coordinate is the median coordinate, the initial target coordinate may be determined based on the maximum coordinate, or may be determined based on the minimum coordinate. Specifically, when the maximum value coordinate is taken as a reference, firstly, the maximum value coordinate is determined from the initial region to be measured, the propagation distance of the maximum value coordinate is obtained, then, the coordinate of the minimum variable value corresponding to the propagation distance is determined, so that the minimum value coordinate corresponding to the propagation distance is determined, and then, the intermediate value of the horizontal and vertical coordinates of the maximum value coordinate and the minimum value coordinate is calculated, so that the intermediate value coordinate is obtained; when the minimum value coordinate is taken as a reference, the minimum value coordinate is determined from the initial region to be measured, the propagation distance of the minimum value coordinate is obtained, the coordinate of the maximum variable value corresponding to the propagation distance is determined, the maximum value coordinate is obtained, and then the intermediate value of the horizontal and vertical coordinates of the maximum value coordinate and the minimum value coordinate is calculated, so that the intermediate value coordinate is obtained.
Optionally, as an implementation manner, step 203 of this embodiment may specifically include:
taking the current target coordinate as a starting point, and acquiring M coordinates along the time direction to obtain a group of reference data points;
taking the corresponding position of the current target coordinate in the next propagation distance as a next starting point, moving the next starting point for N times by a preset step length, and acquiring M coordinates along the time direction every time of moving to obtain N groups of data points to be determined, wherein each group of data points to be determined comprises the M coordinates;
based on the block matching method, calculating matching coefficients between the group of reference data points and the N groups of data points to be determined respectively to obtain N matching coefficients;
determining a target coordinate corresponding to the next propagation distance according to the N matching coefficients;
and updating the current target coordinate to be a target coordinate corresponding to the next propagation distance, and repeatedly executing the step of taking the current target coordinate as a starting point and obtaining M coordinates along the time direction to obtain a group of reference data points until the target coordinate corresponding to each propagation distance in the signal propagation diagram is found.
The block matching method may be a cross-correlation algorithm, and accordingly, the obtained matching coefficient is the cross-correlation coefficient.
Optionally, the destination coordinate corresponding to the next propagation distance is determined according to the N matching coefficients, which may be specifically implemented by determining a matching coefficient representing the best matching degree among the N matching coefficients, and taking the start coordinate in the group of to-be-determined data points corresponding to the matching coefficient representing the best matching degree as the destination coordinate corresponding to the next propagation distance.
Further optionally, the optimal matching degree may be set according to requirements, for example, the correlation coefficient is maximum, the variance value is minimum, and the like.
For example, taking the horizontal axis of the signal propagation diagram as the time axis, the vertical axis as the depth axis representing the propagation distance, and the block matching method as the cross-correlation algorithm for example, a preset region in the diagram is used as an initial region to be measured, and an adaptive technique may be used to find an initial target coordinate in the region to be measured, where the adaptive technique refers to a processing method that is automatically adjusted according to the data characteristics of the processed data during the processing and analysis process, and is used here to adaptively find a first target coordinate point (i.e., an initial target coordinate), where the first target coordinate may be a starting point of a stripe, or may be any coordinate point in the stripe, for example, a 2 nd data point is used as a target coordinate, or a 5 th data point is used as a target coordinate. Assuming that the first target coordinate is a starting point of the stripe, a certain amount of data is selected with the first target coordinate as a starting position, for example, M data points are selected, the M data points are used as reference data points, and then the position of a data point corresponding to the starting point of the depth position (i.e., propagation distance) of the first target coordinate in the next designated depth is used as a starting point, the starting point is moved in a certain search range on a time axis to obtain M data points, wherein the search range can be moved by a preset step size of 1, that is, two groups of data points obtained after each movement are displaced by one data bit position from each other. If the shift is performed N times, N sets of M data points (i.e., data points to be determined) can be obtained, i.e., a total of Nx M data points to be determined is obtained. The values of the N normalized cross-correlation coefficients can be obtained by calculating the normalized cross-correlation coefficients between the M reference data points and the N groups of M data points to be determined, and the maximum value is selected from the N values, so that the starting point of the group of M data points corresponding to the maximum value is the target coordinate corresponding to the specified depth. That is, the point most correlated with the target coordinate of the streak at the previous depth is found in the next specified depth, and a straight line is fitted and the slope parameter is obtained from the obtained target coordinates of the streaks corresponding to the specified depths. Preferably, the precision of the cross-correlation coefficient can be improved by adopting a cosine difference method. And then the next depth data repeats the above process with the corresponding position of the destination coordinate determined just now as the starting point until all the destination coordinates of the specified depth are obtained. By adopting the normalized cross-correlation technique to calculate the organization slope, when the abscissa representing the time change is calculated, because more data are adopted in the normalized cross-correlation calculation (for example, a group of M reference data points and a group of M data points to be determined are needed), and the data are not dependent on one point or a plurality of points, the interference caused by noise in clinic can be effectively weakened, meanwhile, more data participate in the calculation, and the result is more accurate.
Optionally, in order to further improve the accuracy of the result, the target coordinate may be further filtered. Optionally, on the basis of any embodiment, the method may further include:
and screening the target coordinates corresponding to the propagation distances according to a preset screening scheme.
In practical application, various screening schemes can be adopted to optimize the extraction result. The screening scheme is described below by way of example in several embodiments:
in an embodiment, the screening the destination coordinates corresponding to each propagation distance according to a preset screening scheme may specifically include: sequentially fitting P adjacent target coordinates in the target coordinates corresponding to each propagation distance to obtain a corresponding fitting straight line set, wherein P is a positive integer and is smaller than the number of the target coordinates corresponding to each propagation distance; determining a target fitting straight line in the fitting straight line set, wherein the sum of the residual errors from the target coordinates corresponding to each propagation distance to the target fitting straight line is minimum; and removing the target coordinates of which the distance to the target fitting straight line exceeds a preset distance threshold from the target coordinates corresponding to each propagation distance except the P adjacent target coordinates for fitting the target fitting straight line.
Specifically, in the screening scheme, target coordinates with large differences are eliminated, for example, in an actual scene, assuming that 10 target coordinates are obtained altogether, and P is 5, then adjacent 5 target coordinate points (first to fifth target coordinates, second to sixth target coordinates, …, and sixth to tenth target coordinates) are sequentially selected to perform linear fitting based on least square, so as to obtain a corresponding fitting straight line set, a straight line is selected as a target fitting straight line, and the target fitting straight line is a residual error from all target coordinate points (10 target coordinate points) to the fitting straight line and a minimum fitting straight line. The method comprises the steps of solving the distances from all other target coordinate points (the other 5 points) to a target fitting straight line except for 5 points for fitting the target straight line, eliminating the target coordinate points with the distances larger than a preset distance threshold value to finally obtain a series of target coordinate values, namely eliminating the target coordinate points with larger differences from 10 target coordinates, and then fitting the straight line according to the remaining series of target coordinates which are smaller than or equal to 10 target coordinate points to obtain the slope parameter of the position to be measured.
In another embodiment, the screening the destination coordinates corresponding to each propagation distance according to a preset screening scheme may specifically include: calculating the difference value between the X-th target coordinate and the X-1 th target coordinate of the abscissa in the target coordinate corresponding to each propagation distance in turn, wherein X is a positive integer which is more than 2 and less than the number of each propagation distance; if the difference value between the X-th target coordinate and the X-1 th target coordinate is larger than a preset first threshold value, subtracting the corresponding 1 st to X-2 th target coordinate from the X-2 nd to X-1 th target coordinate, and if the average value of the absolute values of the obtained difference values is smaller than a preset second threshold value and X is smaller than a first preset clearing quantity threshold value, clearing the target coordinate before the X-th target coordinate; and if the X is larger than a second preset clearing number threshold value, clearing the target coordinate after the X-th target coordinate.
Specifically, taking an actual scene as an example: after each destination coordinate is obtained, starting from the third destination coordinate, when the abscissa of the X-th destination coordinate minus the abscissa of the X-th destination coordinate is greater than the first threshold (e.g., 5.5), and the abscissa of the 2 … X-1 destination coordinate minus the abscissa of the corresponding 1 … X-2 destination coordinate, the mean of the absolute values of all differences is less than the second threshold (e.g., 4), and X is less than the first preset clearing number threshold, e.g., when the current propagation distance is less than the specified propagation distance (e.g., row 25), then discarding all destination coordinate values of the previous row; similarly, if X is greater than a second predetermined number of erasures threshold, e.g., the current propagation distance is greater than the specified propagation distance (e.g., line 20), then all destination coordinate values of the following lines are discarded.
In another embodiment, the screening the destination coordinates corresponding to each propagation distance according to a preset screening scheme may specifically include: if the difference value between the abscissa of the most-valued coordinates corresponding to the K continuous adjacent propagation distances is smaller than a preset third threshold value, only the K' th most-valued coordinate is reserved; the relationship between K 'and K is subject to practical experience, for example, when K is 2K, K' ═ K; when K is 2K +1, K' is K + 1; wherein k is a positive integer.
Specifically, the target coordinates are traversed, and if the difference value between the abscissa of every two adjacent target coordinates in the K adjacent target coordinates is smaller than a preset third threshold, the K target coordinates are screened. For example, when the abscissa values of consecutive adjacent K (e.g. 4) coordinates differ by less than a third threshold (e.g. 1), and since K ═ 4 is an even number of destination coordinates, only the 2 nd destination coordinate of the 4 destination coordinates is retained; when the difference between the abscissa values of the consecutive adjacent K (e.g. 5) coordinates is smaller than a third threshold (e.g. 1), and K ═ 5 is an odd number of destination coordinates, only the 3 rd destination coordinate of the 5 destination coordinates is reserved; that is, only the coordinates of the intermediate point (K +1) are selected when the number K is an odd number (K2K +1), and only the coordinates of the K-th point are selected when the number K is an even number (K2K).
It is understood that, the order of the above screening schemes executed sequentially may include, but is not limited to: the execution is performed simultaneously or first, and this embodiment is only an exemplary embodiment and does not limit the execution order. The screening protocols can be carried out either individually or in combination.
According to the slope parameter extraction method provided by the embodiment, after the target coordinates corresponding to different propagation distances are obtained, filtering and screening are performed on each target coordinate, and then the screened target coordinates are fitted to obtain the slope parameter value of the tissue, so that the target coordinates with larger errors can be removed, and the accuracy and reliability of slope parameter extraction are further improved.
Optionally, in order to further improve the accuracy of the slope parameter and avoid the randomness and inaccuracy of a single measurement result, in practical application, the accuracy of the measurement result may be improved through multiple measurements. Accordingly, step 102 may be performed multiple times to obtain multiple sets of destination coordinates, based on any of the embodiments; correspondingly, step 103 may specifically include:
respectively fitting according to the multiple groups of target coordinates to obtain multiple slope parameters;
and calculating the median of the slope parameters to obtain the slope parameters of the position to be measured.
For example, the step 102 is repeatedly executed T times, that is, the step of obtaining the target coordinates corresponding to each propagation distance in the preset direction is repeatedly executed T times, and each target coordinate measured T times is obtained; wherein T is a positive integer. For example, T is taken to be 10, i.e., 10 measurements are taken. A set of destination coordinates can be obtained for each measurement and finally 10 sets of destination coordinates can be obtained. Correspondingly, the fitting process of the target coordinate is also executed for T times, namely the target coordinate corresponding to each propagation distance is respectively fitted and processed according to each target coordinate corresponding to the T times of measurement, and the slope parameter set { A ] of the position to be measured is obtained1,..,Ai,..,ATIn which the variable i characterizes the ith measurement, AiA slope parameter obtained according to each target coordinate measured at the ith time; therefore, a slope parameter is determined from a set consisting of T target coordinates obtained after T times of fitting and used for calculating the hardness value of the tissue. Alternatively, one could set the slope parameter { A }1,..,Ai,..,ATDetermining a median value, and determining a slope parameter of the position to be detected according to the median value.
For example, after obtaining the target coordinates of T measurements, taking the previously obtained 10 sets of target coordinates as an example, a set { a ] composed of slope parameters of 10 positions to be measured can be obtained by linear fitting according to the 10 sets of target coordinates1,..,Ai,..,A10Correspondingly, the median in the set can be used as the final slope parameter of the position to be measured.
Optionally, in order to further improve the accuracy of the slope parameter, the obtained variable of the test signal propagating along the preset direction at the position to be measured along with time may be filtered, and then the respective slope parameters of the filtered variable and the filtered variable obtained by T times of measurement are respectively obtained (i.e. 2T slope parameters may be obtained), which may specifically be:
performing band-pass filtering on the variable;
before and after filtering, repeating the step 102 for T times respectively to obtain target coordinates of each propagation distance before and after filtering;
correspondingly, the fitting process of the target coordinate corresponding to each propagation distance to obtain the slope parameter of the position to be measured includes:
according to the target coordinates corresponding to the T times of measurement before filtering, respectively fitting and processing the target coordinates corresponding to the propagation distances to obtain a pre-filtering slope set { B ] of the position to be measured1,..,Bi,..,BT}；
According to the target coordinates corresponding to the T times of measurement after filtering, respectively fitting and processing the target coordinates corresponding to each propagation distance to obtain a filtered slope set { C ] of the position to be measured1,..,Ci,..,CT}；
Set of pre-filtering slopes at position under test { B1,..,Bi,..,BTDetermining a median value as a slope parameter before filtering; set of filtered slopes at position under test { C1,..,Ci,..,CTAnd determining a median value as a slope parameter after filtering.
Optionally, the method of performing the band-pass filtering may be a least-square based band-pass filtering. Optionally, after obtaining the pre-filtering slope parameter and the post-filtering slope parameter, the method may further include: taking the average value of the pre-filtering slope parameter and the post-filtering slope parameter as the slope parameter of the position to be detected; or selecting a target slope from the pre-filtering slope parameter and the post-filtering slope parameter as a slope parameter of the position to be measured.
Optionally, in a tissue elastography scene, after the pre-filtering slope parameter and the post-filtering slope parameter are obtained, a more accurate result can be further selected by referring to the positions of the respective corresponding stripes. Specifically, the slope parameter of the position to be measured is determined according to one or more of the following selection conditions in the pre-filtering slope parameter and the post-filtering slope parameter:
if the tissue elastic hardness value corresponding to the before-filtering slope parameter and the tissue elastic hardness value corresponding to the after-filtering slope parameter are in the same pathological stage, determining a slope parameter corresponding to the tissue elastic hardness value with smaller relative deviation as a slope parameter of the position to be detected; wherein the relative deviation may be a ratio of the quartile to the median of the slope set.
If the tissue elasticity hardness value corresponding to the slope parameter before filtering and the tissue elasticity hardness value corresponding to the slope parameter after filtering are not in the same pathological stage, judging whether the relative deviation of the tissue elasticity hardness values corresponding to the slope parameter after filtering is smaller than a fourth threshold value or not, and if the relative deviation is smaller than the fourth threshold value, determining the slope parameter after filtering as the slope parameter of the position to be detected; and if not, judging whether the relative deviation of the tissue elastic hardness value corresponding to the slope parameter before filtering is smaller than a fifth threshold, and if so, determining the slope parameter before filtering as the slope parameter of the position to be detected.
If the tissue elastic hardness value corresponding to the pre-filtering slope parameter and the tissue elastic hardness value corresponding to the post-filtering slope parameter do not meet the selection condition, discarding each target coordinate corresponding to the current T times of measurement, and re-executing the step 101.
The fourth threshold and the fifth threshold may be equal or unequal, and preferably, the values of the fourth threshold and the fifth threshold may be 0.33; the histopathological stages correspond to different hardness value ranges, and the actual application can refer to the tissue hardness value and a control reference standard of the histopathological stages. And determining whether the slope parameter before filtering and the slope parameter after filtering correspond to the tissue elastic hardness value in the same pathological stage, and determining whether the hardness value obtained by calculating the slope parameter before filtering is used as a final result or the hardness value obtained by calculating the slope parameter after filtering is used as the final result according to various selection conditions. If not, the error of the current T times of measurement is considered to be larger, the target coordinates of the T times of measurement are abandoned, the T times of measurement is carried out again, and the target data of each T group before filtering and after filtering are obtained.
According to the slope parameter extraction method provided by the embodiment, the final slope parameter is obtained based on the multiple measurement results performed before and after filtering through the multiple measurement, filtering and other processing, and the accuracy and reliability of slope parameter extraction are further improved.
Fig. 3 is a schematic structural diagram of a slope parameter extracting apparatus according to the present application, shown in fig. 3, the slope parameter extracting apparatus according to the present application includes:
the acquiring module 31 is configured to acquire a variable of a test signal propagating along a preset direction at a to-be-tested position along with time;
a determining module 32, configured to obtain a destination coordinate corresponding to each propagation distance in the preset direction;
and the processing module 33 is configured to fit and process the destination coordinates corresponding to the propagation distances to obtain a slope parameter of the position to be measured.
The test signal may be a shear wave signal, the direction may be set in combination with an actual scene, and the variable may be any variable reflecting a fluctuation characteristic, such as deformation data. The slope parameter extracting apparatus of this embodiment may be implemented in various forms, for example, as a computer program stored in a storage medium, or as software installed in a physical device, or as a physical device, where the physical device includes, but is not limited to, an ultrasound imaging device, a computer, a terminal, a chip, and the like.
For example, the slope parameter extraction method of this embodiment may use an ultrasound scanning system to perform ultrasound elastography on a tissue, excite a shear wave in the tissue, for example, excite the shear wave in the tissue by pressing, then acquire ultrasound echo data including propagation information of the shear wave propagating in the tissue along a depth, and obtain deformation estimation data using the ultrasound echo data, where the data includes strain data of the shear wave at different depths at different times at a position to be measured. The strain data includes variables reflecting propagation of shear waves in the depth direction at the to-be-measured position of the to-be-measured tissue along with time, and the expression form of the variables may be a signal propagation map, for example, a two-dimensional data map including deformation information at different times and at different depths. The strain data is in a stripe shape in the two-dimensional data map, the slope of the stripe is the slope parameter of the position to be measured and is also the propagation speed of the shear wave in the tissue, and then the selection is carried out by utilizing the formula of Young modulus, wherein E is 3 rho VS 2Wherein E is a hardness value, rho is the density of the tissue to be measured, and VSAnd (4) solving the elastic hardness value of the tissue for the speed value of the shear wave at the position to be measured, and further judging the pathology of the tissue.
The fitting method may be various. As an example, the processing module may be specifically configured to: performing linear fitting on the screened target coordinate based on a least square method to obtain a second fitting straight line; and taking the slope of the second fitting straight line as a slope parameter of the position to be measured.
The target coordinate is a coordinate point determined by a preset algorithm, for example, a coordinate point on a stripe corresponding to different propagation distances (e.g., tissue depths) obtained by a normalized cross-correlation algorithm based on a block matching method. The scheme of solving the tissue slope by the target coordinate can effectively weaken the error caused by noise interference, so that the result is more accurate.
The slope parameter extraction device of this embodiment obtains the target coordinates corresponding to each propagation distance in the preset direction by obtaining the variable of the propagation of the test signal in the position to be measured along the preset direction with time, and performs fitting processing on the target coordinates corresponding to each propagation distance, thereby obtaining the slope parameter of the position to be measured. According to the scheme, the target coordinate is adopted to extract the slope parameter, so that the error caused by noise interference can be effectively weakened, and the accuracy and reliability of slope parameter extraction are improved.
Another exemplary embodiment shows a schematic structural diagram of a slope parameter extraction apparatus according to the present application, and on the basis of the above embodiment, the determining module 32 includes:
the generating unit is used for generating a signal propagation diagram after acquiring a variable of the test signal propagating along a preset direction along with time at a position to be tested;
the determining unit is used for determining an initial target coordinate from the initial region to be measured of the signal propagation diagram and taking the initial target coordinate as the current target coordinate;
and the searching unit is used for determining the target coordinates corresponding to the propagation distances in the signal propagation diagram based on a block matching method and the current target coordinates.
Optionally, as an implementation manner, the searching unit may be specifically configured to obtain M coordinates along a time direction with the current destination coordinate as a starting point, so as to obtain a group of reference data points; the searching unit is further specifically configured to move the next starting point N times with a preset step length by using a corresponding position of the current destination coordinate in a next propagation distance as a next starting point, and acquire M coordinates along a time direction every time the next starting point is moved, so as to obtain N groups of data points to be determined, where each group of data points to be determined includes the M coordinates; the searching unit is further specifically configured to calculate matching coefficients between the group of reference data points and the N groups of data points to be determined, respectively, based on the block matching method, to obtain N matching coefficients; the searching unit is further specifically configured to determine a destination coordinate corresponding to the next propagation distance according to the N matching coefficients; the searching unit is further specifically configured to update the current destination coordinate to be a destination coordinate corresponding to the next propagation distance, and repeatedly perform the step of obtaining M coordinates along the time direction with the current destination coordinate as a starting point to obtain a group of reference data points until the destination coordinate corresponding to each propagation distance in the signal propagation map is searched.
Optionally, the searching unit is specifically configured to determine a matching coefficient representing the best matching degree from the N matching coefficients, and use an initial coordinate in the set of data points to be determined corresponding to the matching coefficient as a destination coordinate corresponding to the next propagation distance.
Optionally, in order to further improve the accuracy of the result, the target coordinate may be further filtered. Optionally, on the basis of any embodiment, the apparatus further includes:
and the screening module is used for screening the target coordinates corresponding to the propagation distances according to a preset screening scheme.
In one embodiment, the screening module is specifically configured to: sequentially carrying out linear fitting based on a least square method on P adjacent target coordinates in the target coordinates corresponding to each propagation distance to obtain a corresponding fitting straight line set, wherein P is a positive integer and is smaller than the number of the target coordinates corresponding to each propagation distance; selecting a first fitting straight line in the fitting straight line set as a target fitting straight line, wherein the sum of the residual errors from the target coordinates corresponding to each propagation distance to the first fitting straight line is minimum; and removing the target coordinates of which the distance to the target fitting straight line exceeds a preset distance threshold from the target coordinates corresponding to each propagation distance except the P adjacent target coordinates for fitting the target fitting straight line.
In another embodiment, the screening module is specifically configured to: calculating the difference value between the X-th target coordinate and the X-1 th target coordinate of the abscissa in the target coordinate corresponding to each propagation distance in turn, wherein X is a positive integer which is more than 2 and less than the number of each propagation distance; if the difference value between the X-th target coordinate and the X-1 th target coordinate is larger than a preset first threshold value, subtracting the corresponding 1 st to X-2 th target coordinate from the X-2 nd to X-1 th target coordinate, and if the average value of the absolute values of the obtained difference values is smaller than a preset second threshold value and X is smaller than a first preset clearing quantity threshold value, clearing the target coordinate before the X-th target coordinate; and if the X is larger than a second preset clearing number threshold value, clearing the target coordinate after the X-th target coordinate.
In another embodiment, the screening module is specifically configured to: if the difference value between the abscissa of the most-valued coordinates corresponding to the K continuous adjacent propagation distances is smaller than a preset third threshold value, only the K' th most-valued coordinate is reserved; the relationship between K 'and K is subject to practical experience, for example, when K is 2K, K' ═ K; when K is 2K +1, K' is K + 1; wherein k is a positive integer.
It is understood that, the order of the above screening schemes executed sequentially may include, but is not limited to: the execution is performed simultaneously or first, and this embodiment is only an exemplary embodiment and does not limit the execution order. The screening protocols can be carried out either individually or in combination.
The slope parameter extraction device provided by this embodiment filters and screens each target coordinate after obtaining target coordinates corresponding to different propagation distances, and then obtains a slope parameter value of a tissue by adopting a least square line fitting, so that a target coordinate with a large error can be removed, and the accuracy and reliability of slope parameter extraction are further improved.
Optionally, in order to further improve the accuracy of the slope parameter, on the basis of any embodiment, the step of obtaining the destination coordinates corresponding to each propagation distance in the preset direction is performed by the determining module for multiple times to obtain multiple sets of destination coordinates; correspondingly, the processing module comprises: the fitting unit is used for respectively fitting to obtain a plurality of slope parameters according to the plurality of sets of target coordinates; and the calculating unit is used for calculating the median of the slope parameters to obtain the slope parameters of the position to be measured.
Optionally, in order to further improve the accuracy of the slope parameter, the apparatus may further include: and the filtering module is used for carrying out band-pass filtering on the variable.
Correspondingly, the determining module 32 respectively repeats the step 102T times before filtering and after filtering, and obtains the target coordinates of each propagation distance before filtering and after filtering;
the processing module 33 is specifically configured to: according to the target coordinates corresponding to the T times of measurement before filtering, respectively fitting and processing the target coordinates corresponding to the propagation distances to obtain a pre-filtering slope set { B ] of the position to be measured1,..,Bi,..,BT}; according to the target coordinates corresponding to the T times of measurement after filtering, respectively fitting and processing the target coordinates corresponding to each propagation distance to obtain a filtered slope set { C ] of the position to be measured1,..,Ci,..,CT}; set of pre-filtering slopes at position under test { B1,..,Bi,..,BTDetermining a median value as a slope parameter before filtering; set of filtered slopes at position under test { C1,..,Ci,..,CTAnd determining a median value as a slope parameter after filtering.
Optionally, after obtaining the pre-filtering slope parameter and the post-filtering slope parameter, the processing module may be further configured to use an average value of the pre-filtering slope parameter and the post-filtering slope parameter as the slope parameter of the position to be measured; or selecting a target slope from the pre-filtering slope parameter and the post-filtering slope parameter as a slope parameter of the position to be measured.
Optionally, the processing module may be further configured to determine, according to one or more of the following selection conditions, a slope parameter of the to-be-measured position in the pre-filtering slope parameter and the post-filtering slope parameter:
if the tissue elastic hardness value corresponding to the before-filtering slope parameter and the tissue elastic hardness value corresponding to the after-filtering slope parameter are in the same pathological stage, determining a slope parameter corresponding to the tissue elastic hardness value with smaller relative deviation as a slope parameter of the position to be detected; wherein the relative deviation may be a ratio of the quartile to the median of the slope set.
If the tissue elasticity hardness value corresponding to the slope parameter before filtering and the tissue elasticity hardness value corresponding to the slope parameter after filtering are not in the same pathological stage, judging whether the relative deviation of the tissue elasticity hardness values corresponding to the slope parameter after filtering is smaller than a fourth threshold value or not, and if the relative deviation is smaller than the fourth threshold value, determining the slope parameter after filtering as the slope parameter of the position to be detected; and if not, judging whether the relative deviation of the tissue elastic hardness value corresponding to the slope parameter before filtering is smaller than a fifth threshold, and if so, determining the slope parameter before filtering as the slope parameter of the position to be detected.
Optionally, the processing module is further configured to discard target coordinates corresponding to the current T times of measurement and instruct the obtaining module to perform the step of 101 again if the tissue elastic hardness value corresponding to the pre-filtering slope parameter and the tissue elastic hardness value corresponding to the post-filtering slope parameter do not satisfy the selection condition.
The slope parameter extraction device provided by this embodiment obtains the final slope parameter based on the multiple measurement results performed before and after filtering through multiple measurements, filtering, and the like, and further improves the accuracy and reliability of slope parameter extraction.
The embodiment of the present invention further provides a slope parameter extraction device, which is characterized by including: at least one processor and memory; the memory stores computer-executable instructions; the at least one processor executes computer-executable instructions stored by the memory to perform the method as previously described.
An embodiment of the present invention further provides a computer-readable storage medium, where a computer program is stored in the computer-readable storage medium, and when the computer program is executed, the computer program implements the method described above.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working process of the apparatus described above may refer to the corresponding process in the foregoing method embodiment, and is not described herein again.
Those of ordinary skill in the art will understand that: all or a portion of the steps of implementing the above-described method embodiments may be performed by hardware associated with program instructions. The program may be stored in a computer-readable storage medium. When executed, the program performs steps comprising the method embodiments described above; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
Finally, it should be noted that: the above embodiments are only used for illustrating the technical solutions of the present application, and not for limiting the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present application.
Claims (11)
1. A slope parameter extraction method is characterized by comprising the following steps:
acquiring a variable of a test signal propagating along a preset direction at a position to be tested along with time;
obtaining target coordinates corresponding to each propagation distance in the preset direction;
fitting and processing the target coordinates corresponding to each propagation distance to obtain the slope parameter of the position to be measured,
wherein, the obtaining of the destination coordinate corresponding to each propagation distance in the preset direction includes:
acquiring a variable of the test signal propagating along a preset direction along with time at a position to be tested, and generating a signal propagation diagram;
determining an initial target coordinate from the initial region to be measured of the signal propagation diagram, and taking the initial target coordinate as a current target coordinate;
taking the current target coordinate as a starting point, and acquiring M coordinates along the time direction to obtain a group of reference data points;
taking the corresponding position of the current target coordinate in the next propagation distance as a next starting point, moving the next starting point for N times by a preset step length, and acquiring M coordinates along the time direction every time of moving to obtain N groups of data points to be determined, wherein each group of data points to be determined comprises the M coordinates;
based on a block matching method, calculating matching coefficients between the group of reference data points and the N groups of data points to be determined respectively to obtain N matching coefficients;
determining a target coordinate corresponding to the next propagation distance according to the N matching coefficients;
and updating the current target coordinate to be a target coordinate corresponding to the next propagation distance, and repeatedly executing the step of obtaining M coordinates along the time direction by taking the current target coordinate as a starting point to obtain a group of reference data points until the target coordinate corresponding to each propagation distance in the signal propagation diagram is found.
2. The method according to claim 1, wherein the determining the destination coordinate corresponding to the next propagation distance according to the N matching coefficients comprises:
and determining a matching coefficient representing the best matching degree in the N matching coefficients, and taking the initial coordinate in the group of data points to be determined corresponding to the matching coefficient representing the best matching degree as the target coordinate corresponding to the next propagation distance.
3. The method according to claim 1 or 2, wherein before the fitting process the destination coordinates corresponding to the propagation distances, the method further comprises:
and screening the target coordinates corresponding to the propagation distances according to a preset screening scheme.
4. The method according to claim 3, wherein the screening the destination coordinates corresponding to each propagation distance according to a preset screening scheme includes:
sequentially fitting P adjacent target coordinates in the target coordinates corresponding to each propagation distance to obtain a corresponding fitting straight line set, wherein P is a positive integer and is smaller than the number of the target coordinates corresponding to each propagation distance; determining a target fitting straight line in the fitting straight line set, wherein the sum of the residual errors from the target coordinates corresponding to each propagation distance to the target fitting straight line is minimum; removing the target coordinates of which the distance to the target fitting straight line exceeds a preset distance threshold from the target coordinates corresponding to each propagation distance except P adjacent target coordinates for fitting the target fitting straight line; alternatively, the first and second electrodes may be,
calculating the difference value between the X-th target coordinate and the X-1 th target coordinate of the abscissa in the target coordinate corresponding to each propagation distance in turn, wherein X is a positive integer which is more than 2 and less than the number of each propagation distance; if the difference value between the X-th target coordinate and the X-1 th target coordinate is larger than a preset first threshold value, subtracting the corresponding 1 st to X-2 th target coordinate from the X-2 nd to X-1 th target coordinate, and if the average value of the absolute values of the obtained difference values is smaller than a preset second threshold value and X is smaller than a first preset clearing quantity threshold value, clearing the target coordinate before the X-th target coordinate; if X is larger than a second preset clearing quantity threshold value, clearing the target coordinate behind the X-th target coordinate; alternatively, the first and second electrodes may be,
if the difference value between the abscissa of the most-valued coordinates corresponding to the K continuous adjacent propagation distances is smaller than a preset third threshold value, only the K' th most-valued coordinate is reserved; the relation between K 'and K is subject to practical experience, and when K is 2K, K' is K; when K is 2K +1, K' is K + 1; wherein k is a positive integer.
5. The method according to claim 1 or 2, wherein the step of obtaining destination coordinates corresponding to each propagation distance in the preset direction is performed a plurality of times to obtain a plurality of sets of destination coordinates;
the fitting process of the destination coordinates corresponding to the propagation distances to obtain slope parameters of the position to be measured includes:
respectively fitting according to the multiple groups of target coordinates to obtain multiple slope parameters;
and calculating the median of the slope parameters to obtain the slope parameter of the position to be measured.
6. The method according to claim 1 or 2, characterized in that the method further comprises:
performing band-pass filtering on the variable;
the obtaining of the destination coordinates corresponding to each propagation distance in the preset direction includes:
before and after filtering, respectively repeating the step of obtaining target coordinates corresponding to each propagation distance in the preset direction for T times to obtain the target coordinates of each propagation distance before and after filtering;
the fitting process of the destination coordinates corresponding to the propagation distances to obtain slope parameters of the position to be measured includes:
according to the target coordinates corresponding to the T times of measurement before filtering, respectively fitting and processing the target coordinates corresponding to the propagation distances to obtain a pre-filtering slope set { B ] of the position to be measured1,..,Bi,..,BT}；
Respectively fitting and processing the target coordinates corresponding to each propagation distance according to the target coordinates corresponding to the T-time measurement after filtering to obtain the position to be measuredSet of filtered slopes { C1,..,Ci,..,CT}；
Set of pre-filtering slopes at position under test { B1,..,Bi,..,BTDetermining a median value as a slope parameter before filtering; set of filtered slopes at position under test { C1,..,Ci,..,CTDetermining a median value as a slope parameter after filtering;
taking the average value of the pre-filtering slope parameter and the post-filtering slope parameter as the slope parameter of the position to be detected; or selecting one of the pre-filtering slope parameter and the post-filtering slope parameter as the slope parameter of the position to be measured.
7. The method of claim 6, wherein selecting one of the pre-filter slope parameter and the post-filter slope parameter as the slope parameter of the location to be measured comprises:
determining the slope parameter of the position to be measured according to one or more of the following selection conditions in the pre-filtering slope parameter and the post-filtering slope parameter:
if the tissue elastic hardness value corresponding to the before-filtering slope parameter and the tissue elastic hardness value corresponding to the after-filtering slope parameter are in the same pathological stage, determining a slope parameter corresponding to the tissue elastic hardness value with smaller relative deviation as a slope parameter of the position to be detected; wherein the relative deviation is a ratio of the quartile range to the median of the slope set,
if the tissue elasticity hardness value corresponding to the slope parameter before filtering and the tissue elasticity hardness value corresponding to the slope parameter after filtering are not in the same pathological stage, judging whether the relative deviation of the tissue elasticity hardness values corresponding to the slope parameter after filtering is smaller than a fourth threshold value or not, and if the relative deviation is smaller than the fourth threshold value, determining the slope parameter after filtering as the slope parameter of the position to be detected; if not, judging whether the relative deviation of the tissue elastic hardness value corresponding to the slope parameter before filtering is smaller than a fifth threshold value, if so, determining the slope parameter before filtering as the slope parameter of the position to be detected,
if the tissue elastic hardness value corresponding to the pre-filtering slope parameter and the tissue elastic hardness value corresponding to the post-filtering slope parameter do not meet the selection condition, discarding each target coordinate corresponding to the current T times of measurement, and re-executing the step of obtaining the variable of the test signal propagating along the preset direction at the position to be tested along with time.
8. The method according to claim 1, wherein the fitting process of the destination coordinates corresponding to the propagation distances to obtain the slope parameter of the position to be measured includes:
fitting the target coordinates corresponding to the propagation distances to obtain a fitted straight line;
and taking the slope of the fitting straight line as a slope parameter of the position to be measured.
9. A slope parameter extraction device, characterized by comprising:
the acquisition module is used for acquiring a variable of the propagation of the test signal along a preset direction at the position to be tested along with time;
the determining module is used for obtaining a target coordinate corresponding to each propagation distance in the preset direction;
a processing module for fitting the target coordinates corresponding to the propagation distances to obtain slope parameters of the position to be measured,
wherein the determining module comprises:
the generating unit is used for generating a signal propagation diagram after acquiring a variable of the test signal propagating along a preset direction along with time at a position to be tested;
the determining unit is used for determining an initial target coordinate from the initial region to be measured of the signal propagation diagram and taking the initial target coordinate as the current target coordinate;
a lookup unit to:
taking the current target coordinate as a starting point, and acquiring M coordinates along the time direction to obtain a group of reference data points;
taking the corresponding position of the current target coordinate in the next propagation distance as a next starting point, moving the next starting point for N times by a preset step length, and acquiring M coordinates along the time direction every time of moving to obtain N groups of data points to be determined, wherein each group of data points to be determined comprises the M coordinates;
based on a block matching method, calculating matching coefficients between the group of reference data points and the N groups of data points to be determined respectively to obtain N matching coefficients;
determining a target coordinate corresponding to the next propagation distance according to the N matching coefficients;
and updating the current target coordinate to be a target coordinate corresponding to the next propagation distance, and repeatedly executing the step of obtaining M coordinates along the time direction by taking the current target coordinate as a starting point to obtain a group of reference data points until the target coordinate corresponding to each propagation distance in the signal propagation diagram is found.
10. A slope parameter extraction device, characterized by comprising: at least one processor and memory;
the memory stores computer-executable instructions; the at least one processor executes computer-executable instructions stored by the memory to perform the method of any of claims 1-8.
11. A computer-readable storage medium, characterized in that a computer program is stored in the computer-readable storage medium, which computer program, when executed, implements the method of any one of claims 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910002361.2A CN111401108B (en) | 2019-01-02 | 2019-01-02 | Slope parameter extraction method and device and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910002361.2A CN111401108B (en) | 2019-01-02 | 2019-01-02 | Slope parameter extraction method and device and computer readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111401108A CN111401108A (en) | 2020-07-10 |
| CN111401108B true CN111401108B (en) | 2021-04-06 |
Family
ID=71428230
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910002361.2A Active CN111401108B (en) | 2019-01-02 | 2019-01-02 | Slope parameter extraction method and device and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111401108B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111735526B (en) * | 2020-07-20 | 2020-11-17 | 深圳中科乐普医疗技术有限公司 | Ultrasonic elastography device and shear wave velocity measurement method for elasticity measurement |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011028973A2 (en) * | 2009-09-04 | 2011-03-10 | Cornell University | Methods and systems for functional imaging of cardiac tissue |
| CN104546014A (en) * | 2014-12-25 | 2015-04-29 | 中国科学院深圳先进技术研究院 | Shear wave velocity estimation method for biological tissue elasticity measurement |
| CN105266851A (en) * | 2015-10-29 | 2016-01-27 | 无锡海斯凯尔医学技术有限公司 | Method and device for detecting viscosity parameter of viscoelastic medium |
| WO2016033752A1 (en) * | 2014-09-03 | 2016-03-10 | 深圳迈瑞生物医疗电子股份有限公司 | Elasticity measurement and detection method and system |
| CN106594529A (en) * | 2017-01-20 | 2017-04-26 | 合肥工业大学 | Leakage signal propagation speed measurement method based on short-time energy and linear fitting |
| CN107616814A (en) * | 2017-08-25 | 2018-01-23 | 深圳中科乐普医疗技术有限公司 | A kind of biological tissue's shear-wave velocity measuring method and medical supersonic wave device |
| CN108158610A (en) * | 2018-01-16 | 2018-06-15 | 苏州国科昂卓医疗科技有限公司 | A kind of elastograph imaging method, device, equipment and ultrasound imaging probe |
-
2019
- 2019-01-02 CN CN201910002361.2A patent/CN111401108B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011028973A2 (en) * | 2009-09-04 | 2011-03-10 | Cornell University | Methods and systems for functional imaging of cardiac tissue |
| WO2016033752A1 (en) * | 2014-09-03 | 2016-03-10 | 深圳迈瑞生物医疗电子股份有限公司 | Elasticity measurement and detection method and system |
| CN104546014A (en) * | 2014-12-25 | 2015-04-29 | 中国科学院深圳先进技术研究院 | Shear wave velocity estimation method for biological tissue elasticity measurement |
| CN105266851A (en) * | 2015-10-29 | 2016-01-27 | 无锡海斯凯尔医学技术有限公司 | Method and device for detecting viscosity parameter of viscoelastic medium |
| CN106594529A (en) * | 2017-01-20 | 2017-04-26 | 合肥工业大学 | Leakage signal propagation speed measurement method based on short-time energy and linear fitting |
| CN107616814A (en) * | 2017-08-25 | 2018-01-23 | 深圳中科乐普医疗技术有限公司 | A kind of biological tissue's shear-wave velocity measuring method and medical supersonic wave device |
| CN108158610A (en) * | 2018-01-16 | 2018-06-15 | 苏州国科昂卓医疗科技有限公司 | A kind of elastograph imaging method, device, equipment and ultrasound imaging probe |
Non-Patent Citations (2)
| Title |
|---|
| Viscoelasticity Measurement in Rat Rivers Using Shear-Wave US Elastography;Katsutoshi Sugimoto 等;《Ultrasound in Medicine&Biology》;20180930;第44卷(第9期);2018-2024 * |
| 超声实时剪切波弹性成像杨氏模量最大值在甲状腺实性结节良恶性鉴别诊断中的应用;薛海燕;《中国优秀硕士学位论文全文数据库 医药卫生科技辑》;20180115(第01期);E060-47 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111401108A (en) | 2020-07-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2549213C2 (en) | Processing seismic data | |
| KR100893286B1 (en) | Image processing system and method | |
| US20140204700A1 (en) | Automatic Processing of Ultrasonic Data | |
| US7632231B2 (en) | Ultrasonic strain imaging device and method providing parallel displacement processing | |
| CN103429163B (en) | Device and method for determining actual tissue layer boundaries of a body | |
| CN104880172B (en) | Measuring flatness of road surface method and device based on Kalman filtering | |
| JP5692079B2 (en) | Displacement estimation method and displacement estimation apparatus | |
| JP2008142519A (en) | Ultrasound diagnostic apparatus and volume data processing method | |
| CN111401108B (en) | Slope parameter extraction method and device and computer readable storage medium | |
| EP2903530A1 (en) | Shear wave attenuation from k-space analysis system | |
| CN110400340B (en) | Depth map processing method and device | |
| CN102204820A (en) | Recording apparatus and recording method thereof | |
| CN111399038B (en) | Slope parameter extraction method and device and computer readable storage medium | |
| CN116908853B (en) | High coherence point selection method, device and equipment | |
| CN111399039B (en) | Slope parameter extraction method, device and storage medium | |
| CN111388012B (en) | Method, device and system for detecting tissue hardness | |
| CN112578435A (en) | Rock ultrasonic testing first arrival picking method and system | |
| CN114972567B (en) | Medical ultrasonic CT multi-parameter image reconstruction method based on wave equation | |
| CN110702786A (en) | Beam structure damage identification method based on multi-scale singular attractor prediction error | |
| CN107169978B (en) | Ultrasonic image edge detection method and system | |
| JP7458641B2 (en) | Method for estimating electrode placement in biological tissue | |
| CN111189912B (en) | Emission reference ultrasonic detection method, device and storage medium | |
| CN111399051B (en) | Data processing method, apparatus and medium | |
| JP2669204B2 (en) | Exploration equipment | |
| RU2577561C1 (en) | Method of measuring pulse response function structure in time in heterogeneous environment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |