A digital dummy intensity calculation method using illumination light feedback, electronic device, and medium

Abstract

The application provides a digital phantom intensity calculation method using illumination light feedback, an electronic device and a medium, and comprises the following steps: acquiring a first excitation power P1 and a first excitation wavelength A(λ) of a fluorescence acquisition system; acquiring a fluorescence molecule concentration, a fluorescence spectrum curve, an excitation efficiency curve E(λ) and a fluorescence emission intensity of a fluorescence molecule sample; acquiring a second excitation power P2 and a second excitation wavelength B(λ) of a to-be-tested device; setting a fluorescence concentration value c, obtaining a fluorescence emission spectrum curve corresponding to the to-be-tested device from a relationship W(c) between a fluorescence center wavelength and the fluorescence concentration; calculating a fluorescence emission ratio η(c) based on the second excitation wavelength B(λ) of the to-be-tested device, the excitation efficiency curve E(λ), the first excitation wavelength A(λ) of the fluorescence acquisition system, the second excitation power P2 of the to-be-tested device and the first excitation power P1 of the fluorescence acquisition system; and multiplying the fluorescence emission ratio η(c) and the fluorescence emission intensity I(c) at any point in space to calculate the fluorescence emission intensity of the to-be-tested device at the position.

Description

Technical Field

[0001] This invention relates to the field of biofluorescence imaging technology, and in particular to a digital phantom intensity calculation method, electronic device, and medium utilizing illumination light feedback. Background Technology

[0002] In the current technological field, accurate reproduction of fluorescence properties is crucial for numerous applications, including biomedical imaging, materials science, environmental monitoring, and industrial inspection. However, traditional digital phantom reproduction methods are often limited by the fixed parameters of the fluorescence acquisition system and the inherent characteristics of fluorescent molecules, which can lead to deviations in measurement results under different devices or conditions. Furthermore, due to the complexity of fluorescence signals, quantitative calculation of fluorescence intensity remains a challenge. For example, existing phantom reproduction methods only consider illumination distribution, neglecting factors such as wavelength and excitation efficiency, thus failing to accurately reproduce the light field information of digital phantoms. Summary of the Invention

[0003] The purpose of this invention is to address the shortcomings of existing technologies by providing a digital phantom intensity calculation method, electronic device, and medium that utilizes illumination light feedback.

[0004] In a first aspect, embodiments of the present invention provide a method for calculating the intensity of a digital phantom using illumination light feedback, the method comprising:

[0005] Obtain the first excitation power P1 and the first excitation wavelength A(λ) of the fluorescence acquisition system;

[0006] Obtain the concentration of fluorescent molecules, fluorescence spectrum curve, excitation efficiency curve E(λ), and fluorescence emission intensity of the fluorescent molecule sample;

[0007] Obtain the second excitation power P2 and the second excitation wavelength B(λ) of the device under test;

[0008] Set the fluorescence concentration value c, and obtain the fluorescence emission spectrum curve of the device under test from the relationship between the fluorescence center wavelength and the fluorescence concentration W(c);

[0009] The fluorescence emission ratio η(c) is calculated based on the second excitation wavelength B(λ) of the device under test, the excitation efficiency curve E(λ), the first excitation wavelength A(λ) of the fluorescence acquisition system, the second excitation power P2 of the device under test, and the first excitation power P1 of the fluorescence acquisition system.

[0010] The fluorescence emission ratio η(c) is multiplied by the fluorescence emission intensity I(c) at any point in space to calculate the fluorescence emission intensity of the device under test at that location.

[0011] Secondly, embodiments of the present invention provide an electronic device, including a memory and a processor, wherein the memory is coupled to the processor; wherein the memory is used to store program data, and the processor is used to execute the program data to implement the above-described digital phantom intensity calculation method using illumination light feedback.

[0012] Thirdly, embodiments of the present invention provide a computer-readable storage medium having a computer program stored thereon, wherein the program, when executed by a processor, implements the above-described method for calculating the intensity of a digital phantom using illumination feedback.

[0013] Fourthly, embodiments of the present invention provide a computer program product, including a computer program / instruction, which, when executed by a processor, implements the above-described method for calculating the intensity of a digital phantom using illumination feedback.

[0014] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0015] This invention provides a method for calculating the intensity of a digital phantom using illumination feedback. This invention not only considers the first excitation power, the first two-dimensional planar distribution of the laser, and the first excitation wavelength of the fluorescence acquisition system, but also integrates the characteristics of the fluorescent molecule sample, including the concentration of fluorescent molecules, the fluorescence spectrum curve, the excitation efficiency curve, and the fluorescence emission intensity. By detecting the second excitation power, the second two-dimensional planar distribution of the laser, and the second excitation wavelength corresponding to the device under test, this invention calculates the fluorescence emission intensity of the device under test based on the principle of fluorescence excitation and emission. This invention makes the mapping from physical space to digital space possible, thereby reproducing fluorescence characteristics in the digital phantom, significantly improving its accuracy and reliability. Attached Figure Description

[0016] To more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0017] Figure 1 This is a flowchart of a digital phantom intensity calculation method using illumination light feedback provided in an embodiment of the present invention;

[0018] Figure 2 This is a schematic diagram of the fluorescence acquisition system provided in an embodiment of the present invention;

[0019] Figure 3 This is a schematic diagram of the image calibration process of the fluorescence acquisition system provided in an embodiment of the present invention;

[0020] Figure 4 This is a schematic diagram of the first two-dimensional planar distribution of the laser and the first excitation wavelength corresponding to the fluorescence acquisition system provided in the embodiments of the present invention;

[0021] Figure 5 This is a schematic diagram of the fluorescence spectrum curve and excitation efficiency curve of the fluorescent molecule provided in the embodiments of the present invention;

[0022] Figure 6 This is a schematic diagram of the fluorescence spectrum curve fitting of the fluorescent molecule provided in the embodiment of the present invention;

[0023] Figure 7 This is a schematic diagram of the fluorescence emission intensity fitting of the fluorescent molecule provided in the embodiments of the present invention;

[0024] Figure 8 This is a schematic diagram of the second laser two-dimensional planar distribution and the second excitation wavelength corresponding to the device under test provided in the embodiment of the present invention;

[0025] Figure 9 This is a schematic diagram of the fluorescence emission intensity corresponding to the device under test provided in the embodiments of the present invention;

[0026] Figure 10 This is a schematic diagram of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0027] Exemplary embodiments will now be described in detail, examples of which are illustrated in the accompanying drawings. When the following description relates to the drawings, unless otherwise indicated, the same numerals in different drawings denote the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the present invention. Rather, they are merely examples of apparatuses and methods consistent with some aspects of the invention as detailed in the appended claims.

[0028] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The singular forms “a,” “the,” and “the” used in this invention and the appended claims are also intended to include the plural forms unless the context clearly indicates otherwise. It should also be understood that the term “and / or” as used herein refers to and includes any or all possible combinations of one or more of the associated listed items.

[0029] It should be understood that although the terms first, second, third, etc., may be used in this invention to describe various information, this information should not be limited to these terms. These terms are only used to distinguish information of the same type from one another. For example, first information may also be referred to as second information without departing from the scope of this invention, and similarly, second information may also be referred to as first information. Depending on the context, the word "if" as used herein may be interpreted as "when," "when," or "in response to a determination."

[0030] The present invention will now be described in detail with reference to the accompanying drawings. Unless otherwise specified, the features of the following embodiments and implementations can be combined with each other.

[0031] like Figure 1 As shown, this embodiment of the invention provides a method for calculating the intensity of a digital phantom using illumination light feedback, the method comprising the following steps:

[0032] Step S1: Obtain the first excitation power P1, the first two-dimensional plane distribution of the laser, and the first excitation wavelength A(λ) of the fluorescence acquisition system.

[0033] Specifically, in this example, such as Figure 2 As shown, a monochrome camera and lens assembly are placed and assembled. Within the imaging distance range, the first laser two-dimensional plane distribution of the illumination light is detected by a two-dimensional plane detector.

[0034] At a fixed central position on the imaging surface, the first excitation power P1 of the illumination light is recorded by an optical power meter;

[0035] At a fixed central position on the imaging plane, the first excitation spectrum of the illumination light is recorded by the fiber optic probe of the spectrometer, thereby obtaining the first excitation wavelength A(λ).

[0036] Step S2: Obtain the fluorescent molecule concentration, fluorescence spectrum curve, excitation efficiency curve E(λ), and fluorescence emission intensity of the fluorescent molecule sample.

[0037] The process of obtaining the excitation efficiency curve E(λ) includes:

[0038] The fluorescence acquisition system irradiates fluorescent molecules with a series of first excitation wavelengths A(λ), records the fluorescence emission power, and fits the fluorescence emission power with the first excitation wavelength A(λ) to obtain the relationship function between the fluorescence emission power and the first excitation wavelength A(λ), i.e., the excitation efficiency curve E(λ).

[0039] Step S3: Obtain the second excitation power P2, the second laser two-dimensional planar distribution, and the second excitation wavelength B(λ) of the device under test.

[0040] Specifically, in this example, a monochrome camera and a 780nm bandwidth color filter lens group are placed and assembled. Within the imaging distance range, the second laser two-dimensional plane distribution of the illumination light is detected by a two-dimensional plane detector.

[0041] At a fixed central position on the imaging surface, the second excitation power P2 of the illumination light is recorded by an optical power meter;

[0042] At a fixed central position on the imaging plane, the second excitation spectrum of the illumination light is recorded by the fiber optic probe of the spectrometer, thereby obtaining the second excitation wavelength B(λ).

[0043] Step S4: Set the fluorescence concentration value c, and obtain the fluorescence emission spectrum curve corresponding to the device under test from the relationship between the fluorescence center wavelength and the fluorescence concentration W(c).

[0044] The process of obtaining the relationship between the fluorescence center wavelength and the fluorescence concentration, W(c), includes:

[0045] A series of fluorescent molecule samples with varying concentrations were prepared. Using a spectrometer with a fixed excitation wavelength of 780 nm, the series of fluorescent molecule solutions with varying concentrations were scanned sequentially to obtain fluorescence spectral curves. The fluorescence center wavelength of the fluorescence spectral curves was extracted to obtain discrete sampling points of the fluorescence spectral curve center wavelength as a function of the fluorescent molecule concentration.

[0046] The relationship between the fluorescence center wavelength and fluorescence concentration, W(c), was obtained by fitting the fluorescence center wavelength with the fluorescence concentration.

[0047] Step S5: Calculate the fluorescence emissivity η(c) based on the second excitation wavelength B(λ) of the device under test, the excitation efficiency curve E(λ), the first excitation wavelength A(λ) of the fluorescence acquisition system, the second excitation power P2 of the device under test, and the first excitation power P1 of the fluorescence acquisition system.

[0048] The relative fluorescence excitation amount of the fluorescence acquisition system is obtained by calculating the integral of the first excitation wavelength A(λ) and the excitation efficiency curve E(λ) of the fluorescence acquisition system.

[0049] The relative fluorescence excitation amount of the device under test is obtained by calculating the integral of the second excitation wavelength B(λ) and the excitation efficiency curve E(λ).

[0050] The ratio of the relative fluorescence excitation amount corresponding to the fluorescence acquisition system to the relative fluorescence excitation amount corresponding to the device under test is taken as the relative fluorescence emission ratio.

[0051] To obtain the absolute fluorescence emission ratio, the relative fluorescence emission ratio needs to be normalized to the excitation power. The ratio of the second excitation power P2 of the device under test to the first excitation power P1 of the fluorescence acquisition system is calculated, and this ratio is multiplied by the relative fluorescence emission ratio to obtain the absolute fluorescence emission ratio η(c).

[0052] Furthermore, the expression is as follows:

[0053]

[0054] Step S6: Multiply the fluorescence emission ratio η(c) by the fluorescence emission intensity I(c) at any point in space to calculate the fluorescence emission intensity of the device under test at that location, which is I(c)η(c).

[0055] Alternatively, the change in fluorescence emission intensity in the two-dimensional plane can be obtained based on the two-dimensional plane distribution of the first laser and the two-dimensional plane distribution of the second laser.

[0056] The process of obtaining the fluorescence emission intensity I(c) at any point in space includes:

[0057] Under the conditions of the first excitation power P1 and the first excitation wavelength A(λ) of the fluorescence acquisition system, the fluorescence emission intensity of a series of fluorescent molecule samples with different concentration gradients is detected, and the fluorescence emission intensity is obtained as a discrete sampling point with varying fluorescent molecule concentration. The fluorescence emission intensity is then fitted to the fluorescent molecule concentration to obtain the fluorescence emission intensity I(c) at any point in space.

[0058] The process of obtaining the change in two-dimensional plane fluorescence emission intensity based on the two-dimensional plane distribution of the first laser and the two-dimensional plane distribution of the second laser includes:

[0059] The fluorescence emission intensity in two-dimensional space is obtained from the fluorescence emission intensity I(c) at any point in space and is denoted as I1(c, x, y).

[0060] Let the first two-dimensional plane distribution of the fluorescence acquisition system be denoted as D1(x, y), and the second two-dimensional plane distribution of the laser from the device under test be denoted as D2(x, y); then the change in fluorescence emission intensity in the two-dimensional plane is: The digital phantom strength I2(c,x,y) = η(c,x,t)I1(c,x,y).

[0061] Example 1

[0062] This embodiment provides a digital phantom calculation case based on the above-described method for quantitative calculation of digital phantom intensity using illumination feedback.

[0063] Step S1: Set up the fluorescence acquisition system.

[0064] The excitation of the fluorescence acquisition system was characterized. In the fluorescence acquisition system, an upright fluorescence imaging setup was used, with the first grayscale camera 1 placed 60 mm away from the sample. An 850 nm filter was placed in front of the camera to detect the fluorescence emission intensity of indocyanine green. Simultaneously, a spectrometer was used to detect the fluorescence spectrum curve, excitation efficiency curve, and excitation wavelength. Figure 3 As shown, the second grayscale camera 2 is placed near the first grayscale camera 1, and the spatial position of the image is calibrated to the first grayscale camera 1; a 780nm color filter is added in front of the second grayscale camera 2 to detect the first laser two-dimensional planar distribution of the excitation light; the first excitation power of the excitation light is detected using a 780nm calibrated optical power meter. Figure 4 The diagram shows that the first excitation power of the fluorescence acquisition system is 5.69W, and also provides a schematic diagram of the first laser two-dimensional planar distribution and the first excitation wavelength.

[0065] Step S2: Prepare a series of gradient indocyanine green solutions.

[0066] Indocyanine green aqueous solutions were prepared at concentrations of 1, 1 / 2, 1 / 4, 1 / 8, and 1 / 16 mg / ml, with each concentration spaced two times apart. A fluorescence acquisition system was used to detect the fluorescence in the samples. Under excitation by a 780 nm laser, a first grayscale camera recorded discrete sample points showing the change in fluorescence emission intensity with concentration. These discrete sample points were then fitted with a function to obtain the functional relationship between the fluorescence emission intensity of indocyanine green and its concentration, as shown below. Figure 7 As shown.

[0067] A fluorescence acquisition system was used to detect the fluorescence spectrum curves of the samples. Under excitation by a 780nm laser, the spectrometer recorded discrete sample points showing the change in fluorescence emission spectra with concentration. Then, the center wavelength of each discrete sample point was extracted, and a function was fitted to the center wavelength to obtain the functional relationship between the fluorescence emission spectrum of indocyanine green and concentration, as shown below. Figure 6 As shown.

[0068] The excitation efficiency curve of the sample was detected using a fluorescence acquisition system. Under excitation light of 700 nm to 800 nm, the spectrometer recorded discrete sample points showing the change in fluorescence emission spectrum with excitation wavelength. Then, the discrete sample points were integrated, and the integration result was normalized to obtain the functional relationship between indocyanine green and excitation wavelength, as shown below. Figure 5 As shown.

[0069] Step S3, characterize the excitation of the device under test, such as... Figure 8As shown. A white diffuse reflection cardboard is placed on the sample plane, the excitation illumination of the device under test is turned on, the second grayscale camera 2 is used to detect the second laser two-dimensional planar distribution of the excitation light, the second excitation wavelength of the excitation light is detected using a spectrometer, and the second excitation power of the excitation light is detected using an optical power meter.

[0070] Step S4: Calculate the data obtained from the above detection, specifically as follows:

[0071] Set the fluorescence concentration value c, and obtain the fluorescence emission spectrum curve of the device under test from the relationship between the fluorescence center wavelength and the fluorescence concentration W(c).

[0072] Based on the second excitation wavelength B(λ) of the device under test (DUT), the excitation efficiency curve E(λ), the first excitation wavelength A(λ) of the fluorescence acquisition system, the second excitation power P2 of the DUT, and the first excitation power P1 of the fluorescence acquisition system, the fluorescence emissivity η(c) is calculated. The fluorescence emissivity η(c) is then multiplied by the fluorescence emission intensity I(c) at any point in space to obtain the fluorescence emission intensity of the DUT at that location, which is I(c)η(c). Figure 9 As shown.

[0073] Accordingly, this application also provides an electronic device, comprising: one or more processors; a memory for storing one or more programs; and, when the one or more programs are executed by the one or more processors, causing the one or more processors to implement the digital phantom intensity calculation method using illumination light feedback as described above. Figure 10 The diagram shown is a hardware structure diagram of any device with data processing capabilities for calculating the intensity of a digital phantom body using illumination feedback, as provided in an embodiment of the present invention. Except for... Figure 10 In addition to the processor, memory, and network interface shown, any data processing device in the embodiment may also include other hardware depending on the actual function of the data processing device, which will not be described in detail here.

[0074] Accordingly, this application also provides a computer-readable storage medium storing computer instructions that, when executed by a processor, implement the digital phantom intensity calculation method using illumination light feedback as described above. The computer-readable storage medium can be an internal storage unit of any data-processing device as described in any of the foregoing embodiments, such as a hard disk or memory. The computer-readable storage medium can also be an external storage device, such as a plug-in hard disk, smart media card (SMC), SD card, flash card, etc., equipped on the device. Furthermore, the computer-readable storage medium can include both internal storage units of any data-processing device and external storage devices. The computer-readable storage medium is used to store the computer program and other programs and data required by the data-processing device, and can also be used to temporarily store data that has been output or will be output.

[0075] The above embodiments are only used to illustrate the design concept and features of the present invention, and their purpose is to enable those skilled in the art to understand the content of the present invention and implement it accordingly. The protection scope of the present invention is not limited to the above embodiments. Therefore, all equivalent changes or modifications made based on the principles and design ideas disclosed in the present invention are within the protection scope of the present invention.

Claims

1. A digital phantom intensity calculation method using illumination light feedback, characterized by, The method specifically comprises the following steps: Acquiring a first excitation power of a fluorescence acquisition system and a first excitation wavelength ; Obtaining the fluorescence molecule concentration, fluorescence spectrum curve, and fluorescence emission intensity of the fluorescence molecule sample and fluorescence emission intensity acquiring a second excitation power of the device under test and a second excitation wavelength ; Setting the fluorescence concentration value c, from the relationship between the fluorescence center wavelength and the fluorescence concentration Obtaining the fluorescence emission spectrum curve corresponding to the device under test; a second excitation wavelength of the device under test , an excitation efficiency curve , a first excitation wavelength of the fluorescence collection system , a second excitation power of the device under test , and a first excitation power of the fluorescence collection system calculating a fluorescence emission ratio ; comprising: A first excitation wavelength of a fluorescence collection system is calculated with an integral of an excitation efficiency curve to obtain a corresponding relative fluorescence excitation amount of the fluorescence collection system; Computing a second excitation wavelength of the device under test with an integral of the excitation efficiency curve to obtain a corresponding relative fluorescence excitation of the device under test; The ratio of the relative fluorescence excitation of the fluorescence acquisition system to the relative fluorescence excitation of the to-be-tested equipment is taken as the relative fluorescence emission ratio; calculating a ratio of the second excitation power of the device under test to the first excitation power of the fluorescence collection system and multiplying the ratio by the relative fluorescence emission ratio to obtain a fluorescence emission ratio ​​ fluorescence emission ratio Fluorescence emission intensity at any point in space The fluorescence emission intensity of the device under test at that point is obtained by multiplication.

2. The digital phantom intensity calculation method using illumination light feedback according to claim 1, wherein, Relationship between fluorescence center wavelength and fluorescence concentration The acquisition process includes: A series of concentration gradient fluorescent molecule samples are configured, a fluorescence spectrum curve is detected, and a fluorescence center wavelength of the fluorescence spectrum curve is extracted; The fluorescence center wavelength and the fluorescence concentration are fitted to obtain the relationship between the fluorescence center wavelength and the fluorescence concentration .

3. The method of claim 1, wherein the method further comprises: Firing efficiency curve The acquisition process comprises: The fluorescence acquisition system acquires fluorescence data at a series of first excitation wavelengths The fluorescence acquisition system acquires fluorescence data at a series of first excitation wavelengths The fluorescence acquisition system acquires fluorescence data at a series of first excitation wavelengths The fluorescence acquisition system acquires fluorescence data at a series of first excitation wavelengths .

4. The method of claim 1, wherein the method is characterized by, fluorescence emission intensity at any point in space The acquisition process comprises: Under the condition of the first excitation power , the first excitation wavelength of the fluorescence acquisition system, the fluorescence emission intensity of a series of concentration gradient fluorescent molecule samples is detected, and the fluorescence molecule concentration and the fluorescence emission intensity are fitted, so as to obtain the fluorescence emission intensity at an arbitrary point in space .

5. The method of claim 1, wherein the method further comprises: The method further comprises: According to the obtained fluorescence emission intensity at any point in space The fluorescence emission intensity in two-dimensional space is denoted as ; Acquiring a first two-dimensional laser plane distribution of a fluorescence acquisition system , acquiring a second two-dimensional laser plane distribution of the device under test ; According to the first laser two-dimensional plane distribution , the second laser two-dimensional plane distribution Obtaining two-dimensional plane fluorescence emission intensity changes , thereby obtaining digital phantom intensity .

6. The digital phantom intensity calculation method using illumination light feedback according to claim 1 or 5, characterized in that, Acquiring a first excitation power of a fluorescence acquisition system and a first excitation wavelength The process comprises: Placing and assembling a black-and-white camera and a lens group, and detecting a first two-dimensional plane distribution of the illumination light through a two-dimensional plane detector within an imageable distance range; at a fixed central position of the imaging plane, a first excitation power of the illumination light is recorded by a power meter ; At a fixed central position of the imaging plane, a first excitation spectrum of the illumination light is recorded by the spectrometer fiber probe, resulting in a first excitation wavelength ; acquiring a second excitation power of the device under test and a second excitation wavelength of the device under test, the process comprising: Placing and assembling a black-and-white camera and a lens group, and detecting a second two-dimensional plane distribution of the illumination light through a two-dimensional plane detector within an imageable distance range; at a fixed central position of the imaging plane, the second excitation power of the illumination light is recorded by a power meter ; At a fixed central position of the imaging plane, a second excitation spectrum of the illumination light is recorded by the spectrometer fiber probe, resulting in a second excitation wavelength .

7. An electronic device comprising a memory and a processor, characterized in that The memory is coupled with the processor; wherein the memory is configured to store program data, and the processor is configured to execute the program data to implement the digital dummy intensity calculation method with illumination light feedback according to any one of claims 1-6.

8. A computer readable storage medium having stored thereon a computer program, characterized in that, The program is executed by the processor to implement the digital dummy intensity calculation method with illumination light feedback according to any one of claims 1-6.

9. A computer program product comprising computer programs / instructions, characterized in that, The computer program / instruction is executed by the processor to implement the digital dummy intensity calculation method with illumination light feedback according to any one of claims 1-6.

Images