WO2021103493A1 - Shear wave-based imaging method, system and apparatus - Google Patents

Shear wave-based imaging method, system and apparatus Download PDF

Info

Publication number
WO2021103493A1
WO2021103493A1 PCT/CN2020/096261 CN2020096261W WO2021103493A1 WO 2021103493 A1 WO2021103493 A1 WO 2021103493A1 CN 2020096261 W CN2020096261 W CN 2020096261W WO 2021103493 A1 WO2021103493 A1 WO 2021103493A1
Authority
WO
WIPO (PCT)
Prior art keywords
shear wave
propagation
detection area
shear
energy
Prior art date
Application number
PCT/CN2020/096261
Other languages
French (fr)
Chinese (zh)
Inventor
刘德清
朱超超
冯乃章
Original Assignee
深圳开立生物医疗科技股份有限公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 深圳开立生物医疗科技股份有限公司 filed Critical 深圳开立生物医疗科技股份有限公司
Publication of WO2021103493A1 publication Critical patent/WO2021103493A1/en

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties

Definitions

  • the present invention relates to the field of tissue lesion analysis, in particular to a shear wave-based imaging method, system and device.
  • Shear wave elastography technology based on acoustic radiation force is an ultrasonic elastography technology for evaluating tissue elasticity, and it is widely used in tissue lesion analysis.
  • the principle of shear wave elastography technology based on acoustic radiation force is: after the probe emits high-energy ultrasound to the tissue, under the action of the acoustic radiation force and the shear stress of the tissue, the tissue in a specific area will propagate to the surroundings. Because of the correlation between tissue hardness and shear wave speed, the tissue hardness can be analyzed by detecting the shear wave speed to determine whether there is tissue disease.
  • the principle of shear wave velocity detection is: at several observation points with a known distance of the forward propagation path of the shear wave, observe the time when the shear wave propagates to several observation points, and then according to The time difference between the shear wave propagation to several observation points and the space distance of several observation points, calculate the average velocity of the shear wave at these observation points. It can be seen that the prior art detects the forward propagation velocity of the shear wave from the excitation source of the shear wave. However, in the actual propagation process of the shear wave, due to the difference in the hardness of different tissues, the shear wave will be based on the hardness of the different tissues. The magnitude of the difference, different degrees of reflection occur. Since the prior art detects the shear wave velocity without considering the reflection that occurs during the propagation of the shear wave, the uniformity of the tissue hardness cannot be accurately analyzed, which is not conducive to the analysis of tissue lesions.
  • the purpose of the present invention is to provide a shear wave-based imaging method, system and device to obtain the shear wave propagation reflection coefficient of the shear wave detection area, so that the shear wave detection area can be accurately determined according to the shear wave propagation reflection coefficient
  • the uniformity of the tissue hardness is conducive to the analysis of regional tissue lesions.
  • the present invention provides a shear wave-based imaging method, including:
  • the shear wave propagation reflection coefficient According to the ratio of the reflected wave energy to the incident wave energy, obtain the shear wave propagation reflection coefficient, and return to the step of selecting a lateral local area in the shear wave detection area until the entire shear wave detection area is traversed, To obtain the shear wave propagation reflection coefficient of the entire shear wave detection area;
  • a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area is generated.
  • the process of selecting a lateral local area in the shear wave detection area and obtaining the incident wave energy and reflected wave energy of the shear wave corresponding to the lateral local area includes:
  • the process of periodically acquiring the tissue displacement of each detection point in the shear wave detection area under the action of the shear wave propagation includes:
  • the tissue displacement of each detection point under the propagation action of the shear wave is obtained according to the periodically obtained position information.
  • the process of obtaining the forward propagation information matrix and the backward propagation information matrix of the shear wave according to the two-dimensional time-space propagation matrix includes:
  • the forward propagation information matrix and the backward propagation information matrix are respectively extracted from the two-dimensional wk frequency domain matrix.
  • the process of correspondingly obtaining the incident wave energy and the reflected wave energy of the shear wave according to the forward propagation information matrix and the backward propagation information matrix includes:
  • the two-dimensional integration of the time frequency and the space frequency is performed on the absolute value of the back propagation information matrix to obtain the reflected wave energy of the shear wave.
  • the process of generating a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area according to the shear wave propagation reflection coefficient includes:
  • the imaging method further includes:
  • the shear wave propagation attenuation is corrected according to the reflected wave energy corresponding to the two adjacent lateral local regions.
  • the present invention also provides a shear wave-based imaging system, including:
  • An energy acquisition module configured to select a lateral local area in the shear wave detection area, and obtain the incident wave energy and the reflected wave energy of the shear wave corresponding to the lateral local area;
  • the coefficient obtaining module is used to obtain the shear wave propagation reflection coefficient according to the ratio of the reflected wave energy to the incident wave energy, and return to execute the energy obtaining module until the entire shear wave detection area is traversed to obtain The shear wave propagation reflection coefficient of the entire shear wave detection area;
  • the tissue imaging module is configured to generate a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area according to the shear wave propagation reflection coefficient.
  • the present invention also provides a shear wave-based imaging device, including:
  • Memory used to store computer programs
  • the processor is configured to implement the steps of any of the foregoing shear wave-based imaging methods when the computer program is executed.
  • the present invention provides an imaging method based on shear waves. Considering that in the actual propagation process of shear waves, due to the difference in the hardness of different tissues, the shear waves will have different degrees of hardness according to the difference in hardness of different tissues. Therefore, the reflection coefficient of shear wave propagation can directly and accurately reflect the uniformity of tissue hardness. Based on this, the present application obtains the shear wave propagation reflection coefficient of the shear wave detection area, so that the uniformity of the tissue hardness of the shear wave detection area can be accurately determined according to the shear wave propagation reflection coefficient, which is beneficial to the analysis of regional tissue lesions.
  • the present invention also provides a shear wave-based imaging system and device, which have the same beneficial effects as the foregoing imaging method.
  • FIG. 1 is a flowchart of a shear wave-based imaging method provided by an embodiment of the present invention
  • FIG. 2 is a diagram of a shear wave propagation waveform corresponding to a two-dimensional time-space propagation matrix provided by an embodiment of the present invention
  • FIG. 3 is a schematic diagram of a two-dimensional wk frequency domain matrix provided by an embodiment of the present invention.
  • FIG. 4 is a schematic diagram of a forward propagation information matrix provided by an embodiment of the present invention.
  • FIG. 5 is a schematic diagram of a backpropagation information matrix provided by an embodiment of the present invention.
  • FIG. 6 is a schematic structural diagram of a shear wave-based imaging system provided by an embodiment of the present invention.
  • the core of the present invention is to provide a shear wave-based imaging method, system and device to obtain the shear wave propagation reflection coefficient of the shear wave detection area, so that the shear wave detection area can be accurately determined according to the shear wave propagation reflection coefficient
  • the uniformity of the tissue hardness is beneficial to the analysis of regional tissue lesions.
  • FIG. 1 is a flowchart of a shear wave-based imaging method according to an embodiment of the present invention.
  • the shear wave-based imaging method includes:
  • Step S1 Select a lateral local area in the shear wave detection area, and obtain the incident wave energy and the reflected wave energy of the shear wave corresponding to the lateral local area.
  • the shear waves will be reflected to different degrees according to the difference in hardness of different tissues, so the shear waves occur during the propagation process.
  • the reflection situation can directly and accurately reflect the uniformity of the tissue hardness. More specifically, the greater the difference in hardness between two adjacent tissues, the greater the degree of reflection of the shear wave at the junction of the two tissues, that is, under the same incident wave energy, the greater the difference in hardness of the shear wave is. The greater the energy of the reflected wave at the junction of the two tissues.
  • the ratio of the reflected wave energy to the incident wave energy can be used, that is, the shear wave propagation reflection coefficient represents the reflection of the shear wave during the propagation process. Therefore, in this embodiment, in order to obtain the shear wave propagation reflection coefficient of a local tissue area, the incident wave energy and the reflected wave energy of the shear wave in the local tissue area should be obtained first.
  • shear wave detection area combining the shear wave propagation reflection coefficients of multiple lateral local areas of the shear wave detection area can directly and accurately determine the organization of the entire shear wave detection area Uniformity of hardness. More specifically, in this embodiment, in the shear wave detection area, a small section of the lateral local area is first selected, and then the incident wave energy and the reflected wave energy of the shear wave corresponding to the selected lateral local area are obtained. The shear wave propagation reflection coefficient of the selected lateral local area.
  • Step S2 According to the ratio of the reflected wave energy to the incident wave energy, obtain the shear wave propagation reflection coefficient, and determine whether to traverse the entire shear wave detection area, if not, return to the execution and select the lateral part in the shear wave detection area Step of the area; if yes, go to step S3.
  • the ratio of the reflected wave energy to the incident wave energy is calculated, thereby obtaining the shear wave energy of the selected lateral local area.
  • the present embodiment needs to traverse the entire shear wave detection area. If the entire shear wave detection area is not traversed, then continue to detect the shear wave that has not been traversed. Select a section of lateral local area in the area, and obtain the shear wave propagation reflection coefficient of this lateral local area according to the above steps; if the entire shear wave detection area has been traversed, the shear wave of the entire shear wave detection area can be obtained Propagation reflection coefficient.
  • the shear wave propagation reflection coefficient of each lateral local area may be different, and the shear wave propagation reflection coefficient of each lateral local area is determined by the tissue hardness of the lateral local area, that is, it has tissue elasticity. The uniqueness and referability of the detection, so as to detect the elasticity of the tissue more accurately.
  • Step S3 Generate a tissue image characterizing the uniformity of the tissue hardness in the shear wave detection area according to the shear wave propagation reflection coefficient.
  • this embodiment can generate a tissue image that clearly characterizes the uniformity of the tissue hardness of the shear wave detection area according to the shear wave propagation reflection coefficient of the entire shear wave detection area, so that the focus tissue can be found in time by analyzing the tissue image.
  • shear wave can be generated by acoustic radiation force, by external vibration, or by the internal motion tissue of the body.
  • the embodiment of the present invention provides a shear wave-based imaging method. Considering that in the actual propagation process of the shear wave, due to the difference in the hardness of different tissues, the shear wave will be different according to the difference in the hardness of the different tissues. Degree of reflection, so the shear wave propagation reflection coefficient can directly and accurately reflect the uniformity of tissue hardness. Based on this, the present application obtains the shear wave propagation reflection coefficient of the shear wave detection area, so that the uniformity of the tissue hardness of the shear wave detection area can be accurately determined according to the shear wave propagation reflection coefficient, which is beneficial to the analysis of regional tissue lesions.
  • the process of selecting a lateral local area in the shear wave detection area and obtaining the incident wave energy and reflected wave energy of the shear wave corresponding to the lateral local area includes:
  • the tissue displacement of each detection point in the shear wave detection area can reflect the shear wave in the shear wave detection area.
  • the energy transmission situation Therefore, in this embodiment, in order to understand the energy propagation of the shear wave in the shear wave detection area, periodically obtain the tissue displacement of each detection point in the shear wave detection area under the action of the shear wave propagation, and according to each detection point The one-dimensional time-displacement waveform of each detection point is obtained one by one of the tissue displacement.
  • the energy propagation of the shear wave in the entire shear wave detection area is specifically obtained by integrating the energy propagation of the shear wave in multiple lateral local areas within it.
  • the transverse direction is selected in the shear wave detection area.
  • the one-dimensional time displacement waveforms of the detection points in the selected lateral local area are combined to obtain the two-dimensional time-space propagation matrix corresponding to the lateral local area. Since the two-dimensional time-space propagation matrix corresponding to the lateral local area contains the forward propagation information matrix and the back-propagation information matrix of the shear wave in the lateral local area, this embodiment can be based on the two-dimensional time and space propagation corresponding to the lateral local area.
  • Matrix obtain the forward propagation information matrix and the reverse propagation information matrix of the shear wave in the lateral local area, and obtain the incident wave of the shear wave in the lateral local area according to the forward propagation information matrix of the shear wave in the lateral local area Energy, in the same way, the reflected wave energy of the shear wave in the lateral local area is obtained according to the back propagation information matrix of the shear wave in the lateral local area. Then use this technical means to traverse the entire shear wave detection area, so as to obtain the energy propagation of the shear wave in the entire shear wave detection area.
  • this embodiment analyzes the shear wave propagation waveform diagram corresponding to the two-dimensional time-space propagation matrix in Fig. 2 (the data at each position on the shear wave propagation waveform forms a two-dimensional time-space propagation matrix), and the shear wave reflection can be obtained.
  • the process of periodically acquiring the tissue displacement of each detection point in the shear wave detection area under the action of the shear wave propagation includes:
  • the tissue displacement of each detection point under the propagation of shear waves is obtained.
  • the technical means of this embodiment to obtain the tissue displacement of each detection point in the shear wave detection area under the action of the shear wave propagation is: first, according to the preset repetition frequency, repeatedly transmit to the shear wave detection area for detection.
  • the detection beam of the shear wave signal that is, the detection beam is periodically transmitted to the shear wave detection area.
  • the detection echo signal is received, and the multi-beam technology (a kind of radar overall technology for measuring target coordinates) is used to obtain the position information of each detection point in the shear wave detection area according to the returned detection echo signal.
  • the tissue displacement of the detection point in this period of time is based on the periodicity of each detection point. Time position information, the tissue displacement of each detection point that changes with time under the action of shear wave propagation is obtained.
  • the process of obtaining the forward propagation information matrix and the backward propagation information matrix of the shear wave according to the two-dimensional time-space propagation matrix includes:
  • the forward propagation information matrix and the backward propagation information matrix are respectively extracted from the two-dimensional wk frequency domain matrix.
  • the two-dimensional time-space propagation matrix is processed by the directional filter, that is, the two-dimensional time-space propagation matrix is subjected to 2D Fourier transform.
  • the two-dimensional wk frequency domain matrix that is, the two-dimensional time frequency-spatial frequency matrix, as shown in Figure 3.
  • the information matrix AMat of the second and fourth quadrants of the two-dimensional wk frequency domain matrix represents the forward propagation of the shear wave.
  • the information matrix BMat of the first quadrant and the third quadrant represents the information of the shear wave back propagation, so this embodiment can extract the forward propagation information matrix AMat from the two-dimensional wk frequency domain matrix respectively ( Figure 4) And the back-propagation information matrix BMat ( Figure 5).
  • this embodiment uses different masks to respectively extract the forward propagation information matrix AMat and the back propagation information matrix BMat in the two-dimensional wk frequency domain matrix, where the value of the mask: extract the forward propagation information matrix AMat
  • the first and third quadrants are set to 0 and the second and fourth quadrants are set to 1
  • the mask is multiplied by the two-dimensional wk frequency domain matrix to obtain the forward propagation information matrix AMat
  • the back propagation information matrix BMat is extracted, the first and third quadrants
  • the value is 1 and the value of the second and fourth quadrants is 0, then the mask is multiplied by the two-dimensional wk frequency domain matrix to obtain the back propagation information matrix BMat.
  • the process of correspondingly obtaining the incident wave energy and the reflected wave energy of the shear wave according to the forward propagation information matrix and the backward propagation information matrix includes:
  • the two-dimensional integration of time frequency and space frequency is performed on the absolute value of the backpropagation information matrix to obtain the reflected wave energy of the shear wave.
  • the calculation method for obtaining the incident wave energy Apower of the shear wave according to the forward propagation information matrix AMat is:
  • the calculation method for obtaining the reflected wave energy Bpower of the shear wave according to the backpropagation information matrix BMat is:
  • the process of generating a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area according to the shear wave propagation reflection coefficient includes:
  • the preset in this embodiment is set in advance, and only needs to be set once, unless it needs to be modified according to the actual situation, it does not need to be reset.
  • a range is set in advance: [min, max], so that the shear wave propagation reflection coefficients in the shear wave detection area are all normalized to the set range according to a certain rule, so that one One gets the normalized value of the reflection coefficient of shear wave propagation.
  • the correspondence relationship between the normalized value of the reflection coefficient and the color of the shear wave detection area is also set in advance, referred to as the coefficient color correspondence relationship.
  • the maximum value max corresponds to the color level 255
  • the minimum value min corresponds to the color level 0
  • the color level Rflct corresponding to the normalized value of the reflection coefficient RflctRatio is: In order to determine the color corresponding to the normalized value of the reflection coefficient.
  • the colors corresponding to the normalized values of the reflection coefficients are determined one-to-one according to the coefficient color correspondence relationship, and then the colors of the shear wave detection area are filled correspondingly according to each color. And control the display screen to display the shear wave detection area filled with the corresponding color for the staff to view, so as to more intuitively show the uniformity of the tissue hardness of the shear wave detection area.
  • gray values may represent different color levels; different colors may also represent different color levels.
  • the imaging method further includes:
  • the shear wave propagation attenuation is corrected according to the reflected wave energy corresponding to two adjacent lateral local areas.
  • ; the absolute value of the shear wave propagation attenuation is sequentially subtracted from the reflected wave energy of the adjacent transverse local areas A and B to obtain the corrected shear Shear wave propagation attenuation, that is, the corrected shear wave propagation attenuation
  • FIG. 6 is a schematic structural diagram of a shear wave-based imaging system provided by an embodiment of the present invention.
  • the shear wave-based imaging system includes:
  • the energy acquisition module 1 is used to select a lateral local area in the shear wave detection area, and obtain the incident wave energy and the reflected wave energy of the shear wave corresponding to the lateral local area;
  • the coefficient calculation module 2 is used to obtain the shear wave propagation reflection coefficient according to the ratio of the reflected wave energy to the incident wave energy, and return to the execution of the energy acquisition module until the entire shear wave detection area is traversed to obtain the entire shear wave detection The shear wave propagation reflection coefficient of the area;
  • the tissue imaging module 3 is used to generate a tissue image characterizing the uniformity of the tissue hardness in the shear wave detection area according to the shear wave propagation reflection coefficient.
  • An embodiment of the present application also provides a shear wave-based imaging device, including:
  • Memory used to store computer programs
  • the processor is used to implement the steps of any of the foregoing shear wave-based imaging methods when the computer program is executed.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

A shear wave-based imaging method. Considering that in the actual propagation process of shear waves, due to the difference in hardness of different tissues, the shear waves will be reflected in different degrees according to the difference in hardness of different tissues, and therefore, the shear wave propagation reflection coefficient can directly and accurately reflect the uniformity of tissue hardness. On this basis, according to the shear wave-based imaging method, the shear wave propagation reflection coefficient of a shear wave detection area is obtained, so that the uniformity of the tissue hardness of the shear wave detection area can be accurately determined according to the shear wave propagation reflection coefficient, facilitating regional tissue lesions analysis. In addition, also disclosed are a shear wave-based imaging system and apparatus, which have the same beneficial effects as the shear wave-based imaging method.

Description

一种基于剪切波的成像方法、系统及装置An imaging method, system and device based on shear wave
本申请要求于2019年11月27日提交至中国专利局、申请号为201911184432.1、发明名称为“一种基于剪切波的成像方法、系统及装置”的中国专利申请的优先权，其全部内容通过引用结合在本申请中。This application claims the priority of a Chinese patent application filed with the Chinese Patent Office on November 27, 2019, the application number is 201911184432.1, and the invention title is "a shear wave-based imaging method, system, and device", and its entire contents Incorporated in this application by reference.
技术领域Technical field
本发明涉及组织病变分析领域，特别是涉及一种基于剪切波的成像方法、系统及装置。The present invention relates to the field of tissue lesion analysis, in particular to a shear wave-based imaging method, system and device.
背景技术Background technique
基于声辐射力的剪切波弹性成像技术是一种评估组织弹性的超声弹性成像技术，其被广泛应用于组织病变分析。基于声辐射力的剪切波弹性成像技术的原理为：由探头向组织发射高能量的超声波后，在声辐射力和组织的剪切应力的作用下，特定区域内的组织会产生向四周传播的振动，从而产生剪切波，而由于组织硬度与剪切波的速度存在着关联关系，因此可以通过检测剪切波的速度分析组织硬度，进而确定是否存在组织病变。Shear wave elastography technology based on acoustic radiation force is an ultrasonic elastography technology for evaluating tissue elasticity, and it is widely used in tissue lesion analysis. The principle of shear wave elastography technology based on acoustic radiation force is: after the probe emits high-energy ultrasound to the tissue, under the action of the acoustic radiation force and the shear stress of the tissue, the tissue in a specific area will propagate to the surroundings. Because of the correlation between tissue hardness and shear wave speed, the tissue hardness can be analyzed by detecting the shear wave speed to determine whether there is tissue disease.
现有技术中，剪切波的速度检测原理为：在剪切波的正向传播路径的已知距离的几个观测点上，观察剪切波分别传播到达几个观测点的时间，然后根据剪切波传播到达几个观测点的时间差和几个观测点的空间距离，计算出剪切波在这几个观测点的平均速度。可见，现有技术检测的是剪切波从剪切波激励源正向传播的速度，但是在剪切波的实际传播过程中，由于不同组织硬度的差异，会导致剪切波按照不同组织硬度差异的大小，发生不同程度的反射。由于现有技术检测剪切波的速度时并未考虑剪切波传播过程中发生的反射情况，所以无法准确分析出组织硬度的均匀性，不利于组织病变分析。In the prior art, the principle of shear wave velocity detection is: at several observation points with a known distance of the forward propagation path of the shear wave, observe the time when the shear wave propagates to several observation points, and then according to The time difference between the shear wave propagation to several observation points and the space distance of several observation points, calculate the average velocity of the shear wave at these observation points. It can be seen that the prior art detects the forward propagation velocity of the shear wave from the excitation source of the shear wave. However, in the actual propagation process of the shear wave, due to the difference in the hardness of different tissues, the shear wave will be based on the hardness of the different tissues. The magnitude of the difference, different degrees of reflection occur. Since the prior art detects the shear wave velocity without considering the reflection that occurs during the propagation of the shear wave, the uniformity of the tissue hardness cannot be accurately analyzed, which is not conducive to the analysis of tissue lesions.
因此，如何提供一种解决上述技术问题的方案是本领域的技术人员目前需要解决的问题。Therefore, how to provide a solution to the above-mentioned technical problems is a problem that needs to be solved by those skilled in the art at present.
发明内容Summary of the invention
本发明的目的是提供一种基于剪切波的成像方法、系统及装置，获取剪切波检测区域的剪切波传播反射系数，从而可根据剪切波传播反射系数准确确定剪切波检测区域的组织硬度均匀性，利于区域组织病变分析。The purpose of the present invention is to provide a shear wave-based imaging method, system and device to obtain the shear wave propagation reflection coefficient of the shear wave detection area, so that the shear wave detection area can be accurately determined according to the shear wave propagation reflection coefficient The uniformity of the tissue hardness is conducive to the analysis of regional tissue lesions.
为解决上述技术问题，本发明提供了一种基于剪切波的成像方法，包括：In order to solve the above technical problems, the present invention provides a shear wave-based imaging method, including:
在剪切波检测区域内选取横向局部区域，并获取所述横向局部区域对应的剪切波的入射波能量和反射波能量；Selecting a lateral local area in the shear wave detection area, and obtaining the incident wave energy and the reflected wave energy of the shear wave corresponding to the lateral local area;
根据所述反射波能量和所述入射波能量的比值，得到剪切波传播反射系数，并返回在所述剪切波检测区域内选取横向局部区域的步骤，直至遍历整个剪切波检测区域，以得到整个剪切波检测区域的剪切波传播反射系数；According to the ratio of the reflected wave energy to the incident wave energy, obtain the shear wave propagation reflection coefficient, and return to the step of selecting a lateral local area in the shear wave detection area until the entire shear wave detection area is traversed, To obtain the shear wave propagation reflection coefficient of the entire shear wave detection area;
根据所述剪切波传播反射系数生成表征所述剪切波检测区域的组织硬度均匀性的组织图像。According to the shear wave propagation reflection coefficient, a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area is generated.
优选地，所述在剪切波检测区域内选取横向局部区域，并获取所述横向局部区域对应的剪切波的入射波能量和反射波能量的过程，包括：Preferably, the process of selecting a lateral local area in the shear wave detection area and obtaining the incident wave energy and reflected wave energy of the shear wave corresponding to the lateral local area includes:
周期性获取剪切波检测区域内各检测点在剪切波的传播作用下的组织位移，并根据所述组织位移得到所述各检测点的一维时间位移波形；Periodically acquiring the tissue displacement of each detection point in the shear wave detection area under the propagation of the shear wave, and obtaining a one-dimensional time displacement waveform of each detection point according to the tissue displacement;
在所述剪切波检测区域内选取横向局部区域，并将所述横向局部区域内各检测点的一维时间位移波形进行组合，得到所述横向局部区域对应的二维时间空间传播矩阵；Selecting a lateral local area in the shear wave detection area, and combining the one-dimensional time displacement waveforms of each detection point in the lateral local area to obtain a two-dimensional time-space propagation matrix corresponding to the lateral local area;
根据所述二维时间空间传播矩阵得到剪切波的正向传播信息矩阵和反向传播信息矩阵，并根据所述正向传播信息矩阵和所述反向传播信息矩阵相应得到剪切波的入射波能量和反射波能量。Obtain the forward propagation information matrix and the backward propagation information matrix of the shear wave according to the two-dimensional time and space propagation matrix, and obtain the incident shear wave according to the forward propagation information matrix and the backward propagation information matrix. Wave energy and reflected wave energy.
优选地，所述周期性获取剪切波检测区域内各检测点在剪切波的传播作用下的组织位移的过程，包括：Preferably, the process of periodically acquiring the tissue displacement of each detection point in the shear wave detection area under the action of the shear wave propagation includes:
向剪切波检测区域周期性发射用于检测剪切波信号的检测波束，并将返回的检测回波信号采用多波束技术得到在所述剪切波检测区域内各检测点的位置信息；Periodically transmit a detection beam for detecting the shear wave signal to the shear wave detection area, and use multi-beam technology to obtain the position information of each detection point in the shear wave detection area for the returned detection echo signal;
根据周期性得到的位置信息得到所述各检测点在剪切波的传播作用下的组织位移。The tissue displacement of each detection point under the propagation action of the shear wave is obtained according to the periodically obtained position information.
优选地，所述根据所述二维时间空间传播矩阵得到剪切波的正向传播信息矩阵和反向传播信息矩阵的过程，包括：Preferably, the process of obtaining the forward propagation information matrix and the backward propagation information matrix of the shear wave according to the two-dimensional time-space propagation matrix includes:
对所述二维时间空间传播矩阵进行2D傅里叶变换，得到包括剪切波的正向传播信息矩阵和反向传播信息矩阵的二维wk频域矩阵；Performing a 2D Fourier transform on the two-dimensional time-space propagation matrix to obtain a two-dimensional wk frequency domain matrix including a forward propagation information matrix of the shear wave and a backward propagation information matrix;
从所述二维wk频域矩阵中分别提取出所述正向传播信息矩阵和所述反向传播信息矩阵。The forward propagation information matrix and the backward propagation information matrix are respectively extracted from the two-dimensional wk frequency domain matrix.
优选地，所述根据所述正向传播信息矩阵和所述反向传播信息矩阵相应得到剪切波的入射波能量和反射波能量的过程，包括：Preferably, the process of correspondingly obtaining the incident wave energy and the reflected wave energy of the shear wave according to the forward propagation information matrix and the backward propagation information matrix includes:
分别对所述正向传播信息矩阵和所述反向传播信息矩阵取绝对值；Taking absolute values for the forward propagation information matrix and the backward propagation information matrix respectively;
对所述正向传播信息矩阵的绝对值进行时间频率和空间频率的二维积分，得到剪切波的入射波能量；Performing two-dimensional integration of time frequency and space frequency on the absolute value of the forward propagation information matrix to obtain the incident wave energy of the shear wave;
对所述反向传播信息矩阵的绝对值进行时间频率和空间频率的二维积分，得到剪切波的反射波能量。The two-dimensional integration of the time frequency and the space frequency is performed on the absolute value of the back propagation information matrix to obtain the reflected wave energy of the shear wave.
优选地，所述根据所述剪切波传播反射系数生成表征所述剪切波检测区域的组织硬度均匀性的组织图像的过程，包括：Preferably, the process of generating a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area according to the shear wave propagation reflection coefficient includes:
将剪切波检测区域的剪切波传播反射系数均归一化到预设范围内，得到剪切波传播的反射系数归一值；Normalize the shear wave propagation reflection coefficients in the shear wave detection area to a preset range to obtain the normalized value of the shear wave propagation reflection coefficient;
根据预设系数颜色对应关系确定与所述反射系数归一值对应的颜色；Determining the color corresponding to the normalized value of the reflection coefficient according to the preset coefficient color correspondence relationship;
根据所述颜色相应填充所述剪切波检测区域的颜色，并控制显示屏显示填充有相应颜色的剪切波检测区域。Fill the color of the shear wave detection area correspondingly according to the color, and control the display screen to display the shear wave detection area filled with the corresponding color.
优选地，所述成像方法还包括：Preferably, the imaging method further includes:
在剪切波正向传播方向上，将相邻的横向局部区域对应的入射波能量相减，得到剪切波传播衰减量；In the forward propagation direction of the shear wave, subtract the incident wave energy corresponding to the adjacent lateral local area to obtain the shear wave propagation attenuation;
根据相邻的两个所述横向局部区域对应的反射波能量修正所述剪切波传播衰减量。The shear wave propagation attenuation is corrected according to the reflected wave energy corresponding to the two adjacent lateral local regions.
为解决上述技术问题，本发明还提供了一种基于剪切波的成像系统，包括：To solve the above technical problems, the present invention also provides a shear wave-based imaging system, including:
能量获取模块，用于在剪切波检测区域内选取横向局部区域，并获取所述横向局部区域对应的剪切波的入射波能量和反射波能量；An energy acquisition module, configured to select a lateral local area in the shear wave detection area, and obtain the incident wave energy and the reflected wave energy of the shear wave corresponding to the lateral local area;
系数求取模块，用于根据所述反射波能量和所述入射波能量的比值，得到剪切波传播反射系数，并返回执行所述能量获取模块，直至遍历整个剪切波检测区域，以得到整个剪切波检测区域的剪切波传播反射系数；The coefficient obtaining module is used to obtain the shear wave propagation reflection coefficient according to the ratio of the reflected wave energy to the incident wave energy, and return to execute the energy obtaining module until the entire shear wave detection area is traversed to obtain The shear wave propagation reflection coefficient of the entire shear wave detection area;
组织成像模块，用于根据所述剪切波传播反射系数生成表征所述剪切波检测区域的组织硬度均匀性的组织图像。The tissue imaging module is configured to generate a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area according to the shear wave propagation reflection coefficient.
为解决上述技术问题，本发明还提供了一种基于剪切波的成像装置，包括：To solve the above technical problems, the present invention also provides a shear wave-based imaging device, including:
存储器，用于存储计算机程序；Memory, used to store computer programs;
处理器，用于在执行所述计算机程序时实现上述任一种基于剪切波的成像方法的步骤。The processor is configured to implement the steps of any of the foregoing shear wave-based imaging methods when the computer program is executed.
本发明提供了一种基于剪切波的成像方法，考虑到在剪切波的实际传播过程中，由于不同组织硬度的差异，会导致剪切波按照不同组织硬度差异的大小，发生不同程度的反射，所以剪切波传播反射系数可直接准确地反映出组织硬度的均匀性。基于此，本申请获取剪切波检测区域的剪切波传播反射系数，从而可根据剪切波传播反射系数准确确定剪切波检测区域的组织硬度均匀性，利于区域组织病变分析。The present invention provides an imaging method based on shear waves. Considering that in the actual propagation process of shear waves, due to the difference in the hardness of different tissues, the shear waves will have different degrees of hardness according to the difference in hardness of different tissues. Therefore, the reflection coefficient of shear wave propagation can directly and accurately reflect the uniformity of tissue hardness. Based on this, the present application obtains the shear wave propagation reflection coefficient of the shear wave detection area, so that the uniformity of the tissue hardness of the shear wave detection area can be accurately determined according to the shear wave propagation reflection coefficient, which is beneficial to the analysis of regional tissue lesions.
本发明还提供了一种基于剪切波的成像系统及装置，与上述成像方法具有相同的有益效果。The present invention also provides a shear wave-based imaging system and device, which have the same beneficial effects as the foregoing imaging method.
附图说明Description of the drawings
为了更清楚地说明本发明实施例中的技术方案，下面将对现有技术和实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。In order to more clearly describe the technical solutions in the embodiments of the present invention, the following will briefly introduce the prior art and the drawings needed in the embodiments. Obviously, the drawings in the following description are only some of the present invention. Embodiments, for those of ordinary skill in the art, without creative work, other drawings can be obtained based on these drawings.
图1为本发明实施例提供的一种基于剪切波的成像方法的流程图；FIG. 1 is a flowchart of a shear wave-based imaging method provided by an embodiment of the present invention;
图2为本发明实施例提供的一种与二维时间空间传播矩阵对应的剪切波传播波形图；FIG. 2 is a diagram of a shear wave propagation waveform corresponding to a two-dimensional time-space propagation matrix provided by an embodiment of the present invention;
图3为本发明实施例提供的一种二维wk频域矩阵示意图；3 is a schematic diagram of a two-dimensional wk frequency domain matrix provided by an embodiment of the present invention;
图4为本发明实施例提供的一种正向传播信息矩阵示意图；4 is a schematic diagram of a forward propagation information matrix provided by an embodiment of the present invention;
图5为本发明实施例提供的一种反向传播信息矩阵示意图；FIG. 5 is a schematic diagram of a backpropagation information matrix provided by an embodiment of the present invention;
图6为本发明实施例提供的一种基于剪切波的成像系统的结构示意图。FIG. 6 is a schematic structural diagram of a shear wave-based imaging system provided by an embodiment of the present invention.
具体实施方式Detailed ways
本发明的核心是提供一种基于剪切波的成像方法、系统及装置，获取剪切波检测区域的剪切波传播反射系数，从而可根据剪切波传播反射系数准确确定剪切波检测区域的组织硬度均匀性，利于区域组织病变分析。The core of the present invention is to provide a shear wave-based imaging method, system and device to obtain the shear wave propagation reflection coefficient of the shear wave detection area, so that the shear wave detection area can be accurately determined according to the shear wave propagation reflection coefficient The uniformity of the tissue hardness is beneficial to the analysis of regional tissue lesions.
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。In order to make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be described clearly and completely in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments These are a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative work shall fall within the protection scope of the present invention.
请参照图1，图1为本发明实施例提供的一种基于剪切波的成像方法的流程图。Please refer to FIG. 1, which is a flowchart of a shear wave-based imaging method according to an embodiment of the present invention.
该基于剪切波的成像方法包括：The shear wave-based imaging method includes:
步骤S1：在剪切波检测区域内选取横向局部区域，并获取横向局部区域对应的剪切波的入射波能量和反射波能量。Step S1: Select a lateral local area in the shear wave detection area, and obtain the incident wave energy and the reflected wave energy of the shear wave corresponding to the lateral local area.
具体地，考虑到在剪切波的实际传播过程中，由于不同组织硬度的差异，会导致剪切波按照不同组织硬度差异的大小，发生不同程度的反射，所以剪切波在传播过程中发生的反射情况可以直接准确地反映出组织硬度的均匀性。更具体地，两个相邻组织的硬度差异越大，剪切波在两个组织交界处发生的反射程度越大，即在同一入射波能量的情况下，剪切波在硬度差异越大的两个组织交界处的反射波能量越大。所以，本实施例可采用反射波能量和入射波能量的比值，即剪切波传播反射系数表示剪切波在传播过程中发生的反射情况。因此，本实施例为了求取一局部组织区域的剪 切波传播反射系数，首先应获取该局部组织区域的剪切波的入射波能量和反射波能量。Specifically, considering that in the actual propagation process of shear waves, due to the difference in hardness of different tissues, the shear waves will be reflected to different degrees according to the difference in hardness of different tissues, so the shear waves occur during the propagation process. The reflection situation can directly and accurately reflect the uniformity of the tissue hardness. More specifically, the greater the difference in hardness between two adjacent tissues, the greater the degree of reflection of the shear wave at the junction of the two tissues, that is, under the same incident wave energy, the greater the difference in hardness of the shear wave is. The greater the energy of the reflected wave at the junction of the two tissues. Therefore, in this embodiment, the ratio of the reflected wave energy to the incident wave energy can be used, that is, the shear wave propagation reflection coefficient represents the reflection of the shear wave during the propagation process. Therefore, in this embodiment, in order to obtain the shear wave propagation reflection coefficient of a local tissue area, the incident wave energy and the reflected wave energy of the shear wave in the local tissue area should be obtained first.
基于此，可以理解的是，在剪切波检测区域内，综合剪切波检测区域的多个横向局部区域的剪切波传播反射系数，可以直接准确地确定出整个剪切波检测区域的组织硬度均匀性。更具体地，本实施例在剪切波检测区域内，首先选取一小段横向局部区域，然后获取所选取的横向局部区域对应的剪切波的入射波能量和反射波能量，目的是后续求取所选取的横向局部区域的剪切波传播反射系数。Based on this, it can be understood that in the shear wave detection area, combining the shear wave propagation reflection coefficients of multiple lateral local areas of the shear wave detection area can directly and accurately determine the organization of the entire shear wave detection area Uniformity of hardness. More specifically, in this embodiment, in the shear wave detection area, a small section of the lateral local area is first selected, and then the incident wave energy and the reflected wave energy of the shear wave corresponding to the selected lateral local area are obtained. The shear wave propagation reflection coefficient of the selected lateral local area.
步骤S2：根据反射波能量和入射波能量的比值，得到剪切波传播反射系数，并判断是否遍历完整个剪切波检测区域，若否，则返回执行在剪切波检测区域内选取横向局部区域的步骤；若是，则执行步骤S3。Step S2: According to the ratio of the reflected wave energy to the incident wave energy, obtain the shear wave propagation reflection coefficient, and determine whether to traverse the entire shear wave detection area, if not, return to the execution and select the lateral part in the shear wave detection area Step of the area; if yes, go to step S3.
具体地，本实施例在获取所选取的横向局部区域对应的剪切波的入射波能量和反射波能量后，计算反射波能量和入射波能量的比值，从而得到所选取的横向局部区域的剪切波传播反射系数。比如，横向局部区域的剪切波传播反射系数＝横向局部区域对应的反射波能量除以其对应入射波能量，则剪切波传播反射系数越大，说明剪切波在横向局部区域内的反射程度越大。Specifically, in this embodiment, after obtaining the incident wave energy and the reflected wave energy of the shear wave corresponding to the selected lateral local area, the ratio of the reflected wave energy to the incident wave energy is calculated, thereby obtaining the shear wave energy of the selected lateral local area. Cut-wave propagation reflection coefficient. For example, the shear wave propagation reflection coefficient of the lateral local area = the reflected wave energy corresponding to the lateral local area divided by its corresponding incident wave energy, the larger the shear wave propagation reflection coefficient, the reflection of the shear wave in the lateral local area The greater the degree.
为了确定出整个剪切波检测区域的组织硬度均匀性，本实施例需遍历完整个剪切波检测区域，如果未遍历完整个剪切波检测区域，则继续在未被遍历的剪切波检测区域内选择一段横向局部区域，并按照上述步骤求取此横向局部区域的剪切波传播反射系数；如果已遍历完整个剪切波检测区域，则可得到整个剪切波检测区域的剪切波传播反射系数。In order to determine the uniformity of the tissue hardness of the entire shear wave detection area, the present embodiment needs to traverse the entire shear wave detection area. If the entire shear wave detection area is not traversed, then continue to detect the shear wave that has not been traversed. Select a section of lateral local area in the area, and obtain the shear wave propagation reflection coefficient of this lateral local area according to the above steps; if the entire shear wave detection area has been traversed, the shear wave of the entire shear wave detection area can be obtained Propagation reflection coefficient.
需要说明的是，各横向局部区域的剪切波传播反射系数可能都是不一样的，每个横向局部区域的剪切波传播反射系数均由该横向局部区域的组织硬度决定，即具有组织弹性检测的唯一性和可参照性，从而更准确地检测出组织弹性情况。It should be noted that the shear wave propagation reflection coefficient of each lateral local area may be different, and the shear wave propagation reflection coefficient of each lateral local area is determined by the tissue hardness of the lateral local area, that is, it has tissue elasticity. The uniqueness and referability of the detection, so as to detect the elasticity of the tissue more accurately.
步骤S3：根据剪切波传播反射系数生成表征剪切波检测区域的组织硬度均匀性的组织图像。Step S3: Generate a tissue image characterizing the uniformity of the tissue hardness in the shear wave detection area according to the shear wave propagation reflection coefficient.
具体地，已知相邻两种组织硬度差异越大，剪切波传播到这两种组织边界时发生的反射越大，即剪切波传播反射系数与组织硬度均匀性具有一定对应关系。基于此，本实施例可根据整个剪切波检测区域的剪切波传播反射系数，生成清晰表征剪切波检测区域的组织硬度均匀性的组织图像，从而可通过分析组织图像及时发现病灶组织。Specifically, it is known that the greater the hardness difference between two adjacent tissues is, the greater the reflection that occurs when the shear wave propagates to the boundary of the two tissues, that is, the shear wave propagation reflection coefficient has a certain corresponding relationship with the uniformity of the tissue hardness. Based on this, this embodiment can generate a tissue image that clearly characterizes the uniformity of the tissue hardness of the shear wave detection area according to the shear wave propagation reflection coefficient of the entire shear wave detection area, so that the focus tissue can be found in time by analyzing the tissue image.
需要说明的是，剪切波可以是通过声辐射力产生，也可以是通过外部振动产生，也可以是身体内部自身运动组织产生。It should be noted that the shear wave can be generated by acoustic radiation force, by external vibration, or by the internal motion tissue of the body.
本发明实施例提供了一种基于剪切波的成像方法，考虑到在剪切波的实际传播过程中，由于不同组织硬度的差异，会导致剪切波按照不同组织硬度差异的大小，发生不同程度的反射，所以剪切波传播反射系数可直接准确地反映出组织硬度的均匀性。基于此，本申请获取剪切波检测区域的剪切波传播反射系数，从而可根据剪切波传播反射系数准确确定剪切波检测区域的组织硬度均匀性，利于区域组织病变分析。The embodiment of the present invention provides a shear wave-based imaging method. Considering that in the actual propagation process of the shear wave, due to the difference in the hardness of different tissues, the shear wave will be different according to the difference in the hardness of the different tissues. Degree of reflection, so the shear wave propagation reflection coefficient can directly and accurately reflect the uniformity of tissue hardness. Based on this, the present application obtains the shear wave propagation reflection coefficient of the shear wave detection area, so that the uniformity of the tissue hardness of the shear wave detection area can be accurately determined according to the shear wave propagation reflection coefficient, which is beneficial to the analysis of regional tissue lesions.
在上述实施例的基础上：On the basis of the above embodiment:
作为一种可选的实施例，在剪切波检测区域内选取横向局部区域，并获取横向局部区域对应的剪切波的入射波能量和反射波能量的过程，包括：As an optional embodiment, the process of selecting a lateral local area in the shear wave detection area and obtaining the incident wave energy and reflected wave energy of the shear wave corresponding to the lateral local area includes:
周期性获取剪切波检测区域内各检测点在剪切波的传播作用下的组织位移，并根据组织位移得到各检测点的一维时间位移波形；Periodically obtain the tissue displacement of each detection point in the shear wave detection area under the action of shear wave propagation, and obtain the one-dimensional time displacement waveform of each detection point according to the tissue displacement;
在剪切波检测区域内选取横向局部区域，并将横向局部区域内各检测点的一维时间位移波形进行组合，得到横向局部区域对应的二维时间空间传播矩阵；Select a lateral local area in the shear wave detection area, and combine the one-dimensional time displacement waveforms of each detection point in the lateral local area to obtain a two-dimensional time-space propagation matrix corresponding to the lateral local area;
根据二维时间空间传播矩阵得到剪切波的正向传播信息矩阵和反向传播信息矩阵，并根据正向传播信息矩阵和反向传播信息矩阵相应得到剪切波的入射波能量和反射波能量。Obtain the forward propagation information matrix and the backward propagation information matrix of the shear wave according to the two-dimensional time and space propagation matrix, and obtain the incident wave energy and the reflected wave energy of the shear wave according to the forward propagation information matrix and the back propagation information matrix. .
具体地，当剪切波传播到某一位置的组织时，此位置的组织会发生位移，所以剪切波检测区域内各检测点的组织位移情况可以体现出剪切波检测区域内剪切波的能量传播情况。因此，本实施例为了了解剪切波检测区域内剪切波的能量传播情况，周期性获取剪切波检测区域内各检测点在剪 切波的传播作用下的组织位移，并根据各检测点的组织位移一一得到各检测点的一维时间位移波形。Specifically, when the shear wave propagates to the tissue at a certain location, the tissue at this location will be displaced. Therefore, the tissue displacement of each detection point in the shear wave detection area can reflect the shear wave in the shear wave detection area. The energy transmission situation. Therefore, in this embodiment, in order to understand the energy propagation of the shear wave in the shear wave detection area, periodically obtain the tissue displacement of each detection point in the shear wave detection area under the action of the shear wave propagation, and according to each detection point The one-dimensional time-displacement waveform of each detection point is obtained one by one of the tissue displacement.
整个剪切波检测区域内剪切波的能量传播情况具体是综合其内多个横向局部区域内剪切波的能量传播情况得到，在一具体实施例中：在剪切波检测区域内选取横向局部区域后，将所选取的横向局部区域内各检测点的一维时间位移波形进行组合，得到横向局部区域对应的二维时间空间传播矩阵。由于横向局部区域对应的二维时间空间传播矩阵包含了横向局部区域的剪切波的正向传播信息矩阵和反向传播信息矩阵，所以本实施例可根据横向局部区域对应的二维时间空间传播矩阵，得到横向局部区域的剪切波的正向传播信息矩阵和反向传播信息矩阵，并根据横向局部区域的剪切波的正向传播信息矩阵得到剪切波在横向局部区域内的入射波能量，同理根据横向局部区域的剪切波的反向传播信息矩阵得到剪切波在横向局部区域内的反射波能量。然后以此技术手段遍历整个剪切波检测区域，从而得到整个剪切波检测区域内剪切波的能量传播情况。The energy propagation of the shear wave in the entire shear wave detection area is specifically obtained by integrating the energy propagation of the shear wave in multiple lateral local areas within it. In a specific embodiment: the transverse direction is selected in the shear wave detection area. After the local area, the one-dimensional time displacement waveforms of the detection points in the selected lateral local area are combined to obtain the two-dimensional time-space propagation matrix corresponding to the lateral local area. Since the two-dimensional time-space propagation matrix corresponding to the lateral local area contains the forward propagation information matrix and the back-propagation information matrix of the shear wave in the lateral local area, this embodiment can be based on the two-dimensional time and space propagation corresponding to the lateral local area. Matrix, obtain the forward propagation information matrix and the reverse propagation information matrix of the shear wave in the lateral local area, and obtain the incident wave of the shear wave in the lateral local area according to the forward propagation information matrix of the shear wave in the lateral local area Energy, in the same way, the reflected wave energy of the shear wave in the lateral local area is obtained according to the back propagation information matrix of the shear wave in the lateral local area. Then use this technical means to traverse the entire shear wave detection area, so as to obtain the energy propagation of the shear wave in the entire shear wave detection area.
此外，本实施例分析图2的与二维时间空间传播矩阵对应的剪切波传播波形图(剪切波传播波形上各位置的数据组成二维时间空间传播矩阵)，可得剪切波反射发生位置，剪切波反向传播的速度、频率及频散等参数。In addition, this embodiment analyzes the shear wave propagation waveform diagram corresponding to the two-dimensional time-space propagation matrix in Fig. 2 (the data at each position on the shear wave propagation waveform forms a two-dimensional time-space propagation matrix), and the shear wave reflection can be obtained. The location, the speed, frequency and dispersion of the shear wave back propagation.
作为一种可选的实施例，周期性获取剪切波检测区域内各检测点在剪切波的传播作用下的组织位移的过程，包括：As an optional embodiment, the process of periodically acquiring the tissue displacement of each detection point in the shear wave detection area under the action of the shear wave propagation includes:
向剪切波检测区域周期性发射用于检测剪切波信号的检测波束，并将返回的检测回波信号采用多波束技术得到在剪切波检测区域内各检测点的位置信息；Periodically transmit detection beams for detecting shear wave signals to the shear wave detection area, and use multi-beam technology to obtain the position information of each detection point in the shear wave detection area with the returned detection echo signal;
根据周期性得到的位置信息得到各检测点在剪切波的传播作用下的组织位移。According to the periodically obtained position information, the tissue displacement of each detection point under the propagation of shear waves is obtained.
具体地，本实施例获取剪切波检测区域内各检测点在剪切波的传播作用下的组织位移的技术手段是：首先按照预设重复频率，重复向剪切波检测区域发射用于检测剪切波信号的检测波束，即向剪切波检测区域周期性发射检测波束。然后接收检测回波信号，以采用多波束技术(一种测定目标坐标的雷达总体技术)实现根据返回的检测回波信号得到在剪切波检测 区域内各检测点的位置信息。可以理解的是，比较本周期得到的一检测点的位置信息与上一周期得到的此检测点的位置信息，可知此检测点在这一时间段的组织位移，即根据各检测点在周期性时间的位置信息，得到各检测点在剪切波的传播作用下随时间变化的组织位移。Specifically, the technical means of this embodiment to obtain the tissue displacement of each detection point in the shear wave detection area under the action of the shear wave propagation is: first, according to the preset repetition frequency, repeatedly transmit to the shear wave detection area for detection. The detection beam of the shear wave signal, that is, the detection beam is periodically transmitted to the shear wave detection area. Then the detection echo signal is received, and the multi-beam technology (a kind of radar overall technology for measuring target coordinates) is used to obtain the position information of each detection point in the shear wave detection area according to the returned detection echo signal. It is understandable that by comparing the position information of a detection point obtained in this cycle with the position information of the detection point obtained in the previous cycle, it can be known that the tissue displacement of the detection point in this period of time is based on the periodicity of each detection point. Time position information, the tissue displacement of each detection point that changes with time under the action of shear wave propagation is obtained.
作为一种可选的实施例，根据二维时间空间传播矩阵得到剪切波的正向传播信息矩阵和反向传播信息矩阵的过程，包括：As an optional embodiment, the process of obtaining the forward propagation information matrix and the backward propagation information matrix of the shear wave according to the two-dimensional time-space propagation matrix includes:
对二维时间空间传播矩阵进行2D傅里叶变换，得到包括剪切波的正向传播信息矩阵和反向传播信息矩阵的二维wk频域矩阵；Perform a 2D Fourier transform on the two-dimensional time-space propagation matrix to obtain a two-dimensional wk frequency domain matrix including the forward propagation information matrix of the shear wave and the backward propagation information matrix;
从二维wk频域矩阵中分别提取出正向传播信息矩阵和反向传播信息矩阵。The forward propagation information matrix and the backward propagation information matrix are respectively extracted from the two-dimensional wk frequency domain matrix.
具体地，本实施例在得到横向局部区域对应的二维时间空间传播矩阵后，对二维时间空间传播矩阵经方向滤波器进行处理，即对二维时间空间传播矩阵进行2D傅里叶变换，得到二维wk频域矩阵，即二维时间频率-空间频率矩阵，如图3所示，二维wk频域矩阵的第二象限和第四象限的信息矩阵AMat表示剪切波正向传播的信息，第一象限和第三象限的信息矩阵BMat表示剪切波反向传播的信息，所以本实施例可从二维wk频域矩阵中分别提取出正向传播信息矩阵AMat(如图4)和反向传播信息矩阵BMat(如图5)。Specifically, in this embodiment, after obtaining the two-dimensional time-space propagation matrix corresponding to the lateral local area, the two-dimensional time-space propagation matrix is processed by the directional filter, that is, the two-dimensional time-space propagation matrix is subjected to 2D Fourier transform. Obtain the two-dimensional wk frequency domain matrix, that is, the two-dimensional time frequency-spatial frequency matrix, as shown in Figure 3. The information matrix AMat of the second and fourth quadrants of the two-dimensional wk frequency domain matrix represents the forward propagation of the shear wave. Information, the information matrix BMat of the first quadrant and the third quadrant represents the information of the shear wave back propagation, so this embodiment can extract the forward propagation information matrix AMat from the two-dimensional wk frequency domain matrix respectively (Figure 4) And the back-propagation information matrix BMat (Figure 5).
更具体地，本实施例利用不同掩膜分别提取二维wk频域矩阵中的正向传播信息矩阵AMat和反向传播信息矩阵BMat，其中，掩膜的取值：提取正向传播信息矩阵AMat时，一三象限取值为0，二四象限取值为1，则掩膜与二维wk频域矩阵相乘得到正向传播信息矩阵AMat；提取反向传播信息矩阵BMat时，一三象限取值为1，二四象限取值为0，则掩膜与二维wk频域矩阵相乘得到反向传播信息矩阵BMat。More specifically, this embodiment uses different masks to respectively extract the forward propagation information matrix AMat and the back propagation information matrix BMat in the two-dimensional wk frequency domain matrix, where the value of the mask: extract the forward propagation information matrix AMat When the first and third quadrants are set to 0 and the second and fourth quadrants are set to 1, then the mask is multiplied by the two-dimensional wk frequency domain matrix to obtain the forward propagation information matrix AMat; when the back propagation information matrix BMat is extracted, the first and third quadrants The value is 1 and the value of the second and fourth quadrants is 0, then the mask is multiplied by the two-dimensional wk frequency domain matrix to obtain the back propagation information matrix BMat.
作为一种可选的实施例，根据正向传播信息矩阵和反向传播信息矩阵相应得到剪切波的入射波能量和反射波能量的过程，包括：As an optional embodiment, the process of correspondingly obtaining the incident wave energy and the reflected wave energy of the shear wave according to the forward propagation information matrix and the backward propagation information matrix includes:
分别对正向传播信息矩阵和反向传播信息矩阵取绝对值；Take the absolute value of the forward propagation information matrix and the backward propagation information matrix respectively;
对正向传播信息矩阵的绝对值进行时间频率和空间频率的二维积分，得到剪切波的入射波能量；Perform two-dimensional integration of time frequency and space frequency on the absolute value of the forward propagation information matrix to obtain the incident wave energy of the shear wave;
对反向传播信息矩阵的绝对值进行时间频率和空间频率的二维积分，得到剪切波的反射波能量。The two-dimensional integration of time frequency and space frequency is performed on the absolute value of the backpropagation information matrix to obtain the reflected wave energy of the shear wave.
具体地，根据正向传播信息矩阵AMat求取剪切波的入射波能量Apower的计算方法为：
Figure PCTCN2020096261-appb-000001
同样的，根据反向传播信息矩阵BMat求取剪切波的反射波能量Bpower的计算方法为：
Figure PCTCN2020096261-appb-000002
Specifically, the calculation method for obtaining the incident wave energy Apower of the shear wave according to the forward propagation information matrix AMat is:
Figure PCTCN2020096261-appb-000001
Similarly, the calculation method for obtaining the reflected wave energy Bpower of the shear wave according to the backpropagation information matrix BMat is:
Figure PCTCN2020096261-appb-000002
作为一种可选的实施例，根据剪切波传播反射系数生成表征剪切波检测区域的组织硬度均匀性的组织图像的过程，包括：As an optional embodiment, the process of generating a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area according to the shear wave propagation reflection coefficient includes:
将剪切波检测区域的剪切波传播反射系数均归一化到预设范围内，得到剪切波传播的反射系数归一值；Normalize the shear wave propagation reflection coefficients in the shear wave detection area to a preset range to obtain the normalized value of the shear wave propagation reflection coefficient;
根据预设系数颜色对应关系确定与反射系数归一值对应的颜色；Determine the color corresponding to the normalized value of the reflection coefficient according to the preset coefficient color correspondence relationship;
根据颜色相应填充剪切波检测区域的颜色，并控制显示屏显示填充有相应颜色的剪切波检测区域。Fill the color of the shear wave detection area according to the color, and control the display screen to display the shear wave detection area filled with the corresponding color.
需要说明的是，本实施例的预设是提前设置好的，只需要设置一次，除非根据实际情况需要修改，否则不需要重新设置。It should be noted that the preset in this embodiment is set in advance, and only needs to be set once, unless it needs to be modified according to the actual situation, it does not need to be reset.
进一步地，本实施例提前设定一个范围：[min，max]，以将剪切波检测区域的各剪切波传播反射系数均按照一定规则归一化到所设定的范围内，从而一一得到剪切波传播的反射系数归一值。Further, in this embodiment, a range is set in advance: [min, max], so that the shear wave propagation reflection coefficients in the shear wave detection area are all normalized to the set range according to a certain rule, so that one One gets the normalized value of the reflection coefficient of shear wave propagation.
本实施例还提前设定剪切波检测区域的反射系数归一值与颜色之间的对应关系，简称系数颜色对应关系，如最大值max对应颜色等级255，最小值min对应颜色等级0，则反射系数归一值RflctRatio对应的颜色等级Rflct为：
Figure PCTCN2020096261-appb-000003
从而确定反射系数归一值对应的颜色。基于此，本实施例在得到各反射系数归一值后，根据系数颜色对应关系确定与各反射系数归一值一一对应的各颜色，然后根据各颜色相应填充剪切波 检测区域的颜色，并控制显示屏显示填充有相应颜色的剪切波检测区域，供工作人员查看，从而更直观地展示出剪切波检测区域的组织硬度均匀性。 In this embodiment, the correspondence relationship between the normalized value of the reflection coefficient and the color of the shear wave detection area is also set in advance, referred to as the coefficient color correspondence relationship. For example, the maximum value max corresponds to the color level 255, and the minimum value min corresponds to the color level 0, then The color level Rflct corresponding to the normalized value of the reflection coefficient RflctRatio is:
Figure PCTCN2020096261-appb-000003
In order to determine the color corresponding to the normalized value of the reflection coefficient. Based on this, in this embodiment, after obtaining the normalized values of the reflection coefficients, the colors corresponding to the normalized values of the reflection coefficients are determined one-to-one according to the coefficient color correspondence relationship, and then the colors of the shear wave detection area are filled correspondingly according to each color. And control the display screen to display the shear wave detection area filled with the corresponding color for the staff to view, so as to more intuitively show the uniformity of the tissue hardness of the shear wave detection area.
需要说明的是，本实施例可以由不同灰度值表示不同颜色等级；也可以由不同颜色表示不同颜色等级。It should be noted that in this embodiment, different gray values may represent different color levels; different colors may also represent different color levels.
作为一种可选的实施例，成像方法还包括：As an optional embodiment, the imaging method further includes:
在剪切波正向传播方向上，将相邻的横向局部区域对应的入射波能量相减，得到剪切波传播衰减量；In the forward propagation direction of the shear wave, subtract the incident wave energy corresponding to the adjacent lateral local area to obtain the shear wave propagation attenuation;
根据相邻的两横向局部区域对应的反射波能量修正剪切波传播衰减量。The shear wave propagation attenuation is corrected according to the reflected wave energy corresponding to two adjacent lateral local areas.
进一步地，考虑到剪切波在组织传播过程中，剪切波实际的传播衰减量并不包含反射造成的剪切波衰减量，所以本实施例计算剪切波传播衰减量的过程包括：在剪切波正向传播方向上，将相邻的横向局部区域A、B各自对应的入射波能量相减，得到剪切波传播衰减量，即剪切波传播衰减量＝|横向局部区域A对应的入射波能量-横向局部区域B对应的入射波能量|；将剪切波传播衰减量的绝对值依次减去相邻的横向局部区域A、B各自对应的反射波能量，得到修正后的剪切波传播衰减量，即修正后的剪切波传播衰减量＝|横向局部区域A对应的入射波能量-横向局部区域B对应的入射波能量|-横向局部区域A对应的反射波能量-横向局部区域B对应的反射波能量，从而更加准确地计算出剪切波传播过程中的衰减。由于剪切波传播衰减量与组织粘性等物理参数具有一定相关性，所以本实施例同样能够更加准确地得到组织粘性等相关物理参数。Furthermore, considering that the shear wave propagation attenuation in the tissue propagation process does not include the shear wave attenuation caused by reflection, the process of calculating the shear wave propagation attenuation in this embodiment includes: In the forward propagation direction of the shear wave, subtract the respective incident wave energies of the adjacent lateral local areas A and B to obtain the shear wave propagation attenuation, that is, the shear wave propagation attenuation=|corresponding to the lateral local area A The incident wave energy-the incident wave energy corresponding to the transverse local area B|; the absolute value of the shear wave propagation attenuation is sequentially subtracted from the reflected wave energy of the adjacent transverse local areas A and B to obtain the corrected shear Shear wave propagation attenuation, that is, the corrected shear wave propagation attenuation=|incident wave energy corresponding to transverse local area A-incident wave energy corresponding to transverse local area B|-reflected wave energy corresponding to transverse local area A-transverse The reflected wave energy corresponding to the local area B, so as to more accurately calculate the attenuation during the propagation of the shear wave. Since the attenuation of shear wave propagation has a certain correlation with physical parameters such as tissue viscosity, this embodiment can also more accurately obtain related physical parameters such as tissue viscosity.
请参照图6，图6为本发明实施例提供的一种基于剪切波的成像系统的结构示意图。Please refer to FIG. 6, which is a schematic structural diagram of a shear wave-based imaging system provided by an embodiment of the present invention.
该基于剪切波的成像系统包括：The shear wave-based imaging system includes:
能量获取模块1，用于在剪切波检测区域内选取横向局部区域，并获取横向局部区域对应的剪切波的入射波能量和反射波能量；The energy acquisition module 1 is used to select a lateral local area in the shear wave detection area, and obtain the incident wave energy and the reflected wave energy of the shear wave corresponding to the lateral local area;
系数求取模块2，用于根据反射波能量和入射波能量的比值，得到剪切波传播反射系数，并返回执行能量获取模块，直至遍历整个剪切波检测区域，以得到整个剪切波检测区域的剪切波传播反射系数；The coefficient calculation module 2 is used to obtain the shear wave propagation reflection coefficient according to the ratio of the reflected wave energy to the incident wave energy, and return to the execution of the energy acquisition module until the entire shear wave detection area is traversed to obtain the entire shear wave detection The shear wave propagation reflection coefficient of the area;
组织成像模块3，用于根据剪切波传播反射系数生成表征剪切波检测区域的组织硬度均匀性的组织图像。The tissue imaging module 3 is used to generate a tissue image characterizing the uniformity of the tissue hardness in the shear wave detection area according to the shear wave propagation reflection coefficient.
本申请提供的分析系统的介绍请参考上述分析方法的实施例，本申请在此不再赘述。For the introduction of the analysis system provided in this application, please refer to the embodiment of the above analysis method, which will not be repeated in this application.
本申请实施例还提供了一种基于剪切波的成像装置，包括：An embodiment of the present application also provides a shear wave-based imaging device, including:
存储器，用于存储计算机程序；Memory, used to store computer programs;
处理器，用于在执行计算机程序时实现上述任一种基于剪切波的成像方法的步骤。The processor is used to implement the steps of any of the foregoing shear wave-based imaging methods when the computer program is executed.
本申请提供的分析装置的介绍请参考上述分析方法的实施例，本申请在此不再赘述。For the introduction of the analysis device provided in this application, please refer to the embodiment of the above analysis method, which will not be repeated in this application.
还需要说明的是，在本说明书中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。It should also be noted that in this specification, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply these entities or operations. There is any such actual relationship or sequence between operations. Moreover, the terms "include", "include" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements not only includes those elements, but also includes those that are not explicitly listed Other elements of, or also include elements inherent to this process, method, article or equipment. If there are no more restrictions, the element defined by the sentence "including a..." does not exclude the existence of other identical elements in the process, method, article, or equipment that includes the element.
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其他实施例中实现。因此，本发明将不会被限制于本文所示的 这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。The above description of the disclosed embodiments enables those skilled in the art to implement or use the present invention. Various modifications to these embodiments will be obvious to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown in this document, but should conform to the widest scope consistent with the principles and novel features disclosed in this document.

Claims (9)

  1. 一种基于剪切波的成像方法，其特征在于，包括：A shear wave-based imaging method, which is characterized in that it includes:
    在剪切波检测区域内选取横向局部区域，并获取所述横向局部区域对应的剪切波的入射波能量和反射波能量；Selecting a lateral local area in the shear wave detection area, and obtaining the incident wave energy and the reflected wave energy of the shear wave corresponding to the lateral local area;
    根据所述反射波能量和所述入射波能量的比值，得到剪切波传播反射系数，并返回在所述剪切波检测区域内选取横向局部区域的步骤，直至遍历整个剪切波检测区域，以得到整个剪切波检测区域的剪切波传播反射系数；According to the ratio of the reflected wave energy to the incident wave energy, the shear wave propagation reflection coefficient is obtained, and the step of selecting a lateral local area in the shear wave detection area is returned until the entire shear wave detection area is traversed, To obtain the shear wave propagation reflection coefficient of the entire shear wave detection area;
    根据所述剪切波传播反射系数生成表征所述剪切波检测区域的组织硬度均匀性的组织图像。According to the shear wave propagation reflection coefficient, a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area is generated.
  2. 如权利要求1所述的基于剪切波的成像方法，其特征在于，所述在剪切波检测区域内选取横向局部区域，并获取所述横向局部区域对应的剪切波的入射波能量和反射波能量的过程，包括：The shear wave-based imaging method according to claim 1, wherein the transverse local area is selected in the shear wave detection area, and the incident wave energy and the shear wave corresponding to the transverse local area are obtained. The process of reflecting wave energy includes:
    周期性获取剪切波检测区域内各检测点在剪切波的传播作用下的组织位移，并根据所述组织位移得到所述各检测点的一维时间位移波形；Periodically acquiring the tissue displacement of each detection point in the shear wave detection area under the propagation of the shear wave, and obtaining a one-dimensional time displacement waveform of each detection point according to the tissue displacement;
    在所述剪切波检测区域内选取横向局部区域，并将所述横向局部区域内各检测点的一维时间位移波形进行组合，得到所述横向局部区域对应的二维时间空间传播矩阵；Selecting a lateral local area in the shear wave detection area, and combining the one-dimensional time displacement waveforms of each detection point in the lateral local area to obtain a two-dimensional time-space propagation matrix corresponding to the lateral local area;
    根据所述二维时间空间传播矩阵得到剪切波的正向传播信息矩阵和反向传播信息矩阵，并根据所述正向传播信息矩阵和所述反向传播信息矩阵相应得到剪切波的入射波能量和反射波能量。Obtain the forward propagation information matrix and the backward propagation information matrix of the shear wave according to the two-dimensional time and space propagation matrix, and obtain the incident shear wave according to the forward propagation information matrix and the backward propagation information matrix. Wave energy and reflected wave energy.
  3. 如权利要求2所述的基于剪切波的成像方法，其特征在于，所述周期性获取剪切波检测区域内各检测点在剪切波的传播作用下的组织位移的过程，包括：The shear wave-based imaging method according to claim 2, wherein the process of periodically acquiring the tissue displacement of each detection point in the shear wave detection area under the propagation of the shear wave comprises:
    向剪切波检测区域周期性发射用于检测剪切波信号的检测波束，并将返回的检测回波信号采用多波束技术得到在所述剪切波检测区域内各检测点的位置信息；Periodically transmit a detection beam for detecting the shear wave signal to the shear wave detection area, and use multi-beam technology to obtain the position information of each detection point in the shear wave detection area for the returned detection echo signal;
    根据周期性得到的位置信息得到所述各检测点在剪切波的传播作用下的组织位移。The tissue displacement of each detection point under the propagation action of the shear wave is obtained according to the periodically obtained position information.
  4. 如权利要求2所述的基于剪切波的成像方法，其特征在于，所述根据所述二维时间空间传播矩阵得到剪切波的正向传播信息矩阵和反向传播信息矩阵的过程，包括：The shear wave-based imaging method according to claim 2, wherein the process of obtaining the forward propagation information matrix and the backward propagation information matrix of the shear wave according to the two-dimensional time-space propagation matrix includes :
    对所述二维时间空间传播矩阵进行2D傅里叶变换，得到包括剪切波的正向传播信息矩阵和反向传播信息矩阵的二维wk频域矩阵；Performing a 2D Fourier transform on the two-dimensional time-space propagation matrix to obtain a two-dimensional wk frequency domain matrix including a forward propagation information matrix of the shear wave and a backward propagation information matrix;
    从所述二维wk频域矩阵中分别提取出所述正向传播信息矩阵和所述反向传播信息矩阵。The forward propagation information matrix and the backward propagation information matrix are respectively extracted from the two-dimensional wk frequency domain matrix.
  5. 如权利要求4所述的基于剪切波的成像方法，其特征在于，所述根据所述正向传播信息矩阵和所述反向传播信息矩阵相应得到剪切波的入射波能量和反射波能量的过程，包括：The shear wave-based imaging method according to claim 4, wherein the incident wave energy and the reflected wave energy of the shear wave are obtained according to the forward propagation information matrix and the backward propagation information matrix. The process includes:
    分别对所述正向传播信息矩阵和所述反向传播信息矩阵取绝对值；Taking absolute values for the forward propagation information matrix and the backward propagation information matrix respectively;
    对所述正向传播信息矩阵的绝对值进行时间频率和空间频率的二维积分，得到剪切波的入射波能量；Performing two-dimensional integration of time frequency and space frequency on the absolute value of the forward propagation information matrix to obtain the incident wave energy of the shear wave;
    对所述反向传播信息矩阵的绝对值进行时间频率和空间频率的二维积分，得到剪切波的反射波能量。The two-dimensional integration of the time frequency and the space frequency is performed on the absolute value of the back propagation information matrix to obtain the reflected wave energy of the shear wave.
  6. 如权利要求1-5任一项所述的基于剪切波的成像方法，其特征在于，所述根据所述剪切波传播反射系数生成表征所述剪切波检测区域的组织硬度均匀性的组织图像的过程，包括：The shear wave-based imaging method according to any one of claims 1 to 5, wherein the shear wave propagation reflection coefficient is used to generate a signal that characterizes the uniformity of the tissue hardness of the shear wave detection area. The process of organizing images includes:
    将剪切波检测区域的剪切波传播反射系数均归一化到预设范围内，得到剪切波传播的反射系数归一值；Normalize the shear wave propagation reflection coefficients in the shear wave detection area to a preset range to obtain the normalized value of the shear wave propagation reflection coefficient;
    根据预设系数颜色对应关系确定与所述反射系数归一值对应的颜色；Determining the color corresponding to the normalized value of the reflection coefficient according to the preset coefficient color correspondence relationship;
    根据所述颜色相应填充所述剪切波检测区域的颜色，并控制显示屏显示填充有相应颜色的剪切波检测区域。Fill the color of the shear wave detection area correspondingly according to the color, and control the display screen to display the shear wave detection area filled with the corresponding color.
  7. 如权利要求1所述的基于剪切波的成像方法，其特征在于，所述成像方法还包括：The shear wave-based imaging method of claim 1, wherein the imaging method further comprises:
    在剪切波正向传播方向上，将相邻的横向局部区域对应的入射波能量相减，得到剪切波传播衰减量；In the forward propagation direction of the shear wave, subtract the incident wave energy corresponding to the adjacent lateral local area to obtain the shear wave propagation attenuation;
    根据相邻的两个所述横向局部区域对应的反射波能量修正所述剪切波传播衰减量。The shear wave propagation attenuation is corrected according to the reflected wave energy corresponding to the two adjacent lateral local regions.
  8. 一种基于剪切波的成像系统，其特征在于，包括：A shear wave-based imaging system, which is characterized in that it comprises:
    能量获取模块，用于在剪切波检测区域内选取横向局部区域，并获取所述横向局部区域对应的剪切波的入射波能量和反射波能量；An energy acquisition module, configured to select a lateral local area in the shear wave detection area, and obtain the incident wave energy and the reflected wave energy of the shear wave corresponding to the lateral local area;
    系数求取模块，用于根据所述反射波能量和所述入射波能量的比值，得到剪切波传播反射系数，并返回执行所述能量获取模块，直至遍历整个剪切波检测区域，以得到整个剪切波检测区域的剪切波传播反射系数；The coefficient obtaining module is used to obtain the shear wave propagation reflection coefficient according to the ratio of the reflected wave energy to the incident wave energy, and return to execute the energy obtaining module until the entire shear wave detection area is traversed to obtain The shear wave propagation reflection coefficient of the entire shear wave detection area;
    组织成像模块，用于根据所述剪切波传播反射系数生成表征所述剪切波检测区域的组织硬度均匀性的组织图像。The tissue imaging module is configured to generate a tissue image characterizing the uniformity of the tissue hardness of the shear wave detection area according to the shear wave propagation reflection coefficient.
  9. 一种基于剪切波的成像装置，其特征在于，包括：An imaging device based on shear wave, which is characterized in that it comprises:
    存储器，用于存储计算机程序；Memory, used to store computer programs;
    处理器，用于在执行所述计算机程序时实现如权利要求1-7任一项所述的基于剪切波的成像方法的步骤。The processor is configured to implement the steps of the shear wave-based imaging method according to any one of claims 1-7 when the computer program is executed.
PCT/CN2020/096261 2019-11-27 2020-06-16 Shear wave-based imaging method, system and apparatus WO2021103493A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201911184432.1A CN110840488B (en) 2019-11-27 2019-11-27 Imaging method, system and device based on shear wave
CN201911184432.1 2019-11-27

Publications (1)

Publication Number Publication Date
WO2021103493A1 true WO2021103493A1 (en) 2021-06-03

Family

ID=69605458

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/096261 WO2021103493A1 (en) 2019-11-27 2020-06-16 Shear wave-based imaging method, system and apparatus

Country Status (2)

Country Link
CN (1) CN110840488B (en)
WO (1) WO2021103493A1 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110840488B (en) * 2019-11-27 2021-09-28 深圳开立生物医疗科技股份有限公司 Imaging method, system and device based on shear wave
CN111449681B (en) * 2020-04-08 2023-09-08 深圳开立生物医疗科技股份有限公司 Shear wave imaging method, device, equipment and readable storage medium
CN113243889B (en) * 2020-08-10 2022-05-10 北京航空航天大学 Method and apparatus for obtaining information of biological tissue
CN114631848B (en) * 2022-02-14 2024-04-05 逸超医疗科技(北京)有限公司 Method, device, equipment and storage medium for detecting tissue uniformity

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103126721A (en) * 2013-03-05 2013-06-05 飞依诺科技(苏州)有限公司 Method and device for compound ultrasonic imaging based on multiple focuses and multiple ultrasonic velocities
CN104302233A (en) * 2013-03-05 2015-01-21 日立阿洛卡医疗株式会社 Ultrasonic diagnosis device and transmission/reception method
US20150073275A1 (en) * 2013-09-11 2015-03-12 Kabushiki Kaisha Toshiba Ultrasound diagnosis apparatus and ultrasound imaging method
CN105748100A (en) * 2014-12-19 2016-07-13 深圳开立生物医疗科技股份有限公司 Quasi-static ultrasonic elastography displacement calculation method and device
CN105899140A (en) * 2013-08-12 2016-08-24 三星电子株式会社 Method for producing elastic image and ultrasonic diagnostic apparatus
CN105997147A (en) * 2016-07-28 2016-10-12 飞依诺科技(苏州)有限公司 Ultrasonic pulse Doppler imaging method and device
CN109310401A (en) * 2016-06-10 2019-02-05 皇家飞利浦有限公司 The damage growth in heating ablation is monitored using reflection shear-wave
CN110840488A (en) * 2019-11-27 2020-02-28 深圳开立生物医疗科技股份有限公司 Imaging method, system and device based on shear wave

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7611465B2 (en) * 2003-07-15 2009-11-03 Board Of Regents, The University Of Texas System Rapid and accurate detection of bone quality using ultrasound critical angle reflectometry
US7785259B2 (en) * 2003-10-03 2010-08-31 Mayo Foundation For Medical Education And Research Detection of motion in vibro-acoustography
CN102283680B (en) * 2011-05-18 2013-03-06 中国医学科学院生物医学工程研究所 Cataract crystalline lens hardness ultrasonic measurement system
CN103598897A (en) * 2013-11-11 2014-02-26 南京大学 Nanomechanics based ultrasonic measurement and characterization method for pathological changes of skin tissues
CN105300856B (en) * 2015-11-11 2017-09-29 上海理工大学 Based on ultrasonic impedances spectrum to granule density and the measuring method of size

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103126721A (en) * 2013-03-05 2013-06-05 飞依诺科技(苏州)有限公司 Method and device for compound ultrasonic imaging based on multiple focuses and multiple ultrasonic velocities
CN104302233A (en) * 2013-03-05 2015-01-21 日立阿洛卡医疗株式会社 Ultrasonic diagnosis device and transmission/reception method
CN105899140A (en) * 2013-08-12 2016-08-24 三星电子株式会社 Method for producing elastic image and ultrasonic diagnostic apparatus
US20150073275A1 (en) * 2013-09-11 2015-03-12 Kabushiki Kaisha Toshiba Ultrasound diagnosis apparatus and ultrasound imaging method
CN105748100A (en) * 2014-12-19 2016-07-13 深圳开立生物医疗科技股份有限公司 Quasi-static ultrasonic elastography displacement calculation method and device
CN109310401A (en) * 2016-06-10 2019-02-05 皇家飞利浦有限公司 The damage growth in heating ablation is monitored using reflection shear-wave
CN105997147A (en) * 2016-07-28 2016-10-12 飞依诺科技(苏州)有限公司 Ultrasonic pulse Doppler imaging method and device
CN110840488A (en) * 2019-11-27 2020-02-28 深圳开立生物医疗科技股份有限公司 Imaging method, system and device based on shear wave

Also Published As

Publication number Publication date
CN110840488B (en) 2021-09-28
CN110840488A (en) 2020-02-28

Similar Documents

Publication Publication Date Title
WO2021103493A1 (en) Shear wave-based imaging method, system and apparatus
JP6053860B2 (en) Computer program for determining shear wave information in medical ultrasound imaging systems
JP6230827B2 (en) Ultrasonic acoustic absorption measurement method, ultrasonic acoustic absorption measurement system, and computer-readable non-convertible storage medium for measuring ultrasonic acoustic attenuation
JP6140953B2 (en) Classification preprocessing in medical ultrasonic shear wave imaging
US9468421B2 (en) Visualization of associated information in ultrasound shear wave imaging
US9907539B2 (en) Sparse tracking in acoustic radiation force impulse imaging
US20060079777A1 (en) Ultrasonic image boundary extracting method, ultrasonic image boundary extracting apparatus, and ultrasonic imaging apparatus
US8801614B2 (en) On-axis shear wave characterization with ultrasound
US20060079780A1 (en) Ultrasonic imaging apparatus
JP2012081269A5 (en)
CN110974296B (en) Method and device for acquiring shear wave elastic image and ultrasonic imaging system
US20060079776A1 (en) Ultrasonic imaging apparatus
US20130343627A1 (en) Suppression of reverberations and/or clutter in ultrasonic imaging systems
US11766297B2 (en) Apparatus and method for detecting an interventional tool
US20180296190A1 (en) Ultrasonic diagnostic device and ultrasonic signal processing method
JP2016512130A (en) System and method for detecting and presenting interventional devices via ultrasound imaging
KR20160140858A (en) Ultrasound imaging system and method for tracking a specular reflector
KR20170094521A (en) Sound speed imaging using shear waves
US20200337679A1 (en) Ultrasonic signal processing apparatus, ultrasonic diagnostic apparatus, ultrasonic signal processing method, and ultrasonic signal processing program
KR102319397B1 (en) Angles for ultrasound-based shear wave imaging
US20190328363A1 (en) Ultrasound diagnostic apparatus and ultrasound signal processing method
US20210128100A1 (en) Analyzing apparatus and ultrasound diagnosis apparatus
US20240138816A1 (en) Ultrasound imaging method and ultrasound imaging system
US20220378403A1 (en) Ultrasound diagnostic apparatus and diagnosis assistance method
JP7302651B2 (en) Ultrasonic signal processing device, ultrasonic diagnostic device, ultrasonic signal processing method, and program

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20893697

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20893697

Country of ref document: EP

Kind code of ref document: A1