クリティカルサイズの骨欠損におけるユーザー定義の関心領域に対する縦断的マイクロコンピュータ断層撮影画像解析

クリティカルサイズの骨欠損は、整形外科治療管理において重大な臨床的課題をもたらします。ASTM F2721によると、臨界サイズの欠陥は、対象骨の直径1.5〜2倍の長さの欠陥として特徴付けられます¹。これらの欠陥の修復は、伝統的に、手続き費用、二次手術の関連リスク、および必要な骨移植量によって制限された自家移植および同種移植の使用を通じて行われてきました²。現在の骨再生技術は、機械的特性、生体適合性、生物活性、血管新生能、および^{分解プロファイル3,4,5}を最適化することにより、骨伝導性と骨誘導性の両方の効果を生み出すように設計された同種および異種性の足場の使用に焦点を当てています。バイオマテリアルは、バイオセラミックスやバイオポリマーから金属などの複合材料まで幅広く研究されています⁶。これらの生体材料のバリエーションは、骨再生の足場としての可能性を調べるために、in vitroとin vivoの両方で試験されています。

μCTは、げっ歯類モデル^7,8,9の骨の形態、構造、微細構造を評価するための非侵襲的で高精度なイメージングのゴールドスタンダードです。このイメージングモダリティは、骨折治癒モデル¹⁰における骨再生の縦断的、生体内進行を評価するために記載されている。μCT^{スキャン9}からの皮質骨および海綿骨の定量化を標準化するための方法が開発されました。半自動セグメンテーションワークフローは、複雑な解剖学的構造を持つ全骨セグメンテーションのための市販の視覚化ソフトウェアを利用して開発されました¹¹。これらの方法により、さまざまな経験レベルのユーザーが、標準化された再現性のある結果を生成するための、シンプルで親しみやすい方法が可能になります。ただし、これらの方法では、ユーザー定義のROIを調査する能力は依然として限られています。

ここでは、可視化ソフトウェアを使用して、縦断的な in vivo ラットモデルの臨界サイズの骨欠損を取り巻くユーザー定義のROI骨量分析を可能にするために、現在の方法を拡張したプロトコルを紹介します。縦断的研究の数週間間で一貫したアライメントとROI選択方法を確立することは、堅牢なプロトコルの開発に不可欠でした。初期時点は、ソリッドモデルの一貫した向きを確保するために、後続の週のアライメントのベースラインとして使用されます。このアライメントにより、重ね合わせたソリッドモデルから対応するμCT画像スライスを選択し、臨界サイズの欠陥を包含することができます。一貫したROIは、スライスの位置だけでなく、領域内のスライス数の比較によっても検証されます。ベースラインモデルから選択したROIは、その後の数週間で再現できるため、比較可能な定量的分析が可能になります。

この研究の縦方向の μCT 画像は、ポリカプロラクトン (PCL) ベースの足場で治療された成体雌チャールズ リバー SASCO-SD ラットの 3 mm の重篤な橈骨欠損から 0 週目と 6 週目に収集されました。すべての動物の使用は、ロチェスター大学の動物資源委員会(UCAR)によって承認されたプロトコルに従って行われました。μCT画像収集は、Scanco Medical VivaCT 40を使用して行いました。

注:このプロトコルの主要なステップは、μCT画像のセグメンテーション、アライメント、ROIの選択とクロッピング、および分析と視覚化に分かれています(図1)。μCT画像セグメンテーションのプロトコルは、Kenney et al. (2022)¹¹から採用されています。

図 1: 概要ワークフロー図。 プロトコルのステップは、主に画像のセグメンテーション、モデルのアライメント、ROIの選択、およびボリューム分析に分かれています。 この図の拡大版を表示するには、ここをクリックしてください。

1. μCT画像セグメンテーション

1. ソフトウェア内で画像を開く
   1. Amiraソフトウェア(以下、可視化ソフトウェア)を起動し、[ Open Data]を選択します。ポップアップウィンドウで、目的のフォルダディレクトリに移動して、ファイルの初期時点.dcmを選択します。 CTRL+A を押してすべての画像を選択し、[ 画像を開く]をクリックします。
   2. ポップアップウィンドウには、開いている画像データセットに関する情報が表示されます。 [OK ] をクリックして続行します。
   3. スケーリングに関する警告が表示されます。「 Convert to Float 」(デフォルトオプション)をクリックして続行します。
   4. データセットのアイコンが 「Project View 」セクションに表示されます。データ・セットの名前を変更するには、データ・セット・アイコンを 1 回クリックして F2 キーを押します。目的の名前を入力したら、[ OK]をクリックします。スキャンを識別するための命名パターンを確立します (例: 'Bone Scaffold Material' - 'WKX')。
2. μCT画像のフィルタリングと閾値処理
   1. [プロジェクト ビュー] セクションに画像を読み込んだ状態で、[オルソ スライス] のデフォルトの 2D ビューアをクリックして、2D の向きとスライスを表示します。2D ビジュアライゼーションが右側の表示ウィンドウに表示されます。
   2. 画像を 3D で視覚化するには、データ セットを右クリックし、ボリューム レンダリングを検索して選択します。ボリュームレンダリングの[プロパティ]セクションで、[カラーマップ]の下限しきい値を2500に調整し、Enterキーを押します。
      注: 2D と 3D は、それぞれ [Ortho Slice] または [Volume Rendering] の横にある青い四角をクリックすることでオンとオフを切り替えることができます。右側の表示ウィンドウでは、バーの上部に沿って複数の表示オプションを使用できます。これらには、画像を移動、回転、またはズームイン/ズームアウトするためのさまざまなツールが含まれます。画像の向きを XYZ、XZ、YZ ビューに変更するためのツールと、簡単な測定のためのツール。
   3. 必要に応じて、一部のデータ セットでは、 ボリューム編集を使用して対象領域外のアーティファクトまたはデータを削除する必要がある場合があります。 Volume Editを追加するには、データセットを右クリックし、 Volume Editを検索して選択します。
   4. Project View セクションの Volume Edit をクリックします。[プロパティ] セクションで、[ツール] の右側にある最初のドロップダウン ウィンドウから [TabBox] を選択します。
   5. 保持する領域を選択するときは、Interactカーソル(表示ウィンドウの上部バーにある尖ったカーソルアイコン)を使用して緑色の角をクリックしてドラッグし、 TabBox を調整します。 TabBox を調整するには、別のビュー (XY、XZ、YZ) を使用して、削除する領域を特定する必要がある場合があります。
   6. 保持する領域が TabBox で囲まれたら、[ 外側を切り取 る] を選択して、 TabBox の外側にあるアーティファクトまたはデータを削除します。新しい変更されたデータセットがプロジェクトビューに生成されます。
   7. ボリューム レンダリング文字列をクリックして、元のデータ セットから変更したデータ セットにドラッグします。これにより、表示ウィンドウに表示される画像が変更されたデータセットに変更されます。ボリュームレンダリングの[プロパティ]セクションで、[カラーマップ]の下限しきい値を2500に調整し、Enterキーを押します。
   8. 変更したデータを右クリックし、 Median Filterを検索して選択します。 [Median Filter] の [Properties] セクションで、[Interpretation] に [3D ] を選択し、[ Apply] をクリックします。これにより、新しいフィルタリングされた(.filtered)データセットが プロジェクトビュー セクションに作成されます。
   9. フィルタリングされたデータを右クリックし、Interactive Thresholdingを検索して、を選択します。Interactive Thresholding の Properties セクションで、Preview Type に 3D を選択し、Intensity Range の下限値を 2500 に調整して、Apply をクリックします。これにより、Project Viewセクションに新しいしきい値設定(.thresholded)データセットが作成されます。
   10. [Volume Rendering] 文字列をクリックして、変更したデータ セットからしきい値が設定されたデータ セットにドラッグします。Interactive Thresholding をオフにし、Volume Rending をオンにします (オフの場合)。
3. μCT画像のセグメンテーション
   1. フィルタリングされたデータセットをクリックし、メニューバーの下にある「セグメンテーション」をクリックして「セグメンテーション」タブに変更します。「セグメンテーション」タブ内のコントロールは、「プロジェクト」タブのコントロールと似ています。
      注: [セグメンテーション ] タブには、表示ウィンドウに関連する追加のツールがいくつかあります。表示ウィンドウの上隅には、画像が表示ウィンドウを埋めるために使用できる +/-ズーム ボタンがあります。画像は、表示ウィンドウの下部にある左/右のスクロールバーを使用して、表示ウィンドウの中央に配置することができます。このスクロールバーの上には、μCTデータスライス内を移動するための2番目のスクロールバーがあります。
   2. 「 セグメンテーション 」タブが適切に設定されていることを確認するために、「 セグメンテーションエディタ 」セクションで、「 画像 」をフィルタリングされたデータセットに設定し、「 ラベル 」フィールドをしきい値設定されたデータセットに設定します。 [Materials ] セクションで [ Material ] という名前のマテリアルをダブルクリックし、名前を 「Bone」に変更します。
   3. セグメンテーションエディタの下の「新規」をクリックして、新しいラベルフィールドを作成します。ポップアップウィンドウで、前と同様の命名規則 (例: 'Bone Scaffold Material' - 'WKX' - Markers) を使用して名前を Markers に変更し、OK をクリックします。
   4. 代替ビュー (XYZ、XZ、YZ) を使用して 2D 表示ウィンドウを調整し、画像の矢状面ビューまたは冠状平面ビューを表示します。
   5. 選択セクション内でブラシ、なげなわ、ワンド、しきい値を使用して、関心のある領域を特定します。このプロトコルでは、ブラシツールを選択します。ブラシツールのサイズを調整します。一般に、中サイズのポインターが適しています。
   6. 画像スタックを、半径ボーンが見え始める最初のスライスに調整します。これには、スライスを複数回移動して、橈骨と尺骨に属するボリュームと、それらがいつ開始するかを判断する必要がある場合があります。
   7. スライスが特定されたら、 ブラシカーソル を使用して、このボーンのセグメントの周りのカーソルをクリックしてドラッグし、描画します。カーソルが引かれた場所に赤い線が表示されます。骨の外側ではなく、骨のセクションの内側に描きます。これにより、視覚化ソフトウェアが 2 つのボーンを区別するのに役立ちます。
   8. スライスを移動し続け、半径の周りを描画します。通常、画像スタックの中央部分を通じて、20〜30枚の画像スライスごとに描画を作成できます。画像スタックの端に沿って、または臨界サイズの欠陥を囲むように、5〜10枚の画像スライスごとに描画を行い、骨がフェードインまたはフェードアウトするときに視覚化ソフトウェアを支援します。
   9. 半径の全長が完成したら、メニューバーのドロップダウンで [選択]、[塗りつぶし]、[ すべてのスライス ]の順にクリックします。これにより、ブラシツールで作成されたすべての描画が塗りつぶされます。これらは、アウトラインではなく赤で網掛けされて表示されるようになりました。
   10. メニュー バーのドロップダウンで [Selection and Interpolate] をクリックします。これにより、作成された描画に基づいて、すべてのスライスの半径ボーンの影付きの領域が埋められます。
       注:これにより、骨の一部が欠落することがあります。続行する前に、すべてのスライスをスクロールして、ボーンが一般的に識別されたことを確認します。ボーンの一部が欠落している場合は、ブラシツールを使用して、欠落しているセクションを囲むいくつかのスライスでこの領域を一周し、上記の2つの手順(1.3.9と1.3.10)を繰り返します。
   11. 半径が完全に特定されたら、[ Materials ] セクションで [ Inside ] という名前のマテリアルをダブルクリックし、名前を 「Radius」に変更します。
   12. [選択] セクションの下にある [追加] (+) 記号をクリックして、これをマテリアルとして追加します。これで、半径が表示ウィンドウに輪郭が引かれます。先に進む前に、[マテリアル]セクションの[半径]マテリアルの[ロック]アイコンをクリックします。これにより、この資料に変更が加えられることはありません。
   13. 尺骨についても同じプロセスを完了します。 [マテリアル ] セクションで [ 追加 ] をクリックして、新しいマテリアルを作成します。新しいマテリアルをダブルクリックし、名前を 「Ulna」に変更します。
   14. 尺骨について、手順1.3.6から1.3.12を繰り返します。
   15. Radius ボーンと Ulna ボーンの両方の手順が完了したら、[Materials] セクションの [Lock] アイコンをクリックして [Radius] のロックを解除します。メニューバーの下にある「プロジェクト」をクリックして、「プロジェクト」タブに変更します。これで、Markers という名前のデータ セットが [Project View] セクションのフィルター済みデータ セットに適用されました。
   16. フィルタリングされたデータ セットを右クリックし、「 Marker Based Watershed Inside Mark (Image Segmentation)」 を検索して選択します。[ プロパティ ] セクションで、ドロップダウン メニューをクリックして、[ データ ] をフィルター処理されたデータ セットとして、[マーカー] データ セットとして、[マーカー] データ セットとして、 バイナリ マスク をしきい値データ セットとして設定し、[ 適用] をクリックします。
   17. .grown ファイルは [Project View ] セクション内に生成されます。しきい値が設定されたデータ セットの [Volume Rendering] 文字列をクリックして .grown データ セットにドラッグし、表示ウィンドウ内の橈骨と尺骨のセグメンテーションを視覚化します。
   18. .grown ファイルを 8 ビットに変換するには、.grown データ セットをクリックし、[ Convert Image Type ] を検索して選択します。.grown ファイルは 16 ビットから 8 ビットに変換されます。

2. アライメント

1. 軸のアライメントと橈骨と尺骨の抽出
   注: 現在のモデルが初期時点の場合は、すべての手順を完了する必要があります。それ以外の場合は、このセクションをステップ 2.1.4 から開始します。半径の横方向のスライスを可能にするには、初期時間点モデルをオルソスライス平面に対して垂直に位置合わせする必要があります。このモデルは、その後の週のアライメントのベースラインとして機能します。
   1. [Project View] セクションで、[Ortho Slice] をオンにし (オフの場合)、.grown データセットをクリックします。[プロパティ] セクションで、[トランスフォーム エディター] アイコンをクリックします。Interactカーソルを使用して、緑色の軸の位置合わせポイントをクリックして半径と尺骨の角度を調整し、Ortho Sliceが半径を通る横方向のスライスを作成するようにします。
   2. これを整合したら、この変換を保存するために新しいデータ セットを作成する必要があります。.grown データセットを右クリックし、 Resample Transformed Image を検索して選択します。
   3. [Resample Transformed Image] をクリックし、[Properties] セクションで [Data] を .grown データ セットに、[Interpolation] を [Nearest Neighbor] に、[Mode] を [extended] に、[Preserve (保存)] を [Voxel Size] に、[Padding (パディング)] の値を 0 に設定します。[適用] をクリックすると、新しい .transformed データ セットが作成されます。
   4. 結合されたセグメンテーションモデルから半径と尺骨を抽出するには、.transformed データセットを右クリックし、 Extract Label を検索して選択します。
   5. [ Extract Label] をクリックし、[ Properties ] セクションで [ Labels ] を .transformed データセット に、[ Label ID ] を 1 に、[ Export to Binary] をオンにします。 [適用] をクリックすると、結果データ セットが作成されます。
   6. ボリュームレンダリング文字列をクリックして、.growedデータセットから結果データセットにドラッグします。Result データセットをクリックし、F2 キーを押してファイルの名前を変更します (例: 'Bone Scaffold Material' - 'WKX' - Radius)。
   7. 尺骨に対して手順 2.1.4 から 2.1.6 を繰り返して、 ラベル ID を 2 に設定し、尺骨の結果ファイルの名前を変更します。
2. 縦方向の画像を保存するファイルを開く
   注: 現在のモデルが初期時点の場合は、すべての手順を完了する必要があります。それ以外の場合は、このセクションをステップ 2.2.2 から開始します。
   1. 初期時点ファイルを保存するには、メニュー バーのドロップダウンで [File]、[Save Project As]、および [Set File Directory Location ] をクリックします。保存タイプをAmiraプロジェクトおよびデータファイル(*.hx)として設定し、ファイルにマスターバージョンとして名前を付けます(例:'Bone Scaffold Material'_'WKX' MASTER)。
   2. 初期時点のマスター ファイルを開き (まだ開いていない場合)、メニュー バーのドロップダウンで [File]、[Save Project As]、 および [Set File Directory Location ] をクリックしてコピーを保存します。保存タイプを Amira プロジェクトおよびデータ ファイル (pack &go) (*.hx) ファイルとして設定し、ファイル名 (例: 'Bone Scaffold Material'_ 'WK0' と 'WKX') を付けます。このファイルは、マスターファイルを上書きせずに、最初の時点とその後の週を比較するために使用されます。
   3. 新しい比較ファイル内で、[Project View] セクションの [Open Data] をクリックして、比較時点の画像を開きます。ポップアップウィンドウで、目的のフォルダディレクトリに移動して、ファイルの初期時点.dcmを選択します。CTRL+Aを押してすべての画像を選択し、[画像を開く]をクリックします。
   4. ポップアップウィンドウには、開いている画像データセットに関する情報が表示されます。 [OK ] をクリックして続行します。
   5. スケーリングに関する警告が表示されます。「 Convert to Float 」(デフォルトオプション)をクリックして続行します。
   6. データセットのアイコンが 「Project View 」セクションに表示されます。データ・セットの名前を変更するには、データ・セット・アイコンを 1 回クリックして F2 キーを押します。目的の名前を入力したら、[ OK]をクリックします。スキャンを識別するための命名パターンを確立できます (例: 'Bone Scaffold Material' - 'WKX')。
3. モデルの位置合わせ
   注:このプロセスでは、橈骨のセグメンテーションのアライメントについて説明します。必要に応じて、尺骨にも同じプロセスを適用できます。両方のボーンを同時に位置合わせしないでください。これにより、適切な位置合わせに問題が発生します。
   1. 初期時点の ボリュームレンダリングのいずれかをオフにします。 表示ウィンドウ内には何も表示されません。
   2. 比較 (後続の週) データ セットを右クリックし、Ortho Slice を検索して選択します。
   3. セクション 1 (ステップ 1.2 から開始)、セクション 2 (ステップ 2.1) を繰り返し、ステップ番号に関する追加の注意事項を確認します。これらのセクションが完了したら、ステップ 2.3.4 から再開します。これらのセクションを完了すると、翌週のデータセットのセグメンテーションが生成され、橈骨と尺骨が抽出されます。
   4. 抽出された橈骨の比較(後続の週)データセットを右クリックし、 Image Registration Wizardを検索して選択します。 [プロパティ ] セクションで、抽出された半径ボーンの比較週データ セットとして [データ ] を設定し、抽出された半径ボーンの初期時点データ セットとして [参照 ] を設定します。
   5. [Image Registration Wizard Actions] セクションで、[Step 1 of 4] の [Skip] をクリックします。ステップ 2 と 3/4 では、Interact カーソルを使用して、TabBox を初期時点と比較週のデータ セットの間の共通の領域に調整し、各ステップの後の [アクション] の下の [適用] をクリックします。ステップ 4/4 で、メトリックを相関に、変換をリジッドに、事前アライメントを主軸の整列として設定し、アクションの下の適用をクリックします。
   6. データセットが整列したら、抽出された半径ボーンの比較(翌週)データセットを右クリックし、Resample Transformed Imageを検索して選択します。[プロパティ] セクションで、抽出された半径ボーンの比較週データ セットとして [データ] を設定し、[Interpolation to Nearest Neighbor] を [Extended] に、[Mode] を [Extended] に、[Preserve to Voxel Size] に、[Padding Value] を 0 に設定し、[Apply] をクリックします。比較データセットに対して、新しい .transformed データセットが作成されます。

3. ROIの選択とクロップ

注:最初にROIクロップを完了して、クリティカルサイズのフラクチャーの周りのスライス番号を決定します。これらのスライス番号が決定されたら、必要に応じて尺骨の同じ時点でこれらを使用することができます。このプロセスでは、モデルに一度適用したトリミングを元に戻すことはできません。

1. 対象作物の領域
   1. をクリックして、初期時点のオルソスライスをオンにし、抽出された半径の初期時点データセットにデータを設定します。平面が半径ボーンを横切るように Orientation を設定します。
   2. [Properties] セクションの [Slice Number] スライドバーを使用して、スライス番号を変更し、クリティカルサイズの欠陥を囲む近位および遠位のスライス位置を決定します。両方のスライスについて、骨折が橈骨の骨幹と出会う最も遠位および近位のスライスを決定します。スライス番号を文書化します。この番号はモデルごとに異なります。
   3. 比較週の [オルソ スライス] をクリックしてオンにし、[ データ ] を抽出された半径の初期 時点データ セット に設定します。平面が半径ボーンを横方向にカットするように、向きを設定します。
   4. [プロパティ] セクションの [スライス番号] スライド バーを使用し、遠位オルソ スライスを示す初期タイム ポイント データ セットを使用して、比較週スライス番号を初期タイム ポイントの遠位スライスに揃えように変更します。視覚化ソフトウェアでは、位置合わせを視覚的にチェックできます。スライスが重なっている場合、表示ウィンドウには 1 つのスライスとして表示されます。
   5. 比較週データセットの遠位スライスのスライス番号に注意してください。近位スライスに対して手順 3.1.4 を繰り返します。
      注:追加のチェックとして、初期週と比較週のデータセットの両方について、近位スライスと遠位スライスの違いを決定します。これらの値は、一貫したサイズの ROI を示すために、同等である必要があります。
   6. 抽出された半径の初期時点をクリックします。 [プロパティ ]セクションで、[ クロップエディター ]ツールをクリックします。
   7. クロップ エディターのポップアップ内で、適切な X、Y、または Z フィールドに最小値と最大値を入力します (これらが入力されると、表示ウィンドウによって ROI が変更されます)。値が入力されたら、[OK] をクリックします。データセットは入力ROIにトリミングされます。
   8. 比較週データセットについて、ステップ 3.1.7 を繰り返します。

4. 分析と可視化

1. ボリューム分析
   1. 初期タイム ポイント データ セットのボリュームを確認するには、抽出された半径の変換された初期タイム ポイント データ セットを右クリックし、[ マテリアル統計]を検索して を選択します。 [プロパティ ] セクションで、抽出された半径の変換された初期時点として [データ ] を設定し、[ 選択] を [マテリアル ] に設定して [適用] をクリックします。新しい .MaterialStatistics データ セットは、初期タイム ポイント データ セットに対して作成されます。
   2. をクリックします。抽出された半径の初期時点のMaterialStatisticsデータセットと、プロパティウィンドウ内でスプレッドシートの表示をクリックします。プロパティウィンドウの上にあるテーブルタブをクリックします。このタブ内には、抽出された半径に設定されたクロップされた初期タイムポイントデータセットのボリュームが表示されます。
   3. 抽出された半径の比較週データセットについて、手順4.1.1と4.1.2を繰り返します。 [テーブル ] タブ内では、最初の時点と比較週に個別のタブがあり、クリックで移動できます。
2. データセットの可視化
   1. ボーンの体積の変化を視覚化するには、抽出された半径の比較週変換データセット を右クリックし 、 Arithmeticを検索して選択します。 [算術演算 ] をクリックし、[ プロパティ ] ウィンドウで、抽出された半径の比較週変換データ セットとして [入力 A ] を設定し、抽出された半径の [ 入力 B ] を初期時点データ セットとして設定し、[ 入力 C ] を [ソースなし] として設定し、[ 結果タイプ] を [入力 A] として設定し、[ オプション ] をオフのままにして、[ 結果チャネル] を入力 A と同様に設定します。 をクリックし 、[Expression] を [A-B] に設定します。
   2. 結果データセットをクリックし、F2 キーを押してファイルの名前を変更します (例: 'Bone Scaffold Material' - 'WKX' - Bone Change)。この結果データセットを右クリックし、[Generate Surface] を検索して選択します。Generate Surface をクリックし、Properties ウィンドウで Apply をクリックし、ポップアップウィンドウで Continue をクリックします。新しい .surf データセットが作成されます。
   3. .surf データセットを右クリックし、 Surface View を検索して選択します。算術結果のサーフェス ビューが [表示] ウィンドウに表示されます。
   4. Surface Viewの色を変更するには、Project ViewウィンドウでSurfaceビューをクリックし、次にPropertiesウィンドウでColorsのドロップダウンをクリックしてConstantを選択します。カラーマップをクリックして、好みの色を割り当てます。
   5. 最初の週のデータ セットのボーン体積の変化を表示するには、.transformed データ セットを右クリックし、 Extract Label を検索して選択します。[ Extract Label] をクリックし、[ Properties ] セクションで [ Labels ] を .transformed データセット に、[ Label ID ] を 1 に、[ Export to Binary] をオンにします。 [適用] をクリックすると、 結果 データ セットが作成されます。 Result データセットをクリックし、 F2 キーを押してファイルの名前を変更します (例: 'Bone Scaffold Material' - 'WKX' - Full - Radius)。
   6. この結果データセットを右クリックし、[ Generate Surface] を検索して選択します。 [Generate Surface] をクリックし、[ Properties ] ウィンドウで [Apply ] をクリックし、ポップアップ ウィンドウで [Continue] をクリックします。新しい .surf データセットが作成されます。
   7. .surf データセットを右クリックし、 Surface View を検索して選択します。算術結果のサーフェス ビューが [表示] ウィンドウに表示されます。
   8. Surface Viewの色を変更するには、Project ViewウィンドウでSurfaceビューをクリックし、次にPropertiesウィンドウでColorsのドロップダウンをクリックしてConstantを選択します。カラーマップをクリックして、好みの色を割り当てます。
   9. さらに数週間続行するには、ステップ 2.2.2 に戻ります。

μ3つのユニークなラットモデルのCT画像を調査し、それぞれがポリカプロラクトン(PCL)足場で処理し、肯定的な結果を示しました。時点をまたいで縦断的研究を分析するには、ROIを選択してクロップする前に、収集されたソリッドモデルを位置合わせする必要があります。この能力を数週間にわたって示すために、0週目と6週目に収集されたソリッドモデルを共通の領域を使用して位置合わせしました(図2)。

図2:週ごとのソリッドモデルのアライメント。 イメージレジストレーションウィザードを使用した、週0をベースライン方向として使用した週間の代表的なソリッドモデルのアライメント。この登録プロセスは半自動化されており、2つのモデル間にユーザー定義の共通領域(点線)を適用して、ソリッドモデルの自動位置合わせを支援します。 この図の拡大版を表示するには、ここをクリックしてください。

画像の位置合わせが成功し、ユーザー定義の ROI が決定された後、0 週目のトリミングされたソリッド モデルが 6 週目のトリミングされたソリッド モデルから差し引かれました。これにより、ROI内の骨量が変化しました(図3)。

図3:骨の体積変化を決定するための算術演算。代表的な算術法は、2週間間の骨量を差し引くことにより、骨量の変化を決定するために使用されます。ROIの骨量の変化は、縦断的研究を通じて傾向を示すことができます。 この図の拡大版を表示するには、ここをクリックしてください。

次に、骨量の変化を使用して、0週目のベースラインからのこれらの変化を視覚的および定量的に示しました。骨の変化のソリッドモデルを 0 週目のベースライン骨ソリッドモデルに投影して、骨量の変化を視覚的に示しました (図 4)。

図4:1人のユーザーの重大な欠損時の骨容積変化の代表的な画像。 PCLスキャフォールドのSurface Viewモデルの視覚化は、経験豊富なユーザー向けに再現されます。各ユーザーは、3つのユニークなラットモデル(PCLスキャフォールド)を解析し、それぞれで3回の試行を行いました。0 週目のボーンは灰色で表示され、0 週目と 6 週目のボーンの変化は青色で表示されます。 この図の拡大版を表示するには、ここをクリックしてください。

さらに、0 週目から 6 週目までのこれらの骨量の変化は、各 PCL 足場について三重に定量化されました。2人の経験豊富なユーザーと1人の初心者ユーザーが分析を実行しました。各ユーザー(n = 3)は、0週目と6週目の時点でこれらのモデルを3回に分けて分析し、重篤な欠損領域におけるこれらの週間の骨量の変化を決定しました(図5)。PCL2サンプルでは、経験豊富なユーザー1と2の間で、初心者ユーザーと比較して統計的な有意性が観察されました(それぞれ0.0485と0.0235のp値)。他のすべての差は統計的に有意ではなかった。各PCLスキャフォールドについて得られたユーザー値を 表1 に示します(データは平均±SDで表されます)。各PCLスキャフォールドのユーザー平均を使用して、二元配置混合モデルと一貫性タイプ解析を使用して、クラス内相関係数(ICC)を実行しました。計算されたICCは0.997で、p値は4.45×10-05であり、優れた信頼性を示しています。

図5:代表的な骨量変化結果。 0 週目から 6 週目までの骨量の変化は、3 人のユーザーによって 3 回に取得された各 PCL 足場について計算され (n = 3 ユーザーあたりの足場あたり、エラーバーは平均 SD を示します)、二元配置 ANOVA とテューキーの多重比較検定 (nsP > 0.05、*P < 0.05) を使用して統計的有意性を分析しました。 この図の拡大版を表示するには、ここをクリックしてください。

表 1: PCL グループ別のユーザー測定の概要。 各 PCL グループのユーザー結果は、平均 ± SD として示されます (n=3/ユーザーあたりスキャフォールドあたり)。

著者には、開示すべき利益相反はありません。