Ⅰ. INTRODUCTION
비만은 우리나라도 산업의 발달로 인한 생활습관과 신체활동 부족 등으로 비만인구가 급증하고 있다. 비만은 남성의 경우 체지방률이 20% 이상, 여성의 경우 30% 이상일 경우로 정의하고 있으며, 같은 정도의 비만수치일지라도 지방의 분포와 종류에 따라 질병의 유병률도 다르게 나타난다. 특히 복부의 지방 축적은 심혈관계 질환, 뇌졸중, 당뇨병 발생 위험과 밀접한 관련이 있는 것으로 알려져 있다. 우리나라 성인 남성의 경우 전체의 25.2%, 성인 여성의 경우 29.4%가 비만이며, 비만여성 중 44.1%가 복부비만으로 판정되었고[1-3] 특히 45세에서 64세 여성의 경우 61.4%가 복부비만인 것으로 나타나 중년여성의 복부 비만이 심각한 수준인 것으로 나타났다. 복부지방은 피하지방 조직과 내장 지방조직으로 나눠지게 되는데 현재까지 보고되고 있는 많은 연구에서 피하복부 비만 중 내장지방의 과다한 축적이 있는 경우 당대사 뿐 아니라 혈중 중성지방 증가, 고밀도 지단백 콜레스테롤 저하 등의 이상지질혈증도 초래하는 것으로 알려져 있고 심혈관 질환의 독립적인 위험인자이면서 대사증후군 발생에 중요한 역할을 하는 것으로 알려져 있다.[4-5]
이러한 복부 내 지방의 분포를 알기 위한 쉽고 간편하게 측정할 수 있는 연구들은 계속되어왔다. BMI검사, 지질대사 지표 분석, 등이 있고 전산화단층영상을 이용한 방법이 있다. 전산화단층영상을 이용해 지방량을 측정하는 방법에서 반자동방식의 프로그램을 사용 시 체지방과 복부지방의 경계를 의사가 손으로 직접 분할하고 경계를 기준으로 지방의 비율을 자동으로 계산하게 된다. 관련 문제를 해결하기 위한 방법들이[6] 제시되고 있으나 본 연구에서는 모폴로지 연산을 이용한 알고리즘 절차가 간단하고 계산량이 적은 새로운 방법으로 문제를 해결하고자 한다. 본 연구에서는 자동으로 체지방 계산을 시도하고 환자에게 정확한 지방의 정보를 제공을 위한 영상분할을 완성 하고자 한다.
Ⅱ. MATERIALS AND METHODS
1. Image Aquisition
전산화단층촬영장치(TOSHIBA, Tokyo, Japan)로 부터 획득되어진 복부 영상을 실험을 위하여 사용하였다. 전산화단층촬영장치 검사 옵션은 Normal-CT이고 이때 획득되어진 영상의 조건은 슬라이스 두께 7.0 mm, 120 kVp, 150 mA의 조건이다. 후처리과정에서 필터로 FC01를 사용하였다. 프로그램 개발은 PC (Intel® Core 2Duo CPU @ 2.93Ghz, 1.96GB RAM) 와 소프트웨어 (MATLAB R2008b, Massachusetts, U.S.A.)를 기반으로 개발하였다.
2. FAT Calculation ALGORITHM
2.1 연산과정
복부 피하지방량을 자동 계산하기 위하여 Fig. 1의 과정으로 진행하였다. 획득된 전산화단층상으로부터 메디언 필터를 적용하고, 모폴로지 연산을 수행한다. 모폴로지 연산을 위해서 확장과 침식 연산법을 적용하여 최종 경계지역을 획득하였다. 다음 경계로부터 좌표를 획득하여 전체 픽셀의 분포와 피하지방의 픽셀의 분포를 정량적으로 획득하였고 정량적 픽셀 분포로부터 복부 지방율을 계산하였다.
2.2 전처리과정
메디언 필터의 연산방법은 입력 영상에서 필터 윈도우에 포함되는 화소의 크기를 비교한 뒤 오름차순으로 정리하고, 그 중간에서 중간 순위에 해당하는 값을 출력 화소로 결정한다. 영상 한 가운데에 임펄스 잡음이 있다면 제거된다. 입력영상에 대해서 정렬을 한 뒤에 그 중간 값으로 재 할당해서 모든 픽셀에 대해서 연산을 한 후에 잡음제거를 끝내게 된다. 전산화단층영상에서 각각 다른 조작들은 항살 각각의 텍스쳐를 갖고 지방과 근육의 경계, 피질과의 경계는 명확하지 않다. 복부 영상에서 지방질과 근육 피질과의 관계는 조영제를 사용해서 대조도를 향상 시킬 수 없기 때문에 각 조직을 자동으로 구분하는 것은 매우 어렵다. 따라서 콘트라스트를 스트레칭하고 영상의 경계가 일정하지 않거나 단절되어 있을 시 이를 일정하게 연결할 때는 형태학의 처리를 적용한다. 이때 침식처리 연산과 팽창연산을 반복적으로 수행해서 목적에 부합하는 결과를 얻게 된다. 침식연산은 다음과 같이 A와 B는 화소의 집합일때 B로 A를 침식시키면 다음의 식과 같다.[7-8] A와 B를 화소의 집합이라고 했을 때, 형태소 B로 A를 팽창한 식은 다음과 같이 정의할 수 있다.
Fig. 1. Flow chart of fat segmentation program
2.3 경계선 검출
경계선 검출은 영상처리 문제에서 매우 유용하게 쓰일 수 있다. 경계정보는 실험 목적영상의 정량적 분석에 반드시 필요한 단계로서 본 연구에서는 LoG(Laplace of Gaussian)을 사용한다. LoG는 가우지안 스무딩을 수행하여 잡음제거 과정을 거친 뒤 에지를 강조하려고 라플라시안을 이용한 것이다.
Eq. (4)의 LoG 함수 필터의 합이 0이며, 최대값이 커져 경사가 급격한 에지를 검출하는 특징을 가지고 있다. 수행하는 방법으로 가우지안 스무딩을 먼저 수행하고 그 결과 값에 라플라시안을 수행한다.[7-8]
Fig. 2. Preprocessing results
Ⅲ. RESULTS AND REVIEW
1. 복부지방 계산을 위한 영상분할
Fig. 2(a)는 필터링 처리를 통해서 경계의 잡음을 제거한 결과를 나타내었다. 이에 대해서 Fig. 2(b)에서 이진 처리를 하였다. 이진 처리 결과 피질 영역 그리고 복부지방영역 근육 및 내장장기 영역이 시각적으로 구분이 되는 것을 볼 수 있다. 그러나 피질 및 근육주변의 경계들은 연속적이지 않기 때문에 영상분할을 통한 지방부분을 분할 할 때 문제를 야기 시킬 수 있다. Fig. 2(b)의 모폴로지 연산을 통해 침식과 팽창을 반복한 결과 영상으로부터 문제점이 해결된 것을 알 수 있어 결과로부터 각 조직 간의 경계를 더욱 정교하게 얻을 수 있다. Fig. 3에서 LoG(Laplace of Gaussian) 함수를 이용해서 각 조직간 경계 부분을 분할하였다. 각 경계는 명확하게 구분이 되었으며, 피하지방을 계산하기 위한 필요조건이 충분한 것을 보여준다. Fig. 4에서 피하지방의 분할결과를 보여준다. Fig. 5는 내장지방분할에 성공한 결과를 보여준다.
Fig. 3. Subcutaneous fat auto segmentation
Fig. 4. Subcutaneous fat auto segmentation process
Fig. 5. Visceral fat splitting result(blue)
2. 연산시간
복부지방을 자동계측하기 위한 알고리즘 개발의 실험은 PC (Intel® Core 2Duo CPU @ 2.93Ghz, 1.96GB RAM) 와 소프트웨어 (MATLAB R2008b, Massachusetts, U.S.A.)를 기반으로 사용되었다. Fig. 6에서는 제안한 자동복부지방 분할과 기존에 사용되고 있는 수동으로 복부지방 분할을 계산하는 방법 간의 시간에 대한 결과를 보여준다. 10개의 서로 다른 복부영사에 대한 계산 시간을 비교하였다. 제안된 방법은 서로 다른 10개의 영상에 대해서 평균 21.84 s, 표준편차 0.37 s 의 연산시간을 나타내었다. 반면 수동으로 복부지방을 분할하는 방법은 평균 88.02 s 표준편차 6.68 s 이었다. 제안하는 방법과 66.18 s의 차이를 보였고 제안하는 방법이 임상에 적용하기에 매우 유리함을 알 수 있었다.
Fig. 6. Comparison of body fat calculation time (x-axis: sample, y-axis: sec)
Ⅳ. CONCLUSION
개발된 알고리즘을 통해서 표피와 피하지방 그리고 내장지방의 경계를 자동분할 하였다. 피하지방 경계 좌표를 자동으로 얻을 수 있다. 개발된 분할결과를 사용하여 자동 지방량 계산을 할 수 있다. 정확한 분할 도구를 제공함으로써 의사에게 편리함을 주고 재검사로 인한 피폭과 검사비용을 줄이는데 도움이 될 것으로 판단한다.
Acknowledgement
본 연구는 2018년 부산가톨릭대학교 교내연구비의 지원에 의해서 수행되었다.
참고문헌
- B. J. Kwon, Y. S. Kim, H. C. Rhim, "Three-phase Dynamic CT Findings of Liver Abscess: Related Factors with Multiple Layering Enhancement Pattern," Journal of the Korean Radiological Society, Vol. 44, No 1, pp. 69-75, 2001. https://doi.org/10.3348/jkrs.2001.44.1.69
- H. M. Yoo, J. S. Cho, K. S. Shin, "Usefulness of Helical CT for the Preoperative Evaluation of Small Advanced Gastric Cancer Mimicking as Early Gastric Cancer at Endoscopy," Journal of the Korean Radiological Society, Vol. 52, No. 6, pp. 385-393, 2005. https://doi.org/10.3348/jkrs.2005.52.6.385
- D. J. Kim, J. S. Yu, S. M. Lee, "Advanced Gastric Cancer: Differentiation of Borrmann Type versus Borrmann Type by Two-Phased Dynamic Multi-Detector Row CT with Use of the Water Filling Method," Journal of the Korean Radiological Society, Vol. 68, No. 2, pp. 117-124, 2013. https://doi.org/10.3348/jksr.2013.68.2.117
- B. C. Yoo, D. C. Kwon, J. S. Lee, K. J. Jeon, S. H. Kim, Y. S. Kim, D. C. Kweon, "Comparison radiation dose of z-axis automatic tube current modulation technique with fixed tube current multi-detector row CT scanning of lower extremity venography," Journal of Radiation. Protection and Research," Vol. 32, No. 3, pp. 123-133, 2007.
- H. J. Kim, J. H. Cho, C. S. Park, "Evaluation of Image Quality in Low Tube-Voltage Chest CT Scan" Journal of Radiation protection, Vol. 35, No. 4, pp. 135-141, 2010.
- Okumura M, Ota T, Tsukagoshi S, Katada K. New Method of Evaluating Edge-preserving Adaptive Filters for Computed Tomography(CT) : Digital Phantom Method. Japanese Journal of Radiological Technology, Vol. 7 No. 1, pp. 971-978, 2006.
- H. S. Yoon, Predict of Fracture Risk Rate According to Morphological Measuring of Proximal Femoral Part Using Dual Energy X-ray Absoptiometry. Journal of Korean Society of radiological technology, Vol. 25, No. 4, pp. 49-53, 2002.
- S. J. Kim, S. W. Yoo, J. C. Kim, S. M. Noh, J. W. Park, A Study on the Segmentation of Kidney from the Non Enhanced CT Image : Journal of KIISE Vol. 27, No. 3, pp. 514-516, 2000.
- M. D. Becker, Y. M. Berkmen, J. H. M. Austin, I. K. Mun, B. M. Romney, A. Rozenshtein, Lung Volumes before and after Lung Volume Reduction Surgery. Vol. 157, No. 6, pp. 1593-1599, 1998. https://doi.org/10.1164/ajrccm.157.5.9706066
- S. S. Hong, J. Huh, " Uses Edge Method, MTF Analyses of CT Equipment which Follows in Channel and Convolution Algorithm, Journal of Korean Society of Computed Tomography Technology, Vol. 13, No. 1, pp. 105-114, 2011.
- H. S. Hwang, K. B. Lee, S. M. Choi, "Evaluation of Image Quality with the use of SAFIRE(Sinogram Affirmed Iterative Reconstruction) and ASIR(Adaptive Statistical Iterative Reconstruction) by changing Radiation Dose," Journal of Korean Society of Computed Tomographic Technology, Vol. 15, No. 1, pp. 215-221, 2013.
- I. S. Lee, M. C. Kim, Y. C. Nam D. K. Han, "Usefulness of Automatic Exposure Control for dose Reduction: in Neck Study Using Multi_detector Computed Tomography," Journal of Korean Society of Computed Tomographic Technology, Vol. 2 No. 1, pp. 7-12, 2010.
- K. J. Park, M. C. Kim, D. K. Han, "A Study on Image Quality and Dose of CT Image Using Adaptive Statistical Iterative Reconstruction Algorithm," Journal of Korean Society of Computed Tomographic Technology, Vol. 12, No. 1, pp. 95-104, 2010.
- M. J. Uhm, S. M. Cheong, H. J. Yoo. "The Evaluation of Comparison between Dose and image Quality using spiral mode and Flash spiral MODE According to Pitch of Head Angiography," Journal of Korean Society of Computed Tomographic Technology, Vol. 14 No. 1 pp. 85-93, 2012.
- D. Y. Lee, Y. S. Jang, Y. B. Kim S. M. Lim, K. W. Lee, Y. G. Jo, "The Variation of the Dose and Noise of HRCT according to Tube Voltage in MDCT," Journal of Korean Society of Computed Tomographic Technology, Vol. 14, No. 1, pp. 103-110, 2012.