この記事では、3D U-Net畳み込みニューラルネットワークを使用して、体重負荷CBCT画像で重力負荷を受けた人間の脊椎をセグメント化するためのプロトコルを紹介します。この方法は、脊椎のアライメントを研究するための貴重なツールを提供し、耐荷重条件下での脊椎の3Dプリンティングと生体力学の研究をサポートします。
方法論記事
* These authors contributed equally
この記事では、3D U-Net畳み込みニューラルネットワークを使用して、体重負荷CBCT画像で重力負荷を受けた人間の脊椎をセグメント化するためのプロトコルを紹介します。この方法は、脊椎のアライメントを研究するための貴重なツールを提供し、耐荷重条件下での脊椎の3Dプリンティングと生体力学の研究をサポートします。
脊柱側弯症や変性椎間板疾患などの筋骨格系疾患の診断と管理には、重力負荷条件下での脊椎のアライメントの正確な分析が不可欠です。2DX線写真や仰臥位のCTやMRIなどの従来のイメージング技術では、荷重を受けた状態の脊椎の完全な3D形状をキャプチャすることはできません。体重を支えるコーンビームコンピュータ断層撮影法(CBCT)は、立った状態で高解像度の3D画像を提供し、脊椎のアライメントをより正確に評価することができます。この研究では、3D 畳み込み層と残留接続を備えた U-Net 畳み込みニューラル ネットワーク (CNN) を使用してこれらの画像をセグメント化するためのディープラーニング ベースのプロトコルを提示します。このプロトコルには、画像取得、手動アノテーション、前処理、およびモデルトレーニングが含まれます。解剖学的構造のセグメンテーションが成功することで、コブ角や椎体の回転など、臨床的に関連する脊椎パラメータを正確に測定することができます。さらに、セグメント化された3Dモデルは、3Dプリント用にSTL形式でエクスポートできるため、手術計画や患者教育が容易になります。さらに、このプロトコルは、下肢など、重力負荷条件下で評価が必要な他の解剖学的構造のセグメンテーションモデルのトレーニングに適合させることができます。この方法は、重力負荷がかかる条件下で脊椎やその他の解剖学的構造を正確にセグメンテーションするための信頼性と適応性のあるツールを提供し、臨床診断と調査研究の両方で実用化できます。
脊柱側弯症、変性椎間板疾患、脊椎変形などの筋骨格系障害の診断と管理には、重力負荷条件下での脊椎のアライメントの正確な分析が不可欠です1,2,3,4。2Dレントゲン写真などの従来のイメージング技術では、脊椎の完全な3D形状をキャプチャすることはできません。従来のCTやMRIでは3D画像を生成することができますが、これらのモダリティでは患者を仰臥位で画像化する必要があります。対照的に、体重を支えるコーンビームコンピュータ断層撮影法(CBCT)は、患者が立った状態で高解像度の3D画像を提供し、脊椎のアライメント5をより詳細に評価することができます。
3Dイメージングによる骨格構造のセグメンテーションは、脊椎変形の診断や手術の計画に必要な、コブ角、椎体回転、骨盤発生率など、臨床的に関連性のある脊椎パラメータを抽出するために重要です6。さらに、セグメンテーションにより、解剖学的構造の3Dモデルの作成が可能になり、患者教育と手術準備がサポートされます7。ディープラーニング技術、特に畳み込みニューラルネットワーク(CNN)は、医用画像解析、特に画像からの複雑な解剖学的構造のセグメンテーションに広く使用されています8,9。ディープラーニングベースのセグメンテーションにより、手作業の負荷が軽減され、オペレーターのばらつきが最小限に抑えられ、脊椎イメージングの精度と一貫性が向上します。
この記事では、体重負荷CBCTを使用して重力負荷スパインをセグメント化するためのディープラーニングベースのプロトコルを紹介します。3D U-Net CNNは、耐荷重条件下で椎骨、骨盤、大腿骨頭を自動的にセグメント化するために使用されます。この方法を日常的な体重負荷CBCTと統合することで、より正確な脊椎評価のための標準化された再現性のあるツールが提供され、それによって臨床的意思決定と患者の転帰の改善をサポートします。このアプローチは、重力負荷を受けた下肢などの評価が必要な他の解剖学的構造のセグメンテーションモデルのトレーニングに複製可能で適応可能です。
この研究は、ヘルシンキ宣言に概説されている倫理基準に従って実施され、四川大学西中国病院の治験審査委員会によって承認されました (承認 834 / 2021)。この研究に参加した参加者は、退行性脊柱側弯症または青年期特発性脊柱側弯症(AIS)と診断されました。除外基準には、脊椎手術の既往歴のある人、脊椎インプラントの存在、または画像診断中に立位を維持できない人が含まれていました。研究に関与したすべての参加者は、画像の取得と研究への参加に先立って、インフォームド コンセントを提供しました。すべての参加者のプライバシーと機密性は、研究全体を通じて厳格に維持されました。この調査で使用した機器とソフトウェアは、 資料表に記載されています。
1.脊椎のイメージング
2. ITK-SNAPによるアノテーション
3. 画像の前処理
注:このプロトコルでは、入力画像と対応するターゲットラベルデータの画像前処理は、Pythonの統合開発環境(IDE)である PyCharmを使用して実装されます。ステップは、読みやすさを向上させ、主要なステップを明確に概説するために 、補足ファイル 1 に擬似コードとして提供されている Python コードを介して実行されます。
4. モデルトレーニング
セグメンテーションの結果
この研究では、脊椎、骨盤、大腿骨頭の3つの解剖学的構造をセグメント化しました。骨盤と大腿骨頭が含まれていたのは、特定の脊椎パラメータが測定のためにこれらの構造からのランドマークを必要とするためです。各構造のデータセットは、 表 2 に示すように、トレーニング セットとテスト セットに分割されました。
セグメンテーションのパフォーマンスは、Dice 係数、ボリューム類似性 (VS)、再現率 (Recall) の 3 つの標準メトリックを使用して評価されました。



ここで、TP は真陽性、FP は偽陽性、FN は偽陰性を示します。
ネットワークを 500 エポックで学習させた後、セグメンテーションの結果はメトリクスに基づいて評価されました。各解剖学的構造の性能を 表3に示します。これらの結果は、提案されたU-Netベースのセグメンテーションモデルが脊椎、骨盤、および大腿骨頭全体にわたって有効性を示しています。 図 2 は、解剖学的構造の代表的なセグメンテーション結果を示しています。
セグメンテーション後のアプリケーション
臨床的に関連性のある3Dパラメータの測定
自動セグメンテーションによって生成された3Dモデルは、コブ角、椎体回転、胸部後弯症、腰椎前弯症、骨盤発生率、骨盤傾斜、仙骨傾斜など、臨床的に重要なさまざまなパラメータを測定するために使用できます。ここでは、セグメント化されたデータを処理し、脊椎と関連構造の正確かつ自動的な測定を提供する自社開発のプログラムを使用して、これらのパラメータの測定値を表示します(図3)。
3Dプリンティング
AIの自動セグメンテーションと後処理で取得した3Dモデルは、3Dプリントソフトで広く対応している一般的な3Dモデルファイル形式であるSTL形式にエクスポートすることができます。本研究では、紫外線(UV)硬化型樹脂を印刷材料として用いた産業用光造形(SLA)プリンターを使用しました。セグメント化されたモデルの分離された部分間で接続が適用され、単一の統一モデルとして印刷できるまとまりのある構造が確保されました。このプロセスにより、患者の解剖学的構造の物理的な複製が可能になり、手術計画、患者とのコミュニケーション、教育ツールに使用できます。 図4 は、AIでセグメント化されたデータから導き出された3Dプリントモデルを示しています。

図1:体重負荷コーンビームコンピュータ断層撮影(CBCT)イメージングのための患者のポジショニング。 患者は立っているプラットフォームに座り、スキャン中の動きや転倒を防ぐためにベルトで固定されます。患者はオーバーヘッドサポートをつかんで安定した位置を維持し、CBCTシステムが重力負荷条件下で脊椎の3D画像をキャプチャします。 この図の拡大版を表示するには、ここをクリックしてください。

図2:脊椎、骨盤、大腿骨頭の軸方向、矢状方向、冠状方向の図でのセグメンテーション結果。 最初の列には CBCT 画像が表示され、その後にグラウンド トゥルース、ディープ ラーニング モデルから取得したセグメンテーション結果が表示されます。2列目と3列目の赤いオーバーレイは、3つのビューすべてで椎体、骨盤、大腿骨頭のセグメント化された領域を示しています。画像は、解剖学的構造をキャプチャするモデルの精度が高く、グラウンドトゥルースと密接に一致することを示しています。 この図の拡大版を表示するには、ここをクリックしてください。

図3:セグメント化された脊椎と骨盤の構造に関する臨床的に関連するパラメータの測定 (A)冠状面の椎骨T11とL3の間のコブ角(8.9°)の測定。(B)矢状面における胸部後弯症(T1-T12:21.2°)の測定。(C)矢状面における腰椎前弯症(L1-L5:26.4°)の測定。(D)軸面の骨盤に対するL2の椎骨回転(VR:3.3°)。(E) 矢状図で測定された骨盤傾斜 (PT: 21.0°)、仙骨傾斜 (SS: 26.1°)、骨盤発生率 (PI: 47.1°) を含む骨盤パラメータ。 この図の拡大版を表示するには、ここをクリックしてください。

図4:脊椎と骨盤のSTLモデルと3Dプリントされた物理モデル。 左の画像は、3Dプリント用に準備されたSTL形式の脊椎と骨盤の3Dモデルを示しています。右の画像は、光造形(SLA)技術を使用して印刷された対応する物理モデルを示しています。 この図の拡大版を表示するには、ここをクリックしてください。
| パラメーター | デフォルト設定 |
| チューブ電圧 | 110キロワット |
| チューブ電流 | 6 mA |
| フレームレート | 12コマ/秒 |
| 再構築FOV | 350ミリメートル |
| スライスの厚さ | 2ミリメートル |
| スライス間の間隔 | 1ミリメートル |
表 1: 露出と再構成のデフォルト パラメーター。脊椎全体の体重を支えるコーンビームコンピューター断層撮影イメージングのデフォルトの露出および再構成パラメータ。
| ラベル | トレーニングデータセット | データセットのテスト |
| 背骨 | 100 | 20 |
| 骨盤 | 57 | 6 |
| 大腿骨頭 | 47 | 5 |
表 2: CBCT 画像モデルのトレーニングとテストに使用されるデータセット。この表は、脊椎、骨盤、大腿骨頭の構造をセグメント化するためのモデルのトレーニングとテストに使用された画像の数を示しています。
| ラベル | サイコロ係数 | ボリュームの類似性 | 思い出す |
| 背骨 | 0.93 | 0.954 | 0.89 |
| 骨盤 | 0.925 | 0.965 | 0.93 |
| 大腿骨頭 | 0.96 | 0.971 | 0.959 |
表3:さまざまな解剖学的構造(脊椎、骨盤、大腿骨頭)に対する深層学習モデルのセグメンテーションパフォーマンス。パフォーマンスメトリクスには、Dice Coefficient、Volume Similarity (VS)、Recallがあり、これらはすべて0から1の範囲で、1に近い値はセグメンテーションの精度が高いことを示します。
補足ファイル1:画像前処理(ステップ3)とモデルトレーニング(ステップ4)の擬似コード。このファイルをダウンロードするには、ここをクリックしてください。
このプロトコルでは、体重負荷CBCTを使用して重力負荷を受けた脊椎画像をセグメント化するためのディープラーニングベースの方法を導入しています。主な手順には、画像取得、手動注釈付け、前処理、モデルトレーニングが含まれます。患者の適切なポジショニングと露出設定は、高品質の画像を得るために重要です。位置ずれや不適切な配置は、最適でない画像につながる可能性があり、ニューラルネットワークが構造を正確にセグメント化することが困難になります。前処理も重要なステップであり、このプロトコルでは、ノイズの追加、明るさの変動、回転、および変換が含まれます。これらの前処理手法は、ニューラルネットワークに入力する画像を準備するように設計されており、さまざまなタイプの画像歪みにさらされることでモデルのロバスト性を高めます。
手動アノテーションは、プロトコルのもう1つの重要な要素です。手動アノテーションの精度は、脊椎体、骨盤、大腿骨頭などの構造を学習してセグメント化するモデルの能力に直接影響します。このプロトコルは、椎体のセグメンテーションに焦点を当てており、コブ角や椎体の回転など、臨床的に最も重要なパラメーターを測定するのに十分です。完全な椎骨 (椎弓板、横突起、棘突起など) を必要とするアプリケーションの場合、モデルは追加のラベル付きデータを使用してトレーニングする必要があります。さらに、手動アノテーションのばらつきも課題となる可能性があり、結果の再現性に影響を与える可能性があります。これを軽減するには、経験豊富な複数のアノテーターを配置し、オブザーバー間の信頼性チェックを実装することをお勧めします。
U-Netアーキテクチャは、3D畳み込み層と残留接続を備えており、3D画像12,13,14の処理に適しています。マルチレベル損失機能は、より深いレイヤーからの詳細を保持し、セグメンテーションの精度を確保するのに役立ちます。ただし、学習率や拡張設定などの学習パラメーターを慎重に調整するには、過適合を防ぎ、モデルのロバスト性を向上させる必要があります。最近の研究では、トランスフォーマーやYOLOアーキテクチャ15,16などのニューラルネットワークを使用して、特定のタスクに合わせた脊椎セグメンテーションまたは検出モデルが提案されていますが、これらは通常、CTまたはMRI脊椎画像に依存しています。対照的に、このプロトコルは、ピクセルレベルのセグメンテーションタスク17でより一般的に使用されるU-Netアーキテクチャを採用している。さらに重要なことに、この方法は、直立CBCTを使用して重力負荷条件下で脊椎をセグメント化するように設計されており、立った状態または体重負荷状態の脊椎の評価を可能にします。このユニークな機能により、脊椎の変形を診断し、治療戦略を計画するために重要な、脊椎のアライメントと機械的ストレスをより現実的に評価できます。特に、脊椎セグメンテーションにおけるダイス係数0.93で示されているように、モデルのパフォーマンスは、脊椎セグメンテーション9,18,19にCTおよびMRIデータを使用する他のディープラーニングベースのモデルのパフォーマンスと同等またはそれを超えています。
この方法にはいくつかの制限があります。この研究で使用されたデータセットは、椎体、骨盤、および大腿骨頭をセグメント化するのに適していますが、まれまたは複雑な脊椎変形をより適切に処理するための拡張の恩恵を受けます。データセットが大規模で多様であればあるほど、モデルをより広範なケースに一般化するのに役立ちます。また、CBCTイメージングは金属アーチファクトに敏感であるため、術後の脊椎を金属インプラントでセグメント化するのが難しいという制約もあります。さらに、この方法は計算量が多いため、強力なGPUへのアクセスが必要であり、すべての臨床および研究環境で利用できるわけではありません。手作業によるアノテーション処理は時間がかかり、そのため、アノテーションの精度をコンセンサスで検証するために複数のアノテーターを巻き込むことができませんでした。さらに、このモデルはまだ外部データセットで検証されていないため、他の臨床環境への一般化可能性は限られています。
脊椎関連のパラメータを正確にセグメント化して測定する能力は、特に脊柱側弯症のような症例において、診断と治療計画を強化します。モデルをSTL形式でエクスポートして3Dプリントする機能は、術前計画にも大きな応用があり、外科医は患者固有のモデルを練習したり参照したりできます。これらの具体的なモデルは、患者教育のための貴重なツールであり、患者が自分の状態と治療をよりよく理解するのに役立ちます20。脊椎イメージング以外にも、このプロトコルは、下肢や関節など、重力負荷を受ける他の解剖学的構造のセグメンテーションにも適応でき、変形性膝関節症や股関節変形などの場合に特に価値があります。さらに、このプロトコルは、荷重負荷条件下での脊椎または他の構造の詳細な3Dモデルが力学および応力分布の分析に重要であるバイオメカニクス研究に潜在的に応用できる可能性がある21,22。これにより、奇形の進行に関する貴重な洞察が得られ、新しい治療戦略の開発を導くことができます。
結論として、このプロトコルは、重力負荷条件下で脊椎画像をセグメント化および分析するための信頼性と効率的な方法を提供します。いくつかの制限がありますが、診断、外科計画、バイオメカニクス研究、教育への応用が期待できるため、臨床と研究の両方のシナリオで貴重なツールとなります。
著者らは、議定書に記載されている作業に関連して競合する利益を宣言しない。
この研究は、四川省の主要な研究開発プロジェクト(2024YFFK0050および2023YFG0126)、四川大学西中国病院の1.3.5プロジェクトピーク規律開発基金人工知能プロジェクト(ZYAI24013)、およびハイエンド人材サポートプログラムの1.3.5プロジェクト(ZYGD23024)の支援を受けました。
| 名前 | 会社 | カタログ番号 | コメント |
|---|---|---|---|
| 3Dプリンター(G1400) | UnionTech、上海、中国 | G1400 | 産業用光造形(SLA)3Dプリンター。 |
| ディープラーニング モデルのトレーニングのための | CUDA バージョン 11.4 | NVIDIA | 11.4 |
| ITK-SNAP | 3D生物医学画像のセグメンテーションとアノテーションのための | オープンソース | 4.0.2 |
| Numpy | 画像 | 1.24.4 | python ライブラリ、Web サイト: https://numpy.org/about/ |
| PyCharm | JetBrains | 2023.2.1 | Python プログラミングに使用される統合開発環境 (IDE)。 |
| Pythonバージョン3.8 | 3.8 | ラーニングモデル開発のためにPyTorchと組み合わせて使用します。 | |
| PyTorch バージョン 1.10 | オープンソース | 1.10 | セグメンテーションモデルのトレーニングに使用されるディープラーニングフレームワーク。 |
| SimpleITK | 画像 | 2.3.1Python | ライブラリ、Webサイト:https://simpleitk.org/ |
| UV-Curable Resin | UnionTech、上海、中国 | UV-001SLA3D | 印刷プロセスで使用される材料。 |
| 重量ベアリングコーンビームCTシステム (WR3Dの) | Angell Technology, Shenzhen, China | DX6290 | 重力負荷状態での背骨全体の画像取得に使用 |
| Intel i7-10700 CPU搭載Windows 10システム | カスタムビルド (Various Manufacturer) | Windows 10 Professional,22H2 | モデルトレーニング用実行環境; 32GB RAM, NVIDIA GeForce RTX 3090 (24GB). |
このJoVE記事のテキストまたは図の再利用許可をリクエスト
許可をリクエスト