Photoacoustic imaging system, method for processing photoacoustic signal and storage medium

Abstract

The embodiment of the invention provides a photoacoustic imaging system, a method for processing a photoacoustic signal and a storage medium. The method comprises the following steps: acquiring an oxyhemoglobin saturation image of an imaging target; a first target detection area is determined on the blood oxygen saturation image, the first target detection area comprises a target blood flow blood oxygen signal area, and the target blood flow blood oxygen signal area is provided with a first position; a target blood flow blood oxygen signal area is detected from the blood oxygen saturation image every other preset time period, and a second position of the target blood flow blood oxygen signal area is obtained; calculating a first position deviation value of the target blood flow and blood oxygen signal area according to the first position and the second position; and adjusting the first target detection area based on the first position deviation value.

Description

Technical Field

[0001] The embodiments of the present invention relate to the field of photoacoustics, and in particular to a photoacoustic imaging system, a method for processing photoacoustic signals, and a storage medium. Background Art

[0002] Blood oxygen saturation parameters are one of the important parameters provided by photoacoustic imaging in clinical diagnosis and treatment. The dual-wavelength photoacoustic signals collected during photoacoustic imaging can obtain a blood oxygen saturation image after formula calculation and post-processing, in which the values ​​of blood oxygen saturation in different positions can be distinguished by color. When quantitatively analyzing the blood oxygen saturation in a specific area of ​​the blood oxygen saturation image, the doctor frames the region of interest (ROI) and confirms it, and then the photoacoustic imaging system measures the mean, maximum or minimum value of the blood oxygen saturation in the area. During the measurement process, due to changes in the content and position of the imaging section caused by the patient's breathing or the shaking of the photoacoustic probe, the target blood flow signal in the initially framed ROI will deviate from the ROI, resulting in inaccurate measurement results. To solve the above problems, doctors need to continuously adjust the position of the target blood flow signal according to the real-time feedback of the blood oxygen saturation image when performing photoacoustic imaging so that it is within the ROI. On the one hand, continuous attention to the blood oxygen saturation image for a long time will affect the doctor's observation of the target blood flow signal; on the other hand, cumbersome manual operations will cause the problem of low intelligence of signal detection. Summary of the invention

[0003] The embodiments of the present invention provide a photoacoustic imaging system and a method and storage medium for processing photoacoustic signals, which can reduce the time doctors pay attention to blood oxygen saturation images to reduce the impact of doctors observing target blood flow signals and improve the intelligence of signal detection.

[0004] The technical solution of the present invention is achieved in this way:

[0005] In a first aspect, an embodiment of the present invention provides a method for processing a photoacoustic signal, the method comprising:

[0006] Acquire a blood oxygen saturation image of an imaging target;

[0007] Determine a first target detection area on the blood oxygen saturation image, wherein the first target detection area includes a target blood flow oxygen signal area, and the target blood flow oxygen signal area has a first position;

[0008] At every preset time period, detecting the target blood flow oxygen signal area from the blood oxygen saturation image, and obtaining a second position of the target blood flow oxygen signal area;

[0009] Calculate a first position offset value of the target blood flow oxygen signal area according to the first position and the second position;

[0010] The first object detection area is adjusted based on the first position offset value.

[0011] In a second aspect, an embodiment of the present application provides a method for processing a photoacoustic signal, the method comprising:

[0012] acquiring a photoacoustic image of an imaging target;

[0013] determining a second target detection region on the photoacoustic image, wherein the second target detection region includes a target photoacoustic signal, and the target photoacoustic signal has a third position;

[0014] At every preset time period, detecting the target photoacoustic signal region from the photoacoustic image, and obtaining a fourth position of the target photoacoustic signal region;

[0015] Calculating a second position offset value of the target photoacoustic signal area according to the third position and the fourth position;

[0016] The second object detection area is adjusted based on the second position offset value.

[0017] In a third aspect, an embodiment of the present invention provides a photoacoustic imaging system, characterized in that the photoacoustic imaging system comprises:

[0018] A laser emitting device, wherein the laser generating device is used to generate laser light for irradiating a target object;

[0019] optical fiber;

[0020] Ultrasonic transducer;

[0021] A scanning control device, wherein the scanning control device controls the laser emitting device to generate the laser;

[0022] a receiving circuit, wherein the receiving circuit receives the photoacoustic signal generated by the target object through the ultrasonic transducer to obtain photoacoustic data;

[0023] a processor, wherein the processor processes the photoacoustic data to obtain a blood oxygen saturation image;

[0024] a display, wherein the display displays the blood oxygen saturation image;

[0025] The processor further performs the following steps:

[0026] Acquire a blood oxygen saturation image of an imaging target; determine a target detection area on the blood oxygen saturation image, wherein the target detection area includes a target blood flow oxygen signal area, and the target blood flow oxygen signal area has a first position; detect the target blood flow oxygen signal area from the blood oxygen saturation image at intervals of a preset time period, and obtain a second position of the target blood flow oxygen signal area; calculate a position offset value of the target blood flow oxygen signal area according to the first position and the second position; and adjust the target detection area based on the position offset value.

[0027] In a fourth aspect, an embodiment of the present invention provides a storage medium having a computer program stored thereon, which is applied to a photoacoustic imaging system. When the computer program is executed by a processor, the above-mentioned method for processing photoacoustic signals is implemented.

[0028] An embodiment of the present invention provides a photoacoustic imaging system and a method for processing photoacoustic signals, and a storage medium, the method comprising: acquiring a blood oxygen saturation image of an imaging target; determining a first target detection area on the blood oxygen saturation image, wherein the first target detection area includes a target blood flow blood oxygen signal area, and the target blood flow blood oxygen signal area has a first position; detecting the target blood flow blood oxygen signal area from the blood oxygen saturation image at intervals of a preset time period, and obtaining a second position of the target blood flow blood oxygen signal area; calculating a first position offset value of the target blood flow blood oxygen signal area according to the first position and the second position; and adjusting the first target detection area based on the first position offset value. With the above implementation scheme, the target blood flow oxygen signal is first located at a first position in the first target detection area. Thereafter, at preset time intervals, the target blood flow oxygen signal area is detected from the oxygen saturation image, and the second position of the target blood flow oxygen signal is determined. The first position offset value of the second position from the first position is determined, and the first target detection area is adjusted based on the first position offset value. The new target detection area can be automatically located according to the deviation of the target blood flow oxygen signal, thereby reducing the time for doctors to pay attention to the oxygen saturation image and reducing the impact of doctors observing the target blood flow signal; at the same time, the intelligence of signal detection is improved. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] Figure 1 A structural diagram of a photoacoustic imaging system provided by an embodiment of the present invention;

[0030] Figure 2 A process for processing a photoacoustic signal provided by an embodiment of the present invention Figure 1 ;

[0031] Figure 3 A schematic diagram of an exemplary target observation area in a blood oxygen saturation image provided by an embodiment of the present invention;

[0032] Figure 4A schematic diagram of an exemplary first target detection area 1 before displacement provided by an embodiment of the present invention;

[0033] Figure 5 A schematic diagram of an exemplary second target detection area 2 after displacement provided by an embodiment of the present invention;

[0034] Figure 6 A block diagram of an exemplary method for automatically detecting and tracking target blood flow oxygen signal and abnormal prompts provided by an embodiment of the present invention;

[0035] Figure 7 A process for processing a photoacoustic signal provided by an embodiment of the present invention Figure 2 ;

[0036] Figure 8 A schematic structural diagram of a photoacoustic imaging system 1 provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0037] In order to enable a more detailed understanding of the features and technical contents of the embodiments of the present invention, the implementation of the embodiments of the present invention is described in detail below in conjunction with the accompanying drawings. The attached drawings are for reference only and are not intended to limit the embodiments of the present invention.

[0038] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art to which the present invention belongs. The terms used herein are only for the purpose of describing the embodiments of the present invention and are not intended to limit the present invention.

[0039] In the following description, reference is made to "some embodiments", which describe a subset of all possible embodiments, but it is understood that "some embodiments" may be the same subset or different subsets of all possible embodiments, and may be combined with each other without conflict. It is also necessary to point out that the terms "first\second\third" involved in the embodiments of the present invention are only used to distinguish similar objects and do not represent a specific order for the objects. It is understood that "first\second\third" may be interchanged in a specific order or sequence where permitted, so that the embodiments of the present invention described herein can be implemented in an order other than that illustrated or described herein.

[0040] like Figure 1 As shown, a photoacoustic imaging system provided by an embodiment of the present invention includes an ultrasonic probe, optical fibers arranged at both ends of the ultrasonic probe, a laser transmitting module, an ultrasonic transmitting module, a scanning control device, an echo receiving module, a processor and a display module, wherein the processor includes an ultrasonic signal processing module, a photoacoustic processing module, an ultrasonic processing module, a laser energy reading module and a fusion processing module. The functions of each module include:

[0041] Scanning control device: controls the laser emitting module to generate laser and perform laser scanning or the ultrasonic emitting module to generate ultrasonic wave and perform ultrasonic scanning, and notifies the ultrasonic signal processing module that the current scanning mode is laser scanning or ultrasonic scanning.

[0042] Laser emission module: generates laser with variable wavelength, which is emitted to the target object through optical fibers located on both sides of the probe.

[0043] Ultrasonic transmitting module: generates ultrasonic transmitting waveform and transmits it to the target object through the ultrasonic probe.

[0044] Echo receiving module: receives the echo signal emitted by the target object.

[0045] Ultrasonic signal processing module: performs general signal processing on echo signals, including filtering, beamforming, signal processing, etc.

[0046] Photoacoustic processing module: performs photoacoustic image processing on the photoacoustic data output by the ultrasonic signal processing module, including image smoothing, enhancement, compensation and correction of photoacoustic images, etc.

[0047] Ultrasonic processing module: performs ultrasonic image processing on the ultrasonic data output by the ultrasonic signal processing module, including image smoothing, enhancement, removal of noise in the image, and improvement of image contrast.

[0048] The laser energy reading module is used to read the energy value of each laser emission.

[0049] Fusion processing module: Fusion processing of the data output by the ultrasonic processing module and the photoacoustic processing module.

[0050] Display module: displays the final image.

[0051] The following embodiments are all described based on the above-mentioned photoacoustic imaging system.

[0052] The embodiment of the present invention provides a method for processing photoacoustic signals, such as Figure 2 As shown, the method may include:

[0053] S101. Acquire a blood oxygen saturation image of an imaging target.

[0054] In an embodiment of the present invention, a blood oxygen saturation image of an imaging target can be acquired by two methods. One method is to transmit a laser signal having a first wavelength to the imaging target, receive a first ultrasonic echo signal returned by the imaging target, generate first photoacoustic image data based on the first ultrasonic echo signal, transmit a laser signal having a second wavelength to the imaging target, receive a second ultrasonic echo signal returned by the imaging target, generate second photoacoustic image data based on the second ultrasonic echo signal, and perform blood oxygen reconstruction on the first photoacoustic image data and the second photoacoustic image data to obtain a blood oxygen saturation image; or, the second method is to acquire a blood oxygen saturation image of the imaging target from a storage medium; or, the data of the photoacoustic image can be acquired from a storage medium, and the data of the photoacoustic image can be processed to acquire a corresponding blood oxygen saturation image. The specific method for acquiring the blood oxygen saturation image of the imaging target can be selected according to actual conditions, and the embodiment of the present invention does not make any specific limitation.

[0055] S102. Determine a first target detection area on the blood oxygen saturation image, wherein the first target detection area includes a target blood flow oxygen signal area, and the target blood flow oxygen signal area has a first position.

[0056] In an optional embodiment, a first target observation area containing a blood flow oxygen signal can be automatically selected from a blood oxygen saturation image obtained by photoacoustic imaging. Specifically, the first target detection area is determined on the blood oxygen saturation image, including: in the blood oxygen saturation image, dividing the blood flow oxygen signal into at least one group of blood flow oxygen signal sets; determining a region boundary of at least one blood flow oxygen signal set in the blood oxygen saturation image, and determining an initial detection area for the image area framed by the region boundary; collecting signal parameters of at least one group of blood flow oxygen signals in the initial detection area within another preset time period; searching for a target blood flow oxygen signal area from at least one group of blood flow oxygen signals according to the signal parameters; determining a first position of the target blood flow oxygen signal area; and adjusting the initial detection area using the first position to obtain a first target detection area.

[0057] Optionally, the first position may be a center position, a centroid position, or any position, and the specific position may be selected according to actual conditions, and the embodiment of the present invention does not make any specific limitation.

[0058] It should be noted that the normal human arterial blood oxygen saturation is 93%-98%, and the venous blood oxygen saturation is 70%-75%. The actual measured value will fluctuate within the above range. Based on this feature, the blood oxygen signal generated by some other human tissues can be distinguished from the blood flow blood oxygen signal. The noise interference can also be filtered out through threshold conditions such as signal strength and signal continuity.

[0059] On the basis of the blood flow oxygen saturation conditions described above, the depth range most suitable for observation in photoacoustic clinical applications can be preferentially traversed and searched, and the blood flow oxygen signals that meet the above blood flow oxygen saturation conditions can be divided into at least one group of blood flow oxygen signal sets by clustering or segmentation methods. After that, image growth or other means can be used to determine the regional boundaries of at least one group of blood flow oxygen signal sets in the blood oxygen saturation image, and the initial detection area containing the blood flow oxygen signal is determined according to the image area framed by the regional boundaries, and the areas are sorted and presented according to fixed standards for user selection. The fixed standards can come from prior indications and accumulated experience in photoacoustic clinical applications.

[0060] Exemplary, reference Figure 3 , within the specified search interval, the initial detection areas corresponding to the four positions 1, 2, 3 and 4 all contain blood flow oxygen signals.

[0061] It should be noted that the search operation for blood flow oxygen signal can also be performed in other depth ranges according to the above operation sequence. If no blood flow oxygen signal that meets the blood flow oxygen saturation conditions can be found in the entire depth range, the user will be prompted that no blood flow oxygen signal is detected.

[0062] In the embodiment of the present invention, at least one set of signal parameters of blood flow oxygen signal in the initial detection area within another preset time period is collected; wherein the signal parameters include at least one of signal strength and signal continuity. The specific parameters can be selected according to the actual situation, and the embodiment of the present invention does not make specific limitations.

[0063] Specifically, according to at least one signal parameter of signal strength and signal continuity, a region where a stable blood flow oxygen signal is located is searched from at least one group of blood flow oxygen signals, and the region where the stable blood flow oxygen signal is located is the target blood flow oxygen signal region.

[0064] In the embodiment of the present invention, the center position of the target blood flow oxygen signal area can be determined based on the centroid weighted method, wherein the centroid weighted method is to use all anchor nodes within the communication range of the unknown node for positioning, obtain the horizontal and vertical coordinates of all anchor nodes, and then, consider the influence of the distance between the anchor node and the unknown node on the unknown node, assign a weight value to each anchor node; finally, based on the weight value and the horizontal and vertical coordinate values ​​of each anchor node, calculate the estimated position of the unknown node, which is the first position of the present invention.

[0065] In the embodiment of the present invention, the initial observation area is adjusted based on the first position of the target blood flow oxygen signal to obtain a first target detection area including the target blood flow oxygen signal.

[0066] For example, see Figure 4, the first position of the target blood flow oxygen signal is located using the centroid weighted method, and the white box is the first target detection area 1.

[0067] In another optional embodiment, a target observation area containing a blood flow oxygen signal can also be manually selected from a blood oxygen saturation image obtained by photoacoustic imaging. Specifically, a first target detection area is determined on the blood oxygen saturation image, including: displaying the blood oxygen saturation image, obtaining a first user input through a user interaction device, and determining the first target detection area based on the user input, wherein the user input includes a first area selected by the user on the blood oxygen saturation image through the user interaction device.

[0068] In the embodiment of the present invention, the blood oxygen saturation image can be displayed on the display module of the photoacoustic imaging system.

[0069] S103. Detect a target blood flow oxygen signal region from the blood oxygen saturation image at every preset time period, and obtain a second position of the target blood flow oxygen signal region.

[0070] In an embodiment of the present invention, the second position of the target blood flow and blood oxygen area can be obtained based on any motion tracking method among the preset area matching method, the preset optical flow method and the preset feature matching method. The specific motion tracking method can be selected according to actual conditions, and the embodiment of the present invention does not make any specific limitation.

[0071] It should be noted that the preset region matching method is one of the most classic methods of speckle tracking in the field of motion estimation. Its core concept is to use the nth frame image as a reference and find the most likely target region position in the n+1th frame image. The iterative process needs to follow a fixed reference standard. The standard value can be positive, such as the normalized cross-correlation coefficient; it can also be negative, such as the minimum square error. The preset region matching method expresses the spatial motion of the target using a mathematical model containing a control equation, such as a Kalman filter.

[0072] It should be noted that optical flow is the instantaneous speed of the pixel motion mapped on the observation imaging plane when the object moves in space. The basic principle of the preset optical flow method is to give each pixel in the image a velocity vector, thereby constructing an image motion field that corresponds to the points on the three-dimensional object one by one. According to the velocity vector of each pixel, the dynamic analysis of the image is realized.

[0073] It should be noted that the preset feature matching method achieves feature matching by continuously reselecting the target point and calculating the motion vector, thereby tracking the second position of the target blood flow oxygen signal after a preset period of time.

[0074] S104. Calculate a first position offset value of the target blood flow oxygen signal region according to the first position and the second position.

[0075] In the embodiment of the present invention, the displacement change amount between the first position and the second position is determined, and the displacement change amount is determined as the first position offset value of the target blood flow oxygen signal.

[0076] Furthermore, after determining the first position offset value of the target blood flow oxygen signal area according to the first position and the second position, the risk degree and the abnormality type can also be determined according to the first position offset value; and an abnormality prompt corresponding to the risk degree and the abnormality type is output.

[0077] In an embodiment of the present invention, a signal displacement velocity value, a signal strength change value, a gradient change value and a signal morphology change of a target blood flow blood oxygen signal in a target blood flow blood oxygen signal area within a preset time period are obtained; an overall offset value of the target blood flow blood oxygen signal area is determined according to the first position offset value, the signal displacement velocity value, the signal strength change value, the gradient change value and the signal morphology change; if the overall offset value is within a preset offset range, the risk level is determined to be low risk or medium risk; if the overall offset value is not within the preset offset range, the risk level is determined to be high risk; and the abnormality type is determined according to the comparison result of the signal position offset value, the signal displacement velocity value, the signal strength change value, the gradient change value and the signal morphology change with the corresponding preset value range.

[0078] In the embodiment of the present invention, the long-distance displacement of the signal center position of the target blood flow oxygen signal, that is, the signal center position offset value is greater than a preset position offset threshold; and the signal center position of the target blood flow oxygen signal is rapidly displaced, that is, the signal displacement speed value is greater than a preset displacement speed threshold; and the signal strength jump, that is, the signal strength change is greater than a preset signal strength change threshold; and the signal morphology similarity is weakened, that is, the signal morphology change is greater than a preset morphology threshold; and the signal stability is weakened, that is, the gradient change value is greater than a preset gradient change threshold; the above multiple manifestations can reflect the change of the current measurement signal state, which can be uniformly described as the degree of deviation from the target blood flow oxygen signal. Therefore, the offset value can be determined according to the center position offset value, signal displacement speed value, signal strength change value, gradient change value and signal morphology change; and the risk level is determined to be low risk, medium risk or high risk according to the comparison between the offset value and the corresponding preset offset range.

[0079] Optionally, the type of the abnormal prompt includes one of the following three types: a pop-up window warning prompt, a sound prompt, and a display identifier prompt. The specific type can be selected according to actual conditions, and the embodiment of the present application does not make specific limitations.

[0080] In an optional embodiment, if the offset value is within a preset offset range, the risk level is determined to be low risk or medium risk, and a low risk or medium risk prompt can be given by displaying an identifier prompt; when the offset value is not within the preset offset range, the risk level is determined to be high risk, and at this time, a high risk prompt can be given by a pop-up warning or a sound alarm.

[0081] It should be noted that, while providing a high-risk prompt, it may also be possible to suggest re-adjusting the position for collecting blood flow and blood oxygen signals; and to determine the abnormality type based on the comparison results of the signal position offset value, signal displacement speed value, signal intensity change value, gradient change value and signal morphology change with the corresponding preset value range, and provide a prompt for the abnormality type, etc. The specific selection can be made according to the actual situation, and the embodiments of the present invention do not make specific limitations.

[0082] It should be noted that factors such as imaging section changes, probe position offset, and output energy intensity changes will affect the accuracy and effectiveness of blood oxygen saturation results and interfere with doctors' observations and judgments. Among them, imaging slice changes correspond to two manifestations: signal intensity jump and weakened signal morphology similarity; probe position offset corresponds to three manifestations: long-distance displacement of the signal center position of the target blood flow oxygen signal, rapid displacement of the signal center position of the target blood flow oxygen signal, and weakened signal stability; output energy intensity changes correspond to one manifestation: weakened signal morphology similarity. Based on this, the abnormal types such as imaging slice changes, probe position offsets, and output energy intensity changes can be determined by comparing the signal position offset value, signal displacement speed value, signal intensity change value, gradient change value, and signal morphology change with the corresponding preset value range.

[0083] S105: Adjust the first target detection area based on the first position offset value.

[0084] In an embodiment of the present invention, the first position offset value is compared with the first preset offset threshold value. If the first position offset value is less than the first preset offset threshold value and the target blood flow oxygen signal area is located within the first target detection area, indicating that the target blood flow oxygen signal does not exceed the first target detection area, the first target detection area is not adjusted; if the first position offset value is greater than or equal to the first preset offset threshold value, or the target blood flow oxygen signal area is not within the first target detection area, the first target detection area is adjusted according to the first position offset value and the target blood flow oxygen signal area to ensure that the target blood flow oxygen signal is always within the detection area, unless the target blood flow oxygen signal is in the blood oxygen saturation image due to imaging section or other abnormal types.

[0085] See also Figure 5 , is the second position of the target blood flow oxygen signal based on the centroid weighted positioning, Figure 4The second target detection area 2 is obtained by adaptively adjusting the first target detection area 1 in .

[0086] It should be noted that after adjusting the first target detection area, the blood oxygen saturation of the imaging target in the first target detection area is quantitatively analyzed.

[0087] Specifically, quantitatively analyzing the blood oxygen saturation of the imaging target in the first target detection area includes: calculating the minimum blood oxygen saturation value, the maximum blood oxygen saturation value and / or the average blood oxygen saturation value of the imaging target in the first target detection area. Specifically, the type of blood oxygen saturation to be analyzed can be selected according to actual conditions, and the embodiment of the present invention does not make specific limitations.

[0088] It can be understood that the target blood flow oxygen signal is first located at the first position of the first target detection area, and then, at preset time intervals, the target blood flow oxygen signal area is detected from the blood oxygen saturation image, and the second position of the target blood flow oxygen signal is determined, the first position offset value of the second position from the first position is determined, and the first target detection area is adjusted based on the first position offset value. The new target detection area can be automatically located according to the deviation of the target blood flow oxygen signal, thereby reducing the time that the doctor pays attention to the blood oxygen saturation image to reduce the impact of the doctor's observation of the target blood flow signal; at the same time, the intelligence of the signal detection is improved.

[0089] Based on the above embodiment, a method for automatically detecting and tracking target blood flow oxygen signal and abnormal prompt is proposed. Figure 6 First, manually or automatically select the first target observation area in the blood oxygen saturation image of the photoacoustic imaging interface, and ensure that there is a stable target blood flow oxygen signal in the first target observation area, and use the centroid weighted positioning method to determine the first position of the target blood flow oxygen signal; during the monitoring process, periodically detect the target blood flow oxygen signal area from the blood oxygen saturation image, and obtain the second position of the target blood flow oxygen signal area, and judge whether the first target detection area deviates from the original position according to the second position of the target blood flow oxygen signal and the first position offset value of the first position; if so, adjust the ROI (first target detection area) according to the first position offset value, and judge the risk level according to the first position offset value, and issue an abnormal prompt when the risk is high, if not, return to continuously detect the target blood flow oxygen signal area, and obtain the second position of the target blood flow oxygen signal area.

[0090] Based on the above embodiments, the present invention further provides a method for processing photoacoustic signals, such as Figure 7 As shown, the method includes:

[0091] S201 , obtaining a photoacoustic image of an imaging target.

[0092] In an embodiment of the present invention, the photoacoustic image of the imaging target can be obtained by emitting a laser signal to the imaging target, receiving an ultrasonic signal returned by the imaging target, and generating the photoacoustic image of the imaging target based on the ultrasonic signal; it can also be obtained from a storage medium. The specific method of obtaining the photoacoustic image of the imaging target can be selected according to actual conditions, and the embodiment of the present invention does not make any specific limitation.

[0093] S202: Determine a second target detection area on the photoacoustic image, wherein the second target detection area includes a target photoacoustic signal, and the target photoacoustic signal has a third position.

[0094] It should be noted that, for S202 , reference may be made to the description of S102 , which will not be repeated here.

[0095] S203 . Detect a target photoacoustic signal region from the photoacoustic image at every preset time period, and obtain a fourth position of the target photoacoustic signal region.

[0096] It should be noted that, for S203 , reference may be made to the description of S103 , which will not be repeated here.

[0097] S204 : Calculate a second position offset value of the target photoacoustic signal region according to the third position and the fourth position.

[0098] It should be noted that, for S204 , reference may be made to the description of S104 , which will not be repeated here.

[0099] Furthermore, after determining the position offset value of the target blood flow oxygen signal area according to the third position and the fourth position, the risk degree and the abnormality type can also be determined according to the second position offset value; and an abnormality prompt corresponding to the risk degree and the abnormality type is output.

[0100] Optionally, the prompt type of the abnormal prompt includes one of the following three types: a pop-up window warning prompt, a sound prompt, and a display identifier prompt. The specific type can be selected according to actual conditions, and the embodiment of the present invention does not make a specific limitation.

[0101] S205: Adjust the second target detection area based on the second position offset value.

[0102] In an embodiment of the present invention, if the second position offset value is less than a preset offset threshold and the target photoacoustic signal area is within the second target detection area, the second target detection area is not adjusted; if the second position offset value is greater than or equal to the preset offset threshold, or the target photoacoustic signal area is not within the second target detection area, the second target detection area is adjusted according to the second position offset value and the target blood flow oxygen signal area.

[0103] It should be noted that, for S205, reference may be made to the description of S105, which will not be repeated here. In the above steps S201-S205, based on the automatic adjustment of the target photoacoustic signal area, not only can the real-time tracking of the blood vessel area be achieved, but also the target photoacoustic signal area on the photoacoustic image can be matched with the target blood flow oxygen signal area on the blood oxygen saturation image, and based on the automatic adjustment of the target photoacoustic signal area, the automatic adjustment of the signal area of ​​the target blood flow oxygen signal can be achieved accordingly.

[0104] The embodiment of the present invention further provides a photoacoustic imaging system 1, such as Figure 8 As shown, the photoacoustic imaging system 1 comprises:

[0105] A laser emitting device 10, wherein the laser emitting device generates a laser to illuminate a target object;

[0106] Optical fiber 11;

[0107] Ultrasonic transducer 12;

[0108] A scanning control device 13, the scanning control device controls the laser emitting device to generate the laser;

[0109] A receiving circuit 14, wherein the receiving circuit receives the photoacoustic signal generated by the target object through the ultrasonic transducer to obtain photoacoustic data;

[0110] A processor 15, which processes the photoacoustic data to obtain a blood oxygen saturation image;

[0111] A display 16, wherein the display displays the blood oxygen saturation image;

[0112] The processor 15 further performs the following steps:

[0113] Acquire a blood oxygen saturation image of an imaging target; determine a target detection area on the blood oxygen saturation image, wherein the target detection area includes a target blood flow oxygen signal area, and the target blood flow oxygen signal area has a first position; detect the target blood flow oxygen signal area from the blood oxygen saturation image at intervals of a preset time period, and obtain a second position of the target blood flow oxygen signal area; calculate a position offset value of the target blood flow oxygen signal area according to the first position and the second position; and adjust the target detection area based on the position offset value.

[0114] Optionally, the processor 15 is also used to transmit a laser signal with a first wavelength to the imaging target, receive a first ultrasonic echo signal returned by the imaging target, generate first photoacoustic image data based on the first ultrasonic echo signal, transmit a laser signal with a second wavelength to the imaging target, receive a second ultrasonic echo signal returned by the imaging target, and generate second photoacoustic image data based on the second ultrasonic echo signal; perform blood oxygen reconstruction on the first photoacoustic image data and the second photoacoustic image data to obtain a blood oxygen saturation image; or, obtain the blood oxygen saturation image of the imaging target from a storage medium.

[0115] Optionally, the processor 15 is also used to not adjust the first target detection area if the first position offset value is less than a first preset offset threshold and the target blood flow oxygen signal area is within the first target detection area; if the first position offset value is greater than or equal to the first preset offset threshold, or the target blood flow oxygen signal area is not within the first target detection area, adjust the first target detection area according to the first position offset value and the target blood flow oxygen signal area.

[0116] Optionally, the processor 15 is further used to quantitatively analyze the blood oxygen saturation of the imaging target in the first target detection area after adjusting the first target detection area.

[0117] Optionally, the processor 15 is further used to calculate the minimum blood oxygen saturation value, the maximum blood oxygen saturation value and / or the average blood oxygen saturation value of the imaging target in the first target detection area.

[0118] Optionally, the processor 15 is also used to display the blood oxygen saturation image, obtain a first user input through a user interaction device, and determine a first target detection area based on the user input, wherein the user input includes a first area selected by the user on the blood oxygen saturation image through the user interaction device.

[0119] Optionally, the processor 15 is further used to divide the blood flow oxygen signal into at least one group of blood flow oxygen signal sets in the blood oxygen saturation image; determine the regional boundary of the at least one blood flow oxygen signal set in the blood oxygen saturation image, and determine the initial detection area with the image area framed by the regional boundary; collect signal parameters of at least one group of blood flow oxygen signals in the initial detection area within another preset time period; search for a target blood flow oxygen signal area from the at least one group of blood flow oxygen signals according to the signal parameters; determine a first position of the target blood flow oxygen signal area; and adjust the initial detection area using the first position to obtain the first target detection area.

[0120] Optionally, the processor 15 is further used to determine the risk level and the abnormality type according to the first position offset value; and output an abnormality prompt corresponding to the risk level and the abnormality type.

[0121] Optionally, the processor 15 is also used to obtain the signal displacement speed value, signal strength change value, gradient change value and signal morphology change of the target blood flow blood oxygen signal in the target blood flow blood oxygen signal area within a preset time period; determine the overall offset value of the target blood flow blood oxygen signal area according to the first position offset value, the signal displacement speed value, the signal strength change value, the gradient change value and the signal morphology change; if the overall offset value is within the preset offset range, determine the risk level as low risk or medium risk; if the overall value is not within the preset offset range, determine the risk level as high risk; and determine the abnormality type according to the comparison result of the signal position offset value, the signal displacement speed value, the signal strength change value, the gradient change value and the signal morphology change with the corresponding preset value range.

[0122] Optionally, the prompt type of the abnormal prompt includes one of the following three types: a pop-up warning, a sound prompt and / or a display identifier prompt.

[0123] Optionally, the processor 15 is further used to obtain the second position of the target blood flow and blood oxygen area based on any one of a preset area matching method, a preset optical flow method and a preset feature matching method.

[0124] An embodiment of the present invention provides a storage medium having a computer program stored thereon. The computer-readable storage medium stores one or more programs. The one or more programs can be executed by one or more processors and are applied to a photoacoustic imaging system. The computer program implements the method for processing photoacoustic signals as described above.

[0125] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element.

[0126] Through the description of the above implementation methods, those skilled in the art can clearly understand that the above-mentioned embodiment methods can be implemented by means of software plus a necessary general hardware platform, and of course by hardware, but in many cases the former is a better implementation method. Based on such an understanding, the technical solution of the present disclosure, or the part that contributes to the relevant technology, can be embodied in the form of a software product, which is stored in a storage medium (such as ROM / RAM, a disk, or an optical disk), and includes a number of instructions for enabling an image display device (which can be a mobile phone, a computer, a server, an air conditioner, or a network device, etc.) to execute the methods described in each embodiment of the present disclosure.

[0127] The above description is only a preferred embodiment of the present invention and is not intended to limit the protection scope of the present invention.

Claims

1. A method for processing photoacoustic signals, characterized in that: The method comprises: Acquire a blood oxygen saturation image of an imaging target; Determine a first target detection area on the blood oxygen saturation image, wherein the first target detection area includes a target blood flow oxygen signal area, and the target blood flow oxygen signal area has a first position; At every preset time period, detecting the target blood flow oxygen signal area from the blood oxygen saturation image, and obtaining a second position of the target blood flow oxygen signal area; Calculate a first position offset value of the target blood flow oxygen signal area according to the first position and the second position; The first object detection area is adjusted based on the first position offset value.

2. The method according to claim 1, characterized in that The step of acquiring a blood oxygen saturation image of an imaging target includes: transmitting a laser signal having a first wavelength to the imaging target, receiving a first ultrasonic echo signal returned by the imaging target, and generating first photoacoustic image data based on the first ultrasonic echo signal, transmitting a laser signal having a second wavelength to the imaging target, receiving a second ultrasonic echo signal returned by the imaging target, and generating second photoacoustic image data based on the second ultrasonic echo signal, Performing blood oxygen reconstruction on the first photoacoustic image data and the second photoacoustic image data to obtain a blood oxygen saturation image; or, Acquire a blood oxygen saturation image of the imaging target from a storage medium.

3. The method according to claim 1, characterized in that: The adjusting the first target detection area based on the first position offset value includes: If the first position offset value is less than a first preset offset threshold, and the target blood flow oxygen signal area is located within the first target detection area, the first target detection area is not adjusted; If the first position offset value is greater than or equal to the first preset offset threshold, or the target blood flow oxygen signal area is not within the first target detection area, the first target detection area is adjusted according to the first position offset value and the target blood flow oxygen signal area.

4. The method according to claim 1 or 3, characterized in that: The method further includes, after adjusting the first target detection area, quantitatively analyzing the blood oxygen saturation of the imaging target within the first target detection area.

5. The method according to claim 4, characterized in that The quantitative analysis of the blood oxygen saturation of the imaging target in the first target detection area includes: calculating the minimum blood oxygen saturation value, the maximum blood oxygen saturation value and / or the average blood oxygen saturation value of the imaging target in the first target detection area.

6. The method according to claim 1, characterized in that The determining a first target detection area on the blood oxygen saturation image includes: The blood oxygen saturation image is displayed, a first user input is obtained through a user interaction device, and a first target detection area is determined based on the user input, wherein the user input includes a first area selected by the user on the blood oxygen saturation image through the user interaction device.

7. The method according to claim 1, characterized in that The determining a first target detection area on the blood oxygen saturation image includes: In the blood oxygen saturation image, dividing the blood flow oxygen signal into at least one group of blood flow oxygen signal sets; Determine a region boundary of the at least one blood flow oxygen signal set in the blood oxygen saturation image, and determine an initial detection region by an image region framed by the region boundary; Collecting at least one set of signal parameters of blood flow oxygen signal in the initial detection area within another preset time period; According to the signal parameters, searching for a target blood flow oxygen signal area from the at least one group of blood flow oxygen signals; Determining a first position of the target blood flow oxygen signal area; The initial detection area is adjusted using the first position to obtain the first target detection area.

8. The method according to claim 1, characterized in that After determining the first position offset value of the target blood flow oxygen signal area according to the first position and the second position, the method further includes: Determining the risk level and the abnormality type according to the first position offset value; Output the risk level and the abnormality prompt corresponding to the abnormality type.

9. The method according to claim 8, characterized in that The method further comprises: Acquire a signal displacement velocity value, a signal intensity change value, a gradient change value, and a signal morphology change of a target blood flow oxygen signal in the target blood flow oxygen signal area within a preset time period; Determine an overall offset value of the target blood flow oxygen signal area according to the first position offset value, the signal displacement speed value, the signal intensity change value, the gradient change value and the signal morphology change; If the overall deviation value is within a preset deviation range, the risk level is determined to be low risk or medium risk; If the overall value is not within the preset deviation range, the risk level is determined to be high risk; and the abnormality type is determined based on the comparison results of the signal position deviation value, the signal displacement speed value, the signal intensity change value, the gradient change value and the signal morphology change with the corresponding preset value range.

10. The method according to claim 8 or 9, characterized in that: The prompt type of the abnormal prompt includes one of the following three types: a pop-up warning prompt, a sound prompt, and a display identifier prompt.

11. The method according to claim 1, characterized in that: The obtaining of the second position of the target blood flow oxygen signal area includes: obtaining the second position of the target blood flow oxygen signal area based on any one of a preset area matching method, a preset optical flow method and a preset feature matching method.

12. A method for processing photoacoustic signals, characterized in that: The method comprises: acquiring a photoacoustic image of an imaging target; determining a second target detection region on the photoacoustic image, wherein the second target detection region includes a target photoacoustic signal region, and the target photoacoustic signal region has a third position; At every preset time period, detecting the target photoacoustic signal region from the photoacoustic image, and obtaining a fourth position of the target photoacoustic signal region; Calculating a second position offset value of the target photoacoustic signal area according to the third position and the fourth position; The second object detection area is adjusted based on the second position offset value.

13. The method according to claim 12, characterized in that The adjusting the second target detection area according to the second position offset value includes: If the second position offset value is less than a preset offset threshold and the target photoacoustic signal area is located within the second target detection area, the second target detection area is not adjusted; If the second position offset value is greater than or equal to the preset offset threshold, or the target photoacoustic signal area is not within the second target detection area, the second target detection area is adjusted according to the second position offset value and the target signal area.

14. The method according to claim 12, characterized in that After determining the position offset value of the target signal area according to the third position and the fourth position, the method further includes: determining a risk level and an abnormality type according to the second position offset value; Output the risk level and the abnormality prompt corresponding to the abnormality type.

15. The method according to claim 14, characterized in that The prompt type of the abnormal prompt includes one of the following three types: a pop-up warning prompt, a sound prompt, and a display identifier prompt.

16. A photoacoustic imaging system, characterized in that: The photoacoustic imaging system comprises: A laser emitting device, wherein the laser generating device generates laser light for irradiating an imaging target; optical fiber; Ultrasonic transducer; A scanning control device, wherein the scanning control device controls the laser emitting device to generate the laser; A receiving circuit, wherein the receiving circuit receives the photoacoustic signal generated by the imaging target through the ultrasonic transducer to obtain photoacoustic data; a processor, wherein the processor processes the photoacoustic data to obtain a blood oxygen saturation image; a display, wherein the display displays the blood oxygen saturation image; The processor further performs the following steps: Acquire a blood oxygen saturation image of an imaging target; determine a target detection area on the blood oxygen saturation image, wherein the target detection area includes a target blood flow oxygen signal area, and the target blood flow oxygen signal area has a first position; detect the target blood flow oxygen signal area from the blood oxygen saturation image at intervals of a preset time period, and obtain a second position of the target blood flow oxygen signal area; calculate a position offset value of the target blood flow oxygen signal area according to the first position and the second position; and adjust the target detection area based on the position offset value.

17. A storage medium having a computer program stored thereon, applied to a photoacoustic imaging system, wherein the computer program, when executed by a processor, implements the method according to any one of claims 1 to 15.

Images