KR100458776B1 - 이미지 처리 장치, 이미지 처리 방법, 기록 매체 및 촬상장치 - Google Patents

이미지 처리 장치, 이미지 처리 방법, 기록 매체 및 촬상장치 Download PDF

Info

Publication number
KR100458776B1
KR100458776B1 KR10-2001-0076051A KR20010076051A KR100458776B1 KR 100458776 B1 KR100458776 B1 KR 100458776B1 KR 20010076051 A KR20010076051 A KR 20010076051A KR 100458776 B1 KR100458776 B1 KR 100458776B1
Authority
KR
South Korea
Prior art keywords
pixel
interest
pixel value
variance
image
Prior art date
Application number
KR10-2001-0076051A
Other languages
English (en)
Other versions
KR20020044074A (ko
Inventor
오기노데츠오
Original Assignee
지이 메디컬 시스템즈 글로발 테크놀러지 캄파니 엘엘씨
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 지이 메디컬 시스템즈 글로발 테크놀러지 캄파니 엘엘씨 filed Critical 지이 메디컬 시스템즈 글로발 테크놀러지 캄파니 엘엘씨
Publication of KR20020044074A publication Critical patent/KR20020044074A/ko
Application granted granted Critical
Publication of KR100458776B1 publication Critical patent/KR100458776B1/ko

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/20—Image enhancement or restoration by the use of local operators

Abstract

이미지내의 실질적인 구조와 관련된 화소와 다른 화소간의 차이를 증대시키기 위해, 관심 화소로서 정의되는 이미지를 구성하는 각 화소를 가지고, 관심 화소가 속하는 로컬 영역내의 화소값 분산이 결정되고(508), 화소값 분산이 잡음 분산보다 훨씬 큰 경우 관심 화소의 화소값이 유지되며, 그렇지 않은 경우, 관심 화소의 화소값은 억제된다(510, 512).

Description

이미지 처리 장치, 이미지 처리 방법, 기록 매체 및 촬상 장치{IMAGE PROCESSING METHOD AND APPARATUS, RECORDING MEDIUM AND IMAGING APPARATUS}
본 발명은 이미지 처리 방법 및 장치, 기록 매체 및 촬상 장치에 관한 것으로, 보다 구체적으로는 이미지를 구성하는 화소값을 조절하는 이미지 처리 방법 및 장치와, 컴퓨터가 그러한 이미지 처리 기능을 구현하도록 하는 프로그램으로 기록된 매체와, 그러한 이미지 처리 장치를 포함하는 촬상 장치에 관한 것이다.
자기 공명 촬상(magnetic resonance imaging; MRI) 장치에서, 촬상될 대상물이 자석 시스템의 내부 공간, 즉 정적 자계가 생성되는 공간으로 운반되고, 그래디언트 자계 및 고주파 자계가 인가되어 대상물내의 스핀(spin)이 자기 공명 신호를 생성하도록 하며, 수신 신호에 근거하여 단층 촬영 이미지(tomographic image)가 생성된다.
혈액 흐름에서의 스핀과 같이 인체내에서 이동하는 스핀과 조직(tissue)내의 스핀과 같이 이동하지 않는 스핀 간에는, 그래디언트 자계 및 고주파 자계가 스핀에 미치는 영향이 상이하다. 이러한 차이점을 이용하여, 인체내에서 이동하는 스핀의 이미지, 예를 들면 혈액 흐름 이미지를 캡처할 수 있다.
혈액 흐름 이미지 캡처시, TOF(time-of-flight) 기법, 위상 대조(phase contrast; PC) 기법 등을 이용한다.
이들 기법 중 하나를 이용하여 3차원 영역에서의 혈액 흐름 투영 이미지를 획득함으로써, 3차원 영역에 대한 멀티슬라이스(multi-slice) 혈액 흐름 단층 촬영 이미지를 캡처하고, 슬라이스 두께 방향에서 멀티슬라이스 혈액 흐름 단층 촬영 이미지에 대해 최대 강도 투영(maximum intensity projection; MIP)을 수행한다.
전술한 바와 같이 혈액 흐름의 투영 이미지가 획득될 때, 희미한 혈액 흐름은 잡음에 의해 불명료해지기 때문에 투영되지 않을 수도 있다. 더욱이, 슬라이스들 사이에 이미지의 평균 신호 강도가 상이할 때, 작은 평균 신호 강도를 갖는 이미지내의 혈액 흐름 이미지는 큰 평균 신호 강도를 갖는 이미지내의 잡음에 의해 불명료해지기 때문에 투영될 수 없다.
따라서, 본 발명의 목적은 이미지내의 실질적인 구조와 관련된 화소와 다른 화소간의 차이를 증대시키는 이미지 처리 방법 및 장치와, 컴퓨터가 그러한 이미지 처리 기능을 구현하도록 하는 프로그램으로 기록된 매체와, 그러한 이미지 처리 장치를 포함하는 촬상 장치에 관한 것이다.
(1) 전술한 문제를 해결하기 위한 하나의 양상에서, 본 발명은 관심 화소(a pixel of interest)로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 것을 특징으로 하는 이미지 처리 방법이다.
(2) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 수단을 포함하는 것을 특징으로 하는 이미지 처리 장치이다.
(3) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 것을 특징으로 하는 기록 매체이다.
(4) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 수단을 포함하는 것을 특징으로 하는, 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치이다.
상기 (1)-(4)에 기술된 바와 같은 양상에서의 본 발명에 따르면, 관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 결정된 분산이 잡음 분산보다 훨씬 큰 경우 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 관심 화소의 화소값을 억제하므로, 이미지내의 실질적인 구조와 관련된 화소와 다른 화소간의 차이를 증대시킬 수 있다.
(5) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 것을 특징으로 하는 이미지 처리 방법이다.
(6) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 수단을 포함하는 것을 특징으로 하는 이미지 처리 장치이다.
(7) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 것을 특징으로 하는 기록 매체이다.
(8) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 수단을 포함하는 것을 특징으로 하는, 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치이다.
상기 (5)-(8)에 기술된 바와 같은 양상에서의 본 발명에 따르면, 관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 결정된 분산이 잡음 분산보다 훨씬 큰 경우 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 관심 화소의 화소값을 유지하므로, 이미지내의 실질적인 구조와 관련된 화소와 다른 화소간의 차이를 증대시킬 수 있다.
더욱이, 상기 (1)-(8)에서 기술된 바와 같은 양상에서, 이미지는 혈액 흐름 이미지일 수 있으므로, 혈액 흐름 이미지와 관련된 화소와 다른 화소간의 차이가 증대된다.
또한, 상기 (4) 및 (8)에 기술된 바와 같은 양상에서의 본 발명에 따르면, 신호는 자기 공명 신호일 수 있으므로, 자기 공명 이미지에 대한 발명이 구현된다.
(9) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 것을 포함하는 화소값 조절을 수행하고, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 것을 특징으로 하는 이미지 처리 방법이다.
(10) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 수단과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 수단을 포함하는 것을 특징으로 하는 이미지 처리 장치이다.
(11) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 것을 포함하는 화소값 조절을 수행하는 기능과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 것을 특징으로 하는 기록 매체이다.
(12) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 수단과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 수단을 포함하는 것을 특징으로 하는, 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치이다.
상기 (9)-(12)에 기술된 바와 같은 양상에서의 본 발명에 따르면, 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 결정된 분산이 잡음 분산보다 훨씬 큰 경우 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 관심 화소의 화소값을 억제하므로, 이미지내의 실질적인 구조와 관련된 화소와 다른 화소간의 차이를 증대시킬 수 있다. 또한, 그와 같이 차이가 증대된 멀티슬라이스 이미지는 최대 강도 투영(maximum-intensity-projected)되므로, 희미한 신호 강도를 가지며 이미지내의 실질적인 구조와 관련된 화소의 투영 이미지를 획득할 수 있다.
(13) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 것을 포함하는 화소값 조절을 수행하고, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 것을 특징으로 하는 이미지 처리 방법이다.
(14) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 수단과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 수단을 포함하는 것을 특징으로 하는 이미지 처리 장치이다.
(15) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 것을 포함하는 화소값 조절을 수행하는 기능과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 것을 특징으로 하는 기록 매체이다.
(16) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 수단과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 수단을 포함하는 것을 특징으로 하는, 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치이다.
상기 (13)-(16)에 기술된 바와 같은 양상에서의 본 발명에 따르면, 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 결정된 분산이 잡음 분산보다 훨씬 큰 경우 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 관심 화소의 화소값을 유지하므로, 이미지내의 실질적인 구조와 관련된 화소와 다른 화소간의 차이를 증대시킬 수 있다. 또한, 그와 같이 차이가 증대된 멀티슬라이스 이미지는 최대 강도 투영되므로, 희미한 신호 강도를 가지며 이미지내의 실질적인 구조와 관련된 화소의 투영 이미지를 획득할 수 있다.
더욱이, 상기 (9)-(16)에서 기술된 바와 같은 양상에서의 본 발명에 따르면, 이미지는 혈액 흐름 이미지일 수 있으므로, 혈액 흐름 이미지와 관련된 화소와 다른 화소간의 차이가 증대되고, 그와 같이 차이가 증대된 멀티슬라이스 이미지를 최대 강도 투영함으로써 희미한 혈액 흐름의 투영 이미지를 획득할 수 있다.
또한, 상기 (12) 및 (16)에 기술된 바와 같은 양상에서의 본 발명에 따르면, 신호는 자기 공명 신호일 수 있으므로, 자기 공명 이미지에 대한 발명이 구현된다.
(17) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하고, 상기 가산된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 것을 포함하는 화소값 조절을 수행하고, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 것을 특징으로 하는 이미지 처리 방법이다.
(18) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값분산을 가산하는 가산 수단과, 상기 가산된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 수단과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 수단을 포함하는 것을 특징으로 하는 이미지 처리 장치이다.
(19) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과, 상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 기능과, 상기 가산된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 것을 포함하는 화소값 조절을 수행하는 기능과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 것을 특징으로 하는 기록 매체이다.
(20) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 가산 수단과, 상기 가산된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 수단과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 수단을 포함하는 것을 특징으로 하는, 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치이다.
상기 (17)-(20)에 기술된 바와 같은 양상에서의 본 발명에 따르면, 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하고, 가산된 분산이 잡음 분산보다 훨씬 큰 경우 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 관심 화소의 화소값을 억제하므로, 이미지내의 실질적인 구조와 관련된 화소와 다른 화소간의 차이를 더 증대시킬 수 있다. 또한, 그와 같이 차이가 증대된 멀티슬라이스 이미지는 최대 강도 투영되므로, 희미한 신호 강도를 가지며 이미지내의 실질적인 구조와 관련된 화소의 보다 우수한 투영 이미지를 획득할 수 있다.
(21) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하고, 상기 가산된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 것을 포함하는 화소값 조절을 수행하고, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 것을 특징으로 하는 이미지 처리 방법이다.
(22) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 가산 수단과, 상기 가산된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 수단과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 수단을 포함하는 것을 특징으로 하는 이미지 처리 장치이다.
(23) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과, 상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 기능과, 상기 가산된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 것을 포함하는 화소값 조절을 수행하는 기능과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 것을 특징으로 하는 기록 매체이다.
(24) 전술한 문제를 해결하기 위한 다른 양상에서, 본 발명은 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 수단과, 상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 가산 수단과, 상기 가산된 분산이 잡음 분산보다 훨씬 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 수단과, 상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 수단을 포함하는 것을 특징으로 하는, 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치이다.
상기 (21)-(24)에 기술된 바와 같은 양상에서의 본 발명에 따르면, 관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하고, 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하고, 가산된 분산이 잡음 분산보다 훨씬 큰 경우 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 관심 화소의 화소값을 유지하므로, 이미지내의 실질적인 구조와 관련된 화소와 다른 화소간의 차이를 더 증대시킬 수 있다. 또한, 그와 같이 차이가 증대된 멀티슬라이스 이미지는 최대 강도 투영되므로, 희미한 신호 강도를 가지며 이미지내의 실질적인 구조와 관련된 화소의 보다 우수한 투영 이미지를 획득할 수 있다.
더욱이, 상기 (17)-(24)에서 기술된 바와 같은 양상에서의 본 발명에 따르면, 이미지는 혈액 흐름 이미지일 수 있으므로, 혈액 흐름 이미지와 관련된 화소와 다른 화소간의 차이가 더 증대되고, 그와 같이 차이가 증대된 멀티슬라이스 이미지를 최대 강도 투영함으로써 희미한 혈액 흐름의 보다 우수한 투영 이미지를 획득할 수 있다.
또한, 상기 (20) 및 (24)에 기술된 바와 같은 양상에서의 본 발명에 따르면, 신호는 자기 공명 신호일 수 있으므로, 자기 공명 이미지에 대한 발명이 구현된다.
전술한 (1)-(24)에 기술된 바와 같은 양상에서의 본 발명에서, 전체 이미지에 걸쳐 정의된 다수의 로컬 영역 각각에 대한 화소값의 오차 제곱의 합(residual sums of squares)을 결정하고, 오차 제곱의 합의 히스토그램을 결정하고, 그 후 잡음 분산은 이미지에 근거하여 직접 획득할 수 있다는 점에서, 히스토그램의 피크를 제공하는 오차 제곱의 합에 근거하여 잡음 분산을 결정하는 것이 바람직하다.
따라서, 본 발명은 이미지내의 실질적인 구조와 관련된 화소와 다른 화소간의 차이를 증대시키는 이미지 처리 방법 및 장치와, 컴퓨터가 그러한 이미지 처리 기능을 수행하도록 하는 프로그램으로 기록된 매체와, 그러한 이미지 처리 장치를 포함하는 촬상 장치를 제공할 수 있다.
본 발명의 다른 목적 및 이점은 첨부 도면에 기술된 바와 같은 본 발명의 바람직한 실시예의 설명으로부터 명확해질 것이다.
도 1은 본 발명의 일실시예에 따른 장치의 블록도,
도 2는 본 발명의 일실시예에 따른 장치의 블록도,
도 3은 도 1 또는 2에 도시된 장치의 동작 흐름도,
도 4는 멀티슬라이스 이미지의 개념도,
도 5는 도 3에 도시된 흐름도의 일부에 대한 상세 흐름도,
도 6은 도 5에 도시된 흐름도의 일부에 대한 상세 흐름도,
도 7은 개념적인 히스토그램도,
도 8은 개념적인 히스토그램도,
도 9는 관심 화소와 로컬 영역간의 관계를 나타내는 도면,
도 10은 도 5에 도시된 흐름도의 일부로서 삽입된 절차의 흐름도,
도 11은 화소값 조절의 영향을 나타내는 이미지 프로파일의 예를 도시한 도면,
도 12는 화소값 조절의 영향을 나타내는 이미지 프로파일의 예를 도시한 도면,
도 13은 도 3에 도시된 흐름도의 일부에 대한 상세 흐름도,
도 14는 화소값 조절의 영향을 나타내는 이미지 프로파일의 예를 도시한 도면,
도 15는 화소값 조절의 영향을 나타내는 이미지 프로파일의 예를 도시한 도면,
도 16은 최대 강도 투영의 개념도.
도면의 주요 부분에 대한 부호의 설명
100 : 자석 시스템 130 : 그래디언트 구동 섹션
140 : RF 구동 섹션 150 : 데이터 수집 섹션
160 : 제어 섹션 170 : 데이터 처리 섹션
180 : 디스플레이 섹션 190 : 조작 섹션
이하, 첨부 도면을 참조하여 본 발명의 실시예를 기술할 것이다. 본 발명은 이들 실시예에 한정되지 않음을 알아야 한다. 도 1은 본 발명의 실시예인, 촬상 장치 또는 자기 공명 촬상(MRI) 장치의 블록도를 도시한다. 본 장치의 구성은 본 발명에 따른 장치의 실시예를 나타낸다. 본 장치의 동작은 본 발명에 따른 방법의 실시예를 나타낸다.
도 1에 도시된 바와 같이, 본 발명의 장치는 자석 시스템(100)을 갖는다. 자석 시스템(100)은 주 자계 코일 섹션(102), 그래디언트 코일 섹션(106) 및 RF(radio frequency) 코일 섹션(108)을 갖는다. 이들 코일 섹션은 전반적으로 원통 형상이며, 동심으로(concentrically) 배치된다. 촬상될 대상물(300)이 장치대(cradle)(500)상에 위치되며, 도시되지 않은 운반 수단에 의해 자석 시스템(100)의 전반적으로 원통형인 내부 공간(보어(bore))으로, 그리고 내부 공간으로부터 운반된다.
주 자계 코일 섹션(102)은 자석 시스템(100)의 내부 공간에 정적 자계를 생성한다. 정적 자계의 방향은 대상물(300)의 몸체 축 방향과 전반적으로 평행하다. 즉, "수평(horizontal)" 자계가 생성된다. 주 자계 코일 섹션(102)은 예를 들면, 초전도성 코일(superconductive coil)을 이용하여 만들어진다. 주 자계 코일 섹션(102)은 초전도성 코일에 한정되지 않으며, 보통의 전도성 코일 등을 이용하여 만들어질 수 있음을 쉽게 알 것이다.
그래디언트 코일 섹션(106)은 그래디언트를 정적 자계 강도에 전달하기 위해그래디언트 자계를 생성한다. 생성될 그래디언트 자계는 3 가지 자계, 즉 슬라이스(slice) 그래디언트 자계, 판독(readout) 그래디언트 자계 및 위상 인코딩(phase encoding) 그래디언트 자계이다. 그래디언트 코일 섹션(106)은 이들 3 개의 그래디언트 자계에 대응하는 3 개의 그래디언트 코일(도시되지 않음)을 갖는다.
RF 코일 섹션(108)은 정적 자계 공간내 대상물(300) 내부의 스핀을 여기하기 위해 고주파 자계를 생성한다. 이하, 고주파 자계의 생성은, 때로는 RF 여기 신호의 전송이라고 지칭될 것이다. 또한, RF 코일 섹션(108)은 여기된 스핀에 의해 생성된 전자기파(electromagnetic wave), 즉 자기 공명 신호를 수신한다.
RF 코일 섹션(108)은 도시되지 않은, 전송 및 수신 코일을 갖는다. 전송 및 수신 코일에 대해, 동일한 코일 또는 개별적인 전용 코일이 사용될 수 있다.
그래디언트 코일 섹션(106)은 그래디언트 구동 섹션(130)과 접속된다. 그래디언트 구동 섹션(130)은 그래디언트 코일 섹션(106)에 구동 신호를 공급하여 그래디언트 자계를 생성한다. 그래디언트 구동 섹션(130)은 그래디언트 코일 섹션(106)내의 3 개의 그래디언트 코일에 대응하는 3 개의 구동 회로(도시되지 않음)를 갖는다.
RF 코일 섹션(108)은 RF 구동 섹션(140)과 접속된다. 그래디언트 구동 섹션(140)은 RF 코일 섹션(108)에 구동 신호를 공급하여 RF 여기 신호를 전송하고, 그로 인해 대상물(300)내의 스핀이 여기된다.
RF 코일 섹션(108)은 데이터 수집 섹션(150)과 접속된다. 데이터 수집 섹션(150)은 RF 코일 섹션(108)에 의해 수신된 수신 신호를 모으고, 그 신호를 뷰 데이터(view data)로서 수집한다.
그래디언트 구동 섹션(130), RF 구동 섹션(140) 및 데이터 수집 섹션(150)은 제어 섹션(160)과 접속된다. 제어 섹션(160)은 그래디언트 구동 섹션(130), RF 구동 섹션(140) 및 데이터 수집 섹션(150)을 제어하여 촬상을 수행한다.
데이터 수집 섹션(150)의 출력은 데이터 처리 섹션(170)에 접속된다. 데이터 처리 섹션(170)은, 예를 들면 컴퓨터를 사용하여 구성된다. 데이터 처리 섹션(170)은 도시되지 않은 메모리를 갖는다. 메모리는 데이터 처리 섹션(170)을 위한 프로그램 및 몇 가지 종류의 데이터를 저장한다. 본 발명의 장치의 기능은 메모리내에 저장된 프로그램을 실행하는 데이터 처리 섹션(170)에 의해 구현된다.
데이터 처리 섹션(170)은 데이터 수집 섹션(150)으로부터 모아진 뷰 데이터를 메모리내에 저장한다. 메모리내에 데이터 공간이 형성된다. 데이터 공간은 2 차원 푸리에(Fourier) 공간을 구성한다. 2 차원 푸리에 공간은, 때로는 k 공간으로 지칭된다. 데이터 처리 섹션(170)은 2 차원 푸리에 공간에서의 데이터에 대해 2 차원 역 푸리에 변환을 수행하여 대상물(300)의 이미지를 생성(재구성)한다.
2 차원 역 푸리에 변환에 의해 재구성된 이미지는 복소수(complex number)의 화소값을 갖는다. 복소수의 절대값을 이용하여 절대값 이미지를 구성한다. 복소수의 실수부(real part)를 이용하여 실수부 이미지를 구성할 수 있다. 복소수의 허수부(imaginary part)를 이용하여 허수부 이미지를 구성할 수 있다. 실수 및 허수부 모두는 양(positive) 및 음(negative)의 값을 가질 수 있다. 그러한 이미지는, 때로는 양-음(positive-negative) 이미지라고 지칭된다.
데이터 처리 섹션(170)은 이미지 처리를 수행하여 재구성 이미지에 대한 화소값의 분산을 결정하는 기능을 갖는다. 또한, 데이터 처리 섹션(170)은 이미지 처리를 수행하여 재구성 이미지에 대한 잡음 분산을 결정하는 기능을 갖는다. 데이터 처리 섹션(170)은 이미지 처리를 수행하여 재구성 이미지에 대한 화소값을 조절하는 기능을 더 갖는다. 또한, 데이터 처리 섹션(170)은 이미지 처리를 수행하여 화소값 조절이 수행된 이미지에 대해 최대 강도 투영(MIP)을 실행하는 기능을 갖는다. 데이터 처리 섹션(170)의 그러한 이미지 처리 기능에 대해서는 이하에 기술될 것이다.
데이터 처리 섹션(170)은 본 발명의 이미지 처리 장치의 실시예이다. 본 장치의 구성은 본 발명에 따른 장치의 실시예를 나타낸다. 본 장치의 동작은 본 발명에 따른 방법의 실시예를 나타낸다.
데이터 처리 섹션(170)은 제어 섹션(160)에 접속된다. 데이터 처리 섹션(170)은 제어 섹션(160)보다 상위에서 제어 섹션(160)을 제어한다. 데이터 처리 섹션(170)은 디스플레이 섹션(180) 및 조작 섹션(190)과 접속된다. 디스플레이 섹션(180)은 그래픽 디스플레이 등을 포함한다. 조작 섹션(190)은 포인팅 장치(pointing device)를 구비한 키보드 등을 포함한다.
디스플레이 섹션(180)은 데이터 처리 섹션(170)으로부터 출력된 재구성 이미지 및 몇 가지 종류의 정보를 디스플레이한다. 조작 섹션(190)은 조작자에 의해 조작되며, 조작 섹션(190)은 몇 가지의 코맨드, 정보 등을 데이터 처리 섹션(170)에 입력한다. 조작자는 디스플레이 섹션(180) 및 조작 섹션(190)을 통해 본 장치를 상호작용적으로 동작시킨다.
도 2는 본 발명의 일실시예인, 다른 타입의 MRI 장치의 블록도이다. 본 장치의 구성은 본 발명에 따른 장치의 실시예를 나타낸다.
도 2에 도시된 장치는 도 1에 도시된 장치에서의 자석 시스템과는 다른 타입의 자석 시스템(100')을 갖는다. 본 장치는 자석 시스템(100')을 제외하고는, 도 1에 도시된 장치와 유사한 구성을 가지므로, 유사한 부분은 유사한 참조 번호로 표시하며, 그 설명은 생략될 것이다.
자석 시스템(100')은 주 자계 자석 섹션(102'), 그래디언트 코일 섹션(106') 및 RF 코일 섹션(108')을 갖는다. 주 자계 자석 섹션(102') 및 코일 섹션 각각은 공간을 가로질러 서로 마주보는 한 쌍의 부재(member)로 구성된다. 이들 섹션은 전반적으로 디스크 형상이며, 공통의 중심 축을 갖도록 배치된다. 대상물(300)은 장치대(500)상에 위치되고, 도시되지 않은 운반 수단에 의해 자석 시스템(100')의 내부 공간(보어)으로, 그리고 내부 공간으로부터 운반된다.
주 자계 자석 섹션(102')은 자석 시스템(100')의 내부 공간에 정적 자계를 생성한다. 정적 자계의 방향은 대상물(300)의 몸체 축 방향과 일반적으로 직교한다. 즉, "수직(vertical)" 자계가 생성된다. 주 자계 자석 섹션(102')은 예를 들면, 영구 자석을 이용하여 만들어진다. 주 자계 자석 섹션(102')은 영구 자석에 한정되지 않으며, 초전도성 또는 보통 전도성의 전자석 등을 이용하여 만들어질 수 있음을 쉽게 알 것이다.
그래디언트 코일 섹션(106')은 그래디언트를 정적 자계 강도에 전달하기 위해 그래디언트 자계를 생성한다. 생성될 그래디언트 자계는 3 가지 자계, 즉 슬라이스 그래디언트 자계, 판독 그래디언트 자계 및 위상 인코딩 그래디언트 자계이다. 그래디언트 코일 섹션(106')은 이들 3 개의 그래디언트 자계에 대응하는 3 개의 그래디언트 코일(도시되지 않음)을 갖는다.
RF 코일 섹션(108')은 정적 자계 공간내 대상물(300) 내부의 스핀을 여기하기 위해 RF 여기 신호를 전송한다. 또한, RF 코일 섹션(108')은 여기된 스핀에 의해 생성된 자기 공명 신호를 수신한다. RF 코일 섹션(108')은 도시되지 않은, 전송 및 수신 코일을 갖는다. 전송 및 수신 코일에 대해, 동일한 코일 또는 개별적인 전용 코일이 사용될 수 있다.
도 3은 본 발명의 장치의 동작 흐름도를 도시한다. 도 1 및 2에 도시된 두 장치는 동일한 방법으로 동작한다. 도 3에 도시된 바와 같이, 단계(302)에서 혈액 흐름 촬상이 수행된다. 혈액 흐름 촬상을 위해, TOF 기법, 위상 대조(PC) 기법 등을 이용한다. 더욱이, 멀티슬라이스에서 촬상이 수행된다. 따라서, 도 4에 개념적으로 도시된 바와 같이, 대상물(300)의 3 차원 영역에 대해 멀티슬라이스 혈액 흐름 단층 촬영 이미지 S1, S2, S3,..., Sm이 캡처된다.
다음, 단계(304)에서, 혈액 흐름 단층 촬영 이미지 S1, S2, S3,..., Sm에 대해 화소값 조절이 수행된다. 화소값 조절은 데이터 처리 섹션(170)의 데이터 처리 기능에 의해 구현된다. 이하, 혈액 흐름 단층 촬영 이미지를 간단히 이미지라고 지칭할 것이다.
도 5는 화소값 조절의 상세 흐름도를 도시한다. 도시된 바와 같이,단계(502)에서 슬라이스 선택이 수행된다. 따라서, 이미지 S1, S2, S3,..., Sm 중 하나, 예를 들면, 이미지 S1이 선택된다.
다음, 단계(504)에서 잡음 분산 계산이 수행된다. 단계(504)에서 잡음 분산을 계산하는 데이터 처리 섹션(170)은 본 발명의 잡음 분산 계산 수단의 실시예이다. 도 6은 잡음 분산 계산의 상세 흐름도를 도시한다. 도시된 바와 같이, 단계(602)에서 이미지내의 로컬 영역이 정의된다. 로컬 영역은 다음 단계에서의 계산을 위해 이용되는 화소값이 속하는 영역이다. 이미지의 중심에서의 로컬 영역은, 예를 들면, 제 1 영역으로 정의된다.
로컬 영역으로서, NxN 화소 매트릭스가 이용된다. 예를 들면, N은 9이다. 매트릭스 크기는 이러한 값에 한정되지 않으며, 임의의 적절한 크기일 수 있음을 알아야 한다. 더욱이, 화소 매트릭스는 정사각형의 매트릭스로 한정되지 않으며, 소정의 화소에 중심을 둔 임의의 적절한 영역일 수 있다. 이하, 로컬 영역을, 때로는 단순하게 영역으로 지칭할 것이다.
다음, 단계(604)에서 그 영역에 속하는 화소값의 오차 제곱의 합 S가 결정된다. 특히, 다음의 수학식과 같으며,
여기서, Pi는 화소값,는 Pi에 중심을 둔 NxN 영역내 화소값의 평균이다. 더욱이, k는 예를 들면, 81이다.
다음, 단계(606)에서, 모든 로컬 영역에 대해 상기 처리가 완료되었는지의 여부에 대한 판정이 수행되고, 완료되지 않았다면, 단계(608)에서 로컬 영역은 시프트된다. 따라서, 예를 들면 인접한 NxN 영역이 새로운 로컬 영역으로서 선택된다.
단계(604)의 처리가 새로운 로컬 영역에 대해 수행되어, 화소값의 오차 제곱의 합이 결정된다. 그 후, 마찬가지 방법으로 이미지내의 모든 로컬 영역에 대해 화소값의 오차 제곱의 합이 결정된다.
그에 따라 획득된 오차 제곱의 합은 x2분포를 가지며, 평균값은 k·σ2이다. k가 클 때, x2분포는 가우시안 분포(Gaussian distribution)에 근사화되며, 그 피크 위치는 대략 k·σ2에 놓이게 된다.
다음, 단계(610)에서, 오차 제곱의 합 S의 히스토그램이 생성된다.
도 7은 이미지가 절대값 이미지인 경우 오차 제곱의 합 S의 히스토그램의 개념을 도시한다. 도시된 바와 같이, 히스토그램은 3 개의 분포 곡선 a, b, c로 이루어진다.
분포 곡선 a는 균일한 구조 부분내 잡음으로부터의 가우시안 분포 곡선이다. 분포 곡선 b는 대상물(300)을 포함하지 않는 FOV(field of view)의 일부내의 잡음, 즉 배경 잡음으로부터의 레일리(Rayleigh) 분포 곡선이다. 이미지는 절대값 이미지이기 때문에, 배경 잡음으로부터의 분포 곡선은 가우시안 분포가 아닌 레일리 분포를 따른다. 분포 곡선 c는 대상물의 미세한(fine) 구조로부터의 곡선이며, 다른 두 곡선과는 달리 중간 분포를 나타낸다.
단계(612)에서, 그러한 히스토그램에 대해 피크 위치 검출이 수행된다. 따라서, 가우시안 분포 곡선 a에 대해 피크 위치 s1이 검출되고, 레일리 분포 곡선 b에 대해 피크 위치 s2가 검출된다.
히스토그램은 실제로는 이산적인 값을 가지므로, 피크 위치가 양호한 정확도로 검출될 수 있다는 점에서, 바람직하게는, 피크 검출 이전에 단계(612)에서 함수에 대한 일치(fitting to a function)가 수행된다. 조화에 이용된 함수는, 예를 들면 각각 가우시안 분포 함수 및 레일리 분포 함수이다. 그러나, 함수는 이것에 한정되지 않으며 임의의 다른 적절한 함수일 수 있다.
다음, 단계(614)에서, 잡음 분산이 계산된다. 잡음 분산 계산은 피크 위치 s1 또는 s2에 근거하여 수행된다.
s1, s2 및 σ는 다음의 수학식과 같은 각각의 관계를 가지므로, 이들 관계로부터 σ의 값이 결정된다.
σ의 값은 그것이 수학식 2로부터 결정되든지, 혹은 수학식 3으로부터 결정되든지 간에 동일하다. 결정된 σ의 값은 잡음 분산 Vn으로서 메모리내에 저장된다.
분포 곡선 c의 일부 조건하에서는, 가우시안 분포 곡선의 피크 위치 s1이 정확하게 검출되지 않을 수도 있다. 이러한 경우, σ의 값은 레일리 분포 곡선 b의 피크 위치 s2에 근거하여 결정된다. 더욱이, 보다 큰 비율(proportion)의 배경 부분 영역을 갖는 이미지에 대해, 양호한 정확도로 잡음 분산을 결정하는 데에는 레일리 분포 곡선 b가 더 적절하다.
전술한 내용은 절대값 이미지의 경우에 대한 것이지만, 처리될 이미지가 양-음 이미지, 즉 실수부 이미지 또는 허수부 이미지인 경우, 배경 부분에서의 잡음은 0에 중심을 둔 양의 값 및 음의 값을 갖는다.
따라서, 단계(610)에서 생성된 히스토그램은 도 8에 예시된 바와 같이 되며, 더 이상은 레일리 분포를 갖지 않는다. 이 경우, 단계(614)에서 가우시안 분포 곡선 a의 피크 위치 s1에 근거하여 잡음 분산이 결정된다.
실제 캡처되는 이미지에 근거하여 잡음 분산값을 직접 획득할 수 있다. 잡음 분산이 이전에 알려진다면, 그 분산을 이용할 수 있어, 계산을 생략할 수 있다.
잡음 분산 Vn이 결정된 후, 도 5의 흐름도의 단계(506)에서, 관심 화소가 이미지내에 정의된다. 제 1 관심 화소는, 예를 들면 이미지 중심에 있는 화소이다.
다음, 단계(508)에서, 관심 화소를 포함하는 로컬 영역내의 화소값 분산 Vi가 계산된다. 관심 화소를 포함하는 로컬 영역은, 예를 들면, 도 9에 도시된 바와 같이 관심 화소 i에 중심이 있는 5x5 매트릭스이다. 매트릭스 크기는 이러한 값에 한정되지 않으며, 임의의 적절한 값일 수 있음을 알아야 한다. 더욱이, 화소 매트릭스는 사각형 매트릭스에 한정되지 않으며, 소정의 화소에 중심이 있는 임의의 적절한 영역일 수 있다. 이하, 로컬 영역은, 때로는 간단하게 영역으로 지칭할 것이다. 단계(508)에서 화소값 분산 Vi를 계산하는 데이터 처리 섹션(170)은 본 발명의 분산 계산 수단의 실시예이다.
다음의 수학식을 이용하여 화소값 분산 Vi를 계산하며, 여기서 k=25이다.
다음, 단계(510)에서 화소값 분산 Vi가 잡음 분산 Vn보다 훨씬 큰지의 여부가 판정된다. 그러한 판정은 다음과 같은 수학식을 이용하여 수행되며, 여기서는 임계값이다.
임계값로는 1보다 큰 적절한 값이 이용된다.
관심 화소를 포함하는 로컬 영역내의 화소값 분산이 잡음 분산보다 훨씬 크지 않으면, 로컬 영역내의 이미지는 아마도 현저한 구조를 갖지 않을 것이며, 화소값 분산은 아마도 잡음으로부터 발생될 것이다.
그러므로, 이러한 경우 단계(512)에서 관심 화소의 화소값이 억제된다. 화소값의 억제는, 예를 들면, 화소값에 계수 α를 승산함으로써 달성된다. 계수 α의 값은 1보다 작은 양수, 예를 들면, 0.8이다. 따라서, 관심 화소의 화소값은,예를 들면 원래 값의 0.8 배로 감소된다. 그러나, 계수 α의 값은 0.8에 한정되지 않으며, 임의의 적절한 값일 수 있다. 더욱이, 화소값의 억제는, 예를 들면 화소값으로부터 사전결정된 소정의 값을 감산함으로써 달성될 수 있다. 상수값은 화소의 최소값을 초과하지 않음을 알아야 한다.
관심 화소를 포함하는 로컬 영역내의 화소값 분산 Vi가 잡음 분산 Vn보다 훨씬 크다면, 로컬 영역내의 이미지는 아마도 에지(edge)와 같은 특정의 구조를 가질 것이며, 화소값 분산은 아마도 이미지의 구조로부터 발생될 것이다. 이러한 경우, 화소값에 대해 특별한 작용이 적용되지 않는다. 따라서, 관심 화소의 화소값은 그 원래 값을 유지한다. 그러한 화소값 조절을 수행하는 데이터 처리 섹션(170)은 본 발명의 화소값 조절 수단의 실시예이다.
다음, 단계(514)에서, 모든 관심 화소에 대해 상기 처리가 완료되었는지의 여부에 대한 판정이 수행되고, 완료되지 않았다면, 단계(516)에서 관심 화소는 예를 들면 인접한 것에 대해 시프트되고, 단계(508)로부터의 처리가 수행된다. 그 후, 이미지 S1내의 모든 화소에 대한 화소값을 조절하기 위해 동일한 처리가 반복된다.
그 후, 단계(518)에서, 모든 슬라이스에 대해 상기 처리가 완료되었는지의 여부에 대한 판정이 수행되고, 완료되지 않았다면, 단계(520)에서 슬라이스는 시프트되고, 해당 슬라이스의 이미지에 대해 동일한 처리가 수행된다. 그 후, 모든 이미지 S1-Sm내의 화소에 대해 화소값 조절을 수행하도록 동일한 처리가 반복된다.
단계들(508 및 510) 사이에, 도 10의 흐름도에 도시된 바와 같은 단계가 부가될 수 있다. 특히, 단계(702)에서 인접 슬라이스내의 대응하는 관심 화소를 포함하는 로컬 영역에 대해 화소값 분산 Vi'이 계산된다.
'인접 슬라이스'라는 용어는 화소값 분산 Vi가 단계(508)에서 결정된 슬라이스에 인접하는 하나 이상의 슬라이스를 의미한다. 그러한 슬라이스를 위해, 예를 들면 전후로 인접하는 슬라이스, 또는 전후로 인접하는 슬라이스들을 이용할 수 있다.
단계(704)에서, 화소값의 분산(들) Vi'이 Vi에 가산되고, 가산된 값은 다음 단계(510)에서의 판정에 이용하기 위한 화소값 분산 Vi로서 정의된다. 가산시에 적절한 가중치가 Vi'에 적용될 수도 있다. 단계(702)에서 화소값 분산을 계산하는 데이터 처리 섹션(170)은 본 발명의 화소값 분산 계산의 실시예이다. 단계(704)에서 화소값 분산을 가산하는 데이터 처리 섹션(170)은 본 발명의 가산 수단의 실시예이다.
따라서, 다수의 슬라이스를 통한 구조가 상기 처리에 의해 획득된 화소값 분산 Vi에 반영된다. 따라서, 예를 들어 이미지가 하나의 이미지상에서 하나의 포인트로서 나타나야 하는, 슬라이스를 통과하는 방향에서 혈액 흐름 이미지가 존재한다면, 그러한 구조를 정확하게 반영하는 화소값 분산을 획득할 수 있고, 분산에 근거하여 보다 정확한 화소값 조절을 달성할 수 있다.
도 11은 그러한 화소값 조절의 영향을 화소값 프로파일에서의 변화로서 도시한다. 도 11에서의 심볼 B는 화소값 조절 이전의 프로파일을 나타내고, 배경 잡음을 통해 뚜렷한 혈액 흐름 이미지 b1 및 희미한 혈액 흐름 이미지 b2가 존재한다.
전술한 화소값 조절의 결과로서, 그러한 프로파일은 혈액 흐름 이미지 b1 및 b2의 화소값을 유지하면서, 예를 들면 0.8 배로 억제된 배경 잡음의 화소값을 가지며, 그 결과 도 11의 A에서 도시된 바와 같은 프로파일이 된다. 프로파일 A에서, 혈액 흐름 이미지 b1은 말할 것도 없이, 원래 이미지에서 희미하였던 혈액 흐름 이미지 b2는 배경 잡음으로부터 증대된 차이를 나타내며 뚜렷해진다. 따라서, 원래 이미지에서 희미하였던 혈액 흐름 이미지 b2의 유도성(elicitability)이 향상될 수 있다.
도 12는 화소값 조절의 다른 영향을 도시한다. 도 12에서의 심볼 P 및 Q는 상이한 슬라이스의 두 이미지 프로파일을 나타내며, 프로파일 Q의 배경 잡음 레벨은 프로파일 P에서의 뚜렷한 혈액 흐름 이미지 b1의 신호 강도보다 크다.
전술한 화소값 조절에 의해, 그러한 프로파일은 예를 들면 0.8 배로 억제된 배경 잡음의 화소값을 가지며, 따라서, 도 12의 Q'에 도시된 바와 같이 혈액 흐름 이미지 b1 및 b2의 신호 강도에 대해 감소된 잡음 레벨을 갖는 프로파일을 획득할 수 있다. 따라서, 슬라이스 q의 이미지 잡음 레벨로부터 혈액 흐름 이미지 b1 및 b2의 차이 또한 뚜렷해지며, 두 이미지를 유도할 수 있다.
도 13은 화소값 조절의 다른 절차에 대한 흐름도를 도시한다. 도 13에서, 도 5에 도시된 단계와 유사한 단계는 유사한 참조 번호에 의해 지정되며, 그 설명은 생략될 것이다. 도 5와 도 13에 도시된 절차간의 차이점은 단계(510)에서의 판정 이후의 화소값 처리에 있다.
특히, 화소값 분산이 관심 화소를 포함하는 로컬 영역내의 잡음 분산보다 훨씬 큰 경우, 단계(512')에서 관심 화소의 화소값이 강화된다. 화소값 강화는 예를 들면, 화소값에 계수 β를 승산함으로써 달성된다. 계수 β의 값은 1보다 큰 양수, 예를 들면 1.2이다. 따라서, 관심 화소의 화소값은 예를 들면 원래 값의 1.2 배로 증가된다. 계수 β의 값은 1.2에 한정되지 않으며, 임의의 적절한 값일 수 있음을 알아야 한다. 더욱이, 계수를 승산하는 대신에, 화소값의 강화는 예를 들면, 사전정의된 소정의 값을 화소값에 가산함으로써 달성될 수 있다.
관심 화소를 포함하는 로컬 영역내의 화소값 분산 Vi가 잡음 분산 Vn보다 훨씬 크지 않은 경우, 화소값에 대해 특별한 작용이 적용되지 않는다. 따라서, 관심 화소의 화소값은 그 원래 값을 유지한다. 그러한 화소값 조절을 수행하는 데이터 처리 섹션(170)은 본 발명의 화소값 조절 수단의 실시예이다.
도 14는 화소값의 프로파일에서의 변화에 의한 그러한 화소값 조절의 영향을 도시한다. 도시된 바와 같이, 화소값 조절 이전의 프로파일에서의 혈액 흐름 이미지 b1 및 b2는 혈액 흐름 이미지 b1' 및 b2'에 의해 도시된 바와 같이, 전술한 화소값 조절의 결과로서 증가된 화소값을 가질 것이다. 따라서, 배경 잡음으로부터의 차이가 증대되고, 유도성이 향상된다.
도 15는 화소값 조절의 다른 영향을 도시한다. 도 15에서의 심볼 P 및 Q는 상이한 슬라이스의 두 이미지의 프로파일을 나타낸다. 프로파일 Q의 배경 잡음 레벨이 프로파일 P에서의 뚜렷한 혈액 흐름 이미지 b1의 신호 강도보다 큰 경우에 조차도, 프로파일 P에서의 혈액 흐름 이미지 b1 및 b2는, 예를 들면 전술한 화소값조절에 의해 1.2 배로 강화되어, 혈액 흐름 이미지 b1' 및 b2'이 된다. 따라서, 혈액 흐름 이미지 b1' 및 b2'은 슬라이스 q의 이미지의 잡음 레벨에 대해 또한 유도될 수 있다.
전술한 바와 같은 화소 조절 이후의 멀티슬라이스 이미지에 대해, 도 3의 흐름도의 단계(306)에서 최대 강도 투영(MIP)이 수행된다. 단계(306)에서 최대 강도 투영을 수행하는 데이터 처리 섹션(170)은 본 발명의 최대 강도 투영 수단의 실시예이다.
전술한 바와 같은 기능을 구현하는 컴퓨터를 위한 프로그램이 기록 매체상에 컴퓨터 판독가능 방법으로 기록된다. 기록 매체를 위해, 예를 들면, 자기 기록 매체, 광학 기록 매체, 자기 광학(magneto-optical) 기록 매체 중 임의의 매체 및 다른 적절한 타입의 기록 매체를 이용할 수 있다. 기록 매체는 반도체 저장 매체일 수 있다. 저장 매체는 본 명세서에서 기록 매체와 동의어이다.
도 16은 최대 강도 투영의 개념적인 도면이다. 도시된 바와 같이, 화소의 최대값은 멀티슬라이스 이미지 S1-Sm을 통과하는 시선(line of sight) E를 따라 추출되고, 추출된 값은 투영 이미지 R에 대한 화소값으로서 이용된다. 투영 이미지 R에서의 화소 수와 동일한 다수의 시선을 시선 E로서 이용한다.
전술한 바와 같은 화소값 조절에 따르면, 혈액 흐름 이미지와 잡음 사이의 차이가 모든 이미지 S1-Sm에 대해 증대되므로, 잡음에 의해 불명료해지지 않으면서, 희미한 혈액 흐름 이미지 조차도 뚜렷이 표시될 수 있다. 따라서, 혈액 흐름 이미지가 희미한 경우에도, 뚜렷한 혈액 흐름 이미지를 갖는 MIP 이미지를 획득할수 있다. 단계(308)에서, 그러한 MIP 이미지가 디스플레이 섹션(180)에 디스플레이된다.
전술한 내용은 이미지 처리가 자기 공명 촬상 장치내의 데이터 처리 섹션에 의해 수행되는 예에 대한 것이지만, 이미지 처리는 자기 공명 촬상 장치와는 분리된, EWS(engineering workstation) 또는 PC와 같은 데이터 처리 장치에 의해서도 수행될 수 있음을 쉽게 알 수 있다.
더욱이, 촬상 장치가 MRI 장치로서 기술되었지만, 촬상 장치는 그것에 한정되지 않으며, X선 CT(computed tomography) 장치, X선 촬상 장치, PET(positron emission tomography) 또는카메라와 같은 임의의 다른 타입의 촬상 장치일 수 있다.
또한, 전술한 내용은 의료용 이미지를 처리하는 예를 참조하여 기술되었지만, 처리 대상물은 의료용 이미지에 한정되지 않으며, 광학 도구에 의해 캡처된 디지털 이미지와 같은 다양한 이미지에 대해 이미지 처리가 수행될 수 있다.
본 발명은 전술한 바와 같은 바람직한 실시예를 참조하여 기술되었지만, 당업자라면 본 발명의 영역을 벗어나지 않고서도 다양한 변경 또는 대체가 가능함을 알 것이다. 따라서, 본 발명의 기술적인 영역은 전술한 실시예 뿐만 아니라, 첨부된 특허 청구 범위의 영역에 속하는 모든 실시예를 포함한다.
본 발명에 의하면, 멀티슬라이스 이미지내의 실질적인 구조와 관련된 화소와다른 화소간의 차이를 더 증대시킬 수 있고, 그와 같이 차이가 증대된 멀티슬라이스 이미지는 최대 강도 투영되므로, 희미한 신호 강도를 가지며 이미지내의 실질적인 구조와 관련된 화소의 보다 우수한 투영 이미지를 획득할 수 있다.

Claims (36)

  1. 이미지 처리 방법에 있어서,
    관심 화소(a pixel of interest)로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 단계와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 단계를 포함하는
    이미지 처리 방법.
  2. 이미지 처리 방법에 있어서,
    관심 화소로서 정의된 멀티슬라이스(multi-slice) 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 단계와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 것을 포함하는 화소값 조절을 수행하는 단계와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영(maximum intensity projection)을 수행하는 단계를 포함하는
    이미지 처리 방법.
  3. 이미지 처리 방법에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 단계와,
    상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 단계와,
    상기 가산된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 것을 포함하는 화소값 조절을 수행하는 단계와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 단계를 포함하는
    이미지 처리 방법.
  4. 제 1 항 내지 제 3 항 중 어느 한 항에 있어서,
    상기 화소값 억제는 1보다 작은 계수를 승산함으로써 수행되는 이미지 처리 방법.
  5. 제 1 항 내지 제 3 항 중 어느 한 항에 있어서,
    상기 화소값 억제는 사전정의된 수치를 감산함으로써 수행되는 이미지 처리 방법.
  6. 이미지 처리 방법에 있어서,
    관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 단계와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 단계를 포함하는
    이미지 처리 방법.
  7. 이미지 처리 방법에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 단계와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 것을 포함하는 화소값 조절을 수행하는 단계와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 단계를 포함하는
    이미지 처리 방법.
  8. 이미지 처리 방법에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 단계와,
    상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 단계와,
    상기 가산된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 것을 포함하는 화소값 조절을 수행하는 단계와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 단계를 포함하는
    이미지 처리 방법.
  9. 제 6 항 내지 제 8 항 중 어느 한 항에 있어서,
    상기 화소값 강화는 1 이상의 계수를 승산함으로써 수행되는 이미지 처리 방법.
  10. 제 6 항 내지 제 8 항 중 어느 한 항에 있어서,
    상기 화소값 강화는 사전정의된 수치를 가산함으로써 수행되는 이미지 처리 방법.
  11. 제 1 항 내지 제 3 항 또는 제 6 항 내지 제 8 항 중 어느 한 항에 있어서,
    전체 이미지에 걸쳐 정의된 다수의 로컬 영역 각각에 대한 화소값의 오차 제곱의 합(residual sums of squares)을 결정하는 단계와, 상기 오차 제곱의 합의 히스토그램을 결정하는 단계와, 이후, 상기 히스토그램의 피크를 제공하는 오차 제곱의 합에 근거하여 상기 잡음 분산을 결정하는 단계를 더 포함하는 이미지 처리 방법.
  12. 제 1 항 내지 제 3 항 또는 제 6 항 내지 제 8 항 중 어느 한 항에 있어서,
    상기 이미지는 혈액 흐름 이미지인 이미지 처리 방법.
  13. 이미지 처리 장치에 있어서,
    관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 장치를 포함하는
    이미지 처리 장치.
  14. 이미지 처리 장치에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 장치와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 장치를 포함하는
    이미지 처리 장치.
  15. 이미지 처리 장치에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 가산 장치와,
    상기 가산된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 장치와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 장치를 포함하는
    이미지 처리 장치.
  16. 제 13 항 내지 제 15 항 중 어느 한 항에 있어서,
    상기 화소값 조절 장치는 1보다 작은 계수를 승산함으로써 상기 화소값 억제를 수행하는 이미지 처리 장치.
  17. 제 13 항 내지 제 15 항 중 어느 한 항에 있어서,
    상기 화소값 조절 장치는 사전정의된 수치를 감산함으로써 상기 화소값 억제를 수행하는 이미지 처리 장치.
  18. 이미지 처리 장치에 있어서,
    관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 장치를 포함하는
    이미지 처리 장치.
  19. 이미지 처리 장치에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 장치와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 장치를 포함하는
    이미지 처리 장치.
  20. 이미지 처리 장치에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 가산 장치와,
    상기 가산된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 장치와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 장치를 포함하는
    이미지 처리 장치.
  21. 제 18 항 내지 제 20 항 중 어느 한 항에 있어서,
    상기 화소값 조절 장치는 1 이상의 계수를 승산함으로써 상기 화소값 강화를 수행하는 이미지 처리 장치.
  22. 제 18 항 내지 제 20 항 중 어느 한 항에 있어서,
    상기 화소값 조절 장치는 사전정의된 수치를 가산함으로써 상기 화소값 강화를 수행하는 이미지 처리 장치.
  23. 제 13 항 내지 제 15 항 또는 제 18 항 내지 제 20 항 중 어느 한 항에 있어서,
    상기 잡음 분산을 결정하는 장치로서, 전체 이미지에 걸쳐 정의된 다수의 로컬 영역 각각에 대한 화소값의 오차 제곱의 합을 결정하고, 상기 오차 제곱의 합의 히스토그램을 결정한 후, 상기 히스토그램의 피크를 제공하는 오차 제곱의 합에 근거하여 상기 잡음 분산을 결정하는 잡음 분산 계산 장치를 더 포함하는 이미지 처리 장치.
  24. 제 13 항 내지 제 15 항 또는 제 18 항 내지 제 20 항 중 어느 한 항에 있어서,
    상기 이미지는 혈액 흐름 이미지인 이미지 처리 장치.
  25. 기록 매체에 있어서,
    관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 기록 매체.
  26. 기록 매체에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 것을 포함하는 화소값 조절을 수행하는 기능과,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 기록 매체.
  27. 기록 매체에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과,
    상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 기능과,
    상기 가산된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 것을 포함하는 화소값 조절을 수행하는 기능과,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 기록 매체.
  28. 기록 매체에 있어서,
    관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 기록 매체.
  29. 기록 매체에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 것을 포함하는 화소값 조절을 수행하는 기능과,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 기록 매체.
  30. 기록 매체에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 기능과,
    상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 기능과,
    상기 가산된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 것을 포함하는 화소값 조절을 수행하는 기능과,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 기능을 컴퓨터가 구현하도록 하는 프로그램으로 컴퓨터 판독가능 방법에 의해 기록되는 기록 매체.
  31. 대상물(object)로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치에 있어서,
    관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 장치를 포함하는
    촬상 장치.
  32. 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 장치와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 장치를 포함하는
    촬상 장치.
  33. 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 가산 장치와,
    상기 가산된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 유지하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 억제하는 화소값 조절 장치와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 장치를 포함하는
    촬상 장치.
  34. 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치에 있어서,
    관심 화소로서 정의된 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 장치를 포함하는
    촬상 장치.
  35. 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 장치와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 장치를 포함하는
    촬상 장치.
  36. 대상물로부터 수집된 신호에 근거하여 이미지를 생성하는 촬상 장치에 있어서,
    관심 화소로서 정의된 멀티슬라이스 이미지를 구성하는 각 화소를 가지고, 상기 관심 화소가 속하는 로컬 영역내의 화소값 분산을 결정하는 분산 계산 장치와,
    상기 결정된 분산에 인접 슬라이스의 이미지내의 대응하는 관심 화소가 속하는 로컬 영역내의 화소값 분산을 가산하는 가산 장치와,
    상기 가산된 분산이 잡음 분산보다 큰 경우 상기 관심 화소의 화소값을 강화하고, 그렇지 않은 경우 상기 관심 화소의 화소값을 유지하는 화소값 조절 장치와,
    상기 화소값 조절이 수행된 멀티슬라이스 이미지에 대해 최대 강도 투영을 수행하는 최대 강도 투영 장치를 포함하는
    촬상 장치.
KR10-2001-0076051A 2000-12-05 2001-12-04 이미지 처리 장치, 이미지 처리 방법, 기록 매체 및 촬상장치 KR100458776B1 (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2000369444A JP4149126B2 (ja) 2000-12-05 2000-12-05 画像処理方法、画像処理装置および画像撮影装置
JPJP-P-2000-00369444 2000-12-05

Publications (2)

Publication Number Publication Date
KR20020044074A KR20020044074A (ko) 2002-06-14
KR100458776B1 true KR100458776B1 (ko) 2004-12-03

Family

ID=18839478

Family Applications (1)

Application Number Title Priority Date Filing Date
KR10-2001-0076051A KR100458776B1 (ko) 2000-12-05 2001-12-04 이미지 처리 장치, 이미지 처리 방법, 기록 매체 및 촬상장치

Country Status (5)

Country Link
US (1) US6985613B2 (ko)
EP (1) EP1225544A3 (ko)
JP (1) JP4149126B2 (ko)
KR (1) KR100458776B1 (ko)
CN (1) CN1287731C (ko)

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1223553A3 (en) * 2000-10-17 2003-09-24 Fuji Photo Film Co., Ltd. Apparatus for suppressing noise by adapting filter characteristics to input image signal based on characteristics of input image signal
JP4045422B2 (ja) * 2002-07-23 2008-02-13 株式会社日立メディコ 画像処理装置
US7417636B2 (en) 2003-05-08 2008-08-26 Siemens Medical Solutions Usa, Inc. Method and apparatus for automatic setting of rendering parameter for virtual endoscopy
US7616792B2 (en) * 2003-11-18 2009-11-10 Siemens Medical Solutions Usa, Inc. Scale selection for anisotropic scale-space: application to volumetric tumor characterization
FR2869723A1 (fr) * 2004-04-29 2005-11-04 Thomson Licensing Sa Element de transition sans contact entre un guide d'ondes et une ligne mocroruban
DE102004025417B4 (de) * 2004-05-24 2008-01-24 Siemens Ag Verfahren zur Verbesserung der Bildhomogenität von Bilddaten aus phasenzyklierten Steady-State-Sequenzen
US7840048B2 (en) * 2004-05-26 2010-11-23 Guardian Technologies International, Inc. System and method for determining whether there is an anomaly in data
US7496218B2 (en) * 2004-05-26 2009-02-24 Ramsay Thomas E System and method for identifying objects of interest in image data
US20090324067A1 (en) 2005-03-15 2009-12-31 Ramsay Thomas E System and method for identifying signatures for features of interest using predetermined color spaces
US20060269135A1 (en) * 2005-03-15 2006-11-30 Ramsay Thomas E System and method for identifying objects of interest in image data
US8538099B2 (en) * 2005-03-23 2013-09-17 General Electric Company Method and system for controlling image reconstruction
DE102005033851A1 (de) * 2005-07-20 2007-02-01 Sci-Worx Gmbh Verfahren zur Reduktion von Rauschartefakten in digitalen Bildern und Einrichtung hierfür
EP1934622A1 (en) * 2005-10-07 2008-06-25 Koninklijke Philips Electronics N.V. Multiple-channel transmit magnetic resonance
JP4859516B2 (ja) * 2006-04-27 2012-01-25 株式会社ソニー・コンピュータエンタテインメント 画像処理装置および画像復元方法
TWI319676B (en) * 2006-10-18 2010-01-11 Quanta Comp Inc Image processing apparatus and method
CN101472058B (zh) * 2007-12-29 2011-04-20 比亚迪股份有限公司 图像噪声去除装置和方法
EP2309285B1 (en) 2008-07-24 2019-08-21 Toshiba Medical Systems Corporation Magnetic resonance imaging apparatus for contrast enhancement of flow images
JP5414335B2 (ja) * 2008-07-24 2014-02-12 株式会社東芝 磁気共鳴イメージング装置
DE102008058740B4 (de) * 2008-11-24 2012-05-31 Siemens Aktiengesellschaft Korrektur von Artefakten in der "Time of Flight"- MR-Angiographie
CN101770636B (zh) * 2008-12-30 2012-09-26 比亚迪股份有限公司 图像噪声去除方法及装置
JP5036785B2 (ja) * 2009-10-23 2012-09-26 キヤノン株式会社 光断層画像生成方法及び光断層画像生成装置
JP5653098B2 (ja) * 2010-07-09 2015-01-14 Hoya株式会社 ノイズ除去装置
IT1403150B1 (it) 2010-11-24 2013-10-04 St Microelectronics Srl Procedimento e dispositivo per depurare dal rumore un segnale video digitale, relativo prodotto informatico.
JP5703068B2 (ja) * 2011-02-28 2015-04-15 綜合警備保障株式会社 顔検出方法および顔検出装置、ならびに、プログラム
JP5438066B2 (ja) 2011-05-25 2014-03-12 日立アロカメディカル株式会社 超音波画像処理装置およびプログラム
JP6073661B2 (ja) * 2012-11-21 2017-02-01 株式会社日立製作所 磁気共鳴イメージング装置および磁気共鳴イメージング方法
JP6647836B2 (ja) * 2014-11-10 2020-02-14 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置、画像処理装置及び画像処理方法
US11300646B2 (en) 2014-11-10 2022-04-12 Canon Medical Systems Corporation Magnetic resonance imaging apparatus, image processing apparatus, and image processing method
JP7455508B2 (ja) 2018-12-26 2024-03-26 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置および医用複素数画像処理装置
CN113506260B (zh) * 2021-07-05 2023-08-29 贝壳找房(北京)科技有限公司 一种人脸图像质量评估方法、装置、电子设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1990000052A1 (en) * 1988-07-01 1990-01-11 G.D. Searle & Co. Poly(aminoalkyl)aminocarbonyl aminodiol amino acid derivatives as anti-hypertensive agents
WO1993001055A2 (en) * 1991-07-04 1993-01-21 John Anthony Cayzer Improved label printing apparatus
US5832134A (en) * 1996-11-27 1998-11-03 General Electric Company Data visualization enhancement through removal of dominating structures

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4761819A (en) * 1987-02-27 1988-08-02 Picker International, Inc. Adaptive noise reduction filter for reconstructed images
JPH0751134B2 (ja) * 1988-12-16 1995-06-05 アロカ株式会社 超音波診断画像処理装置
US5038388A (en) * 1989-05-15 1991-08-06 Polaroid Corporation Method for adaptively sharpening electronic images
JPH0736818B2 (ja) * 1990-05-07 1995-04-26 アロカ株式会社 超音波診断画像処理装置
JP2952687B2 (ja) * 1990-07-11 1999-09-27 ジーイー横河メディカルシステム株式会社 画像処理装置
JPH0736820B2 (ja) * 1990-09-14 1995-04-26 アロカ株式会社 画像輪郭強調処理装置
US5148479A (en) 1991-03-20 1992-09-15 International Business Machines Corp. Authentication protocols in communication networks
US5233299A (en) * 1991-03-25 1993-08-03 General Electric Company Projection methods for producing two-dimensional images from three-dimensional data
JPH04329930A (ja) * 1991-04-30 1992-11-18 Shimadzu Corp 医療用診断画像作成装置
DE69331719T2 (de) * 1992-06-19 2002-10-24 Agfa Gevaert Nv Verfahren und Vorrichtung zur Geräuschunterdrückung
JPH06217958A (ja) * 1993-01-22 1994-08-09 Toshiba Medical Eng Co Ltd 磁気共鳴イメージング装置
EP0610604A1 (en) 1993-02-11 1994-08-17 Agfa-Gevaert N.V. Multiple processing of radiographic images based on a pyramidal image decomposition
EP0610607A1 (en) * 1993-02-11 1994-08-17 Agfa-Gevaert N.V. Method of evaluating image processing performed on a radiographic image
JP3350176B2 (ja) * 1993-09-24 2002-11-25 株式会社東芝 磁気共鳴イメージング装置
US5933540A (en) * 1995-05-11 1999-08-03 General Electric Company Filter system and method for efficiently suppressing noise and improving edge definition in a digitized image
US6043655A (en) * 1997-01-09 2000-03-28 Kabushiki Kaisha Toshiba MR imaging utilizing the time of flight effect
JP3896188B2 (ja) 1997-06-13 2007-03-22 株式会社日立製作所 放射線治療計画のための画像処理装置
FR2772225A1 (fr) * 1997-12-09 1999-06-04 Philips Electronics Nv Procede de traitement d'une image ou le bruit depend du signal
US6718068B1 (en) * 2000-03-10 2004-04-06 Eastman Kodak Company Noise reduction method utilizing statistical weighting, apparatus, and program for digital image processing

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1990000052A1 (en) * 1988-07-01 1990-01-11 G.D. Searle & Co. Poly(aminoalkyl)aminocarbonyl aminodiol amino acid derivatives as anti-hypertensive agents
WO1993001055A2 (en) * 1991-07-04 1993-01-21 John Anthony Cayzer Improved label printing apparatus
US5832134A (en) * 1996-11-27 1998-11-03 General Electric Company Data visualization enhancement through removal of dominating structures

Also Published As

Publication number Publication date
EP1225544A3 (en) 2004-12-15
US6985613B2 (en) 2006-01-10
KR20020044074A (ko) 2002-06-14
CN1287731C (zh) 2006-12-06
US20020094114A1 (en) 2002-07-18
EP1225544A2 (en) 2002-07-24
JP2002183709A (ja) 2002-06-28
CN1357302A (zh) 2002-07-10
JP4149126B2 (ja) 2008-09-10

Similar Documents

Publication Publication Date Title
KR100458776B1 (ko) 이미지 처리 장치, 이미지 처리 방법, 기록 매체 및 촬상장치
JP3976684B2 (ja) 画像における動きの影響を低減する方法および装置
JP5737725B2 (ja) 局所化および高度に限定された画像再構成法
US10386439B2 (en) Density guided attenuation map generation in PET/MR systems
US6842638B1 (en) Angiography method and apparatus
Chung et al. Statistical 3D vessel segmentation using a Rician distribution
US5233299A (en) Projection methods for producing two-dimensional images from three-dimensional data
EP0621546A2 (en) Magnetic resonance angiography method and apparatus employing an integration projection
JP5485663B2 (ja) 対称性検出及び画像位置合わせを用いた自動式走査計画のシステム及び方法
US7848554B2 (en) Sub-voxel motion correction for phase-contrast magnetic resonance imaging
Soltanian-Zadeh et al. A multidimensional nonlinear edge-preserving filter for magnetic resonance image restoration
KR100452645B1 (ko) 화상 처리 방법 및 장치, 기록 매체, 및 화상 촬영 장치
KR100444094B1 (ko) 이미지 처리 장치 및 그 방법, 촬상 장치 및 기록 매체
US7433086B2 (en) Edge detection and correcting system and method
US6192265B1 (en) Diagnostic image processing method
Simmons et al. Improvements to the quality of MRI cluster analysis
JP3865887B2 (ja) 画像補正方法
KR100452644B1 (ko) 영상 처리 방법과 장치, 기록 매체 및 촬상 장치
Park et al. Development of a bias field-based uniformity correction in magnetic resonance imaging with various standard pulse sequences
CN111031909A (zh) 用于校正磁共振数据内的伪影的系统和方法
JPH10165387A (ja) 医用画像診断装置における画像表示方法
US20220308146A1 (en) Identifying radio frequency interference in a magnetic resonance system
WO2022261006A1 (en) Organ segmentation in image
Joshi et al. MRI Denoising for Healthcare
Nair et al. Intensity inhomogeneity correction using modified homomorphic unsharp masking

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20121105

Year of fee payment: 9

FPAY Annual fee payment

Payment date: 20131106

Year of fee payment: 10

LAPS Lapse due to unpaid annual fee