KR100677017B1 - 0차 위상 검출 방법 및 mri 시스템 - Google Patents
0차 위상 검출 방법 및 mri 시스템 Download PDFInfo
- Publication number
- KR100677017B1 KR100677017B1 KR1020020065514A KR20020065514A KR100677017B1 KR 100677017 B1 KR100677017 B1 KR 100677017B1 KR 1020020065514 A KR1020020065514 A KR 1020020065514A KR 20020065514 A KR20020065514 A KR 20020065514A KR 100677017 B1 KR100677017 B1 KR 100677017B1
- Authority
- KR
- South Korea
- Prior art keywords
- phase detection
- zero
- phase
- pulse sequence
- echo
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56554—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Signal Processing (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
본 발명은 MR 신호의 0차 위상을 정확하게 나타내는 0차 위상을 검출하도록 한다. MR 신호의 퓨리에 변환에 의한 모든 샘플링점의 복소 벡터를 이용하여 계산된 합성 벡터는 0차 위상으로 채택된다.
Description
도 1은 본 발명의 실시예에 따른 MRI 시스템의 블럭도,
도 2는 위상 검출 펄스 시퀀스를 도시하는 도면,
도 3은 본 발명의 실시예에 따른 0-차 위상 검출을 나타내는 흐름도,
도 4는 본 발명에 따른 O-차 위상 및 O-차 위상 차이를 나타내는 개념도,
도 5는 멀티-샷 확산 강화 EPI 법에 의한 이미징 펄스 시퀀스를 도시하는 도면,
도 6은 이미징 데이터의 수집 궤적을 나타내는 개념도,
도 7은 위상 검출 시퀀스를 도시하는 도면,
도 8은 종래의 1차 위상 보정과 0차 위상 보정에 관한 설명도,
도 9는 종래의 O-차 위상 및 O-차 위상 차이를 나타내는 개념도.
도면의 주요 부분에 대한 부호의 설명
1 : 자석 어셈블리 3 : 자기장 구배 구동 회로
4 : RF 전력 증폭기 5 : 전치 증폭기
8 : 순차 메모리 회로 9 : 게이트 변조 회로
11 : A/D 변환기 13 : 조작자 콘솔
본 발명은 0-차 위상 검출 방법 및 MRI(Magnetic Resonance Imaging) 시스템에 관한 것으로, 좀 더 구체적으로 MR 신호의 0-차 위상을 정확하게 검출할 수 있는 0차 위상 검출 방법 및 MRI 시스템에 관한 것이다.
도 5는 멀티-샷 확산 강화(multi-shot diffusion enhancement) EPI(Echo Plannar Imaging) 법에 의한 이미징 펄스 시퀀스의 기본예를 도시한다.
이 이미징 펄스 시퀀스에 따르면, 여기 펄스(excitation pulse) RF90과 슬라이스-선택 자기장 구배(slice-selective magnetic field gradient) SG90을 인가한다. 다음에 MPG(Motion Probing Gradient) 펄스 MPG을 인가한다. 반전 RF 펄스 RF180과 반전 슬라이스-선택적 자기장 구배 SG180을 인가한다. 다음에 MPG 펄스 MPG를 인가한다. 다음에 위상-인코딩 자기장 구배 pdn을 인가한다. 교대로 극성을 반전하는 데이터 수집 판독 자기장 구배(data acquisition readout magnetic field gradient)의 펄스 rl, ···, rm을 연속적으로 인가한다. 또한 판독 자기장 구배가 극성을 반전할 때에 위상-인코딩 자기장 구배의 p2, ···, pM을 인가한다. 제 1 에코 e1로부터 제 M 에코 eM이 순서대로 집속하는 때에 샘플링하여, 에코 el, ···, eM에 기반을 두는 이미징 데이터 항목 F(n, 1), ···, F(n, M)를 수집한다. 펄스의 이러한 인가는 위상-인코딩 자기장 구배 pdn을 바꾸면서 n번(여기서 n은 1에서부터 N까지임) 반복된다. k 공간에 기록되는 이미징 데이터 항목 F(1,1) 내지 F(N, M)를 수집한다. 이러한 이미징 기법을 N 샷 M 에코라고 한다. 또한, 시간적으로 순차적인 샷에 대해 순서대로 붙인 번호 n을 샷 번호라고 한다. 또한, 어떤 샷에 응답하여 반송된 후에 시간적으로 순차적으로 집속하는 에코에 대해 시간 순서대로 붙인 번호를 에코 번호라고 한다.
도 6은 k 공간 KS에서의 이미징 데이터 항목 F(1,1) 내지 F(N, M)이 수집되는 궤적(trajectory)을 나타내는 예시도이다. 여기서, N은 4, M은 4를 나타낸다.
k 공간 KS를 위상 인코드축 방향에 제 1 행에서부터 제 N*M 행(도 6으로서는 제 16행)까지 분할한다. 이 경우에 있어서, 제(n+(m-1) N) 행에 기록되는 이미징 데이터 F(n,m)이 제 n 샷에 응답하여 반송되는 제 m 에코에 기반을 두고 수집될 수 있도록 위상-인코딩 자기장 구배 pdn의 펄스 p2,..., pM을 인가한다.
EPI 법은 위상 오류에 민감하다. 그러므로 MR 신호의 위상을 검출하여, 위상 보정을 해야 한다.
도 7은 위상 검출용 펄스 시퀀스의 일례이다.
이 위상 검출 펄스 시퀀스는 위상-인코딩 자기장 구배가 배제된다는 점에서 도 5에 도시된 이미징 펄스 시퀀스와는 다르다.
순서대로 집속하는 제 1 위상 검출 에코 E1로부터 제 M 위상 검출 에코 EM을 기초로 하여, 제1 위상 검출 데이터 D_1로부터 제 M 위상 검출 데이터 D_M을 수집 한다.
다음으로, 복소 벡터 Z (n)를 얻기 위해 제 m 위상 검출 데이터 D_m을 퓨리에 변환한다. 여기서, 1≤m≤M, 1≤n≤N 이다.
다음으로, 제 m 위상 검출 데이터 D_m로 나타내어지는 제 1 차 위상 φ1_m은 다음식에 의해 산출된다.
여기서, arg{}는 복소수의 편각을 나타내는 함수이다. 다음식에 의해 일차 위상이 보정된다.
여기서, exp{}는 지수함수를 나타낸다.
다음으로, 제 m 위상 검출 데이터 D_m으로 나타내는 0-차 위상 φ0는 다음식에 의해 산출된다.
상기 1차 위상 φ1_m 및 O차 위상 φ0_m은 이미징 데이터의 보정이나 그 이외의 여러가지의 처리에 이용된다.
또한, 제 m 위상 검출용 데이터 D_m으로 나타내는 0차 위상 φ0_m과 제 (m+1) 위상 검출용 데이터 D_m+1으로 나타내는 0차 위상 φ0_m+1사이의 0차 위상 차이 △φ0_m은 다음식으로 계산된다.
O 차 위상 차이 △φ0_m도 여러 가지의 처리에 이용된다.
도 8은 복소 벡터 Z(n), Zcorl(n), Zcor(n)의 위상을 도시하는 개념도이다.
1차 위상 및 0차 위상을 보정함으로써, 제 1 위상 검출용 데이터 항목 D_1 내지 제 M 위상 검출 데이터 항목 D_M에 대한 위상 오류의 악영향을 없앨 수 있다.
1차 위상 및 O차 위상의 보정을 연속적으로 하기 위해 식 (3)에서 제공된 0차 위상 φ0_m이 식 (1)에서 제공된 1차 위상 φ1_m과 결합하여 이용되는한, 문제는 발생하지 않는다.
그러나, 식 (3)에 의해 제공되는 O차 위상 φ0_m이 보정 등에 단독으로 이용되는 경우에는, 보정이 올바르게 되지 않는 문제점이 야기된다. 이 이유는, 식 (3)에 의해 제공되는 O차 위상 φ0_m이, 제 1 샘플링점의 복소 벡터 Z(1)의 위상 그 자체이며, 제 1 샘플링점의 복소 벡터 Z(1) 내지 제 N 샘플링점의 복소 벡터 Z(N)의 O차 위상을 정확하게 표현하지 않기 때문이다.
또한, EPI 법이나 GRASE (GRadient And Spin Echo) 법에 있어서, 제 m 위상 검출 데이터 D_m을 제공하도록 적용되는 판독 자기장 구배 펄스의 극성은 제 (m+1) 위상 검출 데이터 D_m+1을 제공하도록 적용되는 판독 자기장 구배 극성과 반대이다. 이 경우에 있어서, 도 9에 도시하는 바와 같이 식 (5)에 따라 산출되는 0차 위상 차이 △φ0_m은 정확한 0차 위상 차이를 나타내지 않는다. 도 9를 참조하면, "양극의 판독 자기장 구배 펄스"는 도 7에 도시된 기수의 판독 자기장 구배 펄스 rl, r3··등을 가리킨다. "음극의 판독 자기장 구배 펄스"는 도 7에 도시된 우수의 판독 자기장 구배 펄스 r2, r4 등을 가리킨다.
그러므로, 본 발명의 제 1 목적은, 제 1 샘플링점 내지 제 N 샘플링점에서 복소벡터 Z(1) 내지 복소벡터 Z(N)의 임의의 O차 위상을 정확하게 나타내는 O차 위상을 검출할 수 있는 0차 위상 검출 방법 및 MRI 시스템을 제공하는 것이다.
본 발명의 제 2 목적은 연속 위상 검출 에코(consecutive phase detection echoes)에 대한 판독 자기장 구배(readout magnetic field gradient)의 극성이 반전하는 경우, 정확한 0차 위상 차이를 검출할 수 있는 0차 위상 검출 방법 및 MRI 장치를 제공하는 것이다.
본 발명의 제 1 측면에 따르면, 0차 위상 φ0을 검출하는 0차 위상 검출 방법이 제공된다. 여기서 위상 검출 데이터는, 이미징 펄스 시퀀스와는 달리 위상-인코딩 자기장 구배 펄스를 포함하지 않는 위상 검출 펄스 시퀀스에 따라 집속(refocuse)되는 위상 검출 에코에 기초하여 구해진다. 위상 검출 데이터는 퓨리에 변환되어 다음과 같이 표현되는 제 n 샘플링점에서의 복소 벡터 Z(n)를 계산한다.
합성 벡터(composite vector) Zsum은 다음 식에 의해 계산된다.
0차 위상 φ0은 다음과 같이 표현되는 합성 벡터 Zsum을 이용하여 검출된다.
제 1 측면에 따른 0차 위상 검출 방법에 있어서, 제 1 샘플링점 내지 제 n 샘플링점에서 복소 벡터 Z(1) 내지 복소 벡터 Z(N)을 사용하여 계산되는 합성 벡터 Zsum의 위상은 MR 신호의 0차 위상으로 채택된다. 종래 기술과 같이 제 1 샘플링점의 복소 벡터 Z(1)의 위상이 MR 신호의 0차 위상으로 채택되는 경우에 비교하여, 제 1 샘플링점 내지 제 N 샘플링점에서 복소 벡터 Z(1) 내지 복소 벡터 Z(N)의 임의의 O차 위상(any of the 0-order phases)을 정확하게 나타내는 0차 위상을 검출할 수 있다.
본 발명의 제 2 측면에 따르면, 상기 0차 위상 검출 방법에 있어서, 제 1 시간 간격 동안에 집속하는 위상 검출 에코에 기초하여 계산된 합성 벡터를 Zsum_1로 가정하고 제 2 시간 간격 동안에 집속하는 위상 검출 에코에 기초하여 계산된 합성 벡터를 Zsum_2로 가정하면, 0차 위상 차이 △φ0가 다음과 같이 표현된다.
제 2 측면에 따른 0차 위상 검출 방법에 있어서, 제 1 위상 검출 에코에 기초하여 구해진 복소 벡터를 이용하여 계산된 합성 벡터 Zsum의 위상 arg{Zsum_1}와 제 2 위상 검출 에코에 기초하여 구해진 복소 벡터를 이용하여 계산된 합성 벡터 Zsum의 위상 arg{Zsum 2}의 차이는 MR 신호의 0차 위상 차이 △φ0로 채택된다. 제 1 위상 검출 에코에 대해 적용된 판독 자기장 구배 펄스의 극성이 제 2 위상 검출 에코에 대해 적용된 판독 자기장 구배 펄스의 위상과 반전하는 경우라도, 정확한 0차 위상 차이가 계산될 수 있다.
본 발명의 제 3 측면에 따르면, 전술한 0차 위상 검출 방법에 있어서, 제 1 시간 간격 동안에 집속하는 위상 검출 에코에 기초하여 계산된 합성 벡터를 Zsum_1 로 가정하고 제 2 시간 간격 동안에 집속하는 위상 검출 에코에 기초하여 계산된 합성 벡터를 Zsum_2로 가정하면, 0차 위상 차이 △φ0가 다음과 같이 표현된다.
제 3 측면에 따른 0차 위상 검출 방법은 제 2 측면에 따른 0차 위상 검출 방법과 균등(equivalent)하여 정확한 0차 위상 차이가 계산될 수 있다.
본 발명의 제 4 측면에 따르면, 전술한 0차 위상 검출 방법에 있어서, 이미징 펄스 시퀀스는 판독 자기장 구배의 극성을 반전시킴으로써 에코를 집속하려 한다.
제 4 측면에 따른 0차 위상 검출 방법에 있어서, 이미징 펄스 시퀀스는 판독 자기장 구배의 극성을 반전시킴으로써 에코를 집속하려 한다. 펄스 시퀀스는 위상 에러에 민감하기 때문에, 0차 위상 및 0차 위상 차이를 정확하게 검출하는 것이 필수적이다.
본 발명의 제 5 측면에 따르면, 전술한 0차 위상 검출 방법에 있어서, 이미징 펄스 시퀀스는 EPI법 또는 GRASE법에 적합하다.
제 5 측면에 따른 0차 위상 검출 방법에 있어서, 펄스 시퀀스는 EPI법 또는 GRASE법에 적합하다. 펄스 시퀀스는 위상 에러에 민감하기 때문에, 0차 위상 및 0차 위상 차이를 정확하게 검출하는 것이 필수적이다.
본 발명의 제 6 측면에 따르면, 전술한 0차 위상 검출 방법에 있어서, 제 1 시간 간격 동안에 집속하는 위상 검출 에코 및 제 2 시간 간격 동안에 집속하는 위상 검출 에코는 연속 에코(consecutive echoes)이다.
제 6 측면에 따른 0차 위상 검출 방법에 있어서, 연속 에코가 다루어지기 때문에, 판독 자기장 구배의 극성은 반전한다. 그럼에도 불구하고, 0차 위상 및 0차 위상 차이는 정확하게 검출될 수 있다.
본 발명의 제 7 측면에 따르면, 전술한 0차 위상 검출 방법에 있어서, 위상 검출 데이터는, 이미징 데이터가 이미징 펄스 시퀀스에 의해 수집되는 스캔과는 다른 기준 스캔(reference scan) 동안 위상 검출 펄스 시퀀스에 의해 수집된다.
제 7 측면에 따른 0차 위상 검출 방법에 있어서, 위상 검출 데이터는 이미징 펄스 시퀀스가 이용되는 스캔과는 다른 기준 스캔 동안에 수집된다. 데이터 수집은 시간적 제약이 적어진다.
본 발명의 제 8 측면에 따르면, 전술한 0차 위상 검출 방법에 있어서, 펄스는, 펄스가 이미징 펄스 시퀀스와 일치하기 전에 위상 검출 펄스 시퀀스에 따라 인가된다.
제 8 측면에 따른 0차 위상 검출 방법에 있어서, 펄스는, 펄스가 이미징 펄스 시퀀스와 일치하기 전에 위상 검출 펄스 시퀀스에 따라 인가된다. 그러므로, 위상 검출 펄스 시퀀스를 이용하여 검출된 0차 위상 및 0차 위상 차이는 이미징 펄스 시퀀스를 이용하여 수행되는 데이터 수집 동안에 이용될 수 있다.
본 발명의 제 9 측면에 따르면, 주로 라디오-주파수(RF) 전송 수단, 구배 펄스 인가 수단, MR 신호 수신 수단, 위상 검출 데이터 수집 수단, 퓨리에 변환 수단 및 0차 위상 계산 수단으로 구성되는 MRI 시스템이 제공된다. 위상 검출 데이터 수집 수단은 RF 펄스 전송 수단, 구배 펄스 인가 수단 및 MR 신호 수신 수단을 제 어한다. 위상 검출 데이터 수집 수단은, 이미징 펄스 시퀀스와 유사한 그리고 위상-인코딩 자기장 구배를 포함하는 펄스 시퀀스를 이용하여 수신되는 위상 검출 에코에 기초하여 위상 검출 데이터를 수집한다. 퓨리에 변환 수단은 위상 검출 데이터를 퓨리에 변환하여 복소 벡터를 계산한다. 제 n 샘플링점의 복소 벡터 Z(n)을 다음의 식과 같다고 가정하면:
합성 벡터 Zsum은 다음과 같이 계산된다:
0차 위상 φ0는 다음과 같이 표현되는 합성 벡터 Zsum을 이용하여 검출된다.
제 9 측면에 따른 MRI 시스템에 있어서, 0차 위상 검출 방법은 제 1 측면에 따라 구현되는 것이 바람직하다.
본 발명의 제 10 측면에 따르면, 전술한 MRI 시스템은 0차 위상 차이 계산 수단을 더 포함한다. 제 1 시간 간격 동안에 집속하는 위상 검출 에코를 이용하여 계산된 합성 벡터를 Zsum_1로 가정하고 제 2 시간 간격 동안에 집속하는 위상 검출 에코를 이용하여 계산된 합성 벡터를 Zsum_2로 가정하여, 0차 위상 차이 계산 수단은 다음과 같이 0차 위상 차이 △φ0를 계산한다:
제 10 측면에 따른 MRI 시스템에 있어서, 0차 위상 검출 방법은 제 2 측면에 따라 구현되는 것이 바람직하다.
본 발명의 제 11 측면에 따르면, 전술한 MRI 시스템은 0차 위상 차이 계산 수단을 더 포함한다. 제 1 시간 간격 동안에 집속하는 위상 검출 에코를 이용하여 계산된 합성 벡터를 Zsum_1로 가정하고 제 2 시간 간격 동안에 집속하는 위상 검출 에코를 이용하여 계산된 합성 벡터를 Zsum_2로 가정하면, 0차 위상 차이 계산 수단은 다음과 같이 0차 위상 차이 △φ0를 계산한다:
제 11 측면에 따른 MRI 시스템에 있어서, 0차 위상 검출 방법은 제 3 측면에 따라 구현되는 것이 바람직하다.
본 발명의 제 12 측면에 따르면, 전술한 MRI 시스템에 있어서, 이미징 펄스 시퀀스는 판독 자기장 구배의 극성을 반전시킴으로써 에코를 집속하려 한다.
제 12 측면에 따른 MRI 시스템에 있어서, 0차 위상 검출 방법은 제 4 측면에 따라 구현되는 것이 바람직하다.
본 발명의 제 13 측면에 따르면, 전술한 MRI 시스템에 있어서, 이미징 펄스 시퀀스는 EPI 법 또는 GRASE 법에 적합하다.
제 13 측면에 따른 MRI 시스템에 있어서, 0차 위상 검출 방법은 제 5 측면에 따라 구현되는 것이 바람직하다.
본 발명의 제 14 측면에 따르면, 전술한 MRI 시스템에 있어서, 제 1 시간 간격 동안에 집속하는 위상 검출 에코 및 제 2 시간 간격 동안에 집속하는 위상 검출 에코는 연속 에코이다.
제 14 측면에 따른 MRI 시스템에 있어서, 0차 위상 검출 방법은 제 6 측면 에 따라 구현되는 것이 바람직하다.
본 발명의 제 15 측면에 따르면, 전술한 0차 위상 검출 방법에 있어서, 위상 검출 데이터는, 이미징 데이터가 이미징 펄스 시퀀스에 의해 수집되는 스캔과는 다른 기준 스캔 동안 위상 검출 펄스 시퀀스에 의해 수집된다.
제 15 측면에 따른 MRI 시스템에 있어서, 0차 위상 검출 방법은 제 7 측면에 따라 구현되는 것이 바람직하다.
본 발명의 제 16 측면에 따르면, 전술한 MRI 시스템에 있어서, 펄스는, 펄스가 이미징 펄스 시퀀스와 일치하기 전에 위상 검출 펄스 시퀀스에 따라 인가된다.
제 16 측면에 따른 MRI 시스템에 있어서, 0차 위상 검출 방법은 제 8 측면에 따라 구현되는 것이 바람직하다.
본 발명이 구현하는 0차 위상 검출 방법 및 MRI 시스템에 따르면, MR 신호의 퓨리에 변환에 의한 모든 샘플링점의 복소 벡터를 이용하여 계산된 합성 벡터의 위상은 0차 위상으로 채택된다. 제 1 샘플링점의 복소 벡터의 위상을 0차 위상으로 하는 종래 기술에 비교하여, 검출된 0차 위상은 MR 신호의 0차 위상을 정확하게 나타낸다.
본 발명의 또 다른 목적 및 장점은 첨부한 도면에 예시된 바와 같이 후속하는 본 발명의 바람직한 실시예의 설명으로부터 분명해질 것이다.
본 발명을 도면에 도시된 실시예와 관련하여 자세히 설명한다.
도 1은 본 발명의 실시예에 따른 MRI 시스템을 도시하는 블럭도이다.
이 MRI 장치(100)에 있어서, 자석 어셈블리(1)는, 내부에 피검체를 삽입하기 위한 구멍(bore)을 갖는다. 이 구멍 주위에 영구 자석(1p)과, 구배 코일(1g) 및 수신기 코일(1r)이 배치된다. 영구 자석(1p)은 일정한 주 자장을 피검체에 인가한다. 구배 코일(1g)은 x축, y축, z축 코일(슬라이스-선택 축(slice-selective axis), 판독 축 및 위상-인코딩 축은 코일이 결합되는 방법에 따라 결정됨)을 포함하고, 자기장 구배를 발생시킨다. 송신 코일(1t)은 피검체내의 원자핵의 스핀을 여기하는 RF 펄스를 송신한다. 수신 코일(1r)은 피검체에서 야기된 MR 신호를 수신한다. 구배 코일(1g), 송신 코일(1t) 및 수신 코일(1r)은 각각 자기장 구배 구동 회로(3), RF 전력 증폭기(4) 및 전치증폭기(front-end amplifier)(5)에 접속되어 있다. 또, 초전도 자석이 영구 자석(1p)을 대신할 수도 있다.
컴퓨터(7)는, 펄스 시퀀스를 구성하여, 순차 기억 회로(sequent memory circuit)(8)에 건네 준다.
순차 기억 회로(8)는 펄스 시퀀스를 기억한다. 그 펄스 시퀀스에 근거하여 자기장 구배 구동 회로(3)를 작동한다. 자석 어셈블리(1)에 포함된 구배 코일(1g)은 자기장 구배를 발생시킨다. 또한, 게이트 변조 회로(9)는 RF 발진기 회로(10)의 캐리어 출력 신호를 사전 결정된 시간, 사전 결정된 포락선을 갖는 펄스 형성 신호로 변조하여, 펄스 형성 신호를 RF 펄스로서 RF 전력 증폭기(4)에 인가한다. RF 전력 증폭기(4)가 RF 펄스를 증폭하고 그 결과 펄스를 자석 어셈블리(1)에 포함된 송신 코일(1t)에 인가한다.
전치 증폭기(5)는 자석 어셈블리(1)에 포함된 수신 코일(1r)에서 수신한 MR 신호를 증폭하여, 그 결과 신호를 위상 검출기(12)에 전송한다. 위상 검출기(12)는 RF 발진 회로(10)의 캐리어 출력 신호를 참조 신호로서 사용하여, MR 신호의 위상을 검출하고, 그 MR 신호를 AD 변환기(11)에 전송한다. AD 변환기(11)는 아날로그 신호를 디지털 신호로 변환하여 그 디지털 신호를 컴퓨터(7)에 전송한다.
컴퓨터(7)는 AD 변환기(11)에서 생성된 데이터를 판독하고, 위상 검출, 위상 보정, 화상 재구성 등을 실행하여 이미지를 생성한다. 이 이미지는 디스플레이 장치(6)에서 디스플레이된다.
또한, 컴퓨터(7)는 조작자 콘솔(13)로부터 입력된 정보를 수신하고 전체적인 제어를 담당한다.
도 2는 본 발명의 실시예에 이용되는 위상 검출 펄스 시퀀스를 도시한다.
이 위상 검출 펄스 시퀀스는 도 5에 도시된 멀티-샷 확산 강화 EPI 법에 적합한 이미징 펄스 시퀀스와는 달리 위상-인코딩 자기장 구배를 포함하지 않는다.
특히, 여기 펄스 RF90와 슬라이스-선택 자기장 구배 SG90를 인가한 다음에 MPG(Motion Probing Gradient) 펄스 MPG을 인가한다. 다음에 반전 RF 펄스 RF180와 반전 슬라이스-선택 자기장 구배 SG180을 인가한 다음에 다른 MPG 펄스 MPG를 인가한다. 다음에 교대로 극성을 반전하는 데이터 수집 판독 자기장 구배의 펄스 rl, ···, rM을 연속적으로 인가하나 위상-인코딩 자기장 구배의 펄스는 인가하지 않는다. 순서대로 집속하는 제 1 위상 검출 에코 E1 내지 제 M 위상 검출 에코 EM을 기초로 하여 제 1 위상 검출 데이터 D_1 내지 제 M 위상 검출 데이터 D_M을 수집한다.
이 위상 검출 펄스 시퀀스는 참조 스캔으로 채택된다. 이 참조 스캔을 완료한 뒤, 이미징 펄스 시퀀스에 따라 이미징 데이터를 수집하기 위해 스캔이 실행된다.
도 3은 본 발명의 실시예에서 수행되는 0차 위상 검출을 나타내는 흐름도이다.
단계 S1에서, 에코 번호 카운터 m을 1로 초기화한다.
단계 S2에서, 제 m 위상 검출 데이터 D_m은 복소 벡터 Z(n)_m을 구하기 위해 판독축 방향에 1차원 퓨리에 변환한다. 여기서, n은 샘플링점 번호를 나타내며, 1≤n≤N 이다.
단계 S3에서, 제 n 샘플링점의 복소벡터 Z (n)_m을 다음과 같이 표현된다:
합성 벡터 Zsum_m은 다음과 같이 계산된다:
단계 S4에서, 다음과 같이 표현되는 합성벡터 Zsum_m을 이용하여 0차 위상 φ0_m이 검출된다.
단계 S5에서, 에코 번호 카운터 m을 1씩 증가시킨다. 단계 S6에서, 카운터 값 m이 2이면 단계 S2에 되돌아간다. 카운터 값 m이 3 이상이면 단계 S7로 진행한다.
단계 S7에서, 0차 위상 차이 △φ_m-2는 다음과 같이 계산된다:
단계 S8에서, 카운터 값 m이 3 내지 M 이면 단계 S2에 되돌아간다. 카운터 값 m이 M+1이면 처리를 종료한다.
또, 단계 S7에서, 0차 위상 차이 △φ_m-2는 다음과 같이 계산된다:
도 4에 도시하는 바와 같이, 검출된 0차 위상 φ0은 모든 샘플링점의 복소벡터의 임의의 0차 위상을 나타낸다. 따라서, 연속 위상 검출 에코에 대해 인가되는 판독 자기장 구배의 극성이 반전하더라도, 검출된 O차 위상 차이는 정확하다. 따라서, 이미징 데이터의 보정 및 0차 위상 차이 △φ0가 이용되는 다른 여러 가지 처리는 정확한 데이터를 구할 수 있다.
본 발명의 사상과 범주를 벗어나지 않으면서 본 발명의 상당히 다른 다수의 실시예를 형성할 수도 있다. 본 발명은 첨부한 청구항 범위에서 정의된 것을 제외하고는, 명세서에 설명된 특정 실시예로 제한받지 않는다.
본 발명은 0차 위상 검출 방법 및 MRI 시스템에 관한 것으로, MR 신호를 퓨리에 변환하여 얻은 모든 샘플링점의 복소벡터의 합성 벡터의 위상을 0차 위상으로 하는 종래 기술과 같이 제 1 샘플링점의 복소벡터의 위상을 O 차 위상으로 하는 경우에 비교하여, MR 신호의 0차 위상을 정확하게 나타내는 0차 위상을 구할 수 있 다.
Claims (16)
- 0차 위상 φ0을 검출하는 0차 위상 검출 방법에 있어서,위상-인코딩 자기장 구배(phase-encoding magnetic field gradient)를 포함하지 않는 위상 검출 펄스 시퀀스에 따라 집속(refocusing)되는 위상 검출 에코에 기초하여 수집한 위상 검출 데이터가, 제 n 샘플링점의 복소벡터 Z(n)를 계산하기 위해 퓨리에 변환되고,상기 제 n 샘플링점의 복소 벡터 Z(n)이 다음과 같이 표현된다고 가정(assume)하고,합성 벡터 Zsum은 아래의 식에 따라 계산되며,아래와 같이 상기 합성 벡터 Zsum을 이용하여 0차 위상 φ0를 검출하는O차 위상 검출 방법.
- 제 1 항에 있어서,상기 이미징 펄스 시퀀스가 판독 자기장 구배의 극성을 반전시켜 에코를 집속시키려 하는O 차 위상 검출 방법.
- 제 1 항에 있어서,상기 이미징 펄스 시퀀스는 EPI(echo planar imaging) 법 또는 GRASE(gradient and spin echo) 이미징 법에 적합한0차 위상 검출 방법.
- 제 1 항에 있어서,상기 제 1 시간 간격 동안에 집속하는 위상 검출 에코와 상기 제 2 시간 간격 동안에 집속하는 위상 검출 에코가 연속 에코(consecutive echoes)인O차 위상 검출 방법.
- 제 1 항에 있어서,이미징 펄스 시퀀스에 의해 이미징 데이터를 수집하는 스캔과는 다른 참조 스캔(reference scan)동안에 위상 검출 펄스 시퀀스에 의해 위상 검출 데이터를 수집하는0차 위상 검출 방법.
- 제 1 항에 있어서,펄스는 이미징 펄스 시퀀스에 일치하기 전에 위상 검출 펄스 시퀀스에 따라 인가되는0차 위상 검출 방법.
- RF 펄스 송신 장치와,구배 펄스 인가 장치와,MR 신호 수신 장치와,상기 RF 펄스 송신 장치, 구배 펄스 인가 장치 및 MR 신호 수신 장치를 제어하여, 위상-인코딩 자기장 구배를 포함하지 않는 위상 검출 펄스 시퀀스에 의해 수신된 위상 검출 에코에 기초하여 위상 검출 데이터를 수집하는 위상 검출 데이터 수집 장치와,상기 위상 검출 데이터를 퓨리에 변환하여 복소 벡터를 계산하는 퓨리에 변환 장치와,제 n 샘플링점의 상기 복소 벡터 Z(n)를 다음과 같이 가정하고,이하의 식에 따라 합성 벡터 Zsum을 계산하여아래와 같이 상기 합성 벡터 Zsum에 의해 0차 위상 φ0을 검출하는 0차 위상 계산 장치를 포함하는 MRI 시스템.
- 제 9 항에 있어서,상기 이미징 펄스 시퀀스가 판독 자기장 구배의 극성을 반전시켜 에코를 집속시키려 하는MRI 시스템.
- 제 9 항에 있어서,상기 이미징 펄스 시퀀스가 EPI 법 또는 GRASE 이미징 법에 적합한 MRI 시스템.
- 제 9 항에 있어서,상기 제 1 시간 간격 동안에 집속하는 위상 검출 에코와 상기 제 2 시간 간격 동안에 집속하는 위상 검출 에코가 연속 에코인MRI 시스템.
- 제 9 항에 있어서,이미징 펄스 시퀀스에 의해서 이미징 데이터를 수집하는 스캔과는 다른 참조 스캔 동안에 위상 검출 펄스 시퀀스에 의해서 위상 검출 데이터를 수집하는MRI 시스템.
- 제 9 항에 있어서,펄스는 이미징 펄스 시퀀스에 일치하기 전에 위상 검출 펄스 시퀀스에 따라 인가되는MRI 시스템.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001328447A JP3720752B2 (ja) | 2001-10-26 | 2001-10-26 | 0次位相検出方法およびmri装置 |
JPJP-P-2001-00328447 | 2001-10-26 |
Publications (2)
Publication Number | Publication Date |
---|---|
KR20030035963A KR20030035963A (ko) | 2003-05-09 |
KR100677017B1 true KR100677017B1 (ko) | 2007-01-31 |
Family
ID=19144505
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
KR1020020065514A KR100677017B1 (ko) | 2001-10-26 | 2002-10-25 | 0차 위상 검출 방법 및 mri 시스템 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20030083571A1 (ko) |
JP (1) | JP3720752B2 (ko) |
KR (1) | KR100677017B1 (ko) |
CN (1) | CN1231180C (ko) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006075380A (ja) * | 2004-09-10 | 2006-03-23 | Ge Medical Systems Global Technology Co Llc | Mr画像生成方法およびmri装置 |
JP4192139B2 (ja) * | 2004-10-27 | 2008-12-03 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ssfpパルスシーケンスにおけるフェーズサイクリング法および磁気共鳴撮影装置 |
US20120274322A1 (en) * | 2011-04-27 | 2012-11-01 | Sangwoo Lee | Magnetic resonance imaging apparatus |
CN111239657B (zh) * | 2020-01-20 | 2022-05-06 | 上海东软医疗科技有限公司 | 谱图的相位校正方法、装置及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05245123A (ja) * | 1992-03-03 | 1993-09-24 | Toshiba Corp | 磁気共鳴イメージング装置 |
JPH0795972A (ja) * | 1993-06-01 | 1995-04-11 | Siemens Ag | 核スピン断層撮影装置 |
KR960033403A (ko) * | 1995-03-28 | 1996-10-22 | 제이 엘. 채스킨 | 자기 공명 영상 장치 |
KR980012936A (ko) * | 1996-07-11 | 1998-04-30 | 제이 엘. 채스킨 | 위상 시프트 측정 방법과 위상 시프트 보정 방법과 자기 공명 촬상 장치 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL9002842A (nl) * | 1990-12-21 | 1992-07-16 | Philips Nv | Magnetische resonantie werkwijze en inrichting ter reductie van beeldfouten in een magnetisch resonantie beeld. |
US5825185A (en) * | 1996-11-27 | 1998-10-20 | Picker International, Inc. | Method for magnetic resonance spin echo scan calibration and reconstruction |
US5923168A (en) * | 1997-06-17 | 1999-07-13 | General Electric Company | Correction of artifacts caused by Maxwell terms in slice offset echo planar imaging |
US6275038B1 (en) * | 1999-03-10 | 2001-08-14 | Paul R. Harvey | Real time magnetic field mapping using MRI |
US6259250B1 (en) * | 1999-04-28 | 2001-07-10 | General Electric Company | Method and apparatus for reducing artifacts in echo planar imaging |
US6285187B1 (en) * | 1999-04-28 | 2001-09-04 | General Electric Company | Method and apparatus for reducing artifacts in echo planar imaging |
DE10015068C2 (de) * | 2000-03-25 | 2002-06-27 | Bruker Medical Gmbh | Verfahren zur Erzeugung von Bildern der magnetischen Resonanz |
US6700374B1 (en) * | 2000-03-29 | 2004-03-02 | Koninklijke Philips Electronics, N.V. | EPI calibration method to minimize ghosting in reconstructed images |
US6320380B1 (en) * | 2000-10-03 | 2001-11-20 | Marconi Medical Systems, Inc. | MRI method and apparatus for increasing the efficiency of echo lanar imaging and other late echo techniques |
US6933720B2 (en) * | 2001-12-11 | 2005-08-23 | Toshiba America Mri, Inc. | Sequence preconditioning for ultra-fast magnetic resonance imaging |
-
2001
- 2001-10-26 JP JP2001328447A patent/JP3720752B2/ja not_active Expired - Fee Related
-
2002
- 2002-10-24 US US10/279,803 patent/US20030083571A1/en not_active Abandoned
- 2002-10-25 CN CNB021470553A patent/CN1231180C/zh not_active Expired - Fee Related
- 2002-10-25 KR KR1020020065514A patent/KR100677017B1/ko not_active IP Right Cessation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05245123A (ja) * | 1992-03-03 | 1993-09-24 | Toshiba Corp | 磁気共鳴イメージング装置 |
JPH0795972A (ja) * | 1993-06-01 | 1995-04-11 | Siemens Ag | 核スピン断層撮影装置 |
KR960033403A (ko) * | 1995-03-28 | 1996-10-22 | 제이 엘. 채스킨 | 자기 공명 영상 장치 |
KR980012936A (ko) * | 1996-07-11 | 1998-04-30 | 제이 엘. 채스킨 | 위상 시프트 측정 방법과 위상 시프트 보정 방법과 자기 공명 촬상 장치 |
Also Published As
Publication number | Publication date |
---|---|
KR20030035963A (ko) | 2003-05-09 |
JP2003135418A (ja) | 2003-05-13 |
JP3720752B2 (ja) | 2005-11-30 |
CN1231180C (zh) | 2005-12-14 |
US20030083571A1 (en) | 2003-05-01 |
CN1413555A (zh) | 2003-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1635184B1 (en) | MR image production method and MRI apparatus | |
US6489765B2 (en) | Magnetic field variation measuring method and magnetic field variation compensating method for MRI apparatus, and MRI apparatus | |
EP1391746B1 (en) | Parallel magnetic resonance imaging using navigator echos | |
US6512372B1 (en) | MRI apparatus which has reduced ghost artifacts due to Maxwell term phase error | |
JPH07323021A (ja) | Mrイメージング装置 | |
US6836113B2 (en) | Measurement and correction of gradient-induced cross-term magnetic fields in an EPI sequence | |
US6737865B2 (en) | MRI apparatus | |
EP1083438A2 (en) | MR imaging method and MRI apparatus | |
KR19990037474A (ko) | 멀티 슬라이스 자기 공명 투영 방법 및 장치 | |
EP0835454B1 (en) | Echo planar imaging involving echo time shifting and motion compensation | |
US4954779A (en) | Correction for eddy current caused phase degradation | |
US6470203B2 (en) | MR imaging method, phase error measuring method and MRI system | |
KR100677017B1 (ko) | 0차 위상 검출 방법 및 mri 시스템 | |
US6573719B2 (en) | MR method for generating MR signal corresponding to k-space excitations along mutually offset trajectories | |
US5685304A (en) | Diffusion sensitizing imaging method and MRI apparatus | |
US4760339A (en) | NMR imaging method | |
US20020050816A1 (en) | MR imaging method, phase error measuring method, and MRI apparatus | |
US20050036944A1 (en) | Diffusion-weighted parallel imaging with navigator-signal-based phase correction | |
EP1126286A1 (en) | Compensation for the influence of eddy currents and residual magnetization in fast spin-echo MR imaging | |
JP3576641B2 (ja) | Mri装置 | |
JP4795565B2 (ja) | Mrデータ収集方法およびmri装置 | |
JP2899821B2 (ja) | Mriのrf自動調整法 | |
JP3542650B2 (ja) | Mri装置 | |
JP2950165B2 (ja) | 核磁気共鳴イメージング装置 | |
JP2006116216A (ja) | 磁気共鳴イメージング装置および画像補正評価方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A201 | Request for examination | ||
E902 | Notification of reason for refusal | ||
E601 | Decision to refuse application | ||
J201 | Request for trial against refusal decision | ||
J301 | Trial decision |
Free format text: TRIAL DECISION FOR APPEAL AGAINST DECISION TO DECLINE REFUSAL REQUESTED 20050511 Effective date: 20061130 |
|
S901 | Examination by remand of revocation | ||
GRNO | Decision to grant (after opposition) | ||
GRNT | Written decision to grant | ||
FPAY | Annual fee payment |
Payment date: 20120109 Year of fee payment: 6 |
|
FPAY | Annual fee payment |
Payment date: 20130109 Year of fee payment: 7 |
|
LAPS | Lapse due to unpaid annual fee |