A method, device, storage medium and electronic device for photoacoustic imaging demixing

Abstract

This specification discloses a method, device, storage medium and electronic device for photoacoustic imaging unmixing. When using original lasers of multiple wavelengths to irradiate the sample to be unmixed, for each wavelength of original laser, a first light energy distribution of the sample to be unmixed is determined when the laser of the wavelength is emitted from a plurality of light outlets to irradiate the sample to be unmixed. The first light energy distribution characterizes the light energy distribution of the points inside the sample to be unmixed. Then, based on the first light energy distribution of the sample to be unmixed at the wavelength and the first photoacoustic image, the second photoacoustic image at the wavelength can be determined. Thus, the sample to be unmixed can be unmixed according to the second photoacoustic image corresponding to the original laser of each wavelength. This method takes into account the different light energy distributions inside the sample to be unmixed under the irradiation of lasers of different wavelengths, thereby improving the accuracy of photoacoustic imaging unmixing based on multiple wavelengths.

Description

Technical Field

[0001] The present application relates to the field of biomedical imaging, and in particular to a method, device, storage medium and electronic device for photoacoustic imaging unmixing. Background Art

[0002] With the development of science and technology, photoacoustic imaging technology has received widespread attention. Photoacoustic imaging technology is a non-invasive biomedical imaging technology. Its principle is: use pulsed laser to irradiate the sample, the sample absorbs photons and releases energy, thereby generating photoacoustic signals, and reconstructing the photoacoustic signals to obtain the light absorption image of the sample. Among them, the components in the sample have different absorption coefficients of light, so the proportion of each component in the sample can be calculated through the light absorption image of the sample. The above process of solving the proportion of each component in the sample based on the light absorption image of the sample is called unmixing.

[0003] In practical applications, multi-wavelength pulsed lasers are generally used to irradiate samples. However, the absorption coefficient and scattering coefficient of samples for pulsed lasers of different wavelengths are different. As a result, the different light energy distributions of lasers of different wavelengths in the sample reduce the accuracy of multi-wavelength photoacoustic imaging unmixing. How to use multi-wavelength lasers to irradiate samples for accurate photoacoustic imaging unmixing is an important issue.

[0004] Based on this, the present application specification provides a method for photoacoustic imaging unmixing. Summary of the invention

[0005] The present specification provides a method, device, storage medium and electronic device for photoacoustic imaging demixing to at least partially solve the above-mentioned problems existing in the prior art.

[0006] This manual adopts the following technical solutions:

[0007] This specification provides a method for photoacoustic imaging unmixing, the method comprising:

[0008] Obtaining the positional relationship between the sample to be unmixed and a plurality of light outlets;

[0009] For each of the original lasers of several wavelengths, according to the positional relationship between the sample to be unmixed and the several light outlets, determine the first light energy distribution of the sample to be unmixed when the several light outlets emit the laser of the wavelength to irradiate the sample to be unmixed, wherein the first light energy distribution represents the light energy distribution of points inside the sample to be unmixed; and, according to the photoacoustic signal generated when the sample to be unmixed is irradiated by the original laser of the wavelength, reconstruct the first photoacoustic image corresponding to the sample to be unmixed when it is irradiated by the original laser of the wavelength;

[0010] Determine, according to the first light energy distribution and the first photoacoustic image, a second photoacoustic image when the sample to be unmixed is irradiated by the original laser of the wavelength;

[0011] The sample to be unmixed is unmixed according to the determined second photoacoustic image when the sample to be unmixed is irradiated by the original laser of each wavelength.

[0012] Optionally, determining the first light energy distribution of the sample to be unmixed when the plurality of light outlets emit lasers of the wavelength to irradiate the sample to be unmixed specifically includes:

[0013] Determine the light energy of the laser of the wavelength emitted from each light outlet; and, according to the positional relationship between the sample to be unmixed and each light outlet, determine the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet;

[0014] Determine, according to the light energy of the laser of the wavelength emitted from each light outlet and the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet, the second light energy distribution of the sample to be unmixed corresponding to each light outlet when each light outlet emits the laser of the wavelength to irradiate the sample to be unmixed;

[0015] Based on the determined second light energy distributions, the first light energy distribution of the sample to be unmixed is determined; wherein the second light energy distribution represents the light energy distribution of a point inside the sample to be unmixed when each light outlet emits a laser of this wavelength to irradiate the sample to be unmixed.

[0016] Optionally, determining the light energy of the laser light of the wavelength emitted from each light outlet specifically includes:

[0017] Obtaining the optical energy of the original laser of the wavelength;

[0018] For each of the several splitting ratios, under the splitting ratio, according to the light energy of the original laser of the wavelength, the light energy of the laser of the wavelength emitted from each light outlet is determined, wherein the splitting ratio is used to divide the original laser according to a preset energy ratio, and the number of lasers obtained by division is consistent with the number of light outlets;

[0019] Determining the first light energy distribution of the sample to be unmixed based on the determined second light energy distributions specifically includes:

[0020] Based on the determined second light energy distributions, the first light energy distribution of the sample to be unmixed is determined under the light splitting ratio.

[0021] Optionally, determining the second light energy distribution of the sample to be unmixed when the light outlet emits the laser of the wavelength to irradiate the sample to be unmixed specifically includes:

[0022] Determine the light attenuation distribution of a point inside the sample to be unmixed relative to a point on the surface of the sample to be unmixed corresponding to the light outlet;

[0023] According to the light attenuation distribution, the light energy of the laser of the wavelength emitted from the light outlet, and the light intensity distribution on the surface of the sample to be unmixed corresponding to the light outlet, the second light energy distribution of the sample to be unmixed is determined when the laser of the wavelength emitted from the light outlet irradiates the sample to be unmixed.

[0024] Optionally, reconstructing a first photoacoustic image corresponding to the sample to be unmixed when irradiated by the original laser of the wavelength specifically includes:

[0025] For each of the several splitting ratios, under the splitting ratio, a first photoacoustic image corresponding to the sample to be unmixed when irradiated by the original laser of the wavelength is reconstructed.

[0026] Optionally, determining a second photoacoustic image of the sample to be unmixed when irradiated by the original laser of the wavelength specifically includes:

[0027] According to the first light energy distribution and the first photoacoustic image corresponding to each light splitting ratio, a second photoacoustic image is determined when the sample to be unmixed is irradiated by the original laser of the wavelength.

[0028] Optionally, there are two light outlets and two light splitting ratios.

[0029] This specification provides a device for photoacoustic imaging demixing, and the device specifically includes:

[0030] An acquisition module, used for acquiring the positional relationship between the sample to be demixed and a plurality of light outlets;

[0031] A first determination module is used to determine, for each wavelength of original laser light of several wavelengths, a first light energy distribution of the sample to be unmixed when the laser light of the wavelength is emitted from the plurality of light outlets to irradiate the sample to be unmixed according to the positional relationship between the sample to be unmixed and the plurality of light outlets, wherein the first light energy distribution represents the light energy distribution of a point inside the sample to be unmixed; and, based on a photoacoustic signal generated when the sample to be unmixed is irradiated by the original laser light of the wavelength, reconstruct a first photoacoustic image corresponding to the sample to be unmixed when it is irradiated by the original laser light of the wavelength;

[0032] A second determination module is used to determine a second photoacoustic image when the sample to be unmixed is irradiated by the original laser of the wavelength according to the first light energy distribution and the first photoacoustic image;

[0033] The unmixing module is used to unmix the sample to be unmixed according to the determined second photoacoustic image when the sample to be unmixed is irradiated by the original laser of each wavelength.

[0034] Optionally, the first determination module is specifically used to determine the light energy of the laser of the wavelength emitted by each light outlet; and, based on the positional relationship between the sample to be unmixed and each light outlet, determine the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet; based on the light energy of the laser of the wavelength emitted by each light outlet and the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet, determine the second light energy distribution of the sample to be unmixed corresponding to each light outlet when each light outlet emits the laser of the wavelength to irradiate the sample to be unmixed; based on the determined second light energy distributions, determine the first light energy distribution of the sample to be unmixed; wherein, the second light energy distribution characterizes the light energy distribution of a point inside the sample to be unmixed when each light outlet emits the laser of the wavelength to irradiate the sample to be unmixed.

[0035] Optionally, the first determination module is specifically used to obtain the optical energy of the original laser of the wavelength; for each of the plurality of splitting ratios, under the splitting ratio, according to the optical energy of the original laser of the wavelength, determine the optical energy of the laser of the wavelength emitted from each light outlet, wherein the splitting ratio is used to divide the original laser according to a preset energy ratio, and the number of lasers obtained by the division is consistent with the number of light outlets;

[0036] The first determination module is specifically configured to determine the first light energy distribution of the sample to be unmixed under the light splitting ratio based on the determined second light energy distributions.

[0037] Optionally, the first determination module is specifically used to determine the light attenuation distribution of a point inside the sample to be unmixed relative to a point on the surface of the sample to be unmixed corresponding to the light outlet; and determine a second light energy distribution of the sample to be unmixed when the light outlet emits the laser of the wavelength to irradiate the sample to be unmixed based on the light attenuation distribution, the light energy of the laser of the wavelength emitted from the light outlet, and the light intensity distribution on the surface of the sample to be unmixed corresponding to the light outlet.

[0038] Optionally, the first determination module is specifically configured to, for each of several splitting ratios, reconstruct a first photoacoustic image corresponding to the sample to be unmixed when irradiated by the original laser of the wavelength under the splitting ratio.

[0039] Optionally, the second determination module is specifically configured to determine, according to the first light energy distribution corresponding to each splitting ratio and the first photoacoustic image, a second photoacoustic image when the sample to be unmixed is irradiated by the original laser of the wavelength.

[0040] Optionally, there are two light outlets and two light splitting ratios.

[0041] This specification provides a computer-readable storage medium, wherein the storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned photoacoustic imaging unmixing method is implemented.

[0042] The present specification provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the above-mentioned photoacoustic imaging unmixing method when executing the program.

[0043] At least one of the above technical solutions adopted in this specification can achieve the following beneficial effects:

[0044] It can be seen from the method of photoacoustic imaging unmixing provided in this specification that when using original lasers of multiple wavelengths to irradiate the sample to be unmixed, for each wavelength of original laser, a number of light outlets are determined to emit lasers of this wavelength to irradiate the sample to be unmixed. The first light energy distribution of the sample to be unmixed, the first light energy distribution characterizes the light energy distribution of the points inside the sample to be unmixed. Then, based on the first light energy distribution of the sample to be unmixed at this wavelength and the first photoacoustic image, the second photoacoustic image at this wavelength can be determined. Thus, the sample to be unmixed can be unmixed according to the second photoacoustic image corresponding to the original laser of each wavelength. This method takes into account the different light energy distributions inside the sample to be unmixed under the irradiation of lasers of different wavelengths, thereby improving the accuracy of photoacoustic imaging unmixing based on multiple wavelengths. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] The drawings described herein are used to provide a further understanding of this specification and constitute a part of this specification. The illustrative embodiments and descriptions of this specification are used to explain this specification and do not constitute an improper limitation on this specification. In the drawings:

[0046] Figure 1 A schematic diagram of a method for photoacoustic imaging demixing provided in this specification;

[0047] Figure 2 A schematic diagram of a photoacoustic imaging unmixing system provided in this specification;

[0048] Figure 3 A schematic diagram of establishing a three-dimensional Cartesian coordinate system on a sample to be unmixed provided in this specification;

[0049] Figure 4 A schematic diagram of an experimental sample provided for this specification;

[0050] Figure 5The true value diagram of the content of CuSO4 in the experimental samples provided in this manual;

[0051] Figure 6 A schematic diagram of the content of CuSO4 in the experimental sample obtained based on the traditional photoacoustic imaging unmixing method provided in this specification;

[0052] Figure 7 A schematic diagram of the content of CuSO4 in the experimental sample obtained based on the photoacoustic imaging unmixing method provided in this specification;

[0053] Figure 8 A schematic diagram of a device for photoacoustic imaging demixing provided in this specification;

[0054] Fig. 9 The corresponding Figure 1 Schematic diagram of electronic equipment. DETAILED DESCRIPTION

[0055] In order to make the purpose, technical solutions and advantages of this specification more clear, the technical solutions of this specification will be clearly and completely described below in combination with the specific embodiments of this specification and the corresponding drawings. Obviously, the described embodiments are only part of the embodiments of this specification, not all of the embodiments. Based on the embodiments in this specification, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the specification.

[0056] The technical solutions provided by the embodiments of this specification are described in detail below in conjunction with the accompanying drawings.

[0057] When using original lasers of multiple wavelengths to unmix samples, due to the different absorption coefficients and scattering coefficients of the samples for lasers of different wavelengths, the light distribution in the samples is different under the irradiation of original lasers of different wavelengths. At present, the light distribution corresponding to the sample at different wavelengths cannot be directly measured by instruments. Therefore, if the problem that the light energy distribution of the sample to be unmixed is different under the irradiation of original lasers of different wavelengths is not considered, the proportion of the components of the unmixed sample to be unmixed will be inaccurate. Based on this, this specification provides a method for photoacoustic imaging unmixing.

[0058] In this specification, the sample can be a biological tissue (component), such as muscle tissue, blood, or a non-biological component, such as soil, gas. This specification does not impose any specific restrictions, and the specific sample can be determined according to needs.

[0059] like Figure 2 As shown, it is a schematic diagram of a photoacoustic imaging unmixing system provided in this specification. Figure 2It can be seen that the number ① is the light outlet, ② is the non-fixed ratio beam splitter, ③ is the fixed ratio beam splitter, ④ is the multi-wavelength tunable laser, ⑤ is the light energy meter, ⑥ is the main control device for collecting and processing data, ⑥ can be an imaging ultrasonic transducer and ancillary devices, and ⑦ is the sample to be unmixed. Among them, the light outlet is connected to the non-fixed ratio beam splitter, the fixed ratio beam splitter and the multi-wavelength tunable laser through the optical transmitter, and the laser generated by the multi-wavelength tunable laser is emitted from each light outlet after passing through the non-fixed ratio beam splitter and the fixed ratio beam splitter, and then irradiated onto the sample to be unmixed. The non-fixed ratio beam splitter can be set with different splitting ratios. In this specification, there are at least two splitting ratios. For each light splitting ratio, the non-fixed ratio beam splitter can divide the original laser based on the total light energy of the original laser, so that each divided sub-laser is transmitted to each light outlet through each optical transmitter, and then emitted from each light outlet. For example: Assuming there are two light outlets, and assuming the light splitting ratio is 2:1, the light energy of one light outlet is 1 / 3 of the total light energy of the original laser, and the light energy of the other light outlet is 2 / 3 of the total light energy of the original laser. The fixed ratio optical beam splitter is used to divide the original laser into two parts according to a specified ratio, one part is transmitted to the non-fixed ratio optical splitter, and the other part is transmitted to the optical energy meter. The optical energy meter can measure the light energy of a part of the laser split and / or can measure the light intensity of a part of the laser split, and then according to the splitting ratio of the fixed ratio optical beam splitter, the total light energy (total light intensity) of the original laser can be obtained. In other words, the light energy of the original laser can be measured by a fixed ratio optical beam splitter and an optical energy meter. Of course, how to determine the total light energy of the original laser generated by a multi-wavelength tunable laser is a relatively mature technology. The above method is only an example, and this specification does not limit it. The imaging ultrasonic transducer and its accessories include a planar array ultrasonic transducer probe, an acquisition module, a data processing module, etc., which are used to collect and process data, etc. The imaging ultrasonic transducer and its accessories are currently available equipment and will not be described in detail in this manual.

[0060] Above Figure 2 The photoacoustic imaging unmixing system shown for photoacoustic imaging unmixing is only an example. When the total light energy of the original laser is measured by other methods, Figure 2 ③ and ⑤ are not necessary components of the photoacoustic imaging unmixing system. In this specification, the photoacoustic imaging unmixing system for performing photoacoustic imaging unmixing includes at least: at least two light outlets, that is, there are several light outlets, a non-fixed ratio optical beam splitter (including at least two splitting ratios), a multi-wavelength tunable laser (used to generate original lasers of different wavelengths), and a main control device (used to collect and process data).

[0061] The execution subject of the technical solution of this specification can be the main control device in the above-mentioned photoacoustic imaging unmixing system, or any electronic device with computing capabilities, such as a server, a terminal, etc., and this specification does not make specific restrictions. The technical solution of this specification is described in detail below using the main control device in the above-mentioned photoacoustic imaging unmixing system as an example.

[0062] like Figure 1 As shown, Figure 1 A schematic flow chart of a method for photoacoustic imaging demixing provided in this specification may specifically include the following steps:

[0063] S100: Obtaining the positional relationship between the sample to be demixed and a plurality of light outlets.

[0064] S102: for each wavelength of original laser light among several wavelengths of original laser light, determine, according to the positional relationship between the sample to be unmixed and the several light outlets, a first light energy distribution of the sample to be unmixed when the several light outlets emit the laser light of the wavelength to irradiate the sample to be unmixed, wherein the first light energy distribution characterizes the light energy distribution of points inside the sample to be unmixed; and, according to the photoacoustic signal generated when the sample to be unmixed is irradiated by the original laser light of the wavelength, reconstruct a first photoacoustic image corresponding to the sample to be unmixed when it is irradiated by the original laser light of the wavelength.

[0065] In this specification, first, the main control device can obtain the positional relationship between the sample to be unmixed and several light outlets. Then, the main control device can determine the first light energy distribution of the sample to be unmixed when the laser of the wavelength is emitted from several light outlets to irradiate the sample to be unmixed, for each wavelength of the original laser of several wavelengths, according to the positional relationship between the sample to be unmixed and the several light outlets. Among them, the first light energy distribution represents the light energy distribution of the point inside the sample to be unmixed, or in other words, the first light energy distribution represents the overall light energy distribution inside the sample to be unmixed when several light outlets emit lasers of the wavelength to irradiate the sample to be unmixed. When the laser is irradiated on the sample to be unmixed, the laser is emitted from the surface of the sample to be unmixed into the interior of the sample to be unmixed, and the light is absorbed and / or scattered by the sample to be unmixed, and the light presents a certain distribution inside the sample to be unmixed, and the first light energy distribution is the light energy distribution presented by the laser of the wavelength inside the sample to be unmixed.

[0066] Specifically, in this specification, when determining the first light energy distribution of the sample to be unmixed, the main control device may first determine the light energy of the laser of the wavelength emitted from each light outlet. And, the main control device may obtain the light energy of the original laser of the wavelength. In order to determine the light energy distribution of the sample to be unmixed under the irradiation of lasers of different wavelengths, in this specification, there are at least two light outlets and at least two splitting ratios. As mentioned above, the splitting ratio is used to divide the original laser according to a preset energy ratio, and the number of divided lasers is consistent with the number of light outlets, so that each divided laser can be emitted from the corresponding light outlet. Then, the main control device can obtain several splitting ratios for dividing the original laser of the wavelength, so as to determine the light energy of the laser of the wavelength emitted from each light outlet according to the light energy of the original laser of the wavelength for each of the several splitting ratios under the splitting ratio. Of course, the light energy of the laser of this wavelength emitted from each light outlet can also be obtained by other methods. For example, the light energy of the laser emitted from each light outlet can be directly measured by existing instruments, as long as the original laser of each wavelength is divided into at least two splitting ratios to each light outlet. That is to say, in this specification, for each wavelength of original laser, at each splitting ratio, the light outlet can emit the laser of this wavelength after being acted upon by the splitting ratio, and then irradiate the sample to be demixed.

[0067] In this specification, when the main control device determines the light intensity distribution on the surface of the sample to be unmixed according to the positional relationship between the sample to be unmixed and each light outlet, a coordinate system can be established in advance based on the sample to be unmixed. Figure 3 As shown, it is a schematic diagram of establishing a three-dimensional Cartesian coordinate system on the sample to be unmixed provided in the specification of this application, that is, a three-dimensional Cartesian coordinate system can be established in the space occupied by the sample to be unmixed. Furthermore, the distance between the end face of the light outlet of each light outlet and the surface of the sample to be unmixed can be preset to a specified distance. In this specification, the specified distance can be preset to 1mm, and the optical conductor in the photoacoustic imaging unmixing system where the main control device is located can be preset to a quartz fiber bundle with a length of 1m and a through diameter of 6mm. The optical fiber filaments in the quartz fiber bundle are randomly mapped, so that the light intensity distribution of the laser emitted from each light outlet on the surface of the sample to be unmixed can be considered to be uniformly distributed. Figure 3 , assuming that the sample to be unmixed is a cube with a side length of 40 mm, there are two light outlets, and the center of each light outlet coincides with the center of the surface of the sample to be unmixed illuminated by the light outlet, that is, the center of the light outlet end face of each light outlet is on the y-axis, such as Figure 3 In the three-dimensional Cartesian coordinate system, the end face of each light outlet is circular. Then the illumination area formed by the lasers of the same wavelength emitted from the two light outlets on the surface of the sample to be demixed is as follows: Figure 3As shown by the two shadows in the figure, the light intensity distribution on the surface of the sample to be unmixed is:

[0068] ,

[0069] .

[0070] in, Represents the light intensity distribution on the surface of the sample to be unmixed corresponding to a light outlet (light outlet 1), It represents the light intensity distribution on the surface of the sample to be unmixed corresponding to another light outlet (light outlet 2), L is the side length of the cube, and d is the diameter of the light outlet end face of the light outlet.

[0071] Of course, according to the positional relationship between the sample to be unmixed and each light outlet, the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet is determined to be only Figure 3 Taking the example of an exemplary description, this specification does not make specific restrictions on how to determine the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet, and there are many solutions at present, such as using optical simulation software, etc., which will not be repeated in this specification. Of course, you can also directly use a light intensity detector to scan the surface of the sample to be unmixed to obtain the light intensity distribution on the surface of the sample to be unmixed, and then based on the positional relationship between the sample to be unmixed and each light outlet, obtain the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet. In actual applications, the relative position of the sample to be unmixed and each light outlet, the establishment of the coordinate system, etc. can all be determined according to actual needs and specific scenarios, so that the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet can also be determined based on actual needs and specific scenarios.

[0072] Therefore, the main control device obtains the light energy of the laser of the wavelength emitted by each light outlet and the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet when each light outlet emits the laser of the wavelength onto the sample to be unmixed, so that the main control device can determine the second light energy distribution of the sample to be unmixed corresponding to each light outlet based on the light energy of the laser of the wavelength emitted by each light outlet and the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet. Therefore, the main control device can obtain the first light energy distribution of the sample to be unmixed when it is irradiated by the laser of the wavelength emitted by several light outlets based on the determined second light energy distribution of the sample to be unmixed corresponding to each light outlet. Among them, in this specification, as mentioned above, the first light energy distribution represents the total light energy distribution of the points inside the sample to be unmixed, that is, the light energy distribution present inside the sample to be unmixed when the sample to be unmixed is irradiated by the original laser of that wavelength, and the second light energy distribution represents the light energy distribution when the sample to be unmixed is irradiated by the laser of that wavelength emitted by each light outlet, that is, the second light energy distribution refers to, for each light outlet, the light energy distribution present inside the sample to be unmixed when the light outlet emits the laser of that wavelength to irradiate the sample to be unmixed.

[0073] Specifically, when the main control device determines the second light energy distribution of the sample to be unmixed corresponding to each light outlet when each light outlet emits a laser of the wavelength to irradiate the sample to be unmixed, the main control device can determine the light attenuation distribution of the points inside the sample to be unmixed relative to the points on the surface of the sample to be unmixed corresponding to each light outlet for each of the several light splitting ratios under the light splitting ratio. In fact, the light attenuation distribution of the points inside the sample to be unmixed relative to the points on the surface of the sample to be unmixed can be expressed as the light attenuation coefficient with the natural logarithm e as the base. An exponentially decaying distribution raised to the power of the inverse of the product of the paths, where is the wavelength of the original laser, and the path refers to the displacement vector between the point inside the sample to be demixed and the point on the surface of the sample to be demixed. Figure 3 In the coordinate system established in the above example, the side length L of the cube is 40 mm, the point on the surface of the sample to be unmixed corresponding to light outlet 1 is (0,0,0), that is, the center point of the light spot presented by light outlet 1 on the surface of the sample to be unmixed, the point on the surface of the sample to be unmixed corresponding to light outlet 2 is (0,40,0), that is, the center point of the light spot presented by light outlet 2 on the surface of the sample to be unmixed, and the point inside the sample to be unmixed is (x,y,z). Then the light attenuation distribution of the point inside the sample to be unmixed relative to the point on the surface of the sample to be unmixed corresponding to the light outlet 1 is as follows: As shown, the light attenuation distribution of the point inside the sample to be unmixed relative to the point on the surface of the sample to be unmixed corresponding to the light outlet 2 is as follows: As shown:

[0074] ,

[0075] .

[0076] Therefore, the main control device can determine the second light energy distribution of the sample to be unmixed corresponding to each light outlet according to the light attenuation distribution on the surface of the sample to be unmixed corresponding to each light outlet. In this specification, for each light outlet, the second light energy distribution of the sample to be unmixed corresponding to the light outlet can be expressed as the product of the light attenuation distribution, the light energy of the laser of the wavelength emitted by the light outlet, and the light intensity distribution on the surface of the sample to be unmixed corresponding to the light outlet. Continuing with the above example, the light energy distribution corresponding to light outlet 1 can be expressed as the following formula As shown, the light attenuation distribution of the point inside the sample to be unmixed relative to the point on the surface of the sample to be unmixed corresponding to the light outlet 2 is as follows: As shown, Indicates that the wavelength emitted from light outlet 1 is The laser light energy, Indicates that the wavelength of the light emitted from the light outlet 2 is The light energy of the laser.

[0077] ,

[0078] .

[0079] It should be noted that, for each of the several splitting ratios, when the light energy of the laser beam of the wavelength emitted from each light outlet is determined by the light energy of the original laser beam of the wavelength and the splitting ratio, then under the specified splitting ratio, the wavelength emitted from the i-th light outlet is The laser light energy It can be expressed by the following formula:

[0080] .

[0081] in, It refers to the wavelength emitted by the i-th light outlet (or light outlet i) obtained according to the specified splitting ratio. The ratio of the laser light energy to the original laser light energy, The wavelength is The original laser light energy, It refers to the light intensity transmission coefficient determined based on comprehensive factors such as the transmittance of the light conductor corresponding to the i-th light outlet and the error in measuring the light energy of the original laser.

[0082] Of course, as mentioned above, the light energy of the laser of this wavelength emitted from each light outlet can be calculated by other methods, or can be directly measured by other instruments. And, as mentioned above, the light energy of the original laser can be measured by a fixed ratio beam splitter and a light energy meter. When the light energy of the original laser is measured by a fixed ratio beam splitter and a light energy meter, the error in measuring the light energy of the original laser can be the response of the light energy meter. Of course, other methods can also be used to determine the light energy of the original laser, and this specification does not make specific restrictions.

[0083] In addition, since the light energy of the laser emitted from each light outlet is obtained after the original laser of the wavelength is acted on by a specified splitting ratio, and there are at least two splitting ratios, for each splitting ratio, the main control device can determine the second light energy distribution of the sample to be unmixed corresponding to each light outlet when each light outlet emits the laser of the wavelength to irradiate the sample to be unmixed under the splitting ratio. Thus, the main control device can determine the first light energy distribution of the sample to be unmixed based on the second light energy distribution of the sample to be unmixed corresponding to each light outlet under the splitting ratio. Specifically, under the splitting ratio, the first light energy distribution can be expressed as the sum of the second light energy distribution of the sample to be unmixed corresponding to each light outlet. Continuing with the above example and formula, under the specified splitting ratio, the first light energy distribution It can be shown as the following formula:

[0084] .

[0085] in, It means that under the specified splitting ratio, the wavelength of the n light outlets is When the laser irradiates the sample to be unmixed, the first light energy distribution of the sample to be unmixed, n represents the total number of light outlets, i represents light outlet i or the i-th light outlet, as mentioned above, It means that under the specified splitting ratio, the wavelength emitted by the i-th light outlet is When the laser irradiates the sample to be unmixed, the second light energy distribution of the sample to be unmixed corresponding to the i-th light outlet.

[0086] In addition, the main control device can also reconstruct the first photoacoustic image corresponding to the sample to be unmixed when it is irradiated with the original laser of the wavelength according to the photoacoustic signal generated when the sample to be unmixed is irradiated with the original laser of the wavelength. How to reconstruct the photoacoustic image based on the photoacoustic signal is already a relatively mature technology, and this specification will not be described in detail. For example, the main control device can obtain the photoacoustic signal or collect the photoacoustic signal based on the pre-installed sensor, and then use the back projection algorithm or the Fourier transform algorithm, etc., to reconstruct the photoacoustic image.

[0087] As mentioned above, in this specification, the original laser of each wavelength needs to be subjected to the action of at least two splitting ratios, and then emitted through the light outlet to irradiate the sample to be unmixed. Therefore, the main control device can reconstruct the first photoacoustic image corresponding to the sample to be unmixed when being irradiated by the original laser of the wavelength under the splitting ratio for each splitting ratio in the several splitting ratios. Therefore, based on the above steps, the main control device can obtain the first photoacoustic images corresponding to the original laser of each wavelength in the several wavelengths of original lasers, when being irradiated by the original laser of the wavelength under the several splitting ratios. That is, for each wavelength of original laser of the several wavelengths of original lasers, for each splitting ratio in the several splitting ratios, the main control device can reconstruct the first photoacoustic image corresponding to the sample to be unmixed when being irradiated by the original laser of the wavelength under the splitting ratio according to the photoacoustic signal when the laser of the wavelength is emitted from each light outlet to irradiate the sample to be unmixed under the splitting ratio.

[0088] In addition, it should be noted that in this specification, the lasers emitted from each light outlet can be irradiated on the sample to be demixed at the same time, and the lasers emitted from each light outlet irradiate different positions on the sample to be demixed.

[0089] S104: Determine, according to the first light energy distribution and the first photoacoustic image, a second photoacoustic image when the sample to be demixed is irradiated by the original laser of the wavelength.

[0090] S106: Unmixing the sample to be unmixed according to the determined second photoacoustic image when the sample to be unmixed is irradiated by the original laser of each wavelength.

[0091] At a specific wavelength, the photoacoustic image reconstructed based on the photoacoustic signal can be expressed as the product of the real light absorption image and the total light energy distribution inside the sample to be unmixed. Therefore, in this specification, the main control device can determine the second photoacoustic image of the sample to be unmixed when it is irradiated by the original laser of this wavelength based on the first light energy distribution and the first photoacoustic image.

[0092] Based on the above steps S100 to S102, for each wavelength of the original laser light among the several wavelengths of the original laser light, the main control device has obtained the first photoacoustic images corresponding to the sample to be unmixed when irradiated by the original laser light of the wavelength under several split ratios. Therefore, the main control device can determine the second photoacoustic image of the sample to be unmixed when irradiated by the original laser light of the wavelength under several split ratios according to the first photoacoustic images corresponding to the sample to be unmixed when irradiated by the original laser light of the wavelength under the split ratio, and the first light energy distribution when the sample to be unmixed is irradiated by the original laser light of the wavelength under the split ratio. Using the above example and formula, the relationship between the first photoacoustic image and the second photoacoustic image can be expressed by the following formula:

[0093] .

[0094] in, Indicates that when the splitting ratio is specified When the sample to be unmixed is The first photoacoustic image is obtained by reconstructing the photoacoustic signal generated during the original laser irradiation. Indicates that the sample to be unmixed is The second photoacoustic image when the original laser is irradiated with a wavelength of The real light absorption image when the original laser is irradiated. Indicates that when the splitting ratio is specified When the sample to be unmixed is The first light energy distribution during the original laser irradiation.

[0095] It should be noted that the real light absorption image refers to the second photoacoustic image, which reflects the absorption characteristics of the sample to be unmixed to light, which is closer to the actual situation than the first photoacoustic image. That is to say, the second photoacoustic image can more accurately reflect the absorption characteristics and scattering characteristics of the sample to be unmixed to light than the first photoacoustic image. Under ideal conditions, there is an ideal photoacoustic image for the original laser of the specified wavelength, which can fully and accurately reflect the absorption characteristics and scattering characteristics of the sample to be unmixed to the original laser of the specified wavelength. When the light energy distribution inside the sample to be unmixed can be accurately obtained, the ideal photoacoustic image can be obtained. However, it is currently impossible to directly measure or accurately obtain the real light energy distribution of the sample to be unmixed under the irradiation of the laser of the specified wavelength through an instrument. The reason why the real light energy distribution cannot be known is that the light attenuation coefficient (that is, the absorption characteristics and scattering characteristics of light) of the sample to be unmixed to lasers of different wavelengths is different, and the light attenuation coefficient cannot be directly measured or accurately calculated by the instrument. In this specification, based on the above steps S100~S106, the light attenuation coefficient of the sample to be unmixed under the irradiation of each wavelength of laser can be obtained. , so that a relatively accurate light energy distribution inside the sample to be unmixed can be obtained, so that a more realistic second photoacoustic image can be obtained relative to the first photoacoustic image.

[0096] Finally, the main control device can unmix the sample to be unmixed according to the second photoacoustic image when the sample to be unmixed is irradiated by the original laser of each wavelength, that is, determine the proportion of each component in the sample to be unmixed. How to unmix the sample to be unmixed based on the second photoacoustic image is already a relatively mature technology, and this manual will not go into details.

[0097] based on Figure 1 The above-mentioned photoacoustic imaging unmixing method provided in this specification shown in the specification, when using original lasers of multiple wavelengths to irradiate the sample to be unmixed, for each wavelength of original laser, determine the first light energy distribution of the sample to be unmixed when the laser of this wavelength is emitted from several light outlets to irradiate the sample to be unmixed, and the first light energy distribution characterizes the light energy distribution of the points inside the sample to be unmixed. Then, based on the first light energy distribution of the sample to be unmixed at this wavelength and the first photoacoustic image, the second photoacoustic image at this wavelength can be determined. Therefore, the sample to be unmixed can be unmixed according to the second photoacoustic image corresponding to the original laser of each wavelength. This method takes into account the problem that the light energy distribution inside the sample to be unmixed is actually uneven under the irradiation of lasers of different wavelengths, thereby improving the accuracy of photoacoustic imaging unmixing based on multiple wavelengths.

[0098] Furthermore, the above method has been verified through experiments, and the unmixing result obtained by using the second photoacoustic image obtained based on the above method to unmix the sample to be unmixed is more accurate than that obtained by directly using the first photoacoustic image to unmix the sample to be unmixed. Figure 4 As shown, it is a schematic diagram of an experimental sample provided in this specification. It can be seen that the experimental sample includes four capillaries 101~104, and when verified by experiment, the experimental sample is in the shape of a cube with a side length of 40mm, and the experimental sample is prepared by preparing 64ml of agar solution with a mass fraction of 1.5%, and adding 3ml of fresh milk and 0.5ml of ink to solidify, so as to simulate the scattering and absorption of light by biological tissue. In addition, the four capillaries are tetrafluoroethylene capillaries, and their inner diameters are all 50μm. The four capillaries are arranged radially at an angular interval of 45°. One end of each capillary can be connected to a needle and injected with liquid. During the experiment, mixed solutions of CuSO4 and NiSO4 in different proportions (pure water as solvent) are respectively injected into the four capillaries to verify the above-mentioned photoacoustic imaging demixing method. Specifically, NiSO4 and CuSO4 solutions with a molar concentration of 0.5mol / L can be prepared, and in one capillary ( Figure 4 Inject CuSO4 and NiSO4 solution in a volume ratio of 1:9 into a capillary tube (101). Figure 4 Only CuSO4 solution is injected into the 102) in the capillary tube, that is, the true value of CuSO4 content is 100%. Figure 4 103) in the injection of CuSO4 and NiSO4 solution in a volume ratio of 3:1 mixed solution, that is, the true value of CuSO4 content is 75%, in a capillary ( Figure 4 A mixed solution of CuSO4 and NiSO4 solution with a volume ratio of 1:3 is injected into the 104) in the experiment, that is, the true value of CuSO4 content is 25%. In addition, four capillaries can be fixed at the center of the experimental sample by mold casting. During the experiment, after the experimental sample is made in advance, Figure 2 The photoacoustic imaging unmixing system shown in the figure uses the method described in the above steps S100 to S106 to unmix the experimental sample.

[0099] like Figure 5 As shown, Figure 5 The true value diagram of the content of CuSO4 in the experimental samples provided in this manual is Figure 5 In the figure, the value of the maximum value projection along the z direction represents the content of CuSO4. It can be seen that as mentioned above, the true values ​​of the CuSO4 contents corresponding to the four capillaries 401-404 are 10%, 100%, 75%, and 25%, respectively. Figure 6The schematic diagram of the content of CuSO4 in the experimental sample obtained by the traditional photoacoustic imaging unmixing method provided in this specification, and the value of the maximum projection along the z direction represents the content of CuSO4. Figure 6 It can be seen that the unmixing results of the CuSO4 content in the four capillaries 101-104 are almost close to 100%, which is quite different from the true value. Figure 7 A schematic diagram of the content of CuSO4 in the experimental sample obtained based on the photoacoustic imaging unmixing method provided in this specification, that is, Figure 7 The value obtained by the photoacoustic imaging unmixing method based on the above steps S100 to S106, and the value of the maximum projection along the z direction represents the content of CuSO4. Figure 7 It can be seen that the unmixing results of the CuSO4 contents corresponding to the four capillaries 101-104 are 9.6%, 99.5%, 78.3%, and 26.3%, respectively, which are relatively close to the true values. Obviously, the photoacoustic imaging unmixing method described in the above steps S100-S106 improves the accuracy of unmixing.

[0100] Based on the above-described method for photoacoustic imaging unmixing, the present specification also provides a schematic diagram of a device for photoacoustic imaging unmixing. Figure 8 shown.

[0101] Figure 8 A schematic diagram of a device for photoacoustic imaging demixing provided in an embodiment of this specification, the device specifically includes:

[0102] An acquisition module 800 is used to acquire the positional relationship between the sample to be demixed and a plurality of light outlets;

[0103] A first determination module 802 is used to determine, for each wavelength of original laser light of several wavelengths, a first light energy distribution of the sample to be unmixed when the plurality of light outlets emit the laser light of the wavelength to irradiate the sample to be unmixed according to the positional relationship between the sample to be unmixed and the plurality of light outlets, wherein the first light energy distribution represents the light energy distribution of a point inside the sample to be unmixed; and, based on a photoacoustic signal generated when the sample to be unmixed is irradiated by the original laser light of the wavelength, reconstruct a first photoacoustic image corresponding to the sample to be unmixed when the sample to be unmixed is irradiated by the original laser light of the wavelength;

[0104] A second determination module 804 is used to determine a second photoacoustic image when the sample to be demixed is irradiated by the original laser of the wavelength according to the first light energy distribution and the first photoacoustic image;

[0105] The unmixing module 806 is configured to unmix the sample to be unmixed according to the determined second photoacoustic image when the sample to be unmixed is irradiated by the original laser of each wavelength.

[0106] Optionally, the first determination module 802 is specifically used to determine the light energy of the laser of the wavelength emitted by each light outlet; and, based on the positional relationship between the sample to be unmixed and each light outlet, determine the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet; based on the light energy of the laser of the wavelength emitted by each light outlet and the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet, determine the second light energy distribution of the sample to be unmixed corresponding to each light outlet when each light outlet emits the laser of the wavelength to irradiate the sample to be unmixed; based on the determined second light energy distributions, determine the first light energy distribution of the sample to be unmixed; wherein, the second light energy distribution characterizes the light energy distribution of a point inside the sample to be unmixed when each light outlet emits the laser of the wavelength to irradiate the sample to be unmixed.

[0107] Optionally, the first determination module 802 is specifically used to obtain the optical energy of the original laser of the wavelength; for each splitting ratio among several splitting ratios, under the splitting ratio, according to the optical energy of the original laser of the wavelength, determine the optical energy of the laser of the wavelength emitted from each light outlet, wherein the splitting ratio is used to divide the original laser according to a preset energy ratio, and the number of lasers obtained by the division is consistent with the number of light outlets;

[0108] The first determination module 802 is specifically configured to determine the first light energy distribution of the sample to be unmixed under the light splitting ratio based on the determined second light energy distributions.

[0109] Optionally, the first determination module 802 is specifically used to determine the light attenuation distribution of a point inside the sample to be unmixed relative to a point on the surface of the sample to be unmixed corresponding to the light outlet; based on the light attenuation distribution, the light energy of the laser of the wavelength emitted by the light outlet, and the light intensity distribution on the surface of the sample to be unmixed corresponding to the light outlet, determine the second light energy distribution of the sample to be unmixed when the laser of the wavelength is emitted from the light outlet to irradiate the sample to be unmixed.

[0110] Optionally, the first determining module 802 is specifically configured to, for each of several splitting ratios, reconstruct a first photoacoustic image corresponding to when the sample to be unmixed is irradiated by the original laser of the wavelength under the splitting ratio.

[0111] Optionally, the second determination module 804 is specifically configured to determine, according to the first light energy distribution corresponding to each splitting ratio and the first photoacoustic image, a second photoacoustic image when the sample to be unmixed is irradiated by the original laser of the wavelength.

[0112] Optionally, there are two light outlets and two light splitting ratios.

[0113] The embodiments of the present specification also provide a computer-readable storage medium, which stores a computer program. The computer program can be used to execute the photoacoustic imaging unmixing method described above.

[0114] Based on the above-mentioned method for photoacoustic imaging unmixing, the present specification also proposes Fig. 9 The schematic structure diagram of the electronic device shown in FIG. Fig. 9 At the hardware level, the electronic device includes a processor, an internal bus, a network interface, a memory, and a non-volatile memory, and may also include hardware required for other services. The processor reads the corresponding computer program from the non-volatile memory into the memory and then runs it to implement the above-mentioned method of photoacoustic imaging unmixing.

[0115] Of course, in addition to software implementation, this specification does not exclude other implementation methods, such as logic devices or a combination of software and hardware, etc., that is to say, the executor of the following processing flow is not limited to each logic unit, but can also be hardware or logic devices.

[0116] In the 1990s, it was very clear whether the improvement of a technology was hardware improvement (for example, improvement of the circuit structure of diodes, transistors, switches, etc.) or software improvement (improvement of the method flow). However, with the development of technology, many improvements of the method flow today can be regarded as direct improvements of the hardware circuit structure. Designers almost always obtain the corresponding hardware circuit structure by programming the improved method flow into the hardware circuit. Therefore, it cannot be said that the improvement of a method flow cannot be implemented with hardware entity modules. For example, a programmable logic device (PLD) (such as a field programmable gate array (FPGA)) is such an integrated circuit whose logical function is determined by the user's programming of the device. Designers can "integrate" a digital system on a PLD by programming it themselves, without having to ask chip manufacturers to design and make dedicated integrated circuit chips. Moreover, nowadays, instead of manually making integrated circuit chips, this kind of programming is mostly implemented by "logic compiler" software, which is similar to the software compiler used when developing and writing programs, and the original code before compilation must also be written in a specific programming language, which is called hardware description language (HDL). There is not only one kind of HDL, but many kinds, such as ABEL (Advanced Boolean Expression Language), AHDL (Altera Hardware Description Language), Confluence, CUPL (Cornell University Programming Language), HDCal, JHDL (Java Hardware Description Language), Lava, Lola, MyHDL, PALASM, RHDL (Ruby Hardware Description Language), etc. The most commonly used ones are VHDL (Very-High-Speed ​​Integrated Circuit Hardware Description Language) and Verilog. Those skilled in the art should also know that it is only necessary to program the method flow slightly in the above-mentioned hardware description languages ​​and program it into the integrated circuit, and then it is easy to obtain the hardware circuit that implements the logic method flow.

[0117] The controller may be implemented in any suitable manner, for example, the controller may take the form of a microprocessor or processor and a computer-readable medium storing a computer-readable program code (e.g., software or firmware) executable by the (micro)processor, a logic gate, a switch, an application-specific integrated circuit (ASIC), a programmable logic controller, and an embedded microcontroller, examples of which include but are not limited to the following microcontrollers: ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20, and Silicone Labs C8051F320, and the memory controller may also be implemented as part of the control logic of the memory. It is also known to those skilled in the art that, in addition to implementing the controller in a purely computer-readable program code manner, the controller may be implemented in the form of a logic gate, a switch, an application-specific integrated circuit, a programmable logic controller, and an embedded microcontroller by logically programming the method steps. Therefore, such a controller may be considered as a hardware component, and the devices for implementing various functions included therein may also be considered as structures within the hardware component. Or even, the devices for implementing various functions may be considered as both software modules for implementing the method and structures within the hardware component.

[0118] The systems, devices, modules or units described in the above embodiments may be implemented by computer chips or entities, or by products with certain functions. A typical implementation device is a computer. Specifically, the computer may be, for example, a personal computer, a laptop computer, a cellular phone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

[0119] For the convenience of description, the above device is described in various units according to their functions. Of course, when implementing this specification, the functions of each unit can be implemented in the same or multiple software and / or hardware.

[0120] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Furthermore, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0121] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0122] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0123] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0124] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0125] The memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0126] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0127] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity 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, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0128] It should be understood by those skilled in the art that the embodiments of this specification may be provided as methods, systems or computer program products. Therefore, this specification may take the form of a complete hardware embodiment, a complete software embodiment or an embodiment combining software and hardware. Moreover, this specification may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0129] This specification may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. This specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules may be located in local and remote computer storage media, including storage devices.

[0130] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0131] The above description is only an embodiment of this specification and is not intended to limit this specification. For those skilled in the art, this specification may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of this specification should be included in the scope of the claims of this application.

Claims

1. A method for photoacoustic imaging unmixing, characterized in that: The method comprises: Obtaining the positional relationship between the sample to be unmixed and a plurality of light outlets; For each wavelength of original laser light among several wavelengths of original laser light, determine the light energy of the laser light of the wavelength emitted from each light outlet; and, according to the positional relationship between the sample to be unmixed and each light outlet, determine the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet; according to the light energy of the laser light of the wavelength emitted from each light outlet and the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet, determine the second light energy distribution of the sample to be unmixed corresponding to each light outlet when each light outlet emits the laser light of the wavelength to irradiate the sample to be unmixed; based on the determined second light energy distributions, determine the first light energy distribution of the sample to be unmixed; wherein the second light energy distribution represents the light energy distribution of the points inside the sample to be unmixed when each light outlet emits the laser light of the wavelength to irradiate the sample to be unmixed; the first light energy distribution represents the light energy distribution of the points inside the sample to be unmixed; and, according to the photoacoustic signal generated when the sample to be unmixed is irradiated by the original laser light of the wavelength, reconstruct the first photoacoustic image corresponding to the sample to be unmixed when it is irradiated by the original laser light of the wavelength; Determine, according to the first light energy distribution and the first photoacoustic image, a second photoacoustic image when the sample to be unmixed is irradiated by the original laser of the wavelength; Unmixing the sample to be unmixed according to the determined second photoacoustic image when the sample to be unmixed is irradiated by the original laser of each wavelength; Determining the second light energy distribution of the sample to be unmixed when the light outlet emits the laser of the wavelength to irradiate the sample to be unmixed specifically includes: determining the light attenuation distribution of a point inside the sample to be unmixed relative to a point on the surface of the sample to be unmixed corresponding to the light outlet; determining the second light energy distribution of the sample to be unmixed when the light outlet emits the laser of the wavelength to irradiate the sample to be unmixed according to the light attenuation distribution, the light energy of the laser of the wavelength emitted by the light outlet, and the light intensity distribution on the surface of the sample to be unmixed corresponding to the light outlet, wherein the light attenuation distribution is expressed as a light attenuation coefficient with a natural logarithm e as the base. An exponentially decaying distribution raised to the power of the inverse of the product of the paths, is the wavelength of the original laser, and the path refers to the displacement vector between the point inside the sample to be demixed and the point on the surface of the sample to be demixed.

2. The method according to claim 1, characterized in that Determine the light energy of the laser of the wavelength emitted from each light outlet, specifically including: Obtaining the optical energy of the original laser of the wavelength; For each of the several splitting ratios, under the splitting ratio, according to the light energy of the original laser of the wavelength, the light energy of the laser of the wavelength emitted from each light outlet is determined, wherein the splitting ratio is used to divide the original laser according to a preset energy ratio, and the number of lasers obtained by division is consistent with the number of light outlets; Determining the first light energy distribution of the sample to be unmixed based on the determined second light energy distributions specifically includes: Based on the determined second light energy distributions, the first light energy distribution of the sample to be unmixed is determined under the light splitting ratio.

3. The method according to claim 2, characterized in that Reconstructing a first photoacoustic image corresponding to the sample to be unmixed when irradiated by the original laser of the wavelength specifically includes: For each of the several splitting ratios, under the splitting ratio, a first photoacoustic image corresponding to the sample to be unmixed when irradiated by the original laser of the wavelength is reconstructed.

4. The method according to claim 3, characterized in that Determining a second photoacoustic image of the sample to be unmixed when irradiated by the original laser of the wavelength specifically includes: According to the first light energy distribution and the first photoacoustic image corresponding to each light splitting ratio, a second photoacoustic image is determined when the sample to be unmixed is irradiated by the original laser of the wavelength.

5. The method according to claim 2, characterized in that The number of the light outlets is two and the number of the light splitting ratios is two.

6. A device for photoacoustic imaging demixing, characterized in that: The device specifically comprises: An acquisition module, used for acquiring the positional relationship between the sample to be demixed and a plurality of light outlets; The first determination module is used to determine the light energy of the laser light of the wavelength emitted from each light outlet for each of the original laser light of several wavelengths; and determine the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet according to the positional relationship between the sample to be unmixed and each light outlet; determine the second light energy distribution of the sample to be unmixed corresponding to each light outlet when each light outlet emits the laser light of the wavelength to irradiate the sample to be unmixed according to the light energy of the laser light of the wavelength emitted from each light outlet and the light intensity distribution on the surface of the sample to be unmixed corresponding to each light outlet; determine the first light energy distribution of the sample to be unmixed based on the determined second light energy distributions; wherein the second light energy distribution characterizes the light energy distribution of the point inside the sample to be unmixed when each light outlet emits the laser light of the wavelength to irradiate the sample to be unmixed. the first light energy distribution of the sample to be unmixed; and, based on the photoacoustic signal generated when the sample to be unmixed is irradiated by the original laser of the wavelength, reconstructing the first photoacoustic image corresponding to the sample to be unmixed when being irradiated by the original laser of the wavelength; determining the second light energy distribution of the sample to be unmixed when the light outlet emits the laser of the wavelength to irradiate the sample to be unmixed, specifically including: determining the light attenuation distribution of a point inside the sample to be unmixed relative to a point on the surface of the sample to be unmixed corresponding to the light outlet; determining the second light energy distribution of the sample to be unmixed when the light outlet emits the laser of the wavelength to irradiate the sample to be unmixed according to the light attenuation distribution, the light energy of the laser of the wavelength emitted by the light outlet, and the light intensity distribution on the surface of the sample to be unmixed corresponding to the light outlet, wherein the light attenuation distribution is expressed as a light attenuation coefficient with a natural logarithm e as the base. An exponentially decaying distribution raised to the power of the inverse of the product of the paths, is the wavelength of the original laser, and the path refers to the displacement vector between the point inside the sample to be unmixed and the point on the surface of the sample to be unmixed; A second determination module is used to determine a second photoacoustic image when the sample to be unmixed is irradiated by the original laser of the wavelength according to the first light energy distribution and the first photoacoustic image; The unmixing module is used to unmix the sample to be unmixed according to the determined second photoacoustic image when the sample to be unmixed is irradiated by the original laser of each wavelength.

7. A computer-readable storage medium, characterized in that: The storage medium stores a computer program, and when the computer program is executed by a processor, the method described in any one of claims 1 to 5 is implemented.

8. An electronic device, characterized in that: The method comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method described in any one of claims 1 to 5 is implemented.

Images