MRI(Magnetic Resonance Imaging)の高速化

Slides:



Advertisements
Similar presentations
地図の重ね合わせに伴う 位相関係の矛盾訂正手法 萬上 裕 † 阿部光敏* 高倉弘喜 † 上林彌彦 ‡ 京都大学工学研究科 † 京都大学工学部 * 京都大学情報学研究科 ‡
Advertisements

て -form - Making て -form from ます -form -. With て -form, You can say... ~てもいいですか? (= May I do…) ~てください。 (= Please do…) ~ています。 (= am/is/are doing…) Connecting.
第 5 章 2 次元モデル Chapter 5 2-dimensional model. Contents 1.2 次元モデル 2-dimensional model 2. 弱形式 Weak form 3.FEM 近似 FEM approximation 4. まとめ Summary.
Essay writing rules for Japanese!!. * First ・ There are two directions you can write. ・よこがき / 横書き (same as we write English) ・たてがき / 縦書き (from right to.
高度情報演習1A “テーマC” 実践 画像処理プログラミング 〜画像認識とCGによる画像生成〜 第四回 演習課題 画像中からの物体抽出処理(背景情報を手がかりとして) 芝浦工業大学 工学部 情報工学科 青木 義満 2006/05/15.
医用画像機器工学実習 断層画像CT( Computed Tomography) を得る方法 フィルタ重畳逆投影法
MR based CFD による脳動脈瘤WSS分布推定における血管形状抽出法および 血液の非ニュートン性の影響
PLの実験内容(2002~2003) He-Ne 一様励起のみ 魚の地図を作成 Ti-Sa 魚、船、原子的平面の温度依存性 (4K-60K)
ウェーブレットによる 信号処理と画像処理 宮崎大輔 2004年11月24日(水) PBVセミナー.
CSWパラレルイベント報告 ヒューマンライツ・ナウ        後藤 弘子.
2010年7月9日 統計数理研究所 オープンハウス 確率モデル推定パラメータ値を用いた市場木材価格の期間構造変化の探求 Searching for Structural Change in Market-Based Log Price with Regard to the Estimated Parameters.
第32回MR基礎講座 (関西) 京都国際会館 画像法の原理(6) 拡散画像1 荏原病院放射線科 井田正博.
多数の疑似システムを用いた システム同定の統計力学 三好 誠司 岡田 真人 神 戸 高 専 東 大, 理 研
MEG実験2009 陽電子スペクトロメータの性能評価
画像特徴(点、直線、領域)の検出と識別-2 呉海元@和歌山大学 2007年5月14日
SHIGENOBU HIROSE AND JULIAN H. KROLIK
SP0 check.
わかりやすいマルチスライスCTにおける画像再構成
荏原病院放射線科 総合脳卒中センター 井田正博
医療支援診断のためのコンピュータ分散システムの検討
Irradiated Polarized Target
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その1
画像工学 2011年10月6日 担当教員 北川 輝彦.
デジタル信号処理④
Group meeting 2016/5/13 Katsuhiro Umeda.
The Sacred Deer of 奈良(なら)
Who Is Ready to Survive the Next Big Earthquake?
生物科学科(高分子機能学) 生体高分子解析学講座(第3) スタッフ 教授 新田勝利 助教授 出村誠 助手 相沢智康
脳血管 MR診断に必要な脳動脈の解剖 荏原病院放射線科 井田正博
みかけの拡散係数 Apparent Diffusion Coefficient:ADC
SI = N(H) ・ (1 - e -TR/T1) ・ e -TE/T2 ・ e -bD 双極傾斜磁場 bipolar gradient
MRI検査スピード化 (疎明に映すことと併用されれば、多数の救命に)
Kazuhide Okada, Ph.D. Professor
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その2
画像工学 2012年10月3日 担当教員 北川 輝彦.
Multi-Pixel Photon Counter(MPPC)の開発
独立成分分析 1.問題は何か:例:解法:全体の見通し 2007/10/17 名雪 勲.
第8週 高精度GPSの構築 位相測位の原理 通信システムの構築.
Bottom-UpとTop-Down アプローチの統合による 単眼画像からの人体3次元姿勢推定
第32回MR基礎講座 (関西) 京都国際会館 画像法の原理(6) 拡散画像3 荏原病院放射線科 井田正博.
正規分布における ベーテ近似の解析解と数値解 東京工業大学総合理工学研究科 知能システム科学専攻 渡辺研究室    西山 悠, 渡辺澄夫.
決定木とランダムフォレスト 和田 俊和.
et1 et1 et2 et2 信号 T2減衰曲線 Mxy(t) = M0 e-t/T2 T2*減衰曲線
横磁化成分と歳差運動 M0 横磁化Mxy 回転座標系 90°RFパルスにより、縦磁化成分Moはxy平面に倒れる(横磁化生成)
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その7
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その5
混合ガウスモデルによる回帰分析および 逆解析 Gaussian Mixture Regression GMR
Ⅱ 磁気共鳴の基礎 1.磁場中での磁気モーメントの運動 2.磁気共鳴、スピンエコー 3.超微細相互作用、内部磁場 references:
◆3.0T MRIを導入 ◆世界初、画像を信号を全てデジタル化 3.0T MRI検査開始のお知らせ 大生病院のMRI装置が新しくなりました。
荏原病院放射線科・総合脳卒中センター 井田正博
Susceptibility-Weighted Imaging:SWI
第32回MR基礎講座 (関西) 京都国際会館 画像法の原理(6) 拡散画像2 荏原病院放射線科 井田正博.
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その4
予測に用いる数学 2004/05/07 ide.
QRコードを用いたIDカードに 適した電子透かし
従来CTとの検査比較( 頭部検査) 1 Slice 16 Slices 1 Slice 16 Slices ルーチン検査 撮影範囲
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その3
2.4 Continuum transitions Inelastic processes
SWI 強度画像、位相画像、SWI 強度画像 位相画像 S WI.
Dark Matter Search with μTPC(powerd by μPIC)
Bottom-UpとTop-Down アプローチの組み合わせによる 単眼画像からの人体3次元姿勢推定
北大MMCセミナー 第62回 附属社会創造数学センター主催 Date: 2016年11月4日(金) 16:30~18:00
川崎浩司:沿岸域工学,コロナ社 第4章(pp.58-68)
MO装置開発 Core part of RTR-MOI Photograph of core part.
ポッツスピン型隠れ変数による画像領域分割
実験結果速報 目的 装置性能の向上 RF入射実験結果 可動リミター挿入 RFパワー依存性 トロイダル磁場依存性 密度依存性
間欠型一酸化炭素中毒に対する高気圧酸素治療の限界
物理学実験 II ブラウン運動 ー 第2日目 ー 電気力学結合系の特性評価 物理学実験II (ブラウン運動) 説明資料.
Apply sound transmission to soundproofing
アノテーションガイドラインの管理を行う アノテーションシステムの提案
Presentation transcript:

MRI(Magnetic Resonance Imaging)の高速化 平成28年4月16日 (土) 岡田一秀

紹介 脳血管障害の種類 - 脳梗塞の医療統計 - 診断工程 - CT・MRI双方のメリット ‐ CTへのアルゴリズムは?   紹介 脳血管障害の種類 -  脳梗塞の医療統計 -  診断工程 -  CT・MRI双方のメリット ‐  CTへのアルゴリズムは? ‐  MRIでのデータ取得へのメカニズムは? -  MRIでの、電波を切ってからの横緩和・縦緩和について -  MRIで、高速診断を目指すためのアプローチ [Kolonia, Pohnpei, FSM] 九州の方々には、今回の大地震に関して、南太平洋より、お見舞い申し上げます。 http://twitcasting.tv/kehua_and_kazu http://twitcasting.tv/kehua_and_kazu/show/ (過去分) 医工学について、毎月、第1週と第3週は日本語で、第2週と4週は英語で、病種や 医療機械種を変え、日本時間の午後6時(金)か午前9時(土)頃、やっております。 お時間があれば、どうかお聴き下さい。(録音もしておきます)  ー ポンペイ コロニア第2スタジオより 

Disability-adjusted life year for cerebral vascular disease per 100,000 inhabitants in 2004 (unit:×10 PPM) < Annual Statistics > ・Estimated 15 million people  worldwide survive after disease onset ・In US 4.8 millions survive. ・In Japan 1.5 millions survives and 0.13 millions dies

Stroke Thrombus 3 membranes which cover the brain - Ischemic (Infarction) - 60% by interruption of the blood supply to the brain * Lacunar infarction (Φ<1.5cm) * Atherothrombotic infarct * Cardiogenic embolism - Hemorrhagic (Bleeding) – 40% by the rupture of a blood vessel or an abnormal vascular structure * Epidural hematoma (bleeding between the dura mater and the skull) * Subdural hematoma (bleeding in the subdural space) * Subarachnoid hemorrhage <くも膜下出血> (bleeding between the arachnoid mater and pia mater) * Brain hemorrhage      <脳内出血>

Desired Present Diagnosis 1) Brain check by CT - in order to check whether the affected area is bleeding or blocked on the blood vessel 2) Brain check by MRI - Precise Detection of even the small infarctions or the slight bleeding 3) Electrocardiogram / Heart Echo - In order to judge the said symptom is caused by Cardiogenic embolism (The thrombus flew from the heart to the brain) or not 4) Carotid (Artery around the neck) ultrasonography / Vascular imaging - To observe the status of the hardening/blockade/stenosis of the arteries

得意分野 同じ、脳卒中でも。。。

Diagnosis both by CT and MRI in case of Infarction Slight Early sign by CT ⇒ Unclear part on the nucleus lentiformis (the frustule + the pallidum) inside the basal ganglia appears in 1-2 hours after the disease onset. Unclear part on the boundary between the white substance and the cortex appears in 2-3 hours after disease onset. The phenomenon by which the cerebral ditch becomes narrower and smaller over 3 hours after disease onset. Follow by MRI ⇒ On DWI (diffusion weighted image) mode the affected area is reflected white because the blocked vessel doesn’t have the stream of the Hydrogen molecules.

Indirect Cause (High Risk) - Blood pressure - Blood lipids (High cholesterol levels) - Diabetes mellitus - Anticoagulation drugs (The medicine which obstructs coagulation of blood) - Lack of Exercise Amyloid as the body wastes’ type protein is accumulated inside the blood vessel in the cerebral cortex of the brain(especially for aged persons) Ischemic detected by CT (Cerebral right hemisphere) Hemorrhagic detected by CT (Cerebral right hemisphere)

脳血管疾患の現状

MRI (核磁気共鳴画像法(nuclear magnetic resonance imaging, NMR) 人の身体の2/3は水で構成されており、 更に、各種脂肪酸やアミノ酸などの 構造式にも水素(H)が含まれているために、 現行の医療用MRIでは、この水素原子(1H)の 信号で画像が描出されている

問題 体内の様子を画像化するために必要な情報を丁寧に取得するため に、比較的長い時間を要するために、例えば子供の場合、その計測 時間の間おとなしく待つことができないという問題がある。 しかし仮にこのMRIによる体内の画像取得のための計測時間を少な くすることができるのであれば、このような問題はなくなる。またどん な被験者に対しても気軽に診察のためにMRIを使うことができるメ リットも大きいため、喉から手が出るほとに嘱望された技術 立体物を透かしてみることができるのであれば、それこそ前から上から横から眺めて、中身はどうなっているのだろうかと推定すれば3次元ピクロスとなる。そうやってモノの中身を探るピクロスのような技術、それが圧縮センシングである。 MRAはMRI画像の中でもスパース性が高いことから圧縮センシングとの相性が良い 世界各国で圧縮センシングの利用によるMRIの計測手法の改善が試行錯誤されて、MRI機器の中に搭載する動き

画像の再構成 -CTで ・ 単純逆投影法 ・ フィルタ補正逆投影(FBP)法 画像の再構成  -CTで 単純逆投影の結果(左)を観察すると,単純逆投影分布は,元の2次元分布よりも広がりをもった値になっている。単純逆投影分布b(x,y)とX線減衰係数の2次元分布f(x,y)は,次式のように点拡がり関数h(x,y)によって関連付けられる。ここで*は畳み込み演算である。b(x,y)=f(x,y)*h(x,y)。さらに,h(x,y)は次式で与えられる。h(x,y)=√ (x2+y2) 。したがって,単純逆投影分布b(x,y)に対してh(x,y)をキャンセルするようなフィルタリングを行う(計算処理負担の概念からは、2次元フィルタを通すのではなく、1次元フィルタを用いる)ことで,X線減衰係数の2次元分布f(x,y)を求める ・ 単純逆投影法              ・ フィルタ補正逆投影(FBP)法 (信州大学・小関先生より)

データのInからOutまで

パルスエコー法 最初に被検体に90°パルスを印加して磁化を静磁場方向と垂直な平面に 倒す(横磁化)。 スピンエコー法⇒プロトン励起に90°パルスから1/2TE時間後に180° パルスを印加し、1/2TE時間後、つまりTE時間後にMR信号 (スピエコー信号)を得るシーケンス。 90° パルス→励起パルス 180°パルス→再収束パルス プロトン励起時に、撮像断面を決めるためのスライス選択用傾斜磁場を印 加する。傾斜磁場を印加することで選択する断面の磁場強度を共鳴周波数 に合わせて、その他の部位をプロトンの共鳴周波数外に位置させ、特定の スライス断面を選択的に励起する。スライス断面選択の操作は、他のパル ス系列でも同様に行われ、プロトン励起時にスライス選択用傾斜磁場印加 を行う。 この後に行われる位相及び周波数エンコード操作はパルス系列の違いで大 きく異なり、撮像時間に影響を与える。

MRIの撮像法 位相エンコード方向に直交する周波数エンコード方向に傾斜磁場を印加する。(信号を読取る間に印加される傾斜磁場:読取り勾配磁場)。⇒磁場強度に比例したプロトンの共鳴周波数変化を生じ周波数   エンコード方向の位置情報が収集可。 ⇒周波数エンコードは繰り返し行う必要はなく k-space 1行のデータをすべて収集できる。 (実際には、スライス選択の後、180°RFパルスを印加する前に、一度   frequency encode gradientが原因となり周波数エンコード方向で生じる   位相の乱れ(dephase)をあらかじめ補正する目的で使用される。)  

撮像時間 いかなる種類のgradientでも共鳴プロトンの位相の乱れを生じさせ、 画像劣化の原因となためこれを可能な限り補正する必要がある。 180°RFパルスの印加前に付加されるため180°RFパルス後、極性が 反転し、負の磁気勾配と同様の作用を示す(±を入れ替える)共鳴周 波数の収集時する正の磁場勾配である。 SE法の特徴は実在する様々な磁場の不均一性が180°パルスにより 補正されるため、不均一磁場の影響を受けたT2*強調像ではないT2 強調像が得られる点である。 撮像時間 = (TR)×(位相エンコード数)×(加算回数)

画像種 (i) T1強調画像 繰り返し時間(TR)を短くし、エコー時間 (TE)を最小に設定したシーケンスでは、 画像のコントラストはT1に依存するので、 T1強調画像である。T1WIでは水は黒く 低信号で描出 され(脳室は黒色)、CTと よく似た画 像を呈し、大脳皮質と白質等 の解剖 学的な構造が捉えやすい特徴 (ii) T2強調画像 TRおよびTEを長く設定したシーケンスでは、画像のコントラストはT2に依存する ので、T2強調画像である。T2WIでは水は白く高信号で描出 され(脳室は白色) 多くの病巣が高 信号で描出されるため、病変の抽出 に有用とされている。 (iii) プロトン密度強調画像 TRを長く、TEを短く設定したシーケンスでは、プロトン密度強調画像である。水分 子の拡散運動(自由運動度) を画像化したもの。拡散が低下した領 域が高信号 として描出される。急性期 の脳梗塞では、拡散が低下して来るため、 超急性期の脳梗塞の部位判定 (白色に描出される)に有用

By The fast processing of MRI MRI’s basic mechanism Non-invasive Tomography Procession of the equinoxes ω0 ω0 : Larmor Angular Frequency Magnetic flux

Transition of Procession of the equinoxes on MRI Stage Ⅰ  Ⅱ Ⅲ Ⅳ External Magnetic OFF ON ON ON Field Radio Wave OFF OFF ON OFF Unify Magnetized Direction Fall it down and unify Phase Proton is rotated, being positively charged  The relationship between its Flat surface for spinning and the inner magnetic flux is orthogonal

Evaluation of axial Relaxation on MRI

The utilization of CS (Compressed Sensing) - as one of Informatics’ Strategy In accordance with the time-constant of transient on Longitudinal or traverse relaxation, Getting ≒Lamore Frequency and then doing IFFT for its value, we can get Image by which we can distinguish between fat, water, and the hard organ.  We can reduce number of yi Forecasting x1 to xn change Condition  Initial Data is Sparse xi : brightness ( 1≦i≦n ) yj : Frequency ( 1≦j≦m )

なぜ、いま、圧縮センシング? 胸部CT撮影における放射線の被曝量は6.9mSvであり,他のX線診断 装置と比較して被曝量が多い。これらの問題を解決するには,360° の全周囲方向から撮影する必要をなくすであるとか,1000回以上も必 要な撮影回数を100回以下にするなどの画期的な技術開発が必要で ある。2005年頃まで、少ない方向数で撮影した画像から十分に高精 度な断面画像を得ることは困難であると考えられていた2006年に,標 本化定理の仮定のもとでは不十分なサンプル数の測定データから原 信号を完全に復元することができる圧縮センシングという手法が発見 された。圧 縮 セ ン シ ン グ の 手 法 をCT再構成に応用することで,従 来よりもはるかに少ない投影角度数の投影データからでも断面画像 の構成が可能であることが多数報告されている