Ultrasound thermometry method, apparatus, and storage medium

Abstract

This application discloses an ultrasonic temperature measurement method, device, and storage medium, relating to the field of medical imaging technology. The method includes: during ultrasonic temperature measurement, if the ultrasonic probe's scanning at the target scanning position is interrupted, restoring the ultrasonic probe to the target scanning position; detecting whether the interruption time of the ultrasonic probe's scanning at the target scanning position reaches a time threshold; if the interruption time reaches the time threshold, calculating the temperature of the scanned object based on the envelope signal of the ultrasonic echo signal. This solves the problem in the prior art where temperature measurement cannot continue once interrupted, achieving temperature calculation using an envelope signal with slightly lower resolution but a larger coverage area, realizing the continuation of temperature monitoring after the interruption process and basically maintaining the accuracy of the calculated temperature.

Description

Technical Field

[0001] This invention relates to an ultrasonic temperature measurement method, device, and storage medium, belonging to the field of medical imaging technology. Background Technology

[0002] In ablation or hyperthermia procedures monitored by ultrasound, existing methods use the cumulative offset of the ultrasound analytical signal to calculate a two-dimensional temperature image.

[0003] However, there is a risk of ultrasound signal interruption during the procedure due to the patient's heavy breathing, coughing, or movement. This is because the offset time between adjacent frames or several adjacent frames is calculated based on the phase angle method of ultrasound signal analysis. Only when it is less than the ultrasonic carrier period The timing is accurate, so once the process is interrupted, the ultrasonic offset time is... Exceeded (For example, in one instance of high-power microwave ablation, the interruption exceeded...) Second, Reachable If the reading is incorrect, the temperature monitoring cannot be calculated correctly, causing the monitoring to stop. Summary of the Invention

[0004] The purpose of this invention is to provide an ultrasonic temperature measurement method, device, and storage medium to solve the problems existing in the prior art.

[0005] To achieve the above objectives, the present invention provides the following technical solution:

[0006] According to a first aspect, embodiments of the present invention provide an ultrasonic temperature measurement method, the method comprising:

[0007] If the ultrasonic probe is interrupted during ultrasonic temperature measurement at the target scanning position, the ultrasonic probe will be restored to the target scanning position.

[0008] The detection method checks whether the interruption time of the ultrasound probe scanning at the target scanning position reaches a time threshold.

[0009] If the interruption time reaches the time threshold, the temperature of the object being scanned is calculated based on the envelope signal of the ultrasonic echo signal.

[0010] Optionally, detecting whether the interruption time of the ultrasound probe's scanning at the target scanning location reaches a time threshold includes:

[0011] The interruption time is calculated based on the time of interruption and the time of recovery, and it is detected whether the interruption time reaches the time threshold.

[0012] And / or based on n frames of ultrasound images prior to the interruption of scanning at the target scanning location, obtain the predicted interruption time offset; detect whether the interruption time offset reaches the time threshold, where n is a positive integer.

[0013] Optionally, if the detected interruption time offset reaches the time threshold, the step of calculating the temperature of the scanned object based on the envelope signal of the ultrasonic echo signal includes:

[0014] The predicted signal offset pixel number is obtained from n frames of ultrasound images before the ultrasound probe interrupts scanning at the target scanning position, where n is a positive integer.

[0015] The window size of the sliding window is determined based on the number of pixels offset from the predicted signal;

[0016] Based on the i-th frame of ultrasound image before interruption and the j-th frame of ultrasound image after recovery, calculate the signal parameters and the initial signal offset pixel number in each sliding window within the region of interest. The signal parameters include the correlation function value in the depth direction and / or the envelope signal value in the depth direction, where i and j are positive integers.

[0017] Based on the signal parameters, identify the strong echo region within the region of interest;

[0018] The temperature within the region of interest is calculated based on the strong echo region and the initial signal offset pixel count.

[0019] Optionally, if the signal parameters include correlation function values ​​in the depth direction, identifying strong echo regions within the region of interest based on the signal parameters includes:

[0020] Calculate the first threshold based on the relevant function value;

[0021] The region whose correlation function value is greater than the first threshold is identified as the strong echo region.

[0022] Optionally, if the signal parameters include envelope signal values ​​in the depth direction, identifying strong echo regions within the region of interest based on the signal parameters includes:

[0023] Calculate the second threshold based on the envelope signal value;

[0024] Large signals in the sliding window are identified based on the envelope signal value and the second threshold.

[0025] If the width ratio of the identified large signal in the sliding window exceeds the width threshold, then the sliding window is identified as the strong echo region.

[0026] Optionally, if the signal parameters include correlation function values ​​in the depth direction and envelope signal values ​​in the depth direction, the step of identifying strong echo regions within the region of interest based on the signal parameters includes:

[0027] The first strong echo region is identified based on the correlation function value in the depth direction;

[0028] The second strong echo region is identified based on the envelope signal value in the depth direction;

[0029] The intersection of the first strong echo region and the second strong echo region is determined as the strong echo region within the region of interest.

[0030] Optionally, calculating the temperature within the region of interest based on the strong echo region and the initial signal offset pixel count includes:

[0031] The initial signal offset pixel count of the strong echo region is interpolated and fitted to obtain the actual signal offset pixel count of the region of interest.

[0032] Get the cumulative offset pixel count of the signal before the interruption;

[0033] The temperature within the region of interest is calculated based on the cumulative offset pixel count of the signal before the interruption and the actual offset pixel count of the signal.

[0034] Optionally, the method further includes:

[0035] If the interruption time does not reach the time threshold, the temperature of the scanned object is calculated based on the analytical signal of the ultrasonic echo signal.

[0036] In a second aspect, an ultrasonic temperature measuring device is provided, the device including a memory and a processor, the memory storing at least one program instruction, the processor loading and executing the at least one program instruction to implement the method as described in the first aspect.

[0037] Thirdly, a computer storage medium is provided, wherein at least one program instruction is stored therein, the at least one program instruction being loaded and executed by a processor to implement the method as described in the first aspect.

[0038] By resuming the ultrasonic probe to the target scanning position if the ultrasonic scanning is interrupted during ultrasonic temperature measurement, and detecting whether the interruption time reaches a time threshold, the temperature of the scanned object is calculated based on the envelope signal of the ultrasonic echo signal. This solves the problem in existing technologies where temperature measurement cannot continue once interrupted. It achieves temperature calculation by using an envelope signal with lower resolution but a larger coverage area to calculate the number of offset pixels, thus continuing temperature monitoring after interruption and maintaining the accuracy of the calculated temperature.

[0039] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it in accordance with the contents of the specification, the preferred embodiments of the present invention are described in detail below with reference to the accompanying drawings. Attached Figure Description

[0040] Figure 1 This is a flowchart of an ultrasonic temperature measurement method provided in one embodiment of the present invention;

[0041] Figures 2a to 2d This is a possible schematic diagram of identifying a strong echo region based on signal parameters, according to one embodiment of the present invention. Detailed Implementation

[0042] The technical solution of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention. In the description of the present invention, it should be noted that the terms "center," "upper," "lower," "left," "right," "vertical," "horizontal," "inner," and "outer," etc., indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation of the present invention. In addition, the terms "first," "second," and "third" are used for descriptive purposes only and should not be construed as indicating or implying relative importance. In the description of the present invention, it should be noted that unless otherwise expressly specified and limited, the terms "installed," "connected," and "linked" should be interpreted broadly. For example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection or an indirect connection through an intermediate medium; it can be a connection within two components. Those skilled in the art will understand the specific meaning of the above terms in this invention according to the specific circumstances. Furthermore, the technical features involved in the different embodiments of this invention described below can be combined with each other as long as they do not conflict with each other.

[0043] Please refer to Figure 1 It illustrates a flowchart of an ultrasonic temperature measurement method provided in one embodiment of this application, as follows: Figure 1 As shown, the method includes:

[0044] Step 101: During the ultrasonic temperature measurement process, if the ultrasonic probe is interrupted at the target scanning position, the ultrasonic probe is restored to the target scanning position.

[0045] The solution proposed in this application is for scenarios involving temperature monitoring of a scanned object, such as in ultrasound ablation procedures. Taking temperature monitoring during ultrasound ablation surgery as an example, at the start of the procedure, the surgeon can fix the ultrasound probe, pre-store a grayscale (B) image of the object before ablation, and mark the scanning position of the ultrasound probe on the surface of the object. At this point, T-mode (temperature mode) monitoring can be performed. After entering T-mode, the ultrasound device can calculate the offset pixel count between adjacent frames based on the analytical signal of the real-time acquired ultrasound echo signal, i.e., from the start of the ablation surgery to the... At frame rate, the calculated number of offset pixels is: ,in, This represents the offset in pixels between two adjacent frames. When the last ultrasound image before the interruption is obtained, the cumulative offset in pixels before the interruption can be calculated.

[0046] During real-time scanning of the target location by the ultrasound probe, a real-time B-image is stored at regular intervals, such as 10 seconds. The ultrasound scan at the target location is continuously monitored for interruptions. If an interruption occurs, the ultrasound probe is restored to the target location. The ultrasound equipment can detect changes in the scan position by observing changes in the ultrasound images, thus determining whether the ultrasound scan has been interrupted. Specifically, if there is a significant change in the B-image between two adjacent ultrasound frames, the scan position is considered to have changed. In practice, the matching degree between two adjacent ultrasound frames can be calculated. If the matching degree is lower than a matching degree threshold, the scan position is considered to have changed and the scan has been interrupted; otherwise, the scan position has not changed, meaning the scan has not been interrupted. The matching degree threshold can be a preset value or a value calculated based on the B-images stored before and during ablation; no limitation is made here.

[0047] Upon detecting a scan interruption, the ultrasound probe can be restored to the target scan position. Optionally, when the ultrasound device receives a trigger signal from the reset button, it enters reset mode, after which the ultrasound probe is restored to the target scan position. In practice, the ultrasound probe can be moved manually by the doctor or by a robotic arm. The matching degree between the real-time acquired ultrasound image and the pre-stored image before the interruption, especially the last frame of the ultrasound image, is used to detect whether the ultrasound probe has moved to the target scan position. After the ultrasound probe is fixed in the target scan position, upon receiving a trigger signal from the preset button, it enters reconnection mode, at which point real-time scanning resumes.

[0048] It should be noted that the ultrasonic ablation in this application may include radiofrequency needle ablation, microwave needle ablation or high-focus ultrasonic scalpel ablation, and this application does not limit the specific ablation method.

[0049] Step 102: Detect whether the interruption time of the ultrasound probe scanning at the target scanning position reaches the time threshold;

[0050] When an interruption in the ultrasound scan is detected, it can be checked whether the interruption time has reached a time threshold. Optionally, this step includes at least one of the following two possible implementations:

[0051] In one possible implementation, the interruption time is calculated based on the time of interruption and the time of recovery, and it is then checked whether the interruption time has reached a time threshold. Optionally, when an interruption is detected, a timer can be started and stopped when the scan is resumed to the target scanning position, thus obtaining the interruption time. In one possible embodiment, the moment when the ultrasound image changes excessively is... The moment when the ultrasound probe is repositioned to match the previous ultrasound image is... The time when the reconnection command is issued via button is In this application and That is, the start and end points of the interruption period, that is, the interruption time is - .

[0052] The time threshold can be a preset value or a custom value, and its specific implementation is not limited.

[0053] In the second possible implementation, the predicted interruption time offset is obtained based on n frames of ultrasound images before the ultrasound probe interrupts scanning at the target scanning position; it is then detected whether the interruption time offset reaches the time threshold, where n is a positive integer.

[0054] Optionally, the predicted temperature rise during the interruption period is obtained by polynomial fitting based on n frames of ultrasound images and the interruption time, for each strip within the region of interest. The data line is obtained by integrating both sides of formula (1). This yields the predicted interruption time offset. Formula (1) is as follows: In the formula, It is depth Temperature changes at that location It is the initial sound speed within the organization. It is the tissue expansion coefficient. It is the coefficient of sound speed as a function of temperature. It is deep inside Cumulative time offset, It is the digital sampling interval of the ultrasound image in the z-direction. It is determined by the sampling frequency and the speed of sound.

[0055] In the second possible implementation, the time threshold can be the carrier period, that is, , This is the ultrasonic carrier frequency.

[0056] The above example illustrates detection using only one method. In actual implementation, detection can also be performed using the first possible method. If the detection result is greater than the time threshold, further detection is performed using the second possible method. If both detection results are greater than the time threshold, the interruption time is determined to be greater than the time threshold. When two detection methods are used, the setting of the time threshold for each method is as described in the previous sections and will not be repeated here. Furthermore, if the detection result of the first method is not greater than the time threshold, the temperature can be directly calculated from the analyzed signal, and the process ends. Since the computational complexity of predicting the interruption time offset is higher than that of directly calculating the interruption time, the above scheme achieves the effect of reducing computational complexity.

[0057] It should be noted that the region of interest can be an area manually marked by the doctor or an area identified by the ultrasound equipment based on the ultrasound image. The region of interest can be a lesion area in the scanned object, such as a tumor area in the liver, etc., and its specific content is not limited.

[0058] Step 103: If the interruption time reaches the time threshold, calculate the temperature of the object being scanned based on the envelope signal of the ultrasonic echo signal.

[0059] Optionally, this step may include:

[0060] First, the number of predicted signal offset pixels is obtained from n frames of ultrasound images before the ultrasound probe interrupts scanning at the target scanning position, where n is a positive integer.

[0061] Optionally, the predicted temperature rise during the interruption period is obtained by polynomial fitting based on n frames of ultrasound images and the interruption duration, for each strip within the region of interest. The data line is obtained by integrating both sides of formula (1). The corresponding signal offset pixel number can be calculated using formula (2). Among them, subscript It indicates estimated value. Indicates the direction of the probe surface or Directional index (pixel index) Indicates the direction of ultrasound propagation or the depth direction of the image or The direction index (pixel index). Formula (2) is: In the formula, It is depth Temperature changes at that location It is the initial sound speed within the organization. It is the tissue expansion coefficient. It is the coefficient of sound speed as a function of temperature. It is deep inside Cumulative time offset, The digital sampling interval of the ultrasound image in the z-direction is determined by the sampling frequency and the velocity of sound.

[0062] Second, the window size of the sliding window is determined based on the number of pixels offset by the predicted signal;

[0063] The window size of the sliding window is: .in, This is a proportional parameter, and its value can be an empirical value or a user-defined value. Typically, it is... The range of values ​​is .

[0064] Third, based on the i-th frame ultrasound image before the interruption and the j-th frame ultrasound image after recovery, calculate the signal parameters and the actual signal offset pixel number in each sliding window within the region of interest. The signal parameters include the correlation function value in the depth direction and / or the envelope signal value in the depth direction, where i and j are positive integers.

[0065] Here, the i-th frame of ultrasound image before the interruption refers to the last i-th frame before the interruption, such as the last frame before the interruption, or the second to last frame, etc. Similarly, the j-th frame of ultrasound image after recovery refers to the j-th frame of ultrasound image counted from the recovery time, such as the first frame after recovery, or the second frame after recovery, etc. In actual implementation, we will use the calculation based on the last frame before the interruption and the first frame after recovery as an example, where i and j are both 1. And, unless otherwise specified, the following examples will all use i=j=1.

[0066] After calculating the window size of the sliding window, the signal envelope at each x-position is traversed using the sliding window. The cross-correlation between the last ultrasound image before the interruption and the first ultrasound image after the recovery is then performed to obtain all the signals. The actual offset of the signal segment within the center sliding window (in pixels) and related function values That is, for each data line at position x, the offset number of pixels of the signal at the center point of the sliding window is calculated within the sliding window, that is, within a space.

[0067] Fourth, based on the signal parameters, identify strong echo regions within the region of interest;

[0068] Optionally, the specific identification method in this step may vary depending on the signal parameters. Therefore, the following will explain the different situations separately.

[0069] When the signal parameters include correlation function values ​​in the depth direction, this step includes:

[0070] (1) Calculate the first threshold based on the relevant function values;

[0071] Optionally, the first threshold can be ,in, Standard deviation It is the sample size. This is the proportionality coefficient.

[0072] (2) Set the relevant function values Greater than the first threshold The region was identified as a strong echo zone.

[0073] When the correlation function value is greater than the first threshold, the current signal region is considered a valid signal region, i.e., identified as a strong echo region. Please refer to [reference needed]. Figure 2a and Figure 2b It illustrates a schematic diagram for determining strong echo regions when the signal parameters are the correlation function values ​​in the depth direction. Among them, Figure 2a The x-axis represents the pixel value in the depth direction, and the y-axis represents the calculated value of the relevant function. Figure 2b The horizontal axis represents the pixels in the depth direction, and the vertical axis represents the calculated signal offset in pixels.

[0074] When the signal parameters include envelope signal values ​​in the depth direction, this step includes:

[0075] (1) Calculate the second threshold based on the envelope signal value;

[0076] The second threshold can be ,in, Standard deviation It is the sample size. This is the proportionality coefficient.

[0077] (2) Identify large signals in the sliding window based on the envelope signal value and the second threshold;

[0078] In this application, signals with envelope signal values ​​greater than a second threshold are identified as large signals.

[0079] (3) If the width ratio of the large signals in the i-th and j-th frames exceeds the width threshold, the center point of the sliding window is identified as a strong echo area.

[0080] The width threshold is a preset value. In actual implementation, the width threshold can be 1 / 2, 1 / 3, etc.

[0081] For example, please refer to Figure 2c and Figure 2d It illustrates a schematic diagram for identifying strong echo regions based on envelope signals. Among them, Figure 2cThe horizontal axis represents the pixel value in the depth direction, and the vertical axis represents the calculated width ratio. Figure 2d The horizontal axis represents the pixel value in the depth direction, and the vertical axis represents the signal offset pixel value. Figure 2c The solid line represents the width ratio of the large signal in the i-th frame of the ultrasound image within the window width, and the dashed line represents the width ratio of the large signal in the j-th frame of the ultrasound image within the window width. Figure 2d The bolded part indicates the location of the identified strong echo area.

[0082] When the signal parameters include both the correlation function value in the depth direction and the envelope value in the depth direction, this step includes:

[0083] (1) Identify the first strong echo region based on the correlation function value in the depth direction;

[0084] (2) Identify the second strong echo region based on the envelope signal value in the depth direction;

[0085] The identification methods for the first and second strong echo regions are similar to those described above, and will not be repeated here.

[0086] (3) The intersection of the first strong echo region and the second strong echo region is determined as the strong echo region within the region of interest.

[0087] For example, the two-dimensional subscript range of the strong echo region is: Since the calculation is performed separately for each data line, meaning the x-position remains unchanged, this step only requires calculating the intersection of the first and second strong echo regions in the z-direction.

[0088] Fifth, calculate the temperature within the region of interest based on the strong echo zone and the actual number of signal offset pixels.

[0089] In practice, this step may include:

[0090] (1) Interpolate and fit the actual signal offset pixel count in the strong echo region to obtain the actual signal offset pixel count in the region of interest;

[0091] Optional, for Two-dimensional interpolation is performed. Furthermore, in practice, interpolation in the z-direction can be performed first, followed by interpolation in the x-direction, thus obtaining the complete result. .

[0092] (2) Obtain the number of signal offset pixels before the interruption;

[0093] Optionally, the ultrasound device may acquire the cumulative offset pixel count of the signal calculated and stored in step 101.

[0094] (3) Calculate the temperature in the region of interest based on the number of signal offset pixels before the interruption and the actual number of signal offset pixels.

[0095] Optionally, the sum of the signal offset pixels before the interruption and the actual signal offset pixels in the strong echo region is calculated, and the temperature in the region of interest is calculated based on the calculated sum. Specifically, the temperature is calculated based on the calculated sum and formula (2) in step 102. Then, the latest temperature data can be calculated according to formula (1).

[0096] It should be noted that if the interruption time does not reach the time threshold, the temperature of the scanned object is calculated based on the analytical signal of the ultrasonic echo signal. This will not be elaborated further in this application.

[0097] In summary, by restoring the ultrasonic probe to the target scanning position if the ultrasonic scanning is interrupted during ultrasonic temperature measurement, and detecting whether the interruption time reaches a time threshold, the temperature of the scanned object is calculated based on the envelope signal of the ultrasonic echo signal. This solves the problem in existing technologies where temperature measurement cannot continue once interrupted, achieving temperature calculation using an envelope signal with lower resolution but wider coverage, thus continuing temperature monitoring after interruption and essentially maintaining the accuracy of the calculated temperature.

[0098] Furthermore, this application identifies strong echo regions and calculates the temperature based on the envelope signal within the strong echo regions, thereby improving the accuracy of the calculated temperature.

[0099] This application also provides an ultrasonic temperature measuring device, the device including a memory and a processor, the memory storing at least one program instruction, the processor loading and executing the at least one program instruction to implement the method described above.

[0100] This application also provides a computer storage medium storing at least one program instruction, which is loaded and executed by a processor to implement the method described above.

[0101] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0102] The embodiments described above are merely illustrative of several implementations of the present invention, and while the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the invention patent. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these all fall within the protection scope of the present invention. Therefore, the protection scope of this invention patent should be determined by the appended claims.

Claims

1. An ultrasonic thermometry method, characterized by, The method includes: If the ultrasonic probe is interrupted during ultrasonic temperature measurement at the target scanning position, the ultrasonic probe will be restored to the target scanning position. The detection method checks whether the interruption time of the ultrasound probe scanning at the target scanning position reaches a time threshold. If the interruption time reaches the time threshold, the temperature of the object being scanned is calculated based on the envelope signal of the ultrasonic echo signal. The calculation of the temperature of the object being scanned based on the envelope signal of the ultrasonic echo signal includes: Based on the i-th frame of ultrasound image before interruption and the j-th frame of ultrasound image after recovery, calculate the signal parameters and the initial signal offset pixel number in each sliding window within the region of interest. The signal parameters include the correlation function value in the depth direction and / or the envelope signal value in the depth direction, where i and j are positive integers. Based on the signal parameters, identify the strong echo region within the region of interest; The temperature within the region of interest is calculated based on the strong echo region and the initial signal offset pixel count.

2. The method of claim 1, wherein, The detection of whether the interruption time of the ultrasound probe scanning at the target scanning position reaches a time threshold includes: The interruption time is calculated based on the time of interruption and the time of recovery, and it is detected whether the interruption time reaches the time threshold. And / or, Based on n frames of ultrasound images before the ultrasound probe interrupts scanning at the target scanning position, the predicted interruption time offset is obtained; it is then detected whether the interruption time offset reaches the time threshold, where n is a positive integer.

3. The method of claim 1, wherein, The calculation of the temperature of the scanned object based on the envelope signal of the ultrasonic echo signal further includes: Based on the n frames of ultrasound images before the ultrasound probe interrupts scanning at the target scanning position, the number of pixels offset for the predicted signal is obtained, where n is a positive integer. The window size of the sliding window is determined based on the number of pixels offset from the predicted signal.

4. The method according to claim 3, characterized in that, If the signal parameters include correlation function values ​​in the depth direction, the step of identifying strong echo regions within the region of interest based on the signal parameters includes: Calculate the first threshold based on the relevant function value; The region whose correlation function value is greater than the first threshold is identified as the strong echo region.

5. The method according to claim 3, characterized in that, If the signal parameters include envelope signal values ​​in the depth direction, the step of identifying strong echo regions within the region of interest based on the signal parameters includes: Calculate the second threshold based on the envelope signal value; Large signals in the sliding window are identified based on the envelope signal value and the second threshold. If the width ratio of the identified large signal in the sliding window exceeds the width threshold, then the sliding window is identified as the strong echo region.

6. The method according to claim 3, characterized in that, If the signal parameters include correlation function values ​​in the depth direction and envelope signal values ​​in the depth direction, the step of identifying strong echo regions within the region of interest based on the signal parameters includes: The first strong echo region is identified based on the correlation function value in the depth direction; The second strong echo region is identified based on the envelope signal value in the depth direction; The intersection of the first strong echo region and the second strong echo region is determined as the strong echo region within the region of interest.

7. The method according to claim 3, characterized in that, The step of calculating the temperature within the region of interest based on the strong echo region and the initial signal offset pixel count includes: The initial signal offset pixel count of the strong echo region is interpolated and fitted to obtain the actual signal offset pixel count of the region of interest. Get the cumulative offset pixel count of the signal before the interruption; The temperature within the region of interest is calculated based on the cumulative offset pixel count of the signal before the interruption and the actual signal offset pixel count of the region of interest.

8. The method according to any one of claims 1 to 7, characterized in that, The method further includes: If the interruption time does not reach the time threshold, the temperature of the scanned object is calculated based on the analytical signal of the ultrasonic echo signal.

9. An ultrasonic temperature measuring device, characterized in that, The device includes a memory and a processor, the memory storing at least one program instruction, and the processor loading and executing the at least one program instruction to implement the method as described in any one of claims 1 to 8.

10. A computer storage medium, characterized in that, The computer storage medium stores at least one program instruction, which is loaded and executed by a processor to implement the method as described in any one of claims 1 to 8.

Images