KR101330638B1 - 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치 - Google Patents

자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치 Download PDF

Info

Publication number
KR101330638B1
KR101330638B1 KR1020120049779A KR20120049779A KR101330638B1 KR 101330638 B1 KR101330638 B1 KR 101330638B1 KR 1020120049779 A KR1020120049779 A KR 1020120049779A KR 20120049779 A KR20120049779 A KR 20120049779A KR 101330638 B1 KR101330638 B1 KR 101330638B1
Authority
KR
South Korea
Prior art keywords
data
magnetic resonance
generating
resonance image
interpolation
Prior art date
Application number
KR1020120049779A
Other languages
English (en)
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 고려대학교 산학협력단
Priority to KR1020120049779A priority Critical patent/KR101330638B1/ko
Priority to US13/673,629 priority patent/US9196062B2/en
Application granted granted Critical
Publication of KR101330638B1 publication Critical patent/KR101330638B1/ko

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5611Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection

Abstract

복수개의 코일 각각에서 획득된 복수개의 데이터들을 보간하여 복수개의 제1 보간 데이터들을 생성하는 단계, 복수개의 주파수 대역들 별로, 상기 제1 보간 데이터들을 가중시켜 가중 데이터들을 생성하는 단계, 및 상기 가중 데이터들을 이용하여 촬영 객체에 대응되는 상기 자기 공명 영상을 획득하는 단계를 포함하며, 복원 정확도를 높일 수 있는 자기 공명 영상 방법이 개시된다.

Description

자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치 {METHOD FOR GENERATING A MAGNETIC RESONANCE IMAGING AND APPARATUS FOR GENERATING A MAGNETIC RESONANCE IMAGING THEREOF}
본원 발명은 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치에 관한 것이다.
자기 공명 영상(MRI: magnetic resonance imaging) 촬영 장치는 자기장을 이용해 피사체를 촬영하는 장치로, 뼈는 물론 디스크, 관절, 신경 인대 등을 원하는 각도에서 입체적으로 보여주기 때문에 정확한 질병 진단을 위해서 널리 이용된다.
자기 공명 영상 장치는 자기 공명(MR: magnetic resonance) 신호를 획득하고, 획득된 신호를 영상으로 재구성하여 출력한다. 자기 공명 영상 장치는 RF 코일들, 영구자석 및 그래디언트 코일 등을 이용하여 자기 공명 신호를 획득한다. 자기 공명 신호가 획득될 때, RF 코일들의 접합 부위 등에서 미측정 신호가 발생할 수 있으며, 이러한 미측정 신호로 인하여 최종 자기 공명 영상을 재구성할 경우 결함과 노이즈가 발생할 수 있다. 또한, RF 코일들에서 획득된 K 공간 데이터가 자기 공명 영상으로 복원되는 과정에서, K 공간 데이터에 존재하던 노이즈가 증폭되어 나타날 수 있다.
그에 따라서, 전술한 결함과 노이즈가 제거된 자기 공명 영상을 최종적으로 출력하기 위해서, 획득된 자기 공명 신호를 캘리브레이션(calibration) 등의 영상 처리를 수행하여 보정하여야한다.
획득된 자기 공명 신호를 처리하는 자기 공명 영상 방법으로는 이미지에 기반한 센스(SENSE:Sensitivity Encoding) 기법, 및 K 공간(K-space)에 기반한 그라파(GRAPPA, Generalized Auto-calibrating Partially Parallel Acquisition) 기법 등이 있다.
센스(SENSE)와 같은 이미지 기반의 영상 방법은 이미지 공간에서 셀프 캘리브레이션(self calibration)을 통해 개별 코일에 대응되는 이미지를 분리하여 코일 민감도(coil sensitivity) 정보를 연산한다. 여기서 개별 코일 이미지는 나이퀴스트 샘플링된 K 공간의 중앙부 데이터를 역 푸리에 변환(Inverse Fourier transform)을 수행함으로써 얻어진다. 그리고, 코일 민감도 정보를 이용하여 자기 공명 영상을 재구성함에 있어서, 코일 민감도의 역행렬은 픽셀별로 수행되기 때문에 이미지에 기반한 셀프 캘리브레이션은 매우 정확한 코일 민감도 정보를 필요로 하게 된다.
따라서, K 공간 데이터의 중앙부에서 많은 수의 캘리브레이션 신호가 필요하고 이미지를 구성하는 데 걸리는 시간이 길어질 수밖에 없다. 또한, 이미지되는 피사체보다 시야각(FOV: Field Of View)이 더 작은 경우에는 이미지에 기반한 셀프 캘리브레이션은 이미지 재구성 시 잔존 엘리어싱 결함(residual aliasing artifacts)이 많이 나타내기도 한다.
또한, 그라파(GRAPPA)와 같은 K 공간 기반의 영상 방법은 셀프 캘리브레이션을 통해 캘리브레이션 신호와 인접하는 측정된 소스 신호 사이의 공간 상호작용 값(spatial correlations 또는 convolution kernels)을 계산한다. 그라파 기법은 센스와 달리 높은 코일 민감도 정보를 필요로 하지 않고 시야각을 재구성함에 있어서도 제한이 없다. 그러나, 영상 신호의 데이터가 노이즈에 의해 훼손되거나 공간 상호작용 값이 변화하는 경우에는 재구성된 이미지에 잔존 엘리어싱 결함 및 증폭된 노이즈들이 발생하는 문제가 있다.
따라서, 자기 공명 영상을 획득하는데 있어서, 영상 신호의 데이터가 노이즈에 의해 훼손되는 경우, 엘리어싱 결함 및 노이즈를 감소시킬 수 있는 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치를 제공할 필요가 있다.
본원 발명은 엘리어싱 결함 및 노이즈를 감소시킬 수 있는 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치의 제공을 목적으로 한다.
구체적으로, 본원 발명은 복원 및 획득되는 자기 공명 영상의 정확도는 유지하면서, 잔존 엘리어싱 결함 및 노이즈를 감소시킬 수 있는 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치의 제공을 목적으로 한다.
본 발명의 일 실시예에 따른 자기 공명 영상 생성 방법은 복수개의 코일을 포함하는 고주파 멀티 코일을 이용하여 자기 공명 영상을 획득하는 방법이다.
본 발명의 일 실시예에 따른 자기 공명 영상 생성 방법은 복수개의 코일을 포함하는 고주파 멀티 코일을 이용하여 자기 공명 영상을 생성하기 위한 방법에 있어서, 복수개의 코일 각각에서 획득된 복수개의 데이터들을 보간하여 복수개의 제1 보간 데이터들을 생성하는 단계, 복수개의 주파수 대역들 별로, 상기 제1 보간 데이터들을 가중시켜 가중 데이터들을 생성하는 단계, 및 상기 가중 데이터들을 이용하여 촬영 객체에 대응되는 상기 자기 공명 영상을 획득하는 단계를 포함한다.
또한, 상기 가중 데이터들을 생성하는 단계는 상기 복수개의 주파수 대역들 각각 별로, 상기 제1 보간 데이터들에 K 공간에서 가중치 적용(K-space weighting)을 수행하여 상기 가중 데이터들을 생성하는 단계를 포함할 수 있다.
또한, 상기 복수개의 주파수 대역들은 저 주파수 대역 및 고주파수 대역을 포함할 수 있다. 그리고, 상기 가중 데이터들을 생성하는 단계는 상기 제1 보간 데이터들을 로우 패스 필터링하여 상기 저 주파수 대역에 대응되는 상기 가중 데이터들을 생성하는 단계, 및 서로 다른 복수개의 가중치 값을 각각 적용하여, 상기 제1 보간 데이터들을 하이 패스 필터링하여 상기 고 주파수 대역에 대응되는 상기 가중 데이터들을 생성하는 단계를 포함할 수 있다.
또한, 상기 자기 공명 영상을 획득하는 단계는 상기 복수개의 주파수 대역들 중 적어도 하나에서, 상기 가중 데이터들을 보간하여 제2 보간 데이터들을 생성하는 단계, 상기 가중 데이터들 중 적어도 하나 및 상기 제2 보간 데이터들 중 적어도 하나를 가중 조합(weighted combination)하여, 상기 복수개의 코일 각각에 대응되는 복수개의 복원 영상을 획득하는 단계, 및 상기 복수개의 복원 영상을 이용하여 상기 자기 공명 영상을 획득하는 단계를 포함할 수 있다.
또한, 상기 제2 보간 데이터들을 생성하는 단계는 상기 복수개의 주파수 대역들 중 적어도 하나에서, 상기 가중 데이터들을 회선 보간(convolution interpolation)하여 상기 제2 보간 데이터들을 생성하는 단계를 포함할 수 있다.
또한, 상기 자기 공명 영상을 획득하는 단계는 상기 복수개의 복원 영상들을 역 고속 퓨리에 변환(Inverse Fast Fourier Transform)하는 단계, 및 상기 역 고속 퓨리에 변환된 상기 복수개의 복원 영상들에 근거하여 상기 자기 공명 영상을 획득하는 단계를 포함할 수 있다.
또한, 상기 자기 공명 영상을 획득하는 단계는 상기 역 고속 퓨리에 변환된 상기 복수개의 복원 영상을 제곱 합(sum of squares)하여, 상기 자기 공명 영상을 획득하는 단계를 포함할 수 있다.
또한, 상기 제1 보간 데이터들을 생성하는 단계는 상기 복수개의 데이터들을 회선 보간(convolution interpolation)하여 복수개의 제1 보간 데이터들을 생성하는 단계를 포함할 수 있다.
본 발명의 일 실시예에 따른 자기 공명 영상 생성 장치는 복수개의 코일 각각에서 획득된 복수개의 데이터들을 보간하여 복수개의 제1 보간 데이터들을 생성하는 데이터 획득부, 복수개의 주파수 대역들 별로, 상기 제1 보간 데이터들을 가중시켜 가중 데이터들을 생성하는 주파수 대역별 처리부, 및 상기 가중 데이터들을 이용하여 촬영 객체에 대응되는 상기 자기 공명 영상을 획득하는 영상 획득부를 포함한다.
또한, 상기 주파수 대역별 처리부는 상기 복수개의 주파수 대역들 각각 별로, 상기 제1 보간 데이터들에 K 공간에서 가중치 적용(K-space weighting)을 수행하여 상기 가중 데이터들을 생성할 수 있다.
또한, 상기 주파수 대역별 처리부는 상기 제1 보간 데이터들을 로우 패스 필터링하여 상기 가중 데이터들을 생성하는 적어도 하나의 로우 패스 필터, 및 상기 제1 보간 데이터를 입력받고 소정 가중치를 적용하여 하이 패스 필터링하는 하이 패스 필터를 다수개 포함하는 하이 패스 필터 부 포함할 수 있다. 여기서, 상기 복수개의 하이 패스 필터들 각각은 서로 다른 가중치 값을 적용받는다.
또한, 상기 데이터 획득부는 상기 복수개의 데이터들을 수신하는 데이터 수신부, 및 상기 복수개의 코일 각각에 대응되는 상기 복수개의 데이터들을 회선 보간(convolution interpolation)하여 상기 복수개의 제1 보간 데이터들을 생성하는 제1 보간부를 포함할 수 있다.
또한, 상기 영상 획득부는 상기 복수개의 주파수 대역들 중 적어도 하나에서, 상기 가중 데이터들을 보간하여 제2 보간 데이터들을 생성하는 제2 보간부, 상기 가중 데이터들 중 적어도 하나 및 상기 제2 보간 데이터들 중 적어도 하나를 이용하여, 상기 복수개의 코일 각각에 대응되는 복수개의 복원 영상을 획득하는 복원 영상 획득부, 및 상기 복수개의 복원 영상들을 이용하여 상기 자기 공명 영상을 획득하는 자기 공명 영상 획득부를 포함할 수 있다.
또한, 상기 자기 공명 영상 획득부는 상기 복수개의 복원 영상들을 역 고속 퓨리에 변환(Inverse Fast Fourier Transform)하고, 상기 역 고속 퓨리에 변환된 상기 복수개의 복원 영상에 근거하여 상기 자기 공명 영상을 획득할 수 있다.
본 발명의 다른 실시예에 따른 자기 공명 영상 획득 방법은 복수개의 코일을 포함하는 고주파 멀티 코일을 이용하여 자기 공명 영상을 획득하기 위한 방법으로, 상기 복수개의 코일에서 획득된 복수개의 데이터들을 보간하여 복수개의 제1 보간 데이터들을 생성하는 단계, 및 제1 보간 데이터들에 대응하는 데이터들을 보간하여 제2 보간 데이터들을 생성하는 단계를 포함한다.
도 1은 본 발명의 일 또는 다른 실시예에 따른 자기 공명 영상 생성 장치를 포함하는 자기 공명 영상 시스템을 나타내는 도면이다.
도 2는 도 1의 영상 촬영부를 상세히 나타내는 도면이다.
도 3은 본 발명의 일 실시예에 따른 자기 공명 영상 생성 장치를 나타내는 도면이다.
도 4는 본 발명의 일 실시예에 따른 자기 공명 영상 생성 방법을 나타내는 도면이다.
도 5는 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 장치를 나타내는 도면이다.
도 6은 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 방법을 나타내는 도면이다.
도 7은 도 5의 주파수 대역별 처리부(520)의 일 예를 설명하기 위한 도면이다.
도 8은 그라파 기법에 따라 획득된 자기 공명 영상을 나타내는 도면이다.
도 9는 본 발명의 일 또는 다른 실시예에 따른 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치를 통하여 획득된 자기 공명 영상을 나타내는 도면이다.
도 10은 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 방법을 나타내는 도면이다.
이하에서는 첨부된 도면을 참조하여 본 발명의 일 또는 다른 실시예에 따른 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치를 상세히 설명한다.
도 1은 본 발명의 일 또는 다른 실시예에 따른 자기 공명 영상 생성 장치를 포함하는 자기 공명 영상 시스템을 나타내는 도면이다.
도 1을 참조하면, 자기 공명 영상 시스템(100)은 영상 촬영부(110) 및 발명의 일 또는 다른 실시예에 따른 자기 공명 영상 생성 장치(120)를 포함한다.
영상 촬영부(110)는 그래디언트 코일(미도시) 및 고주파 코일(RF coil: Radio Frequency coil)을 다수개 포함하는 멀티 코일(미도시)을 포함하며, 그래디언트 코일을 구동하여 조절되는 K 공간(K-space)에 맞춰, 멀티 코일들을 이용하여 촬영 객체를 촬영한 데이터들을 생성한다. 영상 촬영부(110)는 이하에서 도 2를 참조하여 상세히 설명한다. 영상 촬영부(110)에서 생성된 데이터는 K 공간 데이터가 될 수 있다. 여기서, K 공간 데이터는 영상을 만들 수 있는 러 데이터(raw data)의 집합을 의미한다.
자기 공명 영상 생성 장치(120)는 영상 촬영부(110)에서 전송되는 복수개의 데이터들을 수신한다. 여기서, 수신되는 복수개의 데이터들은 영상 촬영부(110)의 멀티 코일에 포함되는 복수개의 코일 각각에서 획득된 K 공간의 영상 정보들을 포함하는 이미지들이 될 수 있다.
자기 공명 영상 생성 장치(120)는 수신된 복수개의 데이터들을 신호 처리하여 자기 공명 영상을 생성한다. 여기서, 신호 처리는 미측정된 신호를 복원하고 잔존 엘리어싱 결함(residual aliasing artifacts) 및 노이즈를 제거하기 위한 영상 처리 동작을 포함할 수 있다. 자기 공명 영상 생성 장치(120)의 영상 처리 동작은 이하에서 도 3을 참조하여 상세히 설명한다.
도 2는 도 1의 영상 촬영부를 상세히 나타내는 도면이다. 도 2에 도시된 영상 촬영부(200)는 도 1의 영상 촬영부(110)에 대응된다. 따라서, 도 1에서와 중복되는 설명은 생략한다.
도 2를 참조하면, 영상 촬영부(200)는 MRI 보어(MRI bore)(210) 및 제어부(230)를 포함한다.
MRI 보어(210)는 주 자석(main magnet), 2차 자기장(secondary magnetic field) 장치 및 환자 침대(patient bed)를 포함한다.
주 자석(212)은 주 자장을 형성하는 자석으로 영구자석(permanent magnet), 또는 초전도 전자석(superconducting electromagnets) 등을 포함할 수 있으며, 원통형의 보어 자석(bore magnet) 형태를 가질 수 있다.
2차 자기장 장치는 부 자장을 형성하기 위한 장치로, 그래디언트 코일(gradient coil)(214), 고주파 멀티 코일(216, 217, 218) 등을 포함할 수 있다.
그래디언트 코일(214)은 x, y, z 방향으로의 경사 자장을 형성하며, 촬영 객체에서 촬영하고자 하는 단면인 FoV(Field of View)를 설정하는데 이용된다.
고주파 멀티 코일(216, 217, 218)은 복수개의 코일(216, 217, 218)을 포함하며, 공명 현상을 유도해내는데 필요한 전자파(radio wave)를 전이시키기 위해 구동된다. 복수개의 코일(216, 217, 218)이 구동되어 발생하는 고주파 신호(RF signal)를 수신하여 K 공간 이미지를 생성한다.
제어부(230)는 MRI 보어(210)에 포함되는 자석들을 구동시킨다. 그리고, 고주파 신호를 MRI 보어(210)로 출력하고, 출력된 고주파 신호에 대응하여 생성되는 고주파 신호를 수신한다. 그리고, 수신된 고주파 신호에 근거하여, 복수개의 코일(216, 217, 218)에 대응되는 복수개의 데이터들을 생성한다.
구체적으로, 제어부(230)는 전송부(232) 및 수신부(234)를 포함할 수 있다. 전송부(232)는 고주파(RF) 신호를 생성하여 MRI 보어(210)로 출력한다. MRI 보어(210)는 고주파 전송부(232)에서 생성된 고주파 신호에 따라서 자기 공명(MR: Magnetic Resonance) 신호를 생성한다.
수신부(234)는 MRI 보어(210)에서 생성되는 자기 공명 신호를 고주파 신호로써 수신한다. 그리고, 수신된 고주파 신호에 근거하여 러 데이터들을 생성한다. 여기서, 러 데이터는 전술한 K 공간의 데이터가 되며, 고주파 멀티 코일(216, 217, 218)이 n 개의 코일들을 포함할 경우, 수신부(234)는 각각의 코일들에 대응하는 n 개의 데이터를 생성할 수 있다.
수신부(234)에서 생성된 n 개의 데이터들은 자기 공명 영상 생성 장치(도 1의 120)로 전송되며, 최종적인 자기 공명 영상을 생성하는데 이용된다.
도 3은 본 발명의 일 실시예에 따른 자기 공명 영상 생성 장치를 나타내는 도면이다. 도 3의 자기 공명 영상 생성 장치(300)는 도 1의 자기 공명 영상 생성 장치(120)에 대응되므로, 도 1에서와 중복되는 설명은 생략한다.
도 3을 참조하면, 자기 공명 영상 생성 장치(300)는 데이터 획득부(310), 주파수 대역별 처리부(330) 및 영상 획득부(350)를 포함하며, 복수개의 코일을 포함하는 멀티 코일에서 획득된 복수개의 데이터들을 수신하고, 수신된 복수개의 데이터를 이용하여 자기 공명 영상을 생성한다. 여기서, 자기 공명 영상 생성 장치(300)가 수신하는 복수개의 데이터는 전술한 K 공간에서 생성된 러 데이터가 될 수 있다.
데이터 획득부(310)는 복수개의 코일 각각에서 획득된 복수개의 데이터들을 보간(interpolation)하여 복수개의 제1 보간 데이터들을 생성한다. 여기서, 복수개의 코일은 전술한 영상 촬영부(110, 200)에 포함되는 멀티 코일에 포함되는 코일들로, 복수개의 코일 각각은 대응되는 K 공간의 데이터들을 생성한다.
이하에서는, 데이터 획득부(310)가 n 개의 코일들 각각에서 획득된 n 개의 데이터를 수신하고, 그에 대응되는 n 개의 제1 보간 데이터들을 생성하는 경우를 예로 들어 설명한다.
구체적으로, 데이터 획득부(310)는 복수개의 코일 각각에서 획득된 복수개의 데이터들을 수신할 수 있다. 또한, 데이터 획득부(310)는 복수개의 코일 각각에 대응되는 고주파 신호들만을 수신할 수도 있다. 이 경우, 수신된 고주파 신호를 K 공간에서 가변 밀도 샘플링하여 데이터를 획득한다. 예를 들어, 소정 코일에 대응되는 K 공간에서, K 공간의 중심부에는 나이퀴스트 샘플링(Nyquist sampling)을 통하여 데이터를 획득하고 K 공간의 주변부에는 언더 샘플링(under sampling)을 통하여 데이터를 획득하여, 일 코일에 대응되는 일 데이터를 생성할 수 있다.
여기서, 소정 데이터는 K 공간에 분포하는 영상 신호들로 구성될 수 있다. 또한, K 공간은 일 축이 주파수 엔코딩(frequency encoding) 정보를 포함하며 일 축에 직교하는 다른 축이 위상 엔코딩(phase encoding) 정보를 포함하는 공간을 나타낸다.
데이터 획득부(310)는 전술한 바와 같이, n개의 코일들에 대응되는 n 개의 K 공간의 데이터들을 외부의 영상 촬영부(110, 200)에서 수신할 수 도 있고, n 개의 코일들 각각에 대응되는 고조파 신호들을 수신받고, 수신된 고조파 신호에 근거하여 n 개의 K 공간의 데이터들을 자체적으로 생성할 수 도 있다.
또한, 데이터 획득부(310)는 n 개의 데이터들을 보간(interpolation)하여 n 개의 제1 보간 데이터들을 생성한다. 여기서, 보간은 n 개의 데이터들 각각을 셀프 캘리브레이션(self calibration)하여 수행될 수 있으며, 구체적으로, 회선 보간(convolution interpolation) 기법을 이용하여 수행될 수 있다. 그에 따라서, 영상 촬영 시 손실된 영상 정보 및 미측정된 영상 정보가 복원된, n 개의 코일들에 각각 대응되는 n 개의 제1 보간 데이터들이 생성된다.
주파수 대역별 처리부(330)는 복수개의 주파수 대역들 별로, 데이터 획득부(310)에서 획득된 제1 보간 데이터들을 가중(weighting)시켜 가중 데이터들을 생성한다.
구체적으로, 주파수 대역별 처리부(330)는 복수개의 주파수 대역들 별로, 상기 제1 보간 데이터들에 K 공간 가중치 적용(K-space weighting)을 수행하여 가중 데이터들을 생성한다. 구체적으로, 주파수 대역별 처리부(330)는 복수개의 주파수 대역들 별로, n 개의 제1 보간 데이터들에 대응되는 n 개의 가중 데이터들을 생성한다.
영상 획득부(350)는 주파수 대역별 처리부(330)에서 생성된 복수개의 주파수 대역 별 가중 데이터들을 이용하여 촬영 객체에 대응되는 자기 공명 영상을 획득한다.
도 4는 본 발명의 일 실시예에 따른 자기 공명 영상 생성 방법을 나타내는 도면이다.
본 발명의 일 실시예에 따른 자기 공명 영상 생성 방법(400)은 도 1 내지 도 3을 참조하여 설명한 본 발명의 일 실시예에 따른 자기 공명 영상 생성 장치(120, 300)의 상세 동작 구성과 동일한 구성상 특징을 포함한다. 따라서, 도 1 내지 도 3과 중복되는 설명은 생략한다. 이하에서는, 도 3에서 도시한 자기 공명 영상 생성 장치(300)를 참조하여, 본 발명의 일 실시예에 따른 자기 공명 영상 생성 방법(400)을 설명한다.
본 발명의 일 실시예에 따른 자기 공명 영상 생성 방법(400)은 고주파 멀티 코일(RF multi coil)을 이용하여 자기 공명 영상을 획득하는 방법이다.
본 발명의 일 실시예에 따른 자기 공명 영상 생성 방법(400)은 복수개의 코일 각각에서 획득된 복수개의 데이터들을 보간(interpolation)하여 복수개의 제1 보간 데이터들을 생성한다(420 단계). 420 단계는 데이터 획득부(310)에서 수행될 있다.
복수개의 주파수 대역들 별로, 420 단계에서 생성된 제1 보간 데이터들을 가중시켜 가중 데이터들을 생성한다(430 단계). 430 단계는 주파수 대역별 처리부(330)에서 수행될 있다.
구체적으로, 430 단계는 복수개의 주파수 대역들 별로, 제1 보간 데이터들에 K 공간 가중치 적용(K-space weighting)을 수행하여 가중 데이터들을 생성할 수 있다.
430 단계에서 생성된 가중 데이터들을 이용하여 촬영 객체에 대응되는 자기 공명 영상을 획득한다(440 단계). 440 단계는 영상 획득부(350)에서 수행될 수 있다.
본 발명의 일 실시예에 따른 자기 공명 영상 생성 방법(400)은 이하에서 도 6을 참조하여 상세히 설명한다.
도 5는 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 장치를 나타내는 도면이다. 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 장치(500)는 도 3에서 도시한 본 발명의 일 실시예에 따른 자기 공명 영상 생성 장치(300)에 대응된다. 구체적으로, 도 5의 데이터 획득부(510), 주파수 대역별 처리부(520) 및 영상 획득부(530)는 각각 도 3의 데이터 획득부(310), 주파수 대역별 처리부(330) 및 영상 획득부(350)에 대응된다. 따라서, 도 3에서와 중복되는 설명은 생략한다.
도 5를 참조하면, 데이터 획득부(510)는 데이터 수신부(511) 및 제1 보간부(513)를 포함할 수 있다.
데이터 수신부(511)는 복수개의 코일 각각에서 획득된 복수개의 데이터들을 수신할 수 있다. 또한, 데이터 수신부(511)는 복수개의 코일 각각에 대응되는 고주파 신호들만을 수신할 수도 있으며, 수신된 고주파 신호를 K 공간에서 가변 밀도 샘플링하여 데이터를 획득한다. 그에 따라서, 데이터 수신부(511)는 n 개의 코일들에 대응되는 n 개의 데이터들을 제1 보간부(513)로 전송한다.
제1 보간부(513)는 데이터 수신부(511)에서 전송된 n 개의 데이터들을 셀프 캘리브레이션 등을 통하여 보간하여 출력한다.
주파수 대역별 처리부(520)는 적어도 하나의 로우 패스 필터(521) 및 하이 패스 필터부(525)를 포함한다. 도 5에서는 주파수 대역별 처리부(520)가 하나의 로우 패스 필터(521)를 포함하는 경우를 예로 들어 도시하였다.
구체적으로, 로우 패스 필터(521)는 n 개의 코일들에 대응되는 n 개의 제1 보간 데이터들을 로우 패스 필터링(low pass filtering)하여, 저 주파수 대역에 해당하며 n 개의 코일들에 대응되는 n 개의 가중 데이터들(560)을 생성한다. 구체적으로, 로우 패스 필터(521)에는 가중치 값이 1로 적용되어, n 개의 제1 보간 데이터들이 가중되지 않고 로우 패스 필터링만 되어 출력될 수 있다.
도 5의 도시예에 있어서, 가중 데이터들(560)에 있어서, 맨 앞에 도시된 가중 데이터(561)는 제1 코일에서 획득된 데이터에 대응되는 가중 데이터이며, 맨 마지막에 도시된 가중 데이터(562)는 제 n 코일에서 획득된 데이터에 대응되는 가중 데이터를 나타낸다.
하이 패스 필터부(525)는 복수개의 하이 패스 필터들을 포함하며, 도 5에서는 제1 내지 제 m 하이 패스 필터들(526, 527, 528)을 포함하는 경우를 예로 들어 도시하였다. 여기서, 제1 내지 제 m 하이 패스 필터들(526, 527, 528) 각각은 n 개의 제1 보간 데이터들을 고 주파수 대역 필터링하며, 제1 내지 제 m 하이 패스 필터(526, 527, 528) 각각은 서로 다른 가중치 값을 적용하여 n 개의 제1 보간 데이터들을 필터링한다. 여기서, m 개의 하이 패스 필터들(526, 527, 528)에 적용되는 m 개의 가중치 값들은 0 에서 1 사이의 서로 다른 값을 가질 수 있다.
그에 따라서, 하이 패스 필터부(525)는 고 주파수 대역에 해당하며 m 개의 가중치 값들이 각각 적용된 m 개 세트의 가중 데이터들을 출력할 수 있다. 예를 들어, 제1 하이 패스 필터(526)는 제1 가중치 값이 적용되며, 고 주파수 대역 및 제1 가중치에 대응되는 n 개의 가중 데이터들(565)을 생성한다. 그리고, 제2 하이 패스 필터(527)는 제2 가중치 값이 적용되며, 고 주파수 대역 및 제2 가중치에 대응되는 n 개의 가중 데이터들(570)을 생성한다. 또한, 제m 하이 패스 필터(528)는 제m 가중치 값이 적용되며, 고 주파수 대역 및 제m 가중치에 대응되는 n 개의 가중 데이터들(575)을 생성한다. 구체적으로, 하이 패스 필터부(525)에서 생성된 가중 데이터들(565, 570, 575) 각각에 있어서, n 개의 가중 데이터들 각각은 n 개의 코일에서 획득된 데이터들에 대응된다.
로우 패스 필터(521) 및 하이 패스 필터부(525)는 이하에서 도 8을 참조하여 상세히 설명한다.
영상 획득부(530)는 제2 보간부(535), 복원 영상 획득부(540) 및 자기 공명 영상 획득부(550)를 포함할 수 있다.
제2 보간부(535)는 복수개의 주파수 대역들 중 적어도 하나에서, 가중 데이터들을 보간하여 제2 보간 데이터들을 생성한다.
구체적으로, 제2 보간부(535)는 하이 패스 필터부(525)에서 출력되는 m 개 세트의 가중 데이터들을 보간(interpolation)하여 m 개 세트의 제2 보간 데이터들을 출력한다. 여기서, 보간은 회선 보간(convolution interpolation) 기법을 이용하여 수행될 수 있다.
복원 영상 획득부(540)는 복수개의 코일 각각 별로, 가중 데이터들 중 적어도 하나 및 제2 보간 데이터들 중 적어도 하나를 이용하여 복수개의 복원 영상을 획득한다. 구체적으로, 복원 영상 획득부(540)는 가중 데이터들 중 적어도 하나 및 제2 보간 데이터들 중 적어도 하나를 가중 조합(weighted combination)하여, 복수개의 코일 각각에 대응되는 복수개의 복원 영상을 획득할 수 있다.
구체적으로, 복원 영상 획득부(540)는 n 개의 코일에 대응되는 n 개의 복원 영상을 각각 생성하는 제1 코일 조합부(541) 내지 제 n 코일 조합부(542)를 포함한다. 제1 코일 조합부(541) 내지 제 n 코일 조합부(542) 각각은 소정 코일에 대응되는 가중 데이터 및 제2 보간 데이터들을 조합(combination)하여, 소정 코일에 대응되는 복원 영상을 획득한다.
예를 들어, 제1 코일 조합부(541)는 로우 패스 필터(521)에서 출력된 가중 데이터들 중 제1 코일에 대응되는 가중 데이터(561), 제1 하이패스 필터(526) 및 제2 보간부(535)를 통과하여 생성된 제2 보간 데이터들(565) 중 제1 코일에 대응되는 제2 보간 데이터(566), 및 제2 하이패스 필터(527) 및 제2 보간부(535)를 통과하여 생성된 제2 보간 데이터들(570) 중 제1 코일에 대응되는 제2 보간 데이터(571) 내지 제n 하이패스 필터(528) 및 제2 보간부(535)를 통과하여 생성된 제2 보간 데이터들(575) 중 제1 코일에 대응되는 제2 보간 데이터(576)를 조합하여, 제1 코일에서 획득된 영상에 대응하는 복원 영상을 획득한다.
또한, 예를 들어, 제n 코일 조합부(542)는 로우 패스 필터(521)에서 출력된 가중 데이터들 중 제n 코일에 대응되는 가중 데이터(562), 제1 하이패스 필터(526) 및 제2 보간부(535)를 통과하여 생성된 제2 보간 데이터들(565) 중 제n 코일에 대응되는 제2 보간 데이터(567), 및 제2 하이패스 필터(527) 및 제2 보간부(535)를 통과하여 생성된 제2 보간 데이터들(570) 중 제n 코일에 대응되는 제2 보간 데이터(572) 내지 제n 하이패스 필터(528) 및 제2 보간부(535)를 통과하여 생성된 제2 보간 데이터들(575) 중 제n 코일에 대응되는 제2 보간 데이터(577)를 조합하여, 제n 코일에서 획득된 영상에 대응하는 복원 영상을 획득한다.
자기 공명 영상 획득부(550)는 복수개의 코일에 대응하는 복수개의 복원 영상들을 이용하여 자기 공명 영상을 획득한다.
구체적으로, 자기 공명 영상 획득부(550)는 n 개의 코일들에 대응하는 n 개의 복원 영상을 역 고속 퓨리에 변환(Inverse Fast Fourier Transform)하고, 역 고속 퓨리에 변환된 n 개의 복원 영상에 근거하여 자기 공명 영상을 획득한다. 구체적으로, 자기 공명 영상 획득부(550)는 n 개의 코일들에 대응하는 n 개의 복원 영상을 역 고속 퓨리에 변환(Inverse Fast Fourier Transform)하고, 역 고속 퓨리에 변환된 n 개의 복원 영상을 제곱 합(sum of squares)하여 자기 공명 영상을 획득할 수 있다.
도 6은 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 방법을 나타내는 도면이다. 또한, 도 6의 620, 630 및 640 단계 각각은 도 4의 420, 430 및 440 단계에 대응된다. 따라서, 도 4에서와 중복되는 설명은 생략한다. 또한, 도 6에 도시된 자기 공명 영상 생성 방법(600)은 도 5를 참조하여 설명한 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 장치(500)의 상세 동작 구성과 동일한 구성상 특징을 포함한다. 따라서, 도 5와 중복되는 설명은 생략한다. 이하에서는, 도 5에서 도시한 자기 공명 영상 생성 장치(500)를 참조하여, 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 방법(600)을 설명한다.
자기 공명 영상 생성 방법(600)은 복수개의 코일 각각에서 획득된 복수개의 데이터들(610)을 외부적으로 수신하거나 자체적으로 생성할 수 있다. 구체적으로, 데이터들(610)들은 전술한 바와 같이 K 공간의 영상 정보들로 구성되며, n 개의 코일들 각각에서 획득된 n 개의 K 공간의 데이터들을 포함한다.
데이터들(610)은 데이터 수신부(511)에서 수신 또는 획득될 수 있다.
예를 들어, 데이터들(610) 각각은, x 축에 주파수 엔코딩(frequency encoding) 정보를 포함하며, y 축에 위상 엔코딩(phase encoding) 정보를 포함한다. 도 6에서는 맨 앞에 도시된 데이터(611)가 제1 코일에서 획득된 K 공간의 데이터이며, 맨 뒤에 도시된 데이터(612)가 제 n 코일에서 획득된 K 공간의 데이터인 경우를 예로 들어 도시하였다.
620 단계는 420 단계와 동일하다. 620 단계는 제1 보간부(513)에서 수행될 수 있다.
복수개의 주파수 대역들 별로, 620 단계에서 생성된 제1 보간 데이터들을 가중시켜 가중 데이터들을 생성한다(630 단계). 여기서, 복수개의 주파수 대역들은 저 주파수 대역 및 고주파수 대역을 포함할 수 있다.
구체적으로, 630 단계에 있어서, 제1 보간 데이터들을 로우 패스 필터링하여 저 주파수 대역에 대응되는 가중 데이터들을 생성한다(631 단계). 또한, 서로 다른 복수개의 가중치 값을 각각 적용하여, 제1 보간 데이터들을 하이 패스 필터링하여 고 주파수 대역에 대응되는 가중 데이터들을 생성한다(632 단계). 631 단계는 로우 패스 필터(521)에서 수행될 수 있으며, 632 단계는 하이 패스 필터부(525)에서 수행될 수 있다.
640 단계는 440 단계에 대응되며, 630 단계에서 생성된 가중 데이터들을 이용하여, 촬영 객체에 대응되는 자기 공명 영상을 획득한다.
구체적으로, 640 단계에 있어서, 복수개의 주파수 대역들 중 적어도 하나에서, 가중 데이터들을 보간하여 제2 보간 데이터들을 생성한다(641 단계). 641 단계의 동작은 제2 보간부(535)에서 수행될 수 있다.
도 6에 도시된 가중 데이터들(651) 및 제2 보간 데이터들(653, 654, 655) 각각은 도 5에 도시된 가중 데이터들(560) 및 제2 보간 데이터들(565, 570, 575)에 동일 대응된다.
복수개의 코일 각각 별로, 가중 데이터들 중 적어도 하나 및 상기 제2 보간 데이터들을 가중 조합(weighted combination)하여 복수개의 복원 영상을 획득한다(660 단계). 660 단계의 동작은 복원 영상 획득부(540)에서 수행될 수 있다. 구체적으로, 제1 코일에 대응되는 제1 복원 영상(662)은 제1 코일 조합부(541)에서 생성될 수 있으며(661 단계), 제n 코일에 대응되는 제n 복원 영상(664)은 제n 코일 조합부(542)에서 생성될 수 있다(663 단계).
복수개의 코일에 대응하는 복수개의 복원 영상을 이용하여 자기 공명 영상(680)을 획득한다(670 단계). 670 단계의 동작은 자기 공명 영상 획득부(550)에서 수행될 수 있다.
도 7은 도 5의 주파수 대역별 처리부(520)의 일 예를 설명하기 위한 도면이다. 도 7에 도시된 그래프에 있어서, x 축은 주파수를 나타내며, y 축은 필터링된 신호의 진폭(filtered amplitude) 값을 나타낸다. y 축에 도시된 필터링된 진폭(filtered amplitude) 값에 있어서, 원 신호의 진폭이 100% 유지되어 필터링 된 경우를 1 값으로 나타내고, 원 신호의 진폭이 0% 유지되어 필터링 된 경우를 0 값으로 나타낸다. 또한, 도 7에서는, 하이 패스 필터부(525)에 m=3 개의 하이 패스 필터들이 포함되는 경우를 예로 들어 도시하였다. 따라서, 이하에서는 제 m 하이 패스 필터(528)는 제3 하이 패스 필터(528)인 경우를 예로 들어 설명한다.
도 7을 참조하면, 710 곡선은 로우 패스 필터(521)의 필터 특성 곡선을 나타내고, 720 곡선은 제1 하이 패스 필터(526)의 필터 특성 곡선을 나타낸다. 그리고, 730 곡선은 제2 하이 패스 필터(527)의 필터 특성 곡선을 나타내며, 740 곡선은 제3 하이 패스 필터(528)의 필터 특성 곡선을 나타낸다.
일예로서, 로우 패스 필터(521)에서 출력되는 신호값은 이하의 [수학식 1]에 의하여 산출될 수 있으며, 하이 패스 필터(526, 527, 또는 528)에서 출력되는 신호값은 이하의 [수학식 2]에 의하여 산출될 수 있다.
[수학식 1]
Figure 112012037524615-pat00001
[수학식 2]
Figure 112012037524615-pat00002
여기서, Kx는 K 공간 데이터의 x 축 좌표의 신호값을 나타내고, Ky는 K 공간 데이터의 y축 좌표의 신호값을 나타낸다. 그리고, c 는 필터의 컷 오프 주파수(cutoff frequency)를 조절하는 계수 값이고, w 는 필터 특성 곡선의 천이 구간에서의 기울기(smoothness of transition period)를 조절하는 계수 값이다. 그리고, [수학식 1]의 'LPF' 는 로우 패스 필터(521)에서 필터링되어 출력되는 신호 값이고, [수학식 2]의 'HPF' 는 하이 패스 필터(526, 527, 또는 528)에서 필터링되어 출력되는 신호 값이다.
도시된 [수학식 1] 및 [수학식 2]는 주파수 대역별 처리부(520)에 포함되는 로우 패스 필터 및 하이 패스 필터의 특성에 따라서 다르게 설계될 수 도 있다.
도 7 및 [수학식 1]을 참조하면, 로우 패스 필터(521)는 저 주파수 대역의 신호 성분을 통과시켜 출력하며, 고 주파수 대역의 신호를 차단시킨다.
도 7 및 [수학식 2]를 참조하면, 하이 패스 필터(526, 527, 또는 528)는 저 주파수 대역의 신호를 차단시키고, 고 주파수 대역의 신호를 통과시켜 출력한다. 하이 패스 필터(526, 527, 또는 528)에 있어서, 저 주파수 대역에서의 신호를 차단하는 정도는 하이 패스 필터(526, 527, 또는 528)에 적용되는 가중치 값(Vw)에 따라서 달라질 수 있다. 또한, [수학식 1] 및 [수학식 2]의 계수 c 및 w 값에 따라서 저 주파수 대역의 신호 차단 정도가 달라질 있으며, 계수 c 및 w 값을 조절하여 하이 패스 필터(526, 527, 또는 528)에 적용되는 가중치 값(Vw)을 조절할 있다.
구체적으로, 저 주파수 대역의 신호 차단 정도가 클수록 하이 패스 필터(526, 527, 또는 528)에 적용되는 가중치 값(Vw)이 커질 수 있다.
예를 들어, 전술한 430 단계에서, 소정 주파수 대역, 예를 들어, 고 주파수 대역,에서 제1 보간 데이터들을 가중시키는 데 있어서, 제1 보간 데이터들에 더 큰 가중치 값을 적용하여 제1 보간 데이터들을 가중시키면, 고 주파수 대역의 신호 성분이 더 많이 포함되어 출력된다.
예를 들어, 저 주파수 대역의 신호를 대부분 차단하는 제3 하이 패스 필터(528)에 적용되는 가중치 값(Vw)은 1 값을 가질 수 있고, 제2 하이 패스 필터(527)에 적용되는 가중치 값(Vw)은 0.8, 제1 하이 패스 필터(526)에 적용되는 가중치 값(Vw)은 0.6 값을 가질 수 있다.
또한, 각각의 하이 패스 필터(526, 527, 528)에 적용되는 계수 c 및 w 값, 또는 가중 값(Vw)은 제2 보간부(535)에서 출력되는 신호들의 정확도 또는 보간 정도, 또는 자기 공명 영상 획득부(550)에서 생성되는 자기 공명 영상의 화질 등에 따라서 실험적으로 최적화 될 수 있다.
또한, 로우 패스 필터(521)가 필터링하는 저 주파수 대역의 범위는 도 7에 도시된 로우 패스 필터의 특성 곡선(710)에서 나타나는 필터링 대역인 f3 내지 f4 의 주파수 값, 또는 f1 내지 f2의 주파수 값에 따라서 달라질 수 있다. 또한, 하이 패스 필터(526, 527, 또는 528) 필터링하는 고 주파수 대역의 범위는 도 7에 도시된 하이 패스 필터의 특성 곡선(720, 730, 740)에서 나타내는 필터링 대역인 f5 이상이 되는 주파수 대역 또는 f6 이상이 되는 주파수 대역으로 결정될 수 있다.
도 8은 그라파 기법에 따라 획득된 자기 공명 영상을 나타내는 도면이다.
도 8의 (a)를 참조하면, 종래의 그라파(GRAPPA) 기법에 따라서 멀티 코일 데이터들을 이용하여 자기 공명 영상(810)을 획득한 경우, 획득된 자기 공명 영상(810)을 나타낸다. 도 8의 (b)는 자기 공명 영상(810)의 일 부분(820)을 확대한 영상이다. 도 8의(b)에 도시된 자기 공명 영상(750)은 자기 공명 영상(820)에 동일 대응된다.
도 8의 (b)를 참조하면, 그라파(GRAPPA) 기법에 따라서 생성된 자기 공명 영상(750)에는 노이즈 및 엘리어싱 결함(aliasing artifacts)(851, 852)이 존재한다.
도 9는 본 발명의 일 또는 다른 실시예에 따른 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치를 통하여 획득된 자기 공명 영상을 나타내는 도면이다.
도 9의 (a)를 참조하면, 본 발명의 일 또는 다른 실시예에 따른 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치를 통하여 획득된 자기 공명 영상(910)이 도시된다. 도 9의 (b)는 자기 공명 영상(910)의 일 부분(920)을 확대하여 도시한 도면이다.
도 9의 (b)를 참조하면, 자기 공명 영상(950)의 951 부분과 도 8의 (b)의 노이즈 및 엘리어싱 결함(851, 852)이 존재하는 부분을 비교하면, 자기 공명 영상(950)에는 노이즈 및 엘리어싱 결함이 거의 존재하지 않는다.
도 10은 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 방법을 나타내는 도면이다. 본 발명의 다른 실시예에 따른 자기 공명 영상 생성 방법(1000)은 복수개의 코일을 포함하는 고주파 멀티 코일을 이용하여 자기 공명 영상을 획득하기 위한 방법이다.
도 10을 참조하면, 자기 공명 영상 생성 방법(1000)은 복수개의 코일에서 획득된 복수개의 데이터들을 보간(interpolation)하여 복수개의 제1 보간 데이터들을 생성한다(1010 단계). 1010 단계의 동작은 제1 보간부(513)에서 수행될 수 있다. 또한, 1010 단계는 도 6에서 설명한 620 단계와 동일 대응되므로, 상세한 설명은 생략한다.
그리고, 1010 단계에서 생성된 제1 보간 데이터들에 대응하는 데이터들을 보간(interpolation)하여 제2 보간 데이터들을 생성한다(1020 단계). 1020 단계의 동작은 제2 보간부(535)에서 수행될 수 있다. 또한, 1020 단계는 도 6에서 설명한 641 단계와 동일 대응되므로, 상세한 설명은 생략한다.
전술한 바와 같이, 본 발명의 일 또는 실시예에 따른 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치는 복수개의 주파수 대역들 별로 보간 데이터를 처리하고, 복수개의 주파수 대역들 각각에서 처리된 데이터들을 이용하여 자기 공명 영상을 획득함으로써, 자기 공명 영상을 복원하는 과정에서 발생하는 노이즈 및 엘리어싱 결함을 감소시킬 수 있다. 그에 따라서, 자기 공명 영상을 더욱 정확하게 복원해 낼 수 있다.
이상의 설명은 본 발명의 일 실시예에 불과할 뿐, 본 발명이 속하는 기술 분야에서 통상의 지식을 가진자는 본 발명의 본질적 특성에서 벗어나지 않는 범위에서 변형된 형태로 구현할 수 있을 것이다. 따라서, 본 발명의 범위는 전술한 실시예에 한정되지 않고 특허 청구 범위에 기재된 내용과 동등한 범위내에 있는 다양한 실시 형태가 포함되도록 해석되어야 할 것이다.
110, 200: 영상 촬영부
120, 300, 500: 자기 공명 영상 생성 장치
210: 보어
212: 주 자석
214: 그래디언트 코일
230: 제어부
232: 전송부
234: 수신부
310, 510: 데이터 획득부
330, 520: 주파수 대역별 처리부
350, 530: 영상 획득부
511: 데이터 수신부
513: 제1 보간부
521: 로우 패스 필터
525: 하이 패스 필터부
535: 제2 보간부
540: 복원 영상 획득부
550: 자기 공명 영상 획득부

Claims (19)

  1. 복수개의 코일을 포함하는 고주파 멀티 코일을 이용하여 자기 공명 영상을 생성하기 위한 방법에 있어서,
    상기 복수개의 코일 각각에서 획득된 복수개의 데이터들을 보간하여 복수개의 제1 보간 데이터들을 생성하는 단계;
    복수개의 주파수 대역들 별로, 상기 제1 보간 데이터들을 가중시켜 가중 데이터들을 생성하는 단계; 및
    상기 가중 데이터들을 이용하여 촬영 객체에 대응되는 상기 자기 공명 영상을 획득하는 단계를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
  2. 제1항에 있어서, 상기 가중 데이터들을 생성하는 단계는
    상기 복수개의 주파수 대역들 별로, 상기 제1 보간 데이터들에 K 공간에서 가중치 적용(K-space weighting)을 수행하여 상기 가중 데이터들을 생성하는 단계를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
  3. 제2항에 있어서,
    상기 복수개의 주파수 대역들은 저 주파수 대역 및 고주파수 대역을 포함하며,
    상기 가중 데이터들을 생성하는 단계는
    상기 제1 보간 데이터들을 로우 패스 필터링하여 상기 저 주파수 대역에 대응되는 상기 가중 데이터들을 생성하는 단계; 및
    서로 다른 복수개의 가중치 값을 적용하여, 상기 제1 보간 데이터들을 하이 패스 필터링하여 상기 고 주파수 대역에 대응되는 상기 가중 데이터들을 생성하는 단계를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
  4. 제1항에 있어서, 상기 자기 공명 영상을 획득하는 단계는
    상기 복수개의 주파수 대역들 중 적어도 하나에서, 상기 가중 데이터들을 보간하여 제2 보간 데이터들을 생성하는 단계를 더 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
  5. 제1항에 있어서, 상기 자기 공명 영상을 획득하는 단계는
    상기 복수개의 주파수 대역들 중 적어도 하나에서, 상기 가중 데이터들을 보간하여 제2 보간 데이터들을 생성하는 단계;
    상기 가중 데이터들 중 적어도 하나 및 상기 제2 보간 데이터들 중 적어도 하나를 가중 조합(weighted combination)하여, 상기 복수개의 코일들에 각각 대응되는 복수개의 복원 영상을 획득하는 단계; 및
    상기 복수개의 복원 영상을 이용하여 상기 자기 공명 영상을 획득하는 단계를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
  6. 제5항에 있어서, 상기 제2 보간 데이터들을 생성하는 단계는
    상기 복수개의 주파수 대역들 중 적어도 하나에서, 상기 가중 데이터들을 회선 보간(convolution interpolation)하여 상기 제2 보간 데이터들을 생성하는 단계를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
  7. 제5항에 있어서, 상기 자기 공명 영상을 획득하는 단계는
    상기 복수개의 복원 영상들을 역 고속 퓨리에 변환(Inverse Fast Fourier Transform)하는 단계; 및
    상기 역 고속 퓨리에 변환된 상기 복수개의 복원 영상들에 근거하여 상기 자기 공명 영상을 획득하는 단계를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
  8. 제7항에 있어서, 상기 자기 공명 영상을 획득하는 단계는
    상기 역 고속 퓨리에 변환된 상기 복수개의 복원 영상을 제곱 합(sum of squares)하여, 상기 자기 공명 영상을 획득하는 단계를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
  9. 제1항에 있어서, 상기 제1 보간 데이터들을 생성하는 단계는
    상기 복수개의 데이터들을 회선 보간(convolution interpolation)하여 복수개의 제1 보간 데이터들을 생성하는 단계를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
  10. 복수개의 코일 각각에서 획득된 복수개의 데이터들을 보간하여 복수개의 제1 보간 데이터들을 생성하는 데이터 획득부;
    복수개의 주파수 대역들 별로, 상기 제1 보간 데이터들을 가중시켜 가중 데이터들을 생성하는 주파수 대역별 처리부; 및
    상기 가중 데이터들을 이용하여 촬영 객체에 대응되는 자기 공명 영상을 획득하는 영상 획득부를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 장치.
  11. 제10항에 있어서, 상기 주파수 대역별 처리부는
    상기 복수개의 주파수 대역들 별로, 상기 제1 보간 데이터들에 K 공간에서 가중치 적용(K-space weighting)을 수행하여 상기 가중 데이터들을 생성하는 것을 특징으로 하는 자기 공명 영상 생성 장치.
  12. 제10항에 있어서, 상기 주파수 대역별 처리부는
    상기 제1 보간 데이터들을 로우 패스 필터링하여 상기 가중 데이터들을 생성하는 적어도 하나의 로우 패스 필터; 및
    상기 제1 보간 데이터를 입력받고 소정 가중치를 적용하여 하이 패스 필터링하는 하이 패스 필터를 다수개 포함하는 하이 패스 필터 부 포함하며,
    상기 복수개의 하이 패스 필터들 각각은 서로 다른 가중치 값을 적용받는 것을 특징으로 하는 자기 공명 영상 생성 장치.
  13. 제10항에 있어서, 상기 데이터 획득부는
    상기 복수개의 데이터들을 수신하는 데이터 수신부; 및
    상기 복수개의 코일 각각에 대응되는 상기 복수개의 데이터들을 회선 보간(convolution interpolation)하여 상기 복수개의 제1 보간 데이터들을 생성하는 제1 보간부를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 장치.
  14. 제10항에 있어서, 상기 영상 획득부는
    상기 복수개의 주파수 대역들 중 적어도 하나에서, 상기 가중 데이터들을 보간하여 제2 보간 데이터들을 생성하는 제2 보간부;
    상기 가중 데이터들 중 적어도 하나 및 상기 제2 보간 데이터들 중 적어도 하나를 이용하여, 상기 복수개의 코일에 대응되는 복수개의 복원 영상을 획득하는 복원 영상 획득부; 및
    상기 복수개의 복원 영상들을 이용하여 상기 자기 공명 영상을 획득하는 자기 공명 영상 획득부를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 장치.
  15. 제14항에 있어서, 상기 자기 공명 영상 획득부는
    상기 복수개의 복원 영상들을 역 고속 퓨리에 변환(Inverse Fast Fourier Transform)하고, 상기 역 고속 퓨리에 변환된 상기 복수개의 복원 영상에 근거하여 상기 자기 공명 영상을 획득하는 것을 특징으로 하는 자기 공명 영상 생성 장치.
  16. 제15항에 있어서, 상기 자기 공명 영상 획득부는
    상기 역 고속 퓨리에 변환된 상기 복수개의 복원 영상을 제곱 합(sum of squares)하여, 상기 자기 공명 영상을 획득하는 것을 특징으로 하는 자기 공명 영상 생성 장치.
  17. 제14항에 있어서, 상기 복원 영상 획득부는
    상기 가중 데이터들 중 적어도 하나 및 상기 제2 보간 데이터들 중 적어도 하나를 가중 조합(weighted combination)하여, 상기 복수개의 코일 각각에 대응되는 상기 복수개의 복원 영상을 획득하는 것을 특징으로 하는 자기 공명 영상 생성 장치.
  18. 제12항에 있어서, 상기 영상 획득부는
    상기 복수개의 하이 패스 필터에서 각각 출력되는 상기 가중 데이터들을 회선 보간(convolution interpolation)하여, 상기 복수개의 하이 패스 필터에 각각 대응되는 복수개의 제2 보간 데이터들을 생성하는 제2 보간부;
    상기 복수개의 코일 각각 별로, 상기 로우 패스 필터에서 출력되는 가중 데이터 및 상기 제2 보간 데이터들을 가중 조합(weighted combination)하여 복원 영상을 생성하는 복원 영상 획득부; 및
    상기 복수개의 코일에 대응하는 복수개의 상기 복원 영상을 이용하여 상기 자기 공명 영상을 획득하는 자기 공명 영상 획득부를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 장치.
  19. 복수개의 코일을 포함하는 고주파 멀티 코일을 이용하여 자기 공명 영상을 획득하기 위한 방법에 있어서,
    상기 복수개의 코일에서 획득된 복수개의 데이터들을 보간하여 복수개의 제1 보간 데이터들을 생성하는 단계; 및
    제1 보간 데이터들에 대응하는 데이터들을 보간하여 제2 보간 데이터들을 생성하는 단계를 포함하는 것을 특징으로 하는 자기 공명 영상 생성 방법.
KR1020120049779A 2012-05-10 2012-05-10 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치 KR101330638B1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020120049779A KR101330638B1 (ko) 2012-05-10 2012-05-10 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치
US13/673,629 US9196062B2 (en) 2012-05-10 2012-11-09 Method and apparatus for generating magnetic resonance image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020120049779A KR101330638B1 (ko) 2012-05-10 2012-05-10 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치

Publications (1)

Publication Number Publication Date
KR101330638B1 true KR101330638B1 (ko) 2013-11-18

Family

ID=49548647

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020120049779A KR101330638B1 (ko) 2012-05-10 2012-05-10 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치

Country Status (2)

Country Link
US (1) US9196062B2 (ko)
KR (1) KR101330638B1 (ko)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017007279A1 (ko) * 2015-07-09 2017-01-12 성균관대학교산학협력단 자기 공명 영상 장치 및 방법
WO2017039163A1 (ko) * 2015-09-04 2017-03-09 삼성전자 주식회사 자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치
KR101806902B1 (ko) 2016-02-16 2017-12-12 삼성전자주식회사 자기 공명 영상 복원 방법 및 장치
KR102144320B1 (ko) * 2020-02-13 2020-08-13 주식회사 에어스 메디컬 자기 공명 영상 처리 장치 및 그 방법
KR102429284B1 (ko) * 2021-08-04 2022-08-04 주식회사 에어스메디컬 컴바인이 적용되는 자기 공명 영상 처리 장치 및 그 방법

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101310825B1 (ko) * 2012-05-10 2013-10-14 고려대학교 산학협력단 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치
KR101517606B1 (ko) * 2013-10-29 2015-05-04 한국과학기술원 K-공간의 저주파 정보를 기반으로 설계된 필터뱅크를 이용한 병렬영상기법
KR101667141B1 (ko) 2015-03-11 2016-10-25 한국과학기술원 소멸필터를 이용한 고속 mr 영상 복원 알고리듬 개발
WO2018137191A1 (en) * 2017-01-25 2018-08-02 Shanghai United Imaging Healthcare Co., Ltd. System and method for magnetic resonance imaging

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09238920A (ja) * 1996-03-12 1997-09-16 Toshiba Corp 磁気共鳴イメージング装置
JP2007510488A (ja) 2003-11-12 2007-04-26 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ パラレルmrイメージング方法

Family Cites Families (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2814700B2 (ja) 1990-06-04 1998-10-27 松下電器産業株式会社 自動焼物器
JP3064361B2 (ja) * 1990-09-07 2000-07-12 株式会社日立製作所 核磁気共鳴装置
US5910728A (en) 1996-11-12 1999-06-08 Beth Israel Deaconess Medical Center Simultaneous acquisition of spatial harmonics (SMASH): ultra-fast imaging with radiofrequency coil arrays
US6289232B1 (en) 1998-03-30 2001-09-11 Beth Israel Deaconess Medical Center, Inc. Coil array autocalibration MR imaging
US6208139B1 (en) * 1998-09-18 2001-03-27 General Electric Company System and method of phase sensitive MRI reconstruction using partial k-space data
US6198283B1 (en) * 1998-09-18 2001-03-06 Ge Medical Systems Global Technology Llc System and method of phase sensitive MRI reconstruction using partial k-space data and including a network
US6841998B1 (en) 2001-04-06 2005-01-11 Mark Griswold Magnetic resonance imaging method and apparatus employing partial parallel acquisition, wherein each coil produces a complete k-space datasheet
US7081749B2 (en) * 2004-08-02 2006-07-25 Albert Macovski Magnetic resonance imaging using demodulated k-space segments
US6833700B2 (en) * 2002-09-13 2004-12-21 General Electric Company Method and apparatus for reconstruction of images in parallel MRI systems
US20050203982A1 (en) * 2004-02-13 2005-09-15 Joseph Kolibal Method and apparatus for approximating, deconvolving and interpolating data using berstein functions
KR100611179B1 (ko) * 2004-06-23 2006-08-10 삼성전자주식회사 영상처리장치
US7298144B2 (en) * 2005-05-06 2007-11-20 The Board Of Trustee Of The Leland Stanford Junior University Homodyne reconstruction of water and fat images based on iterative decomposition of MRI signals
US7916144B2 (en) * 2005-07-13 2011-03-29 Siemens Medical Solutions Usa, Inc. High speed image reconstruction for k-space trajectory data using graphic processing unit (GPU)
WO2007085342A1 (en) * 2006-01-27 2007-08-02 Julius-Maximilians-Universität Würzburg Parallel mri involving density weighting or acquisition weighting
US7282917B1 (en) * 2006-03-30 2007-10-16 General Electric Company Method and apparatus of multi-coil MR imaging with hybrid space calibration
US20100142823A1 (en) * 2007-03-07 2010-06-10 Ze Wang 2d partially parallel imaging with k-space surrounding neighbors based data reconstruction
US20080292167A1 (en) * 2007-05-24 2008-11-27 Nick Todd Method and system for constrained reconstruction applied to magnetic resonance temperature mapping
US7881510B2 (en) * 2007-06-08 2011-02-01 Allegheny-Singer Research Institute Method and apparatus for forming an image with dynamic projective data
US7619410B2 (en) * 2007-06-08 2009-11-17 General Electric Company System and method for accelerated magnetic resonance parallel imaging
US8472688B2 (en) * 2008-04-17 2013-06-25 Wisconsin Alumni Research Foundation Method for image reconstruction employing sparsity-constrained iterative correction
CN102804207B (zh) * 2009-06-19 2016-08-03 微雷公司 用于执行断层图像获取和重构的系统和方法
US8686725B2 (en) * 2009-06-29 2014-04-01 General Electric Company System and apparatus for frequency translation of magnetic resonance (MR) signals
US8831318B2 (en) * 2009-07-06 2014-09-09 The Board Of Trustees Of The University Of Illinois Auto-calibrating parallel MRI technique with distortion-optimal image reconstruction
US8354844B2 (en) * 2009-12-23 2013-01-15 Universitaetsklinikum Freiburg Method for data acquisition acceleration in magnetic resonance imaging (MRI) with N-dimensional spatial encoding using two or more receiver coil arrays and non-linear phase distributions
US8638096B2 (en) * 2010-10-19 2014-01-28 The Board Of Trustees Of The Leland Stanford Junior University Method of autocalibrating parallel imaging interpolation from arbitrary K-space sampling with noise correlations weighted to reduce noise of reconstructed images
DE102010061977B4 (de) * 2010-11-25 2013-10-17 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Ergänzung mit partiell parallelen Aufnahmeverfahren aufgenommen Magnetresonanz-Daten
US9390521B2 (en) * 2010-12-22 2016-07-12 Koninklijke Philips N.V. Rapid parallel reconstruction for arbitrary k-space trajectories
US9588207B2 (en) * 2011-10-06 2017-03-07 National Institutes of Health (NIH), U.S. Dept. of Health and Human Services (DHHS), The United States of America NIH Division of Extramural Inventions and Technology Resources (DEITR) System for reconstructing MRI images acquired in parallel
EP2773985A4 (en) * 2011-11-06 2015-05-06 Mayo Foundation PROCESS FOR CALIBRATION-FREE LOCAL RANK LOW ACTIVATION OF THE RECONSTRUCTION OF MAGNETIC RESONANCE IMAGES
US9594141B2 (en) * 2012-02-15 2017-03-14 National Institutes of Health (NIH), U.S. Dept. of Health and Human Services (DHHS), The United States of America NIH Division of Extramural Inventions and Technology Resources (DEITR) System for accelerated magnetic resonance imaging using parallel coils
KR101310825B1 (ko) * 2012-05-10 2013-10-14 고려대학교 산학협력단 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치
US9097780B2 (en) * 2012-11-09 2015-08-04 Siemens Aktiengesellschaft Multi-stage magnetic resonance reconstruction for parallel imaging applications
US9329250B2 (en) * 2013-03-14 2016-05-03 Wisconsin Alumni Resesarch Foundation System and method for combined chemical species separation and high resolution R2* mapping with magnetic resonance imaging

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09238920A (ja) * 1996-03-12 1997-09-16 Toshiba Corp 磁気共鳴イメージング装置
JP2007510488A (ja) 2003-11-12 2007-04-26 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ パラレルmrイメージング方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017007279A1 (ko) * 2015-07-09 2017-01-12 성균관대학교산학협력단 자기 공명 영상 장치 및 방법
WO2017039163A1 (ko) * 2015-09-04 2017-03-09 삼성전자 주식회사 자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치
KR20170028829A (ko) * 2015-09-04 2017-03-14 삼성전자주식회사 자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치
KR102527812B1 (ko) * 2015-09-04 2023-05-02 삼성전자주식회사 자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치
KR101806902B1 (ko) 2016-02-16 2017-12-12 삼성전자주식회사 자기 공명 영상 복원 방법 및 장치
KR102144320B1 (ko) * 2020-02-13 2020-08-13 주식회사 에어스 메디컬 자기 공명 영상 처리 장치 및 그 방법
WO2021162176A1 (ko) * 2020-02-13 2021-08-19 주식회사 에어스 메디컬 자기 공명 영상 처리 장치 및 그 방법
KR20210103363A (ko) 2020-02-13 2021-08-23 주식회사 에어스 메디컬 자기 공명 영상 처리 장치 및 그 방법
KR102588863B1 (ko) 2020-02-13 2023-10-16 주식회사 에어스메디컬 자기 공명 영상 처리 장치 및 그 방법
KR102429284B1 (ko) * 2021-08-04 2022-08-04 주식회사 에어스메디컬 컴바인이 적용되는 자기 공명 영상 처리 장치 및 그 방법

Also Published As

Publication number Publication date
US20130301891A1 (en) 2013-11-14
US9196062B2 (en) 2015-11-24

Similar Documents

Publication Publication Date Title
KR101330638B1 (ko) 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치
KR101310825B1 (ko) 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치
US9482732B2 (en) MRI reconstruction with motion-dependent regularization
JP6073627B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
JP6018401B2 (ja) 拡散強調エコープラナー撮像法において高次渦電流に誘発された歪みを予測補正するためのシステムおよび方法
JP5349673B2 (ja) データ補正装置
RU2672151C2 (ru) Коррекция посторонних эхосигналов epi
JP5138043B2 (ja) 磁気共鳴イメージング装置
JP6367056B2 (ja) 磁気共鳴イメージング装置
JP2014508023A (ja) 画像劣化なしのより速いチャネル毎の再構成のmri方法
JP2007244848A (ja) データ補正装置、データ補正方法、磁気共鳴イメージング装置およびx線ct装置
JPH025927A (ja) 磁気共鳴画像化方法および装置
JP2016022302A (ja) 磁気共鳴信号処理方法、磁気共鳴信号処理装置及び磁気共鳴装置並びにプログラム
US10175328B2 (en) System and method for reconstructing ghost-free images from data acquired using simultaneous multislice magnetic resonance imaging
JP2019535434A (ja) ディクソンタイプの水/脂肪分離を用いたmrイメージング
JP6986659B2 (ja) 磁気共鳴イメージング装置、ナイキストゴースト補正方法、及び、ナイキストゴースト補正用プログラム
JP6820876B2 (ja) 磁気共鳴撮影装置、感度分布算出方法、及び、感度分布算出プログラム
US10228434B2 (en) Multi-shot echo planar imaging using reordered segments and RF excitation pulse phase and slice profiles matched across interleaves
JP2017136170A (ja) 磁気共鳴イメージング装置
JP5819675B2 (ja) 磁気共鳴イメージング装置、およびプログラム
JP3983792B2 (ja) 核磁気共鳴撮影装置
JP4675936B2 (ja) 核磁気共鳴撮影装置
EP2392935A1 (en) EPI distortion correction using non-phase encoded reference echoes
JP5566587B2 (ja) 磁気共鳴イメージング装置およびリンギング低減方法
JP3889770B2 (ja) 磁気共鳴イメージング装置

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: 20161028

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20171030

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20181030

Year of fee payment: 6