I.서 론
X선속 금속 흡수체를 인위적으로 삽입하는 것을 부가여 과(added filtration)라 한다1). 이러한 filter의 두께를 증가 시켜 가면 spectrum의 형태는 고 에너지 쪽으로 이동하지 만, 전체적인 선량률의 감소를 피할 수 없다. Filtration는 환자의 피폭량을 높이며, 영상의 대조도를 낮추는 방사선 광자에 대하여 영상 형성에 유용한 광자의 비율을 높이는 과정이다. 이것은 주로 광전효과가 많이 발생할 때 감약이 심하고, 콤프톤 반응이 많이 일어날 때에는 감약이 감소하 는 현상을 이용하는 것이다2). Filter는 환자선량을 어느 정 도 감소할 수 있으나, 일정 에너지 영역의 광자만 흡수하는 것이 아니라 모든 광자의 에너지를 흡수하기 때문에, 이 손 실을 보상하기 위하여 선량을 증가시켜야 한다. 중금속 필 터나 K-edge filter를 사용하면, 저에너지와 고에너지를 감 소하여 비교적 에너지 spectrum을 좁게 할 수 있다는 장점 이 있다1). X선속에서 저 에너지를 제거하면 환자선량을 감 소할 수 있고, 고 에너지를 제거하면 영상대조도를 향상할 수 있는 장점이 있다3).
Filter 연구의 중요한 조건은 에너지 스펙트럼의 변화와 물질의 특성에 관한 것이다. 진단용 X선 스펙트럼 예측은 1923년 Kramers에 의해 최초로 시도되었으며, 오늘날에 이르러 다양한 방법들이 연구되고 있다4). 에너지 스펙트럼 측 정은 경험적 모형, 반경험적 모형, 몬테칼로 모형이 있다. 경 험적 모형은 특별한 장치가 요구되며, 선속밀도(flux density) 가 높으면 정확한 측정이 어려운 것으로 알려져 있다5). 검 출기 반응도(response function)가 에너지에 의존하고, 대 부분의 검출기에 대한 반응도가 정확하게 알려져 있지 않는 문제점도 있다6). 몬테칼로 모형은 다양한 변수와 한계를 극 복하여 적용할 수 있으며, 경험적 모형이나 반경험적 모형 에 비해 계산시간이 증가하는 문제를 가지고 있지만, 정확 한 물리 모델과 상호작용과정을 적용하기 때문에 정확한 스 펙트럼을 계산할 수 있다.
복합 filter는 두 개 또는 여러 개의 다른 금속 층으로 이 루어졌으며, Al(알루미늄)과 Cu(구리)를 이용한 복합 filter 를 이용할 경우, 항상 X선 관측으로부터 원자번호가 높은 금속을 놓는데, 그 이유는 방사선의 여과작용은 Cu가 담당 하고, Al은 Cu에서 발생한 특성방사선을 흡수할 목적으로 사용하기 때문이다. Cu의 특성방사선 에너지는 환자까지 도달할 수 있는 충분한 에너지(8keV)인 까닭에, 피부선량을 증가시킬 염려가 있으므로 Al을 삽입하여, 이 특성방사선을 흡수해야 한다. Al에서도 특성방사선이 발생하지만, 이 에 너지는 작아서(1.5 keV) 환자에 도달하기 전에 공기 내에서 흡수 된다7).
Rymond et al.은 농도 1.0로 동일하게 나타내는데 표면선량 과 환자의 피폭의 원인이 되는 흡수선량이 Cu< Al< non-filter 로 나타났다고 했다8). 또한 125 kV의 tube voltage와 0.3 mm Cu filter를 사용했던 팬텀 실험에서 두께를 달리한 팬텀에 서, 각각의 표면입사선량 측정 결과 filter를 사용했을 경우 small 팬텀에서 44, mid 46, large 47% 감소율의 차이를 나타 냈다9). Lee et al.은 머리부 전후방향촬영 시 Cu와 Al filter의 16가지 종류의 조합을 통하여 입사표면선량을 측정 결과, 같은 두께에 Cu에서 복합 filter를 사용한 경우가 모두 선량이 작게 측정되었다10). 또한 유방촬영장치의 부가필터에 따른 선량변 화 실험에서, 단일 필터보다 복합filter를 사용 시 선량 감소효 과가 있다고 주장했다1). 선행 연구들에서 Cu filter에서 복합 filter의 역할이 Cu의 특성 X선 제거를 통한 표면선량 감소에 효과가 증거 되었고, 이에 몬테카를로 시뮬레이션에 의한 에너 지 스펙트럼 측정하여 분석하고, 화질 평가를 통하여 복합 filter 사용 시 Al filter의 사용 유, 무에 따른 효율성을 분석하 여 환자에게 유리한 피폭선량 관리를 통하여 환자 안전을 도모 하고자한다 .
II.연구재료 및 방법
본 연구에서는 장비는 Philips digital diagnost value amorphous silicon flat panel detector Pixium 4600를 이용하여 여러 평가 항목들을 측정한다. Trixel사의 Pixium 4600과 DR의 특성은 다음과 같다(Table 1).
1.연구 재료
1)Program
화질 평가를 위해 DICOM image processing software (Image J v.1.43u; National Institutes of Health Bethesda, MD, USA), 통계 사용을 위하여 SPSS(Statistical Package for Social Science v.18), 에너지 스펙트럼 측정을 위하여 SRS-78, GATE(Geant4 Application for Tomography Emission v.6.2) Simulation을 사용하였다.
2)Anthropomorphic Phantom
인체모형 팬텀은 인체부위별 X선 팬텀으로 남성으로 175cm, 74kg 이다(Figure 1, 2, 3).
3)Added Filter
부가필터는 99.9%순도에 크기는 가로, 세로 각각 10 cm 이며 두께는 Cu는 0.1, 0.2, 0.3 mm, Al는 1 mm이다 (Figure 4, 5).
4)Dosimeter
선량계는 Unfors사의 Xi 모델을 이용하였다(Figure 6).
2.연구 방법
1)에너지 스펙트럼 측정
Simulation energy spectrum을 얻기 위해 Genat4(geometry and tracking) Toolkit를 기반으로 한 GATE를 이용하였다. GATE를 이용한 energy spectrum을 얻기 위해선 X-ray source 와 detector에 대한 geometry code 작성해야 한다. Detector는 실제와 같이 작성하며 가장 중요한 것은 scintillator의 두께가 중요하다. 본 연구에선 600 μm를 입력하였다. GATE 내에서는 X-ray source의 spectrum에 대한 정보가 없으므로, SRS-78 이라는 프로그램을 이용하여 spectrum 정보를 획득하여 X-ray source code로 구현하였다. 에너지는 50, 70, 90, 120 kV로 필터물질은 Cu 0.1, 0.2, 0.3 mm, Cu 0.1, 0.2, 0.3 mm에 각각 Al 1 mm를 합친 물질로 실험을 하였다.
2)에너지 스펙트럼의 평균 비교
각 에너지별 같은 두께의 Cu에 Al filter 사용 전, 후 spectrum의 평균 비교를 실시하여 통계적 유의성을 살펴보 았다. 먼저 정규성 분포 test 후 정규성 분포를 만족하지 않 아서 비모수적 통계방법을 하였다. Al filter 사용 전, 후를 대응 Wilcoxon sign ranks test를 실시하였고 통계량(Z), p-value를 구하였다.
3)에너지 스펙트럼의 저에너지 비교 분석
각 에너지에 Al filter 사용 전, 후의 스펙트럼 그래프가 가우스 분포를 이루기 위해 상승하기 전 저에너지 부분의 광자수를 총 광자수로 나누어 백분율로 환산하여 정량화 하 였다.
4)에너지 스펙트럼의 중심 에너지 비교 분석
같은 에너지에 Al filter 사용 전, 후의 광자양의 차이를 사용 전, 후의 광자양 합산 평균으로 나누어 백분율로 환산 하여 정량화 하였다.
5)에너지 스펙트럼의 고에너지 비교분석
양질의 고에너지의 위치차를 알아보기 위하여 기술통계량 의 왜도 값을 이용하였다. Al filter 사용 전, 후의 평균값이 다르고, 그냥 왜도 값을 구했을 시에는 정규 분포 값의 차이 밖에 알 수 없어서 각 에너지에서 Al filter 사용 전, 후 값을 표준화 하여 같은 조건으로 만들어서 왜도를 구하여. 어떤 실험 조건이 더 고 에너지의 위치하는 지를 살펴보았다.
6)SNR, PSNR, RMSE, MAE
Chest 팬텀을 이용하여 필터 장착 전 기준 영상을 촬영하 고 각각의 필터를 장착 한 실험영상을 촬영한다. 촬영 조건 은 AEC(Auto Exposure Control, 이하 AEC)를 사용하여 얻어진 입사 X-ray 평균에너지는 125 kV에 4 mAs를 고정 하고 SID는 180 cm, 12:1격자를 사용하였다. 이렇게 얻어 진 영상을 Image J를 이용하여 비교 분석하였다. PSNR은 서로 다른 두 영상에 대한 차이를 정량적으로 느끼고 알 수 있도록 숫자로 나타낸 값이며, 이는 다음과 같이 식1을 사용 하여 산출한다.
MAE는 식2를 사용하여 산출한다.
I(i,j) : 좌표(i,j)의 비교영상 픽셀값 K(i,j) : 좌표(i,j)의 원래영상 픽셀값
7)CNR, SNR 측정
CNR는 Pelvis 팬텀을 이용하여 촬영 조건은 AEC를 사용 하여 얻어진 입사 X-ray 평균에너지는 77 kV에 12.5 mAs를 고정하고 SID는 110 cm, 10:1격자를 사용하였다. Right iliac wing, femur head, sacrum을 Image J 프로그램을 이용하여 0.3 X 0.3 cm2의 ROI를 설정하였다. ROI의 pixel 신호 강도 의 표준편차를 측정한 후, ROI를 제외한 Background의 신호 강도 평균값과 표준편차를 구하여 식3에 의하여 계산한다.
B: Background
R; ROI
또한 Skull 팬텀을 이용하여 조건은 AEC를 사용하여 얻 어진 77 kV에 12.5 mAs를 고정하고 SID는 110 cm, 10:1격 자를 사용하였다. Frontal bone, Right orbit, mandible을 위의 조건과 동일하게 영상을 획득하였다. SNR은 CNR 측정 과 동일하게 Image J 프로그램을 이용하여 ROI의 신호강도 의 평균과 표준편차를 측정한 후, ROI를 제외한 Background 의 신호강도 평균과 표준편차를 구하여 식4를 이용하여 계 산하였다(Figure 7).
B: Background
R; ROI
8)정성적 평가
Chest phantom에 120 kV에서 AEC mode 사용하여 EI 200값을 설정 후 조사하여 얻은 4 mAs값을 이용하였다. 10 년 이상의 5명의 방사선사가 흉부 X선 검사 임상영상평가표 (국가암검진 임상평가표)를 기준으로 하여 5점 리커트 척도 를 이용하였다. 평가항목은 평가표에서 투과 상태, 해상도 및 대조도 평가 문항을 이용하였다. 평가항목은 총 7가지 항 목으로 구성하였다(Table 2).
III.결 과
1.에너지 스펙트럼
1)에너지 스펙트럼 측정
GATE simulation을 이용하여 각 filter의 에너지별 광자 량을 지수값으로 측정하였다. 측정 후 아래와 같은 에너지 스펙트럼을 작성하였다. 입사 관전압이 50, 70, 90, 120 kV 로 높아질수록 에너지 스펙트럼이 높은 keV 영역으로 이동 하여 분포가 넓어지고 최대 분포값도 넓게 펴져 평활화 됨 을 알 수 있다(Figure. 8). 에너지를 고정으로 놓고 filter 두께를 변화시키면, 두꺼워 질수록 낮은 keV의 스펙트럼 분 포는 적어지며 높은 keV로 옮겨가지만 전체적의 광자의 수 도 줄어든다(Figure 9).
2)에너지 스펙트럼의 평균 비교
스펙트럼 곡선의 정점도는 관전압이 증가하거나 필터의 물질이 두꺼워지면 높은 쪽으로 옮겨갔다. 같은 두께에서 Al filter 사용 전, 후 평균 keV의 이동을 살펴보면 0.1 mm 두께에서는 50 kV에서 46에서 47 keV로 약간의 변화가 있 었으나, 120 kV에서는 51, 69 keV로 동일했다. 0.3 mm 에 서는 70 kV를 제외하고는 평균 keV가 변화가 없었다. 120 kV에서는 두개의 곡선의 정점을 보였는데 두 번째 정점은 모두 69 keV로 동일하였다. 필터의 두께가 증가되면 에너 지도 오른쪽으로 옮겨가게 되는데 90 kV에서 단일 filter 시 49~56 keV이고, 복합 filter 사용 시는 49~56 keV로 차이가 없었다(Table 3).
Wil coxon sign r anks t est를 실시하였고 통계량(Z), p-value를 구하였다. Al filter 사용 전, 후 전체적인 광자 량 값의 p<0.05 이므로 통계적으로 유의한 차이가 있었다 (Table 4).
3)에너지 스펙트럼의 저 에너지 비교 분석
50 kV의 0.1 mm에서 0.39%, 120 kV의 0.3 mm에서 3.3%이였으며, Al filter를 사용한 복합 filter에서는 50 kV 의 0.1 mm에서 0.3%, 120 kV의 0.3 mm에서 3.4%이였다. 대부분 많은 차이를 보이지 않았으며, 90 kV의 0.1 mm에 서 0.2 %차이가 가장 컸다(Table 5).
4)에너지 스펙트럼의 중심 에너지 비교 분석
전반적인 에너지 스펙트럼의 중심 에너지 비교는 Cu filter를 단일로 사용한 경우가 우수하였다. 두 물질의 동일 한 에너지 대역대의 총 광자수의 평균치에서 차이를 나누어 peak치 평균에서 얼마나 차지하는가 하는 양을 구하였다. 전체적으로 관전압의 크기기 클수록, 같은 관전압에서 물질 의 두께 두꺼워 질수록 그 차는 작아졌다. 같은 물질 두께의 차라 하더라도 관전압이 커질수록 물질의 두께의 차에 대한 차이는 작아졌다. 가장 큰 차는 50 kV의 0.1 mm에서 16.9%, 가장 작은 차는 120 kV의 8.6% 이였다(Table 6).
5)에너지 스펙트럼의 고 에너지 비교분석
고 에너지 부분의 왜도를 구함으로써 얼마나 높은 쪽에 에너지가 형성되는가 하는 실험에서, 관전압이 작고, 같은 관전압에서 물질의 두께가 두꺼울수록, 왜도 값은 더 크게 나타났다. 가장 왜도 값이 큰 것은 120 kV에서 단일 filter Cu0.1 mm, 가장 작은 것은 50 kV에서 복합 filter Cu0.3 mm이였다. 모든 값에서 단일 filter가 복합 filter에 비해 같 은 조건에서 더 우향화 되어 있음을 알 수 있었다(Table 7).
2.Image quality
1)SNR, PSNR, RMSE, MAE
필터의 두께가 증가됨에 따라 SNR, PSNR은 감소되며, RMSE, MAE는 증가된다. 복합 filter 실험에서도 같은 양 상으로 나타나며, 단일 filter에 비해 SNR, PSNR은 감소하 고 RMSE, MAE는 증가된다(Table 8).
2)CNR, SNR
Skull 팬텀을 이용한 실험에서 필터 물질의 두께가 두꺼 워 지면 관찰영역 3군데서 모두 SNR, CNR은 낮아진다. 대 부분의 필터에서 관찰영역 모두에서 단일 filter 경우에 CNR, SNR이 더 높았다(Table 9).
Pelvis 팬텀을 이용한 실험에서 필터 물질이 두꺼워지면 CNR, SNR은 모든 부분에서 감소하였고, 복합 filter 시에 도 모든 부분에서 감소하였고 특히, Sacrum 부분에 CNR, SNR값의 큰 차이를 보였다(Table 10).
3)정성적 평가
흉부 X선 검사 임상영상평가표(국가암검진 임상평가표) 를 기준으로 투과 상태, 해상도 및 대조도 평가 문항을 이용 하였다. 연구방법에서 항목은 7가지로 표의 표기를 쉽게 하 기 위하여 1~7번으로 명명하였다.
10년 이상 경력의 방사선사 5명이 평가한 결과 같은 두께 의 두 경우에서 결과 값의 차이가 크게 나타나지 않았다. Pair t-test 결과 모든 결과 값에 p>0.05로 나타나 두 값의 평균 차는 통계학적으로 의미가 없었다. 횡격막에 관련된 부분은 전체가 5점 만점을 기록하였고 전체적인 값들은 크 게 차이가 나지 않았다. 폐야의 전체적인 부분에서 폐혈관 의 관찰 부분과 심장 뒤의 폐혈관과 하행대동맥 관찰 부분 에서 Cu filter단독으로 사용한 경우에 약간 우수 하였다.
복합 filter 시 기관과 기관지 부분의 관찰에서 다소 평가 점수가 떨어지는 경향을 볼 수 있다. 또한 흉부에서 가장 두 꺼운 부분에 해당하는 흉추 부분에서 관찰 또한 점수가 전 체적으로 떨어졌다(Table 11).
IV.고 찰
X선관에서 방출되는 X선 에너지는 연속 spectrum으로 저에너지, 고에너지 영역이 혼재되어 있다. 일반적으로 부 가 filter를 사용하면 저에너지 영역을 많이 흡수하므로 상 대적으로 X선의 실효에너지를 올리는 데는 효과적이다. 그 러나 기존 Conventional radiography system는 얇은 부가 필터를 사용하거나, 관전압을 기존보다 감소시킬 수밖에 없 어, 선량을 충분히 감소시키지 못하는 한계가 있다8). 이는 X선질 경화에 의한 영상의 대조도가 감소되는 현상이 나타 나기 때문이다. 그러나 Flat panel detector system의 개 발과 발전을 통하여, 기존의 장비에 비하여 넓은 관용도 (latitude)와 높은 대조도(contrast resolution)의 특성이 제공되므로, 영상화질의 저하 없이 환자 선량을 감소시킬 수 있게 되었다. 디지털 방사선 장비는 아날로그 방식과 달 리 넓은 노출지수를 갖는다. 노출지수가 넓다는 것은 bit depth가 넓다는 것이고, 필름과 비교하여 검사조건에 대해 kV선택이 넓다는 것을 의미한다. 실제로 DR 장비에서 관전 압을 높이고 filter를 추가하면 촬영거리가 늘어날수록 표면 선량이 감소하는 선행 연구도 있다11).
에너지 스펙트럼을 살펴보면 저에너지 영역에서는 단일, 복합 fil ter의 비율을 살펴보면 단일은 0.39~3.4, 복합은 0.3~3.4 % 이었다. 그러나 광자양은 50 kV 0.1 mm에서는, 단일은 7,850, 복합은 4,870이며, 120 kV에서는 32,880, 30,660 으로 전체적으로 복합이 양적으로 적었으나 전체 광 자양 비율에서는 큰 차이를 나타내지 못했다.
중심 에너지 영역에서는 16.9~8.6 %로 단일 filter를 사용 한 경우가 효율이 더 놓았다. 고에너지 영역에서 더 높은 쪽 에 위치하여, 치우친 정도를 나타내는 왜도를 살펴보면 단일 filter에서는 0.323~2.166, 복합 filter에서는 0.295~2.157 로 단일 조건이 더 높은 에너지 영역에 위치함을 알게 되었 다. 이렇듯 에너지 영역을 나누어 살펴보았을 때 관전압과 filter 물질의 두께 따라서 선택적으로 사용이 가능할 것으 로 사료된다.
SNR은 DR에서 중요한 인자이고, 영상의 화질은 높은 SNR에 비례하여 향상된다. 또한 X선의 저선량에 다양한 변 수, 형광체, 그리고 DR 시스템의 전자적 요소에서 발생하는 양자noise는 SNR을 감소시킨다. DR 시스템은 검출기의 선 량률에 의존한다12). 그러나 SNR이 높더라도 다른 조직과 조직을 구분할 수 있는 CNR이 충분히 높지 않다면, 우수한 화질의 영상을 획득하는데 실패할 것이다. 즉, 본래부터 원 시 데이터나 측정된 신호 값들 사이의 차가 존재하지 않는 다면 각 물체의 조직 간의 대조도가 형성되지 않고, 대조도 를 증강시키는 방법 또한 어렵기 때문이다13).
PSNR은 대게 11 dB 이상이 되면 상당히 유사한 영상으로 보기는 하지만 사실상 30 dB을 기준으로 그 이하의 수치는 상대적으로 화질 수준이 떨어진다고 본다. 영상의 화질은 인간 의 육안 적으로 보고 판단하는 것이기 때문에, 굳이 PSNR이 높다고 좋은 영상이라 말하기는 어렵다14). Cu 0.1~0.3 mm에 서 PSNR은 45.54, 39.82, 36.57이며 RMSE는 43.54, 83.69, 121.55로 필터의 두께가 증가하면, PSNR은 낮아지고 RMSE 은 증가되었다. 이는 복합 filter 실험도 마찬가지였다. Cu 0.1~0.3는 PSNR은 42.54, 37.84, 35.24, RMSE는 43.28, 105.08, 141.58이었다. Al filter를 사용 전이 더 우수한 값을 나타내었으나, 그 차이는 크다고 할 수 없다. 그러나 단일 조건으로 만 실시한 실험으로 좀 더 다양한 관전압과 mAs에 대한 추가 실험이 필요할 것으로 사료되고, 또한 .기존에 임상 에서 잘 사용되어지지 않았던 0.2, 0.3 mm 또한 의미 있는 영상의 범주에 수치를 나타냄으로 제한적 경우, 즉 비만 환자, 밀도와 두께가 두꺼운 부분에서 사용을 고려가 가능하다.
CNR, SNR 실험에는 Frontal bone, Mandible, Orbit 모두 에서 Al filter 사용 전이 비해서 더 높았으며, 이는 Pelvis에서 도 마찬가지로 Iliac wing, Femur head, Sacrum관찰부위에 서는 모든 두께에서 Al filter사용 전이 우수하였다. 특이점은 Sacrum에서만 값의 차이가 규칙적으로 컸는데, filter 사용 시 해부학적 위치에 따라 결과 값이 달라질 수 있으므로 고려해 보아야 하는 점이라고 사료된다.
정성적 평가 부분에서는 통계적으로 유의성이 없는 것으 로 나왔는데, 같은 물질의 두께에서 투과도의 차이가 크지 않는 이유로 사료된다. 그러나 의미 있는 수치의 차는 아니 지만 폐혈관부의 관찰에서 Al filter 사용 전에 약간의 우수 성을 나타냈다. 이는 전체적인 에너지 스펙트럼에 양질의 광자가 많이 발생함에 따라, 전체적인 관용도의 향상으로 기인한다. 또한 두꺼운 부분 흉추등과 공기를 많이 포함하 고 있는 기관과 기관지 부분에서는 입상성이 거칠어져 두께 에 대한 고찰이 필요할 것으로 사료된다.
V.결 론
에너지 스펙트럼을 영역을 구분하여 분석한 결과 저에너 지 영역에서 광자양은 복합 filter가 적었으나, 전체 광자양 비율에서는 차이가 미미하였다. 또한 중심 에너지와 고에너 지 영역에서 단일 filter가 더 우수하였다.
PSNR, RMSE에서는 전체적으로 Al filter 사용 전 PSNR 은 더 높고, RMSE는 더 낮아 우수했다. 모든 결과 값에서 PSNR이 30 dB이상으로 수치화됨으로 영상 분석적으로, 영 상의 차이가 없는 것으로 나타났다. CNR, SNR은 부위마다 차이가 있었으나, Al filter 사용 전 우수하였다. 정성적 평 가에서는 전체적으로 Cu filter 단독 사용 시에 우수한 평가 는 받았으나, 통계학적으로 p>0.05로 유의하지는 않았다.
위의 결과 값들에서 알 수 있듯이 Al를 사용한 복합 부가 filter는 저에너지 영역에서 발생하는 광자양은 작으나 중심, 고에너지 영역에서는 효과적이 못했다. DR system의 보상 능력을 이용하면 적은 선량으로도 Al 복합 filter 보다 우수 한 에너지 스펙트럼의 이용이 가능하고 화질 평가에서 질 저 하 없이 효과적인 선량 관리가 가능할 것으로 사료된다.