Graduation thesis Diffusion model and LVAD

卒業研究でのプレゼン資料

1. ⾳響解析による左室補助⼈⼯⼼臓の異常検知 知能データ科学講座 櫻井研究室 岡本卓也 2024.2.13
2. 1. 研究の背景・⽬的 (1/3) 【社会的な背景】 重症⼼不全に対する植込み型左室補助⼈⼯⼼臓(LVAD)治療 ・⼼臓移植までの橋渡しの役割 ・Destination Therapy(DT)への適応拡⼤ ⽶国ではDTでの植え込みLVAD患者数が増加 ⽇本でも2021年に保険適応 LVADの⻑期在宅管理の必要性 © 2024 Sakurai & Matsubara Lab. @ Osaka U DT症例数 (Miller LW, NEJM 2007) (Pinney SP, JACC 2017)
3. 2. 研究の背景・⽬的 (2/3) 【⽬的】 LVADに起因する合併症をLVADの駆動⾳を⽤いて検知したい。 合併症の例︓ © 2024 Sakurai & Matsubara Lab. @ Osaka U ポンプ⾎栓 (Scandroglio AM, JACC 2016 ⼤動脈弁逆流 (Grinstein J, JACC 2016)
4. 3. 研究の背景・⽬的 (3/3) 【⽬的】 LVADに起因する合併症をLVADの駆動⾳を⽤いて検知したい。 ⽣活⾳が⼊る在宅環境下でこれらの異常を検知したい →LVAD駆動⾳を⼊⼒とした拡散確率モデルを⽤いた分類を提案する © 2024 Sakurai & Matsubara Lab. @ Osaka U LVAD駆動音 人の声 家電
5. 4. 前提知識 拡散確率モデルとは ・深層⽣成モデルの⼀種 ・データにノイズを徐々に加えていく拡散過程を逆向きに辿る逆拡散過程によって、 ⽣成モデルを定義する。 ・多くのタスクで⾼性能なデータを⽣成可能(e.g.画像、動画、⾳声) (Ho et al.2020) © 2024 Sakurai & Matsubara Lab. @ Osaka U
6. 5. 仮説 © 2024 Sakurai & Matsubara Lab. @ Osaka U 6 【仮説】 正常クラスの波形データだけを⼊⼒として拡散モデルを学習する。 分類したい波形データにノイズを載せる。 その波形を学習したモデルを⽤いたノイズ除去の是⾮によって 正常か異常か分類できるのではないか。
7. 6. 類似研究との差異 【電⼒消費での検知】 電⼒消費をモニタリングし，ポンプ⾎栓の検出 →重篤な症状に発展する軽微な症状の、早期検出に⼗分な感度・特異度ではない 電⼒消費のみでの検知には限界あり (Jorde UP. JACC.2015) (Grabska J. ASAIO J.2020) 【⾳響データでの検知】 LVAD駆動⾳を対象に、機械学習による⼤動脈弁逆流の検出 →良い性能の分類モデルが得られた ポンプ⾎栓を同様に分類した結果、分類性能は⼗分でなかった 特徴ベクトルを正しく選択しないと性能が出ない VAD駆動⾳と無関係なノイズが発⽣. (Yusuke Misumi. J Artif Organs. 2021) © 2024 Sakurai & Matsubara Lab. @ Osaka U 周波数 (Hz) 時間 (0.05sec) 強さ 正常時 周波数 (Hz) 時間 (0.05sec) 強さ ポンプ内⾎栓 あり VAD駆動⾳と無 関係なノイズが 出現
8. 7. 問題定義 Given 電⼦聴診器で録⾳した実際のLVAD駆動⾳(44.1kHzサンプリング、25.3 sec≦ len ≦ 79.5 sec ) ・正常クラスの波形データ108個 ・異常クラスの波形データ188個 Output 波形データのクラス分類 © 2024 Sakurai & Matsubara Lab. @ Osaka U
9. 8. 提案⼿法 1.拡散確率モデル(DDPM)の学習 正常クラスの波形データのみを⼊⼒とし、DDPMベースのニューラルモデルの学習 ・本研究では波形データを⼊⼒とするDiffWaveモデルを採⽤(Z.Kong et al.2021) ・学習⽤・テスト⽤に正常の波形データを8:2に分割 2.波形データのノイズ除去 テスト⽤波形データにノイズを載せ、モデルに⼊⼒しノイズ除去を⾏う。 ・正常テストデータと、同じ数の異常テストデータを⽤いる。 3.出⼒波形の評価 出⼒波形と、ノイズを載せる前の元の波形との絶対平均誤差をスコアとして分類を⾏う。 ⽐較⼿法として、Local Outliner Factor(LOF) と Isolation Forest(iForest)を採⽤した。 © 2024 Sakurai & Matsubara Lab. @ Osaka U
10. 9. 実験結果 ・分類結果 LOFやiForestは極端な分類結果 提案⼿法では他の⼿法より⾼い精度 ・ROC曲線とAUC ROC曲線が左上に寄るほど AUCが⾼いほど良い分類 提案⼿法は、ランダムを含め他の⼿法よりも性能が良い →他の⼿法と⽐較して優位性がみられた。 © 2024 Sakurai & Matsubara Lab. @ Osaka U
11. (a)正常クラス (b)異常クラス 同じ程度のスコアとなり、分類できなかった正常クラスと異常クラスの 波形データの周波数成分を⽐較する。 →異常クラス波形では低周波数帯で、 振幅の強い⾳が帯状に観測された。 周波数成分も学習に⽤いることで分類性能が向上すると予想 10. 考察 © 2024 Sakurai & Matsubara Lab. @ Osaka U
12. 11. まとめ • LVAD駆動⾳を⼊⼒とする拡散確率モデルを⽤いた 異常検知⼿法の提案 • 実際の医療データを⽤いた実験による提案⼿法の評価 • 今後の課題 周波数成分を⽤いたモデルの学習・分類 分類スコアの⼯夫 © 2024 Sakurai & Matsubara Lab. @ Osaka U
13. DiffWave 構造 Input Conv1×1 = Conv1×1 Bi-DilConv-2$ = Bi-directional Dilated Conv (dilation = 2$) FC = Fully connected = Broadcast over length Conditioner ↔ = Element-wise addition + = Element-wise multiplication Diffusion-step embedding Residual layer ' = 0 Residual layer ' = 1 Residual layer ' = * − 1 ⋯ Skip connections Output ⋯ FC FC FC Conv1×1 Conv1×1 Conv1×1 Conv1×1 Conv1×1 Bi-DilConv-2- ./0 1 ↔ + + + + ⋅ ⋅ tanh 7 = Connect to next residual layer = Input of each residual layer Conv1×1 ReLU swish ReLU swish © 2023 Sakurai & Matsubara Lab. @ Osaka U 13
14. 詳細な式 © 2023 Sakurai & Matsubara Lab. @ Osaka U 14 <latexit sha1_base64="Bd3UUUFuC1QzBTek4tRT9ByB8Io=">AAADI3ichZFLb9QwEMcn4VWWRxe4IHGJWIpawa6cigJCQqrggji13W5bqS6Rkzobq86jiXdFa/kL8AU4cAKpQojvwIULEjcEhx74AIhjkbhwYJINj1IBE8Uez/x/47HtZ1IUipBdyz50+MjRY2PHGydOnjo93jxzdqlIB3nAe0Eq03zFZwWXIuE9JZTkK1nOWexLvuxv3C3zy0OeFyJNFtVWxtdi1k9EKAKmMOQ1X9GYqcgP9UPjadV2zW2HhjkLtGs0LTZzpSmTWcQ8ZYxDJQ+VM/kbopx2rac+V6Wqptw29Vlew6gzDuI8K4TEXTVVEarNvkpXHTVFc9GP1JRzxaFKyHU+qorJB/pnV9PG/MC2jddskQ6pzDnouLXTgtrm0uZzoLAOKQQwgBg4JKDQl8CgwG8VXCCQYWwNNMZy9ESV52CggewAVRwVDKMbOPZxtVpHE1yXNYuKDnAXiX+OpAMT5AN5QfbIG/KSfCLf/lpLVzXKXrZw9kcsz7zxR+e7X/9LxTgriH5R/+xZQQg3q14F9p5VkfIUwYgfbj/e695amNCXyTPyGft/SnbJazxBMvwS7MzzhSfQwAdw/7zug87SdMe93pmZv9aavVM/xRhcgIswifd9A2bhHsxBDwLrknXf6lqL9o791n5nvx9JbatmzsE+sz9+B3C81c0=</latexit> xt 1 = 1 p ↵t ✓ xt t p 1 ¯ ↵t ✏✓(xt, t) ◆ + ˜1 2 t z <latexit sha1_base64="df/tZCmQx/qZq6OOi/iQvsQmhnY=">AAAC4HichVHPS9xAFH5JrT/WqrG9FLyELpYV7PJW/FEKBVGEHjz4o6uCsSFJZ3cHk0lIZre1MXfxUnoS6klBivhneOk/0IPHHsWjhV568G02alVqX8jM975535tvZuzA5ZFEPFHUB20P2zs6u3Ldj3p6+7T+x0uRXw8dVnZ81w9XbCtiLhesLLl02UoQMsuzXbZsr08315cbLIy4L97KjYCteVZV8Ap3LEmUqcWGx4UZG7LGpJXMmnFdfGC8WpNJocUN6a8Nz5I1245nEipMcSX+mJg4bLAg4q4vkpyxqV8m+otLdNW28JdKDsshXdeNzXcjppbHIqah3wWlDOQhizlf+wYGvAcfHKiDBwwESMIuWBDRtwolQAiIW4OYuJAQT9cZJJAjbZ2qGFVYxK7TWKVsNWMF5c2eUap2aBeX/pCUOgziDzzEc/yOR3iKf/7ZK057NL1s0Gy3tCww+7afLv7+r8qjWULtWnWvZwkVeJl65eQ9SJnmKZyWvvFp53zx1cJg/Bz38Yz87+EJHtMJROOXczDPFnYhRw9Qun3dd8HSSLE0XhybH81PTmVP0QkD8AwKdN8TMAlvYA7KtO9PpUPRlH7VVrfUz+qXVqmqZJoncCPUrxccrblV</latexit> min ✓ Lunweight(✓) = E x0,✏ k✏ ✏✓(xt, t)k2
15. 実験結果 混合⾏列 © 2023 Sakurai & Matsubara Lab. @ Osaka U 15