CN110269590A - 药代动力学参数确定方法、装置、计算机设备和存储介质 - Google Patents
药代动力学参数确定方法、装置、计算机设备和存储介质 Download PDFInfo
- Publication number
- CN110269590A CN110269590A CN201910537164.0A CN201910537164A CN110269590A CN 110269590 A CN110269590 A CN 110269590A CN 201910537164 A CN201910537164 A CN 201910537164A CN 110269590 A CN110269590 A CN 110269590A
- Authority
- CN
- China
- Prior art keywords
- drug
- activity time
- tissue
- image
- frame
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 239000003814 drug Substances 0.000 claims abstract description 206
- 229940079593 drug Drugs 0.000 claims abstract description 156
- 230000000694 effects Effects 0.000 claims abstract description 123
- 238000005070 sampling Methods 0.000 claims abstract description 48
- 238000012905 input function Methods 0.000 claims abstract description 47
- 210000002381 plasma Anatomy 0.000 claims description 89
- 210000001519 tissue Anatomy 0.000 claims description 78
- 230000006870 function Effects 0.000 claims description 38
- 238000004590 computer program Methods 0.000 claims description 28
- 210000000709 aorta Anatomy 0.000 claims description 17
- 238000002591 computed tomography Methods 0.000 claims description 11
- 238000003384 imaging method Methods 0.000 claims description 9
- 230000008520 organization Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 13
- 239000000700 radioactive tracer Substances 0.000 description 16
- 210000004369 blood Anatomy 0.000 description 11
- 239000008280 blood Substances 0.000 description 11
- 230000008859 change Effects 0.000 description 10
- 239000007924 injection Substances 0.000 description 6
- 238000002347 injection Methods 0.000 description 6
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 5
- 239000008103 glucose Substances 0.000 description 5
- 230000026731 phosphorylation Effects 0.000 description 5
- 238000006366 phosphorylation reaction Methods 0.000 description 5
- 239000000243 solution Substances 0.000 description 5
- 229920002472 Starch Polymers 0.000 description 4
- 239000002872 contrast media Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000002708 enhancing effect Effects 0.000 description 4
- 230000008289 pathophysiological mechanism Effects 0.000 description 4
- 235000019698 starch Nutrition 0.000 description 4
- 239000008107 starch Substances 0.000 description 4
- 238000012800 visualization Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000003902 lesion Effects 0.000 description 3
- 239000002002 slurry Substances 0.000 description 3
- 238000010408 sweeping Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 2
- 238000002583 angiography Methods 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- 230000002427 irreversible effect Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 230000004060 metabolic process Effects 0.000 description 2
- 230000037323 metabolic rate Effects 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000000163 radioactive labelling Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- PXGOKWXKJXAPGV-UHFFFAOYSA-N Fluorine Chemical compound FF PXGOKWXKJXAPGV-UHFFFAOYSA-N 0.000 description 1
- XUHXFSYUBXNTHU-UHFFFAOYSA-N Iotrolan Chemical compound IC=1C(C(=O)NC(CO)C(O)CO)=C(I)C(C(=O)NC(CO)C(O)CO)=C(I)C=1N(C)C(=O)CC(=O)N(C)C1=C(I)C(C(=O)NC(CO)C(O)CO)=C(I)C(C(=O)NC(CO)C(O)CO)=C1I XUHXFSYUBXNTHU-UHFFFAOYSA-N 0.000 description 1
- 239000002616 MRI contrast agent Substances 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 235000009754 Vitis X bourquina Nutrition 0.000 description 1
- 235000012333 Vitis X labruscana Nutrition 0.000 description 1
- 240000006365 Vitis vinifera Species 0.000 description 1
- 235000014787 Vitis vinifera Nutrition 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 235000014113 dietary fatty acids Nutrition 0.000 description 1
- 230000036267 drug metabolism Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 229930195729 fatty acid Natural products 0.000 description 1
- 239000000194 fatty acid Substances 0.000 description 1
- 150000004665 fatty acids Chemical class 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 239000011737 fluorine Substances 0.000 description 1
- 125000001153 fluoro group Chemical group F* 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000010253 intravenous injection Methods 0.000 description 1
- 229960003182 iotrolan Drugs 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- MIKKOBKEXMRYFQ-WZTVWXICSA-N meglumine amidotrizoate Chemical compound C[NH2+]C[C@H](O)[C@@H](O)[C@H](O)[C@H](O)CO.CC(=O)NC1=C(I)C(NC(C)=O)=C(I)C(C([O-])=O)=C1I MIKKOBKEXMRYFQ-WZTVWXICSA-N 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- VIKNJXKGJWUCNN-XGXHKTLJSA-N norethisterone Chemical compound O=C1CC[C@@H]2[C@H]3CC[C@](C)([C@](CC4)(O)C#C)[C@@H]4[C@@H]3CCC2=C1 VIKNJXKGJWUCNN-XGXHKTLJSA-N 0.000 description 1
- 238000009206 nuclear medicine Methods 0.000 description 1
- 102000039446 nucleic acids Human genes 0.000 description 1
- 108020004707 nucleic acids Proteins 0.000 description 1
- 150000007523 nucleic acids Chemical class 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000008288 physiological mechanism Effects 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 210000004881 tumor cell Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4836—Diagnosis combined with treatment in closed-loop systems or methods
- A61B5/4839—Diagnosis combined with treatment in closed-loop systems or methods combined with drug delivery
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Radiology & Medical Imaging (AREA)
- Pharmacology & Pharmacy (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medicinal Chemistry (AREA)
- Vascular Medicine (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Nuclear Medicine (AREA)
Abstract
本申请涉及一种药代动力学参数确定方法、装置、计算机设备和存储介质。所述方法包括：获取组织药物活度时间曲线以及血浆药物活度时间采样点，将血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数，并根据组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值，再根据得到的至少一个参数值获得药物代谢动力学参数。采用本方法能够简化确定药物代谢动力学参数的流程，使提高效率。
Description
技术领域
本申请涉及医学成像技术领域，特别是涉及一种药代动力学参数确定方法、装置、计算机设备和存储介质。
背景技术
药代动力学(Pharmackinetices)简称药动学，主要研究药物进入机体的动态变化，包括药物在机体内的吸收、分布、生化转换(或称代谢)及排泄的过程。血液中药物浓度自始至终都处于动态变化之中，且药物在机体内处置过程较为复杂，受到体内外因素的影响。药物在作用部位的浓度受机体内的影响从而产生的动态变化。为了揭示药物在机体内的动态变化规律性，常常要借助数据的方法来阐明体内药物浓度随时间而变化的规律性，根据机体内药物浓度和时间的数据，建立一定的数学模型，再求得相应的药动学参数，通过这些参数来描述药物机体内的动态变化规律性。
然而，目前药代动力学参数确定方法，存在计算复杂耗时，效率低下等问题。
发明内容
基于此，有必要针对上述技术问题，提供一种的药代动力学参数确定方法、装置、计算机设备和存储介质。
一种药代动力学参数确定方法，所述方法包括：
获取组织药物活度时间曲线以及血浆药物活度时间采样点；
将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数；
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值；
根据至少一个所述参数值获得药物代谢动力学参数。
在其中一个实施例中，所述获取组织药物活度时间曲线包括：
获取多帧扫描图像；
在各帧所述扫描图像中查找组织相关位置的扫描图像，并将组织相关位置的扫描图像作为第一图像；
提取各帧所述第一图像像素点的像素值，并根据所述像素点的像素值计算各帧组织平均像素值；
将各帧所述组织平均像素值根据扫描时间进行排序，得到所述组织药物活度时间曲线。
在其中一个实施例中，所述获取血浆药物活度时间采样点包括：
获取多帧扫描图像；
在各帧所述扫描图像中查找血浆相关的位置的扫描图像，并将血浆相关位置的扫描图像作为第二图像；
提取各帧所述第二图像像素点的像素值，并根据所述像素点的像素值计算各帧血浆平均像素值；
将得到各帧主动脉位置扫描图像的主动脉平均像素值；
将各帧所述血浆平均像素值根据扫描时间进行排序，得到所述血浆药物活度时间采样点。
在其中一个实施例中，所述根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到多个药代动力模型的参数值包括：
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型建立药代动力等式；
对所述药代动力等式进行求解，得到多个药代动力模型的参数值。
在其中一个实施例中，所述获取组织药物活度时间曲线以及血浆药物活度时间采样点之后包括：
根据所述血浆药物活度时间采样点，选择相应的预设函数；
根据药物类型以及组织类型，选择相应的预设药代动力模型。
在其中一个实施例中，所述预设函数包括三指数函数。
在其中一个实施例中，多帧所述扫描图像由动态图像扫描设备扫描获取；
所述动态图像设备包括：正电子发射型计算机断层显像、磁共振成像设备以及电子计算机断层扫描设备。
一种药代动力学参数确定装置，所述装置包括：
获取模块，用于获取组织药物活度时间曲线以及血浆药物活度时间采样点；
拟合模块，用于将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数；
参数值得到模块，用于根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到多个药代动力模型的参数值；
药物代谢动力学参数获取模块，用于根据多个所述参数值获得药物代谢动力学参数。
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现以下步骤:
获取组织药物活度时间曲线以及血浆药物活度时间采样点；
将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数；
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值；
根据至少一个所述参数值获得药物代谢动力学参数。
一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现以下步骤:
获取组织药物活度时间曲线以及血浆药物活度时间采样点；
将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数；
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值；
根据至少一个所述参数值获得药物代谢动力学参数。
上述药代动力学参数确定法、装置、计算机设备和存储介质，通过获取组织药物活度时间曲线以及血浆药物活度时间采样点，将血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数，并根据组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值，再根据得到的至少一个参数值获得药物代谢动力学参数。这样通过利用预设函数拟合的方式优化药代动力学的输入函数，以简化得到药物代谢动力学参数的流程，使提高效率。
附图说明
图1为一个实施例中2Ti房室模型示意图图；
图2为一个实施例中药代动力学参数确定方法的流程示意图；
图3为一个实施例中药代动力学参数确定装置的结构框图；
图4为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
在医学成像技术领域中，在利用示踪剂造影成像时，可根据注入机体内的示踪剂结合药代动力学，以研究机体的病理生理机制，并且使其可视化。首先，将示踪剂注入机体内，之后跟踪示踪剂的运动，再通过药代动力学分析计算出药代动力学各参数，基于个参数得到失踪剂在机体内的运动状态以研究机体的病理生理机制。所述药代动力学分析包括药代动力学模型、输入函数(机体主动脉中血浆药物活度时间曲线)以及组织中药物活度时间曲线。通过这三个部分，可以计算出药代动力学模型中，具有各种生理学意义的临床参数，例如葡萄糖代谢率、脑氧摄取分数以及心肌血流量等。
其中，药代动力学模型为房室模型，而房室模型是基于特殊类型的数学模型。在药动学中，为了定量地描述药物体内过程的动态变化规律性，常常要借助数学的原理和方法来系统地阐明，机体内药量随时间变化的规律。房室模型理论从速度论的角度出发，建立一个数据模型来模拟机体，它将整个机体视为一个系统，并将该系统按动力学特性划分为若干个房室。因此将机体看成是由若干个房室组成的一个完整的系统，称之为房室模型。
将输入函数以数组的形式输入房室模型，得到组织中药物活度时间曲线，并建立等式。通过求解可得到房室模型中的各参数值，再利用各参数值进一步计算出药代动力学的参数值。
通常正电子发射型计算机断层显像技术结合药代动力学，利用正电子发射型计算机断层显像设备采用反射性标记生物分子作为示踪剂以研究生物体中的病理的生理机制，并且使其可视化。其中正电子发射型计算机断层显像(Positron Emission Tomography，PET)，是核医学领域比较先进的临床检查影像技术。是将某种物质，一般是生物生命代谢中必须的物质，如：葡萄糖、蛋白质、核酸、脂肪酸，标记上短寿命的放射性核素(如18F，11C等)，注入人体后，放射性核素在衰变过程中释放出正电子，一个正电子在行进十分之几毫米到几毫米后遇到一个电子后发生湮灭，从而产生方向相反的一对能量为511KeV的光子。这对光子，通过高度灵敏的照相机捕捉，并经计算机进行散射和随机信息的校正。经过对不同的正电子进行相同的分析处理，我们可以得到在生物体内聚集情况的三维图像，从而达到诊断的目的。
如图1所示，以常见的药代动力学模型，2Ti房室模型(两组织不可逆房室模型)以及在PET扫描成像中采用放射性标记生物分子(FDG，氟代葡萄糖)作为示踪剂为例进一步阐述求解房室模型中各参数值过程。
模块100为中央室，中央室为主动脉血浆，其输入函数为血浆中的FDG浓度CP(t)(主动脉血浆中的示踪剂活度时间曲线)；模块102以及模块104为同一外周室，模块102为目标组织内的FDG浓度Ce(t)，模块104为该组织内由部分FDG磷酸化转化为FDG-6-P(磷酸化氟代葡萄糖)浓度Cm(t)，则该组织中示踪剂活度时间曲线为CT(t),CT(t)为Ce(t)与Cm(t)之和。将2Ti房室模型通过数学公式表达为以下方程：
其中
其中药代动力学模型中的K1、k2、k3以及k4为需要求解的参数值。在2Ti房室模型中，参数值K1为药物(FDG)从血浆到组织中的渗透速率；参数值k2描述的是药物(FDG)从组织中的到血浆中的反向渗透速率；参数值k3为FDG磷酸化的过程，或者其他药物特异性结合到细胞受体上的转化率；参数值k4为FDG逆磷酸化的转化率，一般可以认为肿瘤细胞的不可逆磷酸化(k4＝0)。在其他的药代动力学模型中，参数不一定都有明确的临床含义，但是临床用到的肯定是模型中都具有明显临床意义的参数值或参数值组合。
在利用传统方法在求解各参数值过程中，会将主动脉血浆中的示踪剂活度时间曲线进行预处理，拟合出可以使用的输入函数CP(t)，此时输入函数CP(t)为数组形式。之后再将输入函数插值到组织中示踪剂活度时间曲线对应采集时刻的位置，然后再求解方程(1)中的各参数值。在求解过程中，由于方程(1)中有卷积(积分)项，把输入函数相邻点之间视为直线段，并进行多个直线段的拼接，从而求解出模型中各参数值。但由于模型求解中通常包含有积分项，或者拟合后的输入函数与动态造影图像时间不匹配的情况。此时需要对数组形式的输入函数进行分段累积积分，插值补点，这样的求解方式复杂耗时，效率低下。并且在模型中出现时间延迟校正参数等于时间相关的修正项参数需要计算时，输入函数分段累积积分和插值补点带来的求解复杂缺点进一步放大。
在一个实施例中，如图2所示，有必要针对上述技术问题，提供了一种药代动力学参数确定方法，所述方法包括以下步骤：
步骤202，获取组织药物活度时间曲线以及血浆药物活度时间采样点。
在本实施例中，药物通过静脉注射进入血管后，通过人体循环系统进入中央室，所述中央室为主动脉血浆。然后药物由中央室转移扩散到外周室中，也就是组织中。在这个过程中，主动脉血浆中的药物浓度逐渐降低，组织中药物浓度逐渐升高，直至主动脉血浆中的药物浓度以及组织中的药物浓度逐渐达到动态平衡。血浆药物活度时间采样点为：静脉注射以后，根据不同时间点采集主动脉中血浆的药物活度。而组织药物活度时间曲线为：静脉注射以后，组织中药物活度根据时间变化而获取的曲线。
在本实施例中，以PET设备扫描过程中采用放射性标记生物分子(FDG，氟代葡萄糖)作为示踪剂为例，获取的血浆药物活度时间采样点为主动脉血浆中FDG活度，而组织药物活度时间曲线为组织中FDG以及FDG-6-P的活度时间曲线。
其中，获取组织药物活度时间曲线包括：获取多帧扫描图像，在各帧所述扫描图像中查找组织相关位置的扫描图像，并将组织相关位置的扫描图像作为第一图像，提取各帧所述第一图像像素点的像素值，并根据所述像素点的像素值计算各帧组织平均像素值，将各帧所述组织平均像素值根据扫描时间进行排序，得到所述组织药物活度时间曲线。
在本实施例中，由PET设备获取多帧扫描图像，其中多帧扫描图像为PET设备对待测物体进行连续扫描后得到的原始扫描图像。在PET设备对人体进行扫描时，扫描不是瞬间结束的，而是会持续一段时间，在这段时间内人体各部位的示踪剂浓度随着时间变化。通过将PET设备获取的多帧扫描图像按时间排序将会得到动态的扫描图像，而动态的扫描图像中也会显示示踪剂浓度随时间的变化。
进一步的，在各帧扫描图像中查找与组织相关位置的扫描图像，此时组织为药代动力学模型中的周边室，根据感兴趣的部位不同，周边室可为器官组织，或者是血流量供应较少的组织。在每一帧中，将查找的与感兴趣组织相关的图像区域作为第一图像。
进一步的，在各帧第一图像中的提取像素点的像素值，并且根据各帧中的多个像素值计算平均像素值，得到各第一图像中的组织平均像素值。再将各帧第一图像中的组织平均像素值根据扫描时间进行排序，得到组织中的药物活度时间曲线。
在本实施例中，可以在各帧第一图像中提取一个像素点的像素值，作为各第一图像的组织平均像素值。也可以提取多个像素点的像素值，再根据各帧中的多个像素值计算平均像素值，得到各第一图像中的组织平均像素值。
其中，获取血浆药物活度时间采样点包括：获取多帧扫描图像，在各帧所述扫描图像中查找血浆相关的位置的扫描图像，并将血浆相关位置的扫描图像作为第二图像。提取各帧所述第二图像像素点的像素值，并根据所述像素点的像素值计算各帧血浆平均像素值。将得到各帧主动脉位置扫描图像的主动脉平均像素值。将各帧所述血浆平均像素值根据扫描时间进行排序，得到所述血浆药物活度时间采样点。
在本实施例中，可以在各帧第二图像中提取一个像素点的像素值，作为各第二图像的血浆平均像素值。也可以提取多个像素点的像素值，再根据各帧中的多个像素值计算平均像素值，得到各第二图像中的血浆平均像素值。
在本实施例中，由PET设备获取多帧扫描图像，其中多帧扫描图像为PET设备对待测物体进行连续扫描后得到的原始扫描图像。在PET设备对人体进行扫描时，扫描不是瞬间结束的，而是会持续一段时间，在这段时间内人体各部位的示踪剂浓度随着时间变化。通过将PET设备获取的多帧扫描图像按时间排序将会得到动态的扫描图像，而动态的扫描图像中也会显示示踪剂浓度随时间的变化。
进一步的，在各帧扫描图像中查找与主动脉相关位置的扫描图像，此时主动脉为药代动力学模型中的中央室。在每一帧中，将查找到的与主动脉相关的图像区域作为第二图像。
进一步的，在各帧第二图像中的提取像素点的像素值，并且根据各帧中的多个像素值计算平均像素值，得到各帧第二图像中的主动脉平均像素值。再将各帧第二图像中的主动脉平均像素值根据扫描时间进行排序，得到主动脉中的血浆药物活度时间采样点。
在其他实施例中，还可通过静动脉采血的方式获取血浆药物活度时间采样点。
步骤204，将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数。
在本实施例中，此时得到的血浆药物活度时间采样点为数组的形式，需要将数组形式的血浆药物活度时间采样点与预设函数进行拟合，进而得到血浆活度关于时间的模型化输入函数。
在步骤202之前还包括：根据所述血浆药物活度时间采样点，选择相应的预设函数。
具体的，根据血浆药物活度时间采样点的形状，选择形状相适合的预设函数。在PET动态扫描中，一般采用弹丸注射，这种方式产生的输入函数通常为一个快速上升尖峰伴随着持续衰减的形式。由此可以选用指数函数为预设函数。在本实施例中，可选择三指函数作为预设函数。在其他实施中，如果药物注射时比较缓慢，产生的输入函数的形状也相应发生变化，此时需要选用的拟合函数形状尽量与输入函数的形状相似，以使拟合效果更好。如果注射药物的方式让血浆药物活度时间曲线类似一条k<0的直线，则可以用y＝kx+b，作为预设函数，与血浆药物活度时间采样点进行拟合。通过将血浆药物活度时间采样点与预设的函数进行拟合解决了在使用传统方法时，将数组形式的采样点，时间序列不对应的问题，并且减少模型参数值求解的难度。
进一步的，将血浆药物活度时间采样点与三指数函数进行拟合后得到方程如下：
其中，w1、w2、w3、w4、w5以及w6是已知方程参数。
步骤206，根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值。
步骤206还包括：根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型建立药代动力等式，对所述药代动力等式进行求解，得到多个药代动力模型的参数值。
在步骤202之前还包括：根据药物类型以及组织类型，选择相应的预设药代动力模型。
具体的，根据注射到体内的药物不同，以及感兴趣的组织类型不同，可相应的选择预设药代动力模型。以注射的药物为FDG的肿瘤检测为例，则选择数学方程表达如方程(1)所示的2Ti模型。
在本实施例中，根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型建立药代动力等式具体为，将方程(2)带入方程(1)中的，得到方程(3)：
其中，由于FDG-6-P不可在逆转化为FDG,则将参数值k4＝0带入方程(3)，简化得到方程(4)：
通过求解方程(4)，可以确定2Ti模型中的参数值K1、k2、k3。
步骤208，根据至少一个所述参数值获得药物代谢动力学参数。
在本实施例中，根据步骤206中确定的至少一个参数值再计算得到药物代谢动力学参数。通过一个参数，或者多个参数的计算可获得葡萄糖代谢率等药物代谢动力学参数。
在其他实施例中，多帧所述扫描图像由动态图像扫描设备扫描获取，所述动态图像设备包括：正电子发射型计算机断层显像、磁共振成像设备以及电子计算机断层扫描设备。
在本实施例中，除了将本方法运用在PET动态扫描中，还可以运用当磁共振成像设备以及电子计算机断层扫描设备需要造影成像过程中。
其中，在利用磁共振成像设备做增强扫描成像以及对人体血管扫描成像过程中，需要往人体内注入造影剂。通过注射MRI造影剂，缩短组织在外磁场作用下的共振时间、增大对比信号的差异、提高成像对比度和清晰度的一类诊断试剂。它能有效改变生物体内组织中局部的水质子弛豫速率，缩短水分子中质子的弛豫时间，准确地检测出正常组织与患病部位之间的差异的一种检查方式。此时，可以运用上述方法，将注入人体的造影剂作为示踪剂以研究人体中的病理生理机制，并且使其过程可视化。
其中，在利用电子计算机断层扫描设备在做增强扫描或者造影扫描过程中，需要往人体感兴趣位置注入造影剂。在做增强扫描是，需要用高压注射器经静脉注入水溶性有机碘剂，如60％～76％泛影葡胺60ml后再行扫描的方法。血内碘浓度增高后，器官与病变内碘的浓度可产生差别，形成密度差，可能使病变显影更为清楚。在做造影扫描时，先作器官或结构的造影，然后再行扫描的方法。例如向脑池内注入碘曲仑8～10ml或注入空气4～6ml进行脑池造影再行扫描，称之为脑池造影CT扫描，可清楚显示脑池及其中的小肿瘤。此时，可以运用上述方法，将注入人体的造影剂作为示踪剂以研究人体中的病理生理机制，并且使其过程可视化。
上述药代动力学参数确定方法中，通过获取组织药物活度时间曲线以及血浆药物活度时间采样点，将血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数，并根据组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值，再根据得到的至少一个参数值获得药物代谢动力学参数。这样通过利用预设函数拟合的方式优化药代动力学的输入函数，以简化得到药物代谢动力学参数的流程，使提高效率。
应该理解的是，虽然图2流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中，如图3所示，提供了一种药代动力学参数确定装置，包括：获取模块300、拟合模块302、参数值得到模块304和药物代谢动力学参数获取模块306，其中：
获取模块300，用于获取组织药物活度时间曲线以及血浆药物活度时间采样点。
拟合模块302，用于将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数。
参数值得到模块304，用于根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到多个药代动力模型的参数值。
药物代谢动力学参数获取模块306，用于根据多个所述参数值获得药物代谢动力学参数。
关于药代动力学参数确定装置的具体限定可以参见上文中对于药代动力学参数确定方法的限定，在此不再赘述。上述药代动力学参数确定装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图4所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种药代动力学参数确定方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解，图4中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
获取组织药物活度时间曲线以及血浆药物活度时间采样点；
将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数；
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值；
根据至少一个所述参数值获得药物代谢动力学参数。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
获取多帧扫描图像；
在各帧所述扫描图像中查找组织相关位置的扫描图像，并将组织相关位置的扫描图像作为第一图像；
提取各帧所述第一图像像素点的像素值，并根据所述像素点的像素值计算各帧组织平均像素值；
将各帧所述组织平均像素值根据扫描时间进行排序，得到所述组织药物活度时间曲线。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
获取多帧扫描图像；
在各帧所述扫描图像中查找血浆相关的位置的扫描图像，并将血浆相关位置的扫描图像作为第二图像；
提取各帧所述第二图像像素点的像素值，并根据所述像素点的像素值计算各帧血浆平均像素值；
将得到各帧主动脉位置扫描图像的主动脉平均像素值；
将各帧所述血浆平均像素值根据扫描时间进行排序，得到所述血浆药物活度时间采样点。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型建立药代动力等式；
对所述药代动力等式进行求解，得到多个药代动力模型的参数值。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
根据所述血浆药物活度时间采样点，选择相应的预设函数；
根据药物类型以及组织类型，选择相应的预设药代动力模型。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
所述预设函数包括三指数函数。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
多帧所述扫描图像由动态图像扫描设备扫描获取；
所述动态图像设备包括：正电子发射型计算机断层显像、磁共振成像设备以及电子计算机断层扫描设备。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取组织药物活度时间曲线以及血浆药物活度时间采样点；
将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数；
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值；
根据至少一个所述参数值获得药物代谢动力学参数。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
获取多帧扫描图像；
在各帧所述扫描图像中查找组织相关位置的扫描图像，并将组织相关位置的扫描图像作为第一图像；
提取各帧所述第一图像像素点的像素值，并根据所述像素点的像素值计算各帧组织平均像素值；
将各帧所述组织平均像素值根据扫描时间进行排序，得到所述组织药物活度时间曲线。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
获取多帧扫描图像；
在各帧所述扫描图像中查找血浆相关的位置的扫描图像，并将血浆相关位置的扫描图像作为第二图像；
提取各帧所述第二图像像素点的像素值，并根据所述像素点的像素值计算各帧血浆平均像素值；
将得到各帧主动脉位置扫描图像的主动脉平均像素值；
将各帧所述血浆平均像素值根据扫描时间进行排序，得到所述血浆药物活度时间采样点。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型建立药代动力等式；
对所述药代动力等式进行求解，得到多个药代动力模型的参数值。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
根据所述血浆药物活度时间采样点，选择相应的预设函数；
根据药物类型以及组织类型，选择相应的预设药代动力模型。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
所述预设函数包括三指数函数。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
多帧所述扫描图像由动态图像扫描设备扫描获取；
所述动态图像设备包括：正电子发射型计算机断层显像、磁共振成像设备以及电子计算机断层扫描设备。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种药代动力学参数确定方法，其特征在于，所述方法包括：
获取组织药物活度时间曲线以及血浆药物活度时间采样点；
将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数；
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到至少一个药代动力模型的参数值；
根据至少一个所述参数值获得药物代谢动力学参数。
2.根据权利要求1所述的方法，其特征在于，所述获取组织药物活度时间曲线包括：
获取多帧扫描图像；
在各帧所述扫描图像中查找组织相关位置的扫描图像，并将组织相关位置的扫描图像作为第一图像；
提取各帧所述第一图像像素点的像素值，并根据所述像素点的像素值计算各帧组织平均像素值；
将各帧所述组织平均像素值根据扫描时间进行排序，得到所述组织药物活度时间曲线。
3.根据权利要求1所述的方法，其特征在于，所述获取血浆药物活度时间采样点包括：
获取多帧扫描图像；
在各帧所述扫描图像中查找血浆相关的位置的扫描图像，并将血浆相关位置的扫描图像作为第二图像；
提取各帧所述第二图像像素点的像素值，并根据所述像素点的像素值计算各帧血浆平均像素值；
将得到各帧主动脉位置扫描图像的主动脉平均像素值；
将各帧所述血浆平均像素值根据扫描时间进行排序，得到所述血浆药物活度时间采样点。
4.根据权利要求1所述的方法，其特征在于，所述根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到多个药代动力模型的参数值包括：
根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型建立药代动力等式；
对所述药代动力等式进行求解，得到多个药代动力模型的参数值。
5.根据权利要求1所述的方法，其特征在于，所述获取组织药物活度时间曲线以及血浆药物活度时间采样点之后包括：
根据所述血浆药物活度时间采样点，选择相应的预设函数；
根据药物类型以及组织类型，选择相应的预设药代动力模型。
6.根据权利要求5所述的方法，其特征在于，
所述预设函数包括三指数函数。
7.根据权利要求2或3所述的方法，其特征在于，
多帧所述扫描图像由动态图像扫描设备扫描获取；
所述动态图像设备包括：正电子发射型计算机断层显像、磁共振成像设备以及电子计算机断层扫描设备。
8.一种药代动力学参数确定装置，其特征在于，所述装置包括：
获取模块，用于获取组织药物活度时间曲线以及血浆药物活度时间采样点；
拟合模块，用于将所述血浆药物活度时间采样点与预设函数进行拟合，得到血浆药物活度关于时间的模型化输入函数；
参数值得到模块，用于根据所述组织药物活度时间曲线、输入函数以及预设药代动力模型，得到多个药代动力模型的参数值；
药物代谢动力学参数获取模块，用于根据多个所述参数值获得药物代谢动力学参数。
9.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至7中任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910537164.0A CN110269590A (zh) | 2019-06-20 | 2019-06-20 | 药代动力学参数确定方法、装置、计算机设备和存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910537164.0A CN110269590A (zh) | 2019-06-20 | 2019-06-20 | 药代动力学参数确定方法、装置、计算机设备和存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110269590A true CN110269590A (zh) | 2019-09-24 |
Family
ID=67961157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910537164.0A Pending CN110269590A (zh) | 2019-06-20 | 2019-06-20 | 药代动力学参数确定方法、装置、计算机设备和存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110269590A (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110911016A (zh) * | 2019-12-12 | 2020-03-24 | 南京航空航天大学 | 基于非房室模型的血浆及组织的药代动力学参数计算方法 |
| CN113397581A (zh) * | 2021-08-19 | 2021-09-17 | 浙江太美医疗科技股份有限公司 | 重建医学动态影像的方法及装置 |
| WO2022028322A1 (zh) * | 2020-08-03 | 2022-02-10 | 四川海思科制药有限公司 | 确定丙泊酚及其衍生物群体药代动力学模型的方法及系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1916923A (zh) * | 2006-09-01 | 2007-02-21 | 清华大学 | 人体内物质代谢功能信息可视化的计算机模拟方法 |
| CN101223530A (zh) * | 2005-07-21 | 2008-07-16 | 皇家飞利浦电子股份有限公司 | 用于药物代谢动力学建模的自动输入函数估计 |
| US20100209348A1 (en) * | 2009-02-18 | 2010-08-19 | Arondo Pharma, Inc. | Methods for determining liposome bioequivalence |
| CN102013108A (zh) * | 2010-11-23 | 2011-04-13 | 南方医科大学 | 基于区域时空先验的动态pet重建方法 |
| CN104574385A (zh) * | 2014-12-26 | 2015-04-29 | 北京理工大学 | 一种基于动态pet影像的心肌脂肪酸代谢量化检测方法 |
| CN104833748A (zh) * | 2014-02-11 | 2015-08-12 | 沈阳药科大学 | 一种用于中药药代动力学评价的多组分高通量定量分析技术 |
-
2019
- 2019-06-20 CN CN201910537164.0A patent/CN110269590A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101223530A (zh) * | 2005-07-21 | 2008-07-16 | 皇家飞利浦电子股份有限公司 | 用于药物代谢动力学建模的自动输入函数估计 |
| CN1916923A (zh) * | 2006-09-01 | 2007-02-21 | 清华大学 | 人体内物质代谢功能信息可视化的计算机模拟方法 |
| US20100209348A1 (en) * | 2009-02-18 | 2010-08-19 | Arondo Pharma, Inc. | Methods for determining liposome bioequivalence |
| CN102013108A (zh) * | 2010-11-23 | 2011-04-13 | 南方医科大学 | 基于区域时空先验的动态pet重建方法 |
| CN104833748A (zh) * | 2014-02-11 | 2015-08-12 | 沈阳药科大学 | 一种用于中药药代动力学评价的多组分高通量定量分析技术 |
| CN104574385A (zh) * | 2014-12-26 | 2015-04-29 | 北京理工大学 | 一种基于动态pet影像的心肌脂肪酸代谢量化检测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 范丽亚: "《数学建模方法》", 31 December 2015 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110911016A (zh) * | 2019-12-12 | 2020-03-24 | 南京航空航天大学 | 基于非房室模型的血浆及组织的药代动力学参数计算方法 |
| WO2022028322A1 (zh) * | 2020-08-03 | 2022-02-10 | 四川海思科制药有限公司 | 确定丙泊酚及其衍生物群体药代动力学模型的方法及系统 |
| CN113397581A (zh) * | 2021-08-19 | 2021-09-17 | 浙江太美医疗科技股份有限公司 | 重建医学动态影像的方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11020077B2 (en) | Simultaneous CT-MRI image reconstruction | |
| CN110269590A (zh) | 药代动力学参数确定方法、装置、计算机设备和存储介质 | |
| CN110996800B (zh) | 用于确定pet成像动力学参数的系统、方法及非暂时性计算机可读介质 | |
| CN110960241A (zh) | 医学影像扫描的扫描参数确定方法、装置和计算机设备 | |
| CN111080584A (zh) | 医学图像的质控方法、计算机设备和可读存储介质 | |
| CN111887878B (zh) | Pet扫描方法 | |
| JP4446049B2 (ja) | 心筋血流の定量化装置 | |
| WO2008062366A2 (en) | Image generation based on limited data set | |
| CN113329688B (zh) | 用于确定组织体素的特征的方法、系统和计算机可读介质 | |
| CN110223247B (zh) | 图像衰减校正方法、装置、计算机设备和存储介质 | |
| CN113989231A (zh) | 动力学参数的确定方法、装置、计算机设备和存储介质 | |
| CN111402356B (zh) | 参数成像输入函数提取方法、装置和计算机设备 | |
| CN108873047B (zh) | 检测放射源活度的方法、系统、计算机设备和存储介质 | |
| CN112686967A (zh) | 图像重建系统和方法 | |
| WO2022257154A1 (zh) | 参数成像的系统和方法 | |
| CN110507343B (zh) | 运动场获取方法及pet图像校正方法、装置和计算机设备 | |
| EP2174292A1 (en) | A method, apparatus, computer-readable medium and use for pharmacokinetic modeling | |
| CN114677415A (zh) | 图像配准方法、装置、计算机设备及可读存储介质 | |
| CN113712583A (zh) | 死时间校正方法、装置、计算机设备和存储介质 | |
| CN110215203A (zh) | 心电信号获取方法、装置、计算机设备和存储介质 | |
| CN110443868B (zh) | Pet图像获取方法、装置、计算机设备和存储介质 | |
| Klein | Kinetic model based factor analysis of cardiac 82Rb PET images for improved accuracy of quantitative myocardial blood flow measurement | |
| US20230069017A1 (en) | Systems and methods for image reconstruction | |
| CN110992280B (zh) | Pet图像衰减校正方法、装置、计算机设备 | |
| CN118236077A (zh) | 心功能信息确定方法和装置、心肌信息确定方法和装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CB02 | Change of applicant information |