WO2014196732A1 - 측부순환의 영상화 및 평가를 위한 새로운 자기공명영상기법 - Google Patents

측부순환의 영상화 및 평가를 위한 새로운 자기공명영상기법 Download PDF

Info

Publication number
WO2014196732A1
WO2014196732A1 PCT/KR2014/003346 KR2014003346W WO2014196732A1 WO 2014196732 A1 WO2014196732 A1 WO 2014196732A1 KR 2014003346 W KR2014003346 W KR 2014003346W WO 2014196732 A1 WO2014196732 A1 WO 2014196732A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
circulation
dsc
pwi
phase
Prior art date
Application number
PCT/KR2014/003346
Other languages
English (en)
French (fr)
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 US14/895,597 priority Critical patent/US20160125597A1/en
Priority to EP14807402.4A priority patent/EP3005943B1/en
Priority to JP2016518256A priority patent/JP6240319B2/ja
Publication of WO2014196732A1 publication Critical patent/WO2014196732A1/ko

Links

Classifications

    • 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
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0263—Measuring blood flow using NMR
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • 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/5601—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
    • 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/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/56366—Perfusion imaging
    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00—Medical imaging apparatus involving image processing or analysis
    • A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
    • A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
    • 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]
    • G06T2207/10096—Dynamic contrast-enhanced magnetic resonance imaging [DCE-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/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
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20212—Image combination
    • G06T2207/20224—Image subtraction
    • 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/30016—Brain
    • 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/30048—Heart; Cardiac
    • 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
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing

Definitions

  • the present invention relates to a novel magnetic resonance imaging technique for imaging and evaluation of the side circulation.
  • Cerebral infarction is one of the leading causes of death worldwide. The single cause of death among Koreans is the highest mortality rate and high disability rate. Cerebral infarction is a symptom that is caused by the blockage of the blood supply to the brain tissue as the blood vessels are blocked and the infarction of the brain tissue that received blood. Treatment for cerebral infarction is almost the only treatment by reopening blocked blood vessels early and restoring cerebral blood flow before brain tissue is completely infarcted. Cerebral blood flow is maintained by lateral circulation from the surrounding blood vessels until the blocked cerebrovascular vessels reopen. Thus, the degree of lateral circulation prior to reopening treatment plays a critical role in the fate of brain tissue infarction [Bang OY et al.
  • Representative techniques for cerebral infarction using magnetic resonance imaging include diffusion and per fus i on-we ght ed imaging (FLAIR) and Fluid Attenuation Inversion Recovery (FLAIR). These images and maps created through post-processing of images are actively used as image biomarkers in clinical practice. However, despite the development of these various techniques, The method of directly imaging the lateral circulation using magnetic resonance imaging and evaluating the degree has not been developed yet.
  • DSA digital substraction angiography
  • the DSA technique is invasive to the human body, takes a relatively long time to implement, requires an expert, and carries the risk of developing thromboembolism.
  • the side circulation evaluation may be limited in cases where the delayed venous phase is not sufficiently performed or the evaluation of blood vessels opposite the lesion is missed. Therefore, there is a need for the development of a new technique that is non-invasive to the human body and easy and quick to know the degree of side circulation.
  • the present inventors have made a thorough research to solve the above problems of the prior art, and have completed the present invention by developing a new method for imaging and evaluating the side circulation using magnetic resonance imaging techniques. Accordingly, it is an object of the present invention to provide a method of imaging the lateral circulation in the body.
  • Another object of the present invention is to provide a computer readable recording medium having recorded thereon a program for imaging the side circulation.
  • the present invention provides a method for imaging lateral circulation in the body comprising the following steps:
  • step (b) imaging the side circulation by comparing the image obtained through post-processing in step (a) with an image of digital subtraction angiography (DSA).
  • DSA digital subtraction angiography
  • hypothalastomosis refers to a phenomenon in which blood flows bypass to other blood vessels when normal blood flow is blocked due to narrowing or occlusion of blood vessels.
  • cerebral blood flow In the case of cerebral blood flow, it consists of anterior circulat ion by bilateral carotid artery and posterior circulat ion by posterior vertebral basal artery. After that, it is branched into anterior cerebral artery, middle cerebral artery and posterior cerebral artery. Each of these vessels is ultimately responsible for blood supply to certain areas of the brain, but the anastomosis between these vessels occurs in several places, resulting in lateral circulation in the lesion of one vessel.
  • the side circulation in the present invention includes the side circulation in the brain, heart, extremities, kidneys, and intestine, and more preferably in the present invention, the side circulation is cerebral infarction or Lateral circulation that occurs during myocardial infarction.
  • Cerebral infarction refers to a disease in which part of the brain dies because blood vessels of the brain are blocked. Cerebral infarction is caused by blood clots caused by arteriosclerosis of cerebrovascular vessels and occlusion of cerebrovascular vessels by embolism derived from the heart.
  • magnetic resonance imaging refers to a diagnostic technique that uses a magnetic field to create a picture of an internal structure.
  • the machine used for magnetic resonance imaging forms a very powerful magnetic field in the body by using the pulse of electromagnetic waves, and the structure of the body is imaged using the signals generated in the process of formation and collapse of the magnetic field.
  • dynamic susceptibility contrast enhancement is used herein.
  • Dynamic susceptibility contrast-enhanced perfusion weighted imaging DSC-PWI refers to a technique for magnetic resonance imaging that detects changes in the magnetic field caused by an increase in susceptibility to magnetization.
  • DTPA gadolinium diethylene triamine pent a acetic acid
  • image post-processing refers to a series of operations collectively outputting a new image after correction and reconstruction using information of the image as source data after the acquisition of the image.
  • a medical image may be defined by a three-dimensional function f (x, y, z), where X, y, and z represent spatial coordinates, and the magnitude of f at the coordinates is determined by the signal intensity of the image. It is called.
  • medical images are digital images, they consist of finite elements, each of which has a specific position and value, which is called a voxel. Signal strength values in each voxel are corrected and reconstructed through the image post-processing process, and output as an image composed of new information.
  • DSA digital substraction angiography
  • the image post-processing of the DSC-PWI source data of step (a) is sequentially performed over time by subtracting the signal of the image obtained before the contrast agent injection from the signals of the images obtained after the contrast agent injection.
  • Step (a) -l obtaining an image.
  • the tissue to be imaged is excluded in order to exclude the signal due to the anatomical structure and leave only the signal by the contrast agent. Can be visualized.
  • the sequential images over time of the DSC-PWI mean a sequential image having a predetermined time interval, the predetermined time interval is preferably 0.1 ⁇ 10 seconds, more preferably 0.5 ⁇ 5 seconds, Even more preferably 1.0-3 seconds, most preferably 1.5-2 seconds.
  • the image post-processing of the DSC-PWI source data of step (a) overlaps adjacent continuous image cuts in the sequential images obtained in step (a) -l.
  • Step (a) -2 by superimposing adjacent continuous image cuts in the sequential images of the DSC-PWI, the spatial resolution of the image is improved and the signal intensity change over time is maximized.
  • the number of overlapping continuous image cuts may be appropriately selected within a range capable of maximizing a change in signal strength over time while improving spatial resolution.
  • the superimposition of adjacent continuous image cuts may be performed by overlapping two to ten adjacent continuous image cuts, preferably 2 to 7 pieces, more preferably 2 to 5 pieces, and even more preferably. Superimposes two to three, most preferably two, adjacent consecutive image cuts.
  • the contrast between the post-processed DSC-PWI image and DSA image in the arterial phase (arterial phase, capillary phase, venous phase) , Or images corresponding to the late venous phase.
  • the term “arterial phase” refers to when contrast staining starts from the insular region of the normal brain parenchyma
  • the term “capillary vessel” refers to the time when the whole normal brain parenchyma is stained with the highest signal intensity.
  • Vein refers to the time when contrast medium is washed out in the normal brain parenchyma
  • late venous refers to the time when the contrast brain signal remains in the lesion brain parenchyma, while the contrast brain is lost all contrast medium.
  • the DSC-PWI image in the arterial, capillary, venous, and late venous phases is a continuous image of the above-described "overlapping continuous image cut image". It is composed. More preferably, in the arterial, venous, and late venous phases, the DSC-PWI image consists of two consecutive images of the above-described "overlapping continuous image cuts", wherein the capillary The DSC-PWI image consists of three consecutive images of the above-described "image in which adjacent continuous image cuts are overlapped".
  • the DSC-PWI image constituting the arterial group is overlapped to form an early phase image of the side circulation
  • the DSC-PWI image constituting the capillary group is lateral circulation by overlapping A mid phase image of the late vein and late vein of the DSC-PWI image constituting the overlap (late phase) image (late phase) of the late phase (late phase).
  • the present invention provides a method for producing a side circulation image comprising the following steps:
  • the invention provides a computer readable recording medium having recorded thereon a program for producing a side circulation image, comprising the following steps:
  • the present invention relates to a novel magnetic resonance imaging technique for rapidly and accurately imaging the lateral circulation in the body.
  • the present invention relates to dynamic susceptibility contrast enhancement perfusion-enhanced imaging (DSC-PWI) After post-processing the source data, the post-processing image is contrasted with digital subtraction angiography (DSA) to image the lateral circulation.
  • DSC-PWI dynamic susceptibility contrast enhancement perfusion-enhanced imaging
  • non-invasive, rapid and easy imaging and evaluation of collateral circulation can provide important information in understanding diseases caused by vascular narrowing or occlusion, determining treatment policy, and prognostic prediction. have.
  • the present invention relates to a novel magnetic resonance imaging technique for imaging the lateral circulation in the body.
  • the present invention post-processes source data of dynamic susceptibility contrast-enhanced perfusion-enhanced perfusion-enhanced imaging (DSC-PWI), and then compares the image obtained through post-processing with digital subtraction angiography (DSA) to image the lateral circulation.
  • DSC-PWI dynamic susceptibility contrast-enhanced perfusion-enhanced perfusion-enhanced imaging
  • DSA digital subtraction angiography
  • non-invasive rapid and easy imaging and evaluation of the lateral circulation can provide important information for understanding the disease caused by vascular narrowing or occlusion, determining treatment policy, and predicting prognosis.
  • DSC-PWI dynamic susceptibility contrast-enhanced perfusion weighted imaging
  • FIG. 2 shows the result of superimposing DSC-PWI image cuts in order to obtain the best temporal / spatial resolution during side cycle imaging.
  • the best temporal / spatial resolution can be obtained by superimposing two images of two adjacent images. Therefore In contrast, non-overlapping images have a lower spatial resolution, and when more overlaps, the temporal resolution drops.
  • the first line is an image obtained by digital subtraction angiography performed on the opposite side of the lesion and consists of an arterial phase, a capi 1 lary phase, and a venous phase, respectively.
  • the second row is an image obtained by digital subtraction angiography performed on the lesion side.
  • the occlusion of the proximal middle cerebral artery can be confirmed (thin white arrow).
  • lateral circulation from the capillary vessel to the right forearm cerebral artery can be observed (thick white arrows). Lateral circulation can be observed up to the late venous phase.
  • the third row shows a two-cut superimposed image of the DSC—PWI image corresponding to each period of the digital subtracted angiography image. Staining and washout of contrast media by lateral circulation in the lesion side hemisphere can be identified (black arrow heads).
  • the fourth row shows the generation of side circulation images. Two-cut overlapping images of DSC-PWIs corresponding to arterial, capillary, and venous-late venous phases were once again superimposed to generate images capable of evaluating lateral circulation. Each phase corresponds to the early phase, mid phase and late phase.
  • Figure 4 shows the process of creating a side circulation image using the automation software FAST-C0LL.
  • the patient shows the overall flow from the occurrence of cerebral infarction to the hospital and the generation of lateral circulatory images.
  • the image is sent to the workstation, and the side circulation image is generated using FAST-C0LL. Based on the result, the degree of side circulation is evaluated.
  • the generated circulatory images are sent to a picture archiving system (PACS) to provide a variety of information about the patient, along with information from other multimodal MRIs.
  • PACS picture archiving system
  • Panel 5 shows the results of evaluating the degree of collateral circulation in the collateral circulation image.
  • Panel (1) is an image evaluated with poor collaterals, where the contrast agent stain on the lesion side is completely absent until the late phase or only in some areas. When stained with slow side circulation.
  • Panel 2 is an image evaluated with intermediate col laterals in which some areas of the lesion side are stained with fast lateral circulation while some of the lesions are completely absent.
  • Panel 3 is an image evaluated for good collaterals, where the contrast agent stain on the lesion side is completely stained with slow lateral circulation.
  • Panel 4 is an image evaluated with excellent collaterals, where the contrast agent stain is completely stained with fast collaterals.
  • FIG. 6 is a graph showing the change over time of the combined R 2 * values by adding the R 2 * values of all voxels in each cut of DSC-PWI source data.
  • Figures 7a to 7d is a 'graph showing the outcome after cerebral infarction with and without rating and revascularization opening of collateral circulation that is observed in the collateral image.
  • ER early recanalizat ion
  • G grade
  • Each grade Gl, G2, G3 and G4 of the lateral circulation has poor collaterals, intermediate collaterals, good collaterals, and superior lent col laterals. (See Example 5 and FIG. 5), respectively.
  • FIG. 7A the higher the degree of side circulation, the higher the probability of having a good functional prognosis.
  • FIG. 7B the higher the degree of lateral circulation, the higher the probability of showing good neurological improvement.
  • FIG. 7C the higher the grade of lateral circulation, the lower the probability of developing symptomatic cerebral hemorrhage.
  • FIG. 7D the grade of lateral circulation is higher and the specific mortality rate is lower.
  • Example 2 Image Superposition for Obtaining the Best Temporal and Spatial Resolution
  • images obtained from the DSC ⁇ PWI were gradually superimposed to compare the temporal and spatial resolutions to each other (Fig.
  • Example 2 In the process of Example 2 as a final step for imaging the side circulation The determined two-cut fold images were matched with the images obtained by digital subtraction angiography (FIG. 3). To this end, visual contrast of staining and washout of the normal and lesion brain parenchyma is visually determined in the two-cut superimposed images to determine the arterial phase and capillary phase of digital subtraction angiography. The images corresponding to the venous phase were selected. Specifically, when the arterial phase begins staining of contrast agent from the insular region of the normal brain parenchyma, the capillary vessel is stained with the highest signal intensity, and the venous phase is the normal brain parenchyma.
  • Each of these phases typically consists of three consecutive images of capillary, two consecutive arterial venous and late venous two consecutive superimposed images.
  • Lateral circulatory images were generated by superimposing two superimposed images of two arteries, three superimposed images of capillary vessels, and four superimposed images of veins and late veins.
  • Each phase of the generated side circulation image corresponds to an early phasae, a mid phase and a late ate phase.
  • the degree of side circulation can be evaluated from the generated image.
  • the degree of side circulation can be largely divided into no side circulation, slow side circulation and fast side circulation.
  • No collateral is if the show does end visible staining of the contrast agent in the (late phase) end and (late phase) when the contrast dye completely invisible to slow, collateral is middle (m id phase), Quick collateral circulation It refers to the case where the contrast agent is seen from the mid phase in the lesion side to the extent similar to the normal side. Based on this, the side circulation can be evaluated. As shown in the side circulation image of FIG.
  • the prognosis according to vascular reopening after vascular revascularization can be predicted in patients with acute cerebral infarction.
  • Figures 7a to 7d regardless of whether the vascular reopening in acute cerebral infarction, the better the lateral circulation, the higher the probability of having a good prognosis, and the lower the probability of symptomatic bleeding and death.
  • the vessels are reopened in the case of high lateral circulation observed in the lateral circulation images, a good prognosis can be expected without the risk of symptomatic cerebral hemorrhage.

Landscapes

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

Abstract

본 발명은 신체내 측부순환을 영상화하는 새로운 자기공명영상 기법에 관한 것이다. 본 발명은 역동적 자화율대조조영증강 관류강조영상(DSC-PWI)의 소스 데이터를 영상 후처리한 후 후처리를 통해 얻은 영상을 디지털 감산 혈관 조영술(DSA)의 영상과 대조하여 측부순환을 영상화한다. 본 발명의 방법에 의하면 측부순환을 비침습적으로 신속하면서도 용이하게 영상화하여 평가할 수 있으므로 혈관 협착이나 폐색에 의해 발생되는 질환의 이해, 치료 방침 결정 및 예후 예측에 중요한 정보를 제공할 수 있다.

Description

【명세서】
【발명의 명칭】
측부순환의 영상화 및 평가를 위한 새로운 자기공명영상기법 【기술 분야】
본 발명은 측부순환의 영상화 및 평가를 위한 새로운 자기공명영상기법에 관한 것이다.
【배경 기술】
뇌경색은 전 세계적으로 중요한 사망 원인 중 하나이며, 우리나라 국민의 사망원인 중 단일 질환으로는 가장 높은 사망률을 보이고 높은 장애율을 나타낸다. 뇌경색은 뇌혈관이 막히면서 뇌조직으로 가는 혈액의 공급이 차단되어 혈액을 공급받던 뇌조직이 경색되어 나타나는 증상이다. 뇌경색의 치료 방법은 막힌 혈관을 조기에 재개통하여 뇌조직이 완전히 경색에 빠지기 전에 뇌혈류를 회복시키는 것이 거의 유일한 치료 방법이다. 뇌혈류는 막힌 뇌혈관이 재개통 되기 전까지 주변 혈관에서 공급되는 측부순환에 의해 유지된다. 따라서， 재개통 치료 전 측부순환의 정도는 뇌조직이 경색에 빠질 지에 대한 운명에 결정적인 역할을 한다 [Bang OY et al. , J Neurol Neurosurg Psychiatry, 2008 ;79： 625] . 또한， 측부순환은 막힌 혈관의 재개통률 및 치료 후 뇌조직의 출혈성 변화 여부에도 지대한 영향을 미친다 [Bang OY et al. , Stroke, 2011 ;42： 693； Bang 0Y et al. , Stroke, 2011； 42： 2235] . 즉， 뇌경색 환자의 치료방침 결정 및 예후 예축에 있어서 측부순환의 정도를 조기에 아는 것은 매우 중요하다. 일반적으로 뇌경색은 발생 원인과 치료 반웅이 다양하기 때문에 다양한 자기공명영상 기법을 이용한 영상 바이오마커의 개발이 매우 중요하다. 자기공명영상올 이용한 뇌경색의 연구 기법에는 대표적으로 확산 및 관류 강조 영상 (diffusion and per f us i on-we i ght ed imaging) , 액체감약반전회복 (FLAIR; Fluid Attenuation Inversion Recovery) 기법 등이 있으며， 이러한 영상들 및 영상의 후처리 (post-processing)를 통해 만들어진 지도 (map)들은 임상에서 영상 바이오마커로 활발하게 사용되고 있다. 그러나 이러한 다양한 기법들이 개발되었음에도 불구하고， 현재까지 자기공명영상을 이용하여 직접적으로 측부순환을 영상화하고 그 정도를 평가할 수 있는 방법은 아직 개발되지 않은 실정이다.
현재 임상 진료에서는 통상적으로 디지털 감산 혈관조영술 (digital substraction angiography, DSA)이 측부순환의 영상화 및 정도 측정의 황금 기준 (gold standard)으로 사용되고 있다. 그러나 DSA 기법은 인체에 침습적이고， 시행을 위해 상대적으로 많은 시간이 소요되고 이를 위한 전문가가 필요하며， 혈색전증 발생의 위험이 따른다. 또한， DSA 를 이용한 측부순환의 정도를 평가할 때 지연된 정맥상 (late venous phase)까지 층분히 시행되지 않거나， 병변 반대편 혈관에 대한 평가가 누락되는 등의 경우에는 측부순환 평가가 제한적일 수 있다. 따라서， 인체에 비침습적이고 좀 더 빠른 시간 안에 쉬우면서도 객관적으로 측부순환의 정도를 알 수 있는 새로운 기법의 개발이 필요하다. 본 명세서 전체에 걸쳐 다수의 논문 및 특허문헌이 참조되고 그 인용이 표시되어 있다. 인용된 논문 및 특허문헌의 개시 내용은 그 전체로서 본 명세서에 참조로 삽입되어 본 발명이 속하는 기술 분야의 수준 및 본 발명의 내용이 보다 명확하게 설명된다.
【발명의 상세한 설명】
【기술적 과제】
본 발명자들은 상기의 종래기술의 문제점을 해결하기 위하여 예의 연구 노력한 결과， 자기공명영상 기법을 이용하여 측부순환를 영상화하고 그 정도를 평가할 수 있는 새로운 방법을 개발하여 본 발명올 완성하였다. 따라서, 본 발명의 목적은 신체내 측부순환을 영상화하는 방법을 제공하는데 있다.
본 발명의 다른 목적은 측부순환을 영상화하는 프로그램이 기록된 컴퓨터 판독 가능한 기록매체를 제공하는데 있다. 본 발명의 목적 및 장점은 하기의 발명의 상세한 설명， 청구의 범위 및 도면에 의해 보다 명확하게 된다 . 【기술적 해결방법】
본 발명의 일 양태에 따르면， 본 발명은 다음의 단계를 포함하는 신체 내의 측부순환을 영상화하는 방법을 제공한다:
(a) 신체내의 자기공명영상으로서 역동적 자화율대조조영증강 관류강조영상 (dynamic susceptibility contrast -enhanced perfusion weighted imaging, DSC-PWI)의 소스 데이터를 영상 후처리하는 단계; 및
(b) 상기 단계 (a)에서 후처리를 통해 얻은 영상을 디지털 감산 혈관 조영술 (digital subtraction angiography, DSA)의 영상과 대조하여 측부순환을 영상화하는 단계.
본 명세서에서 용어 "측부순환 (col lateral circulation)" 은 혈관의 협착이나 폐색 때문에 정상적인 혈류의 흐름이 차단되었을 때 혈류가 다른 혈관으로 우회하여 흐르는 현상을 의미한다. 뇌혈류의 경우 양측 경동맥에 의한 전방순환 (anterior circulat ion)과 뒤쪽의 추골 기저동맥에 의한 후부순환 (posterior circulat ion)으로 구성되며， 이들은 뇌 기저부에서 문합하여 월리스 고리 (Circle of Willis)를 이루고, 이후 전뇌동맥， 중뇌동맥， 후뇌동맥으로 분지한다. 이들 각각의 혈관은 궁극적으로 뇌 일정영역의 혈액 공급을 담당하나 또한 이들 혈관간의 문합이 몇몇 곳에서 이루어져 일측혈관의 병변 시 측부순환을 이루게 된다.
본 발명의 바람직한 구현예에 의하면, 본 발명에서 측부순환은 뇌， 심장， 사지말단, 신장， 및 장 (intestine) 부위에서의 측부순환을 포함하며， 보다 바람직하게는 본 발명에서 측부순환은 뇌경색 또는 심근경색시에 발생되는 측부순환을 포함한다.
본 명세서에서 용어 "뇌경색" 은 뇌의 혈관이 막혀 뇌의 일부가 죽는 질병을 의미한다. 뇌경색은 뇌혈관의 동맥경화증에 의한 혈전 및 심장에서 유래된 색전에 의한 뇌혈관의 폐색이 원인이다.
본 명세서에서 용어 "자기공명영상 (magnetic resonance imaging)" 은 체내 구조물의 사진을 생성하기 위해서 자기장을 사용하는 진단기술을 의미한다. 자기공명영상에 사용하는 기계는 전자파의 펄스를 이용하여 신체에 매우 강력한 자기장을 형성하며, 이러한 자기장의 형성과 붕괴의 과정에서 만들어지는 신호를 이용하여 체내 구조물을 영상화한다. 본 명세서에서 용어 "역동적 자화율대조조영증강 관류강조영상 (dynamic susceptibility contrast-enhanced perfusion weighted imaging, DSC-PWI)" 는 자화감수성의 증가에 의해 야기되는 자기장의 변화를 검출하는 자기공명영상의 한 기법을 의미한다. 일반적으로는 상자성 조영제인 GD-DTPA( gadolinium diethylene triamine pent a acetic acid)를 환자 팔의 정맥을 통해 주입한 뒤, 시간에 따른 신호변화를 연속 촬영하여 관찰하고, 이러한 신호변화의 검출을 통하여 뇌혈류， 뇌혈류량， 조영제의 평균이동시간 등의 다양한 뇌혈류 -유래 지도들을 영상화하는 기법이다.
본 명세서에서 용어 "영상 후처리 (image post-processing)" 는 영상의 획득 후에 상기 영상으로 이루어진 정보들을 소스 데이터로 이용하여 보정 및 재구성을 거쳐 새로운 영상으로 출력해내는 일련의 작업들을 통칭하는 의미의 용어이다. 일반적으로 의료영상은 3차원 함수 f(x, y, z)로 정의될 수 있는데, X, y, z는 공간 좌표를 의미하며， 그 좌표에서의 f의 크기를 영상의 신호세기 (signal intensity) 라고 부른다. 의료영상은 디지털 영상이기 때문에 유한한 요소들로 구성되어 있고, 각 요소는 특정한 위치와 값을 가지는데, 이 요소를 복셀 (voxel)이라고 부른다. 각 복셀에서의 신호세기 값들은 영상 후처리 과정을 통하여 보정 및 재구성 되고， 새로운 정보들로 구성된 영상으로 출력된다.
본 명세서에서 용어 "디지털 감산 혈관조영술 ((digital substraction angiography, DSA)" 이란 혈관조영술의 한 방법을 의미한다. 구체적으로 설명하면， 먼저 혈관조영을 할 부분의 단순 X선 촬영을 실시하고 다음에 혈관속에 X선 촬영시 관찰 가능한 조영제를 주입하여 X선 촬영을 시행한다. 이렇게 두 가지의 촬영 정보를 컴퓨터에 입력시켜 두 번째 촬영의 정보에서 첫 번째 촬영의 정보를 빼도록 한다. 이러한 과정을 거치게 되면 혈관만 깨끗이 조영되고 나머지 부분의 영상은 없어진 영상을 얻을 수 있다. 이러한 방식으로 혈관 부분의 영상을 더욱 선명히 얻는 기법을 디지털 감산 혈관조영술이라고 한다.
본 발명의 바람직한 구현예에 의하면， 상기 단계 (a)의 DSC-PWI 소스 데이터의 영상 후처리는 조영제 주입 전에 얻은 영상의 신호를 조영제 주입 후에 얻은 영상들의 신호들로부터 차감하여 시간 경과에 따른 순차적인 영상을 얻는 단계 (a)-l를 포함한다. 상기와 같이 조영제 주입 전에 얻은 영상의 신호를 조영제 주입 후에 얻은 시간에 따른 순차적인 영상들의 신호들로부터 차감함으로써， 해부학적 구조에 의한 신호를 배제하고 조영제에 의한 신호만을 남겨 영상화하고자 하는 조직을 순차적으로 시각화할 수 있다. 상기 DSC-PWI의 시간 경과에 따른 순차적인 영상들은 일정의 시간 간격을 갖는 순차적인 영상을 의미하며, 상기 일정의 시간 간격은 바람직하게는 0.1 ᅳ 10초 이고， 보다 바람직하게는 0.5 ᅳ 5초， 보다 더 바람직하게는 1.0 - 3초, 가장 바람직하게는 1.5 - 2초 이다.
본 발명의 다른 바¾직한 구현예에 의하면， 상기 단계 (a)의 DSC-PWI 소스 데이터의 영상 후처리는 상기 단계 (a)-l 에서 얻은 순차적인 영상들에서 인접하는 연속 영상 컷을 중첩하는 단계 (a)-2 를 포함한다. 본 발명에서 DSC-PWI의 순차적인 영상들에서 인접하는 연속 영상 컷을 중첩시킴으로써 영상의 공간 해상도를 향상시키고 시간에 따른 신호 강도 변화를 최대화시킨다. 상기 중첩시키는 연속 영상 컷의 개수는 공간 해상도를 향상시키면서 시간에 따른 신호 강도 변화를 최대화시킬 수 있는 범위 내에서 적합하게 선택할 수 있다. 예를 들어， 상기 인접하는 연속 영상 컷의 중첩은 2 - 10장의 인접하는 연속 영상 컷을 중첩하여 행할 수 있으며， 바람직하게는 2 - 7장， 보다 바람직하게는 2 - 5장， 보다 더 바람직하게는 2 - 3장, 가장 바람직하게는 2장의 인접하는 연속 영상 컷을 중첩시킨다.
본 발명의 다른 바람직한 구현예에 의하면， 상기 단계 (b)에서 후처리한 DSC-PWI 영상과 DSA 영상의 대조는 동맥기 (arterial phase) , 모세혈관기 (capillary phase), 정맥기 (venous phase), 또는 늦은 정맥기 (late venous phase)에 해당하는 영상에 대해 행한다.
본 명세서에서 용어 "동맥기" 는 정상측 뇌 실질의 뇌섬부위 (insular region)로부터 조영제 염색이 시작되는 시기, 용어 "모세혈관기" 는 정상측 뇌 실질 전체가 최고의 신호 강도로 염색되는 시기, 용어 "정맥기" 는 정상측 뇌 실질에서 조영제가 유실 (washout)되는 시기， 용어 "늦은 정맥기" 는 정상측 뇌 실질은 조영제가 모두 유실된 것에 반하여 병변측 뇌 실질에서 남은 조영제 신호가 보이는 시기를 각각 의미한다. 본 발명의 다른 바람직한 구현예에 의하면， 상기 동맥기， 모세혈관기， 정맥기， 및 늦은 정맥기에서 DSC-PWI 영상은 상기 설명된 "인접하는 연속 영상 컷이 중첩된 이미지" 의 연속된 영상으로 구성된다. 보다 바람직하게는， 상기 동맥기， 정맥기， 늦은 정맥기에서 DSC-PWI 영상은 상기 설명된 "인접하는 연속 영상 컷이 중첩된 이미지" 의 연속된 2장의 영상으로 구성되며， 상기 모세혈관기에서 DSC-PWI 영상은 상기 설명된 "인접하는 연속 영상 컷이 증첩된 이미지" 의 연속된 3장의 영상으로 구성된다.
본 발명의 다른 바람직한 구현예에 의하면， 상기 동맥기를 구성하는 DSC-PWI 영상을 중첩하여 측부순환의 초기 (early phase) 영상을 구성하고， 상기 모세혈관기를 구성하는 DSC— PWI 영상을 중첩하여 측부순환의 증기 (mid phase) 영상을 구성하며， 상기 정맥기 및 늦은 정맥기를 구성하는 DSC-PWI 영상을 중첩하여 측부순환의 말기 (late phase) 영상을 구성한다. 본 발명의 다른 양태에 따르면， 본 발명은 다음의 단계를 포함하는 측부순환영상을 제작하는 방법을 제공한다:
(a) ' 역동적 자화율대조조영증강 관류강조영상 (DSC-PWI)의 소스 데이터에서, 해당 시간에서의 자기공명 신호강도 S(t), 베이스 라인의 자기공명 신호강도 S0, echo time TE를 하기의 수학식에 대입하여 각 복샐 (voxel)의 값을 얻는 단계;
¾*(t) = ᅳ 1/TE log [S(t)/S0]
(b) ' 각 영상 컷에서 모든 복셀의 R2 * 값을 더하고， 더하여 합친 R2* 값의 시간에 따른 변화를 그래프화 하는 단계;
(c) ' 상기 그래프에서 합친 * 값의 최대치를 찾고, 이 최대치에 해당하는 영상 컷을 측부순환영상의 중기 (mid phase) 중간지점으로 정하는 단계; 및
(d) ' 상기 측부순환영상의 증기 (mid phase) 중간지점으로 정한 영상 컷올 기준으로 인접한 영상을 중첩하여 측부순환영상을 작성하는 단계. 본 발명의 또 다른 양태에 따르면， 본 발명은 다음의 단계를 포함하는 측부순환영상을 제작하는 프로그램이 기록된 컴퓨터 판독 가능한 기록매체를 제공한다:
(a) " 역동적 자화율대조조영증강 관류강조영상 (DSC-PWI)의 소스 데이터에서， 해당 시간에서의 자기공명 신호강도 S(t), 베이스 라인의 자기공명 신호강도 S0, echo time TE을 하기의 수학식에 대입하여 각 복셀 (voxel)의 * 값을 얻는 단계;
R2*(t) = - 1/TE log [S(t)/S0]
(b) " 각 영상 컷에서 모든 복샐의 * 값올 더하고, 더하여 합친 * 값의 시간에 따른 변화를 그래프화 하는 단계;
(c) " 상기 그래프에서 합친 * 값의 최대치를 찾고, 이 최대치에 해당하는 영상 컷을 측부순환영상의 중기 (mid phase) 중간지점으로 정하는 단계; 및
(d) " 상기 측부순환영상의 중기 (mid phase) 중간지점으로 정한 영상 컷을 기준으로 인접한 영상을 중첩하여 측부순환영상을 작성하는 단계. 상기 단계 (d)' 및 (d) 에서 인접한 영상을 중첩하여 측부순환영상을 작성하는 단계는 상기 본 발명의 다른 양태인 신체내 측부순환을 영상화하는 방법에서 설명된 방법과 동일하므로 중복하여 설명하지 않는다. 상기 설명된 본 발명의 방법으로부터 얻은 측부순환영상에서 관찰되는 측부순환의 등급을 기초로 혈관 재개통 후의 뇌경색 환자의 예후를 예측할 수 있다 (실시예 6 참조). 본 발명의 특징 및 이점을 요약하면 다음과 같다:
(i) 본 발명은 신체내에서 측부순환을 신속하고 정확하게 영상화하는 새로운 자기공명영상 기법에 관한 것이다.
(ii) 본 발명은 역동적 자화율대조조영증강 관류강조영상 (DSC-PWI)의 소스 데이터를 영상 후처리한 후 이렇게 후처리를 통해 얻은 영상을 디지털 감산 혈관 조영술 (DSA)의 영상과 대조하여 측부순환을 영상화한다.
(iii) 본 발명의 방법에 의하면， 측부순환을 비침습적으로 신속하면서도 용이하게 영상화하여 평가할 수 있으므로 혈관 협착이나 폐색에 의해 발생되는 질환의 이해, 치료 방침 결정 및 예후 예측에 중요한 정보를 제공할 수 있다.
【유리한 효과】
본 발명은 신체내 측부순환을 영상화하는 새로운 자기공명영상 기법에 관한 것이다. 본 발명은 역동적 자화율대조조영증강 관류강조영상 (DSC-PWI)의 소스 데이터를 영상 후처리한 후 후처리를 통해 얻은 영상을 디지털 감산 혈관 조영술 (DSA)의 영상과 대조하여 측부순환을 영상화한다. 본 발명의 방법에 의하면 측부순환을 비침습적으로 신속하면서도 용이하게 영상화하여 평가할 수 있으므로 혈관 협착이나 폐색에 의해 발생되는 질환의 이해， 치료 방침 결정 및 예후 예측에 중요한 정보를 제공할 수 있다. 【도면의 간단한 설명】
도 1 은 우측 중대뇌등맥의 폐색이 동반된 뇌경색 환자에서 역동적 자화율대조조영증강 관류강조영상 (dynamic susceptibility contrast- enhanced perfusion weighted imaging, DSC-PWI ) 소스 데이터를 이용하여 뇌 실질을 시각화한 예를 보여준다. DSC-PWI 의 소스 데이터는 한 슬라이스 당 50장의 이미지 컷으로 구성되며， 각 컷 사이에 약 1.7초의 시간 간격을 갖는다. 각 이미지 하단의 시간은 소스 데이터를 얻기 시작한 시점으로부터 경과된 시간을 의미한다. 우측 중대뇌동맥이 폐색되어 있기 때문에 좌반구에 비해 우반구의 조영제 염색 (staining)과 유실 (washout)이 상대적으로 느려져 있음을 알 수 있다.
도 2 는 측부순환 영상화 과정에서 최상의 시간 /공간 해상도를 얻기 위해 DSC-PWI 이미지 컷을 중첩한 결과를 보여준다. 서로 인접한 2 컷의 이미지를 서로 중첩할 때 최상의 시간 /공간 해상도를 얻을 수 있다. 이에 비해 중첩하지 않은 이미지는 공간 해상도가 떨어지며 이보다 많은 컷을 중첩하는 경우에는 시간 해상도가 떨어지게 된다.
도 3 은 디지털 감산 혈관조영술 영상과 DSC-PWI 영상을 매칭하여 측부순환 영상을 구성하는 과정을 보여준다. 첫 번째 행 (line)은 병변 반대 측에서 시행한 디지털 감산 혈관조영술에 의해 얻은 영상으로서 각각 동맥기 (arterial phase), 모세혈관기 (capi 1 lary phase), 정맥기 (venous phase)로 구성된다. 두 번째 행은 병변 측에서 시행한 디지털 감산 혈관조영술에 의해 얻은 영상으로서， 우측 중대뇌동맥 근위부의 폐색을 확인할 수 있다 (thin white arrow). 또한 모세혈관기부터 우측 전대뇌동맥에서 내려오는 측부순환을 관찰할 수 있다 (thick white arrows). 측부순환은 늦은 정맥기 (late venous phase)까지 관찰 가능하다. 세 번째 행은 디지털 감산 혈관조영술 영상의 각 시기에 해당하는 DSC— PWI 영상의 2 컷 중첩이미지를 보여준다. 병변 측 반구에서 측부순환에 의한 조영제의 염색 (staining)과 유실 (washout)을 확인할 수 있다 (black arrow heads). 네 번째 행은 측부순환 영상의 생성을 보여준다. 동맥기, 모세혈관기， 및 정맥기 -늦은 정맥기에 해당하는 DSC-PWI 의 2 컷 중첩이미지를 다시 한번 중첩하여 측부순환을 평가할 수 있는 영상을 생성하였다. 각 시기는 초기 (early phase), 중기 (mid phase) 및 말기 (late phase)에 해당한다.
도 4는 자동화 소프트웨어 FAST-C0LL을 이용한 측부순환영상의 작성 과정을 보여준다. 환자가 뇌경색 증상이 발생하여 병원에 내원한 후 측부순환영상을 생성하기까지의 전체적인 흐름을 보여준다. 뇌 자기공명영상 시행 후 영상을 워크스테이션으로 보내면 FAST-C0LL 을 이용하여 측부순환영상을 생성하고 그 결과를 바탕으로 측부순환의 정도를 평가한다. 생성된 측부순환영상은 의료영상저장전송시스템 (Picture Archiving Co麵 unication System, PACS)에 보내져 다른 멀티모달 자기공명영상에서 얻어진 정보들과 함께 환자에 대한 다양한 정보를 제공하게 된다. FAST-C0LL 을 이용하면 자기공명영상 소스데이터를 얻은 후 평가까지 약 5분 정도의 시간이 소요된다.
도 5 는 측부순환영상에서 측부순환 정도를 평가한 결과를 보여준다. 패널 (1)은 빈약한 측부순환 (Poor collaterals)으로 평가되는 영상으로서， 병변측의 조영제 염색이 말기 (late phase)까지 완전히 없거나 일부 지역만 느린 측부순환으로 염색되는 경우이다. 패널 (2)는 중간 정도의 측부순환 (Intermediate col laterals)으로 평가되는 영상으로서 , 병변측의 조영제 염색이 일부에서 완전히 없으면서, 일부 지역은 빠른 측부순환으로 염색되는 경우이다. 패널 (3)은 좋은 측부순환 (Good collaterals)으로 평가되는 영상으로서, 병변측의 조영제 염색이 느린 측부순환으로 완전히 염색되는 경우이다. 패널 (4)는 탁월한 측부순환 (Excellent col laterals)으로 평가되는 영상으로서 , 조영제 염색이 빠른 측부순환으로 완전히 염색되는 경우이다.
도 6 은 DSC-PWI 소스 데이터의 각 컷에서 모든 복셀의 R2* 값을 더하여 합친 R2 * 값의 시간에 따른 변화를 보여주는 그래프이다. 상기 * 값은 R2 *(t) = - 1/TE log [S(t)/S0] 식에 의해 산출되며， 상기 식에서 S(t)는 해당 시간에서의 자기공명 신호강도이고， S0 는 베이스 라인의 자기공명 신호강도이며， TE는 echo time이다.
도 7a 내지 도 7d 는 측부순환영상에서 관찰되는 측부순환의 등급 및 혈관 재개통 여부에 따른 뇌경색 후 예후를 보여주는 '그래프이다. ER (early recanalizat ion)은 혈관 재개통, Coll (col laterals)은 측부순환, G (grade)는 등급을 의미한다. 측부순환의 각 등급 Gl, G2, G3 및 G4 는 빈약한 측부 순환 (poor collaterals), 중간 정도의 측부 순환 ( intermediate collaterals), 좋은 측부 순환 (good collaterals), 및 탁월한 측부 순환 (excel lent col laterals) (실시예 5 및 도 5 참조)에 각각 대웅된다. 도 7a 에서 보여지는 바와 같이 측부순환의 등급이 높아질수록 좋은 기능적 예후를 가질 확률이 높아진다. 도 7b 에서 보여지는 바와 같이 측부순환의 등급이 높아질수록 좋은 신경학적 호전을 보일 확률이 높아진다. 도 7c에서 보여지는 바와 같이 측부순환의 등급이 높아질수록 증상성 뇌출혈이 발생할 확률은 낮아진다. 도 7d 에서 보여지는 바와 같이 측부순환의 등급이 높아질수특 사망률은 낮아진다.
【발명의 실시를 위한 형태】
이하， 실시예를 통하여 본 발명을 더욱 상세히 설명하고자 한다. 이들 실시예는 오로지 본 발명을 보다 구체적으로 설명하기 위한 것으로, 본 발명의 요지에 따라 본 발명의 범위가 이들 실시예에 의해 제한되지 않는다는 것은 당업계에서 통상의 지식을 가진 자에 있어서 자명할 것이다. 실시예 실시예 1: DSC-PWI 에서 가돌리늄 조영제에 의한 뇌 실질의 순차적 시각화
측부순환을 영상화하기 위한 첫 번째 단계로 DSC-PWI의 소스 데이터를 이용하여 뇌 실질이 조영제에 의해 염색 (staining) 되었다가, 유실 (washout) 되는 양상을 순차적으로 시각화 하였다 (도 1). 이를 위하여 DSC-PWI의 소스 데이터의 각 슬라이스에서 얻어지는 50장의 이미지 중에서 첫 번째 컷의 신호 (signal)를 나머지 49장의 신호에서 각각 빼주었다. 이렇게 함으로써 해부학적 구조에 의한 신호를 배제하고 가돌리늄 조영제의 효과만을 남기게 되어 뇌 실질을 순차적으로 시각화하는 것이 가능하였다. 실시예 2: 최상의 시간 및 공간 해상도를 얻기 위한 이미지 중첩 측부순환을 영상화하기 위한 두 번째 단계로 DSCᅳ PWI에서 얻은 이미지를 단계적으로 중첩하여 시간 및 공간 해상도를 서로 비교하였다 (도
2). 이를 위하여 상기 실시예 1의 과정을 통해 얻어진 이미지의 연속된 컷을 2장， 3장, 5장， 7장， 10장씩 증첩해 보았다. 증첩하기 전의 이미지가 가장 뛰어난 시간 해상도를 제공하지만 한 컷 안에서 양쪽 반구 및 인접한 컷 사이의 차이가 중첩 이미지에 비해 덜 선명한 단점이 있다. 이에 반해
5장 이상의 컷을 중첩한 이미지는 뛰어난 공간 해상도를 제공하지만 시간에 따른 신호 강도의 변화를 정확히 감지하지 못한다. 남은 중첩 이미지 중에서 2장을 중첩한 이미지가 3장을 증첩한 이미지에 비해 더 좋은 시간 해상도를 가지면서 두 이미지의 공간 해상도는 비슷하였다. 따라서, 최종적으로 측부순환 영상화에는 2장의 인접한 컷을 중첩한 이미지를 이용하기로 결정하였다. 실시예 3: 측부순환 영상 생성
측부순환을 영상화하기 위한 마지막 단계로 실시예 2의 과정에서 결정된 2컷 증첩 이미지를 디지털 감산 혈관조영술에서 얻어진 이미지와 매칭시켰다 (도 3). 이를 위하여 2컷 중첩 이미지에서 정상측 및 병변측 뇌 실질의 조영제 염색 (staining)과 유실 (washout) 정도를 시각적으로 판별하여 디지털 감산 혈관조영술의 동맥기 (arterial phase), 모세혈관기 (capillary phase), 정맥기 (venous phase)에 해당하는 이미지를 선정하였다. 구체적으로 동맥기는 정상측 뇌 실질의 뇌섬부위 (insular region)로부터 조영제 염색 (staining)이 시작되는 시기， 모세혈관기는 정상측 뇌 실질 전체가 최고의 신호 강도로 염색되는 시기， 정맥기는 정상측 뇌 실질에서 조영제가 유실 (washout)되는 시기이며， 늦은 정맥기 (late venous phase)는 정상측 뇌 실질은 조영제가 모두 유실된 것에 반하여, 병변측 뇌 실질에서 남은 조영제 신호가 보이는 시기이다. 전형적으로는 상기 각 시기는 모세혈관기에 3장， 나머지의 동맥기 정맥기 및 늦은 정맥기에는 2장의 연속된 2컷의 중첩 이미지로 구성된다. 측부순환 영상은 동맥기의 2컷 중첩이미지 2장, 모세혈관기의 2컷 중첩이미지 3장， 정맥 및 늦은 정맥기의 2컷 중첩이미지 4장을 각각 다시 한번 중첩하여 이미지를 생성하였다. 이렇게 생성된 측부순환 영상의 각 시기는 초기 (early phasae), 중기 (mid phase) 및 말기 ate phase)에 해당한다. 실시예 4: 자동화 소프트웨어를 이용한 측부순환영상의 작성
상기 실시예 3에서 설명된 측부순환영상을 얻는 방법과 동일한 논리적 알고리즘을 적용하여 자동적으로 측부순환영상을 생성하였다. 이러한 자동화 소프트웨어를 소프트웨어 FAST-C0LL(Fast Analysis System for Collaterals)로 명명하였다. FAST-COLL를 이용한 측부순환영상 생성의 전체적인 흐름은 도 4에 나타내었다. FAST-C0LL은 DSC— PWI의 소스 데이터에서 각 복셀의 ¾ * 값을 계산하며， 그 계산식은 아래와 같다.
^¾ (^)
Figure imgf000013_0001
여기에서 S(t)는 해당되는 시간에서 자기공명 신호강도를 의미하며, So는 베이스라인 자기공명 신호강도, TE는 echo time을 의미한다. 측부순환영상의 생성을 위해서 각 컷에서 모든 복셀의 * 값을 더하고， 더하여 합친 * 값의 시간에 따른 변화를 그래프화 하였다 (도 6). 아 그래프에서 합친 R2 * 값의 최대치를 자동으로 찾아내고， 최대치에 해당하는 컷을 측부순환영상의 중기 (mid phase) 중간지점으로 설정하였다. 이 컷을 기준으로 실시예 3에서와 같이 인접한 이미지를 중첩하여 측부순환영상을 생성하였다. 실시예 5: 측부순환영상에서 측부순환 정도의 평가
측부순환영상을 생성한 후 생성한 영상으로부터 측부순환의 정도를 평가할 수 있다. 측부순환의 정도는 크게 측부순환 없음， 느린 측부순환 및 빠른 측부순환으로 나눌 수 있다. 측부순환 없음은 말기 (late phase)까지 조영제 염색이 전혀 보이지 않는 경우이고， 느린 측부순환은 중기 (mid phase)에서는 보이지 않지만 말기 (late phase)에서 조영제의 염색이 보이는 경우이며， 빠른 측부순환은 정상측과 비슷한 정도로 병변측에서도 중기 (mid phase)에서부터 조영제의 염색이 보이는 경우를 지칭한다. 이를 바탕으로 측부 순환을 평가할 수 있는데， 도 5의 측부순환영상에서 나타나는 바와 같이 , 병변측의 조영제 염색이 말기 (late phase) 까지 완전히 없거나 일부 지역만 느린 측부순환으로 염색되는 경우를 빈약한 측부 순환 (poor collaterals), 병변측의 조영제 염색이 일부에서 완전히 없으면서， 일부 지역은 빠른 측부순환으로 염색되는 경우를 중간 정도의 측부 순환 (intermediate collaterals), 병변측의 조영제 염색이 느린 측부순환으로 완전히 염색되는 경우를 좋은 측부 순환 (good collaterals), 조영제 염색이 빠른 측부순환으로 완전히 염색되는 경우를 탁월한 측부 순환 (excel lent col laterals)으로 평가할 수 있다. 실시예 6: 측부순환영상에서 얻어진 측부순환 등급에 따라 급성 뇌경색 환자의 혈관 재개통 치료 후 예후 예측
측부순환영상에서 얻어지는 측분순환의 등급에 따라 급성 뇌경색 환자에 [서 혈관 재개통 치료 후 혈관 재개통 여부에 따른 예후를 예측할 수 있다. 도 7a 내지 도 7d에서 나타나는 바와 같이， 급성 뇌경색에서 혈관 재개통의 여부와 상관없이， 측부 순환의 정도가 좋을수록 좋은 예후를 보일 확률은 높아지고, 증상성 출혈 및 사망이 발생할 확률은 낮아진다. 즉， 측부순환영상에서 관찰되는 측부순환의 등급이 높은 경우에 혈관이 재개통 되면 증상성 뇌출혈의 위험 없이 좋은 예후를 기대해 볼 수 있다. 이상으로 본 발명의 특정한 부분을 상세히 기술하였는 바， 당업계의 통상의 지식을 가진 자에게 있어서 이러한 구체적인 기술은 단지 바람직한 구현 예일 뿐이며， 이에 본 발명의 범위가 제한되는 것이 아닌 점은 명백하다. 따라서， 본 발명의 실질적인 범위는 첨부된 청구항과 그의 등가물에 의하여ᅵ정의된다고 할 것이다.

Claims

【청구의 범위】
【청구항 1】
다음의 단계를 포함하는 신체 내 측부순환을 영상화하는 방법:
(a) 신체내의 자기공명영상으로서 역동적 자화율대조조영증강 관류강조영상 (dynamic susceptibility contrast— enhanced perfusion weighted imaging, DSOPWI)의 소스 데이터를 영상 후처리하는 단계; 및
(b) 상기 단계 (a)에서 영상 후처리를 통해 얻은 영상을 디지털 감산 혈관 조영술 (digital subtraction angiography, DSA)의 영상과 대조하여 측부순환을 영상화하는 단계.
【청구항 2】
제 1 항에 있어서， 상기 단계 (a)의 DSC-PWI 소스 데이터의 영상 후처리는 조영제 주입 전에 얻은 영상의 신호를 조영제 주입 후에 얻은 영상들의 신호들로부터 차감하여 시간 경과에 따른 순차적인 영상을 얻는 단계 (a)-l를 포함하는 것을 특징으로 하는 방법 .
【청구항 3】
제 2 항에 있어서， 상기 단계 (a)의 DSC-PWI 소스 데이터의 후처리는 상기 단계 (a)-l 에서 얻은 순차적인 영상들에서 인접하는 영상 컷을 중첩하는 단계 (a)-2를 포함하는 것을 특징으로 하는 방법 .
【청구항 4】
제 3 항에 있어서， 상기 인접하는 연속 영상 컷의 중첩은 2 - 10장의 영상컷을 중첩하여 행하는 것을 특징으로 하는 방법 .
【청구항 5】
제 1 항에 있어서， 상기 단계 (b)에서 후처리한 DSC-PWI 영상과 DSA 영상의 대조는 동맥기， 모세혈관기, 정맥기, 또는 늦은 정맥기에 해당하는 영상에 대해 행하는 것을 특징으로 하는 방법.
【청구항 6】 제 5 항에 있어서， 상기 동맥기, 모세혈관기， 정맥기， 및 늦은 정맥기에서 DSC-PWI 영상은 인접하는 연속 영상 컷이 중첩된 이미지의 연속된 영상으로 구성되는 것을 특징으로 하는 방법.
【청구항 7】
제 6 항에 있어서， 상기 동맥기를 구성하는 DSC-PWI 영상올 중첩하여 측부순환의 초기 (early phase) 영상을 구성하고， 상기 모세혈관기를 구성하는 DSC-PWI 영상을 중첩하여 측부순환의 중기 (mid phase) 영상을 구성하며， 상기 정맥기 및 늦은 정맥기를 구성하는 DSC-PWI 영상을 중첩하여 측부순환의 말기 (late phase) 영상을 구성하는 것을 특징으로 하는 방법 .
【청구항 8】
제 1 항에 있어서, 상기 측부순환은 뇌 또는 심장 부위에서의 측부순환인 것을 특징으로 하는 방법.
【청구항 9】
다음의 단계를 포함하는 측부순환영상을 제작하는 방법:
(a) ' 역동적 자화율대조조영증강 관류강조영상 (DSC-PWI)의 소스 데이터에서， 해당 시간에서의 자기공명 신호강도 S(t), 베이스 라인의 자기공명 신호강도 S0, echo time TE 를 하기의 수학식에 대입하여 각 복셀 (voxel)의 * 값을 얻는 단계;
R2 *(t) = - 1/TE log [S(t)/S0]
(b) ' 각 영상 컷에서 모든 복셀의 R2 * 값을 더하고， 더하여 합친 * 값의 시간에 따른 변화를 그래프화 하는 단계 ;
(c) ' 상기 그래프에서 합친 R2 * 값의 최대치를 찾고， 이 최대치에 해당하는 영상 컷을 측부순환영상의 중기 (mid phase) 중간지점으로 정하는 단계; 및
(d) ' 상기 측부순환영상의 중기 (mid phase) 중간지점으로 정한 영상 컷을 기준으로 인접한 영상을 중첩하여 측부순환영상을 작성하는 단계. 【청구항 10】
다음의 단계를 포함하는 측부순환영상을 제작하는 프로그램이 기록된 컴퓨터 판독 가능한 기록매체:
(a) " 역동적 자화율대조조영증강 관류강조영상 (DSC-PWI)의 소스 데이터에서, 해당 시간에서의 자기공명 신호강도 S(t), 베이스 라인의 자기공명 신호강도 S0, echo time TE 을 하기의 수학식에 대입하여 각 복셀 (voxel)의 값을 얻는 단계;
R2 *(t) = - 1/TE log [S(t)/S0]
(b) " 각 영상 컷에서 모든 복셀의 R2 * 값을 더하고， 더하여 합친 R2 * 값의 시간에 따른 변화를 그래프화 하는 단계 ;
(c) " 상기 그래프에서 합친 R2 * 값의 최대치를 찾고, 이 최대치에 해당하는 영상 컷을 측부순환영상의 중기 (mid phase) 중간지점으로 정하는 단계 ; 및
(d) " 상기 측부순환영상의 중기 (mid phase) 중간지점으로 정한 영상 컷을 기준으로 인접한 영상을 중첩하여 측부순환영상을 작성하는 단계.
PCT/KR2014/003346 2013-06-03 2014-04-17 측부순환의 영상화 및 평가를 위한 새로운 자기공명영상기법 WO2014196732A1 (ko)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US14/895,597 US20160125597A1 (en) 2013-06-03 2014-04-17 Novel magnetic resonance image technique for imaging and evaluating collateral circulation
EP14807402.4A EP3005943B1 (en) 2013-06-03 2014-04-17 Novel magnetic resonance image technique for imaging and evaluating collateral circulation
JP2016518256A JP6240319B2 (ja) 2013-06-03 2014-04-17 側副循環の画像化及び評価のための新しい磁気共鳴画像法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR10-2013-0063423 2013-06-03
KR20130063423 2013-06-03

Publications (1)

Publication Number Publication Date
WO2014196732A1 true WO2014196732A1 (ko) 2014-12-11

Family

ID=52008332

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/KR2014/003346 WO2014196732A1 (ko) 2013-06-03 2014-04-17 측부순환의 영상화 및 평가를 위한 새로운 자기공명영상기법

Country Status (5)

Country Link
US (1) US20160125597A1 (ko)
EP (1) EP3005943B1 (ko)
JP (2) JP6240319B2 (ko)
KR (1) KR101568826B1 (ko)
WO (1) WO2014196732A1 (ko)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6240319B2 (ja) * 2013-06-03 2017-11-29 サムスン ライフ パブリック ウェルフェア ファウンデーションＳａｍｓｕｎｇ Ｌｉｆｅ Ｐｕｂｌｉｃ Ｗｅｌｆａｒｅ Ｆｏｕｎｄａｔｉｏｎ 側副循環の画像化及び評価のための新しい磁気共鳴画像法
KR101797141B1 (ko) * 2016-04-18 2017-11-13 성균관대학교산학협력단 자기 공명 영상 처리 장치 및 그 영상 처리 방법
CN106821364B (zh) * 2017-02-08 2020-07-03 中国科学院心理研究所 无创检测颈内外侧支循环的设备、非诊断方法及应用
KR102001790B1 (ko) * 2018-12-24 2019-07-23 (주)제이엘케이인스펙션 인공지능 기반 혈류 구간 분류 방법 및 시스템
WO2020207597A1 (en) * 2019-04-12 2020-10-15 Brainlab Ag Frameless anatomy-based 2d/3d image registration
KR102272741B1 (ko) * 2019-07-11 2021-07-02 가톨릭대학교 산학협력단 4차원 자기공명 혈관조영술의 영상정보 후처리를 통한 3차원 감산 동맥조영술과 3차원 감산 정맥조영술 및 4차원 컬러 혈관조영술의 동시 구현 방법과 의료영상 시스템
US11145057B2 (en) * 2019-11-05 2021-10-12 Siemens Healthcare Gmbh Assessment of collateral coronary arteries
KR102613857B1 (ko) * 2021-11-03 2023-12-14 사회복지법인 삼성생명공익재단 딥러닝 모델 기반의 동적 관류 mri 영상을 이용한 측부 순환 평가 방법 및 분석장치

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005095340A (ja) * 2003-09-24 2005-04-14 Toshiba Corp 血流解析装置及び血流解析方法
JP5171251B2 (ja) * 2004-05-04 2013-03-27 スティフテルセン ウニヴァズィテッフォースクニング ベルゲン 灌流磁気共鳴画像化装置の作動方法、灌流磁気共鳴画像化データ操作方法及びこれらの方法に用いるデータ処理装置、コンピュータソフトウェア製品
JP2013085652A (ja) * 2011-10-17 2013-05-13 Toshiba Corp 医用画像処理システム

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4379559B2 (ja) * 1999-04-30 2009-12-09 株式会社日立メディコ 磁気共鳴イメージング装置
CN103635136B (zh) * 2011-04-29 2016-05-25 磁共振创新有限公司 组织相似度图
JP6240319B2 (ja) * 2013-06-03 2017-11-29 サムスン ライフ パブリック ウェルフェア ファウンデーションＳａｍｓｕｎｇ Ｌｉｆｅ Ｐｕｂｌｉｃ Ｗｅｌｆａｒｅ Ｆｏｕｎｄａｔｉｏｎ 側副循環の画像化及び評価のための新しい磁気共鳴画像法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005095340A (ja) * 2003-09-24 2005-04-14 Toshiba Corp 血流解析装置及び血流解析方法
JP5171251B2 (ja) * 2004-05-04 2013-03-27 スティフテルセン ウニヴァズィテッフォースクニング ベルゲン 灌流磁気共鳴画像化装置の作動方法、灌流磁気共鳴画像化データ操作方法及びこれらの方法に用いるデータ処理装置、コンピュータソフトウェア製品
JP2013085652A (ja) * 2011-10-17 2013-05-13 Toshiba Corp 医用画像処理システム

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
BANG OY ET AL., J NEUROL NEUROSURG PSYCHIATRY, vol. 79, 2008, pages 625
BANG OY ET AL., STROKE, vol. 42, 2011, pages 2235
BANG OY ET AL., STROKE, vol. 42, 2011, pages 693
See also references of EP3005943A4

Also Published As

Publication number Publication date
US20160125597A1 (en) 2016-05-05
JP6240319B2 (ja) 2017-11-29
JP2016520405A (ja) 2016-07-14
JP2017213424A (ja) 2017-12-07
KR101568826B1 (ko) 2015-11-16
EP3005943A4 (en) 2017-05-10
KR20140142652A (ko) 2014-12-12
EP3005943B1 (en) 2020-09-23
EP3005943A1 (en) 2016-04-13

Similar Documents

Publication Publication Date Title
WO2014196732A1 (ko) 측부순환의 영상화 및 평가를 위한 새로운 자기공명영상기법
US6377835B1 (en) Method for separating arteries and veins in 3D MR angiographic images using correlation analysis
Barkhausen et al. MR evaluation of ventricular function: true fast imaging with steady-state precession versus fast low-angle shot cine MR imaging: feasibility study
JP5788574B2 (ja) 磁気共鳴イメージング装置
AU2009224632B2 (en) A computer-based method and system for imaging-based dynamic function evaluation of an organ
Loubeyre et al. Dynamic contrast-enhanced MR angiography of pulmonary embolism: comparison with pulmonary angiography.
JP2002263082A (ja) Ｍｒイメージングにおいて造影強調後の定常自由歳差運動を使用する方法及び装置
Hamamoto et al. Feasibility of non-contrast-enhanced MR angiography using the time-SLIP technique for the assessment of pulmonary arteriovenous malformation
Lethimonnier et al. Three-dimensional coronary artery MR imaging using prospective real-time respiratory navigator and linear phase shift processing: comparison with conventional coronary angiography
Ley et al. MR imaging/magnetic resonance angiography of the pulmonary arteries and pulmonary thromboembolic disease
Tanabe et al. State-of-the-art pulmonary arterial imaging–Part 1
ZHANG et al. Contrast enhanced MR angiography in pulmonary sequestration
Tanaka et al. Noninvasive identification of peripheral vessels of oral and maxillofacial regions by using electrocardiography-triggered three-dimensional fast asymmetric spin-echo sequences
Crean Technical considerations for ACHD imaging
Balci et al. Assessment of the anomalous pulmonary circulation by dynamic contrast-enhanced MR angiography in under four seconds
Zeitler et al. Magnetic resonance imaging of aneurysms and thrombi
Oda et al. Magnetic resonance angiography using fresh blood imaging in oral and maxillofacial regions
Oda et al. Magnetic resonance angiography with fresh blood imaging for identification of hemangiomas and blood vessels around hemangiomas in oral and maxillofacial regions
Murakami et al. Time-of-flight MR angiography of portal system and collaterals in portal hypertension using a 2-DFT fast spoiled gradient recalled steady-state precession technique
Hamamoto et al. Non–contrast-enhanced magnetic resonance angiography using time-spatial labeling inversion pulse technique for differentiation between pulmonary varix and arteriovenous malformation
Haibo et al. Pulmonary sequestration: three dimensional dynamic contrast-enhanced MR angiography and MRI
Raptis et al. Emergency thoracic vascular magnetic resonance imaging: protocols and clinical considerations
Ristow et al. Current Imaging Strategies in Patients with Abdominal Aortic Aneurysms
Tan et al. Coronary artery disease: comprehensive evaluation by cardiovascular magnetic resonance imaging
Glockner Renal MR angiography

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14807402

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2016518256

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 14895597

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2014807402

Country of ref document: EP