Arterial plaque characteristic detection method, device and storage medium based on elastic imaging

Abstract

The present invention provides a method, device, and storable medium for detecting arterial plaque characteristics based on elastic imaging. The method comprises: acquiring a tissue region of interest containing plaque characteristics within an arterial vessel; transmitting and receiving ultrasonic signals within the tissue region of interest based on the tissue motion frequency; calculating and obtaining coarse and fine strain image frames from the ultrasonic signals; fusing and matching the coarse and fine strain image frames to obtain precise strain image frames; counting the number of times each pixel in the precise strain image frame has a strain value, performing image processing on the statistical result to obtain a feature image; performing feature analysis on the feature image to obtain a feature analysis result, and outputting and displaying the feature analysis result. The present invention improves the detection accuracy of arterial plaque while minimizing the amount of data. It also calculates various plaque characteristic indicators based on the plaque strain information, conveniently and quickly transmitting the values ​​to a backend for processing and displaying them on a display for user diagnosis.

Description

Technical Field

[0001] The present invention relates to the technical field of ultrasonic medical equipment, and in particular to a method, device and storable medium for detecting arterial vascular plaque characteristics based on elastic imaging. Background Art

[0002] Carotid artery plaques, a manifestation of carotid atherosclerosis, are common at the bifurcation of the common carotid artery and are closely associated with ischemic stroke in the elderly. Medical color Doppler ultrasound (CDU) is the preferred noninvasive method for carotid artery examination and is widely used for screening and follow-up of carotid atherosclerosis. It not only demonstrates the location and size of plaques, as well as the location and severity of luminal stenosis, but also allows for hemodynamic assessment and morphological evaluation of plaques. Based on their morphology and echogenicity, plaques can be categorized as hypoechoic, lipid-rich soft plaques; moderately echogenic, collagen-rich, flat fibrous plaques; hyperechoic, calcified, hard plaques with acoustic shadowing; and mixed ulcerative plaques with varying echogenicity. Existing ultrasound methods for plaque assessment using CDU include IMT (intima-media thickness) measurement, which can quantitatively determine the presence of carotid atherosclerosis. However, assessment of carotid plaque characteristics still relies on subjective predictions based on hemodynamics, grayscale images, or morphology. In summary, in medical color Doppler ultrasound systems, the diagnosis of intra-arterial plaques is generally based on the echo strength, shape and size, assisted by hemodynamic characteristics for a subjective comprehensive diagnosis.

[0003] Therefore, a solution is urgently needed. Summary of the Invention

[0004] One of the purposes of the present invention is to provide an arterial plaque characteristic detection method based on elastic imaging, which can be conveniently and quickly deployed in a medical color ultrasound equipment system.

[0005] The method for detecting arterial plaque characteristics based on elastic imaging provided by an embodiment of the present invention is applied to medical color Doppler ultrasound equipment, including:

[0006] Acquiring a tissue region of interest containing plaque features within an arterial vessel;

[0007] Transmit and receive ultrasonic signals in the tissue area of ​​interest according to the tissue movement frequency;

[0008] Obtaining coarse and fine strain force image frames by calculating ultrasonic signals;

[0009] fusing and matching the coarse and fine strain image frames to obtain precise strain image frames;

[0010] Counting the number of times each pixel in the precise strain image frame has a strain value, and performing image processing on the statistical results to obtain a feature image;

[0011] Perform feature analysis on the feature image, obtain feature analysis results, and output them for display.

[0012] Preferably, when obtaining a tissue region of interest containing plaque features within an arterial vessel, the ultrasound software is controlled to automatically locate the arterial plaque region of the current image as the tissue region of interest, or the user can manually select the region of interest.

[0013] Preferably, the scanning sampling frequency of transmitting and receiving ultrasonic signals is more than twice the frequency of tissue movement, and the scanning sampling frequency of transmitting and receiving is adjustable by the user.

[0014] Preferably, the ultrasonic signal is received and two or more different methods are used to calculate and obtain coarse and fine strain image frames, wherein the strain image method calculated by the high-precision slow method is a fine image, and the strain image method calculated by the low-precision fast method is a coarse image.

[0015] Preferably, the coarse and fine strain image frames are fused and matched, and the fusion and matching includes matching of positions of pixels between the coarse and fine image frames and superposition of strain forces.

[0016] Preferably, the temporal variation pattern of any pixel in the precise strain image frame can be obtained. The precise strain image frame is a temporally continuous frame. The system can automatically or manually select quantitative measurement points on the continuous frame. The ultrasound software can display and draw the numerical value, curve, or chart of the quantitative measurement point changing over time.

[0017] Preferably, the number of times a certain pixel has a certain strain value is counted for the precise strain image containing the plaque, and the statistical process can be performed in a frozen state of the image or in a real-time state.

[0018] Preferably, the strain corresponding to a certain pixel includes a pixel point collective strain consisting of one or two or more pixel points.

[0019] Preferably, data processing is performed on the statistical results to form a characteristic image, wherein the data processing is a process of normalizing the statistical times or recalculating the probability distribution.

[0020] Preferably, feature analysis of the feature image can obtain various types of mathematical or image feature indicators. The feature analysis is based on the feature image after normalization or re-probability distribution. The mathematical or image feature indicators include but are not limited to brightness, energy, contrast, uniformity, non-uniformity, homogeneity, non-homogeneity, correlation, non-correlation and other Chinese or English abbreviations.

[0021] Preferably, the final strain image and characteristic index are displayed comprehensively and output to a display for user diagnosis, wherein the characteristic index can be displayed dynamically or statically at any position on the display in the form of numerical values, charts, curves but not limited to such forms.

[0022] The arterial plaque characteristic detection device based on elastic imaging provided by the present invention is applied to medical color Doppler ultrasound equipment, comprising:

[0023] Select the tissue region of interest module to automatically obtain the manually obtained carotid artery plaque area;

[0024] Transmitter and receiver modules are used to obtain ultrasonic signals of tissue deformation;

[0025] A coarse and fine strain image frame calculation module is used to calculate multiple frames of coarse and fine strain images;

[0026] A precise strain image frame calculation module, used to calculate precise strain image frames;

[0027] A feature image calculation module is used to calculate the number of times a certain pixel has a certain strain value in response to the strain;

[0028] Parameter acquisition module, used to extract mathematical or image feature indicators from feature images;

[0029] The display module is used to display the final strain image and characteristic index values.

[0030] The present invention provides a computer-storable medium, which is used to execute the above-mentioned method for detecting arterial vascular plaque characteristics based on elastic imaging.

[0031] This application has achieved the following beneficial effects:

[0032] This patent enables quantitative measurement and qualitative analysis of plaque within arterial vessels. This patent utilizes elastic imaging technology from medical ultrasound equipment to improve the accuracy of plaque detection within arterial vessels while minimizing data volume. The patent also calculates various plaque characteristic indicators based on plaque strain information, conveniently and quickly transmitting this information to the backend for processing and display on a monitor for user diagnosis.

[0033] Other features and advantages of the present invention will be described in the following description, and in part will become apparent from the description, or will be understood by practicing the present invention. The purposes and other advantages of the present invention can be realized and obtained by the structures particularly pointed out in the written description, claims, and drawings.

[0034] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] The accompanying drawings are used to provide a further understanding of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the present invention and do not constitute a limitation of the present invention. In the accompanying drawings:

[0036] Figure 1 Schematic diagram of a method for detecting arterial plaque characteristics based on elastography in an embodiment of the present invention;

[0037] Figure 2 Schematic diagram of an arterial plaque characteristic detection device based on elastography in an embodiment of the present invention;

[0038] Figure 3 Schematic diagram of calculating coarse and fine strain image frames in an embodiment of the present invention;

[0039] Figure 4 A schematic diagram of the definition of coarse and fine strain image frames in an embodiment of the present invention;

[0040] Figure 5 Schematic diagram of a plaque fusion method in coarse and fine strain images according to an embodiment of the present invention;

[0041] Figure 6 A schematic diagram of real-time display of parameter indicators in an embodiment of the present invention;

[0042] Figure 7 This is a schematic diagram of the final display in an embodiment of the present invention. DETAILED DESCRIPTION

[0043] The preferred embodiments of the present invention are described below with reference to the accompanying drawings. It should be understood that the preferred embodiments described herein are only used to illustrate and explain the present invention, and are not used to limit the present invention.

[0044] An embodiment of the present invention provides an arterial plaque characteristic detection method based on elastography, which is applied to medical color Doppler ultrasound equipment, including:

[0045] Acquiring a tissue region of interest containing plaque features within an arterial vessel;

[0046] Transmit and receive ultrasonic signals in the tissue area of ​​interest according to the tissue movement frequency;

[0047] Obtaining coarse and fine strain force image frames by calculating ultrasonic signals;

[0048] fusing and matching the coarse and fine strain image frames to obtain precise strain image frames;

[0049] Counting the number of times each pixel in the precise strain image frame has a strain value, and performing image processing on the statistical results to obtain a feature image;

[0050] Perform feature analysis on the feature image, obtain feature analysis results, and output them for display.

[0051] When acquiring a tissue region of interest containing plaque features within an arterial vessel, the ultrasound software is controlled to automatically locate the arterial plaque region in the current image as the tissue region of interest. The user can also manually select the region of interest.

[0052] The scanning sampling frequency of transmitting and receiving ultrasonic signals is more than twice the frequency of tissue movement, and the scanning sampling frequency of transmitting and receiving is adjustable by the user.

[0053] The ultrasonic signal can be received by using two or more different methods to calculate and obtain coarse and fine strain image frames, wherein the strain image method calculated by the high-precision slow method is the fine image, and the strain image method calculated by the low-precision fast method is the coarse image.

[0054] The coarse and fine strain image frames are fused and matched, and the fusion matching includes the matching of each pixel position between the coarse and fine image frames and the superposition of strain.

[0055] The temporal variation patterns of any pixel in the precise strain image frame can be obtained. The precise strain image frame is a continuous frame in time. The system can automatically or manually select quantitative measurement points on the continuous frame. The ultrasound software can display and draw the numerical values, curves or charts of the quantitative measurement points changing over time.

[0056] The number of times a certain pixel has a certain strain value for the precise strain image containing the plaque is counted, and the statistical process can be calculated in a frozen state or in a real-time state.

[0057] The strain corresponding to a certain pixel includes the strain of a pixel set consisting of one or two or more pixels.

[0058] The statistical results are processed to form a characteristic image, wherein the data processing is a process of normalizing the statistical times or recalculating the probability distribution.

[0059] Feature analysis of feature images can obtain various types of mathematical or image feature indicators. Feature analysis is based on the feature images after normalization or re-probability distribution. Mathematical or image feature indicators include but are not limited to brightness, energy, contrast, uniformity, non-uniformity, homogeneity, non-homogeneity, correlation, non-correlation and other Chinese or English abbreviations.

[0060] The final strain image and characteristic index are comprehensively displayed and output to a display for user diagnosis, wherein the characteristic index can be displayed dynamically or statically at any position on the display in the form of numerical values, charts, curves but not limited to such forms.

[0061] First, let’s explain the proper nouns involved, as follows:

[0062] Beamforming: Combining multiple ultrasonic signals into one signal is called beamforming;

[0063] Delay parameter: When multiple ultrasonic waves arrive at the target point, there is a sequence of arrival. This time difference is the delay parameter.

[0064] Dynamic reception: In medical color Doppler ultrasound equipment, the direction of ultrasound transmission is artificially defined as the axial direction. Ultrasound is recovered at short intervals along the axial direction. This reception method is called dynamic reception.

[0065] Focusing: Multiple ultrasonic signals arrive at a certain axial depth at the same time, which is called focusing;

[0066] Beam: The ultrasonic signal after beam synthesis is called a beam, which can also be called a beam line or line;

[0067] RF data: the final data after beamforming.

[0068] Specifically, the implementation case includes the following steps:

[0069] like Figure 1 In step S11, the input data of the tissue region of interest of the intra-arterial plaque feature is defined as Pre-Frame (t-1, i, j) or Post-Frame (t, i, j) to represent the previous frame or the next frame, such as Figure 3 In S30 or S31, t represents the frame number at time t; i represents the horizontal pixel label; j represents the vertical pixel label. The left and right boundaries of the region of interest are defined as [-L:1:L], i.e., with a step size of 1, the left boundary is -L, and the right boundary is L; the upper and lower boundaries are defined as [-H:1:H], i.e., with a step size of 1, the upper boundary is defined as -H, and the lower boundary is defined as H.

[0070] like Figure 1 In step S12, the front-end transmitting and receiving modules perform direct scanning in a sequential scanning manner. Specifically, after scanning the grayscale image, they scan the tissue within the regions of interest [-L:1:L] and [-H:1:H], and this cycle repeats. The repetition frequency of the transmitted and received ultrasonic signals is adjustable and is set to the FRF, which is generally at least twice the tissue deformation frame motion frequency.

[0071] like Figure 1 In step S13, the strain image of the tissue can be calculated by the displacement difference between the two frames of images, such as Figure 3In the image, the plaque S32 moves to S33. The position of the previous frame S30 is S35, which is represented by S0. S0 can be regarded as the center of mass of the plaque. In the next frame S31, the position is S2. When the plaque is deformed, the center of mass of the plaque moves from S0 to S2. The process can be considered that the center of mass S0 of the plaque first moves to S1, that is, the overall position of the plaque moves from S32 to S37, and then the center of mass moves from S1 to S2, that is, the overall position of the plaque moves from S37 to S36. Therefore, the displacement vector S34 of the real plaque deformation is expressed as follows:

[0072]

[0073] The vector addition diagram is as follows Figure 5 shown

[0074] The coarse strain image frame in the coarse and fine strain image frames can correspond to the longitudinal strain image frame, that is, the displacement from S0 to S1 The fine strain image frame can correspond to the lateral strain image frame, that is, the displacement from S1 to S2

[0075] like Figure 4 Where S40 is the decomposed longitudinal strain image frame, defined as PH(t, i, j), S41 is the decomposed transverse strain image frame, defined as PL(t, i, j), S42 is the longitudinal direction, and S43 is the transverse direction;

[0076] like Figure 1 In step S13, the coarse and fine strain image frames of the tissue are calculated, wherein the coarse strain image is obtained by using a phase iterative difference method. A calculation formula is as follows:

[0077]

[0078] in:

[0079] w0 is the angular frequency of the deformation; arg: represents the phase angle; H0∈[-H:1:H]; p(t,i,j) represents the displacement of the pixel in the i-th row and j-th column in the frame at time t;

[0080] like Figure 1 In step S13, the coarse and fine strain image frames of the tissue are calculated, wherein the fine strain image is calculated using the following formula:

[0081] pl(t,i,j)=max j {R(t,i,j)}

[0082]

[0083] Where: L0∈[-L:1:L], I: is a fixed row on the image, max j {} is the maximum value function in the j direction, l represents the variable in the interval [-L0, L0]; dl represents the integral variable;

[0084] like Figure 1 In step S14, the coarse strain image frame PH(t, i, j) and the fine strain image frame PL(t, i, j) are fused and matched to obtain the precise strain image frame. A fusion method can directly use the displacement information of pl(t, i, j) to match the p in step 5. t-1 (t, i, j) displacement frame vector calibration, as in step 3 The formula shown is calculated Later In the vector direction, the formula for calculating PH(t, i, j) in step 5 is used to obtain the precise strain force image frame P(t, i, j);

[0085] For the precise strain image frame in step 7, one or more sampling points can be selected in the region of interest, that is, i and j can be arbitrarily selected to obtain the precise strain value P(t) that changes with time, such as Figure 6 As shown in S64, the horizontal axis is the time axis S63, the vertical axis is parameter 0, that is, S60, parameter 0 can be defined as the precise strain value P(t), and the displayed curve is as S64. S62 is the curve 0 in the icon drawn with a solid line or the curve 1 drawn with a dotted line;

[0086] like Figure 1 In step S15, the number of times a certain strain pair PV(t, c(v), i, j) has a certain strain value co(v, n) is counted for the precise strain image P(t, i, j), and the statistical results are processed to form a feature image FI(t, n, n). The specific formula is as follows:

[0087] FI(t,n,n)=SF{PV(t,c(v),i,j)==co(v,n)} / Factor

[0088] Where: v represents the pixel sequence in the pixel pair, v = 1, 2, 3...; N is the maximum strain force, co(v, n) is the prior input value; SF{} is the statistical frequency function of strain force co(v, n); n∈[0, N-1]; Factor is the normalization factor.

[0089] like Figure 1 As shown, S16 performs feature analysis on the feature image FI(t, n, n) to obtain various types of mathematical or image feature indicators, such as Figure 6Parameters 1, 2, 3, 4, etc. shown in S61 in this patent are generally applicable. Based on FI(t, n, n), various parameters such as the energy Powet of the plaque tissue can be calculated. One of the energy Power formulas is as follows:

[0090] Power=sum{FI(t,n,n)*FI(t,n,n)},n∈[0,N-1]

[0091] Where sum{} is the sum function;

[0092] like Figure 1 As shown in step S17, the final strain image and characteristic index are comprehensively displayed and output to the display for user diagnosis, such as Figure 7 As shown, S70 and S71 are dual real-time display images, S70 is the grayscale image area, S71 is the composite image of the grayscale image S73 and the precise strain image S72, where the image of S73 is consistent with the grayscale image of S70, and S74 is the characteristic highlighting area for displaying various parameters such as the indicator and chart. Figure 6 shown.

[0093] This patent enables quantitative measurement and qualitative analysis of plaque within arterial vessels. This patent utilizes elastic imaging technology from medical ultrasound equipment to improve the accuracy of plaque detection within arterial vessels while minimizing data volume. The patent also calculates various plaque characteristic indicators based on plaque strain information, conveniently and quickly transmitting this information to the backend for processing and display on a monitor for user diagnosis.

[0094] The arterial plaque characteristic detection device based on elastic imaging provided by the present invention is applied to medical color Doppler ultrasound equipment, comprising:

[0095] Select the tissue region of interest module to automatically obtain the manually obtained carotid artery plaque area;

[0096] Transmitter and receiver modules are used to obtain ultrasonic signals of tissue deformation;

[0097] A coarse and fine strain image frame calculation module is used to calculate multiple frames of coarse and fine strain images;

[0098] A precise strain image frame calculation module, used to calculate precise strain image frames;

[0099] A feature image calculation module is used to calculate the number of times a certain pixel has a certain strain value in response to the strain;

[0100] Parameter acquisition module, used to extract mathematical or image feature indicators from feature images;

[0101] The display module is used to display the final strain image and characteristic index values.

[0102] like Figure 2 As shown, in this embodiment, a device for detecting arterial vascular plaque characteristics based on elastic imaging of a medical color ultrasound device includes:

[0103] A tissue region of interest selection module S20 is used to select an effective plaque diagnosis area;

[0104] Transmitting and receiving module S21, used to obtain ultrasonic echo data including tissue deformation;

[0105] A coarse and fine strain image frame calculation module S22 is used to calculate longitudinal and transverse strain image frame data of the user's area of ​​interest;

[0106] The precise strain image frame calculation module S23 is used to fuse the coarse and fine strain images to obtain a precise strain image;

[0107] A characteristic image calculation module S24 is used to obtain a data frame for calculating plaque characteristic parameters;

[0108] Parameter acquisition module S25, used to calculate the index value or chart of plaque characteristics;

[0109] Display module S26, used to display strain images and various parameter indicators of plaque characteristics;

[0110] The present invention provides a computer-storable medium, which is used to execute the above-mentioned arterial vascular plaque characteristic detection method based on elastic imaging.

[0111] Obviously, those skilled in the art may make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if such changes and modifications fall within the scope of the claims and their equivalents, the present invention is intended to include such changes and modifications.

Claims

1. A method for detecting arterial plaque characteristics based on elastic imaging, applied to medical color Doppler ultrasound equipment, characterized in that: include: Acquire a tissue region of interest containing arterial plaque features; Transmit and receive ultrasonic signals in the tissue area of ​​interest according to the tissue movement frequency; the scanning sampling frequency of the transmitted and received ultrasonic signals is more than twice the tissue movement frequency; Obtaining a coarse strain image frame and a fine strain image frame by calculating the ultrasonic signal; wherein the strain image frame calculated by a high-precision slow method is the fine strain image frame, and the strain image frame calculated by a low-precision fast method is the coarse strain image frame; The coarse strain image frame and the fine strain image frame are fused and matched to obtain a precise strain image frame; the fusion matching includes matching the positions of each pixel in the coarse strain image frame and the fine strain image frame and superposition of strain; Count the number of times each pixel in the precise strain image frame has a strain value; process the counted number of times and form a feature image, wherein the data processing is to perform normalization calculation on the counted number of times; perform statistics in a frozen image state or in real time; a pixel-to-strain force includes a pixel-point aggregate strain force composed of multiple pixels; Perform feature analysis on the feature image to obtain feature analysis results, which are then output and displayed; the feature analysis results include image feature indicators, which include brightness, contrast, uniformity, and correlation.

2. The method for detecting arterial plaque characteristics based on elastic imaging according to claim 1, wherein: When acquiring a tissue region of interest containing arterial vascular plaque features, the ultrasound software is controlled to automatically locate the arterial vascular plaque region of the current image as the tissue region of interest.

3. The method for detecting arterial plaque characteristics based on elastic imaging according to claim 1, wherein: The scanning sampling frequency of transmitting and receiving ultrasonic signals can be adjusted.

4. The method for detecting arterial plaque characteristics based on elastic imaging according to claim 1, wherein: The temporal variation pattern of any pixel in the precise strain image frame can be obtained. The precise strain image frame is a continuous frame in time. Quantitative measurement points can be automatically selected on the continuous frame, and the ultrasound software can draw a curve showing the variation of the quantitative measurement points over time.

5. The method for detecting arterial plaque characteristics based on elastic imaging according to claim 1, wherein: The image characteristic index is displayed on the monitor as a curve.

6. An arterial plaque characteristic detection device based on elastic imaging, applied to medical color Doppler ultrasound equipment, characterized in that: include: Select the tissue region of interest module to obtain the tissue region of interest containing arterial plaque features; The transmitting and receiving module is used to transmit and receive ultrasonic signals in the tissue area of ​​interest according to the tissue movement frequency; the scanning sampling frequency of the transmitting and receiving ultrasonic signals is more than twice the tissue movement frequency; a coarse and fine strain image frame calculation module, configured to calculate a coarse strain image frame and a fine strain image frame from an ultrasonic signal; wherein the strain image frame calculated using a high-precision slow method is a fine strain image frame, and the strain image frame calculated using a low-precision fast method is a coarse strain image frame; A precise strain image frame calculation module is used to fuse and match the coarse strain image frame and the fine strain image frame to obtain a precise strain image frame; the fusion matching includes matching the positions of each pixel in the coarse strain image frame and the fine strain image frame and superposition of strain forces; The feature image calculation module is used to count the number of times each pixel in the precise strain image frame has a strain value; the statistical number is processed to form a feature image, wherein the data processing is to perform normalization calculation on the statistical number; the statistics are performed in the image frozen state or in real time; a pixel corresponding to the strain includes the pixel point collective strain composed of multiple pixels; The parameter acquisition module is used to perform feature analysis on the feature image and obtain feature analysis results; the feature analysis results include image feature indicators, and the image feature indicators include brightness, contrast, uniformity, and correlation; The display module is used to display image feature indicators.

Images