Susceptibility-Weighted Imaging:SWI

Slides:



Advertisements
Similar presentations
大学院物理システム工学専攻 2004 年度 固体材料物性第 8 回 -光と磁気の現象論 (3) - 佐藤勝昭ナノ未来科学研究拠点.
Advertisements

静脈画像を鍵とする暗号化手 法に関する研究 大山研究室 安藤のぞみ. 研究の背景、目的 近年、バイオメトリクス認証が注目されて いる 静脈は身体内部の情報 → 偽造に強い 環境に左右されることが少ない 利用者の心理的抵抗が軽減される オープンなネットワークへのバイオメトリ クス認証の適用 : Double.
生体情報を利用したオンライン認証システムに関する研 究 情報工学科 大山・山口・小尾研究室 学士課程4年田中 丈登.
三角関数演習問題 r b a [ 三角関数 ] θ 信号理論 (金田) 1演-1 (答は別紙の解答用紙に記入する)
コンソールの利用 零点・極と時間応答の関係 安定性 過渡応答の特性
MR based CFD による脳動脈瘤WSS分布推定における血管形状抽出法および 血液の非ニュートン性の影響
情253 「ディジタルシステム設計 」 (3)Constellation3
プロセス制御工学 3.伝達関数と過渡応答 京都大学  加納 学.
集積回路工学研究室 岩淵 勇樹 秋田 純一 北川 章夫
デジタル信号処理①
画像特徴(点、直線、領域)の検出と識別-2 呉海元@和歌山大学 2007年5月14日
5.アンテナの基礎 線状アンテナからの電波の放射 アンテナの諸定数
DARTs: Efficient scale-space extraction of DAISY keypoints
前回の内容 結晶工学特論 第4回目 格子欠陥 ミラー指数 3次元成長 積層欠陥 転位(刃状転位、らせん転位、バーガーズベクトル)
(質問)  体軸分解能を評価するための「SSPの測定」
透視投影(中心射影)とは  ○ 3次元空間上の点を2次元平面へ投影する方法の一つ  ○ 投影方法   1.投影中心を定義する   2.投影平面を定義する
ディジタル信号処理 Digital Signal Processing
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その1
ガウス誤差関数を利用した 収束の速いヒルベルト変換ディジタルフィルタ
羽佐田葉子 2007年3月24日 アクロス研究会@静岡大学
3次元での回転表示について.
SI = N(H) ・ (1 - e -TR/T1) ・ e -TE/T2 ・ e -bD 双極傾斜磁場 bipolar gradient
2.伝送線路の基礎 2.1 分布定数線路 2.1.1 伝送線路と分布定数線路 集中定数回路:fが低い場合に適用
MRI検査スピード化 (疎明に映すことと併用されれば、多数の救命に)
線形フィルタと畳み込み積分 マスクによる画像のフィルタリング 1.入力画像中の関心の画素のまわりの画素値
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その2
画像処理 基礎.
第32回MR基礎講座 (関西) 京都国際会館 画像法の原理(6) 拡散画像3 荏原病院放射線科 井田正博.
電気回路第1 第11回 ー電力の計算と演習ー 電気回路第1スライド11-1 目次(クリックすると移動します。) 2先週の復習 3電力の復習
et1 et1 et2 et2 信号 T2減衰曲線 Mxy(t) = M0 e-t/T2 T2*減衰曲線
横磁化成分と歳差運動 M0 横磁化Mxy 回転座標系 90°RFパルスにより、縦磁化成分Moはxy平面に倒れる(横磁化生成)
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その5
電力 P ( Power ) 単位 ワット W = J / sec
人工知能特論 9.パーセプトロン 北陸先端科学技術大学院大学 鶴岡 慶雅.
ミンコフスキー時空間 双曲筒 s2= x2+y2ーc2t2 双曲皿 s2= c2t2+y2ーx2 世界円筒 静止 s2 =x2+y2
位相カメラの進捗状況 京都大学修士1回 横山 洋海.
荏原病院放射線科・総合脳卒中センター 井田正博
5章  3次元形状を2次元面に投影する 3次元空間内に定義した形状を,2次元面上(ディスプレイのスクリーン面,プリンタの紙面など)に投影して表示するために必要になる変換について説明する.
構造情報に基づく特徴量を用いた グラフマッチングによる物体識別 情報工学科 藤吉研究室  EP02086 永橋知行.
ベクトル線図 周波数応答 G(jw) (– < w < ) を複素平面内に描いたものが、ベクトル線図である。
ディジタル信号処理 Digital Signal Processing
第32回MR基礎講座 (関西) 京都国際会館 画像法の原理(6) 拡散画像2 荏原病院放射線科 井田正博.
モデルに基づいた PID コントローラの設計 MBD とは モータ駆動系のモデリング モデルマッチング 5.1 節 出力を角速度とした場合
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その4
独立成分分析 (ICA:Independent Component Analysis )
入力付きシステムとラプラス変換 微分方程式がラプラス変換で解けるのなら、 システム (状態変数表現): の解もラプラス変換で求まるはず。
3次元での回転表示について.
東京農業大学 東京情報大学 附属第一高等学校・中等部 附属第二高等学校 附属第三高等学校・中等部
-画像処理(空間フィルタリング)- 画像処理(空間フィルタリング)のモデルとその基本操作 雑音除去・平滑化への適用
臨床診断総論 画像診断(3) 磁気共鳴画像 Magnetic Resonance Imaging: MRI その3
第 6 章 :フィードバック制御系の安定性 6.1 フィードバック系の内部安定性 6.2 ナイキストの安定定理
SWI 強度画像、位相画像、SWI 強度画像 位相画像 S WI.
電磁気学Ⅱ Electromagnetics Ⅱ 6/9講義分 電磁場の波動方程式 山田 博仁.
Mini-RT装置における 強磁場側からの異常波入射による 電子バーンシュタイン波の励起実験
Interferometry and Synthesis in Radio Astronomy 3.1章
平面波 ・・・ 平面状に一様な電磁界が一群となって伝搬する波
Fourier 変換 Mellin変換 演習課題
9. ナイキスト線図と安定余裕 教科書 7.2, 7.3.
偏光X線の発生過程と その検出法 2004年7月28日 コロキウム 小野健一.
4. システムの安定性.
8方向補間ブロックマッチングの実装 福永研究室 数理科学コース 学部4年 能城 真幸.
誘導起電力は 巻数と 磁束の時間変化 に比例する.
1ーQー18 音声特徴量抽出のための音素部分空間統合法の検討
第 5 章 :周波数応答 5.1 周波数応答と伝達関数 周波数伝達関数,ゲイン,位相 キーワード : 5.2 ベクトル軌跡 ベクトル軌跡
電磁気学C Electromagnetics C 5/20講義分 電磁場の波動方程式 山田 博仁.
ソースフィルタモデル.
情報通信システム(2) plala. or 情報通信システム(2) 年4月23日 火曜日  午後4時10分~5時40分 NTT-TX Corp. 加藤 洋一.
電磁気学C Electromagnetics C 7/10講義分 電気双極子による電磁波の放射 山田 博仁.
電気回路学Ⅱ 通信工学コース 5セメ 山田 博仁.
Fourier 変換 Mellin変換 演習課題
Presentation transcript:

Susceptibility-Weighted Imaging:SWI 3D FLASH 高い空間分解能 0.8x0.7x1.6 mm voxel内のdephasingを最小限 大きな位相分散による信号低下を避け、位相差情報を検出するため Rauscher A. Magn Reson Med. 2005 Jul;54(1):87-95.

SWIの理解に虚数が必要? Nombre imaginaire M sinf 虚数 imaginary number 2乗してマイナスになる数 虚数単位:i 2乗して-1になる数 i =√ -1 オイラーの等式 e ip + 1 = 0 横磁化 M f 実成分 Re M cos f 図2 実成分と虚成分,強度画像と位相画像 MR信号は位相差による実成分(real、同位相成分)と虚成分(imaginary、直交位相成分) からなる.実成分と虚成分は直交配列するquadrature型コイルを考えると理解しやすく,直交座標系として表現することができる . Re: 実成分 Im: 虚成分 M: 横磁化成分の絶対値 f : 位相角 Xy平面上の横磁化成分もベクトルとして考えることができる。

複素平面(ガウス平面、アルガン平面) M M f f 複素数=複数の要素をもった数 ベクトル平面=複素平面 虚数の数直線(虚数軸) 虚成分 Im M sinf 5+4i 横磁化 4i M M f f - 5 5 実成分 Re M cos f 実数の数直線(実数軸) Re: 実成分 Im: 虚成分 M: 複素数の絶対値 f : 虚数の偏角 Re: 実成分 Im: 虚成分 M: 横磁化成分の絶対値 f : 位相角 - 4i 複素数=複数の要素をもった数 ベクトル平面=複素平面 Xy平面上の横磁化成分もベクトルとして考えることができる。

実チャネルと虚チャネル、実成分と虚成分 M 実数成分(cos)→同位相 虚数成分(sin)→直角位相 f 虚チャネル 虚成分Im 横磁化 M sinf 横磁化 M 実数成分(cos)→同位相 虚数成分(sin)→直角位相 実チャネル f 実成分 Re M cos f 信号強度と方向 MR信号は位相差による実成分(real、同位相成分)と虚成分(imaginary、直交位相成分) からなる.実成分と虚成分は直交配列するquadrature型コイルを考えると理解しやすく,直交座標系として表現することができる . Re: 実成分 Im: 虚成分 f : 位相角

実成分と虚成分,強度画像と位相画像 強度画像 magnitude image M = √Re2 + Im2 位相画像 phase image tan f = Im / Re f = arctan [Im / Re] (-p/2<f<p/2) 虚成分Im M sinf 横磁化(MR信号)  M f 実成分 Re M cos f Re: 実成分 Im: 虚成分 M: 強度画像 f : 位相角

実成分と虚成分,強度画像と位相画像 実数部分k-space 虚数部分k-space FT M 実数画像 虚数画像 f 虚成分Im M sinf MR信号 FT M 実数画像 虚数画像 f 実成分 Re M cos f 強度画像magnitude 位相画像phase Re: 実成分 Im: 虚成分 M: 強度画像 f : 位相角 SWIは強度画像と位相画像を用いる

Postprocessing of SWI 撮像データからSWIまでの後処理,再構成までのステップ High-pass filter 磁化率差以外の要因による背景の位相変化を除去 背景に含まれる低周波数成分の位相の乱れを取り除く。 “Phase” mask 位相差を強調するために、信号を落としたい位相の部分が0になるようなマスク画像を作成 位相画像上で位相が進んでいる(あるいは遅れている部分)を0、位相差のない部分は1にしたような画像 1.High-pass filter realigned data 2.Phase mask image Phase image 3.multiplication 4.mIP Magnitude image SWI Rauscher A.Magn Reson Med. 2005 Jul;54(1):87-95.

位相マスク画像の掛け合わせ 位相画像 ⇒ マスク ⇒ SWI p < f < 0 0 < f masked < 1 r = r magnitude*[ fmasked ]n  位相差  - p/2 と - p/4  位相マスク画像を4回掛け合わせるとすると → r = r magnitude*[ fmasked ]4 1 = - p/2 f1 masked = 0.50 r = r1 magnitude*0.0625 2 = - p/4 f2 masked = 0.25 r = r2 magnitude*0.0039 2:1 16:1 SWIでは複数回掛け合わせ(n)により、磁化率変化による位相差を強調 T2*WIに比較して10倍以上の鋭敏さがある.

SWIは磁化率に鋭敏 撮像voxelの1/4程度の大きさまで磁化率効果が描出される 静脈 撮像voxelの1/4程度の大きさまで磁化率効果が描出される 0.50mmのvoxelでは0.25mm径の静脈が描出できる 大きさを過大評価 Partial volume effect 静脈 T2WI SWI

Postprocessing of SWI 撮像データからSWIまでの後処理,再構成までのステップ “Phase” mask multiplication 位相マスク画像を強度画像(magnitude image)に複数回掛け合わせる(multiplication)ことで、画像のCNRを向上させる Minimum intensity projection; mIP 静脈を評価する場合には、ある程度のスライス厚情報をもって画像化すると有用なため,最小値投影法mIPによる再構成画像を作成 1.High-pass filter realigned data 2.Phase mask image Phase image 3.multiplication 4.mIP Magnitude image SWI Rauscher A. Magn Reson Med. 2005 Jul;54(1):87-95.