专利名称：分段锥形束ct图像通过拼接获得完整解剖图像的方法
技术领域：
本发明涉及一种医学影像制作技术，具体地说是一种分段锥形束CT图像通过 拼接获得完整解剖图像的方法。
2背景技术：
锥形束CT图像的概念锥形束CT，即Cone Beam CT，简称CBCT，是医学影 像的一种，与普通CT影像类似，也是反映人体解剖结构的断层图像，广泛应用于 肿瘤放射治疗领域。该图像由医用直线加速器附带的"在线图像采集器"(简称OBI) 拍摄得到。当患者在这种加速器上接受放射治疗时，医技人员可以控制加速器机 头与0BI围绕躺在治疗床上的患者同心共面旋转。旋转一圈360。后，加速器控制 端就存储下了 0BI采集到的一套CBCT图像序列，该序列经三维重建后可得到一定 空间范围内的解剖图像。
锥形束CT图像的作用
第一放疗科医生可以将在线获得的一套CBCT图像片段立即导入加速器工作 站，通过图像的三维重建以及与患者定位CT图像的融合，实时地校正患者的治疗 体位，校正完毕后再启动加速器对患者实施实际的照射。这样显然提高了放射治 疗的精度，这种治疗方式通常称为"影像引导放射治疗。
第二患者完成一次治疗后，放疗科医生可以将在线获得的一套CBCT图像导 入到放疗计划系统，通过放疗计划系统对图像进行三维重建。通过对比患者定位 CT图像和重建后的CBCT图像，医生可以观察肿瘤在这次治疗时与若干天前采集定 位CT图像时相比是否发生了变化，如果变化较大，医生将考虑及时地修改患者的 放疗计划。
第三医生在放疗计划系统上对CBCT图像进行三维重建后，可以将患者已经 实施的治疗计划叠加于这套图像之上，并由放疗计划系统进行放射剂量的计算。 如果CBCT图像空间中完整地再现了患者的肿瘤和肿瘤周围的正常l且织，那么医生 就可以根据这一计算结果来评估患者接受放疗时肿瘤和正常组织实际受到照射的 剂量。如果剂量计算结果与患者己经接受的治疗计划有较大差别，包括肿瘤未得 到预期剂量的照射和肿瘤周围正常器官受到了比预期更高的剂量照射，医生将考 虑及时地修改患者的放疗计划。目前，美国Varian公司提供的医用直线加速器上附带了这种OBI设备，该设 备目前已在国内十余家大型肿瘤治疗中心投入使用。由于其机械设计的局限，它 所获取的CBCT图像在人体头脚方向的长度〈14cra。如果采用0. 25cm层厚来重建 CBCT图像空间，则这一套CBCT图像序列将包含55个DIC0M文件，即55张人体断 层CBCT图像(0. 25cm*55=13. 75<14cm)。
现有技术在实际使用中存在的不足是，OBI设备对获得的CBCT图像只能重建 出不足14cm长度图像片段，使用这些图像片段制定和评估鼻咽癌、食管癌、直肠 癌、宫颈癌等大部分肿瘤放疗计划时，不足以完整显示肿瘤在人体中的位置、形 态等信息，即无法完整显示肿瘤组织的解剖结构，影响放疗计划在CBCT图像上的 制定；另一方面，对于胸部的肿瘤，即使拍摄到了完整的肿瘤组织，通常也无法 完整地再现整个肺的解剖结构，这样将导致医生无法通过CBCT图像来准确地评估 对放疗具有敏感反应的肺组织在接受放疗时实际受到的照射剂量。
发明内容
本发明的目的是克服现有技术存在的不足，提供一种分段锥形束CT图像通过 拼接获得完整解剖图像的方法。
本发明的方法是按以下方式实现的，首先，对患者进行锥形束CT扫描，分段 获取两套以上相互衔接的CBCT图像序列，然后，对两套CBCT图像进行观察、处 理得到一套新的CBCT图像序列；把新的CBCT图像序列进行图像三维重建，得到 一套拓宽了的能够完整展现肿瘤发病区域和周围正常器官的CBCT解剖图像；
分段获取CBCT图像步骤如下
1)设计CBCT分段扫描计划获得模拟定位CT扫描图像并将其传入加速器设 备中的Varian Eclipse治疗计划系统，在计划CT空间上设置左右和前后方向上 坐标完全一致的两个等中心，分别制定CBCT扫描计划P1和P2， OBI单次CBCT扫 描长度为14cm ，因此，头脚方向两等中心间的距离设置在13cm以内，从而保证 两段CBCT图像具有足够的信息冗余；
2)CBCT扫描将上述设计好的CBCT计划传入OBI 4D工作站，按计划PI的等 中心进行摆位，根据扫描部位选定好显像设备及相关参数后进行第一次CBCT扫描， 扫描结束并完成在线配准后，不接受系统给出的自动移床参数，即保持治疗床不 动，然后沿头脚方向手动平移治疗床至计划P2的等中心,移动距离为两等中心的 头脚方向坐标之差，以同样显像参数完成第二次CBCT扫描；
CBCT图像拼接步骤如下
加速器工作站对获取的两套CBCT图像序列分别进行编号，设定最小编号的图像为近脚的一张，最大编号的图像为近头的一张，然后对两套CBCT图像序列分别 进行显像解析；
由于两套CBCT图像是患者在治疗床上不发生体位变动的情况下由OBI先后获 取的，而且两套图像的拍摄范围至少有lcm的重叠，所以，在分别观察两套CBCT 图像时，通过参照骨性标记的解剖结构找到两套图像中所共有的最为相似的一层， 记录下该层分别在各自图像序列中所拥有的序列号；
假定两套CBCT图像分别为down序列和up序列，共有的拼接层面为down序列 的第44层(A)、叩序列的第15层(B)，则需要进行处理的图像序列区间是从 叩序列第16层(C)到up序列最上一层(Z)，如图1所示，处理步骤如下
1) 解析存储CBCT图像的DIC0M文件，获取DICOM文件头中各字段信息；
2) 解析down序列的第A层CBCT图像，读取StudyInstanceUID、 InstanceNumber 和ImagePositionPatient三个字段的值；
3) 依次解析up序列的第C至Z层图像，分别将每层图像对应的DICOM文件 头中的StudylnstanceUID字段值更换成down序列A的相应值；
4) 依次解析up序列的第C至Z层图像，分别将每层图像对应的DICOM文件 头中的InstanceNumber字段值更换成down序列A层InstanceNumber字段值的基 础上依次递增l;
5) 依次解析叩序列的第C至Z层图像，分别更换每层图像对应的DICOM文 件头中ImagePositionPatient字段字节3的值，格式为字节1/字节2/字节3; 因为每套CBCT图像的最下一层图像的Z坐标均为-117. 5,所以更换值通过下面的 公式算得
ImagePositionPatient字节3的值=-117. 5+InstanceNumber值X2. 5 其中，InstanceNumber值已由第四步算得；
6) 保存更改后的每个DICOM文件，并按其InstanceNumber值重新依次命名， 得到一套新的CBCT图像序列。
本发明的优异效果
使用本发明的方法拼接的CBCT图像，用于肿瘤靶区勾画和放射剂量计算，实 现精确计划制定、精确计划评估和精确剂量修正的目的，为放疗科医生及放射物 理师在CBCT图像空间中设计和评估放疗计划创造了更广阔的平台，做到了有的放 矢，大幅度提高了肿瘤的放射治疗效果。
图1.分段CBCT图像拼接过程中关键层面的拼接示意6200810138079.9
说明书第3/6页
像为近脚的一张，最图2.胸部上段CBCT图像冠状截面解剖结构示意图(红线为拼接层)； 图3.胸部下段CBCT图像冠状截面解剖结构示意图(红线为拼接层)； 图4.拼接得到的完整胸部CBCT图像冠状截面解剖结构示意图(红线为拼接 层)；
图5.在一次获取的不完整胸部CBCT图像空间观察肿瘤靶区解剖结构示意图 (冠状面)；
图6.在分段获取并拼接得到的完整胸部CBCT图像空间观察肿瘤靶区解剖结 构示意图(冠状面)；
图7.在一次获取的不完整胸部CBCT图像空间观察肿瘤靶区解剖结构示意图 (矢状面)；
图8.在分段获取并拼接得到的完整胸部CBCT图像空间观察肿瘤耙区解剖结 构示意图(矢状面)；
图9.在分段获取并拼接得到的完整胸部CBCT图像空间制定放疗计划并评估 放射剂量示意图10. CBCT图像分段获取及拼接流程图。
具体实施例方式
分段获取CBCT图像步骤如下
1)设计CBCT分段扫描计划获得模拟定位CT扫描图像并将其传入加速器设 备中的Varian Eclipse治疗计划系统，在计划CT空间上设置左右和前后方向上 坐标完全一致的两个等中心，分别制定CBCT扫描计划P1和P2， OBI单次CBCT扫 描长度为14cm ，因此，头脚方向两等中心间的距离设置在13cm以内，从而保证 两段CBCT图像具有足够的信息冗余；
2)CBCT扫描将上述设计好的CBCT计划传入0BI4D工作站，按计划Pl的等 中心进行摆位，根据扫描部位选定好显像设备及相关参数后进行第一次CBCT扫描， 扫描结束并完成在线配准后，不接受系统给出的自动移床参数，即保持治疗床不 动，然后沿头脚方向手动平移治疗床至计划P2的等中心，移动距离为两等中心的 头脚方向坐标之差，以同样显像参数完成第二次CBCT扫描；
CBCT图像拼接步骤如下 加速器工作站对获取的两套CBCT图像序列分别进行编号，设定最小编号的图像为 近脚的一张，最大编号的图像为近头的一张，然后对两套CBCT图像序列分别进行 显像解析；
由于两套CBCT图像是患者在治疗床上不发生体位变动的情况下由OBI先后获取的，而且两套图像的拍摄范围至少有lcm的重叠，所以，在分别观察两套CBCT 图像时，通过参照骨性标记的解剖结构找到两套图像中所共有的最为相似的一层， 记录下该层分别在各自图像序列中所拥有的序列号；
假定两套CBCT图像分别为down序列和up序列，共有的拼接层面为down序列 的第44层(A) 、 up序列的第15层(B)，则需要进行处理的图像序列区间是从 up序列第16层(C)到up序列最上一层(Z)，如图1所示，处理步骤如下
1) 解析存储CBCT图像的DICOM文件，获取DICOM文件头中各字段信息；
2) 解析down序列的第A层CBCT图像，读取StudylnstanceUID、 InstanceNumber 和ImagePositionPatient三个字段的值；
3) 依次解析up序列的第C至Z层图像，分别将每层图像对应的DIC0M文件 头中的StudylnstanceUID字段值更换成down序列A的相应值；
4) 依次解析up序列的第C至Z层图像，分别将每层图像对应的DICOM文件 头中的InstanceNumber字段值更换成down序列A层InstanceNumber字段值的基 础上依次递增1;
5) 依次解析叩序列的第C至Z层图像，分别更换每层图像对应的DICOM文 件头中ImagePositionPatient字段字节3的值，格式为字节l/字节2/字节3; 因为每套CBCT图像的最下一层图像的Z坐标均为-117. 5，所以更换值通过下面的 公式算得
ImagePositionPatient字节3的值=-117. 5+InstanceNumber值X2. 5 其中，InstanceNumber值已由第四步算得；
6) 保存更改后的每个DICOM文件，并按其InstanceNumber值重新依次命名， 得到一套新的CBCT图像序列。
发明方法准确性及可行性验证
1.对拼接得到的CBCT图像几何学完整性进行体模验证对TOPSLANE QIY00A 型胸段体模进行分段CBCT扫描，完成拼接后导入Varian Eclipse治疗计划系统， 对体模CBCT图像进行三维重建；在拼接后CBCT断层图像上逐层勾画体模及其模 拟的双肺组织外轮廓，分别计算体模及双肺的体积，并将该结果与计划CT下重建 并测得的轮廓体积进行比较。结果显示，拼接后CBCT图像测定体模外轮廓体积与 计划CT测量结果相差0. 26% (28. 34 cm3)，拼接CBCT图像测得体模两侧肺体积 分别与计划CT测得结果相差1. 87% (12.82 cm3)和1. 47% (10.07 cm3);
2.对拼接CBCT图像进行包容度在体应用研究选择10例胸部肿瘤患者，从 肺容量包容度角度对CBCT拼接图像进行验证。于患者第一次实施影像引导放射治疗前，在自由呼吸状态下，通过水平移床获取CBCT分段扫描图像，离线后将分段 扫描获得的两套CBCT图像进行拼接并导入Varian Eclipse进行三维重建。在拼 接后CBCT断层图像上逐层勾画人体肺组织轮廓，评估肺容量，并将其与计划CT 重建后测得的相同参数进行比较。10例患者拼接后CBCT图像与计划CT图像肺容 量比对分析见表l。结果显示，全肺容量平均相差1.97%±0.42%
(64.53士26.07cm3),左肺容量平均相差2. 30%±0. 78% (33. 32± 17. 03cm3)，右 肺容量平均相差1.75%±0. 2% (31. 21 ±12. 51 cm3);
验证结果显示，无论体模还是人体，利用本发明的方法很容易实现分段CBCT 的拼接，而且拼接得到的新的CBCT图像序列经三维重建后，精确地再现了体模和 人体肺组织的解剖结构，证明拼接过程无图像信息丢失，几何学精度高，该发明 方法完全可以用于肿瘤放疗临床应用。
10例患者拼接后CBCT图像与计划CT图像肺容量比对分析(体积单位为cm3)
如表l.所示。
表i
病例VcT扁E!tVgbct-左bAV左肺VcT-右肺VcBCT-右據△V右肺厶v全肺
12187.782230.9143.13(1.97%>2243.672289.6746.00(2.05%)89.13(2.01%>
21939.371971.2531.88(1.64%)1993.262031.5338.27(1.92%)70.15(1.78%)
3945.80977.8332.03(3.39%)1331.331353.7022.37(1.68%)54.40(2.39%)
4707.06701.35-5.71(0.81%)1657.141685.9728.83(1.740/。)23.12(0.98%)
51291.21334.4343.23(3.35%)1674.871705.3530.48(1.82%)73.71(2.49%)
6707.13726.5119.38(2.74%)794.36805.8811.52(1.45%)30.90(2.06%)
2339.042392.8553.81(2.30%)2898.772941.3842.61(1.47%)96.42(1.84%)
81886.811922.1035.29(1.87%)1676.41704.7328.33(1.69%)63.62(1.79%)
91178.871209.48'30.61(2.60%)917.19932.7815.59(1.70%>46.20(2.20%)
102083.662133.1949.53(2.38%)2429.312477.4148.10(1.98%)97.63(2.16%)
平1526.671559.9933.32±17.031761.631792.8431.21±12.5164.53±26.07
均±629.33±642.67(2.30%±0.78%)±656.97±668.75(1.75%±0.20%)(1.97%±0.42%)
权利要求
1、分段锥形束CT图像通过拼接获得完整解剖图像的方法，其特征在于对患者进行锥形束CT扫描，分段获取两套以上相互衔接的CBCT图像序列，然后，对每两套相邻CBCT图像进行观察、处理得到一套新的CBCT图像序列；把新的CBCT图像序列进行图像三维重建，得到一套拓宽了的能够完整展现肿瘤发病区域和周围正常器官的CBCT解剖图像；分段获取CBCT图像步骤如下1)设计CBCT分段扫描计划获得模拟定位CT扫描图像并将其传入加速器设备中的Varian Eclipse治疗计划系统，在计划CT空间上设置左右和前后方向上坐标完全一致的两个等中心，分别制定CBCT扫描计划P1和P2，OBI单次CBCT扫描长度为14cm，因此，头脚方向两等中心间的距离设置在13cm以内，从而保证两段CBCT图像具有足够的信息冗余；2)CBCT扫描将上述设计好的CBCT计划传入OBI 4D工作站，按计划P1的等中心对扫描对象进行摆位，根据扫描部位选定好显像设备及相关参数后进行第一次CBCT扫描，扫描结束并完成在线配准后，不接受系统给出的自动移床参数，即保持治疗床不动，然后沿头脚方向手动平移治疗床至计划P2的等中心，移动距离为两等中心的头脚方向坐标之差，以同样显像参数完成第二次CBCT扫描；CBCT图像拼接步骤如下加速器工作站对获取的两套CBCT图像序列分别进行编号，设定最小编号的图像为近脚的一张，最大编号的图像为近头的一张，然后对两套CBCT图像序列分别进行显像解析；. 由于两套CBCT图像是患者在治疗床上不发生体位变动的情况下由OBI先后获取的，而且两套图像的拍摄范围至少有1cm的重叠，所以，在分别观察两套CBCT图像时，通过参照骨性标记的解剖结构找到两套图像中所共有的最为相似的一层，记录下该层分别在各自图像序列中所拥有的序列号；假定两套CBCT图像分别为down序列和up序列，共有的拼接层面为down序列的第44层(A)、up序列的第15层(B)，则需要进行处理的图像序列区间是从up序列第16层(C)到up序列最上一层(Z)，如图1所示，处理步骤如下1)解析存储CBCT图像的DICOM文件，获取DICOM文件头中各字段信息；2)解析down序列的第A层CBCT图像，读取StudyInstanceUID、InstanceNumber和ImagePositionPat ient三个字段的值；3)依次解析up序列的第C至Z层图像，分别将每层图像对应的DICOM文件头中的StudyInstanceUID字段值更换成down序列A的相应值；4)依次解析up序列的第C至Z层图像，分别将每层图像对应的DICOM文件头中的InstanceNumber字段值更换成down序列A层InstanceNumber字段值的基础上依次递增1；5)依次解析up序列的第C至Z层图像，分别更换每层图像对应的DICOM文件头中ImagePositionPatient字段字节3的值，格式为字节1/字节2/字节3；因为每套CBCT图像的最下一层图像的Z坐标均为-117.5，所以更换值通过下面的公式算得ImagePositionPatient字节3的值＝-117.5+InstanceNumber值×2.5其中，InstanceNumber值已由第四步算得；6)保存更改后的每个DICOM文件，并按其InstanceNumber值重新依次命名，得到一套新的CBCT图像序列。
2、根据权利要求1所述的方法，其特征在于，图像三维重建拓宽后的空间由 公式14XX-(X -l)算出，公式中的14是一段CBCT图的宽度长度，单位是cm， X 表示CBCT图像的段数，1表示两张CBCT图像相互重叠的宽度，因此，CBCT图像 扩展的宽度取决于分段获取的CBCT图像的段数。
全文摘要
本发明提供一种分段锥形束CT图像通过拼接获得完整解剖图像的方法，该方法是在计划CT图像空间分别制定两段CBCT扫描计划；在获得一套CBCT图像后，手动平移治疗床，再次行CBCT拍摄以获取相邻分段的CBCT图像；将两段图像分别进行显像解析，再经三维重建后分别与计划CT进行图像融合；通过对融合后骨性标志的分析来选取拼接层面；将两套CBCT图像依次通过计算机进行拼接处理，获得相邻部位分段扫描所得CBCT图像的拼接图像，扩大了肿瘤靶区及周围正常器官的观察范围，保证了以CBCT制定治疗计划的可行性及评价剂量体积参数的完整性，有广泛的临床应用前景。
文档编号G06T17/00GK101427924SQ200810138079
公开日2009年5月13日 申请日期2008年7月16日 优先权日2008年7月16日
发明者于金明, 敏 刘, 刘同海, 洁 卢, 涛 孙, 勇 尹, 张桂芳, 健 朱, 曈 白, 陈进琥 申请人:山东省肿瘤医院