An Open, Reproducible Gamma-Variate Pipeline for CT-Perfusion Time-Attenuation Curve Analysis, with Standardized (ASIST-Japan) Map Visualization

技術要約：CT灌流のためのオープンで再現可能なガンマ変量パイプライン

問題提起
CT灌流（CTP）は急性期脳卒中や腫瘍学的画像診断において極めて重要であるが、商用ベンダーのソフトウェア間における再現性の欠如により、その定量的信頼性が損なわれている。異なるパッケージで同一のソースデータを処理しても、デコンボリューション・アルゴリズム、実装方法、および閾値の差異により、得られる灌流マップや値が著しく異なる。さらに、これらのマップの可視化は、非標準的なルックアップテーブル（LUT）に起因するサイト間の変動に悩まされている。つまり、同一の定量的値であっても、使用されるカラースケールによって視覚的に異なって見えるのである。この「ブラックボックス」的な不透明さは、臨床医や研究者が基礎となる処理ステップを監査または検証することを妨げている。

手法
著者らは、CTP解析のための透明かつパラメータが公開されたパイプラインを提供する、純粋なPython（>=3.9）ベースのオープンソース・ライブラリである ctp-core (v0.1.0) を提示する。本手法は、主に以下の2つのコンポーネントに焦点を当てている：

1. ガンマ変量（GV）フィッティングと解析的指標の導出：

   • モデル： 第一次通過時－減衰曲線は、ガンマ変量関数 $C(t) = K (t −t_0)^\alpha \exp(-\frac{t - t_0}{\beta})$ を用いてモデリングされる。ここで、$t_0$ は開始時刻、$\alpha$ は形状パラメータ、$\beta$ はスケールパラメータ、そして $K$ は係数である。
   • フィッティング： 曲線は、物理的な制約条件（bounds）を設けた非線形最小二乗法（Trust Region Reflective アルゴリズム）を用いてフィッティングされる。初期値はデータ駆動型である（例：$t_0$ は10%ピーク通過点より前に配置される）。
   • 解析的導出： 再最適化を行う代わりに、灌流指標はフィッティングされたパラメータ（$K, t_0, \alpha, \beta$）から解析的に導出される：
     • ピーク時間（$t_{peak}$）： $t_0 + \alpha\beta$
     • ピーク振幅（$y_{peak}$）： $K (\alpha\beta)^\alpha e^{-\alpha}$
     • 曲線下面積（AUC）： $K \beta^{\alpha+1} \Gamma(\alpha + 1)$
     • ボラス到着時間（BAT）： $t_0$
   • 堅牢性： パイプラインには、フィッティングが失敗した場合に例外を発生させるのではなく、NaNを返す「フィットフリー（fit-free）」の生データパスが含まれている。失敗は、エラーメッセージ（例：形状の不一致、オプティマイザの失敗）と共に明示的にフラグ立てされる。

2. 標準化された可視化（ASIST-Japan a-LUT）：

   • 本パイプラインは、値を黒–青–シアン–緑–黄–赤のスペクトルへとマッピングする256レベルのRGBテーブルである、ASIST-Japan標準ルックアップテーブル（a-LUT）を実装している。
   • スカラー値は、標準的な表示範囲（例：CBF 0–80 mL/100 g/min）に対して線形にウィンドウ処理され、決定論的にインデックス化される。
   • 極めて重要な点として、カラーマッピングは可視化にのみ影響を与え、基礎となる定量的なボクセル値は変更されず、別途保存される。

3. スコープと実装：

   • オープンコア（ctp-core）はスタンドアロンであり、GUI、DICOM I/O、または患者データの取り扱いは含まない。
   • デコンボリューションおよび動脈入力関数（AIF）検出による完全なパラメトリックマッピング（CBF, CBV, MTT, Tmax）は文書化され実装されているが、これらは別のプライベートなアプリケーション内に存在する。オープンコアの再現性に関する主張は、GVフィッティングと合成データによる検証に限定される。

検証と結果
本パイプラインは、機密情報や患者データを使用せず、既知の正解（ground truth）を持つ合成モンテカルロ・データのみを用いて検証された。

• 実験設定： 信号対雑音比（SNR）が5から100の範囲で、200回の独立した実行が行われた。
• 性能：
  • 代表的なSNR 20において、パイプラインはピーク時間を0.13 s以内、ピーク振幅を2.0%以内、BATを0.51 s以内（平均絶対誤差）の精度で回復した。
  • フィットの品質は $R^2$ 0.98を達成した。
  • 誤差はノイズとともに単調に悪化した。SNR 5では、ピーク時間誤差は約0.52 s、振幅誤差は約8.1%であった。
  • フィットから導出されたピークは、一貫して生の最大値（フィットフリー・パス）よりも正解に近い値を示した。
• 再現性： コードは決定論的（固定シード）であり、テストスイート（34個のテスト）はビット安定な指標を確認している。合成データ生成器および評価スクリプトは、コアと共にリリースされている。

主要な貢献

1. 透明性： すべてのパラメータ、境界、および処理ステップが公開されている、完全にオープンでパラメータ透過的なパイプラインであり、これは商用の「ブラックボックス」ソリューションと対照的である。
2. 標準化： ASIST-Japan a-LUTの統合により、定量的な値を変更することなく、パラメトリックマップの可視化をサイト間で比較可能にしている。
3. 再現性： 合成データの生成および検証ベンチマークを含むワークフロー全体がMITライセンスの下でリリースされており、独立した監査と再利用を可能にしている。
4. 検証： 幅広いSNR範囲における既知の正解に対する回復精度を実証し、フィッティング失敗への明示的な対処を行っている。

意義と主張
著者らは、ctp-core を臨床用または規制グレードのツールとしてではなく、研究用リファレンスおよび再現可能なベンチマークとして位置づけている。論文では、透明なガンマ変量フィッティングエンジンと標準化された可視化を提供することで、CTP後処理における再現性の危機に対処できると主張している。本研究は、他者が灌流解析の数学的基礎を監査し、コードを再利用して独立した検証を行い、検証済みのベースラインの上に構築することを可能にするものである。著者らは、このプレプリントにおいて実際の患者データを用いた臨床的検証は主張していないことを明記しており、外部の公開データセット（ISLESなど）を用いた検証は今後の課題としている。