KR20140040023A - 회전 코딩 그래디언트를 이용한 ｍｒ 위상 대조 혈관 조영 검사 - Google Patents

회전 코딩 그래디언트를 이용한 ｍｒ 위상 대조 혈관 조영 검사 Download PDF

Info

Publication number
KR20140040023A
KR20140040023A KR1020130112167A KR20130112167A KR20140040023A KR 20140040023 A KR20140040023 A KR 20140040023A KR 1020130112167 A KR1020130112167 A KR 1020130112167A KR 20130112167 A KR20130112167 A KR 20130112167A KR 20140040023 A KR20140040023 A KR 20140040023A
Authority
KR
South Korea
Prior art keywords
gradient
gradients
read
coding
switched
Prior art date
Application number
KR1020130112167A
Other languages
English (en)
Other versions
KR101595354B1 (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 지멘스 악티엔게젤샤프트
Publication of KR20140040023A publication Critical patent/KR20140040023A/ko
Application granted granted Critical
Publication of KR101595354B1 publication Critical patent/KR101595354B1/ko

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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/5635—Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • G01R33/4826—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory in three dimensions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56308—Characterization of motion or flow; Dynamic imaging
    • G01R33/56316—Characterization of motion or flow; Dynamic imaging involving phase contrast techniques

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Vascular Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

본 발명은 속도-의존 위상 정보가 추가의 바이폴라 코딩 그래디언트들(30)의 스위칭을 통해서 검사 대상 내의 이동 스핀들에 새겨지는 검사 대상의 MR 위상 대조 혈관 조영 검사 영상을 생성하기 위한 방법에 관한 것이며, 여기서, MR 위상 대조 혈관 조영 검사 영상의 생성을 위해, 검사 대상의 MR 신호들은 판독 그래디언트(20) 동안 논-카테시안 획득 패턴을 갖는 미가공 데이터 공간에서 판독되고; 상기 추가의 바이폴라 코딩 그래디언트들(30)은 이들이 논-카테시안 획득 패턴에 대응하는 좌표계를 따라서 스위칭되도록 스위칭되고, 이 좌표계의 한 좌표 축이 판독 그래디언트(20)를 따라서 이동한다.

Description

회전 코딩 그래디언트를 이용한 ＭＲ 위상 대조 혈관 조영 검사{MR PHASE CONTRAST ANGIOGRAPHY WITH ROTATING CODING GRADIENTS}
본 발명은 검사 대상의 MR 위상 대조 혈관 조영 검사 영상을 생성하는 방법 - 이 방법에서는 속도-의존 위상 정보가 추가의 바이폴라 코딩 그래디언트(bipolar coding gradients)를 스위칭함으로써 검사 대상에서의 이동 스핀(moving spins)에 새겨진다 - ; 및 이를 위한 MR 시스템에 관한 것이다.
자기 공명 혈관 조영 검사는 검사 대상의 혈관계의 MR 영상을 생성한다. 특히, 2 가지의 혈관 조영 검사 기술이 여기 공개된다. 첫 번째 기술은 유동 혈(flowing blood)과 정지 조직 간에 존재하는 신호 포화에 있어서의 차이들이 이용되는 타임 오브 플라이트(Time of Flight, TOF) 효과로 알려져 있는 것을 기반으로 한다. MR 혈관 조영 검사의 다른 기술은 안정 스핀(resting spins)의 위상 정보와는 다른 위상 정보가 이동 스핀에 새겨진다는 사실에 기반을 두고 있다. 이를 위해, 일반적으로 2 개의 MR 영상이 이러한 위상 대조 혈관 조영 검사 기술에서 획득되는데, 한 MR 영상은 추가(대부분 바이폴라) 코딩 그래디언트가 없는 영상이고 한 영상은 추가의 바이폴라 코딩 그래디언트의 스위칭을 이용한 영상이다. 두 영상 간의 위상 차를 계산함으로써 또는 두 영상의 콤플렉스 차(complex difference) 계산을 통해서, 추가의 바이폴라 코딩 그래디언트가 스위칭된 방향을 따른 위상 차로부터 속도 정보를 얻을 수 있다.
3개의 서로 다른 공간 방향의 유동 정보가 위상 대조 혈관 조영 검사의 도움으로 생성될 수 있는 경우의 영상 시퀀스가 도 1에 개략적으로 간단히 도시되어 있다. RF 펄스(10)가 슬라이스-선택 그래디언트(slice-selection gradient)(11) 동안 슬라이스 내의 스핀들을 여기시키기 위해 스위칭된다. 다른 영상 시퀀스에서와 같이, 위상 코딩 그래디언트(12) 및 판독 그래디언트(13)가 스위칭되고, 판독 그래디언트 동안 기간(14) 동안 신호 판독이 이루어진다. RF 펄스(10)의 스위칭과 생성된 신호 에코(echo) 사이의 기간은 에코 시간 TE라 불린다. 추가의 바이폴라 코딩 그래디언트(15, 16 및 17)는 이제 공간 방향 X, Y 및 Z를 따라서 유동 정보를 각각 얻기 위해서 스위칭될 수 있다. 도 1을 보면 알 수 있듯이, 3-차원 유동 정보를 얻기 위해서는 통상 4개의 측정이 필요하다: 추가의 바이폴라 코딩 그래디언트의 스위칭이 없는 하나의 레퍼런스 측정, 및 3개의 공간 방향 각각에서의 추가 코딩 그래디언트의 스위칭을 이용한 각각의 측정. 바이폴라 코딩 그래디언트에 의해 생성된 그래디언트 모멘트(gradient moment)는 코딩될 최대 속도에 의해 설정된다.
추가의 코딩 그래디언트가 없는 베이스 시퀀스에 관해서, 이는 초기에, 판독 그래디언트(13)가 위상 코딩 그랜디언트(12) 직후에 스위칭될 수 없기 때문에, 최소 에코 시간을 증가시키며; 오히려 바이폴라 코딩 그랜디언트의 스위칭을 위한 기간이 추가로 제공되어야만 한다. 이들 바이폴라 코딩 그래디언트(Venc 그래디언트라 불리기도 함)는 에코 시간을 최소화하기 위하여 하부 영상 시퀀스에 이용된 그래디언트에 일시적으로 겹쳐질 수 있다. 측정 시간을 단축하는 것 외에도, (예를 들어) 인트라복셀 영위상화(intravoxel dephasing)의 원치 않는 효과(즉, 서로 다른 속도 성분들의 상쇄적인 중첩(destructive superposition))가 감소하고 이에 따라서 T*2 감쇠(decay)에 기인한 신호 손실이 감소하기 때문에 품질에 관해서 추가의 장점이 얻어진다.
영상 그래디언트와 바이폴라 코딩 그래디언트의 중첩은 도 1에 도시된 바와 같이 카테시안(Cartesian) k-공간 샘플링의 분야에서 실현될 수 있다. 언급한 바와 같이, 네 개의 데이터 세트가 통상적으로 생성된다: 추가의 바이폴라 그래디언트 모멘트가 없는 유동-보상 그래디언트 스킴(scheme)을 이용한 한 데이터 세트, 및 공간적으로 서로 직교하는 각 바이폴라 코딩 그래디언트를 이용한 3개의 추가 데이터 세트.
도 2는 2차원 경우에 X-방향 및 Y-방향으로 유동 정보를 어떻게 얻는지를 보여주고 있고, 판독 방향은 항상 X-방향이다. 도 2의 좌측에서 알 수 있듯이, 주어진 카테시안 좌표에서 판독 방향은 항상 X-방향으로 이루어진다. 추가의 바이폴라 위상 코딩 그래디언트(17)는 X-방향으로 1회 스위칭되어 유동 코딩이 이 방향으로 이루어지고, Y-방향으로 1회 스위칭되어 유동 코딩이 이 방향으로 이루어진다. 제3 측정은 추가의 유동 코딩 그래디언트 없이 이루어진다.
논-카테시안 k-공간 궤도로 MR 영상에 속하는 미가공 데이터 공간 또는 각 k-공간을 획득하는 것도 알려져 있다. 특히, 언더샘플링(undersampling)을 통해서 보다 큰 시간 효율을 얻을 수 있다.
3-차원 방사상 k-공간 샘플링은 US 6,188,922 B1에 설명되어 있는 바와 같이 벡터 유동 코딩으로 효과적인 3D 측정을 성취하기 위해 위에 설명된 속도 코딩과 결합될 수 있다.
그러한 획득 스킴이 도 3에 개략적으로 도시되어 있다. 3개의 카테시안 공간 방향을 따른 코딩 그래디언트(15-17)는 유지되고; 그러나, 데이터 판독의 그래디언트(18, 19, 20)는 원하는 (예를 들어, 방사상) k-공간 샘플링이 성취될 수 있게 스위칭된다. 그래디언트(18, 19, 20)에 대해 각각 나타내어지는 개별 그래디언트들의 변형은 각각의 방사상 k-공간 투영을 위한 개별 그래디언트의 변형을 보여주어야 한다. 도 4에서, 이는 2차원 경우에 대해서 도시되어 있다. 궤도(8)(하단 왼쪽에서 상단 오른쪽으로 도시된)를 이용한 방사상 판독의 경우, 그래디언트 스위칭은 투영(8)이 성취될 수 있게 그래디언트(19 및 20)의 x-방향 및 y-방향으로 이루어진다. 그러나, 도 2 및 3에 도시된 바와 같이, 추가의 바이폴라 코딩 그래디언트는 또한 추가의 바이폴라 그래디언트(16 및 17)의 스위칭을 통해서 고정된 공간 방향 x 및 y로 생긴다. 도 3 및 도 4로부터 알 수 있듯이, 추가의 바이폴라 코딩 그래디언트는 공간적으로 고정된, 물리적인 xyz-좌표계에서 스위칭되고, 이와 대조적으로 3-차원 방사상 샘플링의 경우 영상 시퀀스는 매 판독마다 입체각으로 회전하는 판독 그래디언트를 포함한다. 각 투영의 경우, 양쪽 레퍼런스 시스템이 다른 그래디언트 커브는 회전 판독 그래디언트와 공간적으로 고정된 유동 코딩 그래디언트와의 일시적인 중첩으로 나타난다. 이는 완전한 k-공간 샘플링을 위해 이용되는 전체 그래디언트 트레인들 중에서, 단지 하나만이 TE의 최소화에 따라 최적화될 수 있고; 다른 모든 커브들은 시간-최적화 방식으로 실현될 수 없다. 그러나, 에코 시간 TE는 모든 측정된 k-공간 투영을 위한 글로벌 측정 파라미터로서 일정(constant)하다. 이는 필연적으로 시간-최적화 중첩에 비해서 TE 시간의 확장을 이끌어낸다.
이는 도 5 내지 도 8을 참조해서 설명된다. 도 5에서, 추가의 바이폴라 코딩 그래디언트는 쇄선(17)으로 표시되어 있다. 방사상 판독의 경우, Gread의 좌표축을 따라서 마찬가지로 나타나는 k-공간 궤도가 한 번 생긴다. 이러한 판독 그래디언트는 도 5에 20으로 표시되어 있다. 바이폴라 코딩 그래디언트(17) 및 판독 그래디언트는 이제 동일한 극성을 가지며 시간 최소화가 이들 두 그래디언트의 중첩에 의해서 시도되고(이는 도 5에 화살표로 상징된 바와 같다), 그 결과는 도 6의 이미지와 같고, 여기서 바이폴라 코딩 그래디언트(17)와 판독 그래디언트(20)는 중첩되어 있다. 이는, 도 7에 도시된 바와 같이, 그래디언트 스위칭이 그래디언트(21A 및 21B)로 표현되어 있는 바와 같이 필요하다는 것을 의미한다. 그래디언트(21A 및 21B)는 판독 그래디언트(20)의 유동-보상 프리-그래디언트와 바이폴라 코딩 그래디언트(17)의 구조적인 중첩에 대응한다. 특히 그래디언트(21A 및 21B)로부터 명백하듯이, 에코 시간을 단축하기 위해서는 전반적으로 매우 높은 그래디언트 세기는 물론이고 이 그래디언트의 고속 회전율(fast slew rate)이 필요하다. 이것이 기존 그래디언트 시스템 때문에 가능하지 않거나 그러한 그래디언트 스위칭이 다른 이유들로 필요하지 않다면, 도 8에 도시된 바와 같이, 이를 단지 그래디언트의 확장된 스위칭으로 대체하여, 그래디언트(21A 및 21B)를 통해서 도 5에 도시된 것과 동일한 그래디언트 모멘트를 성취할 수 있다. 이는 도 7의 그래디언트(21A 및 21B)와 동일한 그래디언트 모멘트를 생성하기 위해 그래디언트(22)가 보다 긴 기간 동안 스위칭되어야만 하기 때문에 전체 에코 시간 TE가 확장된다는 것을 의미한다. 그러나, 이는 에코 시간을 확장시키는 결과를 낳는다. 이때, 투영에서 나타나는 이러한 확장된 에코 시간은 에코 시간이 모든 투영에 대해 일정하게 유지되어야만 하기 때문에 모든 다른 투영에도 이용되어야만 한다.
이는, 종래 기술에서, 전체 에코 시간을 늘리는 것이 수용되어 왔음을 의미한다.
그러므로, 본 발명의 목적은 k-공간의 논-카테시안 신호 판독의 경우에도 에코 시간의 최소화가 가능한 위상 대조 혈관 조영 검사를 실행할 수 있게 하는 것이다.
이 목적은 독립 청구항들의 특징들로 성취된다. 또한 실시예들은 종속 청구항들에 설명되어 있다.
본 발명의 제1 양태에 따르면, 속도-의존 위상 정보가 추가의 바이폴라 코딩 그래디언트들의 스위칭을 통해서 검사 대상 내의 이동 스핀들에 새겨지는 검사 대상의 MR 위상 대조 혈관 조영 검사 영상을 생성하기 위한 방법이 제공된다. MR 위상 대조 혈관 조영 검사 영상의 생성을 위해, 검사 대상의 MR 신호들은 판독 그래디언트 동안 논-카테시안 획득 패턴으로 미가공 데이터 공간에서 판독된다. 본 발명에 따르면, 추가의 바이폴라 코딩 그래디언트들은 이들이 논-카테시안 획득 패턴에 대응하는 좌표계를 따라서 스위칭되도록 스위칭되고, 이 좌표계의 한 좌표 축은 판독 그래디언트를 따라서 이동한다(travel). 이는 추가의 바이폴라 코딩 그래디언트들이 공간적으로 고정된 xyz-좌표계가 아니라 오히려 판독 방향으로 회전하는 좌표계에서 스위칭된다는 것을 의미한다. 그럼으로써, 에코 시간 TE의 단축 - 및 그에 따른 획득 시간의 단축 - 이 판독 그래디언트와 바이폴라 코딩 그래디언트들의 적절한 중첩을 통해서 성취될 수 있다. 용어 "위상 대조 혈관 조영 검사"는 순수한 혈관 조영 표현 및, 적용가능할 경우, 그에 의해서 나타나는 유동 수량화(flow quantification), 즉 유동의 양적 검출을 포함할 수 있다.
한 실시예에서, 판독 그래디언트를 따라서 이어지는 추가의 바이폴라 코딩 그래디언트들 중 하나는 판독 그래디언트에 관해서 반대 극성을 가지고 이동할 수 있다. 이는 전반적으로 더 작은 그래디언트 모멘트가 판독 그래디언트에 관해서 추가의 바이폴라 코딩 그래디언트들의 대향-극성 스위칭을 통해 스위칭되고, 그럼으로써 그래디언트 상승 시간과 전체 레벨에 부과되는 필요조건이 작아져서 에코 시간이 전반적으로 단축될 수 있다는 것을 의미한다.
판독 그래디언트를 따라 이동하는 추가의 바이폴라 코딩 그래디언트는 유리하게도 판독 그래디언트와 시간상 오버랩되게 스위칭된다. 가능하다면, 판독 그래디언트를 따라 이동하는 추가의 바이폴라 코딩 그래디언트는 근본적으로 판독 그래디언트와 동시에 스위칭된다. 이는 판독 그래디언트와 반대인 극성으로 이동하기 때문에, 동시 스위칭이 이루어지는 경우, 최종 그래디언트 모멘트(이는 스위칭된 그래디언트 아래 있는 영역에 대응함)는 추가의 바이폴라 코딩 그래디언트에 의해 생성된 그래디언트 모멘트와 판독 그래디언트의 프리-그래디언트들에 의해 생성된 그래디언트 모멘트가 상쇄되는 경우 전반적으로 매우 작을 수 있다(심지어는 0). 그럼으로써 에코 시간은 더 최적화된다. 판독 그래디언트로 이동하는 죄표계에서는, 스위칭될 바이폴라 코딩 그래디언트의 상대적인 각 배열의 자유도(degree of freedom)가 더 있다. 2차원 경우에, 이들은 서로 직교인 2개의 바이폴라 코딩 그래디언트이고; 3차원 경우에, 이들은 서로 직교인 3개의 바이폴라 코딩 그래디언트이다. 그럼으로써 이 성분은 판독 그래디언트와 반대의 극성으로 이동하고 유리하게도 판독 그래디언트와 시간상 오버랩되도록 각각 판독 그래디언트를 따라서 스위칭된다.
이동 스핀들의 위상 대조 혈관 조영 검사의 경우, 이동 스핀들의 제1 유동 정보는 논-카테시안(2D 또는 3D) 좌표계에서도 판정될 수 있다. 나중에 카테시안 좌표계를 따른 유동 성분들을 계산하기 위해 논-카테시안 좌표계로부터 판정된 유동 정보를 카테시안 좌표계로 변환하는 것이 실행될 수 있다. 계산된 유동 성분들은 이후 카테시안 좌표계를 따라서 표현될 수 있고, 그럼으로써, 유동 정보는, 예를 들어, 유동 속도의 크기 및 방향을 나타내는 벡터로 벡터 맵을 표현함으로써, 뷰어에게 종래 방식으로 표현될 수 있다.
더욱이, 위상 대조 혈관 조영 검사 영상의 판독 그래디언트 또는 획득이 미가공 데이터 공간을 채워서 판독 그래디언트들의 시퀀스에 의해 표시되는 노출(exposure)들을 줄이거나 또는 최소화시키는 것이 가능하다. 판독 그래디언트들의 시퀀스는 가능한 와전류(eddy current)를 최소화하기 위해 바이폴라 코딩 그래디언트들이 연속 방향 변경 또는 다른 적합한 방향 변경(예를 들어, 쌍으로)을 겪을 수 있게 선택될 수 있다. 예를 들어, 구형의(spherical) 미가공 데이터 공간이 3차원 경우에 방사상 투영들로 스캔(scan)되면, 하나의 투영에서 다른 투영으로의 큰 점프(jump)가 방지되는데, 왜냐하면 이들 큰 점프들은 증가된 와전류를 의미하는 큰 그래디언트 변경들에 대응하기 때문이다.
방사상 미가공 데이터 공간 샘플링 이외에도, 다른 논-카테시안 샘플링들, 예를 들어, 나선형-스캔, 로제트(rosettes) 또는 블레이드(blade) 또는 TWIRL 스캔이 가능하다.
위상 코딩 그래디언트들에 대한 규정된 그래디언트 모멘트로 스위칭되는 2개의 MR 데이터 세트들 간의 위상 차와, 바이폴라 코딩 그래디언트가 전혀 스위칭되지 않은 데이터 세트의 계산 대신에, 다른 바이폴라 그래디언트들을 이용하는 것, 예를 들어, 바이폴라 코딩 그래디언트 Gmax 및 0 대신에 그래디언트 쌍
Figure pat00001
및
Figure pat00002
를 한 번 이용하는 것도 가능하다. 3차원 경우에, 예를 들어, 바이폴라 코딩 그래디언트들에 의해 생성된 그래디언트 모멘트들의 4면체 배열(tetrahedral arrangement)이 이용될 수 있다.
더욱이, 다양한 속도를 코딩하기 위해서 다양한 그래디언트 세기들로 좌표축을 따른 공간 방향으로 추가의 바이폴라 코딩 그래디언트들을 스위칭하는 것도 가능한다. 이러한 기법 - 멀티-벤크 획득 기법(multi-venc acquisition technique)이라고도 알려져 있음 - 에 있어서, -180°로부터 +180°까지의 위상 공간에서 너무 높은 속도들의 에일리어싱(aliasing)을 방지하기 위해서 또는 또한 완벽한 푸리에 코딩을 가능하게 하기 위해서 임의 많은 속도 감도 및 결합이 코딩된다. 다양한 코딩 그래디언트 세기들이 또한 규정된 공간 방향에서의 에일리어싱 아티팩트(aliasing artifacts)를 방지하기 위하여 다양한 공간 방향에서 이용될 수 있다.
본 발명은 또한 이를 위한 MR 시스템에 관한 것이다.
본 발명은 첨부 도면을 참조하여 이하 상세히 설명된다.
도 1은 종래 기술에 따른 3개의 공간 방향으로 추가의 바이폴라 코딩 그래디언트들 갖는 그래디언트 스킴을 개략적으로 보여주고 있다.
도 2는 종래 기술에 따른 2차원 경우에 대한 유동 코딩 및 판독 방향을 보여주고 있다.
도 3은 종래 기술에 따른 미가공 데이터 공간의 논-카테시안 판독의 경우 카테시안 코딩 그래디언트들의 이용을 보여주고 있다.
도 4는 2차원 경우의 도 3의 예이다.
도 5는 도 2의 예시적인 실시예에서 에코 시간을 최소화하기 위한 판독 그래디언트와 코딩 그래디언트의 중첩을 개략 보여주고 있다.
도 6은 도 5의 완전한 중첩이 발생한 것을 보여주고 있다.
도 7은 결과로 생긴 도 6으로부터의 필요한 그래디언트의 스위칭을 보여주고 있다.
도 8은 도 7에서와 같은 동일한 그래디언트 모멘트를 구하기 위한 에코 시간의 최종 확장을 보여주고 있다.
도 9는 본 발명에 따른 추가의 코딩 그래디언트 - 이 추가의 코딩 그래디언트는 동일한 논-카테시안 좌표계에서 판독 그래디언트와 함께 이동한다 - 의 스위칭을 보여주고 있다.
도 10은 본 발명의 2차원 경우의 유동 코딩 및 판독을 보여주고 있다.
도 11은 에코 시간을 최소화하기 위한 코딩 그래디언트와 판독 그래디언트의 대향-극성 스위칭을 보여주고 있다.
도 12는 카테시안 공간 방향의 유동 정보로 개별 유동 영상을 재구성하기 위한 단계들의 흐름도이다.
도 13은 미가공 데이터 공간의 방사상 대칭 3차원 획득의 경우 카테시안 좌표계로의 유동 정보의 변환을 보여주고 있다.
도 14는 본 발명에 따른 추가 실시예를 보여주고 있다.
도 15는 본 발명에 따라 판독 그래디언트와 코딩 그래디언트의 제어가 가능한 본 발명에 따른 MR 시스템을 개략 보여주고 있다.
도 15에는 위상 대조 혈과 조영 검사의 시간-최적화된 획득이 가능한 MR 시스템(100)이 개략적으로 도시되어 있다. MR 시스템은 분극 전계 B0를 생성하기 위한 자석(110)을 갖고 있다. 침대(111)에 누워있는 검사 대상(112)은 결과의 자화가 분극 전계 B0의 방향으로 나타나도록 자석에 배열되어 있다. 자계 그래디언트들의 스위칭에 의해서 자화가 되도록 RF 펄스를 스위칭하여 MR 신호를 검출하고 MR 영상을 생성하는 방법은 이 기술 분야에 숙련된 자에게 공지되어 있으므로 여기서 더 이상 상세히 설명하지 않는다. 위상 값의 차를 계산하여 MR 위상 대조 혈관 조영 검사 영상을 생성하는 방법도 역시 이 기술 분야에 숙련된 자에게 공지되어 있다. MR 시스템은 중앙 제어 시스템(120)을 갖고 있다. 그래디언트 유닛(121)은 위상 대조 MR 영상의 생성을 위해 필요한 그래디언트 자계를 생성한다. RF 유닛(122)은 방사되는 RF 전계를 생성한다. 영상 시퀀스 제어기(123)는 선택된 영상 시퀀스(imaging sequence)에 따라서 스위칭되는 자계 그래디언트 및 RF 펄스들의 시퀀스를 제어한다. 영상 컴퓨터(124)는 MR 영상의 계산 및 표현을 가능케 해준다. 입력 유닛(125) - HMI(Human Machine Interface)라고도 불림 - 은, 예를 들어, 영상 시퀀스의 선택, 영상 파라미터의 입력 등에 의해 MR 시스템의 사용자가 검사 대상의 검사 워크플로우를 제어할 수 있게 해준다. 다음에 상세히 설명되는 바와 같이, 영상 시퀀스 제어기(123)는 나중에 코딩 그래디언트들이 판독 그래디언트와 동일한 논-카테시안 좌표계에서 스위칭되도록 그래디언트 유닛(121)을 제어한다.
기본적으로, MR 시스템은 여기에는 상세히 제시되지 않았으며 명료성과 본 발명의 이해를 도모하기 위해 생략된 추가의 컴포넌트들을 갖고 있다. 더욱이, 유닛(121-125)은 제시된 분포와는 다르게 설계될 수 있으며; 개별 유닛에 의해 제어되는 기능들은 다른 구성에서 제어될 수도 있다. 더욱이, 도 11에 도시된 시스템들은 하드웨어 또는 소프트웨어 또는 하드웨어와 소프트웨어의 결합으로 형성될 수 있다.
본 발명은 도 9를 참조하여 상세히 설명된다. 본 발명에 따르면, 도 5 및 도 6에 제시된 바이폴라 코딩 그래디언트(도 9에는 바이폴라 코딩 그래디언트(30)가 쇄선으로 표시되어 있음)는 공간적으로 고정된 xyz-좌표계에서 스위칭되지 않으며; 오히려 이는 각 판독 그래디언트(20)의 판독 방향으로 이동한다. 도 9에 나타나는 표현은 회전 레퍼런스 시스템에서 나타나는데, 이는 판독 그래디언트 Gread가 미가공 데이터 공간의 방사상 판독의 경우 미가공 데이터 공간의 중간 점(middle point)을 통해서 이어지는 방사상 궤도를 따라서 나타난다는 것을 의미한다. 회전 판독 그래디언트(방사상 신호 판독을 고려함)의 경우에, 바이폴러 코딩 그래디언트도 회전한다. 추가의 바이폴라 코딩 그래디언트(20)는 이제 그것이 먼저 각 판독 그래디언트에 평행하게 상기 판독 그래디언트에 대향-극성, 즉 실제 신호 판독(14) 전에 놓여있는 2개의 그래디언트 모멘트에 대해 대향 극성으로 이동하도록 스위칭된다. 판독 방향은 미가공 데이터 공간을 통한 각 투영마다 변한다. 바이폴라 코딩 그래디언트들은 이제 각 투영에 따라 변하고 판독 방향에 맞춰 얼라인(align)된다. 바이폴라 코딩 그래디언트들은 판독 방향에 맞춰 얼라인되기 때문에, 항상 바이폴라 코딩 그래디언트들을 판독 그래디어트들에 대한 대향-극성으로 스위칭할 수 있다.
이는 미가공 데이터 공간의 2차원 방사상 판독의 다음 예에 제시되어 있다. 원의 중앙 점을 통해 이동하는 각 판독 궤도는 반경 r 및 각
Figure pat00003
에 의해 정의된다. 카테시안 xy-방향을 따른 바이폴러 코딩 그래디언트를 스위칭하는 대신에, 바이폴라 코딩 그래디언트는 이제 그 궤도와 같이 방사상으로 얼라인된다. 2개의 카테시안 공간 방향의 유동 정보를 판정하기 위해서 Gx- 및 Gy-방향으로 각 그래디언트를 스위칭하는 대신에, r-방향 및 p-방향의 2개의 속도 성분을 판정하기 위해서 r-방향과 이에 직교 방향 -
Figure pat00004
- 으로 바이폴라 코딩 그래디언트의 스위칭이 1회 이루어진다. 이들 성분은 이후 좌표 변환을 통해서 x-방향 및 y-방향의 속도 성분으로 변환된다(translate).
도 10에는, 영상 투영에 관해 유동 코딩이 어떻게 동작하는지가 도시되어 있다. 도 9는 회전 레퍼런스 시스템에서의 그래디언트(30)를 제시하고 있지만, 도 10에는, MR 시스템에 존재하는 그래디언트 시스템에 대한 이들 그래디언트의 분포는 x-방향 및 y-방향으로 나타난다. 판독 방향이 궤도(31)를 따라서 나타나면, 유동 코딩은 2차원에서 유동 정보의 전체를 얻기 위해 한 번은 판독 방향에 반대로 평행하게 그리고 두 번째는 이에 직교하여 발생한다. 도 10에서 알 수 있듯이, 유동 코딩은 영상 판독의 방향으로 나타나고 또는 영상 투영은 각각 판독 방향에 반대로 나타난다.
도 11에는 바이폴라 코딩 그래디언트가 판독 그래디언트의 프리-그래디언트들과 거의 동시에 스위칭될 때 최소 에코 시간이 성취되는 방법이 도시되어 있다. 스위칭되는 전반적인 그래디언트 모멘트는 기간(40) 동안 그래디언트들의 합이다. 도시된 예에서, 바이폴라 코딩 크래디언트(30)와 판독 그래디언트에 속하는 그래디언트(20)는 기간(40) 동안 동일한 그래디언트 모멘트를 가지고 있으며, 그러므로 기간(40) 동안 스위칭되는 전반적인 그래디언트 모멘트는 0이 될 것이다. 이는 추가의 바이폴라 코딩 그래디언트를 스위칭함이 없이, 단지 판독 그래디언트(20)가 1회 스위칭되고, 기간(40) 동안 바이폴라 코딩 그래디언트의 추가 스위칭이 있다면 그래디언트는 스위칭되지 않는 것을 의미한다. 당연히, 바이폴라 코딩 그래디언트의 그래디언트 모멘트(이는 상쇄될 속도에 의존함) 및 판독 그래디언트의 그래디언트 모멘트는 상쇄되지 않아야 하며; 오히려, 이들은 서로 다른 크기일 수 있어 상기 기간 동안 스위칭되는 최종 그래디언트는 0으로 되지 않을 것이다.
판독 방향을 따라 스위칭되는 바이폴라 코딩 그래디언트(30) 이외에도, 3차원 속도 정보를 위해 판독 방향에 각각 직교인 2개의 추가 바이폴라 코딩 그래디언트(도시되지 않음)가 판독 그래디언트와 함께 회전하는 좌표계에서 3차원 속도 정보를 얻기 위해 스위칭된다.
도 11에 도시된 바와 같이, 에코 시간은 바이폴라 코딩 그래디언트의 대수 부호의 최적 선택을 통해서 최적화될 수 있다. 미가공 데이터 공간의 샘플링은 방사상으로 이루어질 수 있고; 그러나, 다른 논-카테시안 미가공 데이터 공간 샘플링 스킴, 예를 들어, 나선-형 샘플링, 로제트-형 샘플링, 또는 순수한 방사상 샘플링에 관해 바깥쪽 k-공간에서 증가된 샘플링을 성취하기 위하여 실제 방사상 샘플링의 궤도들이 이들 궤도의 마지막에서 S-형으로 굽어지는 TWIRL 샘플링도 가능할 수 있다. 위상 대조 혈관 조영 검사를 위해, xyz 성분들은 판독 그래디언트와 회전하는 바이폴라 코딩 그래디언트(2차원 또는 3차원) 및 각 의존에 따라 코딩 그래디언트가 없는 레퍼런스 측정으로부터 계산될 수 있다(예로, 도 12의 흐름도에 도시된 바와 같이).
제1 단계 S1에서, 코딩 영상 데이터가 모든 관련 코딩 방향에 대해서 획득된다. 도 13에 도시된 바와 같이 미가공 데이터 공간의 구형 획득이 이루어지는 구형 좌표계에서, 바이폴라 코딩 그래디언트 방향들은, 예를 들어, r, p 및 s(
Figure pat00005
)이고, 지수 i는 획득된 모든 투영 및 방향 전체를 가로질러 이어진다. 마지막으로, 단계 S2에서, 단일 투영마다 구형 좌표계에서 코딩 성분 r, p 및 s로부터 xyz-성분이 계산된다. r,
Figure pat00006
,
Figure pat00007
구형 좌표계로부터 카테시안 좌표계로의 역변환은 도 13에 도시된 구형 좌표계로부터 카테시안 좌표계로의 변환으로 이루어진다. 단계 S1 및 S2는 반드시 차례로 진행되어야 할 필요는 없고; 이들은 인터리브(interleave)될 수도 있으며, 따라서 단계 S1이 단계 S2가 시작되기 전에 결론이 나야할 필요가 없다. 위상 대조 혈관 조영 검사로부터 알려져 있는 대로 xyz에 대한 속도-의존 위상 정보(여기서 코딩 그래디언트들은 xyz-방향으로 스위칭된다)는 이후 단계 S3에서 설명될 수 있다.
판독 투영들의 시퀀스는 가능한 와전류 효과를 최소화하기 위해 바이폴라 코딩 그래디언트들 자체가 연속적인 또는 다른 적합한 방향 변경(예를 들어 쌍으로 방향 변경)을 겪도록 전체가 선택될 수 있다. x-방향으로의 코딩 그래디언트, y-방향으로의 코딩 그래디언트 및 z-방향으로의 코딩 그래디언트가 없는 스킴에 따른 코딩 그래딩언트들의 스위칭 대신에, 이용되는 추가의 코딩 그래디언트 모멘트 없는 측정이 아니라 서로 직교인 3개의 그래디언트 세기를 이용한 측정이 이루어지는 다른 코딩 그래디언트 스킴들(예를 들어 코딩 그래디언트 모멘트의 4면체 배열)을 적용하는 것도 가능하다.
이는 도 14를 이용하여 2차원으로 설명된다. 데이터 판독은 도 10에 도시된 바와 같이 이루어진다. 그러나, 유동 코딩은 2개의 코딩이 서로 직교하고 측정이 코딩 없이 이루어지는 도 10과 같이 이루어지지 않으며; 오히려, 유동 코딩은 3개의 측정으로 이루어지며 여기서 유동 코딩들은 가능한 한 서로 멀리 떨어져 있게 위치해 있는 방향들로 향한다. 2D의 경우, 여기서 개별 유동 코딩들 간에 120° 각이 발생한다. 3D의 경우, 4면체(tetrahedron)가 발생할 것이다.
개별 유동 코딩 방향들에 대해서, 에일리어싱 아티팩트를 방지하기 위해서 또는, 예를 들어, 완전한 푸리에 코딩을 가능케 하기 위하여 임의 많은 코딩 감도 및 결합이 코딩될 수도 있다. 예를 들어, 3개의 서로 다른 속도 범위를 커버하기 위해 하나의 공간 방향에서 3개의 서로 다른 그래디언트 세기가 스위칭될 수 있다. 마찬가지로, 예를 들어, 관련 방향들에서 에일리어싱을 피하기 위해 모든 판독 그래디언트들이 작은 코딩 그래디언트로 측정되고 단지 하나의 서브 그룹이 보다 큰 코딩 그래디언트로 추가로 측정된다는 점에서 판독 그래디언트들의 다양한 서로 다른 서브-그룹들을 다양한 코딩 그래디언트들에 연관시키는 것도 가능하다.
본 발명은 선형 독립 유동 코딩들을 이용한 한 세트의 영상 획득이 본질적으로 전체 유동 정보를 담고 있다는 사실을 기반으로 하고 있다. 논-카테시안 획득 방법에 결합하여, 본 발명은 측정 시간을 최적화하는데, 즉 측정 시간을 최소화하는데 이용할 수 있다. 본 발명의 장점은 에코 시간 TE의 단축에 있다. 더욱이, 영상 재구성 이전에 개별 판독 방향들을 xyz-좌표로 변환하는 것을 통해, 바이폴라 코딩 그래디언트들이 카테시안 공간 방향을 따라서 스위칭되는 종래 기술에 따른 방법에서와 동일한 영상 데이터를 얻을 수가 있다.

Claims (9)

  1. 속도-의존 위상 정보가 추가의 바이폴라 코딩 그래디언트들(30)의 스위칭을 통해서 검사 대상 내의 이동 스핀들(moving spins)에 새겨지는 검사 대상의 MR 위상 대조 혈관 조영 검사 영상을 생성하기 위한 방법에 있어서, 상기 MR 위상 대조 혈관 조영 검사 영상의 생성을 위해, 상기 검사 대상의 MR 신호들은 판독 그래디언트(20) 동안 논-카테시안(non-Cartesian) 획득 패턴으로 미가공 데이터 공간에서 판독되고; 상기 추가의 바이폴라 코딩 그래디언트들(30)은 이들이 상기 논-카테시안 획득 패턴에 대응하는 좌표계를 따라서 스위칭되도록 스위칭되고; 이 좌표계의 좌표 축이 상기 판독 그래디언트(20)를 따라서 이동(travel)하는 방법.
  2. 제1항에 있어서,
    상기 판독 그래디언트를 따라서 이어지는 상기 추가의 바이폴라 코딩 그래디언트(30)는 상기 판독 그래디언트(20)에 대해 반대 극성으로 이동하는 것을 특징으로 하는 방법.
  3. 제1항 또는 제2항에 있어서,
    상기 판독 그래디언트를 따라 이동하는 상기 추가의 바이폴라 코딩 그래디언트(30)는 상기 판독 그래디언트(20)와 시간상 오버랩되게 스위칭되는 것을 특징으로 하는 방법.
  4. 제1항 내지 제3항 중 어느 한 항에 있어서,
    상기 이동 스핀들의 유동 수량화(flow quantification)를 위해,
    - 상기 논-카테시안 좌표계에서 상기 이동 스핀들의 제1 유동 정보를 판정하는 단계,
    - 상기 판정된 유동 정보를 카테시안 좌표계로 변환하여 상기 카테시안 좌표계를 따르는 유동 성분들을 계산하는 단계, 및
    - 상기 계산된 유동 성분들을 제공하는 단계
    가 실행되는 것을 특징으로 하는 방법.
  5. 제1항 내지 제4항 중 어느 한 항에 있어서,
    상기 판독 그래디언트들(20)의 시퀀스에 의해 유도된 와전류가 감소하도록 상기 위상 대조 혈관 조영 검사 영상의 획득을 위한 상기 판독 그래디언트들(20)이 상기 미가공 데이터 공간을 채우는 것을 특징으로 하는 방법.
  6. 제1항 내지 제5항 중 어느 한 항에 있어서,
    3차원 미가공 데이터 공간이 획득되는 것을 특징으로 하는 방법.
  7. 제6항에 있어서,
    상기 추가의 바이폴라 코딩 그래디언트들(30)은 상기 바이폴라 코딩 그래디언트들에 의해 생성된 그래디언트 모멘트들의 4면체 배열(tetrahedral arrangement)로 스위칭되는 것을 특징으로 하는 방법.
  8. 제1항 내지 제7항 중 어느 한 항에 있어서,
    상기 추가의 바이폴라 코딩 그래디언트들(30)은 다양한 속도의 코딩을 위해 다양한 그래디언트 세기들을 가지고 상기 좌표계의 좌표축을 따른 공간 방향에서 스위칭되는 것을 특징으로 하는 방법.
  9. 검사 대상의 MR 위상 대조 혈관 조영 검사 영상을 생성하기 위한 자기 공명 시스템(100)에 있어서, 이 자기 공명 시스템(100)은:
    - 상기 MR 위상 대조 혈관 조영 검사 영상을 생성하는데 필요한 그래디언트 전계를 생성하고, 추가의 바이폴라 코딩 그래디언트들(30) - 이들 그래디언트를 통해서 속도-의존 위상 정보가 상기 검사 대상의 이동 스핀들에 새겨진다 - 를 생성하는 그래디언트 유닛(121),
    상기 그래디언트 유닛을 제어하며 어느 획득 패턴으로 상기 위상 대조 혈관 조영 검사 영상에 속하는 미가공 데이터 공간이 획득되는지를 설정하는 영상 시퀀스 제어기(123)를 포함하고,
    상기 영상 시퀀스 제어 유닛(123)은, 상기 MR 위상 대조 혈관 조영 검사 영상의 생성을 위해, 미가공 데이터 공간 내의 상기 검사 대상의 MR 신호들이 판독 그래디언트 동안 논-카테시안 획득 패턴으로 판독되도록 상기 그래디언트 유닛(121)을 제어하고; 상기 추가의 바이폴라 코딩 그래디언트들은 이들이 상기 논-카테시안 획득 패턴에 대응하는 좌표계를 따라서 스위칭되도록 스위칭되고, 상기 좌표계의 좌표축이 상기 판독 그래디언트를 따라 이동하는 자기 공명 시스템.
KR1020130112167A 2012-09-25 2013-09-17 회전 코딩 그래디언트를 이용한 ｍｒ 위상 대조 혈관 조영 검사 KR101595354B1 (ko)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102012217227.0A DE102012217227B4 (de) 2012-09-25 2012-09-25 MR-Phasenkontrastangiographie mit rotierenden Kodierungsgradienten
DE102012217227.0 2012-09-25

Publications (2)

Publication Number Publication Date
KR20140040023A true KR20140040023A (ko) 2014-04-02
KR101595354B1 KR101595354B1 (ko) 2016-02-19

Family

ID=50235210

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020130112167A KR101595354B1 (ko) 2012-09-25 2013-09-17 회전 코딩 그래디언트를 이용한 ｍｒ 위상 대조 혈관 조영 검사

Country Status (4)

Country Link
US (1) US9500734B2 (ko)
KR (1) KR101595354B1 (ko)
CN (1) CN103675739B (ko)
DE (1) DE102012217227B4 (ko)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5611882B2 (ja) * 2010-05-31 2014-10-22 株式会社東芝 磁気共鳴イメージング装置
JP5925529B2 (ja) * 2011-03-31 2016-05-25 株式会社東芝 磁気共鳴イメージング装置
DE102013221938B4 (de) 2013-10-29 2018-11-08 Siemens Healthcare Gmbh Verfahren zur Aufnahme von Magnetresonanzdaten mit einer diffusionsgewichteten Magnetresonanzsequenz und Magnetresonanzeinrichtung
KR102449249B1 (ko) * 2015-05-27 2022-09-30 삼성전자주식회사 자기 공명 영상 장치 및 그 방법
DE102016201458B4 (de) * 2016-02-01 2018-10-31 Siemens Healthcare Gmbh Radiales Erfassen von MR-Daten zur Vermeidung von Artefakten
CN109752681B (zh) * 2019-01-22 2021-10-29 奥泰医疗系统有限责任公司 倾斜平面回波成像方法及核磁共振成像系统
CN110558986A (zh) * 2019-08-27 2019-12-13 苏州朗润医疗系统有限公司 适用于磁共振成像的相位对比成像方法
CN111257809B (zh) * 2020-01-23 2022-05-17 上海东软医疗科技有限公司 磁共振成像方法、装置、存储介质及医疗设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06121782A (ja) * 1991-02-07 1994-05-06 General Electric Co <Ge> 流体の流れの定量測定方法
JP2011036455A (ja) * 2009-08-12 2011-02-24 Hitachi Medical Corp 磁気共鳴イメージング装置
US20110140698A1 (en) * 2009-12-10 2011-06-16 Andreas Greiser Method and device for phase-sensitive flow measurement by magnetic resonance
US20110175610A1 (en) * 2010-01-18 2011-07-21 Case Western Reserve University Non-cartesian under-sampled multi-echo mri

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6188922B1 (en) 1999-01-08 2001-02-13 Wisconsin Alumni Research Foundation Phase contrast imaging using interleaved projection data
ATE404876T1 (de) 2001-11-12 2008-08-15 Wisconsin Alumni Res Found Dreidimensionale phasenkontrastmagnetresonanzabbildung mit verschachtelten projektions-rekonstruktions-daten
US7049816B2 (en) * 2004-09-30 2006-05-23 Wisconsin Alumni Research Foundation Magnetic resonance imaging with dual velocity encoded projection reconstruction acquisition
DE102005051323B4 (de) 2005-10-26 2008-01-31 Siemens Ag Verfahren und Gerät zur Magnet-Resonanz-Bildgebung auf Basis einer Gradientenechosequenz
US9201129B2 (en) 2006-09-13 2015-12-01 Kabushiki Kaisha Toshiba Magnetic-resonance image diagnostic apparatus and method of controlling the same
DE102007041826B4 (de) * 2007-09-03 2018-02-01 Siemens Healthcare Gmbh Verfahren zur Optimierung von angiographischen MR-Bildern
US7847545B2 (en) * 2008-04-11 2010-12-07 General Electric Company System and method for correcting flow velocity measurements in phase contrast imaging using magnetic field monitoring
CN101647699B (zh) * 2008-08-12 2011-08-10 株式会社东芝 磁共振成像装置及磁共振成像方法
JP2010063871A (ja) 2008-08-12 2010-03-25 Toshiba Corp 磁気共鳴イメージング装置
US8148984B2 (en) * 2008-10-03 2012-04-03 Wisconsin Alumni Research Foundation Method for magnitude constrained phase contrast magnetic resonance imaging
CN103597370B (zh) * 2011-04-21 2016-01-06 马普协会 空间编码的相位对比磁共振成像
US9075121B2 (en) * 2011-07-15 2015-07-07 Wisconsin Alumni Research Foundation System and method for rotating angle velocity encoding, phase contrast magnetic resonance imaging
US8744551B2 (en) * 2012-07-12 2014-06-03 Northshore University Healthsystem Method for non-contrast enhanced magnetic resonance angiography
US9081074B2 (en) * 2012-08-21 2015-07-14 Beth Israel Deaconess Medical Center, Inc. (Bidmc, Inc.) Method and apparatus for accelerated phase contrast magnetic resonance angiography and blood flow imaging
US9964621B2 (en) * 2013-10-01 2018-05-08 Beth Israel Deaconess Medical Center, Inc. Methods and apparatus for reducing scan time of phase contrast MRI
DE102013220012B4 (de) * 2013-10-02 2015-04-30 Siemens Aktiengesellschaft Bildrekonstruktion mit mehreren virtuellen Spulen

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06121782A (ja) * 1991-02-07 1994-05-06 General Electric Co <Ge> 流体の流れの定量測定方法
JP2011036455A (ja) * 2009-08-12 2011-02-24 Hitachi Medical Corp 磁気共鳴イメージング装置
US20110140698A1 (en) * 2009-12-10 2011-06-16 Andreas Greiser Method and device for phase-sensitive flow measurement by magnetic resonance
US20110175610A1 (en) * 2010-01-18 2011-07-21 Case Western Reserve University Non-cartesian under-sampled multi-echo mri

Also Published As

Publication number Publication date
CN103675739A (zh) 2014-03-26
CN103675739B (zh) 2016-08-24
DE102012217227A1 (de) 2014-03-27
US20140084916A1 (en) 2014-03-27
DE102012217227B4 (de) 2014-05-15
KR101595354B1 (ko) 2016-02-19
US9500734B2 (en) 2016-11-22

Similar Documents

Publication Publication Date Title
KR101595354B1 (ko) 회전 코딩 그래디언트를 이용한 ｍｒ 위상 대조 혈관 조영 검사
JP6998218B2 (ja) 動き検出を用いるｍｒ撮像
EP2815248B1 (en) Mri with motion correction using navigators acquired using a dixon technique
JP6513398B2 (ja) 事前情報に制約される正則化を使用するｍｒ画像再構成
EP2893363B1 (en) Propeller with dixon water fat separation
EP2699926B1 (en) Spatially encoded phase-contrast mri
JP6270709B2 (ja) 磁気共鳴生データの再構成方法および装置
US10254370B2 (en) Modified pulse sequence for magnetic resonance imaging using MRI-visible markers for motion correction
JP7075420B2 (ja) 可変コントラストのスタック・オブ・スター収集を使用したｍｒイメージング
CN109073720B (zh) 静音磁共振指纹识别
WO2004093682A1 (ja) 磁気共鳴イメージング方法及び装置
JP2019512083A (ja) 永続磁石配列を使用したｍｒｉ画像化システム
JP7350568B2 (ja) 磁気共鳴データ収集方法及び磁気共鳴イメージング装置
JP7023954B2 (ja) プロペラｍｒイメージング
JP2020522344A (ja) 並列マルチスライスｍｒ撮像
JP7399781B2 (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
CN114902062A (zh) 具有内部运动检测的3d mr成像
US20200300952A1 (en) Mr imaging using motion-dependent radial or spiral k-space sampling
JP2022521277A (ja) ウェーブ符号化を用いたパラレルｍｒイメージング

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
FPAY Annual fee payment

Payment date: 20200113

Year of fee payment: 5