KR101974199B1 - 가변절편 자기공명영상 데이터 획득방법 - Google Patents
가변절편 자기공명영상 데이터 획득방법 Download PDFInfo
- Publication number
- KR101974199B1 KR101974199B1 KR1020170035099A KR20170035099A KR101974199B1 KR 101974199 B1 KR101974199 B1 KR 101974199B1 KR 1020170035099 A KR1020170035099 A KR 1020170035099A KR 20170035099 A KR20170035099 A KR 20170035099A KR 101974199 B1 KR101974199 B1 KR 101974199B1
- Authority
- KR
- South Korea
- Prior art keywords
- echo
- frequency
- image
- low
- lines
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
-
- 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
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/5635—Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/446—Multifrequency selective RF pulses, e.g. multinuclear acquisition mode
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56536—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to magnetic susceptibility variations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
- G01R33/482—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
- G01R33/4822—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory in three dimensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- Vascular Medicine (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
본 발명은 자기공명영상의 동맥조영도와 정맥조영도를 얻기 위한 데이터 획득방법에 관한 것으로서, 한 개 이상의 에코를 이용하며, 절편의 수에 따라 최적화되는 동맥조영도와 정맥조영도를 한 번의 촬영으로 동시에 얻거나 동맥조영도의 절편 경계부분 연결성을 개선하는 방법을 제공한다.
Description
본 발명은 자기공명영상을 위한 데이터 획득방법에 관한 것으로서, 한 개 이상의 에코를 이용하며, 한 번의 촬영으로 동맥조영도와 정맥조영도를 한 번에 얻거나 동맥조영도의 질을 개선할 수 있는 기술에 관한 것이다.
자기공명영상은 영상 기술 중 하나로 핵자기공명 원리를 사용한다. 자기장을 발생하는 자기공명 촬영 장치에 인체를 넣고 고주파를 발생시키면 신체의 수소 원자의 전자가 공명하게 된다. 이때 나오는 신호의 차이를 측정하고 컴퓨터를 통해 재구성하여 영상화시키면 자기공명 영상이 된다.
자기공명영상을 얻기 위하여, 동맥조영술(Time-of-flight MR angiography) 및/또는 정맥조영술(Blood oxygenation level dependent MR venogram)을 이용할 수 있다. 이에 관한 선공개 특허문헌으로서, "Correction of artifacts in time-of-flight MR angiography"라는 발명의 명칭으로 2010-05-27에 공개된 미국특허 공개번호 US20100128952A1의 특허문헌과, "Method to evaluate patients for thoracic outlet syndrome"이라는 발명의 명칭으로 2015-02-24에 등록된 미국특허 등록번호 US8965481의 특허문헌이 있다.
한편, CODEA(compatible dual-echo arteriovenography) 기법은 상기 동맥조영술과 정매조영술의 이미지를 동시에 얻을 수 있도록 해주는 기술이다. 이에 관한 선공개 특허문헌으로서, "Echo-specific K-space reordering approach to compatible dual-echo arteriovenography"라는 발명의 명칭으로 2011-09-01에 공개된 미국특허 공개번호 US20110213237A1의 특허문헌이 있다.
일반적으로 동맥조영술은 여러 개의 절편(슬랩,slab)으로 나누어 촬영할 때 혈관이 더 잘 검출되며, 정맥조영술은 하나의 절편으로 촬영해야 SNR이 좋아진다.
기존의 COEA 기법은 동맥조영도과 정맥조영도를 동시에 얻을 수 있도록 해주지만, 절편의 수를 최적화하는데 문제가 있다. 또한, 기존의 동맥조영술을 이용하는 경우 여러 개의 절편으로 나누어 촬영할 때 절편 경계부분에서 인공물이 발생할 가능성이 있다.
상술한 문제를 해결하기 위해 본 발명에서는 절편의 수에 따라 최적화되는 동맥조영도와 정맥조영도를 한 번의 촬영으로 동시에 얻거나 동맥조영도의 절편 경계부분 연결성을 개선하는 데이터 획득방법 및 영상처리방법을 제공하고자 한다.
본 발명의 일 관점에 따라 한 개 이상의 에코를 이용한 이미지 획득방법(=데이터 획득방법)을 제공할 수 있다. 이때, 상기 데이터 획득방법은, 복수 절편(11, 12) 중 제1절편(11)만 여기하는 제1RF(110)를 복수 회 여기하는 제1프로세스를 포함하며, 상기 제1RF(110)를 각각 여기할 때 마다, 제1에코(111)에서는 미리 결정된 순서에 따라 저주파성분(211)을 얻고, 제2에코(112)에서는 미리 결정된 순서에 따라 고주파성분(221)을 얻을 수 있다. 상기 제1프로세스 이후, 복수 절편 중 제2절편(12)만 여기하는 제2RF(120)를 복수 회 여기하는 제2프로세스를 포함하며, 상기 제2RF(120)를 각각 여기할 때 마다, 제1에코(121)에서는 미리 결정된 순서에 따라 저주파성분(212)을 얻고, 제2에코(122)에서는 미리 결정된 순서에 따라 고주파성분(222)을 얻을 수 있다. 상기 제2프로세스 이후, 복수 절편 전체(=전체절편)(10)를 한 번에 여기하는 제3RF(130)를 복수 회 여기하는 제3프로세스를 포함하며, 상기 제3RF(130)를 각각 여기할 때마다, 제1에코(131)에서는 미리 결정된 순서에 따라 고주파성분(213)을 얻고, 제2에코(132)에서는 미리 결정된 순서에 따라 저주파성분(223)을 얻을 수 있다.
이때, 상기 제1프로세스의 상기 제1에코들에서 획득한 제1 저주파 위상인코딩라인들(211)을 이미지 영역으로 변환하여 제1 저주파 서브영상(511)을 생성하는 단계; 상기 제2프로세스의 상기 제1에코들에서 획득한 제2 저주파 위상인코딩라인들(212)을 이미지 영역으로 변환하여 제2 저주파 서브영상(512)을 생성하는 단계; 상기 제1 저주파 서브영상(511)과 상기 제2 저주파 서브영상(512)을 결합하여 제3 저주파 서브영상(611)을 생성하는 단계; 상기 제3 저주파 서브영상(611)을 주파수 도메인으로 변환하여 제1 저주파 데이터 세트(311)를 획득하는 단계; 상기 제1 저주파 데이터 세트(311)와 상기 제3프로세스의 상기 제1에코들에서 획득한 제3 고주파 위상인코딩라인들(312)을 주파수 영역에서 서로 결합하여 제1 K-스페이스데이터(210)를 획득하는 단계;를 포함할 수 있다.
이때, 상기 제1프로세스의 상기 제2에코들에서 획득한 제1 고주파 위상인코딩라인들(221)을 이미지 영역으로 변환하여 제1 고주파 서브영상(721)을 생성하는 단계; 상기 제2프로세스의 상기 제2에코들에서 획득한 제2 고주파 위상인코딩라인들(222)을 이미지 영역으로 변환하여 제2 고주파 서브영상(722)을 생성하는 단계; 상기 제1 고주파 서브영상(721)과 상기 제2 고주파 서브영상(722)을 결합하여 제3 고주파 서브영상(421)을 생성하는 단계; 상기 제3 고주파 서브영상(421)을 주파수 도메인으로 변환하여 제2 고주파 데이터 세트(321)를 획득하는 단계; 상기 제2 고주파 데이터 세트와 상기 제3프로세스의 상기 제2에코들에서 획득한 제3 저주파 위상인코딩라인들(322)을 주파수 영역에서 서로 결합하여 제2 K-스페이스데이터(220)를 획득하는 단계;를 포함할 수 있다.
또는, 상기 제1프로세스의 상기 제1에코들에서 획득한 제1 저주파 위상인코딩라인들(211), 상기 제2프로세스의 상기 제1에코들에서 획득한 제2 저주파 위상인코딩라인들(212), 및 상기 제3프로세스의 상기 제1에코들에서 획득한 제3 고주파 위상인코딩라인들(213)을 각각 이미지 영역으로 변환하여 이미지 영역에서 서로 결합하는 단계;를 포함할 수 있다.
그리고 상기 제1프로세스의 상기 제2에코들에서 획득한 제1 고주파 위상인코딩라인들(221), 상기 제2프로세스의 상기 제2에코들에서 획득한 제2 고주파 위상인코딩라인들(222), 및 상기 제3프로세스의 상기 제2에코들에서 획득한 제3 저주파 위상인코딩라인들(223)을 각각 이미지 영역으로 변환하여 이미지 영역에서 서로 결합하는 단계;를 포함할 수 있다.
이때, 상기 제1에코에서의 촬영조건과 상기 제2에코에서의 촬영조건은 서로 다를 수 있다.
이때, 상기 제1에코는 상기 제2에코보다 선행하며, 제1에코 시에는 Time of Flight MR angiogram 을 얻고, 제2에코 시에는 Susceptiblitiy weighted imaging 을 얻을 수 있다.
이때, 상기 제1에코 시의 절편의 개수는 제2에코 시의 절편의 개수보다 많을 수 있다.
이때, 상기 제1에코 시에는 동맥 이미지를 획득하고 상기 제2에코 시에는 정맥 이미지를 획득하도록 되어 있을 수 있다.
본 발명의 다른 관점에 따라 한 개 이상의 에코를 이용한 데이터획득방법을 제공할 수 있다. 이때, 상기 데이터획득방법은, 제1절편(11)에 대하여 제1RF 신호(110)를 여기하여 발생한 제1에코(111) 및 제2에코(112) 각각에 대하여 K-스페이스의 라인들 중 저주파 라인들(211) 및 고주파 라인들(221)을 획득하고, 제2절편(12)에 대하여 제2RF 신호(120)를 여기하여 발생한 제3에코(121) 및 제4에코(122) 각각에 대하여 K-스페이스의 라인들 중 저주파 라인들(212) 및 고주파 라인들(222)을 획득하고, 상기 제1절편(11) 및 상기 제2절편(12)을 포함하는 전체절편(10)에 대하여 제3RF 신호(130)를 여기하여 발생한 제5에코(131)에서는 고주파 라인들(213)을 획득하고 제6에코(132)에서는 저주파 라인들(223)을 획득하는 것을 특징으로 할 수 있다.
이때, 상기 제1에코가 제2에코에 대해 선행하거나, 또는 상기 제2에코가 제1에코에 대해 선행할 수 있다.
본 발명의 또 다른 관점에 따라 한 개 이상의 에코를 이용한 데이터획득방법을 제공할 수 있다. 상기 데이터획득방법은, 제1절편(11) 및 제2절편(12) 각각에 대해 제1에코(111, 121)에서는 제1오더링에 따른 일부 라인들(211, 212)을 획득하고 제2에코(112, 122)에서는 제2오더링에 따른 일부 라인들(221, 222)을 획득하고, 전체절편(10)에 대하여 제1에코(131)에서는 상기 제1오더링에 따른 잔여 라인들 중 적어도 일부(213)를 획득하고, 제2에코(132)에서는 상기 제2오더링에 따른 잔여 라인들 중 적어도 일부(223)를 획득하는 것을 특징으로 할 수 있다.
본 발명의 또 다른 관점에 따른 데이터 획득방법은, 제1절편(11) 및 제2절편(12) 각각에 대해 제1오더링에 따른 일부 라인들(211, 212)을 획득하고, 전체절편(10)에 대하여 상기 제1오더링에 따른 잔여 라인들(213)을 획득하고, 상기 일부 라인들과 잔여 라인들에 대한 정보를 서로 결합하는 단계를 포함할 수 있다.
본 발명의 또 다른 관점에 따라 한 개 이상의 에코를 이용한 데이터획득방법을 제공할 수 있다. 상기 데이터획득방법은, 전체절편(10)을 분할하여 정의한 복수 개의 절편들(11, 12) 각각에 대해 제1에코(111, 121)에서는 저주파 라인(211, 212)을 획득하고 제2에코(112, 122)에서는 고주파 라인(221, 222)을 획득하며, 상기 전체절편(10)에 대하여 상기 제1에코(131)에서는 고주파 라인(213)을 획득하고 제2에코(132)에서는 저주파 라인(223)을 획득하는 것을 특징으로 할 수 있다.
본 발명에 따르면, 절편의 수에 따라 최적화되는 동맥조영도와 정맥조영도를 한 번의 촬영으로 동시에 얻거나 동맥조영도의 절편 경계부분 연결성을 개선할 수 있다.
도 1은 본 발명의 제1 실시예에 따라 절편 개수에 따른 영상의 질을 설명하기 위한 도면이다.
도 2는 본 발명의 제1 실시예에 따른 RF 펄스 파형 및 에코 파형을 시간의 흐름에 따라 나타낸 것이다.
도 3은 본 발명의 제1 실시예에 따른 K-스페이스(K-space)에서의 데이터 획득 모식도를 나타낸 것이다.
도 4는 본 발명의 제1 실시예에 따른 K-스페이스에서의 데이터 구성을 나타낸 것이다.
도 5는 본 발명의 제3 실시예에 따라 K-스페이스데이터를 획득하기 위한 일 과정을 설명하기 위한 도면이다.
도 6은 본 발명의 제3 실시예에 따른 K-스페이스(K-space)에서의 데이터 획득 모식도를 나타낸 것이다.
도 7은 본 발명의 제3 실시예에 따른 K-스페이스에서의 데이터 구성을 나타낸 것이다.
도 8은 본 발명의 제5 실시예에 따라 동맥조영 이미지를 획득하기 위한 영상처리방법을 설명하기 위한 도면이다.
도 9는 본 발명의 제5 실시예에 따라 정맥조영 이미지를 획득하기 위한 영상처리방법을 설명하기 위한 도면이다.
도 10은 본 발명의 제6 실시예에 따라 동맥조영 이미지 및 정맥조영 이미지를 획득하기 위한 영상처리방법을 설명하기 위한 도면이다.
도 11a는 관찰대상을 종래의 기술에 따라 3개의 절편으로 나누어 관찰하였을 때에 얻은 동맥조영도이다.
도 11b는 관찰대상을 본 발명의 일 실시예에 따라 3개의 절편으로 나누어 관찰하였을 때에 얻은 동맥조영도이다.
도 2는 본 발명의 제1 실시예에 따른 RF 펄스 파형 및 에코 파형을 시간의 흐름에 따라 나타낸 것이다.
도 3은 본 발명의 제1 실시예에 따른 K-스페이스(K-space)에서의 데이터 획득 모식도를 나타낸 것이다.
도 4는 본 발명의 제1 실시예에 따른 K-스페이스에서의 데이터 구성을 나타낸 것이다.
도 5는 본 발명의 제3 실시예에 따라 K-스페이스데이터를 획득하기 위한 일 과정을 설명하기 위한 도면이다.
도 6은 본 발명의 제3 실시예에 따른 K-스페이스(K-space)에서의 데이터 획득 모식도를 나타낸 것이다.
도 7은 본 발명의 제3 실시예에 따른 K-스페이스에서의 데이터 구성을 나타낸 것이다.
도 8은 본 발명의 제5 실시예에 따라 동맥조영 이미지를 획득하기 위한 영상처리방법을 설명하기 위한 도면이다.
도 9는 본 발명의 제5 실시예에 따라 정맥조영 이미지를 획득하기 위한 영상처리방법을 설명하기 위한 도면이다.
도 10은 본 발명의 제6 실시예에 따라 동맥조영 이미지 및 정맥조영 이미지를 획득하기 위한 영상처리방법을 설명하기 위한 도면이다.
도 11a는 관찰대상을 종래의 기술에 따라 3개의 절편으로 나누어 관찰하였을 때에 얻은 동맥조영도이다.
도 11b는 관찰대상을 본 발명의 일 실시예에 따라 3개의 절편으로 나누어 관찰하였을 때에 얻은 동맥조영도이다.
이하, 본 발명의 실시예를 첨부한 도면을 참고하여 설명한다. 그러나 본 발명은 본 명세서에서 설명하는 실시예에 한정되지 않으며 여러 가지 다른 형태로 구현될 수 있다. 본 명세서에서 사용되는 용어는 실시예의 이해를 돕기 위한 것이며, 본 발명의 범위를 한정하고자 의도된 것이 아니다. 또한, 이하에서 사용되는 단수 형태들은 문구들이 이와 명백히 반대의 의미를 나타내지 않는 한 복수 형태들도 포함한다.
<제1실시예 - K-스페이스데이터의 획득>
도 1은 본 발명의 제1 실시예에 따라 절편 개수에 따른 영상의 질을 설명하기 위한 도면이다.
전체절편(즉, 하나의 절편)(10)으로 촬영할 경우 동맥조영술(Time-of-flight MR angiography, TOF MRA)의 이미지(21)에서는 혈관이 잘 검출되진 않지만, 정맥조영술(Blood oxygenation level dependent MR venogram, BOLD MRV)의 이미지(22)에서는 SNR이 좋아질 수 있다.
반면, 전체절편(10)을 여러 개의 절편(11, 12)으로 나누어서 촬영하는 경우 동맥조영술(TOF MRA)의 이미지(23)에서는 혈관이 더 잘 검출되나, 정맥조영술(BOLD MRV)의 이미지(24)에서는 상기 SNR이 전체절편(10)으로 촬영했을 때 보다 좋지 않다.
기존의 compatible dual-echo arteriovenography (CODEA) 기법은 동맥조영술(TOF MRA)과 정맥조영술(BOLD MRV)을 동시에 얻을 수 있도록 해주는 기술이지만 절편의 수는 최적화하지 못한다.
제1실시예는 상기 절편의 수에 따라 다르게 최적화되는 동맥조영도과 정맥조영도를 한 번의 촬영으로 동시에 얻게 해줄 수 있다.
한 세트의 동맥조영술(TOF MRA) 이미지(21)와 정맥조영술(BOLD MRV) 이미지(22)를 얻기 위하여, MRI 데이터 획득대상(10)을 N개(ex: N=2)의 절편(slab)(ex: 11, 12)으로 나누어 처리한다. 즉, 상기 MRI 데이터 획득대상(10)을 '전체절편'이라고 명명하였을 때에, 상기 전체절편(10)은, 서로 겹치지 않는 제1 절편(11), 및 제2 절편(12)들로 이루어질 수 있다.
도 2는 본 발명의 제1 실시예에 따른 RF 펄스 파형, 위상인코딩라인을 얻기 위한 방식에 대한 파형, 및 에코 파형을 시간의 흐름에 따라 나타낸 것이다.
일반적으로 복수 개의 절편으로 분할되는 경우, 각 절편 당 RF 펄스를 복수 회 발생시킬 수 있다. 본 발명의 도 2에서는 설명의 편의를 위하여, 각 절편 당 복수 회의 RF 펄스 중 1회의 RF 펄스를 발생하는 것으로 나타내었다.
본 발명의 제1실시예에 따른 K-스페이스데이터 획득방법에서는, 분할될 절편의 개수가 N(ex: N=2)인 경우, RF 펄스를 N+1회 또는 그 이상 발생시킬 수 있다. 즉, 제1 RF 펄스 여기(excitation)(110), 제2 RF 펄스 여기(120), 및 제3 RF 펄스 여기(130)가 순차적으로 수행될 수 있다.
이때, 기존에는 각 RF 펄스 여기(예컨대, 110)마다 한 개의 에코(예컨대, 111)가 발생하도록 되어 있을 수 있다. 일반적으로 정교한 이미지를 얻기 위하여 다수의 RF 펄스 여기를 발생시키도록 되어 있을 수 있다. 그러나 이러한 경우 시간이 오래 걸린다는 문제가 있다. 따라서 본 발명에서는 각 RF 펄스 여기(예컨대, 110) 당 두 개의 에코(예컨대, 111, 112)를 발생시킬 수 있다.
이때, 위상인코딩라인을 얻기 위한 방식에 대한 파형(예컨대, GPE1 및 GPE2의 파형)에 따라 각 RF 펄스 여기에 대한 두 개의 에코의 파형이 결정될 수 있다. 예컨대, 제1방식(51)에 따라 제1에코(111)의 파형이 결정될 수 있으며, 제2방식(52)에 따라 제2에코(112)의 파형이 결정될 수 있다. 이때, RF 펄스 여기가 발생한 시점부터 시간이 흐를수록 그 신호의 세기는 감소한다. 따라서 제2에코(112)에서 많은 데이터를 얻기 위하여 위상인코딩라인을 얻기 위한 방식은 예컨대 센트릭 PE 오더 방식과 같을 수 있다.
도 3은 본 발명의 제1 실시예에 따른 K-스페이스(K-space)에서의 데이터 획득 모식도를 나타낸 것이다.
이때, K-스페이스는 3D 형태의 공간이며, 각 K-스페이스의 x 축은 시간, y 축은 주파수, 및 z 축은 절편(예컨대, 11)의 높이를 의미할 수 있다.
도 4는 본 발명의 제1 실시예에 따른 K-스페이스에서의 데이터 구성을 나타낸 것이다. 도 4의 (a) 제1RF 및 제2RF에 대한 획득된 K-스페이스를 나타낸 것이고, 도 4의 (b)는 제3RF에 대한 획득된 K-스페이스를 나타낸 것이며, 도 4의 (c)는 동맥조영술에 대하여 각 에코별로 조합한 것을 나타낸 것이고, 도 4의 (d)는 정맥조영술에 대하여 각 에코별로 조합한 것을 나타낸 것이다.
이하, 도 1 내지 도 4를 함께 참조하여 설명한다.
본 발명의 제1실시예에 따른 K-스페이스데이터 획득방법은 아래의 단계들을 포함할 수 있다.
단계(S110):
제1 RF 펄스 여기(110) 내지 제2 RF 펄스 여기(120) 각각에 대하여, 상기 제1 절편(11) 내지 제2 절편(12)에 대한 K-스페이스데이터를 획득할 수 있다.
이때, 제1 RF 펄스 여기(110) 내지 제2 RF 펄스 여기(120) 각각에 대하여, 복수 개의 에코([111, 112] or [121, 122])에서 서로 다른 주파수 영역에 대응하는 위상인코딩라인들을 획득할 수 있다.
즉, 제1 RF 펄스 여기(110) 내지 제2 RF 펄스 여기(120) 각각에 대하여, 상기 복수 개의 에코들 중 제1에코(111 or 121)에 대하여는 제1 저주파 위상인코딩라인들(211) 내지 제2 저주파 위상인코딩라인들(212)을 획득하고, 제2에코(112 or 122)에 대하여는 제1 고주파 위상인코딩라인들(221) 내지 제2 고주파 위상인코딩라인들(222)을 획득할 수 있다. 이때, 바람직하게는 제1에코(111 or 121)는 제2에코(112 or 122)에 대하여 시간 축에서 선행하는 에코일 수 있다.
여기서, 상기 제1 RF 펄스 여기(110) 내지 제2 RF 펄스 여기(120) 각각에 대하여, 즉, 상기 제1 절편(11) 내지 제2 절편(12)에 대하여, 상기 제1에코(111 및 121)에 대하여 얻은 저주파 위상인코딩라인들의 데이터를 제1 저주파 데이터 세트(311)라고 명명할 수 있다. 즉 상기 제1 저주파 데이터 세트(311)는 상기 제1 저주파 위상인코딩라인들(211) 내지 제2 저주파 위상인코딩라인(212)들로 구성될 수 있다.
그리고 상기 제1 RF 펄스 여기(110) 내지 제2 RF 펄스 여기(120) 각각에 대하여, 즉, 상기 제1 절편(11) 내지 제2 절편(12)에 대하여, 상기 제2에코(112 및 122)에 대하여 얻은 고주파 위상인코딩라인들의 데이터를 제2 고주파 데이터 세트(321)라고 명명할 수 있다. 즉 상기 제2 고주파 데이터 세트(321)는 상기 제1 고주파 위상인코딩라인들(221) 내지 제N 고주파 위상인코딩라인들(222)로 구성될 수 있다.
단계(S120):
제3 RF 펄스 여기(130)에 대하여, 상기 전체절편(10)에 대한 K-스페이스데이터를 획득할 수 있다.
이때, 제3 RF 펄스 여기(130)에 대하여, 복수 개의 에코(131 및 132)에서 서로 다른 주파수 영역에 대응하는 위상인코딩라인들을 획득할 수 있다. 즉, 제3 RF 펄스 여기(130)에 대하여, 상기 복수 개의 에코들(131 및 132) 중 제1에코(131)에 대하여는 제3 고주파 위상인코딩라인들(213)을 획득하고, 제2에코(132)에 대하여 제3 저주파 위상인코딩라인들(223)을 획득할 수 있다. 이때, 바람직하게는 제1에코(131)는 제2에코(132)에 대하여 시간 축에서 선행하는 에코일 수 있다.
여기서, 상기 제3 RF 펄스 여기(130)에 대하여, 즉, 상기 전체절편(10)에 대하여, 상기 제1에코(131)에 대하여 얻은 제3 고주파 위상인코딩라인들(213)의 데이터를 제1 고주파 데이터 세트(312)라고 명명할 수 있다. 그리고 상기 제3 RF 펄스 여기(130)에 대하여, 즉, 상기 전체절편(10)에 대하여, 상기 제2에코(132)에 대하여 얻은 제3 저주파 위상인코딩라인들(223)의 데이터를 제2 저주파 데이터 세트(322)라고 명명할 수 있다.
단계(S130):
상기 제1 저주파 데이터 세트(311) 및 상기 제1 고주파 데이터 세트(312)를 서로 조합하여 제1 K-스페이스데이터(210)를 생성하고, 상기 제2 고주파 데이터 세트(321) 및 상기 제2 저주파 데이터 세트(322)를 서로 조합하여 제2 K-스페이스데이터(220)를 생성할 수 있다. 이렇게 생성된 상기 제1 K-스페이스데이터(210)와 제2 K-스페이스데이터(220)는 각각, 서로 다른 영상을 생성하는 데 사용될 수 있다.
<제2실시예 - 획득한 K-스페이스데이터를 이용한 영상처리방법>
본 발명의 제2실시예에 따른 영상처리방법은, 제1실시예에서 상술한 방법으로 획득한 K-스페이스데이터들을 조합하여 영상을 생성하는 기술에 관한 것이다.
본 발명의 제2실시예에 따른 영상처리방법에 의해 생성하는 제1타입영상은, 아래 방법 1 및 방법 2 중 한 가지 방법에 의해 생성될 수 있다.
방법 1: 상술한 제1 K-스페이스데이터를 IFFT 하여 상기 제1타입영상을 획득하는 방법.
방법 2: (1) 상기 제1 저주파 위상인코딩라인들로 이루어진 제1 저주파 서브 K-스페이스데이터 내지 상기 제N 저주파 위상인코딩라인들로 이루어진 제N 저주파 서브 K-스페이스데이터를 각각 IFFT하여 제1 저주파 서브영상 내지 제N 저주파 서브영상을 생성하는 단계, (2) 상기 제N+1 고주파 위상인코딩라인들로 이루어진 제N+1 고주파 서브 K-스페이스데이터를 IFFT하여 제N+1 고주파 서브영상을 생성하는 단계, (3) 상기 제1 저주파 서브영상 내지 상기 제N 저주파 서브영상과 상기 제N+1 고주파 서브영상을 서로 결합하여 상기 제1타입영상을 생성하는 단계를 포함하는 방법.
본 발명의 다른 실시예에 따른 영상처리방법에 의해 생성하는 제2타입영상은, 아래 방법 3 및 방법 4 중 한 가지 방법에 의해 생성될 수 있다.
방법 3: 상술한 제2 K-스페이스데이터를 IFFT 하여 상기 제2타입영상을 획득하는 방법.
방법 4: (1) 상기 제1 고주파 위상인코딩라인들로 이루어진 제1 고주파 서브 K-스페이스데이터 내지 상기 제N 고주파 위상인코딩라인들로 이루어진 제N 고주파 서브 K-스페이스데이터를 각각 IFFT하여 제1 고주파 서브영상 내지 제N 고주파 서브영상을 생성하는 단계, (2) 상기 제N+1 저주파 위상인코딩라인들로 이루어진 제N+1 저주파 서브 K-스페이스데이터를 IFFT하여 제N+1 저주파 서브영상을 생성하는 단계, (3) 상기 제1 고주파 서브영상 내지 상기 제N 고주파 서브영상 및 상기 제N+1 저주파 서브영상을 서로 결합하여 상기 제2타입영상을 생성하는 단계를 포함하는 방법.
<제3 실시예 - K-스페이스데이터의 획득 - 일반화>
도 5는 본 발명의 제3 실시예에 따라 K-스페이스데이터를 획득하기 위한 일 과정을 설명하기 위한 도면이다.
도 6은 본 발명의 제3 실시예에 따른 K-스페이스(K-space)에서의 데이터 획득 모식도를 나타낸 것이다.
이때, K-스페이스는 3D 형태의 공간이며, 각 K-스페이스의 x 축은 시간, y 축은 주파수, 및 z 축은 절편(예컨대, 11)의 높이를 의미할 수 있다.
도 7은 본 발명의 제3 실시예에 따른 K-스페이스에서의 데이터 구성을 나타낸 것이다. 도 7의 (a) 제1RF, 제2RF, 및 제3RF에 대한 획득된 K-스페이스를 나타낸 것이고, 도 7의 (b)는 제4RF에 대한 획득된 K-스페이스를 나타낸 것이며, 도 7의 (c)는 동맥조영술에 대하여 각 에코별로 조합한 것을 나타낸 것이고, 도 7의 (d)는 정맥조영술에 대하여 각 에코별로 조합한 것을 나타낸 것이다.
이하, 도 5 내지 도 7을 함께 참조하여 설명한다.
한 세트의 TOF MRA 이미지(21)와 BOLD MRV 이미지(22)를 얻기 위하여, MRI 데이터 획득대상(10)을 N개의 절편(slab)(ex: 11, 12, 13)으로 나누어 처리한다. 즉, 상기 MRI 데이터 획득대상(10)을 '전체절편'이라고 명명하였을 때에, 상기 전체절편(10)은, 서로 겹치지 않는 제1 절편(11), 제2 절편(12), 제3 절편(13), 내지 제N 절편들로 이루어질 수 있다.
본 발명의 제3 실시예에 따른 K-스페이스데이터 획득방법에서는, 분할될 절편의 개수가 N(ex: N=2 or 3)인 경우, RF 펄스를 N+1회 또는 그 이상 발생시킬 수 있다. 즉, 제1 RF 펄스 여기(excitation), 제2 RF 펄스 여기, 제3 RF 펄스 여기, 내지 제N+1 RF 펄스 여기가 순차적으로 수행될 수 있다. 예컨대 N=2인 경우 제1 RF 펄스 여기(110), 제2 RF 펄스 여기(120), 및 제3 RF 펄스 여기(130)이 순차적으로 수행될 수 있다.
상술한 내용을 일반화 하면, 본 발명의 제3 실시예에 따른 K-스페이스데이터 획득방법은 아래의 단계들을 포함할 수 있다.
단계(S310):
제1 RF 펄스 여기 내지 제N RF 펄스 여기 각각에 대하여, 상기 제1 절편(ex: 11) 내지 제N 절편(ex: 13)에 대한 K-스페이스데이터를 획득할 수 있다.
이때, 제1 RF 펄스 여기 내지 제N RF 펄스 여기 각각에 대하여, 복수 개의 에코에서 서로 다른 주파수 영역에 대응하는 위상인코딩라인들을 획득할 수 있다. 즉, 제1 RF 펄스 여기 내지 제N RF 펄스 여기(ex: N=3) 각각에 대하여, 상기 복수 개의 에코들 중 제1에코에 대하여는 제1 저주파 위상인코딩라인들(ex: 811) 내지 제N 저주파 위상인코딩라인들(ex: 813)을 획득하고, 제2에코에 대하여는 제1 고주파 위상인코딩라인들(ex: 821) 내지 제N 고주파 위상인코딩라인들(ex: 823)을 획득할 수 있다. 이때, 바람직하게는 제1에코는 제2에코에 대하여 시간 축에서 선행하는 에코일 수 있다.
여기서, 상기 제1 RF 펄스 여기 내지 제N RF 펄스 여기 각각에 대하여, 즉, 상기 제1 절편(ex: 11) 내지 제N 절편(ex: 13)에 대하여, 상기 제1에코에 대하여 얻은 저주파 위상인코딩라인들(ex: 811, 812, 813)의 데이터를 제1 저주파 데이터 세트(ex: 311)라고 명명할 수 있다. 즉 상기 제1 저주파 데이터 세트는 상기 제1 저주파 위상인코딩라인들 내지 제N 저주파 위상인코딩라인들로 구성될 수 있다.
그리고 상기 제1 RF 펄스 여기 내지 제N RF 펄스 여기 각각에 대하여, 즉, 상기 제1 절편(ex: 11) 내지 제N 절편(ex: 13)에 대하여, 상기 제2에코에 대하여 얻은 고주파 위상인코딩라인들(ex: 821, 822, 823)의 데이터를 제2 고주파 데이터 세트(ex: 321)라고 명명할 수 있다. 즉 상기 제2 고주파 데이터 세트는 상기 제1 고주파 위상인코딩라인들 내지 제N 고주파 위상인코딩라인들로 구성될 수 있다.
단계(S320):
제N+1 RF 펄스 여기(ex: 제4 RF 펄스 여기)에 대하여, 상기 전체절편(10)에 대한 K-스페이스데이터를 획득할 수 있다.
이때, 제N+1 RF 펄스 여기에 대하여, 복수 개의 에코에서 서로 다른 주파수 영역에 대응하는 위상인코딩라인들을 획득할 수 있다. 즉, 제N+1 RF 펄스 여기에 대하여, 상기 복수 개의 에코들 중 제1에코에 대하여는 제N+1 고주파 위상인코딩라인들(ex: 213)을 획득하고, 제2에코에 대하여 제N+1 저주파 위상인코딩라인들(ex: 223)을 획득할 수 있다. 이때, 바람직하게는 제1에코는 제2에코에 대하여 시간 축에서 선행하는 에코일 수 있다.
여기서, 상기 제N+1 RF 펄스 여기(ex: 제4 RF 펄스 여기)에 대하여, 즉, 상기 전체절편(10)에 대하여, 상기 제1에코에 대하여 얻은 제N+1 고주파 위상인코딩라인들의 데이터를 제1 고주파 데이터 세트(ex: 312)라고 명명할 수 있다. 그리고 상기 제N+1 RF 펄스 여기에 대하여, 즉, 상기 전체절편에 대하여, 상기 제2에코에 대하여 얻은 제N+1 저주파 위상인코딩라인들의 데이터를 제2 저주파 데이터 세트(ex: 322)라고 명명할 수 있다.
단계(S330):
상기 제1 저주파 데이터 세트(ex: 311) 및 상기 제1 고주파 데이터 세트(ex: 312)를 서로 조합하여 제1 K-스페이스데이터를 생성하고, 상기 제2 고주파 데이터 세트(ex: 321) 및 상기 제2 저주파 데이터 세트(ex: 322)를 서로 조합하여 제2 K-스페이스데이터를 생성할 수 있다. 이렇게 생성된 상기 제1 K-스페이스데이터와 제2 K-스페이스데이터는 각각, 서로 다른 영상을 생성하는 데 사용될 수 있다.
<제4 실시예 - 획득한 K-스페이스데이터를 이용한 영상처리방법>
본 발명의 제4 실시예에 따른 영상처리방법은, 상술한 방법으로 획득한 K-스페이스데이터들을 조합하여 영상을 생성하는 기술에 관한 것이다.
본 발명의 제4 실시예에 따른 영상처리방법에 의해 생성하는 제1타입영상은, 아래 방법 1 및 방법 2 중 한 가지 방법에 의해 생성될 수 있다.
방법 1: 상술한 제1 K-스페이스데이터를 IFFT 하여 상기 제1타입영상을 획득하는 방법.
방법 2: (1) 상기 제1 저주파 위상인코딩라인들로 이루어진 제1 저주파 서브 K-스페이스데이터 내지 상기 제N 저주파 위상인코딩라인들로 이루어진 제N 저주파 서브 K-스페이스데이터를 각각 IFFT하여 제1 저주파 서브영상 내지 제N 저주파 서브영상을 생성하는 단계, (2) 상기 제N+1 고주파 위상인코딩라인들로 이루어진 제N+1 고주파 서브 K-스페이스데이터를 IFFT하여 제N+1 고주파 서브영상을 생성하는 단계, (3) 상기 제1 저주파 서브영상 내지 상기 제N 저주파 서브영상과 상기 제N+1 고주파 서브영상을 서로 결합하여 상기 제1타입영상을 생성하는 단계를 포함하는 방법.
본 발명의 다른 실시예에 따른 영상처리방법에 의해 생성하는 제2타입영상은, 아래 방법 3 및 방법 4 중 한 가지 방법에 의해 생성될 수 있다.
방법 3: 상술한 제2 K-스페이스데이터를 IFFT 하여 상기 제2타입영상을 획득하는 방법.
방법 4: (1) 상기 제1 고주파 위상인코딩라인들로 이루어진 제1 고주파 서브 K-스페이스데이터 내지 상기 제N 고주파 위상인코딩라인들로 이루어진 제N 고주파 서브 K-스페이스데이터를 각각 IFFT하여 제1 고주파 서브영상 내지 제N 고주파 서브영상을 생성하는 단계, (2) 상기 제N+1 저주파 위상인코딩라인들로 이루어진 제N+1 저주파 서브 K-스페이스데이터를 IFFT하여 제N+1 저주파 서브영상을 생성하는 단계, (3) 상기 제1 고주파 서브영상 내지 상기 제N 고주파 서브영상 및 상기 제N+1 저주파 서브영상을 서로 결합하여 상기 제2타입영상을 생성하는 단계를 포함하는 방법.
<제5 실시예>
본 발명의 제5 실시예는 복수 개의 에코를 이용한 이미지(데이터) 획득방법에 관한 것이다.
이하, 도 1 내지 도 4를 함께 참조하여 설명한다.
상기 데이터 획득방법은 아래와 같은 순서를 포함할 수 있다.
복수 절편(11, 12) 중 제1절편(11)만 여기하는 제1RF(110)를 복수 회 여기하는 제1프로세스를 수행할 수 있다. 이때, 상기 제1RF(110)를 각각 여기할 때 마다, 제1에코(111)에서는 미리 결정된 순서에 따라 저주파성분(211)을 얻고, 제2에코(112)에서는 미리 결정된 순서에 따라 고주파성분(221)을 얻을 수 있다.
상기 제1프로세스 이후, 복수 절편 중 제2절편(12)만 여기하는 제2RF(120)를 복수 회 여기하는 제2프로세스를 수행할 수 있다. 이때, 상기 제2RF(120)를 각각 여기할 때 마다, 제1에코(121)에서는 미리 결정된 순서에 따라 저주파성분(212)을 얻고, 제2에코(122)에서는 미리 결정된 순서에 따라 고주파성분(222)을 얻을 수 있다.
상기 제2프로세스 이후, 복수 절편 전체(=전체절편)(10)를 한 번에 여기하는 제3RF(130)를 복수 회 여기하는 제3프로세스를 수행할 수 있다. 이때, 상기 제3RF(130)를 각각 여기할 때마다, 제1에코(131)에서는 미리 결정된 순서에 따라 고주파성분(213)을 얻고, 제2에코(132)에서는 미리 결정된 순서에 따라 저주파성분(223)을 얻을 수 있다.
아래의 [표 1]은 각 여기대상 별 얻을 수 있는 주파수 대역을 나타낸 것이다.
여기대상 | 얻는 주파수 대역 | 얻는 주파수 대역 |
제1절편 | 저주파(211) | 고주파(221) |
제2절편 | 저주파(212) | 고주파(222) |
전체 | 고주파(213) | 저주파(223) |
도 8은 본 발명의 제5 실시예에 따라 동맥조영 이미지를 획득하기 위한 영상처리방법을 설명하기 위한 도면이다.
단계(S10)에서, 상기 제1프로세스의 상기 제1에코들(111)에서 획득한 제1 저주파 위상인코딩라인들(211)을 이미지 영역으로 변환(IFFT)하여 제1 저주파 서브영상(511)을 생성할 수 있다.
단계(S11)에서, 상기 제2프로세스의 상기 제1에코들(121)에서 획득한 제2 저주파 위상인코딩라인들(212)을 이미지 영역으로 변환(IFFT)하여 제2 저주파 서브영상(512)을 생성할 수 있다.
단계(S12)에서, 상기 제1 저주파 서브영상(511)과 상기 제2 저주파 서브영상(512)을 결합하여 제3 저주파 서브영상(611)을 생성할 수 있다.
단계(S13)에서, 상기 제3 저주파 서브영상(611)을 주파수 도메인으로 변환(FFT)하여 제1 저주파 데이터 세트(311)를 획득할 수 있다.
단계(S14)에서, 상기 제1 저주파 데이터 세트(311)와 상기 제3프로세스의 상기 제1에코들(131)에서 획득한 제3 고주파 위상인코딩라인들(312)을 주파수 영역에서 서로 결합하여 제1 K-스페이스데이터(210)를 획득할 수 있다.
도 9는 본 발명의 제5 실시예에 따라 정맥조영 이미지를 획득하기 위한 영상처리방법을 설명하기 위한 도면이다.
단계(S20)에서, 상기 제1프로세스의 상기 제2에코들(112)에서 획득한 제1 고주파 위상인코딩라인들(221)을 이미지 영역으로 변환하여 제1 고주파 서브영상(721)을 생성할 수 있다.
단계(S21)에서, 상기 제2프로세스의 상기 제2에코들(122)에서 획득한 제2 고주파 위상인코딩라인들(222)을 이미지 영역으로 변환하여 제2 고주파 서브영상(722)을 생성할 수 있다.
단계(S22)에서, 상기 제1 고주파 서브영상(721)과 상기 제2 고주파 서브영상(722)을 결합하여 제3 고주파 서브영상(751)을 생성할 수 있다.
단계(S23)에서, 상기 제3 고주파 서브영상(751)을 주파수 도메인으로 변환(FFT)하여 제2 고주파 데이터 세트(321)를 획득할 수 있다.
단계(S24)에서, 상기 제2 고주파 데이터 세트(321)와 상기 제3프로세스의 상기 제2에코들(132)에서 획득한 제3 저주파 위상인코딩라인들(322)을 주파수 영역에서 서로 결합하여 제2 K-스페이스데이터(220)를 획득할 수 있다.
도 8 내지 도 9에서 상술한 동맥조영 및 정맥조영 이미지 획득 방법은 상기 이미지를 고속으로 얻기 위한 방법일 수 있다.
상술한 내용과 달리, 일반 속도로 동맥조영 이미지 및 정맥조영 이미지를 획득하는 방법은 다음의 과정을 포함할 수 있다.
먼저, 동맥조영 이미지의 경우 상기 제1프로세스의 상기 제1에코들에서 획득한 제1 저주파 위상인코딩라인들(211), 상기 제2프로세스의 상기 제1에코들에서 획득한 제2 저주파 위상인코딩라인들(212), 및 상기 제3프로세스의 상기 제1에코들에서 획득한 제3 고주파 위상인코딩라인들(213)을 각각 이미지 영역으로 변환하여 이미지 영역에서 서로 결합하는 단계를 포함할 수 있다.
정맥조영 이미지의 경우, 상기 제1프로세스의 상기 제2에코들에서 획득한 제1 고주파 위상인코딩라인들(221), 상기 제2프로세스의 상기 제2에코들에서 획득한 제2 고주파 위상인코딩라인들(222), 및 상기 제3프로세스의 상기 제2에코들에서 획득한 제3 저주파 위상인코딩라인들(223)을 각각 이미지 영역으로 변환하여 이미지 영역에서 서로 결합하는 단계를 포함할 수 있다.
제5 실시예에서, 상기 제1에코(111, 121, 131)에서의 촬영조건과 상기 제2에코(112, 122, 132)에서의 촬영조건은 서로 다를 수 있다.
또한, 상기 제1에코(111, 121, 131)는 상기 제2에코(112, 122, 132)보다 선행할 수 있다. 이때, 제1에코 시에는 Time of Flight MR angiogram 을 얻고, 제2에코 시에는 Susceptiblitiy weighted imaging 을 얻을 수 있다.
또한, 상기 제1에코(111, 121, 131) 시의 절편의 개수는 제2에코(112, 122, 132) 시의 절편의 개수보다 많을 수 있다.
그리고 도 8 및 도 9에서 상술한 바와 같이, 상기 제1에코(111, 121, 131) 시에는 동맥 이미지를 획득할 수 있고, 상기 제2에코 시에는 정맥 이미지(112, 122, 132)를 획득할 수 있다.
<제6 실시예>
제6 실시예는 동맥조영 이미지 및 정맥조영 이미지를 획득하는 다른 방법에 관한 것이다. 제6 실시예는 상술한 제3 실시예에서, N=3인 경우에 해당할 수 있다.
도 10은 본 발명의 제6 실시예에 따라 동맥조영 이미지 및 정맥조영 이미지를 획득하기 위한 영상처리방법을 설명하기 위한 도면이다.
도 10의 (a)는 동맥조영술에 대한 영상처리순서를 나타낸 것이며, 도 10의 (b)는 정맥조영술에 대한 영상처리순서를 나타낸 것이다.
단계(S31)에서, 제1절편(11), 제2절편(12), 및 제3절편(13)의 각 제1에코들에서 획득한 저주파 위상인코딩라인들을 각각 이미지 영역으로 변환(IFFT)하여 제1 저주파 영상(941)을 생성할 수 있다. 이때, 제1 저주파 영상(941)은 상기 각각 이미지 영역으로 변환하여 얻은 저주파 서브영상들을 상기 각 절편(11, 12, 13)의 위치에 따라 배치하여 결합한 것일 수 있다.
단계(S32)에서, 전체절편(10)의 제1에코들에서 획득한 고주파 위상인코딩라인들(912)을 이미지 영역으로 변환(IFFT)하여 제1 고주파 영상(942)을 생성할 수 있다.
단계(S33)에서, 상기 제1 저주파 영상(941)과 상기 제2 저주파 영상(942)을 결합하여 제1 전체영상(943)을 생성할 수 있다.
단계(S34)에서, 상기 제1 전체영상(943)을 MIP(Maximum Intensity Projection)하여 동맥조영 이미지(944)를 획득할 수 있다.
다음, 정맥조영 이미지를 획득하는 방법은 아래와 같다.
단계(S41)에서, 제1절편(11), 제2절편(12), 및 제3절편(13)의 각 제2에코들에서 획득한 고주파 위상인코딩라인들을 각각 이미지 영역으로 변환(IFFT)하여 제2 고주파 영상(951)을 생성할 수 있다. 이때, 제1 고주파 영상(951)은 상기 각각 이미지 영역으로 변환하여 얻은 저주파 서브영상들을 상기 각 절편(11, 12, 13)의 위치에 따라 배치하여 결합한 것일 수 있다.
단계(S42)에서, 전체절편(10)의 제2에코들에서 획득한 저주파 위상인코딩라인들(922)을 이미지 영역으로 변환(IFFT)하여 제2 저주파 영상(952)을 생성할 수 있다.
단계(S43)에서, 상기 제2 고주파 영상(951)과 상기 제2 저주파 영상(952)을 결합하여 제2 전체영상(953)을 생성할 수 있다.
단계(S44)에서, 상기 제2 전체영상(953)을 위상 마스크 필터링(phase mask filtering)과 minIP(Miimum Intensity Projecton)을 통하여 정맥조영 이미지(954)를 획득할 수 있다.
<제7 실시예>
제7 실시예는 복수 개의 에코를 이용한 다른 데이터 획득방법에 관한 것이다.
상기 데이터 획득방법에 따르면, 먼저 제1절편(11)에 대하여 제1RF 신호(110)를 여기하여 발생한 제1에코(111) 및 제2에코(112) 각각에 대하여 K-스페이스의 라인들 중 저주파 라인들(211) 및 고주파 라인들(221)을 획득할 수 있다.
그리고 제2절편(12)에 대하여 제2RF 신호(120)를 여기하여 발생한 제3에코(121) 및 제4에코(122) 각각에 대하여 K-스페이스의 라인들 중 저주파 라인들(212) 및 고주파 라인들(222)을 획득할 수 있다.
그리고 상기 제1절편(11) 및 상기 제2절편(12)을 포함하는 전체절편(10)에 대하여 제3RF 신호(130)를 여기하여 발생한 제5에코(131)에서는 고주파 라인들(213)을 획득하고 제6에코(132)에서는 저주파 라인들(223)을 획득할 수 있다.
이때, 상기 제1에코가 제2에코에 대해 선행하거나, 또는 상기 제2에코가 제1에코에 대해 선행할 수 있다.
<제8 실시예>
제8 실시예는 복수 개의 에코를 이용한 또 다른 데이터 획득방법에 관한 것이다.
상기 데이터 획득방법에 따르면, 먼저 제1절편(11) 및 제2절편(12) 각각에 대해 제1에코(111, 121)에서는 제1오더링에 따른 일부 라인들(211, 212)을 획득하고 제2에코(112, 122)에서는 제2오더링에 따른 일부 라인들(221, 222)을 획득할 수 있다.
그리고 전체절편(10)에 대하여 제1에코(131)에서는 상기 제1오더링에 따른 잔여 라인들 중 적어도 일부(213)를 획득하고, 제2에코(132)에서는 상기 제2오더링에 따른 잔여 라인들 중 적어도 일부(223)를 획득할 수 있다.
<제9 실시예>
제9 실시예는 복수 개의 에코를 이용한 또 다른 데이터 획득방법에 관한 것이다.
상기 데이터 획득방법에 따르면, 제1절편(11) 및 제2절편(12) 각각에 대해 제1오더링에 따른 일부 라인들(211, 212)을 획득할 수 있다. 그리고 전체절편(10)에 대하여 상기 제1오더링에 따른 잔여 라인들(213)을 획득할 수 있다. 이후, 상기 일부 라인들과 잔여 라인들에 대한 정보를 서로 결합할 수 있다.
<제10 실시예>
제10 실시예는 복수 개의 에코를 이용한 또 다른 데이터 획득방법에 관한 것이다.
상기 데이터 획득방법에 따르면, 전체절편(10)을 분할하여 정의한 복수 개의 절편들(11, 12) 각각에 대해 제1에코(111, 121)에서는 저주파 라인(211, 212)을 획득하고 제2에코(112, 122)에서는 고주파 라인(221, 222)을 획득할 수 있다. 그리고 상기 전체절편(10)에 대하여 상기 제1에코(131)에서는 고주파 라인(213)을 획득하고 제2에코(132)에서는 저주파 라인(223)을 획득할 수 있다.
도 11a 및 도 11b는 본 발명의 일 실시예를 이용한 동맥조영술에 따라 획득한 이미지와 종래의 기술에 따른 동맥조영술에 따라 획득한 이미지의 차이를 예시한 것이다.
도 11a는 관찰대상을 종래의 기술에 따라 3개의 절편으로 나누어 관찰하였을 때에 얻은 동맥조영도이다. 도 11a에 따르면 절편들 간의 경계부분(화살표 표시)에 가로줄 형태의 인공물이 나타난다는 점을 확인할 수 있다.
도 11b는 관찰대상을 본 발명의 일 실시예에 따라 3개의 절편으로 나누어 관찰하였을 때에 얻은 동맥조영도이다. 도 11b에 따르면 절편들 간의 경계부분(화살표 표시)에 가로줄 형태의 인공물이 나타나지 않는다는 점을 확인할 수 있다. 즉, 도 11b의 이미지는 도 11a의 이미지에 비하여 개선되었음을 이해할 수 있다. 따라서 본 발명의 일 실시예에 따른 방법을 사용하는 경우, 종래기술에 비하여 동맥조영도의 품질이 향상됨을 알 수 있다. 이러한 동맥조영도의 품질향상은, 본 발명에 따른 실시예에서 멀티 에코를 사용하는 경우에도 얻을 수 있으며, 싱글 에코를 사용하는 경우에도 얻을 수 있다.
상술한 본 발명의 실시예들을 이용하여, 본 발명의 기술 분야에 속하는 자들은 본 발명의 본질적인 특성에서 벗어나지 않는 범위 내에 다양한 변경 및 수정을 용이하게 실시할 수 있을 것이다. 특허청구범위의 각 청구항의 내용은 본 명세서를 통해 이해할 수 있는 범위 내에서 인용관계가 없는 다른 청구항에 결합될 수 있다.
10 : 전체절편 11, 12 : 절편
110 : 제1RF 111 : 제1에코
112 : 제2에코 120 : 제2RF
121: 제1에코 122 : 제2에코
130 : 제3RF 131 : 제1에코
132 : 제2에코 210 : 제1 K-스페이스데이터
211 : 저주파성분 212 : 저주파성분
213 : 고주파성분 220 : 제2 K-스페이스데이터
221 : 고주파성분 222 : 고주파성분
223 : 저주파성분 311 : 제1 저주파 데이터 세트
312 : 제3 고주파 위상인코딩라인들 321 : 제2 고주파 데이터 세트
322 : 제3 저주파 위상인코딩라인들 421 : 제3 고주파 서브영상
511 : 제1 저주파 서브영상 512 : 제2 저주파 서브영상
611 : 제3 저주파 서브영상 721 : 제1 고주파 서브영상
722 : 제2 고주파 서브영상
110 : 제1RF 111 : 제1에코
112 : 제2에코 120 : 제2RF
121: 제1에코 122 : 제2에코
130 : 제3RF 131 : 제1에코
132 : 제2에코 210 : 제1 K-스페이스데이터
211 : 저주파성분 212 : 저주파성분
213 : 고주파성분 220 : 제2 K-스페이스데이터
221 : 고주파성분 222 : 고주파성분
223 : 저주파성분 311 : 제1 저주파 데이터 세트
312 : 제3 고주파 위상인코딩라인들 321 : 제2 고주파 데이터 세트
322 : 제3 저주파 위상인코딩라인들 421 : 제3 고주파 서브영상
511 : 제1 저주파 서브영상 512 : 제2 저주파 서브영상
611 : 제3 저주파 서브영상 721 : 제1 고주파 서브영상
722 : 제2 고주파 서브영상
Claims (14)
- 한 개 이상의 에코를 이용한 이미지 획득방법으로서,
복수 절편 중 제1절편만 여기하는 제1RF를 복수 회 여기하는 제1프로세스를 포함하며, 상기 제1RF를 각각 여기할 때마다, 제1에코에서는 미리 결정된 순서에 따라 저주파성분을 얻고, 제2에코에서는 미리 결정된 순서에 따라 고주파성분을 얻으며,
상기 제1프로세스 이후,
복수 절편 중 제2절편만 여기하는 제2RF를 복수 회 여기하는 제2프로세스를 포함하며, 상기 제2RF를 각각 여기할 때 마다, 제1에코에서는 미리 결정된 순서에 따라 저주파성분을 얻고, 제2에코에서는 미리 결정된 순서에 따라 고주파성분을 얻으며,
상기 제2프로세스 이후,
복수 절편 전체 전부를 한 번에 여기하는 제3RF를 복수 회 여기하는 제3프로세스를 포함하며, 상기 제3RF를 각각 여기할 때마다, 제1에코에서는 미리 결정된 순서에 따라 고주파성분을 얻고, 제2에코에서는 미리 결정된 순서에 따라 저주파성분을 얻는,
데이터 획득방법. - 제1항에 있어서,
상기 제1프로세스의 상기 제1에코들에서 획득한 제1 저주파 위상인코딩라인들을 이미지 영역으로 변환하여 제1 저주파 서브영상을 생성하는 단계;
상기 제2프로세스의 상기 제1에코들에서 획득한 제2 저주파 위상인코딩라인들을 이미지 영역으로 변환하여 제2 저주파 서브영상을 생성하는 단계;
상기 제1 저주파 서브영상과 상기 제2 저주파 서브영상을 결합하여 제3 저주파 서브영상을 생성하는 단계;
상기 제3 저주파 서브영상을 주파수 도메인으로 변환하여 제1 저주파 데이터 세트를 획득하는 단계;
상기 제1 저주파 데이터 세트와
상기 제3프로세스의 상기 제1에코들에서 획득한 제3 고주파 위상인코딩라인들을 주파수 영역에서 서로 결합하여 제1 K-스페이스데이터를 획득하는 단계;
를 포함하는,
데이터 획득방법. - 제1항에 있어서,
상기 제1프로세스의 상기 제2에코들에서 획득한 제1 고주파 위상인코딩라인들을 이미지 영역으로 변환하여 제1 고주파 서브영상을 생성하는 단계;
상기 제2프로세스의 상기 제2에코들에서 획득한 제2 고주파 위상인코딩라인들을 이미지 영역으로 변환하여 제2 고주파 서브영상을 생성하는 단계;
상기 제1 고주파 서브영상과 상기 제2 고주파 서브영상을 결합하여 제3 고주파 서브영상을 생성하는 단계;
상기 제3 고주파 서브영상을 주파수 도메인으로 변환하여 제2 고주파 데이터 세트를 획득하는 단계;
상기 제2 고주파 데이터 세트와 상기 제3프로세스의 상기 제2에코들에서 획득한 제3 저주파 위상인코딩라인들을 주파수 영역에서 서로 결합하여 제2 K-스페이스데이터를 획득하는 단계;
를 포함하는,
데이터 획득방법. - 제1항에 있어서,
상기 제1프로세스의 상기 제1에코들에서 획득한 제1 저주파 위상인코딩라인들,
상기 제2프로세스의 상기 제1에코들에서 획득한 제2 저주파 위상인코딩라인들, 및
상기 제3프로세스의 상기 제1에코들에서 획득한 제3 고주파 위상인코딩라인들
을 각각 이미지 영역으로 변환하여 이미지 영역에서 서로 결합하는 단계;
를 포함하는,
데이터 획득방법. - 제1항에 있어서,
상기 제1프로세스의 상기 제2에코들에서 획득한 제1 고주파 위상인코딩라인들,
상기 제2프로세스의 상기 제2에코들에서 획득한 제2 고주파 위상인코딩라인들, 및
상기 제3프로세스의 상기 제2에코들에서 획득한 제3 저주파 위상인코딩라인들
을 각각 이미지 영역으로 변환하여 이미지 영역에서 서로 결합하는 단계;
를 포함하는,
데이터 획득방법. - 제1항에 있어서, 상기 제1에코에서의 촬영조건과 상기 제2에코에서의 촬영조건은 서로 다른 것을 특징으로 하는, 데이터 획득방법.
- 제1항에 있어서,
상기 제1에코는 상기 제2에코보다 선행하며,
제1에코 시에는 Time of Flight MR angiogram 을 얻고,
제2에코 시에는 Susceptiblitiy weighted imaging 을 얻는 것을 특징으로 하는,
데이터 획득방법. - 제1항에 있어서,
상기 제1에코는 상기 제2에코보다 선행하며,
상기 제1에코 시의 절편의 개수는 제2에코 시의 절편의 개수보다 많은 것을 특징으로 하는,
데이터 획득방법. - 제8항에 있어서,
상기 제1에코 시에는 동맥 이미지를 획득하고
상기 제2에코 시에는 정맥 이미지를 획득하도록 되어 있는,
데이터 획득방법. - 한 개 이상의 에코를 이용한 데이터획득방법으로서,
제1절편에 대하여 제1RF 신호를 여기하여 발생한 제1에코 및 제2에코 각각에 대하여 K-스페이스의 라인들 중 저주파 라인들 및 고주파 라인들을 획득하고,
제2절편에 대하여 제2RF 신호를 여기하여 발생한 제3에코 및 제4에코 각각에 대하여 K-스페이스의 라인들 중 저주파 라인들 및 고주파 라인들을 획득하고,
상기 제1절편 및 상기 제2절편을 포함하는 전체절편에 대하여 제3RF 신호를 여기하여 발생한 제5에코에서는 고주파 라인들을 획득하고 제6에코에서는 저주파 라인들을 획득하는 것을 특징으로 하는,
데이터 획득방법. - 제10항에 있어서,
상기 제1에코가 제2에코에 대해 선행하거나, 또는
상기 제2에코가 제1에코에 대해 선행하는,
데이터 획득방법. - 한 개 이상의 에코를 이용한 데이터획득방법으로서,
제1절편 및 제2절편 각각에 대해 제1에코에서는 제1오더링에 따른 일부 라인들을 획득하고 제2에코에서는 제2오더링에 따른 일부 라인들을 획득하고,
전체절편에 대하여 제1에코에서는 상기 제1오더링에 따른 잔여 라인들 중 적어도 일부를 획득하고, 제2에코에서는 상기 제2오더링에 따른 잔여 라인들 중 적어도 일부를 획득하는 것을 특징으로 하는,
데이터 획득방법. - 제1절편 및 제2절편 각각에 대해 제1오더링에 따른 일부 라인들을 획득하고,
전체절편에 대하여 상기 제1오더링에 따른 잔여 라인들을 획득하고,
상기 일부 라인들과 잔여 라인들에 대한 정보를 서로 결합하는 단계
를 포함하는,
데이터 획득방법. - 한 개 이상의 에코를 이용한 데이터획득방법으로서,
전체절편을 분할하여 정의한 복수 개의 절편들 각각에 대해 제1에코에서는 저주파 라인을 획득하고 제2에코에서는 고주파 라인을 획득하며,
상기 전체절편에 대하여 상기 제1에코에서는 고주파 라인을 획득하고 제2에코에서는 저주파 라인을 획득하는 것을 특징으로 하는,
데이터 획득방법.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020170035099A KR101974199B1 (ko) | 2017-03-21 | 2017-03-21 | 가변절편 자기공명영상 데이터 획득방법 |
PCT/KR2017/011375 WO2018174360A1 (ko) | 2017-03-21 | 2017-10-16 | 가변절편 자기공명영상 데이터 획득방법 |
US16/496,097 US11209512B2 (en) | 2017-03-21 | 2017-10-16 | Method for acquiring variable slab magnetic resonance imaging data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
KR1020170035099A KR101974199B1 (ko) | 2017-03-21 | 2017-03-21 | 가변절편 자기공명영상 데이터 획득방법 |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20180106567A KR20180106567A (ko) | 2018-10-01 |
KR101974199B1 true KR101974199B1 (ko) | 2019-04-30 |
Family
ID=63586125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020170035099A KR101974199B1 (ko) | 2017-03-21 | 2017-03-21 | 가변절편 자기공명영상 데이터 획득방법 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11209512B2 (ko) |
KR (1) | KR101974199B1 (ko) |
WO (1) | WO2018174360A1 (ko) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102001874B1 (ko) * | 2018-01-05 | 2019-07-19 | 한국과학기술원 | 고정 rf 코일과 자유이동 rf 코일의 조합을 이용하여 mri 이미지의 snr을 실시간으로 향상하는 방법 및 이를 이용한 mri 데이터 처리장치 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090278538A1 (en) | 2008-05-07 | 2009-11-12 | Jyh-Horng Chen | Method and apparatus for simultaneously acquiring multiple slices/slabs in magnetic resonance system |
US20110213237A1 (en) | 2008-09-15 | 2011-09-01 | University of Pittsburgh -Of the Commonwealth [Syst] System of Higher Education | Echo-specific k-space reordering approach to compatible dual-echo arteriovenography |
US20110275926A1 (en) | 2008-01-23 | 2011-11-10 | The Regents Of The University Of Colorado | Susceptibility Weighted Magnetic Resonance Imaging Of Venous Vasculature |
US20130225976A1 (en) | 2011-05-02 | 2013-08-29 | Toshiba Medical Systems Corporation | Efficient multi-station mri |
US20140081123A1 (en) | 2012-09-18 | 2014-03-20 | General Electric Company | Magnetic resonance angiography and venography |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011082009B3 (de) * | 2011-09-01 | 2013-02-28 | Siemens Aktiengesellschaft | Schichtselektive Detektion und Korrektur von fehlerhaften Magnetresonanz-Bilddaten in Schicht-Multiplexing-Messsequenzen |
US10203387B2 (en) * | 2013-06-06 | 2019-02-12 | Koninklijke Philips N.V. | MR imaging with enhanced susceptibility contrast |
JP6637053B2 (ja) * | 2015-02-03 | 2020-01-29 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 磁化率重み付け磁気共鳴イメージングのための方法及びシステム |
-
2017
- 2017-03-21 KR KR1020170035099A patent/KR101974199B1/ko active IP Right Grant
- 2017-10-16 WO PCT/KR2017/011375 patent/WO2018174360A1/ko active Application Filing
- 2017-10-16 US US16/496,097 patent/US11209512B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110275926A1 (en) | 2008-01-23 | 2011-11-10 | The Regents Of The University Of Colorado | Susceptibility Weighted Magnetic Resonance Imaging Of Venous Vasculature |
US20090278538A1 (en) | 2008-05-07 | 2009-11-12 | Jyh-Horng Chen | Method and apparatus for simultaneously acquiring multiple slices/slabs in magnetic resonance system |
US20110213237A1 (en) | 2008-09-15 | 2011-09-01 | University of Pittsburgh -Of the Commonwealth [Syst] System of Higher Education | Echo-specific k-space reordering approach to compatible dual-echo arteriovenography |
US20130225976A1 (en) | 2011-05-02 | 2013-08-29 | Toshiba Medical Systems Corporation | Efficient multi-station mri |
US20140081123A1 (en) | 2012-09-18 | 2014-03-20 | General Electric Company | Magnetic resonance angiography and venography |
Non-Patent Citations (1)
Title |
---|
듀얼에코기법을 이용한 T1과 T2* 영상 동시획득(카이스트 석사학위논문, 도원준, 2015) |
Also Published As
Publication number | Publication date |
---|---|
WO2018174360A1 (ko) | 2018-09-27 |
US20200096585A1 (en) | 2020-03-26 |
KR20180106567A (ko) | 2018-10-01 |
US11209512B2 (en) | 2021-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8023705B2 (en) | Method for reconstructing image from echo planar imaging sequence | |
US8742757B2 (en) | Magnetic resonance imaging with high spatial and temporal resolution | |
US20050203377A1 (en) | Method and apparatus for anatomically tailored k-space sampling and recessed elliptical view ordering for bolus-enhanced 3D MR angiography | |
US7898253B2 (en) | Method and apparatus for removing artifacts during magnetic resonance imaging | |
JP2014147756A (ja) | 磁気共鳴システムの作動方法および制御装置 | |
JP2010042245A (ja) | 磁気共鳴イメージング装置 | |
EP2026084A2 (en) | Single coil parallel imaging | |
JP2016022302A (ja) | 磁気共鳴信号処理方法、磁気共鳴信号処理装置及び磁気共鳴装置並びにプログラム | |
KR101974199B1 (ko) | 가변절편 자기공명영상 데이터 획득방법 | |
KR20100010481A (ko) | 초고자장 mri 장비에서 epi 영상 왜곡 보정 방법 | |
EP3683595A1 (en) | Segment-specific undersampling in readout-segmented magnetic resonance imaging | |
JP5675362B2 (ja) | 磁気共鳴撮像における脂肪/水の分離を強化するためのシステム、方法及びマシン読み取り可能プログラム | |
KR20090075644A (ko) | 부분 일시적-보간을 갖는 분할-에코 | |
US11965950B2 (en) | Slice ordering for MB-EPI ASL imaging | |
Xie et al. | Adaptive online self‐gating (ADIOS) for free‐breathing noncontrast renal MR angiography | |
KR101623690B1 (ko) | 자기 공명 영상의 b1 맵 생성 장치 및 방법 | |
JP2008279102A (ja) | 磁気共鳴イメージング装置 | |
JP6554723B2 (ja) | 磁気共鳴装置およびプログラム | |
JP6576726B2 (ja) | 磁気共鳴装置 | |
US9739858B2 (en) | Magnetic resonance imaging apparatus with segmented data acquisition | |
US20240111009A1 (en) | Magnetic resonance imaging apparatus and control method thereof | |
JP2005522250A (ja) | 複数の基礎データセットから複合オブジェクトデータセットを形成するデータ処理 | |
JPWO2006087934A1 (ja) | 磁気共鳴イメージング方法及び装置 | |
JP2016040024A (ja) | 磁気共鳴イメージング装置 | |
JP2004329269A (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 |