바람직한 실시예의 상세한 설명
먼저 도면을 참조하면, 본 발명의 바람직한 실시예의 상세 구성이 도표 및 선도로서 표현되어있다. 각 도면에서 동일한 구성요소는 동일한 인용번호로 표시하고, 유사한 구성요소는 동일한 인용번호에 서로 다른 종류의 소문자 접미어를 붙여서 표시한다.
도1a 및 2b에 도시되어있는 바와 같이, 본 발명은 고에너지 실시간 직접 방사선 변환 X선 이미징 시스템(10)에 관한 것이다. 좀 더 구체적으로는 본 발명은 카드뮴-텔루르(Cd-Te)계 또는 카드뮴-아연-텔루르(Cd-Zn-Te)계 카메라를 이용하는 X선 이미징 시스템(10)에 관한 것이다. 본 발명의 실시간 X선 이미징 시스템(10)은 통상의 이미징 시스템과 마찬가지로 카메라 모듈(12)과, 이미지 프로세서(14)와, 표시 수단(16)을 포함한다. 본 발명의 실시간 X선 이미징 시스템(10)에 있어서 카메라 모듈(12)은 X선 이미징 장치(28)를 구비하고, 이 X선 이미징 장치(28)는 용도 특정 집적 회로(ASIC)의 판독 기판(32)과 전기적으로 연통하는 Cd-Te계 또는 Cd-Zn-Te계 방사선 검출기 기판(30)을 갖는다. 검출기 기판(30)의 각 활성 픽셀(36)은 ASIC 판독 기판(32) 상의 대응 픽셀 회로(31)에 전기적으로 접속되어있다.
다음으로 도1b를 참조하면, 이미징 시스템(10)은 프레임 그라버(frame grabber)(78) 및 이미징 소프트웨어(80)를 내장한 퍼스널 컴퓨터(PC)(76)를 구비하며, 이 퍼스널 컴퓨터(76)는 카메라 링크(82)를 거쳐서 카메라(37), 전원(86) 및 AC/DC 어댑터(88)를 구비한 X선 유닛(84)에 연결된다. 일반적으로 X선 유닛(84)은 예를 들면 네트워크 내의 퍼스널 컴퓨터나 단자에 접속을 하기 위한 네트워크 접속부(90)를 더 구비한다.
다음으로 도2a를 참조하면, 카메라(37)는 인터페이스 인쇄회로기판(PCB)(92)을 갖추고 있으며, 이 인쇄회로기판(92)은 데이터 버스(94)를 통해서 검출기 인쇄회로기판(96)에 접속되는 것으로, 냉각 요소(98)를 갖는다.
X선 이미징 장치(28)는 다수의 이미지 프레임(44)을 생성할 수 있고, 각 이미지 프레임(44)은 미보정 이미지 픽셀 값의 어레이(45)로 이루어진다.
다음으로 도2b를 참조하면, 본 발명의 이미징 시스템(28)의 카메라 모듈(12)에 유용한 이미징 장치(28)가 개략도로서 도시되어있다. 도2b에 개괄적으로 예시한 이 이미징 장치(28)에 있어서 검출기 반도체 기판(30)은 ASIC 판독 기판(32)과의 전기 접속부(35)(예를 들면, 도시된 바람직한 실시예의 범프 접합부)를 갖는다. 반도체 기판(30)의 검출기 소재(34)는 본 발명에 있어서 카드뮴-텔루르계 또는 카드뮴-아연-텔루르계 조성물로 이루어지며 입사 방사선을 흡수한다. 이러한 흡수에 따 라 방사선 에너지는 검출기 소재(34)의 두께 내에서 직접 전하로 변환된다. 이 전하는 각 활성 픽셀 또는 기능 픽셀(36)의 검출기 픽셀 수집 전극(픽셀 접점)(38)에 모이고 전기 접속부(35)를 통해서 ASIC 판독 기판(32)의 픽셀 회로(31) 상의 픽셀 회로 접점(33)과 전기적으로 연통한다. 이렇게 생성된 전하 신호는 ASIC 판독 기판(32)의 검출기 픽셀의 대응 픽셀 회로(31)에서 저장 및/또는 처리된다. 그 후 ASIC 픽셀 회로(31)는 통상적으로 멀티플렉스되고, 이어서 아날로그 출력이 생성되어 칩 내에서 혹은 칩 외에서 디지털화된다. 본 발명에 따르면, 이미지 프레임(44)의 각 픽셀 값(36)은 디지털화되고 해당 프레임의 픽셀 값에 맞춘 개별적 캘리브레이션 보정치를 더 갖는다. 따라서 각 픽셀 값(36)은 실제로는 보정된 디지털 픽셀 값이다. 각 이미지 픽셀(47)에 맞춰진 보정치는 보정 캘리브레이션 프로세스에 따라 보정된 복수의 개별적 단일 프레임 픽셀 값(36)으로부터 유도된다. 서로 다른 이미지 프레임(44)으로부터의 동일한 특정 이미지 픽셀(47)의 보정된 개별적 디지털 픽셀 값(36)은 수집된 이미지 프레임(44)의 적어도 일부에 대하여 정규화 모듈(24)의 알고리즘에 따라 처리함으로써 최종 이미지에 표시될 픽셀 값을 얻는다. 최종 이미지는 실시간 이미지일 수도 있고 정적인 디지털 축적 이미지일 수도 있다.
다음으로 도2c를 참조하면, 본 발명의 이미징 장치(28)는 다수의 이미지 프레임(44)을 생성할 수 있으며, 각 이미지 프레임은 소정의 비트 깊이(즉, 개별 픽셀의 컬러 스케일 또는 그레이 스케일)를 갖는 프레임 픽셀 값(36)의 어레이(45)를 포함한다. (여기서, 컬러 당 8비트의 픽셀은 24 비트의 이미지를 제공하는데, 그 이유는 8 비트x3 컬러=24 비트이기 때문이다. 예를 들어 24 비트 컬러 해상도는 16,700,000 컬러이다.) 이미징 시스템(10)은 서로 다른 프레임(44)의 픽셀 값(36)으로부터 이미지 픽셀 값(47)을 산출하기 위한 처리 수단(24)을 구비하는 바, 이 때 이미지 픽셀 값(47)의 비트 깊이는 개별 프레임(44)로부터의 픽셀 값(36)의 비트 깊이보다 크다. 예를 들면, 단일 프레임을 디지털화하면 12 비트에 불과하며 0 내지 최대 4096로 된다. 이러한 아날로그-디지털 변환기(ADC)는 매우 흔한 것으로 요즈음에는 가격도 저렴하다. 또한, 아날로그-디지털 변환기는 처리 속도가 상당히 빠르며 5 MHz, 심지어는 10 MHz-20 MHz의 클록 속도(clock rates)로 작동한다. 본 발명의 양수인이 구현한 통상의 CdTe-CMOS 카메라는 10 k 내지 1 M 픽셀을 포함할 수 있다. 이것은 단일의 ADC로 20 fps-300 fps, 심지어는 2,000 fps의 프레임률을 실현할 수 있다는 의미이다. 이미지 프레임(44)의 판독 후, 미보정 픽셀 값은 디지털화된다. 따라서 본 발명에 의하면 픽셀 보정 알고리즘(20)을 이용해서 산출한 보정치를 디지털 픽셀 값에 적용함으로써 단일 프레임으로부터 보정된 프레임의 디지털 픽셀 값(36)을 얻을 수 있다. 이어서, 도3a에 나타내는 바와 같이, 서로 다른 프레임(44)으로부터의 디지털 보정 픽셀 값(36)을 축적함으로써 12 비트보다 훨씬 큰 비트 해상도로 표시될 이미지의 디지털 보정 픽셀 값(47)을 얻을 수 있다. 예를 들면, 12 비트 해상도의 프레임이 17개 있다면, 각 프레임은 디지털 축적 후 16 비트 이상이 될 수 있다 (17x4096=69632>16 비트). 이는 디지털 X선 이미징 분야에서 실로 놀라운 기술 혁신이 아닐 수 없다. 그 이유는 종래에 이 정도의 해상도를 얻으려면 긴 집적 시간을 필요로 할 뿐만 아니라 고가이면서 속도가 느린 16 비트 ADC를 사용해야하는 관계로 이미지를 실시간으로 표시할 수 없었기 때문이다. 또한, 위에서도 설명했고 나중에도 추가로 설명하겠지만, 아날로그 신호의 집적 시간이 길어지면 암 전류가 증가하는 등의 다른 문제점이 나타나고 다른 종류의 노이즈가 발생하게 된다.
다음으로 도3b를 참조하면, 서로 다른 프레임(44)으로부터의 보정된 디지털 픽셀 값(36)은 대응하는 픽셀 값일 수도 있고 (스캐닝의 경우처럼) 프레임(44) 내의 다른 위치에서 온 것일 수도 있다. 본질적으로 본 발명의 이미징 시스템(10)은 서로 다른 이미지 프레임(44)으로부터의 이미지 픽셀 값을 보정하기 위한 방법 (20, 49)과 이미지의 보정된 픽셀 값(47)을 산출하기 위한 처리 방법(24)을 포함하며, 이 처리 방법은 넓은 의미에서 서로 대응하는 여러 개의 프레임으로부터 온 보정된 디지털 픽셀 값을 이용한다.
위에서 설명한 바와 같이, 본 발명의 이미징 장치(10)는 바람직하게는 CdTe계 또는 CdZnTe계 X선/감마선 이미징 장치(28)를 포함하며, 이에 따라 CdTe/CdZnTe로 픽셀화된 검출기 기판(30)이 적어도 하나의 ASIC 판독 기판(32)에 범프 접합된다. 이 CdTe/CdZnTe 검출기 기판은 입사하는 X선 또는 감마선을 전자 신호로 직접 변환하기 위한 것이고, ASIC 판독 기판은 각 픽셀(36)로부터의 전자 신호를 저장 및/또는 처리하고 그 신호를 판독하기 위한 것이다. CdTe/CdZnTe 이미징 장치(28)는 바람직하게는 10 fps 더욱 바람직하게는 25 fps-100 fps, 경우에 따라서는 300 fps 이상의 높은 프레임률로 판독된다. 개별 프레임(44)은 각 프레임이 일련의 픽셀 값(36)으로 표현될 수 있도록 디지털화된다. 여기서 각 픽셀 값은 이미징 장 치(28)에 의해 생성된 프레임 내의 특정 픽셀(36)에 대한 디지털 신호 레벨에 대응한다. 각 프레임(44)에 대하여 디지털화된 픽셀 값은 위에서 이미 설명한 픽셀 보정 알고리즘(49)에 따라 보정된다. 동일한 픽셀(36)에 대응하는 서로 다른 프레임(44)의 보정된 개별 디지털 픽셀 값(36)은 수집된 프레임의 적어도 일부에 걸쳐서 소정의 알고리즘에 따라 가산되거나 평균화되거나 처리됨으로써 최종 이미지에 표시될 픽셀 값(47)으로 된다.
본 발명의 핵심은 CdTe 또는 CdZnTe 결정의 결함, CMOS의 비선형성 및 오프셋, 암 전류, 분극화 그리고 후술하는 기타 영향을 모두 고려해야함에도 불구하고 각 개별 프레임(44)으로부터의 디지털 픽셀 값에 대한 보정을 실제로 행할 수 있다는 점이다. 다음 장에서는 픽셀 보정을 실시간으로 효율 좋게 행하는 것에 관해서 설명한다. 본 발명에서는 발명의 범위를 변화시키거나 그 범위로부터 벗어나지 않으면서도 여러 가지의 보정 및/또는 캘리브레이션 기법을 이용할 수 있다.
카메라 모듈(12)과 고속 프레임 프로세서 모듈(18)은 케이블 링크(60)를 통해서 서로 통신한다. 카메라 모듈(12)은 각 픽셀 (36) 또는 픽셀 셀(pixel cell)(29)의 개별적 미처리 픽셀 회로 출력을 나타내는 처리 및 체계화된 픽셀 데이터를 프레임 프로세서 모듈(18)로 공급한다. 이러한 고속 프레임 프로세서 모듈(18)은 본 기술분야에서 카메라 모듈(12)로부터의 픽셀 회로 데이터를 획득하는데 일반적으로 사용되는 프레임 그라버 (78)용 회로를 구비하며 [경우에 따라 프레임 그라버(78)는 카메라 모듈(12)의 일부일 수도 있음], 픽셀 회로 데이터를 처리함으로써 각 픽셀 셀(29)의 미처리 픽셀 회로 출력을 나타내는 미처리 상태의 타임 스탬프(time stamp) 이미지 프레임을 제공한다. 이어서 프레임 프로세서는 이미징 시스템(10)이 캘리브레이션 모드일 경우에는 미처리 상태의 타임스탬프 이미지 프레임을 프레임 데이터 링크(66)를 통해서 캘리브레이션 모듈(20)로 보내고 캘리브레이션 모드가 아닐 경우에는 정규화 모듈(24)로 보낸다.
캘리브레이션 모듈(20)은 캘리브레이션 프로세스(49)를 제어한다. 캘리브레이션 프로세스(49)는 미처리 상태의 타임스탬프 이미지 프레임 데이터 및 기준 필드 방사 강도 등과 같은 기타 캘리브레이션 파라미터를 분석하고 캘리브레이션 데이터 구조 모듈(22)의 룩업 테이블을 탑재하는데 필요한 데이터를 생성한다. 캘리브레이션 모듈(20)은 데이터베이스 링크(68)를 통해서 데이터 구조에 데이터를 기입한다. 만약 룩업 테이블에 적절한 캘리브레이션 데이터를 탑재하지 않으면 정규화 모듈(24)로부터 표시 모듈로의 모든 이미지 출력은 정확성을 잃게 된다. 따라서 본 발명의 이미징 시스템(10)이 정상적인 이미징 동작에 들어가기에 앞서 캘리브레이션 프로세스(49)를 실행하는 것이 필요하다.
캘리브레이션 모드가 아닐 경우 프레임 프로세서(18)는 이미지 프레임(44)의 타임스탬프 데이터를 정규화 모듈(24)로 보낸다. 정규화 모듈(24)은 제2 데이터베이스 링크(70)를 통해 룩업 테이블로부터 추출한 이미지 픽셀의 대응 보정 요구 조건에 따라 미처리 상태의 타임스탬프 이미지 프레임의 각 이미지 픽셀에 대한 정규화 동작을 실행한다. 이어서 정규화 모듈(24)은 정규화된 이미지 프레임을 표시 데이터 링크(74)를 거쳐 표시 모듈(16)로 공급한다. 정규화된 이미지 프레임의 모든 이미지 픽셀은 각각에 대응하는 미처리 이미지 픽셀 강도 값을 표시하며 이 강도 값은 각각에 대응하는 룩업 테이블로부터의 보정 계수에 의해 보정된 것이다.
고품질의 이미지를 얻기 위해서는 카드뮴-텔루르계 검출기 기판(30)과 관련해서 몇 가지의 장애를 극복해야한다. 예를 들면, 누설 전류(암 전류라고도 불림)가 계속해서 발생하며 이에 대한 보상이 필요하다. Cd-Te계 또는 Cd-Zn-Te계 검출기 소재(34) 중 어떤 것은 누설 전류의 레벨을 조절하기 위한 블로킹 접점(도시하지 않음)을 갖도록 제조된다. 검출기 소재(34) 중 다른 어떤 것은 누설 전류를 억제하도록 Zn이나 기타 도판트(dopant)의 함량을 달리해서 제조한다. 어쨌든 누설 전류는 노이즈를 발생시킬 뿐만 아니라 각 픽셀 회로(31) 상의 전하 수집 게이트(33)를 폐색시키게 된다. 또한 블로킹 접점을 사용하면 분극화나 전하 트래핑의 문제점을 발생시키게 되는데, 이러한 문제점은 장치에 따라 다르기는 하지만 작동을 시작하고 나서 수 초, 예를 들면 5초, 10초 또는 60초 후에 현저하게 나타난다.
본 발명의 검출기 기판(30)에서 카드뮴-텔루르계 조성물(즉, Cd-Te 및 Cd-Zn-Te)을 방사선 흡수 매체(34)로 사용할 때의 장점은 매우 높은 방사선 흡수 효율, 잔광(afterglow)의 최소화 그리고 높은 이미지 해상도의 실현 가능성이다. 이러한 장점은 위에서 언급한 문제점을 완화하거나 해소할 수 있는 이미징 장치를 구현하는데 있어 중요하다. 블로킹 접점을 사용하지 않는 경우라고 하더라도, ASIC 판독 기판(32)의 각 픽셀 회로(31)에 있어서 전하 저장 콘덴서의 용량을 증가시키지 않는 이상, 100 msec 이상의 장시간 노출은 누설 전류 및 결정 결함의 문제 때문에 곤란하다. 그러나 전하 저장 콘덴서의 크기를 증가시키는 것은 감도 면에서 불리하다. 그 이유는 전하 저장 용량이 커지는 만큼 감도가 저하하기 때문이다. 일 예이기는 하지만, 본 발명은 전하를 수용하는 각 ASIC 픽셀 회로의 전하 저장 콘덴서로서 50fF용량의 것을 사용함으로써 성공적으로 구현되었다. 이러한 용량의 콘덴서를 사용하면 실용 가능한 최대 노출 시간을 100 msec 이하로 줄여 Cd-Te 또는 Cd-Zn-Te 매체에 있어서 누설 전류 및 기타 결함의 발생을 억제할 수 있다.
다음으로 도4a 및 도4b를 참조하면, 직접 변환(전하 결합) 방사선 검출기 장치에 있어서 과도한 분극화를 방지하는데 매우 유용한 메커니즘을 이용해서 고전압 바이어스의 온/오프 사이클을 단축시킨다. 이 기법을 검출기 바이어스 전압 스위칭 기법이라 부르는데, 이에 의하면 데이터 수집 사이클의 말미에 검출기 기판의 바이어스 전압을 짧은 주기(100 msec 미만)로 스위칭 오프하게 된다. 데이터 수집 사이클의 시간은 선택 가능하며, 예를 들면 3초 내지 20초 이상으로 할 수 있다. 바이어스 전압을 스위칭하면 검출기 기판(30)에서 분극화 또는 전하 트래핑이 발생하는 것을 방지할 수 있다. 그러나 바이어스 전압 스위칭 기법은 X선 이미징 시스템 분야에서 새로운 것이며, 잘 다루지 않을 경우 이미지 품질에 악영향을 줄 수도 있는 몇 가지 속성을 가지고 있다. 그러한 속성 중 하나는 "데드 타임(Dead-time)"이고 다른 하나는 "픽셀 응답 드리프트"이다. "데드 타임"이라 함은 검출기 바이어스 전압이 오프되어 검출기 전하가 수집되지 않는 데이터 수집 사이클의 기간을 말한다. "픽셀 응답 드리프트"는 검출기 바이어스 전압을 다시 온시켰을 때 나타나는 결과로서, 정적 방사선 필드에 대한 픽셀의 응답성이 아직 안정화되지 않은 데이터 수집 사이클의 초기에 해당한다. 이와 같은 2가지의 제약 사항이 도4a에 도시되어있다. 도4b에는 검출기 기판 바이어스 스위칭 회로(121)가 도시되어있다.
도4a에 도시한 실시예에 있어서, 데이터 수집 사이클 타임 Ct는 검출기 바이어스 전압 온/오프 펄스(50)의 발생 개시 시점간의 간격을 말한다. 데드 타임 Dt는 실제의 고전압 다운 시간 Vo과 고전압이 다시 스위치 온된 후에 필요한 얼마간의 안정화 시간을 합산한 것이다. 데드 타임 Dt는 고전압 다운 시간 Vo보다 작을 수 없으므로 검출기 바이어스 전압 스위칭 방식의 이미징 시스템에 있어서 데드 타임 Dt의 영향을 완전히 없애는 것은 불가능하다. 그러나 분극화(전하 트래핑)를 감소시키고/시키거나 데드타임을 데이터 수집 사이클에 있어서 무시할 수 있을 만큼 작게 유지할 수 있도록 바이어스 전압의 오프 시간을 단축함으로써 데드 타임 Dt을 부분적으로 최소화할 수는 있다.
바이어스 전압 스위칭 방식 검출기의 다른 한 가지 잠재적 제약 사항은 픽셀 응답 드리프트 Rd인데, 이는 정적 방사선 필드 노출 레벨에 대하여 픽셀 회로의 경시적 출력 신호(40)가 비선형성을 갖는 것과 관련이 있다(도4a 참조). 이러한 비선형성은 전압 오프-온 펄스(50)의 전압 온 스텝 직후에 가장 현저하게 나타난다. 이러한 비선형성을 바로 잡지 않으면 실시간 이미지 표시기에 있어서 이미지의 전반적 휘도 레벨을 급 등락시키는 원인이 된다. 바이어스 전압 스위칭 방식의 이미징 장치에 있어서 픽셀 셀의 응답이 비선형인 것은 본 발명의 이미징 시스템의 이미지 프레임 생성 후 캘리브레이션 방법을 적용하여 위와 같은 실시간 X선 이미지 표시기의 강도 왜곡을 제거하기에 아주 좋은 대상이다.
다음으로 도9를 참조하면, 본 발명의 캘리브레이션 방법(49)은 암 전류 강도를 비롯한 여러 가지의 서로 다른 균질의 기준 방사선 필드 강도에서 완전한 데이 터 수집 사이클 동안 캘리브레이션 데이터를 수집한다. 그러므로 이 캘리브레이션 방법(49)은 검출기 바이어스 전압 스위칭 기법을 이용하는 디지털 이미징 시스템에서 특히 실용성이 있다. 검출기 바이어스 전압 스위칭 기법을 이용하는 디지털 이미징 시스템의 카메라 모듈(12)은 일반적으로 수천 개의 픽셀 셀(29)을 갖는 검출기/ASIC 조립체(28)를 포함하며, 각각의 픽셀 셀(29)은 검출기 픽셀(36)과 관련 픽셀 회로(31)를 포함한다. 각 픽셀 회로(31)는 해당 픽셀 회로(31)의 디지털화 픽셀 신호를 생성하는 관련 회로 부품과 픽셀 회로 신호 출력부(도시하지 않음)를 구비한다. 픽셀 회로 출력 신호는 관련 검출기 픽셀(36)에 입사하는 X선/감마선 방사 에너지의 강도를 표시한다(도2b 참조). 또한, 무작위로 발생하는 노이즈를 감소시키기 위하여 각각의 기준 강도에 대하여 사이클을 반복한다.
이렇게 수집된 디지털화 픽셀 신호 출력은 카메라 링크(60)를 거쳐서 이미지 프로세서(14)의 고속 프레임 프로세서 모듈(18)로 보내진다. 이 프레임 프로세서 모듈(18)은 각 픽셀 회로(31)로부터의 개별 픽셀 회로 출력 신호를 수신하는 프레임 그라버 회로를 구비한다. 프레임 프로세서 모듈(18)은 개별 디지털화 픽셀 신호를 체계화하여 이미지 프레임을 만든다. 이 때 이미지 프레임의 각 이미지 픽셀은 카메라 모듈(12)의 이미징 장치(28)에 있어서 픽셀 회로에 대응하는 픽셀 신호를 나타낸다. 이미지 프레임에서 이미지 픽셀의 강도는 대응 픽셀 회로(31)로부터 수신된 픽셀 신호의 세기를 나타낸다. 그러나 각 픽셀 셀(29)을 이루는 개별 성분의 기계적 특성과 전기적 특성은 근본적으로 다르기 때문에 이미지 프레임을 포함하는 각종 픽셀의 필드 강도 응답성은 균일한 X선 필드에 대해서조차도 일률적이지 않 다. 따라서 이미지 프레임으로 표시되는 정보를 사용자가 이용하기 전에 이미징 장치를 캘리브레이션할 필요가 있다.
캘리브레이션 과정
다음으로 도8을 참조하면, 캘리브레이션 과정(49)이 상위 개념의 흐름도로서 도시되어있는 바, 여기서 보정 함수(110)에는 픽셀 값(36)을 입력하고 보정 계수(120)를 적용하여 보정된 픽셀 값(47)을 얻는다. 도9는 본 발명의 이미징 시스템(10)에 있어서 캘리브레이션 프로세스(49)의 단계를 좀 더 상세하게 보여주고 있다. 캘리브레이션 데이터는 완전한 데이터 수집 사이클 동안 암 전류 강도 ID를 비롯한 여러 가지의 서로 다른 균질의 기준 방사선 필드 강도로 수집된다. 제1 단계(49a)에서는 암 전류 강도 ID에 대하여 데이터가 수집된다. 제2 단계(49b)에서는 X선 강도 I1에 대하여 데이터가 수집된다. 제3 단계(49c)에서는 X선 강도 IN에 대하여 데이터가 수집된다. 제4 단계(49d)에서는 보정 계수가 산출되고 룩업 테이블(22)이 작성된다.
다음으로 도10 내지 도12를 참조하면, 캘리브레이션 과정(49)이 더욱 상세하게 설명되어있다. 도10에는 데이터 수집 사이클의 하위 개념 방법(120)이 설명되어있다. 제1 단계(122)에서는 데이터 빈(Data Bin)이 초기화되며, 이 데이터 빈은 다음과 같은 구조를 갖는다.
빈 1: 시간 0....T1
빈 2: 시간 T1....T2
빈 N: 시간 TN-1....TN
제2 단계(124)에서는 방사선 필드 강도가 설정된다. 제3 단계(126)에서는 고전압 펄스가 발생된다. 제4 단계(130)에서는 타이머가 리세트된다. 제5 단계(132)에서는 수집 시간이 이미지 프레임의 시간 TIF과 같게 설정되고 데이터 빈 B이 탐색되고 프레임(44)이 데이터 빈 B에 부가되며 사이클이 활성화되어 있는 한 이 루프는 계속 실행된다. 이어서 제7 단계(134)에서는 노이즈 감소를 위하여 위의 동작을 더 반복할 필요가 있는 경우 예를 들면 위의 루프를 다시 실행한다.
도11에 있어서 보정 계수를 산출하기 위한 하위개념 방법(120)은 다음의 단계를 포함한다. 각각의 빈에 대한 루프(140)와 이 루프에 속하는 각 픽셀(36)에 대한 루프(142)의 제1 단계(144)에서는 모든 강도 값 ID 및 I0....IN에 맞춰 다항식을 만들고 픽셀이 한계치 테스트를 통과하지 못하면 해당 픽셀(36)을 플래그한다. 제2 단계(146)에서는 데이터 구조(룩업 테이블)에 기입한다. 제3 단계(148)에서는 하위개념의 방법(120)을 마스킹 루틴(150)으로 이행하여 계속 실행한다.
도12에 있어서 마스킹 하위개념 방법(150)은 다음의 단계를 포함한다. 제1 단계(152)에서 하위개념 방법(150)은 플래그된 픽셀(36)을 체크한다. 제2 단계(154)에서 하위개념 방법(150)은 각각의 플래그된 픽셀(36)과 관련하여 양호한 인접 픽셀(36)을 탐색한다. 제3 단계(156)에서는 이렇게 탐색된 양호한 인접 픽셀의 위치를 데이터 구조(20)에 기입한다. 플래그된 픽셀(36)이 더 이상 없으면 하위 개념 방법(150)을 종료한다.
도13에는 정규화 과정(160)이 설명되어있는 바, 이 정규화 과정은 다음의 단계를 포함하고, 이들 단계는 프레임 프로세서 모듈로부터의 미처리 이미지 픽셀 데이터에 대하여 실행된다. 제1 단계(162)에서는 이미징 시스템(10)의 동작 중 이미지 프레임(44)이 수신되고 타임스탬프가 "T"로 설정된다. 제2 단계(164)에서는 시간 "T" 동안 빈이 탐색된다. 제3 단계(166)에서는 각 픽셀(36)을 반복적으로 체크하여 플래그된 픽셀 또는 불량 픽셀인지의 여부를 판단한다. 긍정으로 판단되면, 제4 단계(168)에서 픽셀 값(47)을 양호한 인접 픽셀들의 가중 평균치로 대체한다. 제5 단계에서는 보정 다항식을 적용한다.
프레임 프로세서 모듈로부터 미처리 이미지 픽셀 데이터가 공급됨에 따라, 캘리브레이션 프로세스에서는 소프트웨어로 구동되는 캘리브레이션 모듈(20)을 이용해서 "룩업 테이블"을 생성하고 이것을 데이터 구조 모듈(22)에 유지한다. 룩업 테이블은 이미지 프레임의 각 픽셀에 대한 시간 의존형 이미지 픽셀 별 보정(pixel specific correction)계수(54)의 집합이다. 픽셀 별 보정계수(54)는 목표 균일 강도 값(52)을 기준으로 하며(도6 참조), 해당 이미지 픽셀의 미처리강도 값을 정규화된 값으로 보정하는데 이용된다. 따라서 이미지 프레임에 표시된 각 이미지 픽셀은 복수의 기준 X선 필드 강도마다 각각 생성된 데이터 구조 모듈(22)의 룩업 테이블 내에 있어서 시간 의존형 보정계수의 데이터 집합을 갖는다.
보정계수/보정치 집합의 시간 의존성은 고속 프레임 모듈에 의해 처리되는 각 이미지 프레임에 타임스탬프를 적용함에 따라 생겨난다. 타임스탬프라 함은 이 미지 프레임(44)이 생성되는 데이터 수집 사이클 Ct의 개시 후 경과한 시간을 말한다. 도5에 도시된 바람직한 실시예에 있어서 타임스탬프를 적용한 이미지 프레임(44)은 데이터 수집 사이클 Ct 동안 균일한 프레임 간격(46)으로 카메라 모듈(12)로부터 캡처(획득)된다. 따라서 타임스탬프를 적용한 이미지 프레임(44)은 항상 서로에 대하여 동일한 시간차를 갖는다. 검출기 바이어스 전압이 스위칭 온된 후 획득된 제1 프레임에는 타임스탬프 0을 할당하고, 제2 프레임에는 타임스탬프 1을 할당하며, 타임스탬프 N까지 이러한 방식으로 할당을 행한다. 실제의 경우 개별적 캘리브레이션 데이터 집합은 각각의 이미지 픽셀에 대하여 산출되는 것으로, 데이터 수집 사이클 Ct에서의 각 타임 스탬프시마다 해당 이미지 픽셀에 대한 보정치가 포함된다. 이와 달리, 캘리브레이션 데이터는 데이터 수집 사이클 Ct당 하나씩 N개의 서로 다른 캘리브레이션 데이터 집합으로 이루어진 것으로 간주하거나 그렇게 구성할 수도 있으며, 이 때 각각의 프레임 데이터 집합은 프레임 내의 각 이미지 픽셀에 대한 개별적 보정치/보정계수를 포함한다. 최고의 화질을 얻기 위해서는 상기 "N"을 여러 가지의 타임스탬프 중 가능한 가장 높은 수 Nmax. 바꾸어 말하면 가능한 가장 높은 프레임률이 되도록 선택해야한다. 그러나 이는 극도로 많은 데이터를 필요로 하는 조건이며, 현재로서는 기술적 제약, 즉 컴퓨터 메모리 처리 시간의 제약을 받기 때문에 N<Nmax가 되도록 선택할 수밖에 없다.
다음으로 도6을 참조하면, 그래프(180)는 각 이미지 프레임(44)의 특정 이미지 픽셀 출력의 미처리 강도 값에 소정의 시간 의존형 보정계수를 적용함으로써 이 미지 픽셀(47)의 강도 값을 정규화하는 것을 보여주고 있다.
데이터의 수집: 캘리브레이션 방법의 제1 단계는 관련 데이터, 구체적으로는 여러 가지의 기준 방사선 필드 강도에 대한 카메라 이미징 장치(28)의 응답 결과를 수집하는 것이다. 이미징 장치(28)의 각 픽셀 셀(29)의 응답 결과는 데이터 수집 사이클 Ct 중의 모든 타임스탬프에 대하여 수집한다. 도면에 도시한 바람직한 실시예에 있어서, 이 단계는 착신되는 퀀텀(quantum) 노이즈의 영향을 감소시킬 수 있도록 1회 이상(대체적으로 20회 이상) 반복한다. 이러한 방식으로 관련 데이터를 수집하면 검출기 또는 ASIC 부품의 불균일성을 보정할 수 있을 뿐만 아니라 기본적으로 "플랫 필드(flat-field)" 보정이 가능해진다. 본 실시예에 있어서, 캘리브레이션 방법은 카메라 모듈(12)의 이미징 장치(28)가 방사선 공급원과 소정의 기하학적 관계에 있음을 전제로 한다. 이 말은 방사선 공급원이 달라지거나 이미징 장치(28)와 방사선 공급원 간의 배치관계가 변화할 때마다 캘리브레이션을 다시 실시해야한다는 뜻이다. 또한 캘리브레이션은 사용된 방사선 스펙트럼 별로 반복해야한다.
픽셀 별 보정(pixel specific correction)계수/보정치의 산출: 여러 가지의 기준 방사선 필드 강도로 노출을 행하였을 때 시간의 함수로서의 단일 픽셀 셀(29)의 응답 결과는 특징적 모습을 갖는다. 본 발명의 캘리브레이션 방법에 있어서 기본이 되는 개념은 균일성이다. 모든 픽셀 셀(29)은 동일한 강도의 방사선에 노출되었을 때 동일한 픽셀 출력 신호를 발생해야한다. 이는 다음의 캘리브레이션 함수가 픽셀 출력 값 xin 로부터 총 출력 값 yout 으로 맵핑된 것임을 의미한다.
이와 같이 맵핑을 하는 것은 다른 모든 픽셀과 동일한 출력을 내는 각 픽셀에 적합한 함수 f pix ()를 구하기 위함이다.
다항식을 이용하는 이유는 극히 빠른 계산을 할 수 있고 이는 실시간 동작을 위해 반드시 필요하기 때문이다. 다항식은 예상치 못한 보간(interpolation) 및 외삽(extrapolation)의 거동을 나타내기 때문에 이와 같은 회귀 문제의 가장 확실한 해법이 되지는 못한다. 이제 위의 함수 fpix ()는 분명히 다음과 같이 표현할 수 있다.
식 중, a i,pix 는 픽셀 pix에 대한 계수이고 M은 다항식의 차수이다. 통상적으로 사용되는 선형 캘리브레이션(이득 및 오프셋 보정)은 M=1일 때 특히 유효하다. 본 실시예에서는 기본적으로 3차까지의 다항식을 이용하지만, 충분히 많은 수의 시간 의존형 계수 데이터 집합을 이용할 경우에는 선형 보정만으로 충분할 수도 있다.
캘리브레이션 파라미터의 추정: 이와 같은 회귀 문제에서 모델 파라미터를 추정하는 일반적 방법은 최대 우도 (Maximum Likelihood/ML) 추정법이다. 이는 함수 및 노이즈 모델이 주어졌다는 가정 하에 특정 시점에서 하나의 픽셀에 대한 모 든 데이터 포인트의 가능성을 최대화할 수 있음을 의미한다. 정규 분포의 제로 평균 노이즈를 가정하며, 하나의 데이터 샘플 x i 의 확률은 다음과 같다.
그리고 모든 샘플이 통계적으로 독립적이라는 가정 하에 해당 샘플의 총 우도는 다음과 같다.
최대 우도 추정법의 한 가지 문제점은 이미 알고 있는 지식을 정확히 적용하기가 매우 어렵다는 점이다. 이를 극복하기 위해서 최대 사후 (Maximum A Posteriori/MAP) 추정법을 이용한다. 이 최대 사후 추정법에서는 모든 샘플의 사후 분포가 다음 식과 같이 최대화된다.
식 중, A는 독립성을 전제로 추정된 샘플의 공분산 행렬(covariance matrix)이고,
및
x=[x
1...x
Ndata]는 데이터 샘플의 벡터이며,
f=[f(x
1)...f(x
Ndata)]는 이 픽셀에 대한 캘리브레이션 값의 벡터이다. p(x)는 상관성이 없는 스케일링 팩터(factor), 즉 에비던스(evidence)이다. 노이즈 및 함수 파라미터 프라이어에 대한 정규 분포가 다음과 같다고 가정하면,
최종 사후 추정치는 다음의 식으로 표현된다.
상기 식의 자연지수를 취하여 모든 상수 계수를 하나로 묶으면 다음과 같은 비용 함수가 얻어진다.
이 비용 함수는
의 페널티 파라미터를 갖는 가중치 및 제약이 부여된 선형 최소 자승 비용 함수로 해석할 수 있다. 최종 파라미터 값은 위의 방정식을 모든 함수 파라미터
a i 에 대하여 미분하고 이어서 미분계수를 0으로 설정함으로써 풀 수 있다. 가중치가 부여된 최소 자승을 이용하는 이유는 서로 다른 X선 강도를 이용할 경우 가장 큰 강도에 대한 퀀텀 노이즈가 예를 들면 암 전류에 대한 퀀텀 노이즈보다도 훨씬 커지기 때문이다. 이에 의하면 상대적으로 작은 값에 상대적으로 큰 가중치를 부여할 수 있어 정확도를 향상시킬 수 있다.
실시 및 성능에 관한 고려사항: 모든 계산에서는 화질을 최적화하기 위해 32 비트 플로팅 포인트 계산법을 이용한다. 요즈음의 86배속 프로세서는 매우 효율적인 병렬 처리를 가능하게 하는 SIMD(단일 명령 복수 데이터) 지령을 잘 실행한다.
보정 프로토콜에 사용하기에 적합한 타임스탬프 캘리브레이션 이미지 프레임 의 선택: 실용상의 이유에서 데이터 수집 사이클 Ct의 모든 타임스탬프를 사용하는 것은 불가능하다. 그 이유는 생성되는 데이터의 양이 막대하고 경우에 따라서는 처리 시간 및 메모리의 할당이 곤란하기 때문이다. 이는 요즈음의 대(large) 면적 카메라의 경우 이미지 크기가 508x512 픽셀에 이르는데서 연유한다. 단위 픽셀 당 파라미터는 4개에 이르고(3차 다항식을 이용할 경우), 각 파라미터는 4 byte이다. 이는 단위 프레임 당 수집되는 데이터의 크기가 3.97 MB라는 것을 의미한다. 본 실시예에 있어서 카메라는 초당 50 프레임의 속도로 동작하는데, 이는 데이터 수집 속도가 198 MB/sec라는 의미이다. 더욱이 이미지는 PCI 버스에서 16 비트 포맷(24.8 MB/sec)으로 판독되고 메모리에 저장된다(또 다시 24.8 MB/sec). 그러므로 50 fps로 동작할 때 총 데이터 처리 속도는 248 MB/sec가 된다. 또한 프레임 평균화 모드에서는 메모리로부터 상기 이미지 값을 읽어내는데 이 때 속도는 또 다시 24.8 MB/sec이므로 메모리 밴드 폭은 총 273 MB/sec가 된다. 이미지를 스크린에 표시할 경우에는 메모리로부터 16 비트 픽셀 값을 읽어내고, 픽셀 단위의 룩업 테이블로부터 32 비트 컬러 값을 읽어내고, 마지막으로 디스플레이 메모리에 32 비트 값이 저장되는 바, 이 때 처리 속도는 총 273 MB/sec중 124 MB/sec가 된다. 그리고 방사선 필드는 더 큰 카메라로 이동해간다.
1차 모델을 이용할 경우 하나의 픽셀은 적어도 2개의 32 비트 플로팅 포인트 수/프레임을 필요로 한다. 30초의 데이터 수집 사이클 시간 동안 300 fps의 프레임률로 96,000 픽셀 이미지 프레임을 처리한다면 1회의 데이터 수집 사이클에서 생성되는 데이터는 6.4 GB가 된다. 도14a 내지 도14c는 이것을 도시하는 도면이다. 도 14a는 균일 샘플링 법(Uniform Sampling)이라고 알려진 종래의 에러 샘플링 방법을 도시하는 바, 이 때 조건은 300 fps의 프레임률, 30초의 사이클, 100,000개의 픽셀, 4개의 파라미터 그리고 파라미터 당 4 byte이다. 여기서 신호에 대하여는 3차 다항식, 시간에 대해서는 0차 다항식을 이용한다. 그러나 300 fps의 프레임률, 30초의 사이클, 100,000개의 픽셀, 4개의 파라미터 그리고 파라미터 당 4 byte의 조건 하에서는 13 GB에 이르는 데이터를 수집해서 처리해야하며, 이는 실용성이 없다. 도14b는 본 발명의 비선형 에러 샘플링 방법을 도시하는 바, 이 때 조건은 0차 보간법, 300개의 데이터 집합, 4개의 파라미터이다. 이 방법에 의하면 도14a와 유사한 카메라 작동 조건에서 불과 480 MB 정도의 데이터가 생성되고 이것을 수집해서 처리하면 된다. 이에 의하면 저장 공간과 처리 용량이 종래 기술에 비해 30배는 줄어든다. 여기서 사이클의 초기와 말기는 인위적인 수정을 가하였음을 밝힌다. 도14c는 선형 보간법, 10개의 데이터 집합 및 4개의 파라미터를 이용한 바람직한 비선형 에러 샘플링 방법을 도시한다. 이 방법에 의하면 도14a와 유사한 카메라 작동 조건에서 불과 16 MB 정도의 데이터가 생성되고 이것을 수집해서 처리하면 된다. 이에 의하면 저장 공간과 처리 용량이 도14a의 종래 방법에 비해 30배는 줄어든다. 여기서 사이클의 초기와 말기는 약간의 인위적 수정을 가하였음을 밝힌다. 쌍선형 보정은 신호에 대해서도 선형이고 시간에 대해서도 선형이다. 시간에 대하여 선형 보간을 행하면 신호의 비선형성이 현저히 감소하므로 신호에 대한 선형 보정은 그것으로 충분하다.
도7 및 도14c에 도시되어있는 바와 같이, 최적의 하위 집합 이미지 프레임을 이용하는 선택을 할 수도 있으며 본 발명의 캘리브레이션 방법에서는 그렇게 하고 있다. 데이터 수집 사이클 Ct의 초기에는 시간의 경과에 따른 픽셀 셀의 회로 출력 신호(40)의 변화가 현저하다. 이처럼 신호의 변화가 크기 때문에 캘리브레이션 데이터 집합은 시간의 경과에 따른 출력 신호(40)의 변화가 비교적 완만해지는 수집 사이클 Ct의 말기에서보다 수집 사이클 Ct의 초기 부분에서 상대적으로 더 많은 기준 프레임을 포함하고 있어야한다. 바람직한 실시예에 있어서는 사용자가 노출 시간(즉, 프레임률) 및/또는 검출기 바이어스 전압(5)의 오프 시간을 변경할 수 있도록 자동실행 방법을 이용한다. 물론 설정은 수동으로 행한다. 도시된 그래프에서 막대 하나는 하나의 캘리브레이션 값 집합을 나타내는 것이다.
마스크 대상 픽셀의 선택 방법: 이미징 장치(28)의 픽셀 셀(29) 중 일부는 소재의 결함이나 제조상의 결함 때문에 실제로는 쓸모가 없다. 그러므로 이러한 픽셀 셀(29)은 색출해서 마스킹을 실시해야한다. 다시 말해 그러한 픽셀 셀의 출력 값은 인접한 픽셀 셀(29)로부터 산출한 적정한 값으로 대체해야한다. 본 발명의 캘리브레이션 방법에서는 인접한 픽셀 셀의 출력 신호 집합의 국소적 평균치를 산출하고 이 평균치를 개별 픽셀의 출력 신호 값과 비교한다. 이 결과 본 발명의 캘리브레이션 방법을 비정지형(non-stationary) 방사선 필드에도 적용할 수 있게 된다. 바람직한 실시예에 있어서는 적어도 5회의 완전한 데이터 수집 사이클 동안 기준 방사선 필드 강도를 1회 설정하여 평균 프레임을 계산한다. 이렇게 함으로써 이미징 장치(28) 내에 불량 픽셀 셀(29)이 존재하는 지의 여부를 최단 시간 내에 확실하고 신뢰성 있게 판단할 수 있다.
대체 값의 계산: 모든 불량 픽셀 셀(29)을 찾아낸 후 그들의 출력 신호 값을 국소적 산술 평균치로 대체한다. 다시 말해 고립성 불량 픽셀 셀(29)의 신호 출력을 4개의 양호한 인접 픽셀의 출력 신호 평균치로 대체한다. 이러한 계산 과정에서 불량 픽셀 셀(29)로부터의 픽셀 출력 신호는 제외한다. 4개의 양호한 인접 픽셀 셀(29)을 선택함에 있어서는 가능한 모든 방향에 동일한 가중치가 부여되도록 한다. 예를 들면, 제1 불량 픽셀 셀(29)의 상측의 픽셀 셀(29)도 역시 불량일 경우, 제1 불량 픽셀 셀(29)의 픽셀 출력 신호에 대한 대체 값을 계산함에 있어서는 좌상 측 및 우상 측의 픽셀 셀(29)을 대신 이용한다.
형상 보정 및 비활성 영역 메우기: 이상적인 면에서 볼 때 ASIC 하이브리드 간의 상대 위치는 밀접하고 균일하다. 이는 인접한 하이브리드 간에 비활성 영역(데드 스페이스)이 일부 존재하고 서로 다른 인접 하이브리드 간의 상대적 거리가 달라질 수 있음을 의미한다. 이 문제에 대한 해결책은 2가지 단계로 이루어진다. 먼저, 기준 물체의 캘리브레이션 이미지에 의거하여 하이브리드 간의 거리와 하이브리드의 회전 가능 각도를 측정한다. 이러한 측정은 카메라 그 자체를 측정 장치로 이용하되, 매우 정확한 치수를 갖는 캘리브레이션 완료 기준 물체를 이미징 함으로써 행한다. 다음으로, 거리를 측정한 후 기지 값과 실측치를 비교하여 불일치 여부를 확인한다.
불일치의 보정과 메우기: 하이브리드의 정확한 위치를 알아낸 후 보정 알고리즘을 실행한다. 측정된 거리에 의거하여 이미지 내에서 픽셀의 위치를 정확히 표시하는 그리드를 만든다. 이에 기초해서 쌍선형 보간법(또는 기타 보간법)으로 하 위 픽셀을 새로운 픽셀 값을 갖도록 이동 및 회전시킨다.
이상에서 설명한 본 발명의 실시예는 다양하게 변경 및 변형할 수 있다. 여기에서는 본 발명의 몇 가지 실시예를 예시적으로 도시하고 설명하였으나 이상의 내용을 광범위하게 변형, 변경 및 대체할 수 있음은 물론이다. 경우에 따라서는 본 발명의 특징 중 일부를 다른 특징과 별개로 이용할 수 도 있다. 따라서 이상의 설명은 넓은 의미로 해석되어야하고 예시적인 것으로만 받아들여야하며, 본 발명의 사상 및 범위는 첨부된 특허청구범위에 의해서만 정해지는 것이다.