Ⅰ. INTRODUCTION
균일한 물팬텀을 컴퓨터단층촬영(Computed Tomo graphy, CT) 스캐너로 스캔하여 영상화하고 국소적 영역에 대한 CT수(CT number)를 보면, CT수가 모두 동일하지는 않고 평균값 주위에서 변한다는 것을 알 수 있다. 이 변화가 시스템의 노이즈이다. 그리고 노이즈는 정상조직과 병리조직 사이에 자연적으로 발생하는 감쇠계수의 차이를 작게 만들기 때문에 CT 스캐너 성능의 매우 중요한 척도로 사용된다[1]. CT 스캐너의 노이즈 측정값은 일반적으로 균일한 물팬텀을 스캔하여 얻은 영상에서 16 cm2의 관심영역(region of interest, ROI)을 설정하여 얻어지는 표준편차로 표시한다[1-3]. 또한 대조도 척도를 결정할 수 있다면 노이즈 측정값은 물 유효선감쇠계수의 백분율로 표시할 수 있다[1,4]. 표준편차에 의한 노이즈 측정은 아주 쉽지만 물 유효선감쇠계수의 백분율에 의한 노이즈 측정은 아주 어렵고, 발표된 자료가 거의 없는 실정이다. 따라서 물 유효선감쇠 계수의 백분율에 의한 노이즈를 측정할 수 있는 방법에 관한 연구의 필요성이 대두된다.
이 연구의 목적은 임상에서 AAPM CT 성능 팬텀을 이용한 CT영상 노이즈 정도관리에 사용할 수 있는 물 유효선감쇠계수의 백분율에 의한 노이즈 측정의 방법을 제시하는데 있다.
본 연구는 120 kVp CT X-선 빔으로 AAPM CT 성능 팬텀을 스캔하여 얻어진 CT영상에서, CT수직선성 삽입부분의 각 핀과 물에 대하여 평균 CT수를 측정하였다. 유효에너지(유효광자에너지)는 각 핀과 물에 대하여 측정된 평균 CT수와 광자에 너지별 선감쇠계수의 선형회귀분석의 상관계수로 부터 결정하였다. 또한 물 유효선감쇠계수는 결정된 유효에너지에 대응하는 물 선감쇠계수로 결정하였다. 그리고 물과 아크릴에 대하여 측정된 평균 CT수와 유효선감쇠계수로부터 대조도 척도를 구하였다. 구하여진 대조도 척도, 물 유효선감쇠계수, 그리고 AAPM CT 성능 팬텀에서 고정핀 부분의 물에서 측정된 표준편차를 사용하여 물 유효선감쇠계수의 백분율에 의한 노이즈를 측정하는 방법을 제시하였다.
Ⅱ. MATERIAL AND METHODS
실험에는 AAPM CT 성능 팬텀(610, Computerized Imaging Reference Systems Cor., USA, 2022)[5]을 CT 슬라이스 영상으로 스캔하기 위하여 CT 스캐너(LightSpeed16, General Electiric Cor., USA, 2019)를 사용하였다. 스캔된 CT 슬라이스 영상에서 평균 CT수를 측정하기 위하여 의료영상저장전송시스템(INFINITT STARPACS, PiviewSTATTM, Infinity Cor., Korea)을 사용하였다.
1. CT 슬라이스 영상 스캔
AAPM CT 성능 팬텀(AAPM CT Performance Phantom)에서 CT수 직선성 삽입부분(CT number linearity insert Part)과 고정핀 부분(Alignment pin Part)을 스캔하기 위하여, 공통 스캔 조건으로 관전압은 120 kVp, scan FOV는 25 cm, display FOV는 25 cm, 그리고 빔 콜리메이션은 10 mm로 설정하였다.
1.1. CT수 직선성 삽입부분 스캔
공통 스캔 조건에 관전류량 250 mAs를 추가로 설정하여, 5번 single slice conventional scan을 시행하였다.
1.2. 고정핀 부분 스캔
공통 스캔 조건에 관전류량 100, 150, 200, 250, 300 mAs를 추가로 설정하여, 각각 5번 single slice conventional scan을 시행하였다.
2. 평균 CT수 측정
스캔 시행 후 영상의 재구성은 standard reconstruction algorithm을 사용하여 CT 슬라이스 영상으로 구성하였다.
2.1. CT수 직선성 삽입부분의 각 핀과 물의 평균 CT수 측정
Fig. 1에서 보여주는 바와 같이, CT 슬라이스 영상에서 각 핀의 물질은 polyethylene, polystyrene, nylon, acrylic, polycarbonate이고, 바탕물질은 물이다. 각 핀과 물에 관심영역(ROI, 5.6 cm2)을 설정하여 CT수(ROI의 평균 CT수)를 측정하였다. 5번 측정값의 평균값으로 평균 CT수는 산출되었다.
Fig. 1. Measurement of mean CT number for each peg and water.
2.2. 고정핀 부분의 물에서 표준편차 측정
Fig. 2에서 보여주는 바와 같이, CT 슬라이스 영상에서 6시 방향에 관심영역(ROI, 16.02 cm2)을 설정하여 CT수(ROI의 평균 CT수)를 측정하였다. 5번 측정값의 평균값으로 평균 CT수와 표준편차는 산출되었다.
Fig. 2. Measurement of standard deviation at water.
3. 질량감쇠계수로부터 선감쇠계수 유도
AAPM CT 성능 팬텀의 CT수 직선성 삽입부부분의 각 핀과 물 물질은 Table 1에서 보여준다.
Table 1. Chemical formula and density of each peg and water
NIST 웹사이트에 들어가서 XCOM: Photon Cross Sections Database에서 Contents의 6. Database Search Form[6] 양식에 Table 1[7]의 화학식을 입력하여, 광자에너지(단일광자에너지) 60 keV와 80 keV에 대한 각 핀과 물의 질량감쇠계수를 얻었다. 이 자료는 Origin pro 8.0 통계프로그램의 linear 내삽법으로 광자에너지 0.5 keV 간격의 질량감쇠계수를 산출하였다. 선감쇠계수는 산출된 질량감쇠계수에 밀도를 곱하여 구하였다.
4. 유효선감쇠계수 유도
CT 스캐너에서 발생된 X-선 빔의 유효에너지는 CT수 직선성 측정으로부터 얻을 수 있다[4]. AAPM CT 성능 팬텀의 CT수 직선성 삽입부분에서 측정된 각 핀과 물의 평균 CT수와 광자에너지별 각 핀과 물의 선감쇠계수를 선형정합하여 얻은 상관계수들 중 가장 큰 상관계수를 갖는 광자에너지는 유효에너지로 결정된다[7-11]. 이 결정된 유효에너지에 대응하는 각 핀과 물의 선감쇠계수가 유효선감쇠계수에 해당된다[4].
5. 물 유효선감쇠계수의 백분율
물 유효선감쇠계수(μw)의 백분율(%σμw)에 의한 노이즈(N) 측정은 Eq. (1)로 기술된다[4].
\(\begin{aligned}\% \rho \mu_{w}=\frac{\sigma \cdot C S \cdot 100}{\mu_{w}}\end{aligned}\) (1)
여기서 σ는 ROI에서 측정된 픽셀값들의 표준편차(standard deviation), 그리고 CS는 대조도 척도(contrast scale)이다. CS는 두 개의 알려진 물질(material)에서 주어진 CT수의 변화에 대한 유효선 감쇠계수의 변화로 정의된다. 일반적으로 하나의 물질은 물이고, 다른 하나는 CT 팬텀의 일반적인 구성 물질인 아크릴이다. 여기서 다른 하나는 아크릴이 아닌 낮은 원자번호(Z)를 갖는 다른 기준 물질을 사용할 수 있다. 대조도 척도[12-15]는 Eq. (2)로 기술된다[4].
\(\begin{aligned}C S=\frac{\mu_{m}(E)-\mu_{w}(E)}{C T_{m}-C T_{w}}\end{aligned}\) (2)
여기서 μm(E)과 μw(E)는 기준 물질과 물의 유효선감쇠계수이고, CTm과 CTw는 기준 물질과 물에서 측정된 CT수이다. 이 연구에는 물과 아크릴을 사용하여 Eq. (2)에 의해 대조도 척도를 산출하였다. 그리고 물 유효선감쇠계수(μw)의 백분율은 Eq. (2)에 의해 산출된 대조도 척도와 고정핀 부분의 물에서 측정된 표준편차를 Eq. (1)에 대입하여 구하였다.
Ⅲ. RESULT
AAPM CT 성능 팬텀에서 CT수 직선성 삽입부분의 각 핀과 물에서 측정된 평균 CT수는 Table 2, 그리고 고정핀 부분의 물에서 측정된 표준편차는 Table 3에서 보여준다. NIST 웹사이트의 XCOM: Photon Cross Sections Database로부터 광자에너지 60 ~ 80 keV 범위에서 각 핀과 물에 대한 질량감쇠 계수는 Table 4에서 보여준다. Table 4의 질량감쇠 계수에 밀도를 곱하여 산출한 선감쇠계수는 Table 5에 나타내었다.
Table 2. CT numbers measured for each peg and water in the CT number linearity insert part
Table 3. Standard deviation measured in water at the alignment pin part
Table 4. Mass attenuation coefficients(g/cm2)
Table 5. Linear attenuation coefficients(cm-1)
CT스캐너에서 120 kVp CT X-선 빔의 유효에너지는 AAPM CT 성능 팬텀의 CT수 직선성 삽입부분에서 측정된 각 핀과 물의 평균 CT수와 광자에너지별 각 핀과 물의 선감쇠계수를 선형정합하여 얻은 상관계수들 중 가장 큰 상관계수를 갖는 광자에너지를 유효에너지로 결정하였다. 이 결정된 유효에너지는 Table 6에서 광자에너지 68 keV이다. 이 유효에너지에 대응하는 선감쇠계수가 유효선감쇠계수로 된다. 이 유효선감쇠계수는 Table 7에 나타내었다.
Table 6. Determination of effective energy from the largest correlation coefficient
Table 7. Effective linear attenuation coefficients(cm-1) at effective energy
대조도 척도는 AAPM CT 성능 팬텀에서 CT수직선성 삽입부분의 물과 아크릴에 대한 측정된 평균 CT수와 유효선감쇠계수를 Eq. (2)에 대입하여 0.000188 cm-1 · HU-1을 얻었다. 물 유효선감쇠계수의 백분율은 Eq. (2)에 의해 얻어진 대조도 척도와 고정핀 부분의 물에서 측정된 표준편차, 물의 유효선감쇠계수 0.197 cm-1를 Eq. (1)에 대입하여 산출한 결과는 Table 8에서 보여준다. 이 결과로부터 물 유효선감쇠계수의 백분율에 의한 노이즈 측정값은 관전류량 100 ~ 300 mAs의 범위에서 0.31 ~ 0.52%를 얻었다.
Table 8. Percentage of the effective linear attenuation coefficient of water(%σμw).
Ⅳ. DISCUSSION
CT스캐너에서 120 kVp X-선 빔의 유효에너지는 CT수 직선성 측정으로부터 얻을 수 있다[4]. 이것은 광자에너지별 AAPM CT 성능 팬텀에서 CT수 직선성 삽입부분의 각 핀과 물의 선감쇠계수에 대한 120 kVp X-선 빔의 스캔으로 얻어진 CT수 직선성 삽입 부분의 CT영상에서 각 핀과 물에서 측정된 평균 CT수를 선형정합 함으로써 얻어지는 상관계수들 중 가장 큰 상관계수를 갖는 광자에너지를 유효에너지로 결정한다는 의미이다. 이 유효에너지에 대응하는 선감쇠계수는 유효선감쇠계수로 표기된다. 실험으로 부터 120 kVp X-선 빔의 유효에너지는 68 keV로 결정되었다. 이 결정된 유효에너지는 AAPM report No. 39에서 보고한 120 ~ 125 kVp CT X-선 빔의 유효에너지 70 keV 정도와 근사적으로 일치하였다[4]. 그리고 유효선감쇠계수는 유효에너지에 대응하는 각 핀과 물의 선감쇠계수로 결정하였다.
대조도 척도는 AAPM CT 성능 팬텀에서 CT수직선성 삽입부분의 물과 아크릴에 대한 측정된 평균 CT수와 유효선감쇠계수를 Eq. (2)에 대입하여 구하였다. 대조도 척도 산출 시 일반적으로 하나의 물질은 물을 선택하고, 다른 하나의 물질은 아크릴을 사용한다.
관전압은 120 kVp로 고정하고, 관전류량은 100, 150, 200, 250, 300 mAs로 다르게 설정하여 얻어진 AAPM CT 성능 팬텀의 고정핀 부분의 CT영상에서, 물 유효선감쇠계수의 백분율에 의한 노이즈 측정값은 0.31 ~ 0.52%로 얻었다. 이 측정된 노이즈값들은 관전류량이 증가할수록 감소하는 양상을 나타내었다. 또한 표준편차로 측정한 노이즈 측정값들과 같은 양상을 나타내었다. 이 이유는 관전류량을 작게 설정하면 X-선관에서 발생되는 광자수가 적게 발생되어 불균일하게 분포하며, 적게 발생된 광자수가 검출기에 도달하여 적게 검출됨으로써 광자수 부족으로 인한 양자 노이즈 신호를 발생하기 때문이다.
임상에서 특수의료장비의 설치 및 운영에 관한 규칙에 근거하여, 노이즈 정도관리검사는 6개월에 한 번 한국의료영상품질관리원에 신청하여 검사를 받고 있다. 이때 노이즈 측정은 표준편차로 한다. 허용오차는 5 HU 이내이다. 노이즈 측정값이 허용 오차 범위를 벗어날 경우, 그 원인은 검출기 문제, 선원 변동, 조준기 오정렬, X-선관 발생기 오교정 등의 문제로 추정할 수 있다.
Ⅴ. CONCLUSION
120 kVp CT X-선 빔으로 AAPM CT 성능 팬텀을 스캔하여 얻어진 CT영상에서, CT수 직선성 삽입부분의 각 핀과 물에 대하여 평균 CT수를 측정하였다. 유효에너지는 각 핀과 물에 대하여 측정된 평균 CT 수와 광자에너지별 선감쇠계수의 선형회귀분석의 상관계수들로부터 가장 큰 상관계수를 갖는 광자에너지로 결정하였다. 그리고 물과 아크릴에 대하여 측정된 평균 CT수와 유효선감쇠계수로부터 대조도척도는 산출되었다. 산출된 대조도 척도, 물의 유효선감쇠계수, 그리고 AAPM CT 성능 팬텀에서 고정핀 부분의 물에서 측정된 표준편차를 사용하여, 물유효선감쇠계수의 백분율에 의한 노이즈를 측정하였다. 결론적으로 이 연구는 물 유효선감쇠계수의 백분율에 의한 노이즈를 측정한 절차와 방법을 제시한다.
참고문헌
- AAPM report No. 1, "Phantoms for Performance Evaluation and Quality Assurance of CT Scanners", American Association of Physicists in Medicine, pp. 1-23, 1977.
- M. E. Phelps, E. J. Hoffman, M. M. Ter-Pogossian, "Attenuation Coefficients of Various Body Tissues, Fluids and Lessions at Photon Energies of 18 to 136 keV", Radiology, Vol. 117, No. 3, pp. 573-583, 1975. https://doi.org/10.1148/117.3.573
- P. S. Rao, E. C. Gregg, "Attenuation of Monoenergetic Gamma Rays in Tissues", American Journal of Roentgenology, Vol. 123, No. 3, pp. 631-637, 1975. https://doi.org/10.2214/ajr.123.3.631
- AAPM report No. 39, "Specification and Acceptance Testing of Computed Tomography Scanners", American Association of Physicists in Medicine, pp. 1-95, 1993.
- http://www.doza.ru/docs/med/aapm_ct.pdf.
- https://www.nist.gov/pml/xcom-photon-cross-sections-database.
- M. R. Millner, W. H. Payne, R. G. Waggener, W. D. McDavid, M. J. Dennis, V. J. Sank, "Determination of Effective Energies in CT Calibration", Medical Physics, Vol. 5, No. 6, pp. 543-545, 1978. http://dx.doi.org/10.1118/1.594488
- J. E. Kim, "Determination of Effective Energy of CT X-ray Beams", Journal of the Korean Society of Radiology, Vol. 13, No. 4, pp. 517-522, 2019. https://doi.org/10.7742/jksr.2019.13.4.517
- J. E. Kim, S. H. Lee, "Derivation of the Effective Energy Calculation Formula of the X-ray Beam Generated by the CT Simulator", Journal of the Korean Society of Radiology, Vol. 15, No. 6, pp. 869-875, 2021. https://doi.org/10.7742/jksr.2021.15.6.869
- P. F. Judy, "Comparison of Equivalent photon Energy Calibration Methods on Computed Tomography", Medical Physics, Vol. 7, No. 6, pp. 685-691, 1980. http://dx.doi.org/10.1118/1.594722
- S. Okayama, T. Soeda, Y. Takami, R. Kawakami, S. Somekawa, S. Uemura, Y. Saito, "The Influence of Effective Energy on Computed Tomography Number Depends on Tissue Characteristics in Monoenergetic Cardiac Imaging", Radiology Research and Practice, Vol. 2012, pp. 1-7, 2012. https://doi.org/10.1155/2012/150980
- J. Nosil, K. I. Pearce, R. A. Stein, "Linearity and Contrast Scale Control in Computed Tomography", Medical Physics, Vol. 16, No. 1, pp. 110-113, 1989. https://doi.org/10.1118/1.596396
- J. E. Barnes, "Characteristics and Control of Contrast in CT", Radiographics, Vol. 12, No. 4, pp. 825-837, 1992. https://doi.org/10.1148/radiographics.12.4.1636042
- C. D. Njiki, J. N. Manyol, Y. E. Yigbedeck, D. W. Abouou, B. C. Yimele, J. F. Sabouang, "Assessment of Image Quality Parameters for Computed Tomography in the City of Yaounde", Open Journal of Radiology, Vol. 8, No. 1, pp. 37-44, 2018. https://doi.org/10.4236/ojrad.2018.81005
- M. J. Medinaa, F. J. Funesa, A. Palaub, V. M. Sanz, A. R. Osorioa, M. Namiasb, "Acceptance testing of a CT scanner with a knowledge-based iterative reconstruction algorithm", Brazilian Journal of Radiation Sciences, Vol. 7, No. 1, pp. 1-17, 2019. https://doi.org/10.15392/bjrs.v7i1A.456