KR101373563B1 - Ｔｏｆ-ｍｒａ를 이용한 혈류특성 및 ｍｒ-신호강도구배(전단율) 유도방법 - Google Patents

Ｔｏｆ-ｍｒａ를 이용한 혈류특성 및 ｍｒ-신호강도구배(전단율) 유도방법 Download PDF

Info

Publication number
KR101373563B1
KR101373563B1 KR1020120081035A KR20120081035A KR101373563B1 KR 101373563 B1 KR101373563 B1 KR 101373563B1 KR 1020120081035 A KR1020120081035 A KR 1020120081035A KR 20120081035 A KR20120081035 A KR 20120081035A KR 101373563 B1 KR101373563 B1 KR 101373563B1
Authority
KR
South Korea
Prior art keywords
image
blood vessel
blood flow
blood
signal intensity
Prior art date
Application number
KR1020120081035A
Other languages
English (en)
Other versions
KR20140014605A (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 KR1020120081035A priority Critical patent/KR101373563B1/ko
Priority to PCT/KR2013/006684 priority patent/WO2014017853A2/ko
Publication of KR20140014605A publication Critical patent/KR20140014605A/ko
Application granted granted Critical
Publication of KR101373563B1 publication Critical patent/KR101373563B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • 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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • 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/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • 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/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104—Vascular flow; Blood flow; Perfusion

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

TOF-MRA를 이용한 혈류특성 및 MR-신호강도구배 유도방법이 제공된다. 본 발명의 실시예에 따른 혈류특성 유도방법은, MRA 혈관 단면 영상을 획득하고 혈관 경계를 검출하여 혈관을 파악하며, 혈관을 의사-컬러링하여 혈류 영상으로 제공한다. 이에 의해, TOF-MRA를 이용한 비교적 간단한 알고리즘을 통해 혈류나 혈류역학 등의 혈류특성을 정확하고 매우 빠르게 유도하여, 이를 분석하여 혈류특성을 파악할 수 있도록 하고, 이를 혈관질환 진단 및 치료에 응용할 수 있게 된다.

Description

ＴＯＦ-ＭＲＡ를 이용한 혈류특성 및 ＭＲ-신호강도구배(전단율) 유도방법{Method of derivation for hemodynamics and MR-signal intensity gradient(or shear rate) using Time-Of-Flight - Magnetic Resonance Angiography}
본 발명은 혈류특성 유도방법에 관한 것으로, 더욱 상세하게는 TOF-MRA를 이용하여, 혈류나 혈류역학 등의 혈류특성을 유도하기 위한 방법 및 이를 적용한 컴퓨팅 장치에 관한 것이다.
전단 응력은 혈관 내피 세포에 영향을 미쳐 혈관기능을 유지시키는 가장 중요한 기계적 힘의 하나로 이를 분석하기 위한 방법들이 개발되어 왔다. CFD(Computational Fluid Dynamics)와 같은 최첨단 컴퓨터를 이용한 정밀분석이 가장 대표적인 방법이나, 이는 전체적인 과정이 매우 복잡하고 여러 단계를 거치며, 컴퓨터 시뮬레이션의 긴 과정을 거쳐야 한다. 또한, 개개의 혈관을 동시에 처리하기 힘들고, 각각 위의 과정들을 되풀이해야하는 어려움이 있어 연구영역 외에서는 범용 되지 못하고 있다.
최근에는 MRI기법을 다양하게 응용하여 이를 구현하려는 시도들이 계속되고 있으나, 이 역시 연구 단계에 멈춰있고 실제 사용되지는 못하고 있다.
벽 전단 응력을 구하기 위한 가장 큰 핵심과제는 혈관벽 근처에서의 속도구배(velocity gradient) 혹은 전단율(shear rate)을 정확히 얻는 것이다.
TOF-MRA(Time-Of-Flight - Magnetic Resonance Angiography)는, 도 1에 도시된 바와 같이, 자장에 완전 포화(saturation)된 영역에 포화되지 않은 혈액이 통과하면서 나타나는 혈류 흐름에 의한 증강효과를 영상화한 것이다.
이 신호 강도는 영상 용적 내에서 포화된 스핀의 이동 정도에 비례한다. 즉, TOF-MRA는 혈액의 혈류특성을 이용해 영상을 구성하게 된다. 기타 신호 강도에 영향을 미치는 인자들은 radiofrequency pulse의 flip angle과 Repetition time (TR), Echo time (TE) 등이다.
TOF-MRA는 뛰어난 영상해상도로 피검자 개개인의 다양한 혈관구조를 상세히 제시할 수 있어 현재 의료분야에서 가장 많이 활용되는 혈관 영상 기법이다. 또한, 관련기술이 계속 발달하여 와류와 같은 복합혈류현상도 진단할 수 있고, 최근에는 미세혈관 영상도 얻을 수 있는 등 그 발전속도가 빠르다.
phase contrast MRI는 TOF-MRA기법과는 다르게 혈류속도와 방향을 얻을 수 있는 영상기법으로, 혈류속도를 직접 측정할 수 있어, 이를 이용해 동맥경화성 플라그의 성장과 밀접한 관계를 갖는 벽전단응력 분포를 계측하는 연구가 수행되고 있다.
그러나, phase contrast MRI는 시간 및 공간 정밀도의 영향으로 시간과 비용적인 측면에서 큰 단점을 가지고 있다. 즉, 혈관구조를 전체적으로 상세히 제시할 수 없어 특정 영역에서 혈류속도를 얻기 위한 연구 외에는 활용할 수 없어 실제 임상에서는 거의 사용되지 못하고 있다.
TOF-MRA는 혈류의 증강효과에 의한 영상임에도 앞서 말한 radiofrequency pulse flip angle과 TR 등과 같은 영상 촬영 기법 혹은 조건에 따라 완벽한 혈류 특성을 얻지는 못한다.
이러한 이유로, TOF-MRA기법은 혈관의 전체적인 영상을 얻기에는 최적이나, 혈류속도와 같은 혈행 인자를 뽑아내기에는 적절치 않은 것으로 알려져 왔다. 특히, 중심부의 빠른 혈류 즉, laminar flow 중심부 최대속도를 보이는 영역이 TR 값에 따라서는 포화된 자장을 영상을 얻기도 전에 벗어나 버려 그 신호가 상쇄되기 때문이다.
TOF-MRA로 적절한 혈류속도를 얻지 못하는 또 다른 이유는 z축에 직각으로 설정된 자장의 방향으로 인해 촬영되는 혈관 단면의 flip angle이 z축에 직각이 아닌 사선 혹은 평행으로 배열되는 경우도 혈관 내부의 적절한 laminar 혈류 특성을 적절히 구현하기 힘들다.
그러나, 최근 TOF-MRA기법에서는 이러한 MR 영상의 제한점을 최대한 활용하여 촬영한다. TR을 오히려 더 줄여 중심속도보다는 주변 속도를 잘 살려낸 다음 이를 MIP(Maximum Intensity Projection) 시켜 3차원 혈관영상을 얻어내면 혈관의 좁아진 부위나 다양한 이상소견을 훨씬 더 효과적으로 진단할 수 있기 때문이다. 즉, 혈관영상의 초점이 중심속도를 정확히 잡아내는 것보다는 혈관 주변 영상(혈관벽주변의 속도임)을 최적화하여 얻는다.
현재 대부분의 TOF-MRA는 임상 및 연구현장에서 repetition time (TR), echo time (TE) 등 촬영 인자들을 모두 이러한 원리에 준해 셋팅해서 사용되고 있다.
본 발명은 상기와 같은 문제점을 해결하기 위하여 안출된 것으로서, 본 발명의 목적은, TOF-MRA를 이용하여 혈류나 혈류역학 등의 혈류특성을 유도하기 위한 방법을 제공함에 있다.
상기 목적을 달성하기 위한 본 발명의 일 실시예에 따른, 혈류특성 유도방법은, MRA 혈관 단면 영상을 획득하는 단계; 상기 영상에서 혈관 경계를 검출하여, 혈관을 파악하는 단계; 상기 혈관을 의사-컬러링(pseudo-coloring)하는 단계; 및 상기 혈관이 의사-컬러링된 영상을 혈류 영상으로 제공하는 단계;를 포함한다.
그리고, 상기 혈류 영상에서 열점(hotspot)의 위치는, 혈류의 치우침을 나타낼 수 있다.
또한, 본 발명의 일 실시예에 따른 혈류특성 유도방법은, 상기 혈관이 의사-컬러링된 영상을 이용하여, 혈관 내부의 특정 위치에서의 MR-신호강도구배(signal intensity gradient)를 산출하는 단계;를 더 포함할 수 있다.
그리고, 상기 MR-신호강도구배 산출단계는, "(혈관내부 영상값 - 혈관벽의 영상값)/(혈관벽으로부터 혈관내부까지의 거리)"를 상기 MR-신호강도구배로 산출할 수 있다.
또한, 본 발명의 일 실시예에 따른 혈류특성 유도방법은, 상기 혈관 내부에서의 MR-신호강도구배 분포를 산출하는 단계; 및 상기 MR-신호강도구배 분포를 상기 영상의 혈관에 도포하여, MR-신호강도구배 영상으로 제공하는 단계;를 더 포함할 수 있다.
그리고, 상기 MR-신호강도구배 영상은, 상기 혈관에 미치는 전단 응력과 장기 또는 병변과의 연관성 분석에 이용될 수 있다.
한편, 본 발명의 다른 실시예에 따른, 컴퓨터로 읽을 수 있는 기록매체는, MRA(Magnetic Resonance Angiography) 혈관 단면 영상을 획득하는 단계; 상기 영상에서 혈관 경계를 검출하여, 혈관을 파악하는 단계; 상기 혈관을 의사-컬러링(pseudo-coloring)하는 단계; 및 상기 혈관이 의사-컬러링된 영상을 혈류 영상으로 제공하는 단계;를 포함하는 혈류특성 유도방법을 수행할 수 있는 프로그램이 기록된다.
그리고, 상기 혈류 영상에서 열점(hotspot)의 위치는, 혈류의 치우침을 나타낼 수 있다.
또한, 상기 혈관이 의사-컬러링된 영상을 이용하여, 혈관 내부의 특정 위치에서의 MR-신호강도구배(signal intensity gradient)를 산출하는 단계;를 더 포함할 수 있다.
그리고, 상기 MR-신호강도구배 산출단계는, "(혈관내부 영상값 - 혈관벽의 영상값)/(혈관벽으로부터 혈관내부까지의 거리)"를 상기 MR-신호강도구배로 산출할 수 있다.
또한, 상기 혈관 내부에서의 MR-신호강도구배 분포를 산출하는 단계; 및 상기 MR-신호강도구배 분포를 상기 영상의 혈관에 도포하여, MR-신호강도구배 영상으로 제공하는 단계;를 더 포함할 수 있다.
이상 설명한 바와 같이, 본 발명에 따르면, TOF-MRA를 이용한 비교적 간단한 알고리즘을 통해 혈류나 혈류역학 등의 혈류특성을 정확하고 매우 빠르게 유도하여, 이를 분석하여 혈류특성을 파악할 수 있도록 하고, 이를 혈관질환 진단 및 치료에 응용할 수 있게 된다.
도 1은 혈관 단면에서 혈류에 의한 증강효과를 나타낸 도면,
도 2는 본 발명의 바람직한 실시예에 따른 TOF-MRA를 이용한 혈류특성 유도방법의 설명에 제공되는 흐름도,
도 3은 TOF-MRA 단면 영상을 의사-컬러링 처리한 예를 나타낸 사진,
도 4 내지 도 6은 의사-컬러링을 통해 나타나는 hotspot과 주변의 laminar flow 특성이 실제 혈류 흐름과 충분히 일치하는지 확인하기 위해, CFD 기법을 이용한 결과와 비교한 이미지들,
도 7은 기저동맥 절단면에서의 MR-신호강도구배와 CFD로 얻은 전단율을 비교한 이미지,
도 8은 기저동맥 측면에서 본 MR-신호강도구배와 CFD로 얻은 전단율을 비교한 이미지, 그리고,
도 9는 본 발명의 바람직한 실시예에 따른 혈류특성 유도방법이 수행될 수 있는 컴퓨팅 장치의 블럭도이다.
이하에서는 도면을 참조하여 본 발명을 보다 상세하게 설명한다.
중심의 빠른 혈류는 다소 신호가 약화되어 정확한 혈류속도를 얻기 힘들어 TOF-MRA를 통해 적절한 혈류특성을 끌어낼 수 없다고 알려져 왔다. 하지만, 혈관 단면 영상에 숨어있는 혈류 신호는 이를 적절히 재구성하면 혈관 특성 및 혈류 역학의 많은 정보를 뽑아낼 수 있다.
laminar flow 중심부의 신호 강도가 다소 약해지더라도 혈관벽 주변 영상 신호 강도들이 갖는 laminar flow 특성 패턴을 최대한 활용하면 가능하다. 특히, 혈관벽에 작용하는 전단 응력 및 전단율은 혈관벽 주변의 속도가 중요한 인자로, 혈관벽 주변속도를 최적화시킨 TOF-MRA 기법을 통해 혈관벽에 작용하는 기계적 힘을 분석하는 데 활용될 수 있는 바, 이하에서 상세히 설명한다.
도 2는 본 발명의 바람직한 실시예에 따른 TOF-MRA를 이용한 혈류특성 유도방법의 설명에 제공되는 흐름도이다. 본 실시예에 따른 혈류특성 유도방법에서는, 진단 및 치료에 광범위하게 사용되는 TOF-MRA(Time-Of-Flight - Magnetic Resonance Angiography) 기법을 이용하여, 혈류, 전단율, 전단응력 등의 혈류특성을 유도하여, 이를 분석하여 혈류특성을 파악할 수 있도록 하고, 이를 혈관질환 진단 및 치료에 응용할 수 있도록 한다.
도 2에 도시된 바와 같이, 먼저 MRA 혈관 단면 영상을 획득하고(S110), 획득된 MRA 혈관 단면 영상에서 혈관 경계를 검출하여, 혈관을 파악한다(S120). 다음, 파악된 혈관 부분을 의사-컬러링(pseudo-coloring)하여(S130), 혈류 영상으로 제공한다(S140).
S130단계에 의해 혈관 부분이 의사-컬러링된 영상은 혈류가 나타난 혈류 영상으로 제공되는데, 혈류 영상에서 hotspot(열점)의 위치는 혈류의 치우침을 나타낸다. 이하에서, 상세히 설명한다.
도 3은 TOF-MRA 단면 영상을 의사-컬러링 처리한 이미지이다. 도 3의 좌측은 TOF-MRA 단면영상으로 어떠한 특성을 파악하기는 힘들다. 하지만, 이를 도 3의 우측과 같이 의사-컬러링 처리하면, 단면 영상에 혈류가 나타남을 확인할 수 있다.
도 3의 단면 영상에서 혈류는 앞쪽, 우측으로(MRA기준) 치우치며 흐르고 있음을 확인할 수 있다. 또한, 혈관벽에서 보면 우측 및 앞쪽 혈관벽에 비해 좌측 및 뒤쪽 혈관벽의 혈류 속도 구배(velocity gradient)가 작은 것을 알 수 있다.
또한, 도 3에서 짙은 빨간색의 hotspot에서 혈류가 가장 빠르고, 주변으로 갈수록 느린 laminar flow 현상을 보인다.
이와 같이, 의사-컬러링을 통해 나타나는 hotspot과 주변의 laminar flow 특성이 실제 혈류 흐름과 충분히 일치하는지를 컴퓨터 시뮬레이션인 CFD(Computational Fluid Dynamics) 기법을 이용해 그 결과와 비교하였다(도 4 내지 도 6).
도 4 내지 도 6에서, 모든 CFD는 환자의 혈관을 3차원으로 재구성하고, 경계부에서의 조건은 환자로부터 직접 얻은 TCD(TransCranial Doppler) 값을 활용하였다.
도 4는 TOF-MRA 의사-컬러링 영상(Basilar artery)과 CFD를 이용해 얻은 단면의 혈류특성으로, TOF-MRA 의사-컬러링 중심부 영상신호가 약해지며(단면번호 25, 35, 40, 45), 혈류가 앞으로 치우쳐 있는 특성은 혈관 앞쪽에서 관찰되는 강한 신호영상(hotspot)과 일치한다.
도 4는 기저혈관(basilar artery) 혈류속도가 빠른 경우로 우측 추골동맥 평균혈류속도는 25.4 cm/sec, 좌측은 27.3 cm/sec이었다. 혈류 속도가 빨라 slice 25, 35, 40, 45 등에서 가운데 빠른 혈류 신호가 감쇄되는 현상이 관찰된다. CFD에서는 중심 혈류가 잘 살아 있다. 그러나, 전반적으로는 pseudo-coloring hotspot에서 보이듯 혈류가 앞으로 살짝 치우쳐 흐르고 있음을 알 수 있고 이는 CFD결과와 일치함을 알 수 있다.
도 5는 혈관협착이 있는 BA 부위에서의 TOF-MRA color tablation과 CFD를 비교한 것으로, 의사-컬러링 hotspot과 주변 혈류특성이 협착부위 (slice 40, 45)에서도 잘 나타나고 있고, 이는 CFD 결과와 일치함을 알 수 있다.
도 5는 기저혈관 중앙부위에 심한 협착이 있으며, 혈류속도는 도 4의 예보다는 다소 낮다(우측 추골동맥 평균혈류속도; 21.2 cm/sec, 좌측, 25,8 cm/sec). 협착이 심한 slice 35, 40 에서도 laminar 혈류 특성이 잘 살아있고, 전반적으로 혈류흐름 양상이 CFD 결과와 거의 일치한다.
도 6은 기저혈관(basilar artery)이 우측으로 치우친 경우의 TOF-MRA 의사-컬러링 영상과 CFD를 비교한 것으로, TOF-MRA를 의사-컬러링한 경우, 중심혈류의 신호강도가 약화되는 현상이 관찰된다(slice 50, 55, 60). 의사-컬러링 hotspot은 우측으로 치우쳐져 있고, 이는 CFD의 우측으로 치우친 laminar flow 특성과 대체로 일치한다.
도 6의 경우도 그림 4와 비슷한 혈류속도를 보인다(우측 추골동맥 평균혈류속도; 25.8 cm/sec, 좌측 19.6 cm/sec). 중심의 빠른 혈류는 신호강도가 약해지는 현상이 slice 50, 55, 60 등에서 관찰된다. 그러나, 우측편향이 심한 기저동맥으로 혈류가 우측으로 치우친 현상이 TOF-MRA 의사-컬러링 영상과 CFD 모두에서 공통으로 관찰된다.
도 4 내지 도 6에서 살펴보았듯 결과는 일관 되게 나타난다. 앞서 설명한 바와 같이 중심부 혈류 및 laminar 혈류특성은 다소 상쇄될 수 있으나, 혈관구조에 따른 혈류 치우침 현상과 주변부 혈류특성은 CFD와 거의 일치하는 양상임을 알 수 있다.
모든 혈관벽 주변부의 신호 강도는 파란색의 가장 느린 혈류특성을 보이며, 그 주변부는 혈관벽보다는 다양한, 약간 더 빠른 혈류속도를 보이고 있다. 특히, 혈류속도가 빨라 중심부 혈류신호가 많이 상쇄되는 경우라 하더라도(예, 도 4의 slice 35, 40; 도 6의 slice 55, 60), 상쇄된 부위의 속도분포는 비록 원래의 혈류속도 특성을 보이지는 않으나, 혈관벽 및 벽주변부 혈류속도와는 현저한 차이가 있음을 알 수 있다.
도 4 내지 도 6에서 다양한 임상례에서 MRA 의사-컬러링을 통해 나타나는 hotspot과 laminar 혈류 특성 및 혈관벽 주변부 혈류특성을 혈관 CFD의 결과와 비교/종합하면, 1) 혈류 중심 영상신호 강도가 실제보다 약해지는 현상이 나타나고, 특히 혈류속도가 빠른 경우 더욱 심해지며, 2) hotspot은 혈류의 치우침 현상을 잘 보여주고, 3) 혈관벽 근처에서의 laminar 혈류특성은 CFD 결과와도 잘 일치되며, 위 1), 2)의 특성에 영향을 받지 않음을 도출할 수 있다.
다시, 도 2를 참조하여 설명한다.
S130단계에서 혈관이 의사-컬러링된 영상을 이용하여, 혈관에 대한 MR-신호강도구배(MR-signal intensity gradient)[또는 MR-전단율(MR-shear rate)] 영상을 생성하여 제공한다(S150). S150단계에 제공되는 MR-신호강도구배 영상은, 혈관이 의사-컬러링된 혈류 영상에서 혈관 내부에서의 MR-신호강도구배 분포를 산출하고, MR-신호강도구배 분포를 혈관에 도포함으로 생성 가능하다.
MR-신호강도구배 영상은, 혈관에 미치는 힘과 뇌(혹은 장기)의 구조나 병변과의 연관성의 정성/정략 분석에 이용될 수 있는 바, 이하에서 상세히 설명한다.
MRA 혈관벽 주변 신호는 도 4 내지 도 6에서 확인하였듯 MRA 특성에 의한 중심 혈류신호 감쇄 효과에 전혀 영향을 받지 않았으므로 이를 이용해 MR-신호강도구배를 구할 수 있다.
구체적으로, 혈관 내부의 특정 위치에서의 MR-신호강도구배는, 아래의 수학식 1을 통해 산출가능하다.
[수학식 1]
MR-signal intensity gradient = (MR signal B - MR signal A) / D
여기서, "MR-signal intensity gradient"는 혈관 내부의 특정 위치에서의 MR-신호강도구배이고, "MR signal B"은 혈관내부 특정 위치에서의 영상값이며, "MR signal A"은 혈관벽의 영상값이고, "D"는 혈관벽으로부터 혈관내부 특정 위치까지의 거리이다.
의사-컬러링 영상과 혈관벽 근처 영상 데이타를 이용하여 구한 MR-신호강도구배를 CFD 결과와 비교한 결과, 절단면 영상을 통해 분석한 결과 MR-신호강도구배와 CFD 결과는 매우 유사하였다. 절단면은 각 단면 slice 번호를 통해 정확한 레벨을 얻어 비교하였다.
도 7은 기저동맥 절단면에서의 MR-신호강도구배(A)와 CFD로 얻은 전단율(B)을 비교한 이미지이다. 두 이미지 모두에서 좌측 전방으로 높은 전단율 부위가 관찰되고, 후방으로 낮은 전단율 부위가 관찰된다.
도 7의 경우 기저혈관의 형상에 따라 혈류속도가 좌측 전방으로 치우쳐 흐름을 알 수 있었고, MR-신호강도구배도 일치해서 나타났다. CFD 결과는 가장 큰 전단율을 MR-신호강도구배 보다 약간 더 전방으로 치우쳐 있다. 이는 혈관벽을 3차원 재구성할 때 사용하는 혈관벽 threshold가 단면 영상에서 사용하는 것보다는 다소 높아 발생하는 오류로 판단된다. 이러한 오류에도 불구하고 두 전단율의 결과는 매우 유사함을 알 수 있다.
도 8은 기저동맥 측면에서 본 MR-신호강도구배(A)와 CFD로 얻은 전단율(B)을 비교한 이미지이다. 아래측 추골 동맥에서는 전반적으로 전단율이 낮은 반면, 기저동맥 협착부에서는 높은 전단율을 보이며, 이는 A, B 모두 유사함을 확인할 수 있다.
도 8의 경우도 아래 추골동맥(vertebral artery)에서는 매우 낮은 전단율을 보이고, 기저혈관 중앙부에서 매우 높은 전단율이 관찰되는 데, MR-신호강도구배와 CFD 결과가 거의 유사한 결과를 보인다.
도 7 및 도 8에서, 왼쪽에 보인 MR-신호강도구배는 총 산출 소요시간이 불과 수 초에 불과한 반면 CFD를 통해 얻은 결과는 수 시간이 소요되어, 유사한 결과임에도 불구하고 그 결과를 얻기까지 소요된 시간 차이는 대단히 크다.
매우 짧은 시간에 수행된 MR-신호강도구배이지만 이는 TOF-MRA 영상에 들어온 모든 혈관을 동시에 처리할 수 있어 실질적으로 모든 혈관에 대한 CFD 정보를 구하고자 한다면 두 방법의 시간차는 각각의 혈관 수를 곱해서 계산해야 한다.
도 9는 본 발명의 바람직한 실시예에 따른 혈류특성 유도방법이 수행될 수 있는 컴퓨팅 장치의 블럭도이다. 도 9에 도시된 바와 같이, 본 실시예에 따른 혈류특성 유도방법을 수행할 컴퓨팅 장치(200)는, 통신 인터페이스(210), 모니터(220), 프로세서(230), 사용자 인터페이스(240) 및 저장매체(250)를 구비한다.
통신 인터페이스(210)는 외부기기와 통신 연결을 설정하고 유지하는 수단으로, 이때 외부기기는 TOF-MRA 영상장비가 될 수 있다.
사용자 인터페이스(240)는 컴퓨팅 장치(200)를 조작하기 위한 명령을 입력하는 수단으로 키보드, 마우스 등을 포함한다.
저장매체(250)는 도 2에 도시된 혈류특성 유도방법을 수행할 수 있는 프로그램과 이 프로그램을 수행하는데 필요한 데이터가 저장되는 매체이다. 또한, 저장매체(250)에는 TOF-MRA 영상이 저장되어 있을 수도 있다.
모니터(220)는 혈류특성 유도방법을 수행함에 있어 나타나는 과정과 결과가 표시되는 디스플레이이다. 프로세서(230)는 도 2에 도시된 혈류특성 유도방법을 알고리즘으로 수행하기 한다.
지금까지, 혈류특성 유도방법 및 이를 수행할 수 있는 컴퓨팅 장치에 대해 바람직한 실시예들을 들어 상세히 설명하였다. 한편, 위 실시예에 따른 혈류특성 유도방법을 프로그램으로 구현하는 경우도 본 발명의 기술적 범주에 포함됨은 물론이다.
또한, 이상에서는 본 발명의 바람직한 실시예에 대하여 도시하고 설명하였지만, 본 발명은 상술한 특정의 실시예에 한정되지 아니하며, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 당해 발명이 속하는 기술분야에서 통상의 지식을 가진자에 의해 다양한 변형실시가 가능한 것은 물론이고, 이러한 변형실시들은 본 발명의 기술적 사상이나 전망으로부터 개별적으로 이해되어져서는 안될 것이다.
200 : 컴퓨팅 장치 210 : 통신 인터페이스
220 : 모니터 230 : 프로세서
240 : 사용자 인터페이스 250 : 저장매체

Claims (11)

  1. MRA(Magnetic Resonance Angiography) 혈관 단면 영상을 획득하는 단계;
    상기 영상에서 혈관 경계를 검출하여, 혈관을 파악하는 단계;
    상기 혈관을 의사-컬러링(pseudo-coloring)하는 단계;
    상기 혈관이 의사-컬러링된 영상을 혈류 영상으로 제공하는 단계; 및
    상기 혈관이 의사-컬러링된 영상을 이용하여, 혈관 내부의 특정 위치에서의 "(혈관내부 영상값 - 혈관벽의 영상값)/(혈관벽으로부터 혈관내부까지의 거리)"를 MR-신호강도구배(signal intensity gradient)로 산출하는 단계;를 포함하는 것을 특징으로 하는 혈류특성 유도방법.
  2. 제 1항에 있어서,
    상기 혈류 영상에서 열점(hotspot)의 위치는,
    혈류의 치우침을 나타내는 것을 특징으로 하는 혈류특성 유도방법.
  3. 삭제
  4. 삭제
  5. 제 1항에 있어서,
    상기 혈관 내부에서의 MR-신호강도구배 분포를 산출하는 단계; 및
    상기 MR-신호강도구배 분포를 상기 영상의 혈관에 도포하여, MR-신호강도구배 영상으로 제공하는 단계;를 더 포함하는 것을 특징으로 하는 혈류특성 유도방법.
  6. 제 5항에 있어서,
    상기 MR-신호강도구배 영상은,
    상기 혈관에 미치는 전단 응력과 장기 또는 병변과의 연관성 분석에 이용되는 것을 특징으로 하는 혈류특성 유도방법.
  7. MRA(Magnetic Resonance Angiography) 혈관 단면 영상을 획득하는 단계;
    상기 영상에서 혈관 경계를 검출하여, 혈관을 파악하는 단계;
    상기 혈관을 의사-컬러링(pseudo-coloring)하는 단계;
    상기 혈관이 의사-컬러링된 영상을 혈류 영상으로 제공하는 단계; 및
    상기 혈관이 의사-컬러링된 영상을 이용하여, 혈관 내부의 특정 위치에서의 "(혈관내부 영상값 - 혈관벽의 영상값)/(혈관벽으로부터 혈관내부까지의 거리)"를 MR-신호강도구배(signal intensity gradient)로 산출하는 단계;를 포함하는 것을 특징으로 하는 혈류특성 유도방법을 수행할 수 있는 프로그램이 기록된 컴퓨터로 읽을 수 있는 기록매체.
  8. 제 7항에 있어서,
    상기 혈류 영상에서 열점(hotspot)의 위치는,
    혈류의 치우침을 나타내는 것을 특징으로 하는 혈류특성 유도방법을 수행할 수 있는 프로그램이 기록된 컴퓨터로 읽을 수 있는 기록매체.
  9. 삭제
  10. 삭제
  11. 제 7항에 있어서,
    상기 혈관 내부에서의 MR-신호강도구배 분포를 산출하는 단계; 및
    상기 MR-신호강도구배 분포를 상기 영상의 혈관에 도포하여, MR-신호강도구배 영상으로 제공하는 단계;를 더 포함하는 것을 특징으로 하는 혈류특성 유도방법을 수행할 수 있는 프로그램이 기록된 컴퓨터로 읽을 수 있는 기록매체.
KR1020120081035A 2012-07-25 2012-07-25 Ｔｏｆ-ｍｒａ를 이용한 혈류특성 및 ｍｒ-신호강도구배(전단율) 유도방법 KR101373563B1 (ko)

Priority Applications (2)

Application Number Priority Date Filing Date Title
KR1020120081035A KR101373563B1 (ko) 2012-07-25 2012-07-25 Ｔｏｆ-ｍｒａ를 이용한 혈류특성 및 ｍｒ-신호강도구배(전단율) 유도방법
PCT/KR2013/006684 WO2014017853A2 (ko) 2012-07-25 2013-07-25 Ｔｏｆ-ｍｒａ를 이용한 혈류특성 및 ｍｒ-신호강도구배(전단율) 유도방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020120081035A KR101373563B1 (ko) 2012-07-25 2012-07-25 Ｔｏｆ-ｍｒａ를 이용한 혈류특성 및 ｍｒ-신호강도구배(전단율) 유도방법

Publications (2)

Publication Number Publication Date
KR20140014605A KR20140014605A (ko) 2014-02-06
KR101373563B1 true KR101373563B1 (ko) 2014-03-12

Family

ID=49997930

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020120081035A KR101373563B1 (ko) 2012-07-25 2012-07-25 Ｔｏｆ-ｍｒａ를 이용한 혈류특성 및 ｍｒ-신호강도구배(전단율) 유도방법

Country Status (2)

Country Link
KR (1) KR101373563B1 (ko)
WO (1) WO2014017853A2 (ko)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101992057B1 (ko) 2018-08-17 2019-06-24 (주)제이엘케이인스펙션 혈관 투영 영상을 이용한 뇌질환 진단 방법 및 시스템
KR20220063031A (ko) 2020-11-09 2022-05-17 울산대학교 산학협력단 3차원 영상 기반 유체 시뮬레이션 방법,　장치 및 프로그램
US11579226B2 (en) 2020-09-16 2023-02-14 Siemens Healthcare Gmbh System and method for time of flight imaging with a tight sequence diagram pattern

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
PE20181268A1 (es) 2011-06-17 2018-08-03 Hanmi Science Co Ltd Un conjugado que comprende oxintomodulina y un fragmento de inmunoglobulina, y el uso del mismo
AR093387A1 (es) 2012-11-06 2015-06-03 Hanmi Pharm Ind Co Ltd Formulacion liquida de un conjugado de proteina que comprende oxintomodulina y un fragmento de inmunoglobulina
TWI802396B (zh) 2014-09-16 2023-05-11 南韓商韓美藥品股份有限公司 長效glp-1/高血糖素受體雙促效劑治療非酒精性脂肝疾病之用途
KR102418477B1 (ko) 2014-12-30 2022-07-08 한미약품 주식회사 글루카곤 유도체
KR102462075B1 (ko) * 2020-08-05 2022-11-01 정슬기 뇌자기공명혈관조영술의 혈관 신호강도그래디언트를 이용한 혈류산출법 및 이를 이용한 뇌혈관질환 및 뇌졸중 위험도 분석 시스템

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08154917A (ja) * 1994-12-06 1996-06-18 Hitachi Medical Corp 血流情報表示方法
JPH10179545A (ja) * 1996-12-26 1998-07-07 Hitachi Medical Corp 磁気共鳴イメージング装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4139869B2 (ja) * 2005-08-31 2008-08-27 国立大学法人岐阜大学 医用画像処理装置
JP5002214B2 (ja) * 2005-09-22 2012-08-15 株式会社東芝 磁気共鳴イメージング装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08154917A (ja) * 1994-12-06 1996-06-18 Hitachi Medical Corp 血流情報表示方法
JPH10179545A (ja) * 1996-12-26 1998-07-07 Hitachi Medical Corp 磁気共鳴イメージング装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101992057B1 (ko) 2018-08-17 2019-06-24 (주)제이엘케이인스펙션 혈관 투영 영상을 이용한 뇌질환 진단 방법 및 시스템
US11579226B2 (en) 2020-09-16 2023-02-14 Siemens Healthcare Gmbh System and method for time of flight imaging with a tight sequence diagram pattern
KR20220063031A (ko) 2020-11-09 2022-05-17 울산대학교 산학협력단 3차원 영상 기반 유체 시뮬레이션 방법,　장치 및 프로그램

Also Published As

Publication number Publication date
WO2014017853A2 (ko) 2014-01-30
KR20140014605A (ko) 2014-02-06
WO2014017853A3 (ko) 2014-03-13

Similar Documents

Publication Publication Date Title
KR101373563B1 (ko) Ｔｏｆ-ｍｒａ를 이용한 혈류특성 및 ｍｒ-신호강도구배(전단율) 유도방법
Grandin et al. Usefulness of magnetic resonance–derived quantitative measurements of cerebral blood flow and volume in prediction of infarct growth in hyperacute stroke
US10677866B1 (en) Systems and methods for correcting measurement artifacts in MR thermometry
Srinivasan et al. State-of-the-art imaging of acute stroke
JP6679467B2 (ja) 磁気共鳴イメージング装置および酸素摂取率算出方法
Taxt et al. Single-channel blind estimation of arterial input function and tissue impulse response in DCE-MRI
JP2007525250A (ja) 組織血流の遅延−補償型算出方法
JP6240319B2 (ja) 側副循環の画像化及び評価のための新しい磁気共鳴画像法
Lin et al. Characterization of cardiac dynamics in an acute myocardial infarction model by four-dimensional optoacoustic and magnetic resonance imaging
WO2022198719A1 (zh) 血流动力的仿真方法及装置
Meves et al. Comparison between echo contrast agent-specific imaging modes and perfusion-weighted magnetic resonance imaging for the assessment of brain perfusion
Pereira et al. 4D flow MRI in neuroradiology: techniques and applications
Bolcaen et al. In vivo DCE-MRI for the discrimination between glioblastoma and radiation necrosis in rats
Celicanin et al. Hybrid ultrasound‐MR guided HIFU treatment method with 3 D motion compensation
US10213178B2 (en) Calculation of perfusion parameters in medical imaging
Zhang et al. Use of cardiac output to improve measurement of input function in quantitative dynamic contrast‐enhanced MRI
WO2016094461A1 (en) Methods, systems, and computer readable media for estimating the index of cerebral oxygen metabolism (i-com) using magnetic resonance (mr) imaging
Lawrence et al. Longitudinal characterization of local perfusion of the rat placenta using contrast-enhanced ultrasound imaging
Nomura et al. Registration method between phase-contrast magnetic resonance angiography and time-of-flight magnetic resonance angiography—a preliminary study
Pan et al. Color Doppler ultrasonography in the evaluation of compensatory arteries in patients with moyamoya disease: combined with cerebral angiography.
KR101507762B1 (ko) Tof-mra의 신호강도를 이용한 혈류속도의 유도 시스템 및 방법
JP2016131836A (ja) 血流プロファイルの測定方法
Hartono et al. High temporal resolution dynamic contrast-enhanced MRI at 7 Tesla: a feasibility study with mouse liver model
Mamatyukov et al. Comprehensive research of human brain hemodynamics: Clinical monitoring and computer simulations
Smith et al. The effects of geometry on stiffness measurements in high-field magnetic resonance elastography: A study on rodent cardiac phantoms

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
FPAY Annual fee payment

Payment date: 20170221

Year of fee payment: 4

FPAY Annual fee payment

Payment date: 20180208

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20190422

Year of fee payment: 6

FPAY Annual fee payment

Payment date: 20200302

Year of fee payment: 7