Tech #5 — WIP/C2P Perfusion Quantification Dual Sequence 徹底深掘り

qPerfusion = Dual Sequence + Fermi/BTEX deconv + MoCo + Gd Map + MBF/MPR

心筋灌流の **絶対値定量化 (Quantitative Perfusion)** を実現する WIP/C2P sequence の徹底解説。Dual Sequence (LV cavity AIF + myocardium signal を同時取得) → Gd Concentration Map 変換 → Fermi or BTEX deconvolutionMBF (Myocardial Blood Flow) / MPR (Myocardial Perfusion Reserve) 算出。KCL Amedeo Chiribiri + NIH Peter Kellman 双璧の論文徹底解説 + 三重大/東北大 routine。Common Limitation は Common Limit page 参照。

TL;DR

  1. Dual Sequence = LV cavity (AIF) + myocardium 同時取得。AIF saturation 回避が key、絶対値定量化の前提
  2. 2 deconvolution model: Fermi function (簡易、empirical) vs BTEX (Blood Tissue EXchange, mechanistic, Chiribiri 推進)
  3. 双璧 KOL: KCL Chiribiri (BTEX, multi-center clinical) + NIH Peter Kellman (PSIR, MOCO, fully-automated workflow)

01qPerfusion 全体技術 cascade — 5 段 + 双璧の研究者

Dual Sequence → Gd Map → MoCo → Deconvolution → MBF/MPR の 5 段技術。

qPerfusion 5 stage cascade
Dual Seq → Gd Map → MoCo → Deconv (Fermi/BTEX) → MBF/MPR

qPerfusion (Quantitative Myocardial Perfusion) は心筋灌流を 絶対値 (mL/g/min) で定量化する WIP/C2P。従来の semi-quantitative (signal intensity ratio) と異なり、multi-center 比較可能 + cut-off 適用可。Common Limitation (WIP=研究用 / Single-vendor / IPA) は Common Limit page に集約。

#Component役割section
1Dual SequenceLV cavity AIF + myocardium 同時取得本 page §02
2Gd Concentration Mapsignal → Gd 濃度変換 (T1 mapping based)本 page §03
3MoCo (Motion Correction)心拍 + 呼吸 + bulk motion 補正本 page §04
4Deconvolution (Fermi or BTEX)AIF + tissue signal → impulse response → MBF本 page §05-06
5MBF / MPR algorithmStress + Rest = MBF、ratio = MPR本 page §07
6臨床例 + Decision Tree三重大 + 東北大 + NCVC routine本 page §08

📚 用語解説 — qPerfusion: Quantitative Myocardial Perfusion。心筋灌流を mL/g/min 絶対値で測定。stress (アデノシン負荷) + rest (安静) で 2 回撮像、MPR = stress / rest。
📚 用語解説 — semi-quantitative vs quantitative: semi = signal intensity ratio (施設 / 装置依存)、quantitative = mL/g/min (multi-center 比較可)。
📚 用語解説 — 関連 ZK note: [[BTEX_Model]] [[Fermi Deconvolution]] [[Kellman PSIR-MOCO]] [[Myocardial Perfusion Reserve]] [[NCVC_qPerfusion]] [[Compare_StressPerfusion_Vendors]] [[4_Perfusion]] 等 6+ files (10_Knowledge/20_Wiki/10_Concepts/)。

🛠️ 運用方法: (1) 医師に「qPerfusion = 5 段 cascade、各 component が定量化精度に影響」と説明 (2) 失敗事例 (MBF が異常値) を component 単位で trouble shoot (3) Chiribiri (KCL) + Kellman (NIH) 双璧を理解、論文選択時に lineage 区別

⚠️ Common Limitation 必読: WIP/C2P = 研究用 (NOT 医療機器) / VA60A 推奨 / IPA 締結必須 / Single-vendor → 詳細 Common Limit page

02Dual Sequence — LV cavity AIF + myocardium 同時取得

AIF saturation 回避が絶対値定量化の前提、2 種 sequence を交互取得。

Dual sequence AIF myocardium
LV cavity AIF (低 FA) + myocardium (高 FA) を交互取得、saturation 回避
要素AIF sequenceMyocardium sequence
対象LV cavity (Gd 濃度高)Myocardium (Gd 濃度低)
FA (Flip Angle)低 FA (~10°)高 FA (~50°)
TR / TE短い長め
Saturation 回避低 FA で linear regime
撮像 timing各 RR interval で交互同上

📚 用語解説 — AIF (Arterial Input Function): 心筋灌流の入力関数、LV cavity 内の Gd 濃度時間曲線。Deconvolution の input、絶対値定量化に必須。
📚 用語解説 — Saturation 問題: LV cavity の Gd 濃度高 (~5 mM peak) で signal が saturation → AIF 過小評価 → MBF 過大評価。低 FA で linear regime に保つ。
📚 用語解説 — Linear regime: signal ∝ Gd concentration の領域。低 FA + short TR で確保、qPerfusion の前提。
📚 用語解説 — Dual sequence の代替: 単一 sequence で AIF + myocardium 両方取得する手法もあるが、saturation 回避が困難 → multi-center 標準は dual。

🛠️ 運用方法: (1) Stress + Rest 各々で dual sequence 取得 (4 撮像/case) (2) AIF curve の peak が saturation してないか確認 (>3 mM で警告) (3) Myocardium signal の SNR 確保のため high FA + 平均化 (4) 撮像 timing を ECG trigger で精密制御

⚠️ アンチパターン: 単一 sequence で AIF + myocardium → saturation で MBF 不正確。AIF curve を目視で peak 確認なし → saturation 見落とし。Stress + Rest の sequence parameter 異なる → MPR (ratio) 計算崩壊。

03Gd Concentration Map — Signal → Gd 濃度変換 (T1 mapping based)

Pre-contrast T1 + post-contrast signal で Gd 濃度を pixel ごと算出。

Gd concentration map
Pre-contrast T1 map + post-contrast signal → Gd 濃度 pixel-wise

📚 用語解説 — T1 mapping: 各 pixel の T1 value を測定する技術。MOLLI (Modified Look-Locker IR) 標準、cardiac syncing 必須。
📚 用語解説 — Gd relaxivity (r1): Gd 1 mM が R1 を増加させる係数、~4 /s/mM @ 3T (Gadovist), 文献値で施設標準化。
📚 用語解説 — Bloch equation 逆計算: 通常 signal = f(T1, FA, TR) を解いて T1 算出。固定 sequence parameter で signal → R1 の look-up table 構築可。
📚 用語解説 — Pre-contrast T1 map の重要性: 患者個人 T1 で baseline 設定、施設間 / 患者間比較可能性向上。

🛠️ 運用方法: (1) Pre-contrast T1 map を必ず撮像 (MOLLI 5(3)3 等) (2) Gd relaxivity は施設で固定値採用 (Gadovist = 4.4, Magnescope = 5.2 @ 3T 等、文献値) (3) Gd map 計算は recon software で自動化、operator 介入最小 (4) Stress + Rest で別々に Gd map 計算 → MBF 算出

⚠️ アンチパターン: Pre-contrast T1 map skip → 仮定値で誤差大。Gd relaxivity を施設で変更 → 過去症例との比較不能。Signal を直接 deconv → Gd map 経由なしで semi-quantitative に逆戻り。

04MoCo for Perfusion — 心拍 + 呼吸 + bulk motion

Stress 時の頻拍 + アデノシン副作用 (深呼吸/咳) で motion 多発、MoCo 必須。

Perfusion MoCo stress
Stress 時 motion (頻拍 + 副作用) → image-based MoCo + outlier rejection
Motion 種類Stress 時の特徴補正方法
Cardiac (心拍)頻拍 (HR ~100-130)ECG gating + bin
Respiratory (呼吸)深呼吸 (副作用)iNAV + soft-gating
Bulk (体動)咳 / 不快k-space rejection
AIF curve driftSlow driftPolynomial detrending

📚 用語解説 — アデノシン負荷: 冠動脈拡張薬、心筋灌流を最大化。副作用 = 頻拍 / 胸部不快感 / 深呼吸 / 咳 → motion 多発。
📚 用語解説 — Kellman PSIR-MOCO: Peter Kellman (NIH) が Cardiac MR の standard MoCo として確立、PSIR + image-based MoCo の組合せ。qPerfusion にも応用。
📚 用語解説 — Outlier rejection: 動きすぎた frame を recon から除外。AIF curve は連続性保つため interpolation。

🛠️ 運用方法: (1) Stress 撮像前に患者に「副作用ある」事前説明 (2) MoCo 失敗時 (image quality 低下) は AIF curve smoothness 確認 (3) Outlier ratio (>20%) は撮像 quality 警告、再撮像検討 (4) Kellman の MoCo paper を施設に共有

⚠️ アンチパターン: アデノシン副作用を skip 説明 → 患者驚愕で motion 増。MoCo を skip → MBF の variance 大幅増、定量化精度低下。Outlier rejection なし → 動いた frame で MBF 異常値。

05Fermi Deconvolution — 簡易 empirical model

Fermi function fitting、計算高速、Cardiac Suite 等で実装。

Fermi deconvolution
Fermi function = empirical impulse response、簡易 fitting で MBF 算出

Fermi function deconvolution: Tissue impulse response R(t) を Fermi function で近似、AIF と convolve して tissue curve を fit。**簡易 + 計算高速**、商用 software (CMR Tools, Circle CVI42) で標準実装。

📚 用語解説 — Fermi function: 物理学の Fermi-Dirac 分布関数の形 (1 / [exp((x-μ)/T) + 1])。Tissue impulse response の近似として 1980s から使用、Jerosch-Herold 1998 (MRM) で広く採用。
📚 用語解説 — Empirical vs mechanistic: Fermi = empirical (実測 fit、機構不問)、BTEX = mechanistic (生理学 model)。
📚 用語解説 — 商用実装: CMR Tools / Circle CVI42 / cvi42 で標準実装、routine 使用可。

🛠️ 運用方法: (1) Fermi は計算高速、初期評価に最適 (2) 三重大 / 東北大 で routine 使用 (3) 商用 software で operator 介入少 (4) BTEX との比較で Fermi の限界を理解

⚠️ アンチパターン: Fermi の物理意義を過剰説明 → empirical なので深追いしない。BTEX と混同 → mechanistic vs empirical の区別必要。

06BTEX Deconvolution — Mechanistic 2-compartment model (Chiribiri 推進)

Blood-Tissue Exchange、生理学的 model、KCL Chiribiri グループ世界リード。

BTEX 2-compartment model
BTEX = Blood-Tissue Exchange、mechanistic 2-compartment model

BTEX (Blood-Tissue Exchange) は心筋を 血管 + 細胞外 + 細胞内 の compartment に分け、Gd の compartment 間移行を記述する mechanistic model。KCL Amedeo Chiribiri グループが multi-center clinical validation を推進、世界 standard 化。

Parameter意味MBF との関係
F (= MBF)Myocardial Blood Flow (mL/g/min)主 output
VpPlasma volume fraction微小血管
PSPermeability-Surface area product膜透過性
VeExtracellular volume fractionECV と関連
ViIntracellular volume fraction細胞内

📚 用語解説 — BTEX (Blood-Tissue Exchange): Bassingthwaighte 1989 の 2-compartment model、心筋灌流の mechanistic model。Fermi より生理学的に正確。
📚 用語解説 — Amedeo Chiribiri: KCL (King's College London) cardiac MR、qPerfusion + BTEX の世界 leader。Multi-center BTEX validation 主推進、ACDC 等 dataset 公開。
📚 用語解説 — BTEX vs Fermi の違い: BTEX = 5 parameters (F + Vp + PS + Ve + Vi)、 fitting 重い + 数値安定性 issue あり。Fermi = 3 parameters、軽い。BTEX が「正しい」が、Fermi が routine 実用。
📚 用語解説 — Manabe との関係: 三重大 / 東北大 qPerfusion で BTEX 比較研究中、Chiribiri との connection 育成中。

🛠️ 運用方法: (1) BTEX は research grade、商用実装少 (2) Chiribiri グループの open dataset (ACDC) で validation (3) 三重大 + 東北大 qPerfusion で BTEX vs Fermi 比較研究 (4) ISMRM/SCMR で Chiribiri と直接交流推奨 (5) 関連 ZK [[BTEX_Model]]

⚠️ アンチパターン: BTEX を「Fermi より良い」と曖昧 push → 商用実装少、運用負担大。BTEX 5 parameters を全て fit しようとする → 数値不安定、PS と Vi は固定推奨。Chiribiri 引用なしで BTEX 説明 → 信頼性低下。

07MBF / MPR Algorithm + Peter Kellman (NIH) Fully-Automated Workflow

Stress + Rest = MBF、ratio = MPR、Kellman の fully-automated workflow が世界基準。

MBF MPR Kellman workflow
Stress + Rest = MBF maps、ratio = MPR、Kellman fully-automated workflow
指標計算正常値異常閾値
MBF restRest deconvolution 結果0.6-1.0 mL/g/min
MBF stressStress deconvolution 結果2.0-4.0 mL/g/min<2.0 で疑い
MPR (Myocardial Perfusion Reserve)MBF stress / MBF rest2.0+ 正常<1.5 = 異常 (虚血)
Segmental MBFAHA 16 segment 別Segment 異常 = local 虚血

📚 用語解説 — MBF (Myocardial Blood Flow): 心筋単位重量あたり灌流量 (mL/g/min)。Stress (アデノシン負荷時) + Rest (安静時) の 2 値。
📚 用語解説 — MPR (Myocardial Perfusion Reserve): MBF stress / MBF rest。冠動脈疾患の機能的指標、CFR (Coronary Flow Reserve) と類似概念。
📚 用語解説 — Peter Kellman: NIH (NHLBI) cardiac MR、元 Siemens collaborator。Fully-automated qPerfusion workflow 確立 (Kellman 2017 JCMR + 2019 Circulation)、operator 介入最小化 + multi-center 標準化。
📚 用語解説 — AHA 16 segment: American Heart Association 心筋 16 区分。Segmental MBF で local 虚血を検出。
📚 用語解説 — Manabe との関係: 三重大 (qPerfusion 2502 active) + 東北大 (active) で Kellman workflow 採用、Chiribiri BTEX と比較中。

🛠️ 運用方法: (1) Kellman fully-automated workflow を施設に紹介、operator burden 最小化 (2) MPR cut-off (>2.0 正常 / <1.5 異常) で診断 (3) AHA 16 segment 報告で local 虚血検出 (4) Chiribiri (KCL) + Kellman (NIH) 双璧の論文を施設に共有 (5) 関連 ZK [[Kellman PSIR-MOCO]] [[Myocardial Perfusion Reserve]]

⚠️ アンチパターン: MBF 単独で診断 → Stress / Rest 両方必須、MPR が key。Operator manual selection 多用 → Kellman 自動化の reproducibility 損失。AHA 16 segment skip → local 虚血見落とし。

08臨床例 + 訪問前 30 秒 Decision Tree + 双璧の論文 5

三重大 + 東北大 + NCVC routine、訪問対応 30 秒判断、KCL/NIH 双璧論文。

qPerfusion clinical decision papers
三重大/東北大 routine + Decision Tree + Chiribiri/Kellman 双璧
施設qPerfusion status代表 use caseWIP/C2P code
三重大学2502_Quantitative_Perfusion activeStress + Rest MBF + MPRWIP/C2P
東北大学qPerfusion active (継続)Multi-vessel disease 評価WIP/C2P
NCVCqPerfusion 関連研究 activeCardiomyopathy + 灌流WIP/C2P
京都大学Neuro 中心、Cardiac qPerf 候補未契約
他施設Cardiac MRI routine 化次第
#Author + YearJournal貢献
1Kellman 2017JCMRFully-automated qPerfusion workflow
2Chiribiri 2013JCMR (or Circ Imag)BTEX vs Fermi multicenter validation
3Hsu 2018JACC Cardiovasc ImagingqPerfusion clinical outcome
4Engblom 2017JCMRFully-quantitative AIF correction
5Brown 2021 (latest)CirculationMulti-center qPerfusion + AI

🛠️ 運用方法 — 訪問前 30 秒 Decision Tree: (1) MAGNETOM Vida/Skyra/Prisma + VA60A+ ? → No → 終了 (2) Cardiac MRI routine ある? → Yes → (3) Stress test (アデノシン) 施設対応? → Yes → qPerfusion 提案 (4) IPA タイムライン提示 (5) Chiribiri + Kellman 論文添付 (6) Fermi (商用) vs BTEX (research) を選択肢提示

⚠️ Common Limitation 必読 (再掲): WIP/C2P = NOT 医療機器 / Single-vendor / IPA 8-12 weeks / Apps Guide 沈黙箇所 / 1.5T-3T extrapolation / SAR (Stress 時高) / 提供責任分界 — 詳細 explainer-common-limit.pages.dev

結論

  1. qPerfusion = 5 段 cascade: Dual Sequence → Gd Map → MoCo → Deconv (Fermi/BTEX) → MBF/MPR。各 component が定量化精度に影響。
  2. 2 deconvolution model: Fermi (empirical, 商用 routine) vs BTEX (mechanistic, KCL Chiribiri 推進、research grade)。Peter Kellman (NIH) の fully-automated workflow が世界基準。
  3. 三重大 + 東北大 + NCVC routine。訪問前 30 秒 Decision Tree + Chiribiri/Kellman 双璧論文 5 添付。Common Limitation は Common Limit page 参照。