1. 引言
能谱CT成像技术已成为CT高级临床应用的发展方向,能谱CT辐射源特性与常规CT不同,还没有建立精准研究方法和评价标准及评价体系。目前能谱CT的辐射剂量与影像质量评价以及最优化匹配尚没有系统性的研究,国内外研究仍处于起步阶段。电离辐射所致人体和器官剂量可直接表现射线对人体可能造成的损伤,所以对辐射剂量进行准确评估具有重要意义 [1] 。而评估受检者辐射剂量是至关重要的,由于难以通过探测器直接测量器官剂量,因此借助计算机人体模型利用蒙特卡罗方法估算剂量成为准确评估剂量的重要手段。
儿童受检者作为辐射敏感特殊人群正处于生长发育的关键时期,对辐射具有较高的敏感性 [2] 。国际放射防护委员会ICRP102号报告指出,全球CT检查的应用频率在过去二十年迅速增长,随着CT检查在儿童疾病诊断中应用更广泛,其中儿童约占CT检查人群的10%以上 [3] 。目前新冠疫情进入疫情常态化阶段,国内外感染新冠的儿童人数逐步增多,放射学CT扫描技术成为诊断新冠确诊病例与无症状感染者及评价新冠患者病情的重要手段 [4] ,因此儿童将受到更多的医疗照射可能性。CT检查引起的照射会增加发生潜在疾病的风险,尤其是正处于生长发育期的儿童,其对电离辐射也更加敏感,同样的照射剂量下儿童存在更高的患病率与致死率 [5] 。评价受检者辐射剂量所采用的物理表征量是器官剂量与有效剂量,而受检者辐射剂量多采用蒙特卡罗方法结合参考人曲面模型进行仿真模拟估算。
目前,辐射防护用儿童人体模型已经有很大进展。国外已经建立了符合自己本国人群生理参数的儿童参考人曲面模型,美国Choonsik等基于CT影像数据构建包含1岁、5岁、10岁和15岁儿童的NURBS表面模型库 [6] ,UF family。日本Tomoaki在已有的成年人模型基础上,通过变形技术建立身高和体重符合标准日本3岁、5岁和7岁儿童的模型 [7] 。国际辐射防护委员会(ICRP)第103工作组完成了儿科网格型参考计算模型(MRCP)的开发,包括10个模型(新生儿、1岁、5岁、10岁和15岁男性和女性) [8] ,解决了ICRP出版物143的儿科参考计算模型的局限性,这些模型采用体素格式。儿科MRCP是通过将体素型参考模型转换为高质量网格格式来构建的,在小而复杂的器官和组织(如骨骼、淋巴结)的详细解剖结构中有显著增强。但是ICRP报告中参考人的数据是基于高加索人种得到的,不符合中国儿童的生理参数,并不适用于国人,因此有必要建立符合中国儿童体格特征与器官位置的体模用于评估辐射剂量。
国内虽已有儿童参考人曲面模型,但还未有系统的儿童参考人曲面模型的研究。中国科学技术大学将已建立的符合ICRP参考人数据的成年男女面元模型(RIP-AM、PIR-AF)进行分段缩放和调整,建立了符合GBZ/T200《辐射防护用参考人》标准的中国儿童参考人面元模型 [9] 。清华大学根据1岁儿童物理体模(Model 704-D)的CT扫描图像建立了符合中国儿童参考人特征的体素模型(CPP01) [10] ,模型CPP01的器官和骨骼是通过缩放成年人得到的。由于儿童的骨骼发育水平、部分器官形态和相对位置与成人差距较大,因此儿童体模CPP01不能准确表现中国儿童的解剖学特性。又基于中国儿童真实CT扫描图像建立了10岁中国儿童参考人面元体模 [11] ,该模型采用多边形网格模型与NURBS曲面模型相结合的方法构建了85个器官,经过优化和调整,各组织和器官质量与GBZ/T200.2-2007和亚洲参考人的参考值差异小于1%,各器官位置和形态符合中国儿童解剖学特征。接着其工作组开发了中国参考模型系列,包括3个月(CRC3m)、1岁(CRC01)、5岁(CRCO5)、10岁(CRC10)、15岁男性(CRCM15)和15岁女性(CRCF15)模型,图1为其开发的网格型中国儿童参考模型系列 [12] 。
Figure 1. Schematic diagram of the child’s reference human meta-model
图1. 儿童参考人面元模型示意图
现有的儿童参考人模型均是基于较早年代的儿童生理参数建立的,而当代同一年龄段儿童生理参数由于营养结构变化等因素而发生变化 [13] ,根据调研的儿童体格数据来看,近几年儿童的体格数据与GBZ/T200.2-2007的参考值对比有小幅度增长,因此现有的儿童参考人模型已经不能准确反映我国当代儿童的身体特征。本文调研的儿童体格数据来自上海市第一妇婴保健院,其中0.25岁、1岁、5岁男孩、5岁女孩、10岁男孩和10岁女孩数据各30例,0.25岁和1岁儿童未区分性别。图2和图3为调研的儿童身高体重数据和GBZ/T200.2-2007参考值的对比。
Figure 2. Children’s weight and GBZ/T200.2-2007 reference value comparison
图2. 儿童体重和GBZ/T200.2-2007参考值对比
Figure 3. Children’s height and GBZ/T200.2-2007 reference value comparison
图3. 儿童身高和GBZ/T200.2-2007参考值对比
目前中国辐射防护用参考人中10岁以前(包括10岁)的身高体重男女并无差异,但图中体现出当代男女儿童体格数据上较参考值有所增长,并且10岁男童的体重数据较中国辐射防护用参考人有明显增加,因此需要重新构建中国儿童参考人模型。本文基于近几年的中国儿童CT影像构建一套不同年龄的中国儿童参考人曲面模型,包括新生儿0.25岁、1岁、5岁和10岁的男性和女性人体模型,后续可用于能谱CT所致儿童受检者的辐射剂量评估。
2. 体模构建
数字化虚拟参考人模型在为各类人群制定辐射防护方案的过程中,具有不可替代的作用。从发展历史和模型的性质来看,分为数学模型、体素模型和曲面模型。与使用圆柱体、封闭二次曲面建模的数学模型相比,曲面模型可真实反映人体器官形态和相对位置,在剂量评估中更高的精度。体素模型有局限性,任何实际尺寸的体素都无法再现精细结构,其中一些器官对于辐射防护非常重要,如眼睛的晶状体等,这些组织吸收剂量的计算必须依赖于专门的部分人体模型。相比于体素模型,曲面模型能很好地保证人体内部细微结构地高精度和部分结构由生理运动导致地可变性 [14] ,易于通过灵活的变形操作反映人体在辐射场中的不同姿态和不同体格人的身体和器官形态,同时曲面模型也可描述呼吸道壁等使用体素模型难以描述的组织,具有更高的精确性。
本文构建了一系列中国儿童参考人曲面模型,模型的器官质量和位置与中国参考儿童的解剖学特征一致。体模构建过程为CT图像分割、器官优化、NURBS曲面重构、模型装配及微调,流程如图4所示。模型基于逆向工程的三维实体建模技术,对初始重建的器官进行优化并构建NURBS曲面模型,对所有的器官和组织的曲面模型进行调整和装配,从而得到中国儿童曲面模型,为后续能谱CT所致儿童受检者的辐射剂量估算做准备。
Figure 4. The construction process of Chinese children’s curved surface model
图4. 中国儿童曲面模型的构建流程
本研究拟建立一个符合我国当代儿童身体特征的参考人曲面模型组,包含0.25岁、1岁、5岁男性、5岁女性、10岁男性和10岁女性的曲面体模。ICRP103工作组和清华实验室完成的儿科网格型参考计算模型中包括15岁的男性和女性,目前15岁儿童在外观上已接近成年人,故不在本研究对象范围。儿童5岁之前,性别器官不太明显,身体各器官的组成相似,因此,对于5岁之前的年龄组,本文儿童参考人曲面模型的构建与性别无关。
2.1. CT扫描图像与分割
本文所用CT影像为同济大学附属第一妇婴保健院提供的不同年龄儿童不同体段的CT扫描图像。本文建立模型基础CT图像由不同儿童的头部、颈部、胸腹部、盆部和足部扫描图像分割拼接而成。
为使得本文建立的儿童体模具有代表性,在选取CT扫描图像时,所对应的儿童受检者身高、体重与中国辐射防护用参考人中提供的各年龄儿童参考人的数值差异并没有太大。同时,体模建立过程中应对身体轮廓、器官形态和质量进行一定的调整,通过调整优化使体模的体格参数及器官质量尽量符合中国儿童参考人的数值。
建立儿童体模首先进行图像分割,将CT扫描图像导入MIMICS软件,对CT影像文件进行阈值分割、蒙版编辑、区域生长、三维重建等操作。通过使用中值滤波操作以增加CT图像中器官的对比度,这提高了CT图像的质量,便于后续的器官手动分割。由于CT扫描对软组织的分辨能力弱,导致影像中器官边缘不明显,因此分割腹腔等部位器官时需要一定的解剖学知识和经验,本工作在专业医师指导下对器官进行准确分割。图5为腹腔CT图像的分割。将分割得到的器官轮廓利用MIMICS的三维重建功能导出STL格式的几何模型。
Figure 5. Segmentation of CT images of abdominal cavity
图5. 腹腔CT图像的分割
2.2. 器官优化
本文中使用了真实的中国儿童的CT图像并对其进行了三维重建,得到了不同器官和组织的几何模型。对于新生儿以及1岁儿童一些分割不清晰的器官,为了使模型更加精准,减少其对蒙特卡罗模拟计算的影响,主要通过Geomagic Studio 2015软件中多边形处理模块对不同器官和组织的几何模型进行优化。将前面所建立的三维几何模型导入Geomagic Studio 2015软件中,首先删除模型表面的钉状物和光滑三角网格;然后修复模型内部相交的区域、消除重叠的三角形;最后进行简化三角面片数量、松弛多边形以及创建流行等操作。在构造参数化曲面模块下,首先绘制器官模型的轮廓线,并对其进行松弛、收缩、合并、细分或者延伸等处理改变器官的轮廓使其更加接近真实器官;其次,探测器官模型的曲率线,并对其进行升级或约束处理;然后移动和修复曲面片;最后,对曲面片进行格栅构造、编辑以及简化处理,即可获得器官的NURBS曲面模型。图6所示为肺的NURBS曲面构建过程。
(a) 优化前 (b) 曲面构造 (c) 格栅构造 (d) NURBS曲面
Figure 6. NURBS surface construction process of lung
图6. 肺的NURBS曲面构建过程
2.3. 器官装配
本文使用Rhinoceros 6.0软件来装配和整合这些模型,装配时保证人体处于正直站立状态,且使体模身高与参考人数据相同。将多个器官和组织的NURBS模型分别导入犀牛软件,并将其置于不同的图层中。皮肤、骨骼以及其他未建立的微小器官和组织直接利用犀牛软件强大的NURBS建模功能进行建模,并通过各种缩放工具进行调整,使其形状、质量均符合中国辐射防护用参考人的参数。肌肉和脂肪在本研究中不是重点研究的对象,并且CT值较为近似,所以将其统一定义为其余组织。将以上所有器官、组织的相对位置进行调整和装配,并确保各个器官和组织之间没有重叠,进而得到中国儿童曲面模型,以OBJ格式保存导出。
3. 曲面体模构建结果
依照上述方法建立了中国儿童参考人的曲面模型。儿童男性和女性体模除生殖器官外,其他器官的差别较小。本文建立的所有儿童体模包含睾丸、卵巢、肺、胃、膀胱、胸腺、肝、食管、甲状腺、脑、心脏、肾、脾、唾液腺、胆囊、小肠、胰、眼球、晶状体、肾上腺、皮肤等器官或组织在内的人体脏器,同年龄的男女共用一套器官和身体外表面,5岁和10岁儿童体模只在涉及生殖器官处有不同。建立的中国儿童参考人体模的器官质量与GBZ/T200.2-2007确定的参考人器官质量相近,体模器官质量通过体积乘以该器官的密度得到。表1与表2列出了体模主要器官和GBZ/T200.2-2007确定的参考人器官质量对比。
Table 1. Body model and GBZ/T200.2-2007 organ quality comparison
表1. 体模与GBZ/T200.2-2007器官质量的对比
Table 2. Body model and GBZ/T200.2-2007 organ quality comparison
表2. 体模与GBZ/T200.2-2007器官质量的对比
本文所建立体模的大部分器官和组织数据与中国辐射防护用参考人数值和亚洲参考人标准相对差异在3%以内。在器官构建过程中,一些器官的体积和质量相较GBZ/T200.2-2007规定有一定程度增加,以体现当代儿童身体素质所导致的体格变化。
4. 结论
本文基于真实的中国儿童CT扫描图像,建立了体格特征和器官质量符合当代中国儿童的参考人曲面模型。该体模组各年龄的曲面模型的身高、体重等体格数据与选取的儿童受检者一致,将所有器官和组织调整到没有重叠,在形态和相对位置上均符合中国儿童的解剖学特性。该体模的所有器官和组织的数据较中国辐射防护用特殊人群参考人都有不同程度的增加,因此所建的中国儿童体模组更能准确表现当代儿童的体格特征。
儿童解剖特征真实准确的体模可以更精确评估受检者辐射剂量。由于儿童尚处于生长发育时期,故儿童的部分器官在形态、相对大小、相对位置上较成年人有很大差别,两者间的差异以腹腔器官尤为突出。儿童肝脏占据腹腔空间的比例明显大于成年人,肝脏附近脏器大多紧密相连,本文通过直接分割儿童的CT图像建立的模型,在调整时改动较小,可以保持儿童器官的解剖学正确性。本文建立的NURBS曲面儿童模型可通过NURBS变形功能实现对体模器官的平滑微调,且调整过程在一定程度上反映了器官形变特性,为今后调整体模以适应不同的体格和姿态提供了有效的方法。
5. 讨论
本文的实验数据来自于上海市第一妇婴保健院,因此儿童体格数据多来自上海小孩,考虑儿童的身高生长和所在行政区的发展水平、人均可支配收入水平有关,因此在数据的选取上存在局限。
本文建立的中国儿童参考人曲面模型为能谱CT下儿童受检者辐射剂量的评估提供了重要工具,后续利用蒙特卡罗方法,可实现能谱CT扫描所致特殊人群受检者的器官剂量与有效剂量的快速估算。结合NURBS曲面模型具备精确的解剖结构和结构姿势的可变性,可实现对人体器官几何模型结构和轮廓的调整,使其满足中国辐射防护用参考人的参数标准。因此可以通过变形调整得到不同体格、不同姿势的体模,可以更加真实地模拟不同姿势下,不同年龄的儿童接受照射时的剂量。
NOTES
*第一作者。
#通讯作者。