姜立朋 張懷文 谷珊珊

近距離治療中放射源位置檢驗的數(shù)學(xué)形態(tài)學(xué)方法

姜立朋1張懷文2谷珊珊3

1（遼寧醫(yī)學(xué)院附屬第一醫(yī)院 放療科 錦州 121001）
2（江西省腫瘤醫(yī)院 放療技術(shù)科 南昌 330029）
3（解放軍總醫(yī)院 放療科 北京 100853）

利用數(shù)學(xué)形態(tài)學(xué)算法計算EBT3膠片拍攝的近距離治療中放射源的位置準確性。目前驗證192Ir放射源到位精度的方法主要包括膠片、Fricke放射量測定器、平板探測器、簡單的直尺結(jié)合計量儀方法等，數(shù)學(xué)形態(tài)學(xué)算法自動計算曝光膠片放射源的位置未見報道。本研究用EBT3免沖洗膠片分別獲得不同步進長度和已知駐留位置誤差的影像。采用Contourlet變換對膠片圖像進行去噪，求得放射源步進長度的誤差小于±0.2 mm，源駐留位置的誤差在±0.5 mm以內(nèi)，符合臨床要求。結(jié)果表明本方法簡單，且結(jié)果可靠。

192Ir，近距離放射治療，膠片，數(shù)學(xué)形態(tài)學(xué)，質(zhì)量保證

近距離治療中放射源在施源器中通過優(yōu)化駐留位置和時間而形成特定的劑量分布。源位置的準確性直接影響吸收劑量分布[1]。放射源的位置準確性是近距離治療重要的質(zhì)量保證(Quality Assurance, QA)項目[2?3]。膠片是驗證192Ir放射源到位精度的常用方法[4?5]，膠片曝光后放射源位置的讀取一般用直尺測量[6]，存在一定的人為誤差。采用數(shù)學(xué)方法自動計算曝光膠片放射源的位置具有重要的應(yīng)用價值，但還未見相關(guān)報道。本文利用數(shù)學(xué)形態(tài)學(xué)算法估算EBT3膠片上放射源的位置坐標，分析放射源的位置準確性。

1 材料與方法

1.1 近距離治療機

瑞典醫(yī)科達核通(Nucletron)公司生產(chǎn)的micro-Selectron HDR (High Dose-rate)高劑量率近距離治療機由治療計劃系統(tǒng)、控制系統(tǒng)和放射源傳送系統(tǒng)組成。利用計算機控制步進馬達驅(qū)動高活性192Ir放射源，配合優(yōu)化的治療計劃實施治療。192Ir放射源的活性長度為3.5 mm，源的最大活度為0.37TBq。治療過程中，放射源以快于50 cm·s?1的速度移動到治療位置，每個駐留位置的照射時間有兩個獨立計時器控制。

1.2 EBT3免沖洗膠片

GAFCHROMIC EBT3是一種基于輻射交聯(lián)效應(yīng)使塑料變色反映劑量變化的特種透明膠片，厚度為0.28 mm，具有免沖洗、防水等特點。表面有反牛頓環(huán)的涂層，可以避免平板掃描儀產(chǎn)生的牛頓環(huán)偽影。劑量在1?800 cGy，能量響應(yīng)和組織等效性較好，均勻一致性優(yōu)于1.5%，空間分辨率優(yōu)于0.1 mm[7?8]。

1.3 膠片測量

將EBT3免沖洗膠片平放粘貼在RW3固體水平板上，塑料施源管拉直平放在EBT3膠片上。以參考位置為中心，選擇步進長度分別為40 mm、30mm、20 mm、10 mm、5 mm、2.5 mm制定計劃，每個步進長度設(shè)置5個放射源駐留位置，駐留時間為10 s，分別得到5列不同步進長度下192Ir源曝光的圖像。為驗證數(shù)學(xué)形態(tài)學(xué)算法的準確性和精度，在?40 mm、0 mm、+40 mm位置分別設(shè)置已知的位置誤差，第2組為?2 mm、+2 mm、?1 mm，第3組為+2 mm、?2 mm、+1 mm，加上參考位置圖像，共得到三組曝光的影像。

1.4 數(shù)學(xué)形態(tài)學(xué)算法

數(shù)學(xué)形態(tài)學(xué)(Morphology)是以形態(tài)為基礎(chǔ)對圖像進行分析的數(shù)學(xué)工具，基本原理是用具有一定形態(tài)的結(jié)構(gòu)元素去度量和提取圖像中的對應(yīng)形狀，對圖像進行分析和識別。數(shù)學(xué)形態(tài)學(xué)算法主要通過選擇相應(yīng)的結(jié)構(gòu)元素采用膨脹、腐蝕、開啟、閉合4種基本運算的組合來處理圖像，對基于形狀細節(jié)進行圖像處理提供了強有力的手段，具有直觀上的簡明性和數(shù)學(xué)上的嚴謹性[8]。因為獲取的掃描膠片圖像存在噪聲，對于后續(xù)斑點質(zhì)心的提取干擾較大，采用基于Contourlet變換的圖像去噪算法對膠片圖像進行預(yù)處理。去噪目的就是從含噪聲圖像y中恢復(fù)原始圖像，保持圖像x的特征，優(yōu)化均方差。去噪一般采用模型：y=x+δe。其中，x是期望圖像，y是原始膠片圖像，e是噪聲，δ是噪聲方差?；贑ontourlet變換的圖像去噪算法具體實現(xiàn)方法為：

(1) 確定Contourlet分解的層次K，對含噪圖像進行Contourlet變換，得到低頻系數(shù)a0和高頻系數(shù)d0, d1, …, dK?1。

(2) 在Contourlet變化域設(shè)定閾值對系數(shù)進行處理，得到新的Contourlet系數(shù)。其中閾值設(shè)定為：

對于已經(jīng)預(yù)處理的膠片圖像，采用二值化的質(zhì)心計算方法完成質(zhì)心提取。設(shè)m×n維二值圖像(i, j)，其中目標部分為A，背景部分為B。即：

取圖像中任一點作為質(zhì)心搜索過程中的初始點，也就是動態(tài)質(zhì)心的初始點。搜索過程將由此點出發(fā)直到目標質(zhì)心。令(K, L)表示動態(tài)質(zhì)心，函數(shù)D(i-K,j-L)為圖像各點到(K, L)的歐氏距離：

D(K, L)表示圖像與模塊在動態(tài)質(zhì)心(K, L)點的相關(guān)值。搜索質(zhì)心就是內(nèi)積H值沿某個軌跡迅速下降的極小化過程。當搜索到達真實質(zhì)心時，H值達到最小。從式(4)可以看出，計算(K, L)點H值時，模板中心與(K, L)重合。計算(K, L) 8個鄰點的H值，其中最小的H值對應(yīng)的鄰點為下一個動態(tài)質(zhì)心。

利用MATLAB 2012b編程，運用形態(tài)學(xué)算法求膠片曝光黑斑的質(zhì)心，步驟如下：(1) 對圖像通過灰度化和反色后閾值選擇得到光斑特征區(qū)域；(2)基于Contourlet變換模糊去噪，消除熱噪聲以及像素不均勻產(chǎn)生的噪聲；(3) 再次進行閾值選擇，得到更清晰的光斑區(qū)域；(4) 形態(tài)學(xué)處理，得到連通域的規(guī)則形狀圖形；(5) 采用二值化質(zhì)心計算方法提取光斑的質(zhì)心作為源的駐留位置，與源的預(yù)設(shè)駐留位置比較，從而得到源的駐留位置偏差。

2 結(jié)果

2.1 步進長度的檢測

放射源步進長度分別為40 mm、30 mm、20 mm、10 mm、5 mm、2.5 mm，連續(xù)步進5個駐留位置得到192Ir源曝光的膠片見圖1。膠片掃描后用數(shù)學(xué)形態(tài)學(xué)算法處理，得到的放射源位置圖像見圖2。由圖1和2，當步進長度小于5 mm時，因黑斑黑度重迭，數(shù)字形態(tài)學(xué)算法不能分辨單個放射源的位置。

圖1 源步進長度的檢測(a) 曝光的膠片，(b) 計算的放射源位置Fig.1 Detection of source step length. (a) Exposed film, (b) Position of radiation source

每個步進長度設(shè)置連續(xù)5個駐留位置，5個駐留位置之間的4個步進長度的參考值和計算值的結(jié)果見表1。

表1 數(shù)學(xué)形態(tài)學(xué)算法計算的192Ir源步進長度Table 1 Source step length calculated using mathematical morphology.

從表1可以看出，利用數(shù)學(xué)形態(tài)學(xué)算法得到的源位置誤差在±0.2 mm以內(nèi)，準確性比較好。

2.2 放射源位置誤差檢測的準確性

設(shè)置了放射源位置誤差，曝光的膠片和數(shù)學(xué)形態(tài)學(xué)方法得到的源質(zhì)心坐標如圖2所示。第1組曝光黑斑為沒有誤差的參考192Ir源曝光位置，第2組和第3組曝光黑斑分別為人為設(shè)置了誤差的192Ir源曝光位置。

從表2可以看出，數(shù)學(xué)形態(tài)學(xué)算法可以檢測到放射源的位置誤差小于±0.3 mm，利用數(shù)學(xué)形態(tài)學(xué)算法可以自動分析膠片影像的192Ir源位置準確性。

圖2 源位置誤差的檢測(a) 膠片，(b) 源位置Fig.2 Detection of source position error. (a) Film, (b) Source position

表2 數(shù)學(xué)形態(tài)學(xué)算法計算的192Ir源駐留位置誤差Table 2 Error of source dwell position calculated using mathematical morphology.

3 結(jié)語

近距離放射治療中放射源的位置和活度以及照射時間是影響吸收劑量分布的主要因素[9]。近距離治療的放射源緊貼靶區(qū)和人體組織，單次劑量較大，照射時間長，需要嚴格的質(zhì)量保證與控制措施，確?；颊叩陌踩玔10?12]。放射源的位置準確性是質(zhì)量保證的重要內(nèi)容，膠片曝光是放射源位置準確性檢測最直接和準確的方法，但利用直尺等工具的目測方法誤差較大。

數(shù)學(xué)形態(tài)學(xué)算法適用于處理192Ir放射源曝光的膠片，可以檢測放射源的位置偏差，誤差在±0.5 mm以內(nèi)，能滿足臨床需求，但放射源步進的距離不能太小。數(shù)學(xué)形態(tài)學(xué)算法分析膠片影像信息確定放射源的位置，方法簡單易行和結(jié)果準確，為近距離放射治療中放射源位置準確性的驗證提供了一種新方法。

1 曹振, 趙華立, 謝芹, 等.125I-103Pd復(fù)合近距離治療源劑量參數(shù)的蒙特卡羅確定[J]. 核技術(shù), 2013, 36(1): 010203

CAO Zhen, ZHAO Huali, XIE Qin, et al. Determination of dosimetric characteristics of125I-103Pd brachytherpy source with Monte-Carlo method[J]. Nuclear Techniques, 2013, 36(1): 010203

2 中華人民共和國衛(wèi)生行業(yè)標準. WS 262-2006, 后裝γ源治療的患者防護與質(zhì)量控制檢測規(guī)范[S]. 北京: 人民衛(wèi)生出版社, 2007

Health industry standard of the People's Republic of China. WS 262-2006, Specification for radiological protection of patients and quality control in γ-ray source afterloading brachytherapy[S]. Beijing: People's Health Publishing House, 2007

3 中華人民共和國國家職業(yè)衛(wèi)生標準. GBZ 121-2002, 后裝γ源近距離治療衛(wèi)生防護標準[S]. 北京: 人民衛(wèi)生出版社, 2002

National occupational health standards of the People's Republic of China. GBZ 121-2002, Radiologic protection standards for γ-ray afterloading brachytherapy[S]. Beijing: People's Health Publishing House, 2002

4 余耘, 馮寧遠. 微型銥源空間劑量角分布的測定[J]. 中華放射腫瘤學(xué)雜志, 1993, 2: 60?61

YU Yun, FENG Ningyuan. Determination of spatial dose angular distribution for micro iridium source[J]. Chinese Journal of Radiation Oncology, 1993, 2: 60?61

5 盧峰, 陳睿, 尚云英, 等. 免沖洗膠片法測量后裝機192Ir放射源到位精度及步進距離精度[J]. 中華放射醫(yī)學(xué)與防護雜志, 2014, 34(7): 543?546

LU Feng, CHEN Rui, SHANG Yunying, et al. Flushing-free film test of192Ir accuracy of position and step distance for afterloading systems[J]. Chinese Journal of Radiological Medicine and Protection, 2014, 34(7): 543?546

6 單國平, 徐明清, 嚴英師. 后裝機192Ir放射源到位精度及劑量重復(fù)性測量[J]. 中華放射醫(yī)學(xué)與防護雜志, 2003, 23(3): 203

SHAN Guoping, XU Mingqing, YAN Yingshi. Measurement of position quality assurance and dose repeatability for192Ir source[J]. Chinese Journal of Radiological Medicine and Protection, 2003, 23(3): 203

7 Borca V C, Pasquino M, Russo G, et al. Dosimetric characterization and use of Gafchromic EBT3film for IMRT dose verification[J]. Journal of Applied Clinical Medical Physics, 2013, 14(2): 158?171

8 Rao M. Application of mathematical morphology to biomedical image processing[J]. University of Westminster, 2004, 3(5): 564?567

9 李文華, 張書旭, 徐海榮. 富頸癌后裝治療的劑量分布測試[J]. 中國醫(yī)學(xué)物理學(xué)雜志, 2005, 22(3): 501?504

LI Wenhua, ZHANG Shuxu, XU Hairong. The experiment on the dosage distribution for cervical cancer treated by after-loading radiation[J]. Chinese Journal of Medical Physics, 2005, 22(3): 501?504

10 李浦, 狄小云, 余華, 等. 宮頸癌后裝治療中直腸、膀胱劑量的研究[J]. 中國醫(yī)學(xué)物理學(xué)雜志, 2007, 24(4): 242?244

LI Pu, DI Xiaoyun, YU Hua, et al. Research on the in vivo dose of rectum and bladder in intracavitary brachytherapy of cervix cancer[J]. Chinese Journal of Medical Physics, 2007, 24(4): 242?244

11 FeHner C, Putter R, Knocke T H, et al. Comparison of radiography and computed tomography-based treatment planning in cervix cancer in brachytherapy with specific attention to some quality assurance aspects[J]. Radiotherapy and Oncology, 2001, 58(1): 53?56

12 Das R K, Patel R, Shah H, et al. 3D CT-based high-dose-rate breast brachytherapy implants: treatment planning and quality assurance[J]. International Journal of Radiation Oncology Biology Physics, 2004, 59(4): 1224?1228

CLC TL929

Detection of brachytherapy source dwell positions based on mathematical morphology

JIANG Lipeng1ZHANG Huaiwen2GU Shanshan3
1(Department of Radiation Oncology, First Hospital of Liaoning Medical University, Jinzhou 121001, China)
2(Department of Radiation Oncology, Jiangxi Cancer Hospital, Nanchang 330029, China)
3(Department of Radiation Oncology, Chinese People's Liberation Army General Hospital, Beijing 100853, China)

Background: Accuracy of brachytherapy source positions is an important part of radiotherapy quality assurance using two-dimension ion chamber array or film. It has already been confirmed that the mathematical morphology algorithm can be used for medical image processing. Purpose: This study aims to detect the dwell position of high-dose-rate brachytherapy source using EBT3film based on mathematical morphology. Methods: EBT3film was exposed for different step lengths of source and preset dwell position errors. The scanned digital images were denoised using contourlet transform. Black spots of192Ir source were obtained on the basis of threshold selection and morphological operation. The positions of source were acquired by means of binary centroid algorithm. The calculated source position coordinates were compared with the reference data. Results: The error of source step length was less than ±0.2 mm, while the error of source dwell position was within ±0.5 mm. Conclusion: Mathematical morphology can be used in automatic detection of brachytherapy source positions. This method is simple and reliable.

192Ir, Brachytherapy, Film, Mathematical morphology, Quality assurance

TL929

10.11889/j.0253-3219.2015.hjs.38.060503

江西省衛(wèi)生廳科技計劃項目(No.20141125)資助

姜立朋，男，1977年出生，2013年于承德醫(yī)學(xué)院獲碩士學(xué)位，研究方向為肺癌的放射治療

谷珊珊，E-mail: guss1990@163.com

2015-03-17，

2015-05-05