《肿瘤放射物理学》三维剂量计算模型和治疗方案优化.pptVIP

高能X (γ) 射线的三维剂量计算模型 IMRT治疗方案优化 高能X (γ) 射线的三维剂量计算模型 治疗计划的设计步骤 射野剂量分布的数字表达 三维剂量计算模型中，最常用的剂量数据的表达方式是剂量网格矩阵 直角坐标系 极坐标系 扇形线坐标系 等离轴比线坐标系 坐标系的选择取决于采用的剂量计算模型 1、直角坐标系 最直观的网格矩阵表示方法 网格点大小按剂量精度选择 射野边缘，剂量变化梯度较大的地方，应选较小网格点 缺点：需存储的数据量较大 2、极坐标系 特别适用于旋转照射的剂量计算 剂量等中心处归一 3、扇形线坐标系 网格点为源于放射源的扇形线与平行于模体（或患者皮肤表面）的平行线（深度线）的交点 扇形线间距随偏离中心轴而变化，半影区扇形线较密 4、等离轴比线坐标系 等离轴比线为位于同一离轴比线上诸点的剂量率与模体内同一深度处中心轴上剂量率之比为常数 剂量网格点为平行于模体的不同深度线与相应离轴比线的交点 射野边缘等离轴比线靠得比较近，射野中心区域的等离轴比线相距较远 剂量计算模型 剂量计算模型： 基于修正的剂量计算模式 以参考条件下对标准野测量深度剂量、离轴比、散射因子等为基础，加之必要修正后得到实际放射物理条件和放射治疗对象剂量分布的剂量计算方法。 矩阵法 解析法 矩阵法 楔形板、射野挡块、组织补偿器、曲面和不均匀组织对剂量分布的影响，按一维校正因子逐点计算与修正。 等剂量线叠加在矩阵网格上，形成射野剂量矩阵。 多个射野剂量矩阵的叠加，形成多野照射的剂量分布 解析法 剂量分布通常表示为两个函数的乘积。 二维情况下： 矩阵法和解析法依赖于实验条件和实验模型，用于实际剂量计算时，通常要根据实际情况进行校正，如射野形状、源皮距、人体曲面、非均匀组织、楔形板等 主要用于计算二维平面上的剂量分布 剂量计算模型 基于模型的剂量计算模式 原射线和散射线剂量分离法 卷积法 蒙特卡罗法 1、原射线剂量和散射线剂量分离法 在电子平衡的条件下，将体模中任意一点的剂量分为原射线剂量和散射线剂量两部分，即体模内的吸收剂量为： 半经验公式模型 原射线剂量：用零射野的组织空气比(TAR)表示 最大剂量空气比(TMR)，组织模体比(TPR) 散射线剂量：用散射空气比(SAR)表示 散射最大剂量比(SMR)，散射模体比(SPR) 射野内任一点的剂量为： 半经验公式模型 在均匀模体中或不规则野的剂量计算，有较高计算精度 没有考虑计算点周围不均匀组织对剂量的影响，不能用于三维计算 对半经验公式模型的改进 等效组织空气比(ETAR)法 有效深度d′:修正扇形线方向上不均匀组织引起的深度变化 有效射野r′:组织密度和散射点与计算点间相对位置的函数 微分散射空气比(DSAR) 微分散射空气比(DSAR)：散射体积单元△V内的散射空气比 微分散射空气比(DSAR) 柱形环单元内r · △r · △θ内的散射空气比： 微分散射空气比(DSAR) 对不均匀组织，P点的散射剂量 2、卷积法 剂量计算是通过将放射性的光通分布(fluence distribution)与一个点扩展函数(point spread distribution)或卷积核相卷积来实现 卷积核的获取 直接实验测量法（粗糙的方法） 采用适当的限光和准直装置，把准直器散射降到最低，把测量得到的小野剂量分布的二维截面剂量分布作为二维卷积核 有限面积的近似卷积核 高斯函数逼近法 用高斯函数的解析式逼近三（二）维卷积核 小角度散射 Monte Carlo模拟法（高精度方法） 等同于在均匀介质中用Monte Carlo方法进行剂量计算 解卷积法（简单可靠） 从测量数据中用解卷积方法抽取卷积核的方法 卷积核包含了全部光子与物质作用的全部物理信息，获取容易可靠 3. 蒙特卡罗法(Monte Carlo) 原理： 用Monte Carlo方法来模拟大量单个光子在输运过程中与物质的作用过程，通过对作用过程的随机采样，对每一次作用应用放射物理定律来预测和统计作用的结果。 Monte Carlo技术是用随机抽样技术模拟三个过程： 原射线的能谱及其离轴分布； 原射线及散射线光子在介质中的输运过程； 模拟由光子与物质相互作用后产生的次级电子的输运和能量沉积过程。 Monte Carlo是一种适用性最强的三维剂量计算方法，但计算时间太长无法满足临床的实时要求。 IMRT治疗方案优化 IMRT计划优化 优化目标函数和优化搜索算法是IMRT优化的重要内容 目标函数(Objective Function) 目标函数 以数学公式表示期望分布和实际分布差异的函数形式 最简单的目标函数