CN-122025031-A - Biological scouring monitoring method, device, computer equipment and storage medium

CN122025031ACN 122025031 ACN122025031 ACN 122025031ACN-122025031-A

Abstract

The method is used for determining time activity curves of different types of physiological areas by collecting dynamic PET data of radioactive markers in a target object, respectively carrying out sectional fitting on the time activity curves of the different types of physiological areas so as to determine kinetic parameters of the different types of physiological areas, and respectively carrying out linear fitting on the kinetic parameters of the different types of physiological areas according to preset fitting parameters, thereby accurately obtaining the biological flushing rate of the different types of physiological areas, improving the accuracy of subsequent biological flushing effect modeling and range verification, and being suitable for different tissue types and individual differences of users.

Inventors

LIU YANG
QIU AO
XIE QINGGUO

Assignees

苏州瑞派宁科技有限公司

Dates

Publication Date: 20260512
Application Date: 20251226

Claims (20)

1. A method of biological flushing monitoring, the method comprising: Collecting dynamic PET data of a radioactive marker in a target object, and determining time activity curves of different types of physiological areas; Respectively carrying out sectional fitting on the time activity curves of the physiological areas of all types to determine the kinetic parameters of the physiological areas of all types; And respectively carrying out linear fitting on the dynamic parameters of the physiological areas of each type according to preset fitting parameters to obtain the biological flushing rate of the physiological areas of each type, wherein the fitting parameters are used for representing the mapping relation between the dynamic parameters of the physiological areas of each type and the corresponding biological flushing rate.
2. The method of claim 1, wherein after obtaining the biological flush rate for each type of physiological region, the method further comprises: Determining a production rate of the radioactive label based on the monitored data of the proton beam irradiation; determining the concentration change rate of the radioactive marker according to the production rate of the radioactive marker and the biological scouring rate of various types of physiological areas.
3. The method of biological flushing monitoring of claim 2, wherein the determining the rate of change of the concentration of the radiolabel based on the rate of production of the radiolabel and the rate of biological flushing of each type of physiological region comprises: determining the concentration change rate of the radioactive marker based on a multi-chamber biological scour dynamics model according to the production rate of the radioactive marker and the biological scour rate of various types of physiological areas.
4. A method of monitoring biological flushing according to claim 3 wherein the multi-chamber biological flushing dynamics model is implemented using the formula: wherein R represents the production rate of the radiolabel, k' 2,i represents the biological flushing rate of the i-th type of physiological region, C i (t) represents the activity concentration of the radiolabel at time t in the i-th type of physiological region, The rate of change of the concentration of the radiolabel at time t is indicated.
5. The method of biological flushing monitoring of claim 4, further comprising: determining the nuclides generated by proton beam irradiation and the initial activity and half-life of the corresponding nuclides for any type of physiological region; Correcting the activity concentration of the radioactive marker in the physiological region according to the initial activity and half-life of the nuclide, and determining the corrected activity concentration in the physiological region.
6. The method of claim 1, wherein before the linear fitting of the kinetic parameters of each type of physiological region according to the preset fitting parameters, the method further comprises: And establishing a mapping relation between the sample kinetic parameters of different types of physiological areas and the corresponding sample biological flushing rates based on experimental data.
7. The method of claim 6, wherein establishing a mapping relationship between the sample dynamics parameters of different types of physiological regions and the corresponding sample biological flush rates based on the experimental data comprises: acquiring a sample biological flushing rate of a sample tracer based on a gold standard mode for each type of physiological region; Synchronously acquiring experimental PET data of the sample tracer, and fitting to obtain sample kinetic parameters of a corresponding type of physiological region; Performing linear fitting according to the sample biological flushing rate and the sample kinetic parameters corresponding to each type of physiological region to determine fitting parameters of the corresponding physiological region; And establishing a mapping relation between the sample kinetic parameters of the corresponding type physiological region and the sample biological flushing rate according to the fitting parameters.
8. The method for monitoring biological flushing according to claim 1, wherein the step of linearly fitting the kinetic parameters of each type of physiological region according to preset fitting parameters to obtain the biological flushing rate of each type of physiological region includes: Aiming at any type of physiological area, determining target fitting parameters of the corresponding type of physiological area in preset fitting parameters; and performing linear fitting on the dynamic parameters and the target fitting parameters of the physiological region based on the linear fitting function to obtain the biological flushing rate of the corresponding physiological region.
9. The method of claim 8, wherein the linear fitting function is implemented using the following equation: k′ 2,i ＝a i k 2,i +b i , Where k 2,i represents the kinetic parameters of the i-th type of physiological region, a i and b i represent the target fitting parameters of the i-th type of physiological region, and k' 2,i is the biological flush rate of the i-th type of physiological region obtained by fitting.
10. The method of claim 1, wherein the step of fitting the time activity curves of the physiological regions of each type separately to determine the kinetic parameters of the physiological regions of each type comprises: and respectively carrying out sectional fitting on the time activity curves of the physiological areas of all types by adopting a distributed dynamics model so as to determine the dynamics parameters of the physiological areas of all types.
11. The method of claim 10, wherein the distributed kinetic model is implemented using the formula: wherein F i (T) represents the blood flow rate of the ith physiological region, T ci represents the average vascular transit time of the ith physiological region, K 1,i represents the blood transport rate of the ith physiological region to the tissue, K 2,i represents the clearance rate of the tissue of the ith physiological region, and R i (T) represents the activity concentration of the ith physiological region.
12. The method of biological flushing monitoring of any one of claims 1 to 11, wherein determining the time activity profile for different types of physiological regions includes: identifying regions with different pathological features based on the dynamic PET data to obtain different types of physiological regions; Determining, for each type of physiological region, an activity concentration sequence over time of the corresponding type of physiological region radiolabel based on the dynamic PET data; And determining the time activity curve of the corresponding type of physiological region according to the activity concentration sequence of the radioactive marker of the corresponding type of physiological region along with time.
13. The method of any one of claims 1 to 11, wherein the different types of physiological regions are different tissue regions divided based on pathological features.
14. The method of claim 13, wherein the pathological characteristic comprises at least one of a focal activity characteristic, a vascular status characteristic, and a tissue type characteristic.
15. A biological flush monitoring device, the device comprising: the data acquisition module is configured to acquire dynamic PET data of the radioactive marker in the target object and determine time activity curves of different types of physiological areas; The dynamic parameter determining module is configured to perform segment fitting on the time activity curves of the various types of physiological areas respectively so as to determine dynamic parameters of the various types of physiological areas; The biological flushing rate determining module is configured to respectively perform linear fitting on the dynamic parameters of each type of physiological area according to preset fitting parameters to obtain the biological flushing rate of each type of physiological area, wherein the fitting parameters are used for representing the mapping relation between the dynamic parameters of each type of physiological area and the corresponding biological flushing rate.
16. The biological flush monitoring device of claim 15, wherein the device further comprises: a data monitoring module configured to determine a production rate of the radiolabel based on the monitored data of the proton beam irradiation; A concentration change rate determination module configured to determine a concentration change rate of the radiolabel based on a production rate of the radiolabel and a biological flush rate of each type of physiological region.
17. The biological flush monitoring device according to claim 16, wherein the concentration change rate determination module is further configured to determine the concentration change rate of the radioactive marker based on a multi-chamber biological flush kinetic model according to a production rate of the radioactive marker and a biological flush rate of each type of physiological region.
18. The biological flush monitoring device of claim 17, wherein, The multi-chamber biological scouring kinetic model is realized by adopting the following formula: wherein R represents the production rate of the radiolabel, k' 2,i represents the biological flushing rate of the i-th type of physiological region, C i (t) represents the activity concentration of the radiolabel at time t in the i-th type of physiological region, The rate of change of the concentration of the radiolabel at time t is indicated.
19. The biological flush monitoring device of claim 18, further comprising a correction module configured to determine an initial activity and half-life of a nuclear species and a corresponding nuclear species generated by proton beam irradiation for any type of physiological region, and to correct an activity concentration of a radioactive marker of the physiological region based on the initial activity and half-life of the nuclear species, and to determine the corrected activity concentration of the physiological region.
20. The biological flush monitoring device according to claim 15, further comprising a mapping relationship establishing module configured to establish a mapping relationship between sample kinetic parameters of different types of physiological regions and corresponding sample biological flush rates based on experimental data.

Description

Biological scouring monitoring method, device, computer equipment and storage medium Technical Field The present disclosure relates to the field of medical image data processing, and in particular, to a biological flushing monitoring method, apparatus, computer device, computer readable storage medium, and computer program product. Background The particle therapy such as proton therapy and carbon ion can intensively release the maximum dose of radiation at the tail end of a tumor target area by virtue of the unique physical characteristics of Bragg peak, so that the damage to surrounding normal tissues is remarkably reduced, and the particle therapy has become an important development direction of accurate radiotherapy. However, the clinical accuracy of this technique is highly dependent on precise control of the particle range, and range uncertainty (mainly due to scattering of the particle beam in human tissue, individual physiological tissue differences, and dynamic biological effect changes, etc.) remains a critical issue to be addressed currently. To verify the consistency of the actual throw of the particle beam with the plan, positron Emission Tomography (PET) is widely used for range verification after proton treatment, indirectly reflecting the dose deposition location by detecting the positron emitter distribution produced by the proton nuclear reaction. However, the biological scour effect (biological washout effect) severely limits the accuracy of PET imaging. The effect means that the electron emitter migrates or is cleared in the body due to factors such as blood perfusion, metabolic clearance, tissue pathological states and the like, so that the spatial distribution of the electron emitter deviates from an original dose deposition area, and particularly, the electron emitter shows obvious differences in tumor microenvironment with rich blood vessels or necrosis, thereby influencing the distribution of the electron emitter in PET imaging. The traditional modeling method of biological scouring effect in proton treatment mainly comprises the following steps of 1) simplifying a scouring process by adopting a uniform attenuation constant, ignoring tissue heterogeneity, 2) improving resolution by introducing fast, medium and slow scouring components into a multicomponent model, and still lacking adaptation to individual pathophysiological characteristics, and 3) learning complex nonlinear relations from dynamic PET data by using an artificial intelligence-based modeling method, wherein the modeling method is sensitive to data quality and has insufficient stability under clinical noise and data deletion conditions. Thus, there is a need for an accurate determination method that can be adapted to clinical practice and that fuses physiological information of an individual. Disclosure of Invention The present disclosure provides a bio-flush monitoring method, apparatus, computer device, computer-readable storage medium, and computer program product to solve at least the problems in the related art. The technical scheme of the present disclosure is as follows: according to a first aspect of embodiments of the present disclosure, there is provided a biological flush monitoring method comprising: Collecting dynamic PET data of a radioactive marker in a target object, and determining time activity curves of different types of physiological areas; Respectively carrying out sectional fitting on the time activity curves of the physiological areas of all types to determine the kinetic parameters of the physiological areas of all types; And respectively carrying out linear fitting on the dynamic parameters of the physiological areas of each type according to preset fitting parameters to obtain the biological flushing rate of the physiological areas of each type, wherein the fitting parameters are used for representing the mapping relation between the dynamic parameters of the physiological areas of each type and the corresponding biological flushing rate. In one embodiment, after the biological flushing rate of each type of physiological region is obtained, the method further comprises determining a production rate of the radioactive marker based on the monitoring data of the proton beam irradiation, and determining a concentration change rate of the radioactive marker according to the production rate of the radioactive marker and the biological flushing rate of each type of physiological region. In one embodiment, the determining the rate of change of the concentration of the radioactive label based on the rate of production of the radioactive label and the rate of biological flushing of each type of physiological region comprises determining the rate of change of the concentration of the radioactive label based on a multi-chamber biological flushing kinetic model based on the rate of production of the radioactive label and the rate of biological flushing of each type of physiological region. In one embodiment