JP2018075075A5 - - Google Patents
Download PDFInfo
- Publication number
- JP2018075075A5 JP2018075075A5 JP2016217307A JP2016217307A JP2018075075A5 JP 2018075075 A5 JP2018075075 A5 JP 2018075075A5 JP 2016217307 A JP2016217307 A JP 2016217307A JP 2016217307 A JP2016217307 A JP 2016217307A JP 2018075075 A5 JP2018075075 A5 JP 2018075075A5
- Authority
- JP
- Japan
- Prior art keywords
- distribution
- oxygen uptake
- uptake rate
- oxygen
- unit
- Prior art date
- Legal status (The legal status 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 status listed.)
- Granted
Links
- MYMOFIZGZYHOMD-UHFFFAOYSA-N oxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 claims description 47
- 229910052760 oxygen Inorganic materials 0.000 claims description 47
- 239000001301 oxygen Substances 0.000 claims description 47
- 210000001519 tissues Anatomy 0.000 claims description 12
- 210000003462 Veins Anatomy 0.000 claims description 11
- 210000004556 Brain Anatomy 0.000 claims description 4
- 238000002595 magnetic resonance imaging Methods 0.000 claims 2
- 238000002592 echocardiography Methods 0.000 description 10
- 230000000704 physical effect Effects 0.000 description 10
- 238000000926 separation method Methods 0.000 description 7
- 238000006243 chemical reaction Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000000034 method Methods 0.000 description 4
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 4
- 210000004369 Blood Anatomy 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 230000000875 corresponding Effects 0.000 description 2
- 238000005534 hematocrit Methods 0.000 description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006011 modification reaction Methods 0.000 description 2
- 230000001702 transmitter Effects 0.000 description 2
- 210000000601 Blood Cells Anatomy 0.000 description 1
- 102000008857 Ferritin Human genes 0.000 description 1
- 108050000784 Ferritin Proteins 0.000 description 1
- 238000008416 Ferritin Methods 0.000 description 1
- 230000037013 METABOLISM RATE Effects 0.000 description 1
- 210000004940 Nucleus Anatomy 0.000 description 1
- 210000000463 Red Nucleus Anatomy 0.000 description 1
- 210000003523 Substantia Nigra Anatomy 0.000 description 1
- 230000002490 cerebral Effects 0.000 description 1
- 108010002255 deoxyhemoglobin Proteins 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drugs Drugs 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 229920000915 polyvinyl chloride Polymers 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 230000003068 static Effects 0.000 description 1
Images
Description
送信コイル105と送信機107は送信部として機能し、送信機107が生成する高周波信号を送信コイル105に送信することにより送信コイルから高周波磁場が照射される。受信コイル106と受信機108は受信部として機能し、受信コイル106が検出した核磁気共鳴信号は受信機108を通して計算機109に送られる。なお、図2では、送信コイル105と受信コイル106とに別個のものを用いる場合を示しているが、送信コイル105と受信コイル106との機能を兼用する1つのコイルで構成してもよい。
本実施形態のMRI装置では、被検体の所定の物性値の影響を含む核磁気共鳴信号を発生させるパルスシーケンスが格納されている。所定の物性値は、被検体の酸素量特性を反映した物性値であり、例えば、磁化率やR 2 * 或いは「R 2 * とR 2 との差R’」などである。R 2 及びR 2 * は、それぞれ、真の横緩和時間T 2 の逆数、みかけの横緩和時間T 2 * の逆数であり、TEの異なる複数のエコー信号を用いて算出することができる。一般にR 2 * (T 2 * )はグラディエントエコーを用いて、R 2 (T 2 )はスピンエコーを用いて算出される。
これらの処理は、図3に示す酸素摂取率分布算出部330の各部により実現することができる。すなわち、図3(a)に示ように、酸素摂取率分布算出部330は、複素画像Iから酸素摂取率を反映した物性値分布を算出する物性値分布算出部331と、算出した物性値分布から、少なくとも2つ種類以上の生体組織に分離した組織分離画像を算出する組織分離部332と、組織分離画像のいずれかを酸素摂取率に対応する画素値に変換する酸素摂取率変換部333と、組織分離画像において、任意の画素の画素値とその周囲の画素の平均値が略等しいという条件に基づいて酸素摂取率分布を推定する酸素摂取率分布推定部334と、を備える。酸素摂取率分布算出部330は、図3(b)に示すように、酸素摂取率分布推定部334が算出した推定酸素摂取率分布に対して補正を行う補正部を備えていてもよい。図示する例では、酸素摂取率分布算出部330は、任意の静脈の酸素摂取率と脳全体の酸素摂取率の平均値が略等しいという条件に基づいて部分体積により低下した酸素摂取率を補正し、最終的な酸素摂取率分布を得る部分体積効果補正部335を備える。
計測された各エコーはkr、kp、ksを座標軸とする3次元のk空間上(メモリ空間)に配置される。このとき、一つのエコーはk空間上でkr軸に平行な1ラインを占める。このRSSGシーケンス550により得られる絶対値画像は、TEが短いエコーの場合はT1(縦緩和時間)強調画像、TEが長いエコーの場合は画素内の位相分散を反映したT 2 * 強調画像となる。
本実施形態の磁化率分布算出部331は、これらの方法を用いて定量的磁化率(QSM)分布を算出する。本実施形態では、いずれの方法を用いて定量的磁化率分布を算出してもよい。
本実施形態では、粒状成分を分離した後、静脈成分を分離するので、磁化率分布において高磁化率となる淡蒼球や赤核、黒質などの鉄タンパク(フェリチン)沈着組織を取り除いた静脈の磁化率分布を算出することができる。
ここで、任意の静脈の酸素摂取率をOEFSSS、バイナリマスクM内にある画素数をN、画素位置iにおける生体組織の酸素摂取率をOEFt(i)とすると、上述した条件(任意の静脈の酸素摂取率と脳全体の酸素摂取率の平均値が等しいという条件)は、式(8)で表すことができる。PVCは式(8)を満たすための部分体積効果補正係数である。
任意の静脈の酸素摂取率OEFSSSは、例えば、組織分離処理において、静脈として分離された画像(磁化率分布画像)から任意の静脈に相当する1ないし複数の画素を選択し、その変換処理S1103後の酸素摂取率を用いることができる。
以上の処理S1101〜S1105により、図4の酸素摂取率分布算出処理(酸素量特性算出)S1003が完了する。なお、図6に示すフローは一例であり、その一部の処理を省いたり、順番を入れ替えたりすることも本実施形態に含まれる。例えば、図6の部分体積効果補正を省いたり、部分体積効果補正を組織分離画像算出後や酸素摂取率変換後に実施したりすることもあり得る。また部分体積効果補正以外に、適宜フィルタなどを用いた平滑処理などを施してもよい。さらに、算出した酸素摂取率を用いて、その他の酸素量特性、例えば酸素飽和度や酸素代謝率などを算出することも可能である。
[画像表示:S1004]
酸素摂取率分布算出部330が算出した酸素摂取率分布や磁化率分布は、ディスプレイ110(図2)に表示することができる。或いは外部記憶装置111に画像データとして格納し、所望の表示装置で表示されてもよい。表示の形態は、特に限定されず、例えば酸素摂取率分布は、図8(c)に示すような等高線やカラーマップとして表示することができる。
酸素摂取率分布算出部330が算出した酸素摂取率分布や磁化率分布は、ディスプレイ110(図2)に表示することができる。或いは外部記憶装置111に画像データとして格納し、所望の表示装置で表示されてもよい。表示の形態は、特に限定されず、例えば酸素摂取率分布は、図8(c)に示すような等高線やカラーマップとして表示することができる。
本実施形態においても計算機109及び酸素摂取率算出部330の構成は、第一実施形態と同様であるが、計測制御部310は2つ以上の複数のエコー時間のエコーを計測し、エコー時間が異なる複数の複素画像を取得する制御を行う。また物性値分布算出部331、及び酸素摂取率変換部333の処理内容が異なる。具体的には、また物性値分布算出部331では、複数のエコー時間の複素画像から、R2 *分布を算出し、酸素摂取率変換部333では、R2 *分布の各画素値を酸素摂取率に変換する。本実施形態では、図3における物性値分布算出部331は、R2 *分布算出部と読み替える。
[酸素摂取率変換:S1103]
酸素摂取率変換部333は、S1102で分離した静脈のR2 *分布の各画素値(R2 *値)を酸素摂取率に変換する。静脈のR2 *値は、デオキシヘモグロビンの濃度、すなわち酸素摂取率に比例することが知られている(Xiang He他、 ” Quantitative BOLD: Mapping of Human Cerebral Deoxygenated Blood Volume and Oxygen Extraction Fraction: Default State”、 Magnetic Resonance in Medicine、2007年、57巻、115−126頁の式1と式2を参照)。ここで、静脈の横緩和速度をR2、ヘマトクリット(血液中に占める血球の体積の割合を示す数値)をHct、単位ヘマトクリットあたりの動静脈間の磁化率差をΔχdoとすると、酸素摂取率OEFは式(11)で表される。
酸素摂取率変換部333は、S1102で分離した静脈のR2 *分布の各画素値(R2 *値)を酸素摂取率に変換する。静脈のR2 *値は、デオキシヘモグロビンの濃度、すなわち酸素摂取率に比例することが知られている(Xiang He他、 ” Quantitative BOLD: Mapping of Human Cerebral Deoxygenated Blood Volume and Oxygen Extraction Fraction: Default State”、 Magnetic Resonance in Medicine、2007年、57巻、115−126頁の式1と式2を参照)。ここで、静脈の横緩和速度をR2、ヘマトクリット(血液中に占める血球の体積の割合を示す数値)をHct、単位ヘマトクリットあたりの動静脈間の磁化率差をΔχdoとすると、酸素摂取率OEFは式(11)で表される。
なお第一実施形態で説明した変形例、例えばモルフォロジーフィルタバンク以外の組織分離方法を採用することや、図4に示す手順の入れ替え、部分体積効果補正の省略等は、本実施形態においても同様に適用することができる。
R 2 が求められるパルスシーケンスとしては、例えば図10に示すように、反転パルス702−2や702−3を追加したMulti−echo SEシーケンス750を用いることができる。このシーケンスの詳細な説明は省略するが、このシーケンスでは、図5のGrEシーケンスと同様に、複数のTEの異なるエコー710、711が得られる。これらのエコーは静磁場の不均一の影響を含まないスピンエコーであり、図10のシーケンスを、スライス方向及び位相エンコード方向の各エンコード量を異ならせて繰り返すことにより複数のSE画像用データ(k空間データ)が得られる。画像再構成部320はこれら画像用データを再構成し、TEが異なる複数のSE画像を得る。
[R2 *分布(R’分布)算出:S1101]
第二実施形態のR2 *分布(R’分布)算出ステップS1101と同様に、物性値算出部(R2 *分布算出部)331はTEの異なる複数の画像の強度情報を用いて、画素毎に見かけの横緩和速度R2 *を算出し、R2 *分布を得る。計測において、GrE系パルスシーケンスのみを用いた場合にはR2 *分布が得られるが、第二実施形態の変形例のようにSE系パルスシーケンスを併用した場合にはR’(=R2 *−R 2 )分布が得られる。以下、R2 *分布とR’ 分布を代表して、R2 *分布が得られたものとして説明する。
第二実施形態のR2 *分布(R’分布)算出ステップS1101と同様に、物性値算出部(R2 *分布算出部)331はTEの異なる複数の画像の強度情報を用いて、画素毎に見かけの横緩和速度R2 *を算出し、R2 *分布を得る。計測において、GrE系パルスシーケンスのみを用いた場合にはR2 *分布が得られるが、第二実施形態の変形例のようにSE系パルスシーケンスを併用した場合にはR’(=R2 *−R 2 )分布が得られる。以下、R2 *分布とR’ 分布を代表して、R2 *分布が得られたものとして説明する。
Claims (1)
- 請求項4記載の磁気共鳴イメージング装置であって、
前記酸素量特性分布推定部は、
脳組織の酸素量特性の任意の画素の画素値とその周囲の画素の平均値が略等しいという条件と、
脳組織の酸素量特性の任意の画素の周囲の平均値と、静脈の酸素量特性の任意の画素の周囲の平均値が略等しいという条件とに基づいて前記対象領域の酸素量特性分布を推定することを特徴とする磁気共鳴イメージング装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016217307A JP6679467B2 (ja) | 2016-11-07 | 2016-11-07 | 磁気共鳴イメージング装置および酸素摂取率算出方法 |
US16/342,612 US10765349B2 (en) | 2016-11-07 | 2017-10-11 | Magnetic resonance imaging device and method for calculating oxygen extraction fractions |
PCT/JP2017/036826 WO2018083952A1 (ja) | 2016-11-07 | 2017-10-11 | 磁気共鳴イメージング装置および酸素摂取率算出方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016217307A JP6679467B2 (ja) | 2016-11-07 | 2016-11-07 | 磁気共鳴イメージング装置および酸素摂取率算出方法 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2018075075A JP2018075075A (ja) | 2018-05-17 |
JP2018075075A5 true JP2018075075A5 (ja) | 2019-06-27 |
JP6679467B2 JP6679467B2 (ja) | 2020-04-15 |
Family
ID=62076118
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016217307A Active JP6679467B2 (ja) | 2016-11-07 | 2016-11-07 | 磁気共鳴イメージング装置および酸素摂取率算出方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US10765349B2 (ja) |
JP (1) | JP6679467B2 (ja) |
WO (1) | WO2018083952A1 (ja) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10890641B2 (en) * | 2017-04-07 | 2021-01-12 | Cornell University | System and method of robust quantitative susceptibility mapping |
EP3301649B1 (en) * | 2017-09-07 | 2019-10-30 | Siemens Healthcare GmbH | Method for processing medical image data and image processing system for medical image data |
KR20190065941A (ko) * | 2017-12-04 | 2019-06-12 | 신에쓰 가가꾸 고교 가부시끼가이샤 | 탄화탄탈 피복 탄소 재료 및 그 제조 방법, 반도체 단결정 제조 장치용 부재 |
CN108829639B (zh) * | 2018-07-06 | 2023-10-27 | 上海联影医疗科技股份有限公司 | 一种磁共振成像方法和设备 |
JP7280997B2 (ja) * | 2018-08-29 | 2023-05-24 | 富士フイルムヘルスケア株式会社 | 画像処理装置及び磁気共鳴イメージング装置 |
JP7140606B2 (ja) * | 2018-08-29 | 2022-09-21 | 富士フイルムヘルスケア株式会社 | 画像処理装置、画像処理方法、画像処理プログラム及び磁気共鳴イメージング装置 |
CN111096748B (zh) * | 2019-12-19 | 2023-06-02 | 首都医科大学宣武医院 | 一种动态测量脑氧代谢率的方法 |
KR102097738B1 (ko) * | 2020-03-03 | 2020-04-06 | 주식회사 휴런 | 인공지능 기반 파킨슨 병 진단 장치 및 방법 |
CN113571188A (zh) * | 2021-08-02 | 2021-10-29 | 首都医科大学宣武医院 | 一种氧摄取分数的计算方法 |
CN115530820B (zh) * | 2022-11-30 | 2023-01-31 | 脑玺(苏州)智能科技有限公司 | 一种氧摄取分数测量方法、装置、设备及存储介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007149971A2 (en) * | 2006-06-22 | 2007-12-27 | General Hospital Corporation | Cancer detection by optical measurement of compression-induced transients |
US8958866B2 (en) * | 2008-03-31 | 2015-02-17 | The General Hospital Corporation | System and method to analyze blood parameters using magnetic resonance imaging |
CN103764025B (zh) * | 2011-10-12 | 2015-12-09 | 株式会社日立制作所 | 磁共振成像装置以及磁化率强调图像生成方法 |
JP6085545B2 (ja) * | 2013-09-26 | 2017-02-22 | 株式会社日立製作所 | 磁気共鳴イメージング装置、画像処理装置および磁化率画像算出方法 |
JP6533571B2 (ja) * | 2015-04-15 | 2019-06-19 | 株式会社日立製作所 | 磁気共鳴イメージング装置及び画像作成方法 |
-
2016
- 2016-11-07 JP JP2016217307A patent/JP6679467B2/ja active Active
-
2017
- 2017-10-11 US US16/342,612 patent/US10765349B2/en active Active
- 2017-10-11 WO PCT/JP2017/036826 patent/WO2018083952A1/ja active Application Filing
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2018075075A5 (ja) | ||
JP6679467B2 (ja) | 磁気共鳴イメージング装置および酸素摂取率算出方法 | |
JP5843876B2 (ja) | 磁気共鳴イメージング装置および磁化率強調画像生成方法 | |
JP6289664B2 (ja) | 磁気共鳴イメージング装置、定量的磁化率マッピング方法、計算機、磁化率分布計算方法、及び、磁化率分布計算プログラム | |
WO2014076808A1 (ja) | 磁気共鳴イメージング装置および定量的磁化率マッピング法 | |
RU2605524C2 (ru) | Магнитно-резонансная визуализация с подавлением артефактов потока | |
JP6085545B2 (ja) | 磁気共鳴イメージング装置、画像処理装置および磁化率画像算出方法 | |
Barral et al. | Diffusion-weighted MR imaging of the normal pancreas: reproducibility and variations of apparent diffusion coefficient measurement at 1.5-and 3.0-Tesla | |
JP6765972B2 (ja) | 拡散強調磁気共鳴イメージングを用いて組織内の鉄の量を定量化する方法、装置、及びコンピュータ・ソフトウェア | |
JP6533571B2 (ja) | 磁気共鳴イメージング装置及び画像作成方法 | |
Mutsaerts et al. | Gray matter contamination in arterial spin labeling white matter perfusion measurements in patients with dementia | |
CN111096748B (zh) | 一种动态测量脑氧代谢率的方法 | |
Biondetti et al. | Investigating the oxygenation of brain arteriovenous malformations using quantitative susceptibility mapping | |
Liu et al. | Quantification of liver iron concentration using the apparent susceptibility of hepatic vessels | |
US10324150B2 (en) | Method and apparatus for generating a magnetic resonance image data set | |
Opriessnig et al. | In vivo cardiovascular magnetic resonance of 2D vessel wall diffusion anisotropy in carotid arteries | |
US7310548B2 (en) | Method of magnetic resonance perfusion imaging | |
CN111195117A (zh) | 一种基于多延迟多参数非造影磁共振脑血灌注成像的方法 | |
Lindemeyer et al. | Iterative restoration of the fringe phase (REFRASE) for QSM | |
US11169237B2 (en) | Phase sensitive magnetic resonance angiography | |
Goldfarb et al. | Imaging of reperfused intramyocardial hemorrhage with cardiovascular magnetic resonance susceptibility weighted imaging (SWI) | |
US20230236272A1 (en) | Deep learning-based water-fat separation from dual-echo chemical shift encoded imaging | |
EP4246169A1 (en) | Augmenting diffusion-weighted magnetic resonance imaging data using a trained artificial neural network | |
Syeda et al. | A Continuum of T | |
US20150241526A1 (en) | Determining a measuring point-in-time in a cardiac cycle for conducting magnetic resonance diffusion measurements |