CN116092632B - 一种用于放射性药物评价的核医学成像数据分析方法 - Google Patents

一种用于放射性药物评价的核医学成像数据分析方法 Download PDF

Info

Publication number
CN116092632B
CN116092632B CN202310247911.3A CN202310247911A CN116092632B CN 116092632 B CN116092632 B CN 116092632B CN 202310247911 A CN202310247911 A CN 202310247911A CN 116092632 B CN116092632 B CN 116092632B
Authority
CN
China
Prior art keywords
diffusion
lesion
medicine
coefficient
extracting
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.)
Active
Application number
CN202310247911.3A
Other languages
English (en)
Other versions
CN116092632A (zh
Inventor
李新平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Greepharma Inc ltd
Original Assignee
Nanjing Greepharma Inc ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Greepharma Inc ltd filed Critical Nanjing Greepharma Inc ltd
Priority to CN202310247911.3A priority Critical patent/CN116092632B/zh
Publication of CN116092632A publication Critical patent/CN116092632A/zh
Application granted granted Critical
Publication of CN116092632B publication Critical patent/CN116092632B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • G06T7/0014—Biomedical image inspection using an image reference approach
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
    • G16H20/10—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
    • G16H20/17—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients delivered via infusion or injection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10116—X-ray image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20212—Image combination
    • G06T2207/20221—Image fusion; Image merging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30096—Tumor; Lesion
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Primary Health Care (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Public Health (AREA)
  • Epidemiology (AREA)
  • Chemical & Material Sciences (AREA)
  • Urology & Nephrology (AREA)
  • Surgery (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Medicinal Chemistry (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine (AREA)

Abstract

本发明涉及图像数据处理技术领域，公开了一种用于放射性药物评价的核医学成像数据分析方法，对各个位置的药物分布数据情况进行提取，并与对应位置的病变情况进行比照，同时对药物施加准确位置和剂量进行提取，导入图片生成软件中得到药物扩散演化图，代入图像处理模块中计算药物扩散至各个方向上的扩散系数，且病变系数分析模块对该方向上的病变系数进行提取计算，并将扩散系数和病变系数代入计算模型中，得出药物扩散的最大偏移位置的坐标，对最大偏移位置的坐标进行提取，导入药物注射参数中，对注射位置的具体坐标进行补偿，以便于通过药物的扩散和组织的病变调节对应位置的注射量和注射方向，有利于药物的精准施加，具有广阔的应用前景。

Description

一种用于放射性药物评价的核医学成像数据分析方法
技术领域
本发明涉及图像数据处理技术领域，特别涉及用于放射性药物评价的核医学成像数据分析方法。
背景技术
放射性药物指含有放射性核素供医学诊断和治疗用的一类特殊药物，用于机体内进行医学诊断或治疗的含放射性核素标记的化合物或生物制剂，在进行放射性药物使用效果评价的过程中，一般是在放射性半衰期内使用核医学成像的方式对放射性药物的吸收和位置进行指示，但是药物在注射过程中会发生扩散，扩散后的药物可能就不会对病变位置产生效果，这样需要根据药物的扩散和组织的病变调节对应位置的注射量和注射方向，对注射位置的具体坐标进行补偿，现有的用于放射性药物评价的核医学成像数据分析方法显然不具备这些功能。
如授权公开号为CN109254018B的中国专利提供一种放射性药物的药物代谢动力学成像检测系统，包括：正电子发射断层成像PET探头、设置于所述PET探头内的微流控芯片、以及与所述微流控芯片连接的流体操作单元；其中，所述微流控芯片包括至少一个微流道；所述流体操作单元，用于为所述微流道注入检测液体。
又如授权公开号为CN110599472B的中国专利提供本发明公开了计算SPECT定量断层图像中SUV归一化系数的方法及系统，包括以下步骤：步骤A，对受检者的全身平片扫描所获得的前位平片图像和后位平片图像求取平均图像，在平均图像中分割出与受检者体内有效放射性摄取部位对应的区域；步骤B，对平均图像中的有效放射性摄取区域内的每个像素值进行人体衰减系数、探测效率和采集时间校正；步骤C，计算步骤B得出的二维放射性摄取活度定量图像针对所有有效放射性摄取区域的平均值，并进行体积归一化和衰变时间校正，得出受检者体内的摄取放射性药物活度平均浓度值。相对于基于注射剂量和病人体重的常规方法，此方法提升SUV归一化系数计算的准确性，实现全自动读取和操作。
但药物在注射过程中会发生扩散，扩散后的药物可能就不会对病变位置产生效果，这样需要根据药物的扩散和组织的病变调节对应位置的注射量和注射方向，对注射位置的具体坐标进行补偿，现有的用于放射性药物评价的核医学成像数据分析方法显然不具备这些功能，为解决这一问题，本发明提出了一种用于放射性药物评价的核医学成像数据分析方法。
发明内容
针对现有技术的不足，本发明的主要目的是提供一种用于放射性药物评价的核医学成像数据分析方法，能够有效解决背景技术中的问题：药物在注射过程中会发生扩散，扩散后的药物可能就不会对病变位置产生效果，这样需要根据药物的扩散和组织的病变调节对应位置的注射量和注射方向，对注射位置的具体坐标进行补偿，现有的用于放射性药物评价的核医学成像数据分析方法显然不具备这些功能，本发明的具体技术方案如下：
一种用于放射性药物评价的核医学成像数据分析方法，包括以下具体步骤：
S1：按照设置的划分参数，对治疗部位的放射性药物分布的投影图进行采集，并对治疗部位投影图按照病变位置和未病变位置进行划分，得到病变位置的药物分布投影图和未病变位置的药物分布投影图；
S2：按照得到的病变位置的药物分布投影图和未病变位置的药物分布投影图，将投影图代入像素处理模块中得到治疗部位各个位置投影图的像素分布图，对各个位置的药物分布数据情况进行提取，并与对应位置的病变情况进行比照，同时对药物施加准确位置和剂量进行提取，导入图片生成软件中得到药物扩散演化图；
S3：对药物扩散演化图进行提取，以药物注射中心点为坐标原点，代入图像处理模块中计算药物扩散至各个方向上的扩散系数，且病变系数分析模块对该方向上的病变系数进行提取计算，并将扩散系数和病变系数代入计算模型中，得出药物扩散的最大偏移位置的坐标；
S4：对最大偏移位置的坐标进行提取，导入药物注射参数中，对注射位置的具体坐标进行补偿。
本发明进一步的改进在于，所述S1中的具体步骤为：
S101：对患者治疗部位进行X光扫描，对治疗部位X光扫描底片进行提取，找出患者病变最严重的位置，以该位置作为坐标原点
Figure SMS_1
，对患者病变最严重位置注射放射性药物，使用放射性检测设备在放射性半衰期内对治疗位置的放射性药物的投影图进行提取；
S102：使用图片融合软件对放射性药物的投影图和扫描底片进行融合，得到治疗部位的药物分布和治疗部位的融合图片；
S103：将融合后的图片按照病变位置和未病变位置进行划分，以得到病变位置的药物分布投影图和未病变位置的药物分布投影图。
本发明进一步的改进在于，所述S2的具体步骤为：
S201：将病变位置的药物分布投影图和未病变位置的药物分布投影图代入像素处理模块中，得到治疗部位各个位置的投影图像素分布图；
S202：读取治疗部位各个位置的投影图像素分布，对各个位置的药物分布数据情况进行提取，并提取对应位置的病变情况数据，所述病变情况数据以病情恶劣情况打分，在投影图图片中提取药物分布图层和病变情况数据图层；
S203：同时在药物分布图层中对药物施加准确位置和剂量进行标注，从而得到药物扩散演化图。
本发明进一步的改进在于，所述S3的具体步骤为：
S301：对药物扩散演化图进行提取，以药物注射中心点为圆心点，即患者病变最严重的位置，药物注射聚集区域的点至圆心点距离的最大值和最小值的和的一半作为半径画圆，作为基准圆，其数学表达式为，
Figure SMS_2
，其中/>
Figure SMS_3
是坐标原点，即为圆心点，/>
Figure SMS_4
为药物注射聚集区域的点至圆心点距离的最大值，/>
Figure SMS_5
为药物注射聚集区域的点至圆心点距离的最小值；
S302：将提取的药物扩散演化图、药物注射中心点位置和半径代入图像处理模块中计算药物扩散至各个方向上的扩散系数，扩散方向以坐标圆点为零点以30度为一个方向扩散划分；
S303：将扩散方向上的病变数据代入病变系数分析模块中，病变系数分析模块对每个扩散方向上的病变系数进行提取计算；
S304：将计算得到的扩散系数和病变系数代入计算模型中，并进行药物偏移位置偏移量计算，得出药物扩散的最大偏移位置的坐标。
本发明进一步的改进在于，所述S4的具体步骤为：
S401、对S3中得到的最大偏移位置的坐标进行提取，并提取其扩散方向的对应位置，同时提取与其扩散方向呈中心对称的对应扩散方向；
S402、将最大偏移位置的坐标、该扩散方向上的药物扩散量进行提取，导入药物注射参数中；
S403、将药物注射参数中的注射位置按照原注射位置向着对应扩散方向偏转最大偏移位置的坐标进行注射。
本发明进一步的改进在于，所述S302中包括扩散方向上扩散系数的计算策略，所述扩散方向上扩散系数的计算策略的具体步骤为：
S3021：对其中一个扩散方向上的药物扩散距离和药物扩散量进行提取，提取方式为以药物注射中心点为中心，提取药物扩散的距离集合
Figure SMS_6
，其中/>
Figure SMS_7
为n倍基准圆半径减去n-1倍基准圆半径的范围，和扩散距离上的药物扩散量
Figure SMS_8
，其中/>
Figure SMS_9
为/>
Figure SMS_10
范围内的药物扩散量；
S3022：将药物扩散距离和药物扩散量代入扩散系数计算公式中计算该扩散方向上的扩散系数，其中扩散系数计算公式为：
Figure SMS_11
，其中/>
Figure SMS_12
为药物扩散的距离集合中的第i项，/>
Figure SMS_13
为/>
Figure SMS_14
范围内的药物扩散量，X为整体注射量；
S3023：对各个扩散方向上的扩散系数进行计算，得到扩散系数序列，计算中心对称位置上对应两个扩散方向的扩散系数差；
本发明进一步的改进在于，所述S303步骤包括以下具体内容：
S3031:对其中一个扩散方向上的病变数据进行提取，提取病变扩散的距离集合
Figure SMS_15
，和对应距离上的病变像素灰度相对值集合/>
Figure SMS_16
，其中病变像素灰度相对值计算方式为对应位置的病变像素灰度与健康的位置像素灰度的比较值，/>
Figure SMS_17
为/>
Figure SMS_18
范围内的病变像素灰度相对值；
S3032:将病变扩散距离和对应距离上的病变像素灰度相对值代入病变系数计算公式中计算该扩散方向上的病变系数，其中病变系数的计算公式为：
Figure SMS_19
，其中/>
Figure SMS_20
为药物扩散的距离集合中的第i项，/>
Figure SMS_21
为病变像素灰度相对值集合中的第i项，/>
Figure SMS_22
为病灶中心的病变像素灰度相对值；
S3033:对各个扩散方向上的病变系数进行计算，得到病变系数序列，储存。
本发明进一步的改进在于，所述S304步骤包括以下具体内容：
所述S304步骤包括以下具体内容：
S3041:将计算得到的扩散系数和病变系数按照其所对应的扩散方向划分为多个序列集合，划分个数与扩散方向的数量相同，将数据代入计算模型中；
S3042:所述计算模型包括以下具体内容：计算各个方向上扩散系数与病变系数的相差值，保存为序列集的形式，在各个方向上扩散系数与病变系数的相差值序列集中找出最大值，此最大值对应的扩散方向设为风险方向，将风险方向对应的扩散方向到圆点所成的扇形区域按照
Figure SMS_23
为划分长度划分为若干个区域，对各个区域中的药物密度和病变系数进行计算，药物密度的计算公式为对应区域的药物总量除以区域面积，对该扩散方向上的最大偏移位置进行提取，提取方式为，将扩散方向上各个区域药物密度和病变系数比值进行计算，得到比值中的最大值，以比值中的最大值所在区域的中心点相对于圆点坐标的位置作为最大偏移位置；
S3043:在下次注射时对于风险方向的注射，向着风险方向呈中心对称的对应扩散方向偏转最大偏移位置。
与现有技术相比，本发明具有如下有益效果：
本发明能够对治疗部位的放射性药物分布的投影图进行采集，并对治疗部位投影图按照病变位置和未病变位置进行划分，得到病变位置的药物分布投影图和未病变位置的药物分布投影图将投影图代入像素处理模块中得到治疗部位各个位置投影图的像素分布图，对各个位置的药物分布数据情况进行提取，并与对应位置的病变情况进行比照，同时对药物施加准确位置和剂量进行提取，导入图片生成软件中得到药物扩散演化图，代入图像处理模块中计算药物扩散至各个方向上的扩散系数，且病变系数分析模块对该方向上的病变系数进行提取计算，并将扩散系数和病变系数代入计算模型中，得出药物扩散的最大偏移位置的坐标，对最大偏移位置的坐标进行提取，导入药物注射参数中，对注射位置的具体坐标进行补偿，以便于通过药物的扩散和组织的病变调节对应位置的注射量和注射方向，有利于药物的精准施加，具有广阔的应用前景。
附图说明
图1为本发明一种用于放射性药物评价的核医学成像数据分析方法的流程示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合说明书附图对本发明的具体实施方式做详细的说明。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是本发明还可以采用其他不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广，因此本发明不受下面公开的具体实施例的限制。
其次，此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例，也不是单独的或选择性的与其他实施例互相排斥的实施例。
实施例
本实施例通过对治疗部位的放射性药物分布的投影图进行采集，并对治疗部位投影图按照病变位置和未病变位置进行划分，得到病变位置的药物分布投影图和未病变位置的药物分布投影图将投影图代入像素处理模块中得到治疗部位各个位置投影图的像素分布图，对各个位置的药物分布数据情况进行提取，并与对应位置的病变情况进行比照，同时对药物施加准确位置和剂量进行提取，导入图片生成软件中得到药物扩散演化图，代入图像处理模块中计算药物扩散至各个方向上的扩散系数，且病变系数分析模块对该方向上的病变系数进行提取计算，并将扩散系数和病变系数代入计算模型中，得出药物扩散的最大偏移位置的坐标，对最大偏移位置的坐标进行提取，导入药物注射参数中，对注射位置的具体坐标进行补偿，以便于通过药物的扩散和组织的病变调节对应位置的注射量和注射方向，有利于药物的精准施加，具有广阔的应用前景，具体方案为，如图1所示，一种用于放射性药物评价的核医学成像数据分析方法，包括以下具体步骤：
S1：按照设置的划分参数，对治疗部位的放射性药物分布的投影图进行采集，并对治疗部位投影图按照病变位置和未病变位置进行划分，得到病变位置的药物分布投影图和未病变位置的药物分布投影图；
其中，S1中的具体步骤为：
S101：对患者治疗部位进行X光扫描，对治疗部位X光扫描底片进行提取，找出患者病变最严重的位置，以该位置作为坐标原点
Figure SMS_24
，对患者病变最严重位置注射放射性药物，使用放射性检测设备在放射性半衰期内对治疗位置的放射性药物的投影图进行提取；
S102：使用图片融合软件对放射性药物的投影图和扫描底片进行融合，得到治疗部位的药物分布和治疗部位的融合图片；
S103：将融合后的图片按照病变位置和未病变位置进行划分，以得到病变位置的药物分布投影图和未病变位置的药物分布投影图；
S2：按照得到的病变位置的药物分布投影图和未病变位置的药物分布投影图，将投影图代入像素处理模块中得到治疗部位各个位置投影图的像素分布图，对各个位置的药物分布数据情况进行提取，并与对应位置的病变情况进行比照，同时对药物施加准确位置和剂量进行提取，导入图片生成软件中得到药物扩散演化图；
其中，S2的具体步骤为：
S201：将病变位置的药物分布投影图和未病变位置的药物分布投影图代入像素处理模块中，得到治疗部位各个位置的投影图像素分布图；
S202：读取治疗部位各个位置的投影图像素分布，对各个位置的药物分布数据情况进行提取，并提取对应位置的病变情况数据，病变情况数据以病情恶劣情况打分，在投影图图片中提取药物分布图层和病变情况数据图层；
S203：同时在药物分布图层中对药物施加准确位置和剂量进行标注，从而得到药物扩散演化图；
S3：对药物扩散演化图进行提取，以药物注射中心点为坐标原点，代入图像处理模块中计算药物扩散至各个方向上的扩散系数，且病变系数分析模块对该方向上的病变系数进行提取计算，并将扩散系数和病变系数代入计算模型中，得出药物扩散的最大偏移位置的坐标；
其中，S3的具体步骤为：
S301：对药物扩散演化图进行提取，以药物注射中心点为圆心点，即患者病变最严重的位置，药物注射聚集区域的点至圆心点距离的最大值和最小值的和的一半作为半径画圆，作为基准圆，其数学表达式为，
Figure SMS_25
，其中/>
Figure SMS_26
是坐标原点，即为圆心点，/>
Figure SMS_27
为药物注射聚集区域的点至圆心点距离的最大值，/>
Figure SMS_28
为药物注射聚集区域的点至圆心点距离的最小值；
S302：将提取的药物扩散演化图、药物注射中心点位置和半径代入图像处理模块中计算药物扩散至各个方向上的扩散系数，扩散方向以坐标圆点为零点以30度为一个方向扩散划分；
其中，S302中包括扩散方向上扩散系数的计算策略，扩散方向上扩散系数的计算策略的具体步骤为：
S3021：对其中一个扩散方向上的药物扩散距离和药物扩散量进行提取，提取方式为以药物注射中心点为中心，提取药物扩散的距离集合
Figure SMS_29
，其中/>
Figure SMS_30
为n倍基准圆半径减去n-1倍基准圆半径的范围，和扩散距离上的药物扩散量
Figure SMS_31
，其中/>
Figure SMS_32
为/>
Figure SMS_33
范围内的药物扩散量；
S3022：将药物扩散距离和药物扩散量代入扩散系数计算公式中计算该扩散方向上的扩散系数，其中扩散系数计算公式为：
Figure SMS_34
，其中/>
Figure SMS_35
为药物扩散的距离集合中的第i项，/>
Figure SMS_36
为/>
Figure SMS_37
范围内的药物扩散量，X为整体注射量；
S3023：对各个扩散方向上的扩散系数进行计算，得到扩散系数序列，计算中心对称位置上对应两个扩散方向的扩散系数差；
S303：将扩散方向上的病变数据代入病变系数分析模块中，病变系数分析模块对每个扩散方向上的病变系数进行提取计算；
其中，S303步骤包括以下具体内容：
S3031:对其中一个扩散方向上的病变数据进行提取，提取病变扩散的距离集合
Figure SMS_38
，和对应距离上的病变像素灰度相对值集合/>
Figure SMS_39
，其中病变像素灰度相对值计算方式为对应位置的病变像素灰度与健康的位置像素灰度的比较值，/>
Figure SMS_40
为/>
Figure SMS_41
范围内的病变像素灰度相对值；
S3032:将病变扩散距离和对应距离上的病变像素灰度相对值代入病变系数计算公式中计算该扩散方向上的病变系数，其中病变系数的计算公式为：
Figure SMS_42
，其中/>
Figure SMS_43
为药物扩散的距离集合中的第i项，/>
Figure SMS_44
为病变像素灰度相对值集合中的第i项，/>
Figure SMS_45
为病灶中心的病变像素灰度相对值；
S3033:对各个扩散方向上的病变系数进行计算，得到病变系数序列，储存；
S304：将计算得到的扩散系数和病变系数代入计算模型中，并进行药物偏移位置偏移量计算，得出药物扩散的最大偏移位置的坐标；
其中，S304步骤包括以下具体内容：
所述S304步骤包括以下具体内容：
S3041:将计算得到的扩散系数和病变系数按照其所对应的扩散方向划分为多个序列集合，划分个数与扩散方向的数量相同，将数据代入计算模型中；
S3042:计算模型包括以下具体内容：计算各个方向上扩散系数与病变系数的相差值，保存为序列集的形式，在各个方向上扩散系数与病变系数的相差值序列集中找出最大值，此最大值对应的扩散方向设为风险方向，将风险方向对应的扩散方向到圆点所成的扇形区域按照
Figure SMS_46
为划分长度划分为若干个区域，对各个区域中的药物密度和病变系数进行计算，药物密度的计算公式为对应区域的药物总量除以区域面积，对该扩散方向上的最大偏移位置进行提取，提取方式为，将扩散方向上各个区域药物密度和病变系数比值进行计算，得到比值中的最大值，以比值中的最大值所在区域的中心点相对于圆点坐标的位置作为最大偏移位置；
S3043:在下次注射时对于风险方向的注射，向着风险方向呈中心对称的对应扩散方向偏转最大偏移位置；
S4：对最大偏移位置的坐标进行提取，导入药物注射参数中，对注射位置的具体坐标进行补偿；
其中，S4的具体步骤为：
S401、对S3中得到的最大偏移位置的坐标进行提取，并提取其扩散方向的对应位置，同时提取与其扩散方向呈中心对称的对应扩散方向；
S402、将最大偏移位置的坐标、该扩散方向上的药物扩散量进行提取，导入药物注射参数中；
S403、将药物注射参数中的注射位置按照原注射位置向着对应扩散方向偏转最大偏移位置的坐标进行注射。
通过本实施例能够实现：能够对治疗部位的放射性药物分布的投影图进行采集，并对治疗部位投影图按照病变位置和未病变位置进行划分，得到病变位置的药物分布投影图和未病变位置的药物分布投影图将投影图代入像素处理模块中得到治疗部位各个位置投影图的像素分布图，对各个位置的药物分布数据情况进行提取，并与对应位置的病变情况进行比照，同时对药物施加准确位置和剂量进行提取，导入图片生成软件中得到药物扩散演化图，代入图像处理模块中计算药物扩散至各个方向上的扩散系数，且病变系数分析模块对该方向上的病变系数进行提取计算，并将扩散系数和病变系数代入计算模型中，得出药物扩散的最大偏移位置的坐标，对最大偏移位置的坐标进行提取，导入药物注射参数中，对注射位置的具体坐标进行补偿，以便于通过药物的扩散和组织的病变调节对应位置的注射量和注射方向，有利于药物的精准施加，具有广阔的应用前景。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (8)

1.一种用于放射性药物评价的核医学成像数据分析方法，其特征在于：包括以下具体步骤：
S1：按照设置的划分参数，对治疗部位的放射性药物分布的投影图进行采集，并对治疗部位投影图按照病变位置和未病变位置进行划分，得到病变位置的药物分布投影图和未病变位置的药物分布投影图；
S2：按照得到的病变位置的药物分布投影图和未病变位置的药物分布投影图，将投影图代入像素处理模块中得到治疗部位各个位置投影图的像素分布图，对各个位置的药物分布数据情况进行提取，并与对应位置的病变情况进行比照，同时对药物施加准确位置和剂量进行提取，导入图片生成软件中得到药物扩散演化图；
S3：对药物扩散演化图进行提取，以药物注射中心点为坐标原点，代入图像处理模块中计算药物扩散至各个方向上的扩散系数，且病变系数分析模块对该方向上的病变系数进行提取计算，并将扩散系数和病变系数代入计算模型中，得出药物扩散的最大偏移位置的坐标；
S4：对最大偏移位置的坐标进行提取，导入药物注射参数中，对注射位置的具体坐标进行补偿。
2.根据权利要求1所述的一种用于放射性药物评价的核医学成像数据分析方法，其特征在于：所述S1中的具体步骤为：
S101：对患者治疗部位进行X光扫描，对治疗部位X光扫描底片进行提取，找出患者病变最严重的位置，以该位置作为坐标原点，对患者病变最严重位置注射放射性药物，使用放射性检测设备在放射性半衰期内对治疗位置的放射性药物的投影图进行提取；
S102：使用图片融合软件对放射性药物的投影图和扫描底片进行融合，得到治疗部位的药物分布和治疗部位的融合图片；
S103：将融合后的图片按照病变位置和未病变位置进行划分，以得到病变位置的药物分布投影图和未病变位置的药物分布投影图。
3.根据权利要求2所述的一种用于放射性药物评价的核医学成像数据分析方法，其特征在于：所述S2的具体步骤为：
S201：将病变位置的药物分布投影图和未病变位置的药物分布投影图代入像素处理模块中，得到治疗部位各个位置的投影图像素分布图；
S202：读取治疗部位各个位置的投影图像素分布图，对各个位置的药物分布数据情况进行提取，并提取对应位置的病变情况数据，所述病变情况数据以病情恶劣情况打分，在投影图图片中提取药物分布图层和病变情况数据图层；
S203：同时在药物分布图层中对药物施加准确位置和剂量进行标注，从而得到药物扩散演化图。
4.根据权利要求3所述的一种用于放射性药物评价的核医学成像数据分析方法，其特征在于：所述S3的具体步骤为：
S301：对药物扩散演化图进行提取，以药物注射中心点为圆心点，即患者病变最严重的位置，药物注射聚集区域的点至圆心点距离的最大值和最小值的和的一半作为半径画圆，作为基准圆，其数学表达式为，
Figure QLYQS_1
，其中/>
Figure QLYQS_2
是坐标原点，即为圆心点，/>
Figure QLYQS_3
为药物注射聚集区域的点至圆心点距离的最大值，/>
Figure QLYQS_4
为药物注射聚集区域的点至圆心点距离的最小值；
S302：将提取的药物扩散演化图、药物注射中心点位置和半径代入图像处理模块中计算药物扩散至各个方向上的扩散系数，扩散方向以坐标圆点为零点以30度为一个方向扩散划分；
S303：将扩散方向上的病变数据代入病变系数分析模块中，病变系数分析模块对每个扩散方向上的病变系数进行提取计算；
S304：将计算得到的扩散系数和病变系数代入计算模型中，并进行药物偏移位置偏移量计算，得出药物扩散的最大偏移位置的坐标。
5.根据权利要求4所述的一种用于放射性药物评价的核医学成像数据分析方法，其特征在于：所述S4的具体步骤为：
S401、对S3中得到的最大偏移位置的坐标进行提取，并提取其扩散方向的对应位置，同时提取与其扩散方向呈中心对称的对应扩散方向；
S402、将最大偏移位置的坐标、该扩散方向上的药物扩散量进行提取，导入药物注射参数中；
S403、将药物注射参数中的注射位置按照原注射位置向着对应扩散方向偏转最大偏移位置的坐标进行注射。
6.根据权利要求5所述的一种用于放射性药物评价的核医学成像数据分析方法，其特征在于：所述S302中包括扩散方向上扩散系数的计算策略，所述扩散方向上扩散系数的计算策略的具体步骤为：
S3021：对其中一个扩散方向上的药物扩散距离和药物扩散量进行提取，提取方式为以药物注射中心点为中心，提取药物扩散的距离集合
Figure QLYQS_5
，其中/>
Figure QLYQS_6
为n倍基准圆半径减去n-1倍基准圆半径的范围，和扩散距离上的药物扩散量/>
Figure QLYQS_7
，其中/>
Figure QLYQS_8
为/>
Figure QLYQS_9
范围内的药物扩散量；
S3022：将药物扩散距离和药物扩散量代入扩散系数计算公式中计算该扩散方向上的扩散系数，其中扩散系数计算公式为：
Figure QLYQS_10
，其中/>
Figure QLYQS_11
为药物扩散的距离集合中的第i项，/>
Figure QLYQS_12
为/>
Figure QLYQS_13
范围内的药物扩散量，X为整体注射量；
S3023：对各个扩散方向上的扩散系数进行计算，得到扩散系数序列，计算中心对称位置上对应两个扩散方向的扩散系数差。
7.根据权利要求6所述的一种用于放射性药物评价的核医学成像数据分析方法，其特征在于：所述S303步骤包括以下具体内容：
S3031:对其中一个扩散方向上的病变数据进行提取，提取病变扩散的距离集合
Figure QLYQS_14
，和对应距离上的病变像素灰度相对值集合/>
Figure QLYQS_15
，其中病变像素灰度相对值计算方式为对应位置的病变像素灰度与健康的位置像素灰度的比较值，/>
Figure QLYQS_16
为/>
Figure QLYQS_17
范围内的病变像素灰度相对值；
S3032:将病变扩散距离和对应距离上的病变像素灰度相对值代入病变系数计算公式中计算该扩散方向上的病变系数，其中病变系数的计算公式为：
Figure QLYQS_18
，其中/>
Figure QLYQS_19
为药物扩散的距离集合中的第i项，/>
Figure QLYQS_20
为病变像素灰度相对值集合中的第i项，/>
Figure QLYQS_21
为病灶中心的病变像素灰度相对值；
S3033:对各个扩散方向上的病变系数进行计算，得到病变系数序列，储存。
8.根据权利要求7所述的一种用于放射性药物评价的核医学成像数据分析方法，其特征在于：所述S304步骤包括以下具体内容：
S3041:将计算得到的扩散系数和病变系数按照其所对应的扩散方向划分为多个序列集合，划分个数与扩散方向的数量相同，将数据代入计算模型中；
S3042:所述计算模型包括以下具体内容：计算各个方向上扩散系数与病变系数的相差值，保存为序列集的形式，在各个方向上扩散系数与病变系数的相差值序列集中找出最大值，此最大值对应的扩散方向设为风险方向，将风险方向对应的扩散方向到圆点所成的扇形区域按照
Figure QLYQS_22
为划分长度划分为若干个区域，对各个区域中的药物密度和病变系数进行计算，药物密度的计算公式为对应区域的药物总量除以区域面积，对该扩散方向上的最大偏移位置进行提取，提取方式为，将扩散方向上各个区域药物密度和病变系数比值进行计算，得到比值中的最大值，以比值中的最大值所在区域的中心点相对于圆点坐标的位置作为最大偏移位置；
S3043:在下次注射时对于风险方向的注射，向着风险方向呈中心对称的对应扩散方向偏转最大偏移位置。
CN202310247911.3A 2023-03-15 2023-03-15 一种用于放射性药物评价的核医学成像数据分析方法 Active CN116092632B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310247911.3A CN116092632B (zh) 2023-03-15 2023-03-15 一种用于放射性药物评价的核医学成像数据分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310247911.3A CN116092632B (zh) 2023-03-15 2023-03-15 一种用于放射性药物评价的核医学成像数据分析方法

Publications (2)

Publication Number Publication Date
CN116092632A CN116092632A (zh) 2023-05-09
CN116092632B true CN116092632B (zh) 2023-06-30

Family

ID=86204711

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310247911.3A Active CN116092632B (zh) 2023-03-15 2023-03-15 一种用于放射性药物评价的核医学成像数据分析方法

Country Status (1)

Country Link
CN (1) CN116092632B (zh)

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7983735B2 (en) * 2003-04-15 2011-07-19 General Electric Company Simulation of nuclear medical imaging
WO2016064750A1 (en) * 2014-10-20 2016-04-28 New York University Method to select radiation dosage for tumor treatment based on cellular imaging
DE102014221634A1 (de) * 2014-10-24 2016-04-28 Siemens Aktiengesellschaft Verfahren zum Bestimmen einer Strahlendosis eines Radiopharmazeutikums
DK178729B1 (en) * 2015-10-06 2016-12-12 Nordic Healthcare Advisory Aps Positioning Device for determining and registering an injection site
WO2017084004A1 (zh) * 2015-11-16 2017-05-26 数码医疗集团 治疗计划制定方法、装置及放疗系统
CN107153771B (zh) * 2017-05-17 2020-06-30 许昌学院 一种药物分子的同步控制方法及其应用
US20200163642A1 (en) * 2017-08-08 2020-05-28 University Of Vermont And State Agricultural College Methods for computational modeling to guide intratumoral therapy
DE102017009040A1 (de) * 2017-09-27 2019-03-28 Städtisches Klinikum Dessau Verfahren zur EPID-basierten Überprüfung, Korrektur und Minimierung des Isozentrums eines Strahlentherapiegeräts
CN111615365B (zh) * 2018-12-26 2021-11-30 西安大医集团股份有限公司 摆位方法、装置及放射治疗系统
EP4022632A1 (en) * 2019-08-30 2022-07-06 Eli Lilly and Company Methods and apparatus for determining device dosage
CN110575625A (zh) * 2019-09-18 2019-12-17 北京大学第三医院 一种放射治疗组织等效补偿物的制作方法
CN112863641A (zh) * 2019-11-12 2021-05-28 西安大医集团股份有限公司 放射治疗系统及其放射源的偏移量确定方法、装置
CN114880960A (zh) * 2021-07-15 2022-08-09 青岛埃米博创医疗科技有限公司 基于流体动力学的放射栓塞剂量及注射位置评估方法
CN115068842A (zh) * 2022-06-15 2022-09-20 宁波市第一医院 后装近距离放射治疗计划验证方法

Also Published As

Publication number Publication date
CN116092632A (zh) 2023-05-09

Similar Documents

Publication Publication Date Title
Flotats et al. Proposal for standardization of 123 I-metaiodobenzylguanidine (MIBG) cardiac sympathetic imaging by the EANM Cardiovascular Committee and the European Council of Nuclear Cardiology
CN110599472B (zh) 计算spect定量断层图像中suv归一化系数的方法及系统
JP4414410B2 (ja) 画像再構成方法
Wu et al. Absolute quantitative total-body small-animal SPECT with focusing pinholes
Keyes et al. Calculation of viable and infarcted myocardial mass from thallium-201 tomograms
CN105476653A (zh) Spect骨显像的定量分析技术及在骨评估中的用途
CN104605875A (zh) 符合线路spect肿瘤显像的定量分析技术及在肿瘤评估中的用途
Tang et al. Implementation of a combined X-ray CT-scintillation camera imaging system for localizing and measuring radionuclide uptake: experiments in phantoms and patients
Khanzadeh et al. Estimating the absorbed dose of organs in pediatric imaging of 99mTc-DTPATc-DTPA radiopharmaceutical using MIRDOSE software
Michel et al. Online brain attenuation correction in PET: towards a fully automated data handling in a clinical environment
CN116092632B (zh) 一种用于放射性药物评价的核医学成像数据分析方法
Wierts et al. Dosimetry in molecular nuclear therapy
Hays Kinetics of the human thyroid trap: a compartmental model
Elawadi et al. The effect of contrast agents on dose calculations of volumetric modulated arc radiotherapy plans for critical structures
Webb et al. Single photon emission computed tomographic imaging and volume estimation of the thyroid using fan-beam geometry
Asl et al. Prediction of Absorbed Dose to Normal Organs with Endocrine Tumors for I-131 by use of 99mTC Single Photon Emission Computed Tomography/Computed Tomography and Geant4 Application for Tomographic Emission Simulation
Kimiaei et al. Collimator design for improved spatial resolution in SPECT and planar scintigraphy
Strigari et al. Evolution of Portable Sensors for In-Vivo Dose and Time-Activity Curve Monitoring as Tools for Personalized Dosimetry in Molecular Radiotherapy
Ejeh Accuracy of iodine-131 activity quantification and dosimetry for three-dimensional patient-specific models
Slart et al. Sex-based differences in nuclear medicine imaging and therapy
Erojo Quantitative I-131 SPECT for thyroid cancer radionuclide therapy
Lassmann et al. Thyroid Imaging and Dosimetry
Tagesson et al. Transformation of activity distribution in quantitative SPECT to absorbed dose distribution in a radionuclide treatment planning system
Barrack Practical Implementation of Dosimetry for Patients with Thyroid Carcinoma Undergoing I-131 Therapy
Hademenos et al. An improved method for simultaneous dual isotope Tc-99m/Tl-201 emission scintigraphy

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
GR01 Patent grant
GR01 Patent grant