KR101974199B1 - 가변절편 자기공명영상 데이터 획득방법 - Google Patents

가변절편 자기공명영상 데이터 획득방법 Download PDF

Info

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
Application number
KR1020170035099A
Other languages
English (en)
Other versions
KR20180106567A (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 한국과학기술원
Priority to KR1020170035099A priority Critical patent/KR101974199B1/ko
Priority to PCT/KR2017/011375 priority patent/WO2018174360A1/ko
Priority to US16/496,097 priority patent/US11209512B2/en
Publication of KR20180106567A publication Critical patent/KR20180106567A/ko
Application granted granted Critical
Publication of KR101974199B1 publication Critical patent/KR101974199B1/ko

Links

Images

Classifications

    • 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/543Control 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
    • 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/563Image 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/5635Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]
    • 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/446Multifrequency selective RF pulses, e.g. multinuclear acquisition mode
    • 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/4818MR 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
    • 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
    • G01R33/56536Correction of image distortions, e.g. due to magnetic field inhomogeneities due to magnetic susceptibility variations
    • 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/4818MR 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/482MR 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/4822MR 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
    • 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

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

가변절편 자기공명영상 데이터 획득방법{A method for variable-slab MRI data acquisition}
본 발명은 자기공명영상을 위한 데이터 획득방법에 관한 것으로서, 한 개 이상의 에코를 이용하며, 한 번의 촬영으로 동맥조영도와 정맥조영도를 한 번에 얻거나 동맥조영도의 질을 개선할 수 있는 기술에 관한 것이다.
자기공명영상은 영상 기술 중 하나로 핵자기공명 원리를 사용한다. 자기장을 발생하는 자기공명 촬영 장치에 인체를 넣고 고주파를 발생시키면 신체의 수소 원자의 전자가 공명하게 된다. 이때 나오는 신호의 차이를 측정하고 컴퓨터를 통해 재구성하여 영상화시키면 자기공명 영상이 된다.
자기공명영상을 얻기 위하여, 동맥조영술(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개의 절편으로 나누어 관찰하였을 때에 얻은 동맥조영도이다.
이하, 본 발명의 실시예를 첨부한 도면을 참고하여 설명한다. 그러나 본 발명은 본 명세서에서 설명하는 실시예에 한정되지 않으며 여러 가지 다른 형태로 구현될 수 있다. 본 명세서에서 사용되는 용어는 실시예의 이해를 돕기 위한 것이며, 본 발명의 범위를 한정하고자 의도된 것이 아니다. 또한, 이하에서 사용되는 단수 형태들은 문구들이 이와 명백히 반대의 의미를 나타내지 않는 한 복수 형태들도 포함한다.
<제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 고주파 서브영상

Claims (14)

  1. 한 개 이상의 에코를 이용한 이미지 획득방법으로서,
    복수 절편 중 제1절편만 여기하는 제1RF를 복수 회 여기하는 제1프로세스를 포함하며, 상기 제1RF를 각각 여기할 때마다, 제1에코에서는 미리 결정된 순서에 따라 저주파성분을 얻고, 제2에코에서는 미리 결정된 순서에 따라 고주파성분을 얻으며,
    상기 제1프로세스 이후,
    복수 절편 중 제2절편만 여기하는 제2RF를 복수 회 여기하는 제2프로세스를 포함하며, 상기 제2RF를 각각 여기할 때 마다, 제1에코에서는 미리 결정된 순서에 따라 저주파성분을 얻고, 제2에코에서는 미리 결정된 순서에 따라 고주파성분을 얻으며,
    상기 제2프로세스 이후,
    복수 절편 전체 전부를 한 번에 여기하는 제3RF를 복수 회 여기하는 제3프로세스를 포함하며, 상기 제3RF를 각각 여기할 때마다, 제1에코에서는 미리 결정된 순서에 따라 고주파성분을 얻고, 제2에코에서는 미리 결정된 순서에 따라 저주파성분을 얻는,
    데이터 획득방법.
  2. 제1항에 있어서,
    상기 제1프로세스의 상기 제1에코들에서 획득한 제1 저주파 위상인코딩라인들을 이미지 영역으로 변환하여 제1 저주파 서브영상을 생성하는 단계;
    상기 제2프로세스의 상기 제1에코들에서 획득한 제2 저주파 위상인코딩라인들을 이미지 영역으로 변환하여 제2 저주파 서브영상을 생성하는 단계;
    상기 제1 저주파 서브영상과 상기 제2 저주파 서브영상을 결합하여 제3 저주파 서브영상을 생성하는 단계;
    상기 제3 저주파 서브영상을 주파수 도메인으로 변환하여 제1 저주파 데이터 세트를 획득하는 단계;
    상기 제1 저주파 데이터 세트와
    상기 제3프로세스의 상기 제1에코들에서 획득한 제3 고주파 위상인코딩라인들을 주파수 영역에서 서로 결합하여 제1 K-스페이스데이터를 획득하는 단계;
    를 포함하는,
    데이터 획득방법.
  3. 제1항에 있어서,
    상기 제1프로세스의 상기 제2에코들에서 획득한 제1 고주파 위상인코딩라인들을 이미지 영역으로 변환하여 제1 고주파 서브영상을 생성하는 단계;
    상기 제2프로세스의 상기 제2에코들에서 획득한 제2 고주파 위상인코딩라인들을 이미지 영역으로 변환하여 제2 고주파 서브영상을 생성하는 단계;
    상기 제1 고주파 서브영상과 상기 제2 고주파 서브영상을 결합하여 제3 고주파 서브영상을 생성하는 단계;
    상기 제3 고주파 서브영상을 주파수 도메인으로 변환하여 제2 고주파 데이터 세트를 획득하는 단계;
    상기 제2 고주파 데이터 세트와 상기 제3프로세스의 상기 제2에코들에서 획득한 제3 저주파 위상인코딩라인들을 주파수 영역에서 서로 결합하여 제2 K-스페이스데이터를 획득하는 단계;
    를 포함하는,
    데이터 획득방법.
  4. 제1항에 있어서,
    상기 제1프로세스의 상기 제1에코들에서 획득한 제1 저주파 위상인코딩라인들,
    상기 제2프로세스의 상기 제1에코들에서 획득한 제2 저주파 위상인코딩라인들, 및
    상기 제3프로세스의 상기 제1에코들에서 획득한 제3 고주파 위상인코딩라인들
    을 각각 이미지 영역으로 변환하여 이미지 영역에서 서로 결합하는 단계;
    를 포함하는,
    데이터 획득방법.
  5. 제1항에 있어서,
    상기 제1프로세스의 상기 제2에코들에서 획득한 제1 고주파 위상인코딩라인들,
    상기 제2프로세스의 상기 제2에코들에서 획득한 제2 고주파 위상인코딩라인들, 및
    상기 제3프로세스의 상기 제2에코들에서 획득한 제3 저주파 위상인코딩라인들
    을 각각 이미지 영역으로 변환하여 이미지 영역에서 서로 결합하는 단계;
    를 포함하는,
    데이터 획득방법.
  6. 제1항에 있어서, 상기 제1에코에서의 촬영조건과 상기 제2에코에서의 촬영조건은 서로 다른 것을 특징으로 하는, 데이터 획득방법.
  7. 제1항에 있어서,
    상기 제1에코는 상기 제2에코보다 선행하며,
    제1에코 시에는 Time of Flight MR angiogram 을 얻고,
    제2에코 시에는 Susceptiblitiy weighted imaging 을 얻는 것을 특징으로 하는,
    데이터 획득방법.
  8. 제1항에 있어서,
    상기 제1에코는 상기 제2에코보다 선행하며,
    상기 제1에코 시의 절편의 개수는 제2에코 시의 절편의 개수보다 많은 것을 특징으로 하는,
    데이터 획득방법.
  9. 제8항에 있어서,
    상기 제1에코 시에는 동맥 이미지를 획득하고
    상기 제2에코 시에는 정맥 이미지를 획득하도록 되어 있는,
    데이터 획득방법.
  10. 한 개 이상의 에코를 이용한 데이터획득방법으로서,
    제1절편에 대하여 제1RF 신호를 여기하여 발생한 제1에코 및 제2에코 각각에 대하여 K-스페이스의 라인들 중 저주파 라인들 및 고주파 라인들을 획득하고,
    제2절편에 대하여 제2RF 신호를 여기하여 발생한 제3에코 및 제4에코 각각에 대하여 K-스페이스의 라인들 중 저주파 라인들 및 고주파 라인들을 획득하고,
    상기 제1절편 및 상기 제2절편을 포함하는 전체절편에 대하여 제3RF 신호를 여기하여 발생한 제5에코에서는 고주파 라인들을 획득하고 제6에코에서는 저주파 라인들을 획득하는 것을 특징으로 하는,
    데이터 획득방법.
  11. 제10항에 있어서,
    상기 제1에코가 제2에코에 대해 선행하거나, 또는
    상기 제2에코가 제1에코에 대해 선행하는,
    데이터 획득방법.
  12. 한 개 이상의 에코를 이용한 데이터획득방법으로서,
    제1절편 및 제2절편 각각에 대해 제1에코에서는 제1오더링에 따른 일부 라인들을 획득하고 제2에코에서는 제2오더링에 따른 일부 라인들을 획득하고,
    전체절편에 대하여 제1에코에서는 상기 제1오더링에 따른 잔여 라인들 중 적어도 일부를 획득하고, 제2에코에서는 상기 제2오더링에 따른 잔여 라인들 중 적어도 일부를 획득하는 것을 특징으로 하는,
    데이터 획득방법.
  13. 제1절편 및 제2절편 각각에 대해 제1오더링에 따른 일부 라인들을 획득하고,
    전체절편에 대하여 상기 제1오더링에 따른 잔여 라인들을 획득하고,
    상기 일부 라인들과 잔여 라인들에 대한 정보를 서로 결합하는 단계
    를 포함하는,
    데이터 획득방법.
  14. 한 개 이상의 에코를 이용한 데이터획득방법으로서,
    전체절편을 분할하여 정의한 복수 개의 절편들 각각에 대해 제1에코에서는 저주파 라인을 획득하고 제2에코에서는 고주파 라인을 획득하며,
    상기 전체절편에 대하여 상기 제1에코에서는 고주파 라인을 획득하고 제2에코에서는 저주파 라인을 획득하는 것을 특징으로 하는,
    데이터 획득방법.
KR1020170035099A 2017-03-21 2017-03-21 가변절편 자기공명영상 데이터 획득방법 KR101974199B1 (ko)

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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102001874B1 (ko) * 2018-01-05 2019-07-19 한국과학기술원 고정 rf 코일과 자유이동 rf 코일의 조합을 이용하여 mri 이미지의 snr을 실시간으로 향상하는 방법 및 이를 이용한 mri 데이터 처리장치

Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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. 磁化率重み付け磁気共鳴イメージングのための方法及びシステム

Patent Citations (5)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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