拍動する心臓の3D冠状微小血管を鮮明に捉える：新たな超音波画像診断技術

11月 6, 2026 · 論文解説

🎧 この記事を音声で聴く（AI生成Podcast）

解説対象論文: 3D non-rigid motion correction using doppler-regularized optical flow for coronary ultrasound localization microscopy in beating heart
Physics in Medicine and Biology｜被引用数: 0（2026年11月6日時点）

心臓の健康を左右する重要な血管ネットワークである冠状微小血管は、その小ささから従来の画像診断では詳細な観察が困難でした。特に、拍動する心臓の動きは、超高分解能イメージング技術である超音波局在顕微鏡（ULM）の精度を低下させる大きな課題です。このたび発表された革新的な研究は、この心臓の激しい動きを高度なアルゴリズムで補正し、これまで不可能だった3Dでの冠状微小血管ネットワークの鮮明な可視化と精密な血流速度測定を実現しました。この技術は、冠状微小血管疾患の理解と診断に新たな光をもたらす可能性を秘めています。

研究サマリー（FINER & PICO）

Feasible実施可能性	既存の超音波技術とコンピュータビジョンの手法を組み合わせ、GPU等で計算効率を向上可能。
Ethical倫理面	ラットの心臓を用いた研究は、倫理委員会および動物管理委員会の承認を得て実施された。
Interesting面白さ	拍動する心臓の3D微小血管をサブマイクロメートルレベルで可視化・定量化する新技術は非常に興味深い。
Novel新規性	3D冠状ULMにおけるドップラー正則化オプティカルフローを用いた非剛体運動補正は新規性がある。
Relevant切実さ	冠状微小血管疾患の診断と病態理解を深め、将来の治療法開発に貢献する。
Measurable測定可能性	マイクロバブル追跡密度、血管の描出度、速度測定分解能（少なくとも40 µm）で定量的に評価された。
Modifiable派生・改善	オプティカルフローの平滑化パラメータとドップラー重みパラメータを調整可能。
Structured文章構造	論文は序論、材料と方法、結果、考察、結論と明確な構成で、研究内容が体系的に記述されている。
PICO / PECO対象・介入・比較・結果	P:拍動する心臓の3D冠状微小血管のULMイメージング、I:ドップラー正則化オプティカルフローによる3D非剛体運動補正、C:補正なしのULM、O:より正確な微小血管ネットワークの再構成と可視化、微小血管の鮮明化、速度測定。

はじめに

どうも、Beyond the Pixelです。心臓の健康は、私たちの生活の質に深く関わる重要なテーマです。特に、冠状微小血管疾患（CMD）は、直径300 µm未満の血管が関係するもので、心血管疾患による罹患率と死亡率に大きく寄与しています。しかし、現在の臨床画像診断技術では、これらの微小な血管ネットワークを直接観察するのに必要な空間的・時間的分解能が不足しており、診断が困難なのが現状です。この課題を解決する可能性を秘めているのが、回折限界を超えたイメージングを可能にする超音波局在顕微鏡（ULM）です。ULMは冠状微小血管構造を明らかにすることができますが、心臓の大きな動きによって画像が影響を受けてしまうという問題がありました。今回ご紹介する論文は、この3D冠状ULMにおける非剛体運動の課題に対処し、冠状微小血管ネットワークのより正確な再構成と可視化を可能にする画期的なフレームワークを提案しています。

画期的なアプローチ：ドップラー正則化オプティカルフロー

本研究チームは、拍動するラットの心臓の3D ULMにおける心臓の動きを補償するため、ドップラー速度正則化を用いたオプティカルフローに基づく3D非剛体運動補正フレームワークを開発しました。これにより、冠状微小血管ネットワークのより正確な再構成と可視化を目指しました。このフレームワークは、心臓の動きを詳細に推定し、画像に写る微小バブルの位置を正確に補正することで、これまで困難だった微小血管の鮮明な画像化を実現します。

**Figure 1.** Figure 1 shows the ultrasound volumetric image processing pipeline. First the motion inside each cardiac cycle is corrected.

オプティカルフローによるモーション推定

オプティカルフロー法は、Bモード画像（超音波の断面画像）上の動きを推定するために用いられます。これは、画像内のピクセルやボクセルの明るさの変化から、その動きの方向と速度を計算する手法です。しかし、明るさの変化だけでは動きの推定が不十分な場合があるため、この問題の解決には正則化が必要です。この研究では、各運動場の勾配に古典的なL2ノルム正則化が採用されました。

ドップラー速度による正則化

この研究の重要な要素は、オプティカルフロー法にドップラー速度正則化を組み込んだ点です。ドップラー法は、超音波信号の位相シフトから組織の軸方向（プローブに沿った方向）の速度を非常に正確に推定できます。このドップラー速度の情報をオプティカルフローの計算に組み込むことで、軸方向の動き推定の精度と安定性が向上し、全体的な3D運動場の信頼性が高まります。これは、ポレーらが2次元で提案した方法を3次元に拡張したものです。

2段階の動き補正パイプライン

この動き補正フレームワークは、2段階のパイプラインで構成されています。まず、各心臓サイクル内の動きを補正します。これには、IQボリューム（超音波の生データに近い信号）からドップラー軸方向速度を計算し、ドップラー正則化を伴うシングルスケールオプティカルフローを用いて連続するボリューム間の3D変位場を取得します。次に、これらの変位場を統合して局在化されたマイクロバブルの座標を補正します。次に、異なる心臓サイクル間の動きを補正します。これは、左心室圧曲線を同期に用いて、ミッドダイアストール（心臓サイクル中の同じ瞬間）で2つのBモードボリュームを選択し、マルチスケールオプティカルフローを用いて変位場を計算することで行われます。このサイクル間の変位場は、既にサイクル内補正が施されたマイクロバブル座標に適用されます。

**Figure 2.** Fig. 2. Experimental setup. A rat heart was isolated and cannulated at the ascending aorta. A peristaltic pump was used to provide a solution of

驚くべき結果：微小血管の鮮明化と速度測定

この方法は、超高速3D超音波で撮像された、灌流された摘出拍動ラット心臓で検証されました。動き補正はマイクロバブルの追跡密度と血管の描出を向上させ、100 µm未満の微小血管の可視化と、少なくとも40 µmの空間分解能での速度測定を可能にしました。

**Figure 3.** Fig. 3. Estimation of the cardiac motion with the volume-to-volume Optical Flow. Tissue velocities maps are displayed as arrows on the two

例えば、心尖部の累積変位はY方向で500 µmを超え、これはボクセルサイズよりもはるかに大きいことが示されました。心臓の収縮期から弛緩期にかけての組織速度のダイナミクスも詳細に可視化されました。

**Figure 4.** Fig. 4. 3D trajectories of tissue points superimposed on the B-mode volume. The trajectories of tissue points are displayed in 3D over the mean

心筋内の異なる点の3D軌跡は、心臓の領域間で動きが大きく異なることを示しており、この研究で提案された非剛体運動補正の重要性を強調しています。

3D冠状ULMへのモーション補正の適用

動き補正は、3D密度マップを大幅に向上させることが示されました。動き補正なしでは、最大の冠状動脈は見えるものの、全体的に密度マップはぼやけており、直径200 µm未満の血管を分解することは困難でした。しかし、オプティカルフローに基づく動き補正を適用すると、ULM密度マップはより密度の高い血管と鮮明な境界を示し、再構成されたボクセルあたりのマイクロバブル追跡数が増加しました。

**Figure 5.** Fig. 5. Effect of the motion correction on the 3D coronaries ULM density map. ULM density maps without (a) and with motion correction (b). Here the MBs density is the number of MBs per rasterized voxels of 30x30x30 µm3.

この密度増加は、直径100 µmを超える大きな血管と、直径100 µm未満の微小血管の両方で観察されました。

**Figure 6.** Fig. 6. Effect of the motion correction on the MBs density distribution. Histograms showing the number and density of MBs for big vessels (diameter > 100 µm) (a) and small vessels (diameter < 100 µm) (b) for 3 different hearts.

これは、3つのラット心臓すべてで一貫して確認されています。特に微小血管の可視化が大幅に改善され、補正なしでは血管が2〜4倍広く見えることがあったのに対し、補正後はY字型分岐の小さい血管が区別できるようになり、冠状ULMの分解能が向上しました。

**Figure 7.** Fig. 7. ULM reconstruction using motion correction. Close up view of ULM density map (Left) both without (top – left) and with (bottom – left) motion correction. Cross-section of the density is shown for 4 vessels (1) to (4) to quantify the effect of motion correction (right).

さらに、動き補正は速度推定にも影響を与え、3D冠状ULM速度マップは動き補正法の適用後に大幅に強化されました。直径360 ± 20 µmの大きな冠状血管の断面におけるマイクロバブルの速度分布は、隣接するボクセル間で統計的に区別可能であり、少なくとも40 µmの速度測定分解能が示されています。

**Figure 8.** Fig. 8. Effect of the motion correction on 3D coronary ULM velocity map. ULM velocity maps without (left) and with (middle) motion correction with the color representing the velocities of the tracked MBs. The velocity distribution profile of the MBs tracked across the width of a selected (5) large coronary vessel (diameter of 360 ± 20 µm) of the ULM map was obtained after motion correction (right). Between each adjacent voxel statistical significance was inferred using an unpaired Student t-test. P > 0.05: no significance (NS), *P < 0.05, **P < 1e-2, ***P < 1e-3, ****P < 1e-4.

研究の意義と今後の展望

本研究で開発されたドップラー正則化オプティカルフローに基づく3Dモーション補正フレームワークは、拍動するラットの心臓における冠状微小循環のロバストなマッピングを可能にしました。このアプローチは、激しい非剛体運動があるにもかかわらず、血流測定を安定させ、希少な灌流微小血管を保存します。これにより、血管内のマイクロバブル検出が強化され、高分解能の可視化を損なうことなく取得時間を短縮できます。この成果は、3D冠状ULMの生体内および臨床応用への将来的な移行を支持するものです。今回の研究では、灌流された摘出拍動ラット心臓での実験が行われましたが、これは生理学的に関連性が高いものの、生体内での血行動態や自律神経調節、周囲組織の影響を完全に捉えるものではありません。今後の研究では、このフレームワークを完全に生体内のモデルで実装し、生理学的条件下での性能を評価することが課題となります。また、オプティカルフローアプローチはノイズに敏感であるという潜在的な限界も指摘されています。ドップラー正則化は軸方向速度推定を改善しますが、ごく低いSN比や乱流のある血管では効果が限定的になる可能性があります。さらに、人間の心臓のより遅い心拍数や異なる動きのダイナミクスを考慮すると、より大きな心臓や臨床画像診断への適用には、取得ハードウェア、動き推定スケール、および処理パイプラインの適応が必要となるでしょう。これらの技術的課題への対処は、提案されたフレームワークを実験設定からヒトの心臓イメージングへと移行させ、その臨床的適用性を評価するための重要なステップとなります。

結論

この研究は、激しい非剛体運動下にある拍動するラット心臓の冠状微小循環のロバストなマッピングを可能にする、3Dオプティカルフローに基づく動き補正フレームワークを開発しました。ドップラー正則化オプティカルフローは、血流測定を安定させ、希少な灌流血管を鮮明に保ちます。このアプローチにより、血管内のマイクロバブル検出が向上し、高分解能の可視化を損なうことなく取得時間が短縮されます。これらの結果は、3D冠状ULMの生体内および臨床応用への将来的な移行を強力に支持するものです。

用語集

冠状微小血管疾患（CMD）: 心臓を栄養する直径300 µm未満の微小な血管に異常が生じる疾患。
超音波局在顕微鏡（ULM）: 超音波造影剤（マイクロバブル）を一つずつ追跡し、従来の超音波の回折限界を超えて微小血管を可視化する技術。
回折限界: 光学機器や超音波機器で、波の性質によって定められる、それ以上小さくできない空間分解能の理論的な限界。
非剛体運動: 物体の全体的な形状や内部構造が変化しながら動くような動き（例：拍動する心臓）。
オプティカルフロー: 連続する画像フレーム間のピクセルの動き（変位ベクトル）を推定するコンピュータビジョン技術。
ドップラー速度正則化: ドップラー効果を用いて測定された速度情報をオプティカルフローの計算に組み込み、動き推定の精度を向上させる手法。
マイクロバブル（MBs）: 超音波造影剤として血液中に注入される微小なガスの泡。超音波で追跡され、血流を可視化する。
超高速超音波: 通常の超音波よりもはるかに速いフレームレートで画像を生成する技術。
時空間クラッターフィルター: 超音波画像から、組織の動きによる不要な信号（クラッター）を除去するフィルター処理。
特異値分解（SVD）: 行列を3つの行列に分解する手法で、画像処理では不要な信号（クラッター）と目的の信号（マイクロバブル）を分離するために用いられる。
ランゲンドルフ灌流: 単離した動物の心臓を体外で灌流液によって栄養し、拍動状態を維持する実験手法。
IQボリューム: 超音波信号のI（同相）成分とQ（直交）成分から構成される、生データに近い複素ボリュームデータ。
心臓サイクル: 心臓が1回の拍動を行う一連の過程。収縮期と拡張期からなる。
拡張期: 心臓が弛緩し、血液が心室に充満する時期。
収縮期: 心臓が収縮し、血液を大動脈や肺動脈に送り出す時期。