KR20170028829A - 자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치 - Google Patents

자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치 Download PDF

Info

Publication number
KR20170028829A
KR20170028829A KR1020160087383A KR20160087383A KR20170028829A KR 20170028829 A KR20170028829 A KR 20170028829A KR 1020160087383 A KR1020160087383 A KR 1020160087383A KR 20160087383 A KR20160087383 A KR 20160087383A KR 20170028829 A KR20170028829 A KR 20170028829A
Authority
KR
South Korea
Prior art keywords
coil
space
dictionary
image
unit
Prior art date
Application number
KR1020160087383A
Other languages
English (en)
Other versions
KR102527812B1 (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 EP16842129.5A priority Critical patent/EP3346443B1/en
Priority to PCT/KR2016/008494 priority patent/WO2017039163A1/ko
Priority to US15/756,666 priority patent/US20180247436A1/en
Publication of KR20170028829A publication Critical patent/KR20170028829A/ko
Application granted granted Critical
Publication of KR102527812B1 publication Critical patent/KR102527812B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T1/00—General purpose image data processing
    • 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
    • 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • 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/561—Image 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/5611—Parallel 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/60—Editing figures and text; Combining figures or text
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20212—Image combination
    • G06T2207/20221—Image fusion; Image merging
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/41—Medical

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Theoretical Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Signal Processing (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치가 제공된다. 자기 공명 영상 처리 장치는, 프로세서 및 프로세서에 연결된 메모리를 포함하고, 프로세서는, 대상체에 대한 불완전하게 샘플링된 멀티 코일 k-공간을 획득하고, 미리 획득된 제1 사전 및 제1 사전을 이용하여 획득된 제2 사전을 이용하여, 멀티 코일 k-공간을 복원하여 대상체의 복원 영상을 획득할 수 있다.

Description

자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치{Method of Magnetic Resonance Image reconstruction and Apparatus for processing Magnetic Resonance Image}
본 발명은 의료 영상 처리에 관한 것으로, 더 구체적으로는 복수의 수신기 코일에 의해 획득된 불완전하게 샘플링된 k-공간 데이터로부터 대상체의 자기 공명(MR) 영상 복원을 하기 위한 방법 및 시스템과 관련된 발명이다.
자기 공명 영상(MRI)은 해부학적 구조와 생리학적 기능 모두를 시각화 시킬 수 있기 때문에 진단에 널리 이용되는 비침습적인 영상 기술이다.
MRI는 핵 자기 공명 원리에 기초한다. 교번 자기장(alternating magnetic field) 및 무선 주파수(RF) 신호는 인체 일부에서의 양성자 밀도 분포에 관한 정보를 얻기 위한 수소 원자 여기를 위해 사용된다. 여기된 수소 원자로부터의 신호 응답은 소위 "수신기 코일" 세트에 의해 수신된다. 이러한 수신기 코일은 무선 안테나(코일), RF 수신기, 증폭기, 검출기, 및 신호의 수신, 증폭, 동기화 및 디지털 형태로의 전환을 수행할 수 있는 다른 요소들의 조합일 수 있다. 현재, 자기 공명(MR) 스캐너는 상술한 여러 RF 수신기 코일을 포함할 수 있다.
복수의 RF 수신기 코일이 획득하는 것은 가시화할 대상체의 의학 영상이 아니라, 대상체의 공간적 이산 퓨리에 스펙트럼이다. 이러한 스펙트럼을 일반적으로 "k-공간"이라고 한다. 만일 하나의 RF 수신기 코일이 사용되는 경우, 하나의 k-공간이 획득된다. 만일 복수의 RF 수신기 코일이 사용되는 경우, MR 스캐너는 수신기 코일의 공간적인 위치 및 개별적인 매개변수에 따라 대상체로부터 각자의 신호를 획득하는, 서로 다른 RF 수신기 코일에 대응하는 몇몇의 k-공간을 제공한다. 서로 다른 수신기 코일의 k-공간은 이하에서 멀티 코일 k-공간(multi-coil k-space)이라 한다. 이 멀티 코일 k-공간은 완전히 또는 불완전하게 샘플링(언더 샘플링)될 수 있다. 불완전하게 샘플링 되는 것은 멀티 코일 k-공간의 일부가 취득 과정에서 생략된다는 것을 의미한다. k-공간을 불완전하게 샘플링하는 이유는 완전한 샘플링의 경우 스캔 절차에서 시간이 매우 소모되며, 스캔될 멀티 코일 k-공간의 상당한 부분을 생략하는 것으로 의학 영상 획득시간을 현저히 감소시킬 수 있기 때문이다.
그러나, 불완전하게 샘플링된 k-공간으로부터 자기 공명 영상을 복원한 예시들은 일반적으로 좋지 않은 영상 품질을 가지고 있으며, 강한 앨리어싱 아티팩트(aliasing artifacts)에 의해 손상된다. 따라서, 대상체의 복원 품질을 개선하기 위해서, 다양한 사후-처리(post-processing) 절차가 적용된다. 몇몇의 저자들은 복수의 앨리어싱된 측정치들로부터 직접 대상체를 복원하기 위해 RF 수신기 코일의 감도를 추정한다. 다른 방법으로, 멀티 코일 k-공간의 생략된 스펙트럼 요소를 통상적으로 "코일 감도 지도" 라고 하는 코일 간의 기능 관계를 이용하여 직접적으로 재계산하는 방법이 있다. 이 유형의 가장 대표적인 방법은 스매시(SMASH)와 그라파(GRAPPA)이다. 특허 문헌 US 8222900 B2 및 US 2014/0218026 A1에 k-공간 처리를 위해 GRAPPA가 적용된 다양한 버전이 있다. US 8427156 B2 및 US 7884604 B2에서, 코일 감도 지도가 내재적으로 사용되었다. 그러나 코일 감도 지도는 정확하게 추정되는 것이 매우 어려우며, 시간 소모적인 절차일 수 있다. 또 다른 공지기술인 증강된(augmented) 라그랑(Lagrangian) 기반 방법(예를 들어, WO 2013/067546 A1을 보라)은 다루기 어렵고, 계산이 복잡하다. 그러므로, 상기 언급된 방법들을 사용하는 것은 실제로 충분히 합리적이지 않을 수 있다.
또한, 희박(sparse) 근사치에 기초하는 경쟁적인 공지의 MRI 복원 기법들이 있다. 이러한 기법 중 하나는 압축 센싱(Compressed sensing, CS)이며, CS의 주된 아이디어는 실시간 신호들/영상들은 일반적으로 일부 도메인에서 희박하고(sparse), 미지수의 개수보다 상당히 적은 수의 측정 치들을 이용하여 실시간 신호들/영상들을 (허용 가능한 품질로) 정확하게 복구 할 수 있다는 것이다. 요점은 오버컴플리트(overcomplete) 기반으로 객체의 적절한 스파스(sparse) 표현을 찾는 것이다. Donoho [특허 출원 US 2006/0029279] 에 의해 시작된 획기적인 CS 기법은 MR 영상의 정확한 복원을 위해, 사용된 k-공간이 완전히 샘플링 될 필요가 없음을 입증한다. CS는 모든 측정에서 객체의 희박성을 활용하고, 큰 잠재력을 보여 오늘 날의 최첨단 기술이 되었다. 그러나, 최적화 문제에 있어서 일반적으로 너무 많은 변형 공식(variational formulation)으로 인해, 특정 희박화 변환(sparsifying transform)을 지정할 필요가 있다. 사용되는 기저와 희박화 변환(sparsifying transform)에 의해 구별되며, 스파스(sparse) 객체 추정을 기초로 하는 많은 신호/영상 복원 기법들이 있다. 특히, 특허 출원 공보 (WO 2002/031756 A1)에서, 객체 추정은 웨이블릿(wavelet) 기저를 통해 이루어지고, 특허 (US 8699773 B2)에서는 "유사성 클러스터(similarity clusters)"라고 부르는 영상 특성들에 의해, 그리고 특허 출원 (US 2010/0239143 A1)에서는 이전에 수신된 완전히 샘플링된 영상의 희박화된 데이터 세트를 사용하여 객체 추정이 이루어진다. 상기 언급된 특허 문헌들은 MR 영상의 특수성에 초점을 맞추지 않으며, 이러한 접근들은 다른 공지의 기술인 사전 기반(dictionary-based)의 스파스(sparse) 추정 방법과 대조적으로 MR 영상 특성에 적응되지 않는다.
GSENSE 및 관련된 응용에 따라, 공동 희박성(joint sparsity)은 대상체에서의 중복된(redundant) 정보를 활용한다. 그러나 복원 품질을 상당히 개선하기 위해서, 현대 CS 기법([R. Otazo 외 "고도로 가속화된 1차 통과 심장 관류(first-pass cardiac perfusion) MRI를 위한 CS 및 병렬 영상의 조합," Magn Reson Med. 64, 767-776 (2010)을 보라])이 사용된다. 특허 US 8587307 B2에서는, CS가 병렬 MRI(pMRI)와 SENSE에 따라 폴딩된 원본 영상들로부터 결과 영상을 복구하기 위해 사용된다. 특허 출원 (US 2013/0099786 A1)에서는, 공동 희박성(joint sparsity)이 pMRI의 관점에서 SENSE 또는 PARS 알고리즘에 의해 사용된다.
모든 CS 알고리즘의 중요한 포인트는 불완전하게 샘플링된 측정된 멀티 코일 k-공간의 초기치 설정이다. 대부분의 문헌들은 기본적 제로 충전(default zero-filling) 사례에 대한 결과를 개시한다. 다수의 특허문헌들 (예를 들어, WO 2014/075005 A1, WO 2012/144957 A1 and CN 103505207 A을 보라) 은 반복적인 분할 브레그만(Iterative Split Bregman) 알고리즘에 의한 MRI 복원에 대해 제시하고 있으나, 아무도 초기치 설정을 위해 GSENSE와 같은 시나리오에서 분할 브레그만을 사용하지 않는다. 이것은 공동 CS의 이점을 활용하게 하고, 보다 정확하고 적절한 값들로 측정된 멀티 코일 k-공간의 샘플링 되지 않은 k-공간의 위치들을 채울 수 있게 한다. 분할 브레그만 알고리즘의 공동 버전에 기초한 초기치 설정은 제안된 CS MRI 복원 알고리즘의 또 다른 키포인트이다.
병렬 영상의 발전과 함께, 수신기 코일들의 어레이들로부터의 객체의 중복된 정보 때문에 MRI 데이터 획득은 상당히 가속화되었다. 또한, 수신기 어레이들이 있는 MRI는 단일 코일이 있는 MRI 기술에 비해 높은 신호 대 잡음비를 제공할 수 있다. 동시에, 소자들의 증가는 더 많은 데이터세트를 야기하고, 이로써 특히 3D k-공간의 경우, 영상 복원에 있어서 더 많은 계산을 필요로 한다. 코일 압축 알고리즘 (예를 들어, WO 2012/123921 및 [T. Zhang 외 "데카르트 샘플링을 이용한 가속화된 이미징을 위한 코일 압축," Magn Reson Med. 69, 571-582 (2013)]을 보라)은 많은 RF 수신기 코일을 더 적은 수의 소위 가상의 수신기 코일로 데이터 압축을 함으로써 이 문제를 완화하는데 효과적이다. 코일 압축(CC)은 미가공 데이터(raw data)를 스펙트럼의 샘플링에 변화 없이 선형적으로 조합함으로써, 가상의 수신기 코일의 k-공간은, RF 수신기 코일로부터 획득된 k-공간의 가중된 합(weighted sum)으로써 표현된다. 가중치 값을 구하는 계산 방법은 현 출원의 범위를 벗어난다. 코일 압축(CC) 방법의 이점은 다수의 수신기 코일에 대한 k-공간 프로세싱을 가능하게 하는 것이다.
본 발명에 있어서, K-SVD는 (WO 2006/106508 A2, US 2012/0177128, 및 [M. Aharon, M. Elad, and A. Bruckstein, "K-SVD:스파스 표현을 위한 오버컴플리트 사전 설계를 위한 알고리즘," IEEE Trans. Image Process. 54, 4311-4322 (2006) ]을 보라) CS 영상 복원을 위한 사전(dictionary)을 학습하는 데 사용된다. 사전은 완전히 샘플링된 일련의 MR 영상들에 앞서서 미리 학습된다. 사전을 학습하는 절차는 본 출원의 범위를 벗어난다.
상기에 고려된 모든 방법들에 의해, 특허 출원 (US 2013/0099786 A1)이 제안된 공동의 다중-사전 CS MRI 복원(joint multi-dictionary CS MRI reconstruction) 방법에 가장 근접한 것으로 나타난다.
개시된 발명 및 개시된 발명의 실시예들은 불완전하게 샘플링된 멀티 코일 k-공간으로부터 MR 영상을 복원 하기 위한 방법 및 시스템에 관한 것이다. 데이터를 획득하는 총 시간은 샘플들 수에 의존하기 때문에, MR 획득 절차를 가속하기 위한 간단한 방법은 샘플링 된 주파수의 양을 줄이고, 불완전하게 샘플링된 k-공간 스펙트럼을 형성하는 MR 영상을 복원하는 것이다.
개시된 발명에 따른 방법에서, pMRI와 CS 기법의 장점들을 결합하는 것이 제안되었으나, 종래의 방법들과 달리, 균일하지 않은 크기의 패치들을 이용한 미리 학습(pre-learn)된 사전(dictionary)에 기초하여, 최첨단 CS MRI 복원의 최근 개발을 사용하는 것이 제안된다. 소위 단일의 합성 영상에 적용되는 사전 기반 스파스 근사(sparse approximation)는 멀티 코일 k-공간에서 측정된 샘플들의 복원과 번갈아 가며 나타날 수 있다. 또한, 상이한 주파수들의 다양하고, 섬세한 특징들을 표적으로 하는 미리 계산된 여러 사전들이 스파스 근사를 위해 적용될 수 있다.
멀티 코일 k-공간의 불완전한 샘플링은 이진 행렬로 표현된 미리 정의된 샘플링 방식(일반적으로 소위 "샘플링 마스크")에 따라 MRI 스캐너에 의해 수행될 수 있다. 샘플링 마스크의 "1"은 해당 공간 주파수(해당 위치에서의 k-공간 스펙트럼 지점)가 샘플링 되는 것을 의미하고, "0"은 각각의 주파수가 샘플링 되지 않는 것을 의미할 수 있다. 전술한 방식에 따라 샘플링한 결과가 불완전하게 샘플링된 k-공간 이다. 획득된 k-공간에서 불완전하게 샘플링된 값들은 몇몇의 초기 값들로 채워지며, 일반적으로 0으로 채워질 수 있다. 샘플링 마스크는 모든 RF 수신기 코일들에 사용된 것과 동일하며, 상기 방법은 서로 다른 수신기 코일들로부터 획득된 복수의 불완전하게 샘플링된 k-공간 스펙트럼인, 멀티 코일 k-공간을 입력으로 한다. 샘플링 주파수들의 양이 적을 수록, 멀티 코일 k-공간의 언더 샘플링 율은 더 높아질 수 있다.
개시된 실시예에 따른 발명의 목적은 MR 스캐너로부터 획득된 불완전하게 샘플링된 멀티 코일 k-공간의 개시된 처리를 통하여, 허용 가능한 품질을 갖는 대상체의 MR 영상을 제공하는 데에 있다.
개시된 발명에 따른 실시예에서, CS 방식은, MR 영상 복원을 위해, 특히, 비균일한 크기의 패치를 갖는 다양한 미리 계산된 사전들에 기초한 스파스 객체 근사를 하기 위해 사용될 수 있다. 일반적으로 개시된 실시예에 따른 발명은 이하의 내용에 주요한 특징이 있다. 사전 기반 CS 스파스 근사와 병렬 MRI의 공동 희박성(joint sparsity)의 조합은 서로 다른 수신기 코일에서의 효과적인 데이터 합성이 가능하게 할 수 있다. 이는 측정 노이즈를 억제하고, 비교적 빠른 복원을 가능하게 하며, 강한 앨리어싱 아티팩트를 필터링 해낼 수 있다. 더욱이, 사전 기반 CS는 사전(dictionary) 학습과정을 통해 MRI의 세부 사항에 적응되기 때문에 적절한 희박화 변환(sparsifying transform)을 선택할 필요가 없으므로, 특정 고정 기반 함수들(예를 들어 웨이블릿)보다 더 유연할 수 있다. 제안된 새로운 공동 분할 브레그만(Joint Split Bregman) 및 다양한 사전들에 의한 스파스 객체 근사에 의한 데이터 초기치 설정의 효율성은 수렴 속도의 상당한 증가를 야기하며, 고도의 언더 샘플링률 케이스인 경우에도 영상의 세밀한 복원이 가능하게 할 수 있다. 멀티 코일 k-공간의 개개의 주파수 대역 처리를 위해 사용되는 다중 대역 분해(multiband decomposition)는 주파수 대역 내의 큰 희박성(sparsity)으로 인한 영상 복원 품질을 더욱 개선할 수 있다.
상술한 기술적 과제를 달성하기 위한 기술적 수단으로서, 본 개시의 제1 측면은, 프로세서 및 프로세서에 연결된 메모리를 포함하고, 프로세서는, 대상체에 대한 불완전하게 샘플링된 멀티 코일 k-공간을 획득하고, 미리 획득된 제1 사전 및 제1 사전을 이용하여 미리 획득된 제2 사전을 이용하여, 멀티 코일 k-공간을 복원하여 대상체의 복원 영상을 획득하는, 자기 공명 영상 처리 장치를 제공할 수 있다.
일 실시예에 따른 제2 사전은 미리 획득된 대상체의 영상과 그 영상을 제1 사전에 기초하여 근사한 것의 잔차를 기초하여 획득될 수 있다.
일 실시예에 따른 근사는 스파스 근사일 수 있다.
일 실시예에 따른, 프로세서는, 불완전하게 샘플링된 멀티 코일 k-공간에 기초하여 제1 영상을 획득하고, 제1 영상의 균일하지 않은 크기의 패치들, 제1 사전 및 제2 사전에 기초하여 멀티 코일 k-공간을 복원하여 복원 영상을 획득할 수 있다.
일 실시예에 따른, 프로세서는, 불완전하게 샘플링된 멀티 코일 k-공간에 기초하여 획득된 서로 다른 주파수 대역에 대응되는 패치들, 제1 사전 및 제2 사전에 기초하여 멀티 코일 k-공간을 복원하여 복원 영상을 획득할 수 있다.
일 실시예에 따른 자기 공명 영상 처리 장치는 오퍼레이터 콘솔을 더 포함하고, 오퍼레이터 콘솔은 복원 영상을 시각화 하기 위한 표시부를 포함할 수 있다.
일 실시예에 따른 오퍼레이터 콘솔은 제어부를 더 포함하고, 제어부는 자기 공명 영상 처리 장치를 제어하도록 구성될 수 있다.
일 실시예에 따른 자기 공명 영상 처리 장치는 불완전하게 샘플링된 멀티 코일 k-공간 데이터를 획득하도록 구성된 MR 스캐너를 더 포함할 수 있다.
일 실시예에 따른 프로세서는, 획득한 불완전하게 샘플링된 멀티코일 k-공간을 복수 개의 주파수 대역으로 분할하고, 복수 개의 주파수 대역으로 분할된 불완전하게 샘플링된 멀티 코일 k-공간들을 각각 제1 사전 및 제2 사전을 이용하여 복원함으로써 복원 영상을 획득할 수 있다.
일 실시예에 따른 프로세서는, 불완전하게 샘플링된 멀티 코일 k-공간을 전처리 하고, 전처리된 멀티 코일 k-공간을, 제1 사전 및 제2 사전을 이용하여 복원할 수 있다.
일 실시예에 따른 프로세서는, 불완전하게 샘플링된 멀티 코일 k-공간의 샘플링 되지 않은 위치들을 초기 추측(initial guess)으로 채우고, 불완전하게 샘플링된 멀티 코일 k-공간을 코일 압축(coil compressing)함으로써, 불완전하게 샘플링된 멀티 코일 k-공간을 전처리 할 수 있다.
또한, 본 개시의 제2 측면은, 대상체에 대한 불완전하게 샘플링된 멀티 코일 k-공간을 획득하는 단계; 및 미리 획득된 제1 사전 및 제1 사전을 이용하여 미리 획득된 제2 사전을 이용하여, 멀티 코일 k-공간을 복원하여 대상체의 복원 영상을 획득하는 단계를 포함하는 자기 공명 영상을 복원하기 위한 방법을 제공할 수 있다.
또한, 상기 불완전하게 샘플링된 멀티 코일 k-공간을 전처리 하는 단계는 상기 불완전하게 샘플링된 멀티 코일 k-공간의 샘플링 되지 않은 위치들을 초기 추측(initial guess)으로 채우는 단계; 및 상기 불완전하게 샘플링된 멀티 코일 k-공간을 코일 압축(coil compressing)하는 단계를 포함할 수 있다.
또한, 상기 멀티 코일 k-공간의 샘플링 되지 않은 위치들을 초기 추측으로 채우는 단계는, 상기 샘플링 되지 않은 위치들을 상수, 랜덤 복소수 값들 및 공동 분할 브레그만 반복 알고리즘(joint Split Bregman iterative algorithm) 중 적어도 하나에 의해 얻어진 값들로 채우는 단계를 포함할 수 있다.
또한, 상기 전처리 된 멀티 코일 k-공간을 다중 대역 공동 CS 복원을 이용하여 복원하는 단계는, 상기 멀티 코일 k-공간을 다중 대역으로 분해하여 서로 다른 주파수 대역들에 대응하는 복수의 멀티 코일 k-공간들을 획득하는 단계; 상기 서로 다른 주파수 대역들에 대응하는 복수의 멀티 코일 k-공간들에 기초하여 코일 감도 지도들을 획득하는 단계; 상기 서로 다른 주파수 대역들에 대응하는 복수의 멀티 코일 k-공간들의 다중 사전 공동 CS 복원(multi-dictionary joint compressed sensing reconstruction)을 수행하는 단계; 상기 복원된 복수의 멀티 코일 k-공간들에 기초하여 최종 멀티 코일 k-공간을 획득하는 단계; 및 상기 획득된 최종 멀티 코일 k-공간으로부터 자기 공명 영상을 획득하는 단계를 포함할 수 있다.
또한, 상기 서로 다른 주파수 대역들에 대응하는 복수의 멀티 코일 k-공간들의 상기 다중 사전 공동 CS 복원을 수행하는 단계는, 상기 서로 다른 주파수 대역들에 대응하는 복수의 멀티 코일 k-공간들로부터 합성된 영상을 획득하는 단계; 다중 사전 스파스(sparse) 근사에 기초하여, 상기 합성된 영상으로부터 상기 합성된 영상의 근사치를 획득하는 단계; 상기 합성된 영상의 근사치로부터 멀티 코일 k-공간 데이터를 계산하는 단계; 및 상기 계산된 멀티 코일 k-공간 데이터를 복원하는 단계를 포함할 수 있다.
또한, 상기 복수의 멀티 코일 k-공간들로부터 상기 합성된 영상을 획득하는 단계는, 상기 복수의 멀티 코일 k-공간들을 코일-와이즈(coil-wise) 퓨리에 역변환하여 복수 개의 코일 영상들을 획득하는 단계; 및 상기 복수 개의 코일 영상들을 상기 코일 감도 지도들의 각 공간 위치에 투영(projecting)하는 단계를 포함할 수 있다.
또한, 상기 합성된 영상으로부터 상기 합성된 영상의 근사치를 획득하는 단계는, 패치 기반 다중 사전 스파스 코딩을 수행하여 스파스 코드를 획득하는 단계; 및 상기 획득된 스파스 코드에 기초하여 상기 합성된 영상을 근사하여, 상기 합성된 영상의 근사치를 획득하는 단계를 포함할 수 있다.
또한, 상기 합성된 영상의 근사치로부터 상기 멀티 코일 k-공간 데이터를 계산하는 단계는, 상기 코일 감도 지도들을 리스케일링(rescaling)하는 단계; 및 상기 합성된 영상의 근사치를 복수 개의 코일 영상들로 분해하고, 상기 복수 개의 코일 영상들을 퓨리에 변환하는 단계를 포함할 수 있다.
또한, 상기 다중 사전 스파스 근사에 기초하여 상기 합성된 영상으로부터 상기 합성된 영상의 근사치를 획득하는 단계는, 상기 합성된 영상으로부터 중복된(overlapping) 패치들을 추출하는 단계; 상기 추출된 패치들의 사전 기반 스파스 근사를 수행하는 단계; 및 상기 사전 기반 스파스 근사가 수행된 패치들을 상기 중복된 패치들과 조합하여 상기 합성된 영상의 근사치를 획득하는 단계를 포함할 수 있다.
또한, 상기 추출된 패치들의 사전 기반 스파스 근사를 수행하는 단계는, 단일 사전, 원본 지향적 사전 및 잔여 지향적 사전 중 적어도 하나를 이용하여 수행될 수 있다.
도 1은 예시적인 실시예에 따른 자기 공명 영상 처리 장치 도시한다.
도 2는 예시적인 실시예에 따른 공동 CS MRI 복원 방법을 나타내는 흐름도이다.
도 3은 예시적인 실시예에 따른 공동 CS MRI 복원 방법을 나타내는 흐름도이다.
도 4는 예시적인 실시예에 따른 공동 분할 브레그만 반복 알고리즘에 기초한 데이터 초기치 설정에 대한 흐름도이다.
도 5는 예시적인 실시예에 따른 사전 기반 공동 CS 복원 방법에 대해 도시한다.
도 6은 예시적인 실시예에 따른 패치들의 단일 사전 스파스 근사에 대해 도시한다.
도 7은 예시적인 실시예에 따른 패치들의 다중 사전 스파스 근사에 대해 도시한다.
도 8은 예시적인 실시예에 따른 사전 기반 스파스 근사에 대해 도시한다.
도 9는 예시적인 실시예에 따른 원본 및 잔차 지향적인 사전들의 학습 과정을 도시한다.
도 10은 예시적인 실시예에 따른 공동 CS MRI 복원 시스템을 도시한다.
도 11은 예시적인 실시예에 따른 다중 사전 CS MRI 복원 시스템을 도시한다.
도 12는 예시적인 실시예에 따른 원본 지향적 및 잔차 지향적인 사전들을 학습하는 시스템을 도시한다.
도 13은 예시적인 실시예에 따른 자기 공명 영상 처리 시스템을 도시한다.
도 14는 예시적인 실시예에 따른 공동 분할 브레그만 복원을 위한 발명의 시스템을 도시한다.
도 15는 예시적인 실시예에 따른 사용자 시나리오 1을 도시한다.
도 16은 예시적인 실시예에 따른 사용자 시나리오 2를 도시한다.
도 17은 예시적인 실시예에 따른 사용자 시나리오 3을 도시한다.
도 18은 MRI 시스템의 개략도이다.
명세서 전체에 걸쳐 동일 참조 부호는 동일 구성요소를 지칭한다. 본 명세서가 실시예들의 모든 요소들을 설명하는 것은 아니며, 본 발명이 속하는 기술분야에서 일반적인 내용 또는 실시예들 간에 중복되는 내용은 생략한다. 명세서에서 사용되는 ‘부’(part, portion)라는 용어는 소프트웨어 또는 하드웨어로 구현될 수 있으며, 실시예들에 따라 복수의 ‘부’가 하나의 요소(unit, element)로 구현되거나, 하나의 ‘부’가 복수의 요소들을 포함하는 것도 가능하다. 이하 첨부된 도면들을 참고하여 본 발명의 작용 원리 및 실시예들에 대해 설명한다.
본 명세서에서 영상은 자기 공명 영상(MRI) 장치, 컴퓨터 단층 촬영(CT) 장치, 초음파 촬영 장치, 또는 엑스레이 촬영 장치 등의 의료 영상 장치에 의해 획득된 의료 영상을 포함할 수 있다.
본 명세서에서 ‘대상체(object)’는 촬영의 대상이 되는 것으로서, 사람, 동물, 또는 그 일부를 포함할 수 있다. 예를 들어, 대상체는 신체의 일부(장기 또는 기관 등; organ) 또는 팬텀(phantom) 등을 포함할 수 있다.
MRI 시스템은 자기 공명(magnetic resonance, MR) 신호를 획득하고, 획득된 자기 공명 신호를 영상으로 재구성한다. 자기 공명 신호는 대상체로부터 방사되는 RF 신호를 의미한다.
MRI 시스템은 주자석이 정자장(static magnetic field)을 형성하여, 정자장 속에 위치한 대상체의 특정 원자핵의 자기 쌍극자 모멘트 방향을 정자장 방향으로 정렬시킨다. 경사자장 코일은 정자장에 경사 신호를 인가하여, 경사자장을 형성시켜, 대상체의 부위 별로 공명 주파수를 다르게 유도할 수 있다.
RF 코일은 영상 획득을 원하는 부위의 공명 주파수에 맞추어 RF 신호를 조사할 수 있다. 또한, RF 코일은 경사자장이 형성됨에 따라, 대상체의 여러 부위로부터 방사되는 서로 다른 공명 주파수의 MR 신호들을 수신할 수 있다. 이러한 단계를 통해 MRI 시스템은 영상 복원 기법을 이용하여 MR 신호로부터 영상을 획득한다.
스캐닝 되는 대상체의 영상을 형성하기 위해, 여기된 수소 원자들로부터의 신호 응답은 RF 수신기 코일들의 세트에 의해 수신될 수 있다. 무선 주파수(radio-frequency) 수신기 코일은 무선 안테나, RF 수신기, 증폭기, 검출기, 및 신호의 수신, 증폭, 동기화 및, 디지털 형태로 전환을 실현시키는 몇몇의 다른 요소들의 조합일 수 있으며, 수신기 코일들의 세트가 동시에 객체의 응답 신호를 얻기 위해 사용된다. 수신기 코일들은 서로 다른 공간적 위치와 매개변수 설정을 가지고 있다. 이러한 병렬적인 데이터 수집은 공동 CS 복원(joint CS reconstruction) 시 대상체에 대한 매우 중복된 정보를 초래할 수 있다. 획득 시간은 k-공간들의 불완전한 샘플링에 의해 감소될 수 있다. k-공간의 불완전한 샘플링은 모든 수신기 코일들에 대해 동일한 소정의 샘플링 마스크에 따라 수행될 수 있다. 샘플링 되지 않은 위치들의 값들은 일부 초기값들(일반적으로 0)로 채워질 수 있다. 이러한 불완전하게 샘플링된 멀티 코일 k-공간은 개시된 실시예에 따른 CS MRI 복원 알고리즘들의 입력으로써 사용된다.
개시된 실시예에 따른 발명에서, 수신기 코일 안테나의 기술적인 사양과 구성에는 제한이 없다.
도 1은 예시적인 실시예에 따른 공동 CS MRI 복원을 위한 자기 공명 영상 장치를 도시한다.
일 실시예에 따르면, 자기 공명 영상 장치는 도 1에 도시된 바와 같이 복원 서버(100)로 구현될 수 있다. 복원 서버(100)는 대상체를 자기 공명 영상 촬영하여 획득한 MR 신호를 처리하기 위한 장치일 수 있다.
복원 서버(100)는 MR 스캐너(도시되지 않음)로부터 MR 신호를 수신할 수 있다. 또한, 복원 서버(100)는 MR 스캐너로부터 MR 신호에 대응되는 멀티 코일 k-공간 데이터를 수신할 수 있다. 복원 서버(100)는 MR 신호에 대응되는 멀티 코일 k-공간 데이터를 이용하여 자기 공명 영상을 복원하는 서버 장치를 포함할 수 있다. 여기서, 서버 장치는 환자가 자기 공명 영상 촬영을 진행하는 병원 내부의 서버 장치일 수도 있고, 병원 외부에 위치된 의료용 서버 장치가 될 수도 있다.
복원 서버(100)는 프로세서(110)와 메모리(120)를 포함할 수 있다. 다만 이는 예시적인 것으로, 복원 서버(100)는 도 1에 도시된 구성보다 더 많은 구성을 포함할 수 있다.
프로세서(110)는 대상체에 대한 불완전하게 샘플링된 멀티 코일 k-공간을 획득할 수 있다. 일 실시예에 따르면, 프로세서(110)는 불완전하게 샘플링된 멀티 코일 k-공간을 MR 스캐너로부터 수신할 수 있다.
일 실시예에 따르면, 프로세서(110)는 복수개의 사전을 이용하여 멀티 코일 k-공간을 복원하여 대상체의 복원 영상을 획득할 수 있다. 사전은 계산을 단순화하기 위해 미리 계산된 데이터를 의미할 수 있다. 예를 들어 사전은, 테이블 및 룩업 테이블등을 포함할 수 있다. 이 때, 프로세서(110)는 메모리(120)에 저장된 영상 복원 프로그램을 작동할 수 있다. 또한, 프로세서(110)는 불완전하게 샘플링된 자기 공명 영상 데이터를 복원하기 위한 모듈을 포함할 수 있다.
또한, 프로세서(110)는 미리 획득된 제1 사전 및 제1 사전을 이용하여 획득된 제2 사전을 이용하여, 멀티 코일 k-공간을 복원하여 대상체의 복원 영상을 획득할 수 있다. 여기서 복원 영상은 최종적으로 얻고자하는 대상체에 대한 자기 공명 영상을 의미할 수 있다. 일 실시예에 따른 제2 사전은 미리 획득된 대상체의 영상과 미리 획득된 영상을 제1 사전에 기초하여 스파스 근사한 것의 잔차를 기초하여 획득될 수 있다.
또한, 프로세서(110)는 불완전하게 샘플링된 멀티 코일 k-공간에 기초하여 제1 영상을 획득하고, 제1 영상의 균일하지 않은 크기의 패치들, 제1 사전 및 제2 사전에 기초하여 멀티 코일 k-공간을 복원하여 대상체의 복원 영상을 획득할 수 있다. 여기서 제1 영상은 서로 다른 주파수 대역들에 대응하는 복수의 멀티 코일 K-공간들로부터 획득된 영상들을 합성한 합성 영상일 수 있다.
또한, 프로세서(110)는 불완전하게 샘플링된 멀티 코일 k-공간에 기초하여 획득된 서로 다른 주파수 대역에 대응되는 패치들, 제1 사전 및 제2 사전에 기초하여 멀티 코일 k-공간을 복원하여 대상체의 복원 영상을 획득할 수 있다.
또한, 프로세서(110)는 획득한 불완전하게 샘플링된 멀티코일 k-공간을 복수 개의 주파수 대역으로 분할하고, 복수 개의 주파수 대역으로 분할된 불완전하게 샘플링된 멀티 코일 k-공간들을 각각 제1 사전 및 제2 사전을 이용하여 복원함으로써 복원 영상을 획득할 수 있다.
또한, 프로세서(110)는 불완전하게 샘플링된 멀티 코일 k-공간을 전처리 하고, 전처리된 멀티 코일 k-공간을, 제1 사전 및 제2 사전을 이용하여 복원할 수 있다.
또한, 프로세서(110)는 불완전하게 샘플링된 멀티 코일 k-공간의 샘플링 되지 않은 위치들을 초기 추측(initial guess)으로 채우고, 불완전하게 샘플링된 멀티 코일 k-공간을 코일 압축(coil compressing)함으로써, 불완전하게 샘플링된 멀티 코일 k-공간을 전처리 할 수 있다.
예를 들어, 프로세서(110)가 초기치 추측으로 샘플링 되지 않은 위치들을 채우는 값들은, 상수, 랜덤 복소수 값들 또는 공동 브레그만 반복 알고리즘에 의해 획득된 값들을 포함할 수 있다. 예를 들어, 초기치 추측 값은 0 일 수 있다. 예를 들어 프로세서(110)가 코일 압축하는 것은, 다수의 RF 수신기 코일을 더 적은 수의 가상의 수신기 코일로 데이터 압축을 하는 것을 포함할 수 있다. 일 실시예에 따르면, 프로세서(110)는 전처리된 멀티 코일 k-공간을 복원 알고리즘을 이용하여, 복원할 수 있다. 예를 들어, 프로세서(110)는 다중 대역 공동 CS 복원을 이용하여 전처리된 멀티 코일 k-공간을 복원할 수 있다. 프로세서(110)가 다중 대역 공동 CS 복원을 이용하여, 전처리된 멀티 코일 k-공간을 복원하는 과정에는, 불완전하게 샘플링된 멀티 코일 k-공간을 복수 개의 주파수 대역(예를 들어, 낮은 주파수 대역, 높은 주파수 대역)으로 분할하는 단계가 포함될 수 있다.
프로세서(110)는 복수 개의 주파수 대역으로 분할된 멀티 코일 k-공간을 각각 다중 사전 공동 CS 복원할 수 있다. 또한, 프로세서(110)는 복원된 복수 개의 멀티 코일 k-공간 결과를 합성함으로써, 복원된 자기 공명 영상을 획득할 수 있다. 그러나, 복원 서버(100)가 다중 대역 공동 CS 복원을 이용하여 자기 공명 영상을 복원하는 과정이 이에 한정되는 것은 아니며, 더 많은 단계에 의해 수행될 수도 있다.
일 실시예에 따른, 메모리(120)는 불완전하게 샘플링된 멀티 코일 k-공간을 획득하기 전에 미리 획득되는 복수개의 사전들을 저장할 수 있다. 예를 들어, 메모리(120)는 미리 획득된 제1 사전 및 제1 사전을 이용하여 획득된 제2 사전을 저장할 수 있다. 또한, 메모리(120)는 자기 공명 영상을 복원하기 위해 필요한 펄스 시퀀스, 코일 감도 지도 등을 저장할 수 있다.
자기 공명 영상 장치는 오퍼레이터 콘솔(도시되지 않음) 더 포함할 수 있다. 오퍼레이터 콘솔은 자기 공명 영상 처리 시스템을 제어하도록 구성되는 제어부(도시되지 않음) 및 복원된 자기 공명 영상을 디스플레이하기 위한 표시부(도시되지 않음)를 포함할 수 있다.
또한, 자기 공명 영상 장치는 k-공간 획득을 위해 구성된 적어도 하나의 수신기 코일을 포함하는 자기 공명 스캐너(도시되지 않음)를 더 포함할 수 있다.
도 2는 예시적인 실시예에 따른 자기 공명 영상을 복원하는 방법을 나타내는 흐름도이다.
도 2의 자기 공명 영상을 복원하는 방법의 각 단계들은 예를 들어, 도 1에서 설명한 자기 공명 영상 복원 시스템에서 수행될 수 있다.
단계 S110에서 자기 공명 영상 복원을 위한 시스템은, 대상체에 대한 불완전하게 샘플링된 멀티 코일 k-공간을 획득할 수 있다(S110).
또한, 단계 S120에서 자기 공명 영상 복원을 위한 시스템은, 제1 사전과 제2 사전을 이용하여 대상체의 복원 영상을 획득할 수 있다(S120). 구체적으로, 자기 공명 영상 복원을 위한 시스템은 미리 획득된 제1 사전 및 제1 사전을 이용하여 획득된 제2 사전을 이용하여, 멀티 코일 k-공간을 복원하여 대상체의 복원 영상을 획득할 수 있다.
도 3a 및 도 3b는 예시적인 실시예에 따른 공동 CS MRI 복원 방법을 나타내는 흐름도이다.
일 실시예에 따르면, 공동 CS MRI 복원 방법은 도 3a에 도시된 바와 같이, 3가지의 주요 단계들을 포함할 수 있다. 또한, 공동 CS MRI 복원 방법은 도 3b에 도시된 바와 같은 세부 단계들을 포함할 수 있다.
일 실시예에 따르면, 공동 CS MRI 복원 방법은, 복원될 불완전하게 샘플링된 멀티 코일 K- 공간을 획득하는 단계를 포함 할 수 있다(S210). 불완전하게 샘플링된 멀티 코일 k-공간은 MR 스캐너(도시되지 않음)로부터 수신될 수 있다. 불완전하게 샘플링된 멀티 코일 k-공간은 소정의 샘플링 방식(샘플링 마스크)에 따라 획득될 수 있다.
일 실시예에 따르면, 공동 CS MRI 복원 방법은, 획득된 불완전하게 샘플링된 멀티 코일 k-공간을 전처리 하는 단계를 포함할 수 있다(S220, S309). 전처리 단계(S220, S309)에서, 샘플링 되지 않은 위치들을 일부 초기값들로 채우기 위한 멀티 코일 k-공간의 초기치 설정(S302)이 수행될 수 있다. 초기치 설정에는, 0또는 랜덤 복소수 값들을 포함하는 다른 유형의 초기치 설정들이 적용될 수 있다. 일 실시예에 따르면, 공동 분할 브레그만 반복 알고리즘(Joint split Bregman iterative algorithm)이 이러한 초기치 설정 절차에 적용될 수 있다. 일 실시예에 따르면, 전처리 단계(S220, S309)에서, 초기치가 설정된 멀티 코일 k-공간은, 총 k-공간들의 수를 감소시키기 위해, CC(Coil compression, 코일 압축)에 의해 압축될 수 있다(S303). 코일 압축으로 생성된 k-공간들은 RF 수신기 코일들 대신에 사용되는 소위 "가상 코일"과 관련될 수 있다.
일 실시예에 따르면, 공동 CS MRI 복원 방법은, 전처리된 멀티 코일 k-공간을 다중 대역 공동 CS 복원을 이용하여 복원하는 단계를 포함할 수 있다(S230, S310). 멀티 코일 k-공간은, 상기 생성된 각각의 멀티 코일 k-공간을 각각의 주파수 대역에 대응하는 방식으로, 여러 개의 주파수 대역들로 분해될 수 있다(S304). 또한 상기 코일 감도 지도들(Coil sensitivity maps)이 완전히 샘플링된 k-공간 데이터의 일부에 기초하여 획득될 수 있다(S305). 또한, 각각의 멀티 코일 k-공간은 각각의 주파수 대역 내의 샘플링 되지 않은 주파수들을 복원하기 위해, 상기 다중 사전 공동 CS 알고리즘을 사용하여 복원될 수 있다(S306). 또한, 서로 다른 주파수 대역들에 대응되는 모든 멀티 코일 k-공간들은 최종의 멀티 코일 k-공간을 형성하기 위해 합성될 수 있다(S307). 또한, 멀티 코일 k-공간의 결과물로부터 결과 MR 영상이 계산될 수 있다(S308).
아래에서, 사용된 CS 기술 및 병렬 MRI 복원의 기초적인 내용이 기술된다.
SENSE와 비슷한 반복 MRI 복원
병렬 MRI 알고리즘들은 측정된 k-공간들을 합성하기 위해 활용되는 코일 감도 지도들과 불가분하게 관련된다. 측정될 대상체의 양성자 밀도 π 와 i번째 수신기 코일에 의해 획득된 대상체의 불완전하게 샘플링된 퓨리에 스펙트럼 u i 간의 물리적인 관계는 아래와 같이 표현될 수 있다.
Figure pat00001
r 은 공간 좌표를 나타내며,
Figure pat00002
는 불완전하게 샘플링된 퓨리에 변환, s i (r)은 공간 좌표 r 에서의 i번째 수신기 코일의 코일 감도 지도 값을 나타낸다. 방정식 (1)은 대상체 π가 매개변수 설정과 수신기 코일들의 위치에 의한 다수의 k-공간의 차이를 생성하는 방법에 대해 설명한다.
기호°가 엘리먼트 와이즈 하다마르 곱(element-wise Hadamard product) 나타내는 것으로 하면, 방정식 (1)은 벡터의 형태로 아래와 같이 재작성될 수 있다.
Figure pat00003
대상체의 관심 대상(양성자 밀도)를 나타내는 것은 아래와 같이 편리해진다.
Figure pat00004
여기서 위첨자 H 는 헤르미이트 공액(Hermitian conjugate)을 나타내며,
Figure pat00005
은 불완전하게 샘플링된 퓨리에 역변환을 의미한다. 코일 감도 지도들의 복소수 벡터들과 k-공간들이 사용되는 것을 주목하라. 굵은 활자체의 소문자는 열(column)-벡터를 나타낸다.
방정식 (1) 및 (2)는 대상체를 멀티 코일 k-공간으로 분할하는 것을 묘사한다. 즉, π는 개별 코일 감도 지도들 s i 을 통하여 다양한 수신기 코일들에 대응되는 여러 개의 영상들 p i 로 분할되며, 이어서 상기 영상들 p i 는 퓨리에 도메인으로 전송된다. 영상 p i 는 일반적으로 i번째 "코일 영상"으로 언급된다. 방정식 (3)은 대응되는 k-공간들로부터 획득된 코일 영상들의 합성 연산을 정의한다. 상기 언급된 연산들(분할 및 합성)은 제안된 공동 CS MRI 복원과 GSENSE를 통합한다.
공동 분할 브레그만을 이용한 초기치 설정
개시된 실시예에 따른 방법에서, 발명자들은 효율적인 데이터 초기치 설정을 위해 분할 브레그먼 반복 알고리즘(Split Bregman iterative algorithm)의 새로운 변형을 제안한다. [T. Goldstein, and S. Osher, "l1-정규화된 문제를 위한 분할 프레그만 방법(The Split Bregman method for l1-regularized problems)", SIAM J. on Imag. Sciences 2, 323-343 (2009) ]에 따르면, 공동 CS 복원 문제는 다양한 수신기 코일들에 의해 획득되는 측정치들을 위해 디노이징(denoising) 문제를 추구하는 기초로서 아래와 같이 공식화 될 수 있다.
Figure pat00006
이것은 볼록 평가 함수(convex criterion function)의 최소화를 야기한다. 여기서
Figure pat00007
는 정규화 매개변수이고, Ψ은 희박화 변환(sparsifying transform)이며(전체 변환), l1-놈(norm)
Figure pat00008
는 벡터 요소들(items)의 절대값들의 합으로서 정의된다. 방정식 (4)에서, l1- 및 l2-놈은 [Y. Wang외, "전체 변환 정규화로 영상 블러링을 위한 고속 알고리즘", CAAM Technical Report (2007)]에 따라 아래와 같이 분리되며,
Figure pat00009
여기서 d는 스파스 객체 근사치이다. 이것은 아래의 공동 분할 브레그만 알고리즘을 따른다.
Figure pat00010
Figure pat00011
Figure pat00012
여기서 b k 는 k번째 반복에서의 보조 브레그만 매개변수(auxiliary Bregman parameters)의 벡터를 나타낸다. 단일 수신기 코일을 위한 일반적인 분할 브레그만 알고리즘과 단일의 수신기 코일의 측정치들과는 대조적으로, 방정식(6.1)의 공동 공식화(joint formulation)는 병렬 MRI 근사법의 모든 장점들을 가질 수 있다.
도 4는 예시적인 실시예에 따른 공동 분할 브레그만 반복 알고리즘에 기초한 데이터 초기치 설정에 대한 흐름도이다.
[T. Goldstein, and S. Osher, "l1-정규화된 문제를 위한 분할 프레그만 방법(The Split Bregman method for l1-regularized problems)", SIAM J. on Imag. Sciences 2, 323-343 (2009)] 및 GSENSE에 따르면, 공동 분할 브레그만 알고리즘 (6.1)-(6.3)의 구현은 두 개의 중첩된 루프로 구성될 수 있다. 첫 번째 단계(S401)에서, 스파스 객체 근사의 보조 변수들 d 0와 브레그만 매개변수 b 0는 제로 벡터들과 함께 초기화될 수 있다. 내부 루프 동안, 객체 π, 객체 π 의 스파스 근사 d, 및 보조 매개변수 벡터 b는 재계산될 수 있다. 단계(S404)에서, 객체 π는 추정된 영상들의 합성과 함께 서로 다른 수신기 코일들로부터의 k-공간들의 퓨리에 역변환에 의해 방정식(S3)에 따라 추정될 수 있다. 그 후, 방정식 (6.1) 및 (6.2)의 최소화에 의한 객체 재계산 (S405) 및 객체 근사치 d의 갱신 (S406)이 각각 수행될 수 있다. l1놈에 대한 최적화 (S406)는 1/2μ와 같은 임계치와 소프트 임계치 설정을 통해 실현될 수 있다. 단계 (S407)에서, 브레그만 매개변수 벡터 b는 방정식 (6.3)에 따라 갱신될 수 있다. 마지막으로, 생성된 객체는 방정식 (2)에 따라, 코일 영상들에 대응되는 여러 개의 k-공간들로 표현될 수 있다. 내부 루프는 정지 조건 (S403)이 충족될 때 종료되며, 그 후에 내부 루프 독립 변수들은 서로 다른 수신기 코일들의 측정된 k-공간들과 객체 복원에 맞게 갱신될 수 있다(S409).
외부 루프의 조건(S402)이 만족되면, 복원된 객체 추정치는 복원된 퓨리에 스펙트럼 에서 측정된 k-공간 샘플들을 복원(S410)하여 갱신될 수 있다. 이 절차는 공동 분할 알고리즘이 주어진 k-공간 측정치들에게 변화를 주지 않고, 공동 CS 복원을 위한 데이터를 준비하면서, 샘플링 되지 않은 주파수 하에서만 동작하도록 할 수 있다. 만일 상기 연산이 종료되면, 공동 분할 브레그만에 기초하여 초기치 설정이 완료될 수 있다.
사전 기반 스파스 객체 근사
개시된 실시예에 따른 방법의 주요 부분은 미리 계산된 사전을 이용하여 입력 체적 신호(input volumetric signal)의 스파스 근사에 기초하는 것이다. 이 단계에서 다음의 최적화 문제는 아래 식으로 해결될 수 있다.
Figure pat00013
방정식 (4)에서의 l1놈을 대신하여, 여기서는 벡터의 0이 아닌 요소들로서 정의되는 l0-norm
Figure pat00014
이 사용된다. 사전 기반 접근법(approach)에서, 희박화 변환(sparsifying transform)은 사전의 기본 요소들("원자들(atoms)"이라고 함)로 신호를 스파스 분해하는 것에 의해 나타난다. 개시된 실시예에 따른 방법에서, 패치 기반 스파스 분해가 적용된다: 체적 대상체의 비교적 작고 일정하지 않은 크기의 일부가 적절한 사전의 원자들의 선형 조합에 의해 벡터화되고 근사 될 수 있다. 이 선형 조합의 가중치들은 "스파스 코드(sparse codes)"라고 불린다. 패치의 차원수는 패치의 벡터화로 인해 문제 되지 않는다: 2D 또는 3D 패치는 열 벡터(column-vector)로 표현될 수 있다. 만일 패치의 원본 차원수가 원자들의 차원수에 대응된다면, 벡터화된 형태에서 모든 원본 패치(2D 또는 3D)는 사전의 벡터화된 원자들에 의해 근사화 될 수 있다. 수학적으로, 다음의 최적화 문제는 아래 식으로 해결될 수 있다.
Figure pat00015
여기서 D는 사전 변환 행렬이고, δ는 정규화 매개변수이다. 방정식 (8)의 열 벡터 π j 는 체적 MR 대상체의 합성된 영상의 패치이고, 대응하는 스파스 코드 z j 의 j번째 열 벡터는 D에 의한 스파스 분해이다. 스파스 코드의 벡터의 l0-놈의 사용은 근사에서 최소한의 원자들, 예를 들어, z에서의 최소한의 0이 아닌 가중치들만, 사용하는 것이 바람직하다는 것을 의미한다. 일 실시예에서, K-SVD 방법은 각각의 크기, 스케일, 및 진단되는 사람의 기관들의 세부에 대한 양질의 MR 영상들에 있어서 사전을 학습하기 위해 이용될 수 있다.
불균일한 크기의 패치들
개시된 실시예에 따른 방법은 다양한 샘플링 기술들을 이용하여 작동될 수 있다. 대상체의 합성된 영상으로부터 패치들을 추출하는 절차는 이 영상의 체적 부분들의 집합을 생성할 수 있다. 패치의 최적의 종횡비 (x축, y축, 및 z축 크기에 대한 이방성 비율(anisotropy proportion))는 앨리어싱 아티팩트들의 분포 때문에 k-공간 샘플링의 주된 방향에 의존할 수 있다. 만족스러운 근사를 위해서, 패치가 클수록, 사전의 요소들(elements)이 더 많아야 하기 때문에, 사전의 크기가 제한될 때, 패치의 크기는 크지 않아야 한다. 만일 샘플링 레이트가 더 높은 쪽에서, 패치가 더 많은 데이터를 커버하는 경우, 그 패치의 형태는 좋은 선택일 것이다. 패치 당 고정된 양의 정보를 커버하도록 하는 전략은 불균일한 크기의 패치들을 적용함으로서 수행 가능해진다. 이러한 방법은 고정된 크기의 패치들을 사용한 방법보다 더 좋은 품질의 영상 복원을 가능하게 할 수 있다.
최적화 문제 (7)-(8)은 이하에서 설명할 도 5에 도시된 제안된 공동 CS 반복 알고리즘에 의해 해결될 수 있다.
도 5는 예시적인 실시예에 따른 사전 기반 공동 CS 복원 방법에 대해 도시한다.
사전 기반 공동 CS MRI 복원
도 5를 참조하여 보면, 공동 CS 반복 알고리즘의 입력으로서 사용되는 상당히 불완전하게 샘플링된 멀티 코일 k-공간(S501)은 상기에 언급된 공동 분할 브레그만 알고리즘에 의해 초기 전처리 될 수 있다. 솔루션이 수렴(S502)할 때까지, 이하의 단계들이 수행될 수 있다. 먼저, 생성된 코일 영상들의 합성과 멀티 코일 k-공간의 퓨리에 역변환(IFFT) (S504)에 의해 대상체의 합성 영상(S506)이 계산될 수 있다. 이를 위해, 미리 계산된 코일 감도 지도들(S505)이 활용될 수 있다. 중첩된 패치들(S508)의 세트가 합성된 객체(S506)로부터 추출될 수 있(S다(S507). 그리고 나서, (알고리즘의 데이터 및 매개변수들에 의존하는) 각 2차원 또는 3차원의 패치의 벡터화를 위해, 스파스 객체 분해(sparse object decomposition) (S509)가 적합한 사전 D(S510)로 방정식 (9)에 따라 수행될 수 있다. 스파스 객체 분해는 직교 매칭 추적 알고리즘(Orthogonal Matching Pursuit algorithm)에 의해 해결될 수 있다. 스파스 코드들(S511)은 각각의 벡터화된 패치를 위해 획득되지만, 패치들의 근사(S512) 는 사전 요소들(S510)의 선형 조합과 계수들로서 계산된 스파스 코드들(S511)을 통하여 수행될 수 있다. 그 다음, 합성된 객체의 스파스 근사치(S515)는 복원된 패치들을 중첩된 패치들과 조립(S514)함으로써 얻어질 수 있다. 다음 단계(S516)에서, 합성된 객체(S515)는 코일 영상들로 분해되고, 생성된 코일 영상들은 퓨리에 변환(FFT)에 의해 퓨리에 도메인으로 전달될 수 있다. 이러한 연산(S516)은 중간의 멀티 코일 k-공간(intermediate multi-coil k-space)(S517)을 야기할 수 있다. 그 다음, 멀티 코일 k-공간의 샘플링된 위치에서의 스펙트럼 데이터는 자신의 측정된 값으로 복원될 수 있다(S518). 만일 솔루션이 수렴하는 경우(S502), 복원된 멀티 코일 k-공간 (S519)는 공동 CS MRI 복원 알고리즘의 결과로서 출력(S503)될 수 있으며, 수렴하지 않는 경우엔 새로운 반복이 시작될 수 있다.
다중대역 분해
개시된 실시예에 따른 방법에서, 데이터의 다중(multilevel) 구조를 활용하기 위한 방법 중 하나는 서로 다른 주파수 대역에 대응하는 복수의 객체 요소들을 복원하는 것이다. 일 실시예에 따른 방법에서, 멀티 코일 k-공간은, 낮은, 중간 및 높은 주파수 등에 대응하는 여러 주파수 대역으로 분할될 수 있다. 이러한 분할에서, 새로운 멀티 코일 k-공간들이 나타나며, 이 k-공간들은 다른 모든 대역들이 0으로 채워지는 반면, 이들 각각은 측정된 데이터로 채워진 각각의 주파수 대역을 포함할 수 있다. 대역들로 분해되는 동안, 대역 당 신호의 변화가 신호의 총 변화보다 상당히 적기 때문에 신호의 희박성(sparsity)이 개선되며, 이 때문에 더 좋은 품질로 복원될 수 있다. 분해된 k-공간은 제안된 다중 사전 복원의 메인 루프의 입력으로 여겨지며, 병렬로 처리될 수 있다. 모든 주파수 대역에 대해, 모든 멀티 코일 k-공간이 복원되면, 복원된 결과들을 합성함으로써 복원 영상이 얻어질 수 있다. 사전 학습 단계에서 객체는 각 주파수 대역 각각에 대한 사전들을 학습하기 위한 방법과 마찬가지로 각 주파수 대역들로 분해될 수 있다. 복원 단계에서, 각 주파수 대역은 대응하는 사전을 이용하여 복원될 수 있다.
다중 사전 CS MRI 복원
개시된 실시예에 따른 발명의 중요한 포인트는 다중 사전 스파스 객체 근사이다. CS 복원의 사전들 및 관련된 매개변수들(스파스 근사의 허용 오차와 같은)은 복원 절차 동안 변화될 수 있다. 개시된 실시예에 따른 발명에서, CS 복원에 다중 사전들을 적용하는 2가지 접근법이 제안된다.
첫 번째 접근법에서, 사전들은 특정 횟수의 반복 이후에 서로 교체될 수 있다(도 5의 블록 510 참조). 예를 들어, 만일 예비 복원이 이미 수행된 경우라면, 스파스 근사치의 기저를 변경하는 것은 매우 효율적이며, 영상의 세부 디테일들의 복원과 명확화가 요구된다.
또 다른 제안된 접근법은 CS 복원의 한 번의 반복 동안, 동시에 다수의 사전을 이용하는 스파스 객체 근사와 관련된 것이다. 상기 언급한 접근법과 이 접근법의 구별을 위해, 체적 대상체의 영상을 형성하는 패치들의 사전 기반 근사에 대한 면밀한 관찰이 아래에 개시된다.
도 6은 예시적인 실시예에 따른 패치들의 단일 사전 스파스 근사에 대해 도시한다.
구체적으로, 도 6에서, 단일 사전 근사를 목표로 하는 CS MRI 복원의 일부가 도시된다. 대상체의 체적 부분들(volumetric portions of the object)의 벡터화된 패치들(S600)은 사전 기반 스파스 근사의 입력일 수 있다. 이 패치들은 대응하는 스파스 코드들(S603)을 산출하면서, 사전(S602)에 의해 성기게(sparsly) 분해될 수 있다(S601). 그 다음, 요구되는 정확도(스파스 근사 허용 오차 범위)까지의 패치들(S604)의 근사가 수행될 수 있다. 도 6에서, 근사 결과는 (S605)로 도시된다. 패치들의 이러한 근사치들이 이 방식의 출력에 해당하며, 대상체 영상의 조립을 위해 더 사용될 수 있다.
도 7은 예시적인 실시예에 따른 패치들의 다중 사전 스파스 근사에 대해 도시한다.
도 7을 참조하여 보면, 다중 사전 근사 기법에서, 연산의 비교적 간단한 시퀸스(sequence)가 반복되며, 도 7에 도시된 여러 단계 스파스 근사 기법을 야기한다. 체적 대상체로부터 추출된 벡터화된 패치들(S700)이 입력 신호로서 사용될 수 있다. 그 다음, 제1 사전(S702)에 의해 스파스 코드들(S703)이 생성되도록 분해될 수 있다(S701). 여기에 원본 신호가 사용되기 때문에, 제1 사전(S702)을 이하에서 원본 지향적 사전(original-oriented dictionary)이라고 부른다. 그 다음, 패치들의 스파스 근사치(S705)는 직교 매칭 추구 알고리즘(Orthogonal Matching Pursuit (OMP))을 통해 계산될 수 있다(S704). 다음 단계 (S706)에서, 패치들 잔차들(S707)은 입력 패치들(S700)로부터 패치들의 스파스 근사치(S705)를 성분 별로 뺄셈하는 것에 의해 계산될 수 있다. 그 다음, 스파스 근사는 패치 잔차들(S707)에 대하여 반복될 수 있다: 패치 잔차들(S707)은 새로운 스파스 코드들(S710)을 생성하는 제2 사전(S709)에 의해 분해될 수 있다(S708). 제2 사전(S709)은 중간 및 높은 주파수 특징들의 잔차들을 대상으로 하므로, 이하에서 잔차 지향적인 사전(residual-oriented dictionary)이라고 부른다. 마지막으로, 패치들의 근사치들(S705)(제1 사전에 의해 획득된) 및 패치 잔차들의 스파스 근사치들(S712)(제2 사전에 의해 획득된)은 이후 합산될 수 있다(S713). 이것은 원하는 패치들의 다중 사전 스파스 근사치(S714)를 초래할 수 있다. 이러한 패치들의 조립으로, 상기 언급된 단일 사전 접근법에 비해 더 명확한 MRI 복원이 획득될 수 있다: 미세한 디테일이 복원될 수 있으며, 노이즈와 앨리어싱 아티팩트 발생이 비교적 잘 억제될 수 있다.
도 8은 예시적인 실시예에 따른 사전 기반 스파스 근사에 대해 도시한다.
사전 기반 스파스 근사의 예시가 도 8(a)에 도시된다. 또한, 도 8에서, 다중 스케일 사전들(multi-scale dictionaries)(도 8(b)) 및 대상 지향적 사전들(target-oriented dictionaries)(도 8(c))에 의한 서로 다른 종류의 사전 기반 스파스 근사가 도시된다. 다중 스케일 사전들은 4N × 4M 사이즈의 패치 π j (801)는 예를 들어, 2N × 2M, N × M 등의 더 작은 크기의 블록 안에 표현된 정보 요소(informative component)를 가진 패치들(802)에 의해 스파스 근사될 수 있다. 이 패치들(802)의 블록들의 나머지는 0으로 가정되며, 일반적으로, 이 정보 영역들은 메인 4N × 4M 블록 안에 다양한 위치에 포함될 수 있다. 영상의 다양한 영역들은 서로 다른 주파수 대역들에 의한 특징들을 포함하므로, 대상 지향적 사전은, 서로 다른 주파수 대역들에 대응되는 패치들을 포함할 수 있다. 사전들을 낮은, 중간의, 높은 주파수 패치와 같이, 추출된 패치들의 추정된 다양성에 따라, 합성된 영상의 다양한 영역들에 적용하는 것이 더 효율적이라는 것을 알 수 있다. 도 8(c)에서, 서로 다른 컬러들은 다양한 관심영역들을 강조한다: 어두운 컬러의 영역을 위해서, 고주파수 패치들(803)이 사전과 함께 적용될 수 있고, 좀 더 밝은 컬러들에 의해 표시된 두 개의 영역들은 상대적으로 낮은 주파수 패치들(804)이 사전과 함께 적용될 수 있다.
다중 사전 복원: 사전 학습( dictionary learning )
초기 패치들의 스파스 근사 지향적인 사전들의 학습은 K-SVD 알고리즘[M. Aharon, M. Elad, and A. Bruckstein,"K-SVD: 스파스 표현을 위한 오버컴플리트 사전들의 설계를 위한 알고리즘(An Algorithm for Designing Overcomplete Dictionaries for Sparse Representation)," IEEE Trans. Image Process. 54, 4311-4322 (2006)]에 기초한다. 상기 절차의 세부사항은 개시된 실시예에 따른 발명의 범위를 벗어난다. 여기에 여러 사전들로 동작하는 방법이 제안된다. 학습 기술은 영상 패치들(image patches) 및 이른바 패치-잔차(patch-residuals) 즉, 원본 패치들과 원본 패치들을 스파스 근사한 패치들 사이의 차이, 지향적인 다양한 사전들을 위한 것이다. 종래의 사전은 원본의 고 품질의 영상의 세트에서 학습되는 반면, 잔차 지향적인 사전들은 원본 영상과 원본 영상의 원본 지향적인 사전에 의한 스파스 근사들 사이의 차이에서 학습된다. 이하의 용어들, 510(반복적으로 대체된 사전들을 포함한다), 602, 및 702는 원본 지향적인 사전들이다.
도 9는 예시적인 실시예에 따른 원본 및 잔차 지향적인 사전들의 학습 과정을 도시한다.
먼저, 패치 기반 사전을 학습하는 단계(S901)는 대상체의 고 품질 영상의 커다란 데이터 세트(S900)를 사용하여 수행될 수 있다. 이것은 원본 지향적인 사전을 제공할 수 있다. 일 실시예에 따르면, 기존의 K- SVD 알고리즘이 사전 학습을 위해 사용될 수 있다. 다음 단계 903에서, 계산된 원본 지향적인 사전(S902)은 원본 영상들의 패치들과 모든 영상들에 대한 스파스 근사들 사이의 차이를 주어진 데이터 세트(S900)로부터 계산하기 위해 사용될 수 있다. 이것은 잔차들의 데이터 세트(S904)를 발생시킬 수 있다. 그 후, 원하는 잔차 기반 사전(S906)을 계산하기 위해, 패치 기반 사전 학습(S905)이 상기 발생된 패치 잔차들(S904)를 이용하여 수행될 수 있다.
상기 원본 지향적 사전(S902) 및 잔차 지향적 사전(S906)은 원자들(atoms)의 수를 감소시키고, 이로 인해 스파스 근사를 가속화 시키기 위해 그람 슈미트(Gram-Schmidt) 절차에 의해 직교화 될 수도 있다는 것을 주목하라.
도 10은 예시적인 실시예에 따른 공동 CS MRI 복원 시스템을 도시한다.
공동 CSMRI 복원 시스템
제안된 시스템은 CS MRI 의료 영상 복원에 적용가능하며, 도 10에 도시된 시스템을 형성하기 위한 아래의 구성들(units)(예를 들어, 복원 서버)을 포함할 수 있다:
샘플링 방식(언더 샘플링 마스크)을 생성하고, 생성된 샘플링 방식을 데이터 획득부(data-acquisition unit, 1002), 전처리부(preprocessing, 1003), 감도 추정부(sensitivity-estimating unit, 1001), 및 CS MRI 복원부(1005)에게 전송하는 데이터 샘플링부(Data-samling unit, 1000);
데이터 샘플링부 (1000)로부터 수신한 샘플링 방식과 관련하여, 다수의 RF 수신기 코일에 의해 불완전하게 샘플링된 k-공간들을 획득하고 나서, 획득한 멀티 코일 k-공간들을 전처리부(1003) 및 감도 추정부(1001)에게 전송하는 데이터 획득부(Data-acquisition unit, 1002);
샘플링된 멀티 코일 k-공간 데이터에 기초하여, 코일 감도 지도들을 추정하는 감도-추정부(Sensitivity-estimating)(1001). 샘플링된 주파수들의 공간 좌표들은 (데이터 샘플링부(700)에 의해 준비된)샘플링 마스크로부터 획득될 수 있으며, 필요한 멀티 코일 k-공간은 상기 데이터 획득부(1002)로부터 수신될 수 있다;
초기 값들(0 또는 공동 분할 브레그만과 같은 초기치 설정 절차에 의해 제공된 값들)로 멀티 코일 k-공간의 샘플링 되지 않은 위치들을 채우고, 초기치 설정된 멀티 코일 k-공간을 CS MRI 복원부(1005)에게 전송하는 전처리부(1003). 초기치 설정될 k-공간들은 데이터 획득부(1002)로부터 수신되고, 샘플링 되지 않은 위치들에서의 정보는 데이터 샘플링부(1000)로부터 수신되며, 초기치 설정 알고리즘을 위해 필요한 코일 감도 지도들은 감도 추정부(1001)로부터 수신될 수 있다;
스파스 객체 근사를 위한 원본 및/또는 잔차 지향적인 사전들을 학습하고, 상기 사전들을 CS MRI 복원부(1005)에게 전송하는 스파스 근사부(Sparse-approximating unit, 1004);
전처리부(1003)로부터 획득한 전처리된 멀티 코일 k-공간을 다중 사전 공동 CS MRI 알고리즘에 따라 복원하는 CS MRI 복원부(1005). k-공간들에서의 복원될 샘플들의 정보는 데이터 샘플링부(1000)로부터 수신되고, 공동 CS 복원을 위한 사전들과 코일 감도 지도들은 각각 스파스 근사부(1004)와 감도 추정부(1001)로부터 수신될 수 있다.
데이터 샘플링부(1000), 감도 추정부(1001), 및 전처리부(1003)는 멀티 코일 k-공간을 전처리하는 블록을 형성한다. 스파스 근사부(1004) 및 CS 자기 공명 영상(MRI) 복원부(1005)는 다중 사전 공동 CS 복원(multi-dictionary joint compressed sensing reconstruction) 블록을 형성할 수 있다.
도 11은 예시적인 실시예에 따른 다중 사전 CS MRI 복원 시스템을 도시한다.
구체적으로, 도 11는 다중 사전 접근법(이 케이스에서 사전들은 복원이 되는 동안 순차적으로 변경된다) 각각에 대한 복원부의 특정 케이스를 도시한다. 이 시스템은 도 5에 도시된 상술된 방법에 대한 구현을 나타내며, 아래의 구성들(units)을 포함할 수 있다:
다중 RF 수신기 코일들에 의해 대상체의 k-공간들을 획득하고, 몇몇의 적절한 값들(0 또는 공동 분할 브레그만의 초기 추측들)로 샘플링 되지 않은 위치들을 채우는 것을 통하여 획득한 k-공간들의 전처리를 수행하며, 전처리된 k-공간들을 측정치 복원부(1110) 및 작동부(Trigger unit, 1105)에게 전송하는 k-공간 입력부(1100);
코일 감도 지도들을 추정하고, 추정된 코일 감도 지도들을 합성부(1107) 및 분할부(1109)에게 전송하는 감도 입력부(1101);
k-공간 획득의 샘플링 방식(샘플링 마스크)를 생성하고, 생성된 샘플링 방식을 측정치 복원부(1110)에게 전송하는 데이터 샘플링부(1102);
복원될 멀티 코일 k-공간을, 반복된 수 또는 각각의 정지 기준 분석에 의해 생성된 제어부(1103)의 신호에 따라 측정치 복원부(1110)로부터 합성부(1107) 또는 k-공간 출력부(1106)에게 전송하는 작동부(1105);
복원된 멀티 코일 k-공간을 출력하고 저장하는 k-공간 출력부(1106);
작동부(1105)로부터 멀티 코일 k-공간을 수신하고, 수신한 멀티 코일 k공간을 위한 퓨리에 역변환을 수행하며, 감도 입력부(1101)로부터 획득한 코일 감도 지도들을 사용하여 "코일 영상들"을 합성하는 합성부(merging unit, 1107); 그 후 생성된 합성 영상을 다중 사전 근사부(1108)에게 전송된다;
제어부(1103)로부터 현재 반복 수를 수신하고, 사전 기반 객체 근사를 위한 적절한 사전들을 선택하며, 선택된 사전들을 다중 사전 근사부(1108)에게 전송하는 사전 입력부(1104);
합성부(1107)로부터 합성된 영상을 수신하고, 사전 입력부(1104)로부터의 원본 및/또는 잔차 지향적인 사전들에 의해 합성된 영상의 스파스 근사를 수행하며, 결과를 분할부(1109)에게 전송하는 다중 사전 근사부(1108);
감도 입력부(1101)로부터 획득된 코일 감도 지도들을 사용하여 다중 사전 근사부(1108)로부터 획득된 객체 스파스 근사를 분할하는 분할부(Splitting unit, 1109). 또한, 분할부(1109)는 코일들과 관련된 영상들로부터 k-공간들을 획득하고, 획득한 영상들("멀티 코일 k-공간")을 측정치 복원부(1110)에게 전송한다;
분할부(1109)로부터 멀티 코일 K- 공간을 수신하고, k-공간 데이터의 샘플링된 위치들에서의 원본 값을 복원하는 측정치 복원부(1110). 샘플링된 위치들은 데이터 샘플링부(1102)로부터 획득되며, 복원 결과는 작동부(1105)에게 전송된다;
제어부(1103)는 공동 CS MRI 복원의 두 가지 양상을 결정할 수 있다: 이 두 가지 양상은, 각각의 k-공간 전송을 위한 목적지가 되는 구성(unit)을 결정하는 것과 현재 반복에서 사용되는 사전들을 선택하는 것이다. 최초 반복에서, 제어부(1103)는 데이터 처리를 시작하기 위해서, k-공간 입력부(1100)로부터 초기 데이터를 수신하기 위해 작동부(1105)를 스위치할 수 있다. 다른 모든 케이스들에서, 작동부(1105)는 측정치 복원부(1110)로부터 수신한 k-공간을 제공할 수 있다. 더욱이, 이 멀티 코일 k-공간은 합성부(1107)에게 전송될 수 있다. 만일 정지 조건(반복 수 또는 정확성 기준이 정지 조건을 위해 확인됨)이 만족되면, 제어부(1103)는 측정치 복원부(1110)으로부터 합성부(1107)로의 멀티 코일 k-공간 전송이 정지되도록 하기 위해 작동부(1105)를 스위치할 수 있다. 이 케이스에서, 작동부(1105)는 결과 k-공간을 k-공간 출력부(1106)에게 전송할 수 있다. 동시에, 제어부(1103)는 적절한 사전들을 선택하기 위해 현재 반복의 수를 사전 입력부(1104)에게 알릴 수 있다.
도 12는 예시적인 실시예에 따른 원본 지향적 및 잔차 지향적인 사전들을 학습하는 시스템을 도시한다.
도 12은 도 9에 도시된 알고리즘의 구현을 나타낸다. 이 시스템은 아래의 구성들(units)을 포함할 수 있다:
대상체들의 고품질의 MR 영상들의 데이터세트를 수집하고, 수집한 데이터세트를 사전 학습부(1201) 및 스파스 근사부(1202)에게 전송하는 영상 수집부(image-collecting unit, 1200);
(예를 들어, K-SVD에 의해) 사전 학습을 수행하고, 입력 데이터에 따라 원본 및 잔차 지향적인 사전들을 생성하는 사전 학습부(1201). 이러한 모든 사전들은 사전들을 저장하는 사전 출력부(1204)에게 전송될 수 있다. 생성된 원본 지향적인 사전은 사전 학습부(1201)로부터 스파스 근사부(1202)에게 전송되고, 잔차 지향적인 사전은 잔차 근사부(1203)에게 전송될 수 있다;
영상 수집부(1200)으로부터의 MR 영상들의 모든 데이터 세트를 위한 스파스 근사치를 계산하는 스파스 근사부(1202); 상기 계산은 사전 학습부(1201)로부터 획득된 사전과 함께 수행된다. 잔차들은 사전 학습부(1201)와 잔차 근사부(1203)에게 다시 전송될 수 있다;
사전 학습부(1201)로부터 획득된 잔차 지향적인 사전에 의해 스파스 근사부(1202)로부터 수신된 잔차들의 근사를 계산하는 잔차 근사부(1203);
계산된 원본 및 잔차 지향적인 사전들을 수신하고, 수신한 사전들을 저장(저장된 사전들은 스파스 근사를 위해 사용되는 것으로 가정된다)하는 사전 출력부(1204).
도 13은 예시적인 실시예에 따른 자기 공명 영상 처리 시스템을 도시한다.
도 13에 도시된 자기 공명 영상 처리 시스템은 다수의 수신기 코일에 기초하여 불완전하게 샘플링된 멀티 코일 k-공간을 획득하기 위한 MR 스캐너(1301), 복원 서버(1302), 및 오퍼레이터 콘솔(1303)을 포함할 수 있다. 불완전하게 샘플링된 멀티 코일 k-공간이 획득되면, 데이터와 함께 파일은 복원 서버 메모리(1305)에 저장될 수 있다. 오퍼레이터는 오퍼레이터 콘솔(1303)의 컨트롤러(1307)를 이용하여 영상 복원 프로그램을 실행할 수 있다. 이 프로그램은 복원 서버(1302)의 프로세서(1304)에서 작동할 수 있다; 프로세서(1304)는 그것의 서브루틴으로서 불완전하게 샘플링된 MRI 복원을 위한 모듈(1308)을 포함할 수 있다. 상기 프로그램 작동의 결과는 오퍼레이터 콘솔(1303)의 디스플레이(1306)에 표시될 수 있다. 컨트롤러(1307)는 작동의 결과를 제어하고, 복원되고, 디스플레이 될 MR 영상의 매개변수들을 모니터 할 수 있다.
도 14는 예시적인 실시예에 따른 공동 분할 브레그만 복원을 위한 발명의 시스템을 도시한다.
구체적으로 도 14는 공동 분할 브레그만 알고리즘에 의한 데이터 초기치 설정을 위한 시스템(1400)일 수 있다.. 도 14는 도 4에 도시된 알고리즘의 구현을 나타낸다. 이 시스템(1400)은 아래의 구성들(units)을 포함할 수 있다.
시스템(1400)은, 선택된 샘플링 마스크에 따라 수신기 코일에 의해 멀티 코일 k-공간을 획득하고, k-공간 및 샘플링 마스크를 제어 및 갱신부(1401)와 측정치 복원부(1407)에게 전송하는 k-공간 입력부(1409)를 포함할 수 있다.
또한, 시스템(1400)은, 수신기 코일의 추정된 감도 지도들을 수신하고, 수신한 감도 지도들을 합성부(1403) 및 분할부(1406)에게 전송하는 감도 입력부(1402)를 포함할 수 있다.
또한, 시스템(1400)은, 제어 및 갱신부(1401)로부터 멀티 코일 k-공간을 수신하고, k-공간들의 퓨리에 역변환을 수행하고, 감도 지도들(감도 입력부(1402)로부터 획득된)을 이용하여 획득된 영상들을 합성하며, 합성된 결과 영상을 객체 추정부(1404)에게 전송하는 합성부(1403)를 포함할 수 있다.
또한, 시스템(1400)은, 합성부(1403)로부터 필요한 매개변수들과 합성된 영상을 수신하고, 각각의 기준(criterion) 최소화에 의해 객체 스파스 근사를 다시 계산하며, 계산된 결과를 내부 매개변수 갱신부(1405)에게 전송하는 객체 추정부(1404)를 포함할 수 있다. 객체 추정부(1404)는 l1노름(norm)에 포함된 평가 함수를 최소화함으로써 상기 객체의 스파스 근사치를 재계산할 수 있다.
또한, 시스템(1400)은, 객체 추정부(1404)로부터 갱신된 객체 (영상)의 스파스 표현을 수신하고, 관련 매개변수들을 갱신하며, 갱신된 객체 근사를 분할부(1406)에게 전송하는 내부 매개변수 갱신부(1405)를 포함할 수 있다.
또한, 시스템(1400)은, 내부 매개변수 갱신부(1405)로부터 획득된 객체 스파스 근사를 이 목적을 위한 코일 감도 지도들(감도 지도들은 감도 입력부(1402)로부터 획득된다)을 이용하여 분할하는 분할부(1406)를 포함할 수 있다. 상기, 분할부(1406)는 생성된 코일 영상들의 퓨리에 변환을 수행하고, 갱신된 멀티 코일 k-공간을 제어 및 갱신부(1401)에게 전송할 수 있다.
제어 및 갱신부 (1401)는 공동 분할 브레그만 알고리즘의 외부 및 내부 루프를 위한 프로세싱 상태를 제어할 수 있다. 내부 루프에 대한 조건이 충족되는 경우, 제어 및 갱신부는 외부 매개변수들을 갱신할 수 있다. 만일 외부 루프에 대한 조건이 충족되는 경우, 복원된 k-공간들은 측정치 복원부(1407)에게 전송되고, 반면에 예비적으로 재구성된 k-공간들은 합성부(1403)에게 전송될 수 있다.
또한, 시스템(1400)은, 제어 및 갱신부(1401)로부터 복원된 멀티 코일 k-공간을 수신하고, k-공간 데이터의 샘플링된 성분들을 그 성분들의 (측정된) 원본 값으로 복원하는 측정치 복원부(1407)를 포함할 수 있다.
또한, 시스템(1400)은, 측정치 복원부(1407)로부터 수신된 복원된 멀티 코일 k-공간을 출력하고 저장하는 k-공간 출력부(1408)를 포함할 수 있다.
도 15는 예시적인 실시예에 따른 사용자 시나리오 1을 도시한다.
상기 시나리오는 세 가지 역할을 전제로 할 수 있다.
MR 스캐너(1520)와 MR 영상 복원 시스템(1530)을 제어하는 의사 또는 오퍼레이터(1510), 측정 프로세스를 실현하고, 대상체의 스펙트럼 데이터 측정치를 제공하는 MR 스캐너(1520), 및 불완전하게 샘플링된 스펙트럼 데이터에 기초하여 MR 영상 복원 방법을 실현하는 소프트웨어 및 하드웨어를 포함하는 MR 영상 복원 시스템(1530)의 역할을 전제할 수 있다.
사용자 시나리오1의 주요한 단계는:
의사(1510)가 환자의 선택된 기관(organ)에 대하여 조사(investigation)를 위한 방식으로 환자를 MR 스캐너(1520) 안에 위치시키는 단계(S1511), 의사(1510)가 원하는 스캐닝 모드 및 조사할 기관(organ)을 선택하는 단계(S1512), 의사(1510)가 스캐닝 절차의 속도를 높이기 위해, 원하는 샘플링 마스크를 선택하는 단계(S1513), 의사(1510)가 영상 복원 절차를 위해 매개변수 세트를 선택하는 단계(S1514), 의사(1510)가 스캐닝 절차를 시작하는 단계(S1515), MR 스캐너(1520)가 영상 복원 절차의 속도를 높이기 위해서, 몇몇의 데이터를 건너뛰면서 환자를 스캔하는 단계(S1516, 언더 샘플링), MR 영상 복원 시스템(1530)이 획득한 데이터를 처리하고, 최종적으로, 조사할 기관의 MR 영상을 복원하는 단계(S1517) (상기 MR 영상 복원 시스템에 의해 구현되는 알고리즘의 주요 단계는 도 3과 관련된 본문에 설명되어 있다); 및 의사(1510)가 복원된 MR 영상을 조사(investigation)하는 단계(S1518)를 포함할 수 있다.
개시된 실시예에 따른 발명의 주요한 효과는 최소한의 시간으로 좋은 품질의 MR 영상을 제공하는 것이다(스캔된 데이터는 불완전하게 샘플링되며, 스캔 시간의 총량은 "품질/ 소요시간"에 대해 소정의 절충을 함에 따라4-6배 감소할 수 있다).
도 16는 예시적인 실시예에 따른 사용자 시나리오 2를 도시한다.
사용자 시나리오 2에서, 복원 절차는 MR 스캐너(1620)와 함께 제공될 수 있는 추가적인 하드웨어에 의해 수행되는 것이 아니라, 복원 서비스를 제공하는 전용 서버에서 수행될 수 있다. 복원 서버의 주인(owner) 및 복원 알고리즘은 메디컬 센터와의 계약 하에서 접근을 허용할 수 있다. 이 경우, 모든 (의료 센터에 위치한) MR 스캐너(1620)는 복원 알고리즘에 대한 어떠한 참조(reference) 없이도 사용될 수 있다. 복원은 클라우드 컴퓨팅(cloud computing) 기술을 이용하여 수행될 수 있으며, 특히, 분산 컴퓨팅(distributed computing)을 사용하여 수행될 수 있다.
사용자 시나리오2의 주요 단계는:의사(또는 오퍼레이터)(1610)가 환자를 의료 센터의 MR 스캐너 안에 위치 시키는 단계(S1611), 의사(또는 오퍼레이터)(1610)가 바람직한 가속 인자(acceleration factor)를 선택하는 단계(S1612), 의사(또는 오퍼레이터)(1610)가 스캔 절차를 시작하는 단계(S1613),; MR 스캐너(1620)가 환자를 스캔하고, 불완전하게 샘플링된 스펙트럼 데이터를 제공하는 단계(S1614), MR 스캐너(1620)가 불완전하게 샘플링된 데이터를 복원 서비스를 제공하는 클라우드에 있는 복원 서버(1630)에게 전송하는 단계(S1615), 클라우드에 있는 복원 서버(1630)가 수신한 불완전하게 샘플링된 데이터로부터 3D 영상 복원을 수행하는 단계(S1616), 클라우드에 있는 복원 서버(1630)가 복원된 영상을 다시 의료 센터에 있는 MR 스캐너(1620)에게 전송하는 단계(S1617), 의료 센터에 있는 의사(1610)가 수신한 복원된 영상을 조사(investigate)하는 단계(S1618)를 포함 할 수 있다.
도 17은 예시적인 실시예에 따른 사용자 시나리오 3을 도시한다.
도 17에 도시된 사용자 시나리오는 비파괴 분석을 위한 다른 사용자 시나리오일 수 있다. 대상체는 간섭 레이저 빔(coherent laser beam)으로 조사(illuminate)되고, 산란 방사선은 스캐너에 의해 검출될 수 있으며, 대상체의 영상은 컴퓨터의 위상 검색 알고리즘(phase-retrieval algorithm)에 의한 많은 강도(intensity) 측정치들로부터 정확하게 복원될 수 있다. 사용된 강도 측정치들은 서로 다른 거리로의 웨이브필드(wave field) 전파를 묘사하는 전달 함수 때문에 서로 다를 수 있다. 웨이브필드 전파의 정확한 거리들은 종종 특정 회절 요소들: 공간 광 변조기(SLMs)에 의해 인코딩 될 수 있다. 공간 광 변조기들의 충전율(fill factor)은 완벽해지는 것과 거리가 멀기 때문에(공간 광 변조기의 비용에 의존하여, 시판되는 SLM의 충전율은 72~98%가 될 수 있다.), 작동하는 퓨리에 스펙트럼은 일반적으로 불완전하게 샘플링될 수 있다. 복원된 영상들에서 링잉 노이즈(ringing artifacts) 생성을 피하기 위해, 사용된 전달함수들의 보상이 필요할 수 있다. 예를 들어, 이것은 [A. Migukin외, "4f 광학 시스템에 있어서의 위상 검색: 백그라운드 보상 및 이진 진폭 객체의 스파스 정규화(Phase retrieval in 4f optical system: background compensation and sparse regularization of object with binary amplitude)," Appl. Opt. 52, A269-A280 (2013)]에서 비롯된 두 단계의 검색 단계를 통해 실현될 수 있다.
사용자 시나리오3은, 관심 객체에게 레이저 빔을 조사하고, 객체로부터 산란 방사선을 감지하며, (SLM에 의해 인코딩 된) 서로 다른 전달 함수에 대응하는 다수의 강도 측정치들을 획득하는 레이저 스캐너(1720), 영상 복원을 위해 필요한 소프트웨어가 설치된 영상 복원 시스템(컴퓨터)(1730),및 영상 복원 시스템(1730)을 제어하고, 복원된 영상을 조사하는 사용자(오퍼레이터)(1710)를 포함하는 3가지 역할을 가정할 수 있다.
사용자 시나리오3의 주요 단계는:
사용자(오퍼레이터)(1710)가 관심 대상체(또는 그것의 단편) 앞에 레이저 스캐너(1720)를 위치시키는 단계(S1711), 사용자(오퍼레이터)(1710)가 레이저 스캐너(1720)를 이용하여 관심 대상체(또는 그것의 단편)을 스캔하는 단계(S1712), 사용자(오퍼레이터)(1710)가 테스트 대상체(test-object) 앞에 레이저 스캐너(1720)를 위치 시키고, 레이저 스캐너(1720)를 이용하여 테스트 대상체를 스캔하는 단계(S1713), 영상 복원 시스템(컴퓨터)(1730)이 알려진 테스트 대상체 및 관심 대상체의 스캔 데이터를 수신하는 단계(S1714), 영상 복원 시스템(컴퓨터)(1730)이 테스트 대상체의 스캔 데이터에 의해 전달 함수를 수정하는 단계(S1715), 관심 대상체의 영상을 수정된 전달 함수들을 이용하여, 관심 대상체의 스캔 데이터로부터 복원하는 단계(S1716), 및 사용자(오퍼레이터)가 대상체의 복원된 영상을 조사하는 단계(S1717)를 포함할 수 있다.
개시된 실시예에 따른 발명의 기술적 효과는 최소한의 시간으로 양질의 MR 영상을 제공하는 것이며(상기 취득된 데이터가 불완전하게 샘플링 되기 때문에 총 획득시간이 4-6배 정도 비례하여 감소한다), 즉, 공지의 솔루션과 비교하면 본 발명은 비슷하거나 더 좋은 품질의 MR 영상을 얻는 반면, 데이터 획득시간은 상당히 가속화 하는 것을 가능하게 할 수 있다. 본 명세서는 본 발명의 권리범위를 명확히 하고, 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자가 본 발명을 실시할 수 있도록, 본 발명의 원리를 설명하고, 실시예들을 개시한다. 개시된 실시예들은 다양한 형태로 구현될 수 있다.
도 18은 MRI 시스템의 개략도이다.
도 18을 참조하면, MRI 시스템(1)은 오퍼레이팅부(10), 제어부(30) 및 스캐너(50)를 포함할 수 있다. 여기서, 제어부(30)는 도 18에 도시된 바와 같이 독립적으로 구현될 수 있다. 또는, 제어부(30)는 복수 개의 구성 요소로 분리되어, MRI 시스템(1)의 각 구성 요소에 포함될 수도 있다. 이하에서는 각 구성 요소에 대해 구체적으로 살펴보도록 한다.
스캐너(50)는 내부 공간이 비어 있어, 대상체가 삽입될 수 있는 형상(예컨대, 보어(bore) 형상)으로 구현될 수 있다. 스캐너(50)의 내부 공간에는 정자장 및 경사자장이 형성되며, RF 신호가 조사된다.
스캐너(50)는 정자장 형성부(51), 경사자장 형성부(52), RF 코일부(53), 테이블부(55) 및 디스플레이부(56)를 포함할 수 있다. 정자장 형성부(51)는 대상체에 포함된 원자핵들의 자기 쌍극자 모멘트의 방향을 정자장 방향으로 정렬하기 위한 정자장을 형성한다. 정자장 형성부(51)는 영구 자석으로 구현되거나 또는 냉각 코일을 이용한 초전도 자석으로 구현될 수도 있다.
경사자장 형성부(52)는 제어부(30)와 연결된다. 제어부(30)로부터 전송 받은 제어신호에 따라 정자장에 경사를 인가하여, 경사자장을 형성한다. 경사자장 형성부(52)는 서로 직교하는 X축, Y축 및 Z축 방향의 경사자장을 형성하는 X, Y, Z 코일을 포함하며, 대상체의 부위 별로 공명 주파수를 서로 다르게 유도할 수 있도록 촬영 위치에 맞게 경사 신호를 발생 시킨다.
RF 코일부(53)는 제어부(30)와 연결되어, 제어부(30)로부터 전송 받은 제어신호에 따라 대상체에 RF 신호를 조사하고, 대상체로부터 방출되는 MR 신호를 수신할 수 있다. RF 코일부(53)는 세차 운동을 하는 원자핵을 향하여 세차운동의 주파수와 동일한 주파수의 RF 신호를 대상체에게 전송한 후 RF 신호의 전송을 중단하고, 대상체로부터 방출되는 MR 신호를 수신할 수 있다.
RF 코일부(53)는 원자핵의 종류에 대응하는 무선 주파수를 갖는 전자파를 생성하는 송신 RF 코일과, 원자핵으로부터 방사된 전자파를 수신하는 수신 RF 코일로서 각각 구현되거나 또는 송/수신 기능을 함께 갖는 하나의 RF 송수신 코일로서 구현될 수도 있다. 또한, RF 코일부(53)외에, 별도의 코일이 대상체에 장착될 수도 있다. 예를 들어, 촬영 부위 또는 장착 부위에 따라, 헤드 코일(Head coil), 척추 코일(spine coil), 몸통 코일(torso coil), 무릎 코일(knee coil) 등이 별도의 코일로 이용될 수 있다.
스캐너(50)의 외측 및/또는 내측에는 디스플레이부(56)가 마련될 수 있다. 디스플레이부(56)는 제어부(30)에 의해 제어되어, 사용자 또는 대상체에게 의료 영상 촬영과 관련된 정보를 제공할 수 있다.
또한, 스캐너(50)에는 대상체의 상태에 관한 모니터링정보를 획득하여 전달하는 대상체 모니터링정보 획득부가 마련될 수 있다. 예를 들어, 대상체 모니터링정보 획득부(미도시)는 대상체의 움직임, 위치 등을 촬영하는 카메라(미도시), 대상체의 호흡을 측정하기 위한 호흡 측정기(미도시), 대상체의 심전도를 측정하기 위한 ECG 측정기(미도시), 또는 대상체의 체온을 측정하기 위한 체온 측정기(미도시)로부터 대상체에 관한 모니터링정보를 획득하여 제어부(30)로 전달할 수 있다. 이에 따라, 제어부(30)는 대상체에 관한 모니터링정보를 이용하여 스캐너(50)의 동작을 제어할 수 있다. 이하에서는 제어부(30)에 대해 살펴보도록 한다.
제어부(30)는 스캐너(50)의 전반적인 동작을 제어할 수 있다.
제어부(30)는 스캐너(50) 내부에서 형성되는 신호들의 시퀀스를 제어할 수 있다. 제어부(30)는 오퍼레이팅부(10)로부터 수신받은 펄스 시퀀스(pulse sequence) 또는 설계한 펄스 시퀀스에 따라 경사자장 형성부(52) 및 RF 코일부(53)를 제어할 수 있다.
펄스 시퀀스란, 경사자장 형성부(52), 및 RF 코일부(53)를 제어하기 위해 필요한 모든 정보를 포함하며, 예를 들어 경사자장 형성부(52)에 인가하는 펄스(pulse) 신호의 강도, 인가 지속시간, 인가 타이밍 등에 관한 정보 등을 포함할 수 있다.
제어부(30)는 펄스 시퀀스에 따라 경사 파형, 즉 전류 펄스를 발생시키는 파형 발생기(미도시), 및 발생된 전류 펄스를 증폭시켜 경사자장 형성부(52)로 전달하는 경사 증폭기(미도시)를 제어하여, 경사자장 형성부(52)의 경사자장 형성을 제어할 수 있다.
제어부(30)는 RF 코일부(53)의 동작을 제어할 수 있다. 예를 들어, 제어부(30)는 공명 주파수의 RF 펄스를 RF 코일부(53)에 공급하여 RF 신호를 조사할 수 있고, RF 코일부(53)가 수신한 MR 신호를 수신할 수 있다. 이때, 제어부(30)는 제어신호를 통해 송수신 방향을 조절할 수 있는 스위치(예컨대, T/R 스위치)의 동작을 제어하여, 동작 모드에 따라 RF 신호의 조사 및 MR 신호의 수신을 조절할 수 있다.
제어부(30)는 대상체가 위치하는 테이블부(55)의 이동을 제어할 수 있다. 촬영이 수행되기 전에, 제어부(30)는 대상체의 촬영 부위에 맞추어, 테이블부(55)를 미리 이동시킬 수 있다.
제어부(30)는 디스플레이부(56)를 제어할 수 있다. 예를 들어, 제어부(30)는 제어신호를 통해 디스플레이부(56)의 온/오프 또는 디스플레이부(56)를 통해 표시되는 화면 등을 제어할 수 있다.
제어부(30)는 MRI 시스템(1) 내 구성요소들의 동작을 제어하기 위한 알고리즘, 프로그램 형태의 데이터를 저장하는 메모리(미도시), 및 메모리에 저장된 데이터를 이용하여 전술한 동작을 수행하는 프로세서(미도시)로 구현될 수 있다. 이때, 메모리와 프로세서는 각각 별개의 칩으로 구현될 수 있다. 또는, 메모리와 프로세서는 단일 칩으로 구현될 수도 있다.
오퍼레이팅부(10)는 MRI 시스템(1)의 전반적인 동작을 제어할 수 있다. 오퍼레이팅부(10)는 영상 처리부(11), 입력부(12) 및 출력부(13)를 포함할 수 있다.
영상 처리부(11)는 메모리를 이용하여 제어부(30)로부터 수신 받은 MR 신호를 저장하고, 이미지 프로세서를 이용하여 영상 복원 기법을 적용함으로써, 저장한 MR 신호로부터 대상체에 대한 영상 데이터를 생성할 수 있다.
예를 들어, 영상 처리부(11)는 메모리의 k-공간(예컨대, 푸리에(Fourier) 공간 또는 주파수 공간이라고도 지칭됨)에 디지털 데이터를 채워 k-공간 데이터가 완성되면, 이미지 프로세서를 통해 다양한 영상 복원기법을 적용하여(예컨대, k-공간 데이터를 역 푸리에 변환하여) k-공간 데이터를 영상 데이터로 복원할 수 있다.
또한, 영상 처리부(11)가 MR 신호에 대해 적용하는 각종 신호 처리는 병렬적으로 수행될 수 있다. 예를 들어, 다채널 RF 코일에 의해 수신되는 복수의 MR 신호를 병렬적으로 신호 처리하여 영상 데이터로 복원할 수도 있다. 한편, 영상 처리부(11)는 복원한 영상 데이터를 메모리에 저장하거나 또는 후술할 바와 같이 제어부(30)가 통신부(60)를 통해 외부의 서버에 저장할 수 있다.
입력부(12)는 사용자로부터 MRI 시스템(1)의 전반적인 동작에 관한 제어 명령을 입력 받을 수 있다. 예를 들어, 입력부(12)는 사용자로부터 대상체 정보, 파라미터 정보, 스캔 조건, 펄스 시퀀스에 관한 정보 등을 입력 받을 수 있다. 입력부(12)는 키보드, 마우스, 트랙볼, 음성 인식부, 제스처 인식부, 터치 스크린 등으로 구현될 수 있다.
출력부(13)는 영상 처리부(11)에 의해 생성된 영상 데이터를 출력할 수 있다. 또한, 출력부(13)는 사용자가 MRI 시스템(1)에 관한 제어 명령을 입력 받을 수 있도록 구성된 유저 인터페이스(User Interface, UI)를 출력할 수 있다. 출력부(13)는 스피커, 프린터, 디스플레이 등으로 구현될 수 다.
한편, 도 18에서는 오퍼레이팅부(10), 제어부(30)를 서로 분리된 객체로 도시하였으나, 전술한 바와 같이, 하나의 기기에 함께 포함될 수도 있다. 또한, 오퍼레이팅부(10), 및 제어부(30) 각각에 의해 수행되는 프로세스들이 다른 객체에서 수행될 수도 있다. 예를 들어, 영상 처리부(11)는, 제어부(30)에서 수신한 MR 신호를 디지털 신호로 변환하거나 또는, 제어부(30)가 직접 변환할 수도 있다.
MRI 시스템(1)은 통신부(60)를 포함하며, 통신부(60)를 통해 외부 장치(미도시)(예를 들면, 서버, 의료 장치, 휴대 장치(스마트폰, 태블릿 PC, 웨어러블 기기 등))와 연결할 수 있다.
통신부(60)는 외부 장치와 통신을 가능하게 하는 하나 이상의 구성 요소를 포함할 수 있으며, 예를 들어 근거리 통신 모듈(미도시), 유선 통신 모듈(61) 및 무선 통신 모듈(62) 중 적어도 하나를 포함할 수 있다.
통신부(60)가 외부 장치로부터 제어 신호 및 데이터를 수신하고, 수신된 제어 신호를 제어부(30)에 전달하여 제어부(30)로 하여금 수신된 제어 신호에 따라 MRI 시스템(1)을 제어하도록 하는 것도 가능하다.
또는, 제어부(30)가 통신부(60)를 통해 외부 장치에 제어 신호를 송신함으로써, 외부 장치를 제어부의 제어 신호에 따라 제어하는 것도 가능하다.
예를 들어 외부 장치는 통신부(60)를 통해 수신된 제어부(30)의 제어 신호에 따라 외부 장치의 데이터를 처리할 수 있다.
외부 장치에는 MRI 시스템(1)을 제어할 수 있는 프로그램이 설치될 수 있는바, 이 프로그램은 제어부(30)의 동작의 일부 또는 전부를 수행하는 명령어를 포함할 수 있다.
프로그램은 외부 장치에 미리 설치될 수도 있고, 외부장치의 사용자가 어플리케이션을 제공하는 서버로부터 프로그램을 다운로드하여 설치하는 것도 가능하다. 어플리케이션을 제공하는 서버에는 해당 프로그램이 저장된 기록매체가 포함될 수 있다.
한편, 개시된 실시예들은 컴퓨터에 의해 실행 가능한 명령어 및 데이터를 저장하는 컴퓨터로 읽을 수 있는 기록매체의 형태로 구현될 수 있다. 상기 명령어는 프로그램 코드의 형태로 저장될 수 있으며, 프로세서에 의해 실행되었을 때, 소정의 프로그램 모듈을 생성하여 소정의 동작을 수행할 수 있다. 또한, 상기 명령어는 프로세서에 의해 실행되었을 때, 개시된 실시예들의 소정의 동작들을 수행할 수 있다.
이상에서와 같이 첨부된 도면을 참조하여 개시된 실시예들을 설명하였다. 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자는 본 발명의 기술적 사상이나 필수적인 특징을 변경하지 않고도, 개시된 실시예들과 다른 형태로 본 발명이 실시될 수 있음을 이해할 것이다. 개시된 실시예들은 예시적인 것이며, 한정적으로 해석되어서는 안 된다.

Claims (20)

  1. 프로세서 및 상기 프로세서에 연결된 메모리를 포함하고,
    상기 프로세서는,
    대상체에 대한 불완전하게 샘플링된 멀티 코일 k-공간을 획득하고,
    미리 획득된 제1 사전 및 상기 제1 사전을 이용하여 획득된 제2 사전을 이용하여, 상기 멀티 코일 k-공간을 복원하여 상기 대상체의 복원 영상을 획득하는, 자기 공명 영상 처리 장치.
  2. 제1항에 있어서, 상기 제2 사전은
    미리 획득된 대상체의 영상과 상기 제1 사전에 기초하여 상기 영상을 근사한 것의 잔차를 기초하여 획득되는, 자기 공명 영상 처리 장치.
  3. 제1항에 있어서, 상기 근사는 스파스 근사인, 자기 공명 영상 처리 장치.
  4. 제1항에 있어서, 상기 프로세서는,
    상기 불완전하게 샘플링된 멀티 코일 k-공간에 기초하여 제1 영상을 획득하고,
    상기 제1 영상으로부터 획득된 균일하지 않은 크기의 패치들, 상기 제1 사전 및 상기 제2 사전에 기초하여 상기 멀티 코일 k-공간을 복원하여 상기 복원 영상을 획득하는, 자기 공명 영상 처리 장치.
  5. 제1항에 있어서, 상기 프로세서는,
    상기 불완전하게 샘플링된 멀티 코일 k-공간에 기초하여 획득된 서로 다른 주파수 대역에 대응되는 패치들, 상기 제1 사전 및 상기 제2 사전에 기초하여 상기 멀티 코일 k-공간을 복원하여 상기 복원 영상을 획득하는, 자기 공명 영상 처리 장치.
  6. 제1항에 있어서,
    오퍼레이터 콘솔을 더 포함하고,
    상기 오퍼레이터 콘솔은 복원 영상을 시각화 하기 위한 표시부를 포함하는, 자기 공명 영상 처리 장치.
  7. 제6항에 있어서,
    상기 오퍼레이터 콘솔은 제어부를 더 포함하고,
    상기 제어부는 상기 자기 공명 영상 처리 장치를 제어하도록 구성되는, 자기 공명 영상 처리 장치.
  8. 제1항에 있어서,
    상기 불완전하게 샘플링된 멀티 코일 k-공간 데이터를 획득하도록 구성된 MR 스캐너를 더 포함하는, 자기 공명 영상 처리 장치.
  9. 제1항에 있어서,
    상기 프로세서는,
    상기 획득한 불완전하게 샘플링된 멀티 코일 k-공간을 복수 개의 주파수 대역으로 분할하고,
    상기 복수 개의 주파수 대역으로 분할된 불완전하게 샘플링된 멀티 코일 k-공간들을 각각 상기 제1 사전 및 상기 제2 사전을 이용하여 복원함으로써 복원 영상을 획득하는, 자기 공명 영상 처리 장치.
  10. 제1항에 있어서, 상기 프로세서는,
    상기 불완전하게 샘플링된 멀티 코일 k-공간을 전처리 하고, 상기 전처리된 멀티 코일 k-공간을, 상기 제1 사전 및 상기 제2 사전을 이용하여 복원하는, 자기 공명 영상 처리 장치.
  11. 제10항에 있어서, 상기 프로세서는,
    상기 불완전하게 샘플링된 멀티 코일 k-공간의 샘플링 되지 않은 위치들을 초기 추측(initial guess)으로 채우고, 상기 불완전하게 샘플링된 멀티 코일 k-공간을 코일 압축(coil compressing)함으로써, 상기 불완전하게 샘플링된 멀티 코일 k-공간을 전처리 하는, 자기 공명 영상 처리 장치.
  12. 자기 공명 영상을 복원하는 방법에 있어서,
    대상체에 대한 불완전하게 샘플링된 멀티 코일 k-공간을 획득하는 단계; 및
    미리 획득된 제1 사전 및 상기 제1 사전을 이용하여 획득된 제2 사전을 이용하여, 상기 멀티 코일 k-공간을 복원하여 상기 대상체의 복원 영상을 획득하는 단계를 포함하는, 자기 공명 영상 복원 방법.
  13. 제12항에 있어서, 상기 제2 사전은
    미리 획득된 대상체의 영상과 상기 제1 사전에 기초하여 상기 영상을 근사한 것의 잔차를 기초하여 획득되는, 자기 공명 영상 복원 방법.
  14. 제12항에 있어서, 상기 복원 영상을 획득하는 단계는
    상기 불완전하게 샘플링된 멀티 코일 k-공간에 기초하여 제1 영상을 획득하는 단계; 및
    상기 제1 영상으로부터 획득된 균일하지 않은 크기의 패치들, 상기 제1 사전 및 상기 제2 사전에 기초하여 상기 멀티 코일 k-공간을 복원하여 상기 복원 영상을 획득하는 단계를 포함하는, 자기 공명 영상 복원 방법.
  15. 제12항에 있어서, 상기 복원 영상을 획득하는 단계는
    상기 불완전하게 샘플링된 멀티 코일 k-공간에 기초하여 획득된 서로 다른 주파수 대역에 대응되는 패치들, 상기 제1 사전 및 상기 제2 사전에 기초하여 상기 멀티 코일 k-공간을 복원하여 상기 복원 영상을 획득하는 단계를 포함하는, 자기 공명 영상 복원 방법.
  16. 제12항에 있어서,
    상기 복원 영상을 디스플레이하는 단계를 더 포함하는, 자기 공명 영상 복원 방법.
  17. 제12항에 있어서, 상기 대상체에 대한 불완전하게 샘플링된 멀티 코일 k-공간을 획득하는 단계는
    MR 스캐너로부터 상기 불완전하게 샘플링된 멀티 코일 k-공간 데이터를 수신하는 단계를 포함하는, 자기 공명 영상 복원 방법.
  18. 제12항에 있어서, 상기 복원 영상을 획득하는 단계는
    상기 획득한 불완전하게 샘플링된 멀티 코일 k-공간을 복수 개의 주파수 대역으로 분할하는 단계; 및
    상기 복수 개의 주파수 대역으로 분할된 불완전하게 샘플링된 멀티 코일 k-공간들을 각각 상기 제1 사전 및 상기 제2 사전을 이용하여 복원함으로써 복원 영상을 획득하는 단계를 포함하는, 자기 공명 영상 복원 방법.
  19. 제12항에 있어서, 상기 복원 영상을 획득하는 단계는
    상기 불완전하게 샘플링된 멀티 코일 k-공간을 전처리 하는 단계; 및
    상기 전처리된 멀티 코일 k-공간을, 상기 제1 사전 및 상기 제2 사전을 이용하여 복원하는 단계를 포함하는, 자기 공명 영상 복원 방법.
  20. 제12항 내지 제19항에 기재된 자기 공명 영상을 복원하는 방법의 각 단계를 수행하기 위한 프로그램이 저장된, 컴퓨터로 판독 가능한 기록 매체.
KR1020160087383A 2015-09-04 2016-07-11 자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치 KR102527812B1 (ko)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP16842129.5A EP3346443B1 (en) 2015-09-04 2016-08-02 Method for restoring magnetic resonance image and magnetic resonance image processing apparatus
PCT/KR2016/008494 WO2017039163A1 (ko) 2015-09-04 2016-08-02 자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치
US15/756,666 US20180247436A1 (en) 2015-09-04 2016-08-02 Method for restoring magnetic resonance image and magnetic resonance image processing apparatus

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
RU2015137794A RU2626184C2 (ru) 2015-09-04 2015-09-04 Способ, устройство и система для реконструкции магнитно-резонансного изображения
RU2015137794 2015-09-04

Publications (2)

Publication Number Publication Date
KR20170028829A true KR20170028829A (ko) 2017-03-14
KR102527812B1 KR102527812B1 (ko) 2023-05-02

Family

ID=58454163

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020160087383A KR102527812B1 (ko) 2015-09-04 2016-07-11 자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치

Country Status (5)

Country Link
US (1) US20180247436A1 (ko)
EP (1) EP3346443B1 (ko)
KR (1) KR102527812B1 (ko)
RU (1) RU2626184C2 (ko)
WO (1) WO2017039163A1 (ko)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101883806B1 (ko) * 2017-04-07 2018-07-31 울산과학기술원 영상 복원 장치 및 방법
KR20210075826A (ko) * 2019-12-13 2021-06-23 한양대학교 산학협력단 학습 가능한 공간-스펙트럴 변환 기법 기반의 딥러닝 영상 분해 및 압축 센싱 방법

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101806902B1 (ko) * 2016-02-16 2017-12-12 삼성전자주식회사 자기 공명 영상 복원 방법 및 장치
WO2019090447A1 (en) * 2017-11-10 2019-05-16 Shenzhen United Imaging Healthcare Co., Ltd. System and method for image reconstruction
US10573031B2 (en) * 2017-12-06 2020-02-25 Siemens Healthcare Gmbh Magnetic resonance image reconstruction with deep reinforcement learning
US10489943B2 (en) 2018-02-28 2019-11-26 General Electric Company System and method for sparse image reconstruction
DE102018214636A1 (de) * 2018-08-29 2020-03-05 Siemens Healthcare Gmbh Verfahren zum Erzeugen eines Magnetresonanzbildes, Computerprogramm, Datenspeicher und Magnetresonanzanlage
CN109471053B (zh) * 2018-10-18 2020-01-31 电子科技大学 一种基于双约束的介电特性迭代成像方法
US11860256B2 (en) * 2021-04-19 2024-01-02 Arizona Board Of Regents On Behalf Of The University Of Arizona MRI streak-artifact suppression, system and method
WO2023114317A1 (en) * 2021-12-14 2023-06-22 Regents Of The University Of Minnesota Noise-suppressed nonlinear reconstruction of magnetic resonance images

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040165757A1 (en) * 2001-04-13 2004-08-26 Jingfei Ma Adaptive data differentiation and selection from multi-coil receiver to reduce artifacts in reconstruction
US20120099774A1 (en) * 2010-10-21 2012-04-26 Mehmet Akcakaya Method For Image Reconstruction Using Low-Dimensional-Structure Self-Learning and Thresholding
KR101330638B1 (ko) * 2012-05-10 2013-11-18 고려대학교 산학협력단 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8310233B2 (en) * 2009-02-18 2012-11-13 Mayo Foundation For Medical Education And Research Method for image reconstruction from undersampled medical imaging data
GB2495447B (en) * 2010-06-24 2017-04-12 Koninklijke Philips Nv Dynamic contrast enhanced MR imaging with compressed sensing reconstruction
RU2013104364A (ru) * 2010-07-02 2014-08-10 Конинклейке Филипс Электроникс Н.В. Компьютерный программный продукт, осуществляемый компьютером способ и система магнитно-резонансной визуализации для получения магнитно-резонансного изображения
WO2014006633A1 (en) * 2012-07-02 2014-01-09 Arjun Arunachalam A method and system for aliased k-space data restoration

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040165757A1 (en) * 2001-04-13 2004-08-26 Jingfei Ma Adaptive data differentiation and selection from multi-coil receiver to reduce artifacts in reconstruction
US20120099774A1 (en) * 2010-10-21 2012-04-26 Mehmet Akcakaya Method For Image Reconstruction Using Low-Dimensional-Structure Self-Learning and Thresholding
KR101330638B1 (ko) * 2012-05-10 2013-11-18 고려대학교 산학협력단 자기 공명 영상 생성 방법 및 그에 따른 자기 공명 영상 생성 장치

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
SAIPRASAD RAVISHANKAR, 'MR Image Reconstruction From Highly Undersampled k-Space Data by Dictionary Learning’, IEEE TRANSACTIONS ON MEDICAL IMAGING, VOL. 30, NO. 5, MAY 2011(2011.05.31.) 1부.* *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101883806B1 (ko) * 2017-04-07 2018-07-31 울산과학기술원 영상 복원 장치 및 방법
KR20210075826A (ko) * 2019-12-13 2021-06-23 한양대학교 산학협력단 학습 가능한 공간-스펙트럴 변환 기법 기반의 딥러닝 영상 분해 및 압축 센싱 방법

Also Published As

Publication number Publication date
EP3346443A4 (en) 2018-10-17
EP3346443B1 (en) 2021-10-27
US20180247436A1 (en) 2018-08-30
RU2626184C2 (ru) 2017-07-24
RU2015137794A (ru) 2017-03-10
WO2017039163A1 (ko) 2017-03-09
KR102527812B1 (ko) 2023-05-02
EP3346443A1 (en) 2018-07-11

Similar Documents

Publication Publication Date Title
KR102527812B1 (ko) 자기 공명 영상 복원을 위한 방법 및 자기 공명 영상 처리 장치
KR101659578B1 (ko) 자기 공명 영상 처리 방법 및 장치
Fabian et al. Data augmentation for deep learning based accelerated MRI reconstruction with limited data
CN111656392B (zh) 用于合成磁共振图像的系统和方法
US8692549B2 (en) Method for reconstructing images of an imaged subject from a parallel MRI acquisition
RU2557334C2 (ru) Динамическая контрастная улучшенная мр визуализация с реконструкцией сжатого измерения
US9482732B2 (en) MRI reconstruction with motion-dependent regularization
KR20150125543A (ko) 자기 공명 영상 장치 및 자기 공명 영상의 생성 방법
JP5085832B2 (ja) サブサンプリングされた捕捉を用いた磁気共鳴撮像方法
US10705168B2 (en) System and method for low rank approximation of high resolution MRF through dictionary fitting
JP6332891B2 (ja) サイドバンドアーチファクトを抑制する並列マルチスライスｍｒ撮像
KR101605130B1 (ko) 자기 공명 영상 장치 및 그에 따른 자기 공명 영상의 이미징 방법
US10527695B2 (en) Systems and methods for efficient magnetic resonance fingerprinting scheduling
JP6548531B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
KR101525014B1 (ko) 자기 공명 영상 장치 및 그에 따른 자기 공명 영상의 이미징 방법
JP2007503904A (ja) 高磁場磁気共鳴撮像のための適応的画像均一性補正
KR101806902B1 (ko) 자기 공명 영상 복원 방법 및 장치
KR101580532B1 (ko) 자기 공명 영상 처리 장치 및 방법
US20060164087A1 (en) Magnetic resonance method and device
CN113558600A (zh) 磁共振成像装置以及图像处理方法
KR102393288B1 (ko) 자기공명영상장치 및 그 제어방법
KR101826063B1 (ko) 화학 교환 포화 전이 자기 공명 영상 장치 및 방법
JP5156958B2 (ja) 磁気共鳴イメージング装置
CN113567901A (zh) 一种磁共振旋转坐标系下的自旋晶格弛豫成像方法和系统
WO2015170394A1 (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