[发明专利]一种用于观察子宫肿瘤放疗总剂量的分布图像的生成方法

摘要

本发明涉及X射线辐射的测量，具体涉及一种用于观察子宫肿瘤放疗总剂量的分布图像的生成方法，该方法由以下步骤组成：先根据图像中像素点CT值的大小，分割出含施源器HDR CT图像中的施源器区域，并进行二值化处理得到掩膜图像；然后，利用掩膜图像对HDRCT图像进行收缩变形和线性插值处理，再以处理后的图像为浮动图像，以不含施源器的盆腔IMRT CT图像为参考图像进行变形配准；最后，用得到的配准变形场对去除施源器区域剂量的HDR剂量分布图像进行变形，该变形结果与不包含施源器的IMRT剂量分布图像进行叠加后再与不包含施源器的IMRT CT图像进行图像融合。本方法所生成的图像可准确地评价子宫肿瘤病人受照的总剂量。

主权项

一种用于观察子宫肿瘤放疗总剂量的分布图像的生成方法，该方法由以下步骤组成：（1）读取含施源器的盆腔HDR CT图像、盆腔HDR剂量分布图像、不含施源器的盆腔IMRT CT图像和盆腔IMRT剂量分布图像；（2）根据图像中像素点CT值大小的不同，采用区域生长的分割法分割出含施源器盆腔HDR CT图像中的施源器区域，并赋予所分割出的施源器区域像素值为0，得到施源器区域像素值为0的盆腔CT图像；接着，将施源器区域像素值为0的盆腔CT图像中的施源器区域用1表示，其他区域用0表示，得到掩膜图像；同时，将所述的盆腔HDR剂量分布图像中对应的施源器区域的像素值赋予为0，得到施源器区域像素值为0的盆腔HDR剂量分布图像；（3）利用变形场将掩膜图像进行收缩变形，该收缩变形的具体步骤如下：在三维直角坐标系中，令收缩变形力F在x和y方向上的收缩变形力为掩膜图像在x和y方向上的梯度，收缩变形力F在z方向上的收缩变形力为0，然后求解下式（Ⅰ）所示的Navier‑Stokes方程，得到收缩的变形场u和收缩变形的掩膜图像，▽2v+▽(▽·v)+F=0     （Ⅰ）上式（Ⅰ）中，▽为梯度，▽2为拉普拉斯算子，v为变形场的速度，F为收缩变形力；（4）用步骤（3）所得的变形场u分别对施源器区域像素值为0的盆腔HDR CT图像和施源器区域像素值为0的盆腔HDR剂量分布图像进行收缩形变，得到收缩变形的盆腔HDRCT图像和收缩变形的盆腔HDR剂量分布图像；然后分别搜索出收缩变形的盆腔HDR CT图像和收缩变形的盆腔HDR剂量分布图像中对应于收缩变形的掩膜图像中像素值不等于0的点，并分别对所述不等于0的点进行线性插值，得到去除施源器后的盆腔HDR CT图像和去除施源器区域剂量的盆腔HDR剂量分布图像；（5）以去除施源器后的盆腔HDR CT图像为浮动图像，以盆腔IMRT CT图像为参考图像进行变形配准，同时得到配准变形场r和变形配准后的盆腔HDR CT图像；（6）用步骤（5）得到的配准变形场r对去除施源器区域剂量的盆腔HDR剂量分布图像进行变形，得到变形配准后的盆腔HDR剂量分布图像；（7）先将变形配准后的盆腔HDR剂量分布图像与不包含施源器的盆腔IMRT剂量分布图像进行叠加，再把叠加结果与不包含施源器的盆腔IMRT CT图像进行图像融合，得到具有人体解剖结构的辐射总剂量的分布图像。