Photoacoustic imaging method and system based on data analysis

Abstract

The invention relates to the technical field of photoacoustic imaging, in particular to a photoacoustic imaging method and system based on data analysis, and the method comprises the steps: determining the imaging state of an initial photoacoustic image according to an image abnormal coefficient and a detection domain variable value, and determining an analysis mode according to the imaging state; in the detection anomaly analysis, determining a first optimization mode according to the noise modulation coefficient; in tissue influence analysis, a second optimization mode is determined according to a distortion influence coefficient and a regional dislocation coefficient, and the second optimization mode is that an adjustment mode is determined according to a blood flow influence index, or a compensation mode is determined according to a tissue difference coefficient and deviation uniformity; the adjusting mode is that the detection angle range is adjusted according to the blood flow influence index, or the laser repetition frequency is adjusted according to the distortion influence coefficient, the compensation mode is characteristic frequency compensation or related frequency compensation, and the image quality of photoacoustic imaging can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of photoacoustic imaging, and in particular to a photoacoustic imaging method and system based on data analysis. Background Art

[0002] As an emerging biomedical imaging technology, photoacoustic imaging combines the high contrast of optical imaging and the deep penetration of ultrasonic imaging, and can provide structural and functional information inside tissues. However, there are problems such as artifacts, blurring, and insufficient resolution in photoacoustic imaging, resulting in unstable quality of the generated photoacoustic images. Therefore, how to improve the image quality of photoacoustic imaging is a technical problem that needs to be solved urgently by those skilled in the art.

[0003] Chinese Patent Publication No. CN119257543A discloses a photoacoustic imaging method and a photoacoustic imaging system. The method includes: acquiring ultrasonic image data of an imaging target, and determining the type of biological tissue of the imaging target based on the characteristics of the ultrasonic image data; determining a light flux distribution result corresponding to the type of biological tissue, where the light flux distribution result is generated after laser irradiates a simulated target tissue during the simulation imaging process of a simulated target tissue corresponding to the type of biological tissue of the imaging target, and is used to characterize the light flux distribution of the biological tissue relative to the body surface at different imaging depths; acquiring photoacoustic image data of the imaging target; determining a compensation parameter for the photoacoustic image data based on the light flux distribution result; compensating the photoacoustic image data based on the compensation parameter to obtain compensated photoacoustic image data; and generating a photoacoustic image according to the compensated photoacoustic image data. It can be seen that the above technical solution has the following problems: only compensating the photoacoustic image data through the light flux distribution, it is impossible to adaptively adjust the imaging parameters according to the image state, and it is difficult to adapt to complex biological tissue environments and dynamic physiological processes, resulting in poor image quality of photoacoustic imaging. Summary of the Invention

[0004] Therefore, the present invention provides a photoacoustic imaging method and system based on data analysis to overcome the problem in the prior art that it is impossible to adaptively adjust the imaging parameters according to the image state, and it is difficult to adapt to complex biological tissue environments and dynamic physiological processes, resulting in poor image quality of photoacoustic imaging.

[0005] To achieve the above object, the present invention provides a photoacoustic imaging method based on data analysis, including: Determining the imaging state of the initial photoacoustic image according to the image anomaly coefficient and the detection domain variation value, and determining the analysis method according to the imaging state, where the analysis method is abnormal detection analysis or tissue influence analysis; In the abnormal detection analysis, determining a first optimization method according to the noise modulation coefficient, where the first optimization method is to adjust the laser wavelength according to the noise modulation coefficient, or, detection adjustment optimization; In the tissue influence analysis, a second optimization method is determined according to the distortion influence coefficient and the regional dislocation coefficient. The second optimization method is to determine the adjustment method according to the blood flow influence index, or to determine the compensation method according to the tissue difference coefficient and the deviation uniformity. The adjustment method is to adjust according to the blood flow influence index for the detection angle range, or to adjust according to the distortion influence coefficient for the laser repetition frequency. The compensation method is characteristic frequency compensation or related frequency compensation.

[0006] Further, if the imaging state of the initial photoacoustic image is that the image anomaly coefficient is greater than or equal to the preset image anomaly coefficient and the detection domain variation value is less than the preset detection domain variation value, the analysis method is anomaly detection analysis.

[0007] Further, if the imaging state of the initial photoacoustic image is that the image anomaly coefficient is less than the preset image anomaly coefficient or the detection domain variation value is greater than or equal to the preset detection domain variation value, the analysis method is tissue influence analysis.

[0008] Further, if the noise modulation coefficient is greater than or equal to the preset noise modulation coefficient, the first optimization method is to perform a decreasing adjustment on the laser wavelength according to the noise modulation coefficient. The decreasing value of the laser wavelength has a positive correlation with the noise modulation coefficient.

[0009] Further, if the noise modulation coefficient is less than the preset noise modulation coefficient, the first optimization method is detection adjustment optimization. In the detection adjustment optimization, the detection adjustment interval is determined based on the detection anomaly coefficient, and the partitioning method is determined based on the interval comparison coefficient to obtain several divided intervals, and the detection angle range of each divided interval is adjusted according to the interval influence coefficient.

[0010] Further, determining the partitioning method based on the interval comparison coefficient includes: If the interval comparison coefficient is greater than or equal to the preset interval comparison coefficient, the partitioning method is to perform associated partitioning according to the point position correlation degree. If the interval comparison coefficient is less than the preset interval comparison coefficient, the partitioning method is to perform uniform partitioning according to the anomaly coefficient.

[0011] Further, the detection angle range of each divided interval is increased according to the interval influence coefficient; The increasing value of the detection angle range corresponding to a single divided interval has a positive correlation with the interval influence coefficient corresponding to the divided interval.

[0012] Further, if the distortion influence coefficient is less than the preset distortion influence coefficient and the regional dislocation coefficient is less than the preset regional dislocation coefficient, the second optimization method is to determine the compensation method according to the tissue difference coefficient and the deviation uniformity. If the tissue difference coefficient is greater than or equal to the preset tissue difference coefficient or the deviation uniformity is greater than or equal to the preset deviation uniformity, the compensation method is related frequency compensation; If the tissue difference coefficient is less than the preset tissue difference coefficient and the deviation uniformity is less than the preset deviation uniformity, the compensation method is characteristic frequency compensation.

[0013] Furthermore, if the distortion influence coefficient is greater than or equal to the preset distortion influence coefficient or the regional dislocation coefficient is greater than or equal to the preset regional dislocation coefficient, the second optimization method is to determine the adjustment method according to the blood flow influence index; If the blood flow influence index is greater than or equal to the preset blood flow influence index, the adjustment method is to increase the detection angle range according to the blood flow influence index; If the blood flow influence index is less than the preset blood flow influence index, the adjustment method is to increase the laser repetition frequency according to the distortion influence coefficient.

[0014] The present invention also provides a photoacoustic imaging system based on data analysis, including: A state analysis unit for determining the imaging state of the initial photoacoustic image according to the image anomaly coefficient and the detection domain variation value, and determining the analysis method according to the imaging state. The analysis method is anomaly detection analysis or tissue influence analysis; A first optimization unit connected to the state analysis unit for determining a first optimization method according to the noise modulation coefficient in the anomaly detection analysis. The first optimization method is to adjust the laser wavelength according to the noise modulation coefficient, or detection adjustment optimization; A second optimization unit connected to the state analysis unit for determining a second optimization method according to the distortion influence coefficient and the regional dislocation coefficient in the tissue influence analysis. The second optimization method is to determine the adjustment method according to the blood flow influence index, or to determine the compensation method according to the tissue difference coefficient and the deviation uniformity; An adjustment unit connected to the second optimization unit for determining the adjustment method according to the blood flow influence index to adjust the detection angle range according to the blood flow influence index, or to adjust the laser repetition frequency according to the distortion influence coefficient; A compensation unit connected to the second optimization unit for determining the compensation method as characteristic frequency compensation or related frequency compensation according to the tissue difference coefficient and the deviation uniformity.

[0015] Compared with the prior art, the beneficial effect of the present invention lies in that, in the technical solution of the present invention, the imaging state of the initial photoacoustic image is determined according to the image anomaly coefficient and the detection domain variation value, the abnormality degree of the initial photoacoustic image is effectively reflected by the image anomaly coefficient and the detection domain variation value, and then different analysis methods are adaptively selected according to the imaging state, so that the selection of the analysis method is more in line with the actual application scenario, and the abnormal problems in the imaging can be accurately identified and processed, and the artifacts and noise in the image can be significantly reduced, the clarity and resolution of the image are improved, and the uncertainty in the imaging process is reduced, thereby improving the image quality of the photoacoustic imaging.

[0016] Furthermore, in the present invention, the abnormal situation of noise is effectively reflected by the noise modulation coefficient, and then different first optimization methods are adaptively selected according to the noise modulation coefficient, so that the selected first optimization method can improve the signal collection efficiency and quality, thereby enhancing the image contrast, and by adjusting the wavelength according to the noise modulation coefficient, the wavelength that best suits the current imaging target can be selected, thereby improving the adaptability of imaging, reducing signal noise and artifacts, and thus improving image quality.

[0017] Furthermore, in the present invention, by analyzing the detection anomaly coefficient, the abnormal time point in the imaging process can be identified, and the detection adjustment interval can be determined, so that the abnormal area in the imaging can be accurately located and unnecessary adjustments can be reduced. The partitioning method is determined based on the interval comparison coefficient. Through the associated partitioning, the detection parameters can be adjusted more finely to improve the imaging quality. Through the uniform partitioning method, the detection resources can be more evenly distributed to improve the imaging efficiency. By analyzing the interval influence coefficient of each divided interval, the detection angle range can be adjusted in a targeted manner, thereby optimizing the detection parameters and reducing the impact of noise on imaging.

[0018] Furthermore, in the present invention, a second optimization method is determined according to the distortion influence coefficient and the regional dislocation coefficient. The distortion influence coefficient and the regional dislocation coefficient are used to effectively reflect the distortion and regional dislocation problems in the imaging process, and then different second optimization methods are adaptively selected according to the distortion influence coefficient and the regional dislocation coefficient. When the distortion or regional dislocation problem is more serious, the blood flow influence index can quantify the influence of blood flow on the imaging quality. By adjusting the detection angle range or the laser repetition frequency, the imaging parameters can be optimized to reduce the negative impact of blood flow on the imaging quality. When the distortion and regional dislocation problems are not serious, the influence of tissue characteristics on imaging is further analyzed. The tissue difference coefficient and the deviation uniformity can reflect the tissue heterogeneity and the characteristic uniformity of the abnormal area. The imaging quality can be optimized through related frequency compensation or characteristic frequency compensation, and the artifacts and noise in the image can be significantly reduced, thereby improving the clarity and resolution of the image.

[0019] Furthermore, in the present invention, the tissue difference coefficient and the deviation uniformity effectively reflect the significant tissue characteristic differences or the signal distribution uniformity in the imaging region. Then, different compensation methods are adaptively selected according to the tissue difference coefficient and the deviation uniformity. The relevant frequency compensation can effectively improve the frequency characteristics of the signal by adjusting the frequencies with higher frequency coupling degrees, and reduce the image quality problems caused by frequency distortion. The characteristic frequency compensation focuses on optimizing the signals of the characteristic frequencies, further improving the clarity and resolution of the image, and thus significantly improving the quality and efficiency of photoacoustic imaging.

[0020] Furthermore, in the present invention, the adjustment method is determined according to the blood flow influence index. The blood flow influence index effectively reflects the influence degree of blood flow on imaging. When the blood flow influence index is high, it indicates that the influence of blood flow on the imaging quality is significant. By increasing the detection angle range, the artifacts and noises caused by blood flow can be reduced. When the blood flow influence index is small, by increasing the laser repetition frequency, the resolution and sensitivity of imaging can be improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] Figure 1 is a schematic diagram of the photoacoustic imaging method based on data analysis of the present invention; Figure 2 is a flowchart of determining the first optimization method according to the noise modulation coefficient of the present invention; Figure 3 is a flowchart of determining the second optimization method according to the distortion influence coefficient and the regional dislocation coefficient of the present invention; Figure 4 is a unit connection diagram of the photoacoustic imaging system based on data analysis of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0022] In order to make the objectives and advantages of the present invention clearer, the present invention will be further described below in conjunction with embodiments; it should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0023] The preferred embodiments of the present invention will be described below with reference to the accompanying drawings. Those skilled in the art should understand that these embodiments are only used to explain the technical principle of the present invention and do not limit the protection scope of the present invention.

[0024] It should be noted that in the description of the present invention, the terms indicating the direction or position relationship such as "upper", "lower", "left", "right", "inner", "outer", etc. are based on the direction or position relationship shown in the drawings. This is only for convenience of description and does not indicate or imply that the device or element must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation to the present invention.

[0025] In addition, it should be noted that in the description of the present invention, unless otherwise clearly specified and defined, the terms "installation", "connection", and "linkage" should be understood in a broad sense. 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, and it can be the communication inside two components. For those skilled in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0026] Please refer to Figures 1 to 3 As shown, the present invention provides a photoacoustic imaging method based on data analysis, including: Determining the imaging state of the initial photoacoustic image according to the image anomaly coefficient and the detection domain variation value, and determining the analysis method according to the imaging state. The analysis method is abnormal detection analysis or tissue influence analysis; In the abnormal detection analysis, determining the first optimization method according to the noise modulation coefficient. The first optimization method is to adjust the laser wavelength according to the noise modulation coefficient, or, detection adjustment optimization; In the tissue influence analysis, determining the second optimization method according to the distortion influence coefficient and the regional dislocation coefficient. The second optimization method is to determine the adjustment method according to the blood flow influence index, or, to determine the compensation method according to the tissue difference coefficient and the deviation uniformity; The adjustment method is to adjust the detection angle range according to the blood flow influence index, or to adjust the laser repetition frequency according to the distortion influence coefficient. The compensation method is characteristic frequency compensation or related frequency compensation.

[0027] The application scenario of the present invention is to optimize the subsequent imaging process through the initial photoacoustic image to generate a photoacoustic image. In the present invention, a number of historical records are correspondingly set. Any historical record records the sub-image anomaly coefficient, image anomaly coefficient, detection domain variation value, sub-noise modulation coefficient, detection anomaly coefficient, and interval comparison coefficient, etc. in the historical process of generating a photoacoustic image at least once. And each historical record corresponds to a qualified mark. The qualified mark records whether the process of generating a photoacoustic image meets the user's requirements. The qualified mark can be manually recorded. It can be understood that the user can determine whether the process of generating a photoacoustic image meets the requirements according to the self-set index. The self-set index can be, but is not limited to, the artifact area, which will not be elaborated here. Among them, the artifact area is the area of the artifact in the generated photoacoustic image.

[0028] The initial photoacoustic image is the photoacoustic image generated by irradiating the target tissue with pulsed laser for the first time; The photoacoustic imaging steps include: irradiating the target tissue with a pulsed laser, detecting the ultrasonic signals propagating from the tissue surface through an ultrasonic transducer, and using a reconstruction algorithm to invert the light absorption distribution inside the tissue based on the received ultrasonic signals to generate a photoacoustic image. The reconstruction algorithm includes, but is not limited to, the back-projection algorithm and the iterative reconstruction algorithm. The target tissue includes, but is not limited to, tumor tissue, cardiovascular tissue, nervous system tissue, lymphatic system tissue, and skin system tissue. This is easily understood by those skilled in the art and will not be elaborated herein.

[0029] Specifically, if the imaging state of the initial photoacoustic image is that the image anomaly coefficient is greater than or equal to the preset image anomaly coefficient and the detection domain variation value is less than the preset detection domain variation value, the analysis method is anomaly detection analysis.

[0030] Among them, the imaging state includes a first imaging state and a second imaging state. The first imaging state is that the image anomaly coefficient is greater than or equal to the preset image anomaly coefficient and the detection domain variation value is less than the preset detection domain variation value. The second imaging state is that the image anomaly coefficient is less than the preset image anomaly coefficient or the detection domain variation value is greater than or equal to the preset detection domain variation value. The image anomaly coefficient = |the sub-image anomaly coefficient corresponding to the initial photoacoustic image - the average value of the sub-image anomaly coefficients corresponding to the photoacoustic images generated by the historical records that can meet the user's needs|. The sub-image anomaly coefficient = the standard deviation of the pixel values corresponding to each pixel point in the initial photoacoustic image / the average value of the pixel values corresponding to each pixel point in the initial photoacoustic image. After irradiating the tissue with the pulsed laser, the time period during which the tissue generates photoacoustic signals and is received by the detector is recorded as the full interval. The full interval is equally divided into a parts. The value of a is positively correlated with the duration of the full interval. The time point with the latest time among each equal division point and the reference time point is recorded as the time point. The optical power corresponding to a single time point is measured by using an optical power meter or an energy meter. Among them, a is the smallest integer greater than or equal to a1, and a1 = 0.1 × the duration of the full interval. The unit of the duration of the full interval is s. The detection domain variation value = (the optical flux corresponding to the earliest time point in the full interval - the optical flux corresponding to the latest time point in the full interval) / the standard deviation of the optical fluxes corresponding to each time point in the full interval. The values of the preset image anomaly coefficient and the preset detection domain variation value can be determined by the user according to the actual application scenario. The smaller the value of the preset image anomaly coefficient and the larger the value of the preset detection domain variation value, the greater the user's need for detecting anomaly analysis. Provide a set of values for the preset image anomaly coefficient and the preset detection domain variation value, detect the historical records of the user's anomaly analysis detection, and record the average value of the image anomaly coefficients corresponding to the historical records that can meet the user's needs as the preset image anomaly coefficient, and record the average value of the detection domain variation values corresponding to the historical records that can meet the user's needs as the preset detection domain variation value.

[0031] Specifically, if the imaging state of the initial photoacoustic image is that the image anomaly coefficient is less than the preset image anomaly coefficient or the detection domain variation value is greater than or equal to the preset detection domain variation value, the analysis method is tissue influence analysis.

[0032] Specifically, if the noise modulation coefficient is greater than or equal to the preset noise modulation coefficient, the first optimization method is to reduce the laser wavelength according to the noise modulation coefficient; The reduction value of the laser wavelength has a positive correlation with the noise modulation coefficient.

[0033] Among them, the noise modulation coefficient is the maximum value of the sub-noise modulation coefficients corresponding to the initial photoacoustic image and each photoacoustic image generated from the historical records that can meet the user's needs; The calculation formula for the sub-noise modulation coefficient Z corresponding to a single photoacoustic image generated from the historical records that can meet the user's needs and the initial photoacoustic image is: ; Among them, I(α,β) and K(α,β) are the pixel values of a single photoacoustic image generated from the historical records that can meet the user's needs and the initial photoacoustic image at the position (α,β), respectively, and m and ε are the number of rows and columns of the pixels included in the image; The reduction value of the laser wavelength = (K / preset noise modulation coefficient) × noise modulation coefficient, where K is a proportionality coefficient, K = 0.9nm; The value of the preset noise modulation coefficient can be determined by the user according to the actual application scenario. The smaller the value of the preset noise modulation coefficient, the greater the user's need to reduce the laser wavelength according to the noise modulation coefficient. Provide a value for the preset noise modulation coefficient, detect the historical records of the user's reduction adjustment of the laser wavelength according to the noise modulation coefficient, and record the average value of the noise modulation coefficients corresponding to the historical records that can meet the user's needs as the preset noise modulation coefficient; The laser wavelength is the wavelength of the pulsed laser used.

[0034] Specifically, if the noise modulation coefficient is less than the preset noise modulation coefficient, the first optimization method is detection adjustment optimization; During detection adjustment optimization, a detection adjustment interval is determined based on the detection anomaly coefficient, and a partitioning method is determined based on the interval comparison coefficient to obtain several divided intervals, and the detection angle range of each divided interval is adjusted according to the interval influence coefficient.

[0035] Among them, when determining the detection adjustment interval based on the detection anomaly coefficient, the time points with the detection anomaly coefficient greater than the preset detection anomaly coefficient are recorded as abnormal time points, and the smallest time period that can include each abnormal time point is recorded as the detection adjustment interval; The confirmation method of the detection anomaly coefficient is that for a single time point, the detection anomaly coefficient corresponding to this time point = |the luminous flux corresponding to this time point - the average value of the luminous fluxes corresponding to each time point in the full interval|, For the value of the preset detection anomaly coefficient, the user can determine it according to the actual application scenario. The greater the user's demand for improving the imaging quality, the smaller the value of the preset detection anomaly coefficient. Provide a value of the preset detection anomaly coefficient, and record the minimum value of the detection anomaly coefficients corresponding to each abnormal time point in the historical records that can meet the user's needs as the preset detection anomaly coefficient.

[0036] Specifically, determining the partitioning method based on the interval comparison coefficient includes: If the interval comparison coefficient is greater than or equal to the preset interval comparison coefficient, the partitioning method is to perform associated partitioning according to the point position correlation; If the interval comparison coefficient is less than the preset interval comparison coefficient, the partitioning method is to perform uniform partitioning according to the anomaly coefficient.

[0037] Among them, the interval comparison coefficient is the standard deviation of the luminous fluxes corresponding to each abnormal time point in the detection adjustment interval; For the value of the preset interval comparison coefficient, the user can determine it according to the actual application scenario. The greater the user's demand for improving the imaging quality, the smaller the value of the preset interval comparison coefficient. Provide a value of the preset interval comparison coefficient, detect the historical records of the user's uniform partitioning according to the anomaly coefficient, and record the average value of the interval comparison coefficients corresponding to the historical records that can meet the user's needs as the preset interval comparison coefficient; Record each time point in the detection adjustment interval as a range point; When performing associated partitioning according to the point position correlation, the point position correlation between any two range points in a single divided interval is greater than or equal to the preset point position correlation, and there are adjacent time points for any one range point in a single divided interval, and the point position correlation between any edge time point corresponding to a single divided interval and any one range point in this divided interval is less than the preset point position correlation; For a single range point, denote this range point as the target range point. The adjacent time points corresponding to the target range point are the range points that are adjacent to the target range point and in the same division interval as the target range point. The edge time points corresponding to a single division interval are the time points that are adjacent to the range points in this division interval and not in this division interval. When performing uniform partitioning according to the anomaly coefficient, divide the detection adjustment interval into several division intervals with the same time length. The number of division intervals is negatively correlated with the anomaly coefficient, and the anomaly coefficient is the time length of the detection adjustment interval. The number of division intervals is the smallest integer greater than or equal to M, where M = (the number of anomaly time points included in the detection adjustment interval / the anomaly coefficient) × the total number of range points included in the detection adjustment interval × 0.5. For any two time points, the point position correlation degree = 1 - (the absolute value of the difference between the anomaly reference values corresponding to the two time points / the larger value of the anomaly reference values corresponding to the two time points). For a single time point, denote this time point as the target time point, detect the historical time points with the same detection anomaly coefficient as this time point, and denote the average value of the sub - anomaly reference values corresponding to each historical time point as the anomaly reference value corresponding to the target time point. Denote the time points in the detection adjustment interval corresponding to each initial photoacoustic image in the historical records that can meet the user's needs as historical time points. The confirmation method of the sub - anomaly reference value is as follows: for a single historical time point, denote this historical time point as the target historical time point, detect the initial photoacoustic image corresponding to the target historical time point. The sub - anomaly reference value corresponding to the target historical time point = the image anomaly coefficient of the initial photoacoustic image corresponding to the target historical time point × [the order in which the target historical time point is located in the detection adjustment interval corresponding to the target historical time point in ascending order / (the total number of historical time points in the detection adjustment interval corresponding to the target historical time point + 1)]. The value of the preset point position correlation degree can be determined by the user according to the actual application scenario. The greater the user's demand for improving the imaging accuracy, the greater the value of the preset point position correlation degree. Provide a value of the preset point position correlation degree, detect the historical records of the user's association partitioning according to the point position correlation degree, and denote the average value of the reference point position correlation degrees corresponding to each division interval in the historical records that can meet the user's needs as the preset point position correlation degree. The reference point position correlation degree is the point position correlation degree corresponding to any two time points in a single division interval in the historical records that can meet the user's needs.

[0038] Specifically, increase and adjust the detection angle range for each division interval according to the interval influence coefficient. The increase value of the detection angle range corresponding to a single division interval is positively correlated with the interval influence coefficient corresponding to this division interval.

[0039] Among them, the way to confirm the interval influence coefficient is as follows: for a single divided interval, denote this interval as the target divided interval, and the interval influence coefficient corresponding to the target divided interval = [|the first standard deviation - the second standard deviation| / the larger value of the first standard deviation and the second standard deviation] - [the number of frequencies in the target power spectrum that are the same as those in the reference power spectrum / (the number of frequencies appearing in the reference power spectrum + 1)]; Denote the standard deviation of the power corresponding to each frequency in the target power spectrum as the first standard deviation, and denote the standard deviation of the power corresponding to each frequency in the reference power spectrum as the second standard deviation. Denote the power spectrum corresponding to the target divided interval as the target power spectrum, denote the time period when the ultrasonic signal propagated from the tissue surface is detected by the ultrasonic transducer as the full interval, and denote the power spectrum corresponding to the time period in the full interval except the detection adjustment interval as the reference power spectrum; Convert the ultrasonic signal of a single interval into a frequency-domain signal through the Fourier transform method, so as to obtain the power spectrum of the signal, and the power spectrum reflects the power distribution of the signal at different frequencies; The detection angle range corresponding to a single divided interval is the angle interval at which the ultrasonic transducer can receive the photoacoustic signal; The increment value of the detection angle range corresponding to a single divided interval The calculation formula is: ; Among them, C is the interval influence coefficient corresponding to this divided interval, s1 is the first coefficient, s2 is the second coefficient, s1 = 1.353, s2 = 2; When performing an increase adjustment on the detection angle range, the increase amount extends the detection angle range from both sides of the initial detection angle range based on the initial detection angle range and the increase amounts on both sides are the same, and the initial detection angle range is 45°.

[0040] Specifically, if the distortion influence coefficient is less than the preset distortion influence coefficient and the regional dislocation coefficient is less than the preset regional dislocation coefficient, then the second optimization method is to determine the compensation method according to the tissue difference coefficient and the deviation uniformity; If the tissue difference coefficient is greater than or equal to the preset tissue difference coefficient or the deviation uniformity is greater than or equal to the preset deviation uniformity, then the compensation method is related frequency compensation; If the tissue difference coefficient is less than the preset tissue difference coefficient and the deviation uniformity is less than the preset deviation uniformity, then the compensation method is characteristic frequency compensation.

[0041] Among them, the distortion influence coefficient = |the image reference value corresponding to the initial photoacoustic image - the average value of the image reference values corresponding to each photoacoustic image generated by the historical records that can meet the user's needs|; Image reference value corresponding to a single photoacoustic image = Standard deviation of the power corresponding to each frequency in the power spectrum of the entire interval corresponding to the photoacoustic image / Average value of the power corresponding to each frequency in the power spectrum of the entire interval corresponding to the photoacoustic image; The regional dislocation coefficient is the average value of the distance reference values corresponding to each abnormal region. For a single abnormal region, the distance reference value corresponding to the abnormal region is the average value of the shortest distances from the abnormal region to other abnormal regions; Divide the initial photoacoustic image into several rectangular regions with the same and equal area. The number of rectangular regions is positively correlated with the area of the initial photoacoustic image; Number of rectangular regions = (Area of the initial photoacoustic image / S) × K, where S is the correction coefficient, S = 1 square centimeter, the unit of the area of the initial photoacoustic image is square centimeter, and K is the weight coefficient, K = 10; Mark the rectangular regions with characteristic coefficients greater than the preset characteristic coefficient as abnormal regions; Characteristic coefficient of a single rectangular region = |Average value of the pixel values corresponding to each pixel point in the rectangular region - Average value of the pixel values corresponding to each pixel point in the initial photoacoustic image|; For the value of the preset characteristic coefficient, the user can determine it according to the actual application scenario. The greater the user's accuracy requirement for improving the quality of photoacoustic imaging, the smaller the value of the preset characteristic coefficient. Provide a value of the preset characteristic coefficient, and record the average value of the characteristic coefficients corresponding to each abnormal region in the historical records that can meet the user's needs as the preset characteristic coefficient; For the values of the preset distortion influence coefficient and the preset regional dislocation coefficient, the user can determine them according to the actual application scenario. The greater the values of the preset distortion influence coefficient and the preset regional dislocation coefficient, the greater the user's need to determine the compensation method according to the tissue difference coefficient and the deviation uniformity. Provide the values of the preset distortion influence coefficient and the preset regional dislocation coefficient, detect the historical records of the user determining the compensation method according to the tissue difference coefficient and the deviation uniformity, and record the average value of the distortion influence coefficients corresponding to the historical records that can meet the user's needs as the preset distortion influence coefficient, and record the average value of the regional dislocation coefficients corresponding to the historical records that can meet the user's needs as the preset regional dislocation coefficient; The tissue difference coefficient is the average value of the sub - tissue difference degrees corresponding to each rectangular region. The confirmation method of the sub - tissue difference degree corresponding to a single rectangular region is as follows: for a single rectangular region, mark the connected domain that appears in the rectangular region as the reference connected domain. The sub - tissue difference degree corresponding to the rectangular region = 1 / (Average value of the regional emergence coefficients corresponding to each reference connected domain + 1), and the regional emergence coefficient corresponding to a single reference connected domain is the number of rectangular regions where the reference connected domain appears; Each connected region contains a number of pixel points. The pixel values corresponding to the pixel points in a single connected region are all the same, and the pixel values of the adjacent pixel points corresponding to the connected region are all different from the pixel value corresponding to the connected region. The adjacent pixel points corresponding to a single connected region are the pixel points adjacent to the connected region; Deviation uniformity = 1 / standard deviation of the characteristic coefficients corresponding to each abnormal region; The values of the preset tissue difference coefficient and the preset deviation uniformity can be determined by the user according to the actual application scenario. The smaller the values of the preset tissue difference coefficient and the preset deviation uniformity, the greater the user's need for relevant frequency compensation. Provide a set of values for the preset tissue difference coefficient and the preset deviation uniformity, detect the historical records of the user's relevant frequency compensation, and record the average value of the tissue difference coefficients corresponding to the historical records that can meet the user's needs as the preset tissue difference coefficient, and record the average value of the deviation uniformities corresponding to the historical records that can meet the user's needs as the preset deviation uniformity; Each time period obtained by equally dividing the full interval into a parts is recorded as a sub-interval; The power spectrum corresponding to the full interval is recorded as the full power spectrum, and the power spectra corresponding to each sub-interval are recorded as sub-power spectra; In relevant frequency compensation, the frequencies in the full power spectrum with a pre-compensation coefficient greater than the preset pre-compensation coefficient are recorded as pre-compensation frequencies. For a single frequency, this frequency is recorded as the target frequency. The pre-compensation coefficient corresponding to the target frequency = |the compensation reference value corresponding to the target frequency in the initial photoacoustic image - the average value of the compensation reference values corresponding to the photoacoustic images generated in each historical record that can meet the user's needs|. The compensation reference value corresponding to the target frequency in a single photoacoustic image is the standard deviation of the power corresponding to the power spectra of each sub-interval corresponding to the target frequency in this photoacoustic image; The other frequencies except the pre-compensation frequencies are recorded as candidate frequencies. Compensate for the relevant frequencies corresponding to each pre-compensation frequency and each pre-compensation frequency; For a single pre-compensation frequency, this pre-compensation frequency is recorded as the target pre-compensation frequency, and the relevant frequency corresponding to the target pre-compensation frequency is the candidate frequency whose frequency coupling degree with the target pre-compensation frequency is greater than the preset frequency coupling degree; For any two frequencies, the calculation formula for the frequency coupling degree r is: ; where n is the number of sub-intervals corresponding to the initial photoacoustic image; and are respectively the powers corresponding to the two frequencies in the power spectrum of the i-th sub-interval, is the average value of the powers corresponding to the corresponding frequency in the power spectra of each sub-interval, is The average value of the power corresponding to the corresponding frequency in the power spectrum of each sub-interval, i = 1, 2, 3, ……, n.

[0042] For the values of the preset pre-compensation coefficient and the preset frequency coupling degree, the user can determine them according to the actual application scenario. The greater the user's demand for improving the photoacoustic imaging accuracy, the smaller the value of the preset pre-compensation coefficient and the greater the value of the preset frequency coupling degree. Provide a set of values for the preset pre-compensation coefficient and the preset frequency coupling degree. Denote the average value of the pre-compensation coefficients corresponding to each pre-compensation frequency in the historical records that can meet the user's needs as the preset pre-compensation coefficient, and the preset frequency coupling degree is 0.8; When compensating for a single frequency, denote this frequency as the target frequency, and adjust the power corresponding to the target frequency in the full power spectrum according to the frequency comparison deviation value; If the frequency comparison deviation value corresponding to the target frequency is greater than or equal to the preset frequency comparison deviation value, then reduce and adjust the power corresponding to the target frequency in the full power spectrum according to the compensation coefficient; If the frequency comparison deviation value corresponding to the target frequency is less than the preset frequency comparison deviation value, then increase and adjust the power corresponding to the target frequency in the full power spectrum according to the compensation coefficient; Both the increase value and the decrease value of the power corresponding to the target frequency are positively correlated with the compensation coefficient corresponding to the target frequency; the compensation coefficient corresponding to the target frequency = |the pre-compensation coefficient corresponding to the target frequency × (the frequency comparison deviation value corresponding to the target frequency + 1)|; The increase value or the decrease value of the power corresponding to the target frequency = the compensation coefficient corresponding to the target frequency × k1 × the power corresponding to the target frequency in the full power spectrum, K1 = 0.4; The frequency comparison deviation value = (the power corresponding to the target frequency in the full power spectrum corresponding to the initial photoacoustic image / the frequency average value corresponding to the initial photoacoustic image) - (the average value of the powers corresponding to the target frequency in the full power spectra corresponding to the photoacoustic images generated in each historical record that can meet the user's needs / the average value of the frequency average values corresponding to each full power spectrum in the photoacoustic images generated in each historical record that can meet the user's needs); For the value of the preset frequency comparison deviation value, the user can determine it according to the actual application scenario. The greater the value of the preset frequency comparison deviation value, the greater the user's demand for increasing and adjusting the power corresponding to the frequency. Provide a set of values for the preset frequency comparison deviation value, and the preset frequency comparison deviation value is 0; The frequency average value is the average value of the powers corresponding to each frequency in a single full spectrum; In the characteristic frequency compensation: Denote the sub-intervals with the sub-interval turbulence coefficient greater than or equal to the preset sub-interval turbulence coefficient as characteristic sub-intervals, denote the sub-intervals with the sub-interval turbulence coefficient less than the preset sub-interval turbulence coefficient as non-characteristic sub-intervals, and denote each frequency in the characteristic sub-intervals as characteristic frequencies; When compensating for a single characteristic frequency, the characteristic frequency is denoted as the target characteristic frequency, and the power corresponding to the target characteristic frequency in the full power spectrum is adjusted according to the characteristic amplitude; If the characteristic amplitude corresponding to the target characteristic frequency is greater than or equal to the preset characteristic amplitude, the power corresponding to the target characteristic frequency in the full power spectrum is adjusted to decrease according to the characteristic compensation degree; If the characteristic amplitude corresponding to the target characteristic frequency is less than the preset characteristic amplitude, the power corresponding to the target characteristic frequency in the full power spectrum is adjusted to increase according to the characteristic compensation degree; It can be understood that by adjusting the corresponding frequencies in the full power spectrum, the image quality can be improved; after adjusting the power corresponding to the frequencies in the full power spectrum, the inverse Fourier transform is performed on the adjusted spectrum to obtain an optimized ultrasonic signal. Using a reconstruction algorithm, the light absorption distribution inside the tissue is inverted according to the received ultrasonic signal to generate a photoacoustic image, which is easy for those skilled in the art to understand and will not be elaborated here; Characteristic compensation degree = |average value of the turbulence coefficients of each sub-interval in the first interval corresponding to the target characteristic frequency × (characteristic amplitude corresponding to the target characteristic frequency + 1)|; Both the increase value and the decrease value of the power corresponding to the target characteristic frequency are positively correlated with the characteristic compensation degree corresponding to the target characteristic frequency; Increase value or decrease value of the power corresponding to the target characteristic frequency = (characteristic compensation degree corresponding to the target characteristic frequency / power corresponding to the target characteristic frequency) × k2 × power corresponding to the target characteristic frequency in the full power spectrum, k2 = 0.3; The sub-interval turbulence coefficient is the average value of the powers corresponding to each frequency in the power spectrum corresponding to the sub-interval; For the value of the preset sub-interval turbulence coefficient, the user can determine it according to the actual application scenario. The larger the value of the preset sub-interval turbulence coefficient, the greater the user's demand for increasing the power corresponding to the frequency. A value of the preset sub-interval turbulence coefficient is provided, and the preset sub-interval turbulence coefficient is 0; Characteristic amplitude corresponding to the target characteristic frequency = average value of the powers corresponding to the target characteristic frequency in the power spectra corresponding to each first interval - average value of the powers corresponding to the target characteristic frequency in the power spectra corresponding to each second interval; The characteristic sub-interval where the target characteristic frequency appears is denoted as the first interval, and the non-characteristic sub-interval where the target characteristic frequency appears is denoted as the second interval.

[0043] Specifically, if the distortion influence coefficient is greater than or equal to the preset distortion influence coefficient or the region dislocation coefficient is greater than or equal to the preset region dislocation coefficient, the second optimization method is to determine the adjustment method according to the blood flow influence index; If the blood flow influence index is greater than or equal to the preset blood flow influence index, the adjustment method is to increase the detection angle range according to the blood flow influence index; If the blood flow influence index is less than the preset blood flow influence index, the adjustment method is to increase the laser repetition frequency according to the distortion influence coefficient.

[0044] Among them, the blood flow influence index is the maximum value of the broadening reference values corresponding to each sub-interval of the initial photoacoustic image. The broadening reference value corresponding to a single sub-interval is the number of frequencies appearing in the power spectrum corresponding to this sub-interval. It can be understood that the turbulence or uneven velocity distribution of the blood flow will cause changes in the broadening reference value; For the value of the preset blood flow influence index, the user can determine it according to the actual application scenario. The smaller the value of the preset blood flow influence index, the greater the user's need to increase the detection angle according to the blood flow influence index. Provide a value of the preset blood flow influence index, detect the historical record of the user increasing the detection angle according to the blood flow influence index, and record the average value of the blood flow influence indexes corresponding to each historical record that can meet the user's needs as the preset blood flow influence index; Increase the detection angle range of the entire interval according to the blood flow influence index. The increase value of the detection angle range is positively correlated with the blood flow influence index. When increasing the detection angle range, the increase amount expands the detection angle range from both sides based on the initial detection angle range and the increase amounts on both sides are the same. The initial detection angle range is 45°; The increase value of the detection angle range = (blood flow influence index / preset blood flow influence index) × initial detection angle range; Increase the laser repetition frequency according to the distortion influence coefficient. The increase value of the laser repetition frequency is positively correlated with the distortion influence coefficient; The increase value of the laser repetition frequency = initial laser repetition frequency × ln (distortion influence coefficient + 1), and the initial laser repetition frequency is 2 pulses / s; The laser repetition frequency is the number of laser pulses output per unit time.

[0045] Please refer to Figure 4 As shown, it is the unit connection diagram of the photoacoustic imaging system based on data analysis of the present invention. The present invention also provides a photoacoustic imaging system based on data analysis, including: A state analysis unit for determining the imaging state of the initial photoacoustic image according to the image anomaly coefficient and the detection domain variation value, and determining the analysis method according to the imaging state. The analysis method is anomaly detection analysis or tissue influence analysis; The first optimization unit is connected to the state analysis unit and is used to determine a first optimization method according to the noise modulation coefficient in the detection anomaly analysis. The first optimization method is to adjust the laser wavelength according to the noise modulation coefficient, or to perform detection adjustment optimization; The second optimization unit is connected to the state analysis unit and is used to determine a second optimization method according to the distortion influence coefficient and the regional dislocation coefficient in the tissue influence analysis. The second optimization method is to determine the adjustment method according to the blood flow influence index, or to determine the compensation method according to the tissue difference coefficient and the deviation uniformity; The adjustment unit is connected to the second optimization unit and is used to determine that the adjustment method according to the blood flow influence index is to adjust the detection angle range according to the blood flow influence index, or to adjust the laser repetition frequency according to the distortion influence coefficient; The compensation unit is connected to the second optimization unit and is used to determine that the compensation method according to the tissue difference coefficient and the deviation uniformity is characteristic frequency compensation or correlation frequency compensation.

[0046] So far, the technical solution of the present invention has been described in combination with the preferred embodiments shown in the drawings. However, it is easy for those skilled in the art to understand that the protection scope of the present invention is obviously not limited to these specific embodiments. Without departing from the principle of the present invention, those skilled in the art can make equivalent changes or replacements to the relevant technical features, and the technical solutions after these changes or replacements will fall within the protection scope of the present invention.

[0047] The above are only the preferred embodiments of the present invention and are not used to limit the present invention; for those skilled in the art, the present invention can have various changes and modifications. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A photoacoustic imaging method based on data analysis, characterized in that: include: Determine the imaging state of the initial photoacoustic image according to the image abnormality coefficient and the detection domain variation value, and determine the analysis method according to the imaging state, the analysis method being detection abnormality analysis or tissue impact analysis; In the abnormality detection analysis, a first optimization method is determined according to the noise modulation coefficient, and the first optimization method is to adjust the laser wavelength according to the noise modulation coefficient, or to detect the adjustment optimization; In the tissue impact analysis, the second optimization method is determined according to the distortion impact coefficient and the regional dislocation coefficient. The second optimization method is to determine the adjustment method according to the blood flow impact index, or to determine the compensation method according to the tissue difference coefficient and the deviation uniformity. The adjustment method is to adjust the detection angle range according to the blood flow influence index, or to adjust the laser repetition frequency according to the distortion influence coefficient, and the compensation method is characteristic frequency compensation or related frequency compensation.

2. The photoacoustic imaging method based on data analysis according to claim 1, characterized in that: If the imaging state of the initial photoacoustic image is that the image anomaly coefficient is greater than or equal to the preset image anomaly coefficient and the detection domain change value is less than the preset detection domain change value, the analysis method is detection anomaly analysis.

3. The photoacoustic imaging method based on data analysis according to claim 2, characterized in that: If the imaging state of the initial photoacoustic image is that the image abnormality coefficient is less than the preset image abnormality coefficient or the detection domain variation value is greater than or equal to the preset detection domain variation value, the analysis method is tissue impact analysis.

4. The photoacoustic imaging method based on data analysis according to claim 2, characterized in that: If the noise modulation coefficient is greater than or equal to the preset noise modulation coefficient, the first optimization method is to reduce and adjust the laser wavelength according to the noise modulation coefficient; The reduction value of the laser wavelength is positively correlated with the noise modulation coefficient.

5. The photoacoustic imaging method based on data analysis according to claim 4, characterized in that: If the noise modulation coefficient is less than the preset noise modulation coefficient, the first optimization method is detection adjustment optimization; In the detection adjustment optimization, the detection adjustment interval is determined based on the detection anomaly coefficient, and the partitioning method is determined based on the interval comparison coefficient to obtain several divided intervals. The detection angle range of each divided interval is adjusted according to the interval influence coefficient.

6. The photoacoustic imaging method based on data analysis according to claim 5, characterized in that: The partitioning method is determined based on the interval comparison coefficient, including: If the interval comparison coefficient is greater than or equal to the preset interval comparison coefficient, the partitioning method is to perform associated partitioning according to the point correlation degree; If the interval comparison coefficient is less than the preset interval comparison coefficient, the partitioning method is uniform partitioning according to the abnormal coefficient.

7. The photoacoustic imaging method based on data analysis according to claim 5, characterized in that: According to the interval influence coefficient, the detection angle range of each divided interval is increased and adjusted; The increase value of the detection angle range corresponding to a single divided interval is positively correlated with the interval influence coefficient corresponding to the divided interval.

8. The photoacoustic imaging method based on data analysis according to claim 3, characterized in that: If the distortion influence coefficient is less than the preset distortion influence coefficient and the regional displacement coefficient is less than the preset regional displacement coefficient, the second optimization method is to determine the compensation method according to the tissue difference coefficient and the deviation uniformity; If the tissue difference coefficient is greater than or equal to the preset tissue difference coefficient or the deviation uniformity is greater than or equal to the preset deviation uniformity, the compensation method is the related frequency compensation; If the tissue difference coefficient is less than the preset tissue difference coefficient and the deviation uniformity is less than the preset deviation uniformity, the compensation method is characteristic frequency compensation.

9. The photoacoustic imaging method based on data analysis according to claim 8, characterized in that: If the distortion influence coefficient is greater than or equal to the preset distortion influence coefficient or the regional displacement coefficient is greater than or equal to the preset regional displacement coefficient, the second optimization method is to determine the adjustment method according to the blood flow influence index; If the blood flow influence index is greater than or equal to the preset blood flow influence index, the adjustment method is to increase the detection angle range according to the blood flow influence index; If the blood flow influence index is less than the preset blood flow influence index, the adjustment method is to increase the laser repetition frequency according to the distortion influence coefficient.

10. An imaging system using the photoacoustic imaging method based on data analysis according to any one of claims 1 to 9, characterized in that: include: A state analysis unit, used to determine the imaging state of the initial photoacoustic image according to the image abnormality coefficient and the detection domain variation value, and determine the analysis method according to the imaging state, the analysis method being detection abnormality analysis or tissue impact analysis; A first optimization unit, connected to the state analysis unit, is used to determine a first optimization method according to the noise modulation coefficient in the abnormality detection analysis, wherein the first optimization method is to adjust the laser wavelength according to the noise modulation coefficient, or to detect the adjustment optimization; A second optimization unit, which is connected to the state analysis unit, is used to determine a second optimization method according to the distortion influence coefficient and the regional dislocation coefficient in the tissue impact analysis, wherein the second optimization method is to determine an adjustment method according to a blood flow influence index, or to determine a compensation method according to a tissue difference coefficient and a deviation uniformity; an adjustment unit connected to the second optimization unit, and configured to determine, according to the blood flow influence index, an adjustment mode, namely, adjusting the detection angle range according to the blood flow influence index, or adjusting the laser repetition frequency according to the distortion influence coefficient; A compensation unit is connected to the second optimization unit and is used to determine whether the compensation method is characteristic frequency compensation or related frequency compensation according to the tissue difference coefficient and the deviation uniformity.

Images