Method Article
セグメント化および定量化 ImageJ プラグインを使用して骨の外側の境界に基づく 2 D と 3 D 画像の海綿骨のワークフローを提案します。このアプローチはより効率的かつ現在の手動手輪郭アプローチよりも正確な現在の商用ソフトウェアでは利用できないする層によって数量を提供しています。
マイクロ コンピューター断層撮影 (マイクロ CT) は、定期的に骨量と異なる骨損失条件下で小動物の骨梁の微細構造プロパティを評価するために使用されます。マイクロ CT 画像の骨の分析のための標準的なアプローチは、スライスによって半自動手-輪郭、労働集約的、間違いやすい作業であります。骨の外側の境界、海綿骨を識別して少ないオペレーター バイアスとき適切な精度で自動的に分割によると海綿骨の自動セグメンテーションのための効率的な方法は、ここで説明しました。分割パラメーターが設定されます。分割パラメーターのすべての組み合わせが変更された 1 つずつ順番に、関連するパラメーターと領域分割結果が満足のいくセグメンテーション パラメーターをプロファイルに領域分割結果の画像のスタックが表示されます。簡単に視覚的に確認します。プラグインの品質管理機能として理論値と測定量を比較することができます, シミュレートされた標準的なオブジェクトの量を示します。海綿の性質、骨梁の厚さの層によって定量化がこのようなプラグインによって報告され、選択した領域内のような特性の分布を簡単にプロファイリングできます。措置がここで一つだけ現在の商用ソフトウェアの出力から利用できない層によって定量化海綿骨に関するより多くの情報を保持し、構造変化の統計的分析が容易にさらに、各パラメーターの定量化された値は、サンプルごとに報告されます。したがって、記載されているワークフローは、精度と効率による海綿骨の分析により良いアプローチです。
海綿骨のマイクロ CT 解析は、異なる骨損失条件1,2,3、いくつかの変数が関連の下で小動物の骨の形態学的変化を追跡するための標準的なアプローチ、骨の構造は、報告された4です。しかし、このようなパラメーターが5長骨の骨幹端に均等に配布されていない、現在商用マイクロ CT 機6,7 から各サンプルの各構造変数の要約または平均値のみが報告されます。、単一の値は完全に分析領域の測定パラメーターの特性を表すことはできません。海綿骨の層によって定量化だけでなく各変数の詳細を保持、また構造の後の統計的分析を促進する分析地域のような変数の分布のプロファイリングを有効に5さまざまな条件下での変更。したがって、このメソッドの目的は、現在、市販のマイクロ CT 解析パッケージで利用されていません、各スライス レベルでマイクロ CT スキャンの海綿骨を定量化は。
効率的にセグメント海綿骨--スライスごとに自動分割方法が望ましいです。ただし、マイクロ CT 解析のための現在の標準的な技術、労働集約的でエラーを起こしやすい、骨皮質のコンパートメントから分離する半自動補間続いて手動対話型輪郭に基づいて、実質的な演算子バイアス8,9,10に関連付けられています。自動セグメンテーション方法11,12は報告されたが、海綿骨と皮質骨の間良好な分離がある地域で地域の明確な分離なしではなく、このようなメソッドは最適です。さらに、異なるセグメンテーション パラメーターは異なるサンプル12、必須、様々 なパラメーターの組み合わせ12を試みることによって骨のサンプルのグループに適用される満足な分割パラメーターを手動で選択するは面倒です。にもかかわらず、すべての関連するパラメーターが設定されている場合、分割プロセスは自動です。骨の外側の境界がスキャンの背景を持つ最大のコントラストと長骨の骨幹端部皮質シェルは、選択した領域を分析ことができます長い骨の外側輪郭に応じてセグメンテーション方法いくつかの変更を表示確実に、正確には、皮質の殻から海綿骨を区切ります。このような分割方法の利点は、セグメンテーションは海綿骨と皮質骨6,12の違いではなく背景と骨の外側の境界の違いに基づいています。 13, したがって、それは一般的に異なるグループの間の骨梁変化のより信頼性の高い解析を促進する骨のサンプルのグループの満足できる分割パラメーターの組み合わせを見つけるは簡単。
各スライスにはボリューム、サーフェス、および 3 次元 (3 D) の厚さは 3 D 数量で報告される間、レベル、エリア、境界、および二次元 (2 D) 厚さが 2D 解析の報告されます。そのような情報は、一般的に現在画像解析ツールは、そのような情報が必要な一般のイメージに報告された手順を適用できることを示すから報告されません。
動物を対象とする手順ケアおよび実験動物の使用 (NIH パブリケーション、第 8 版、2011) ガイドに従って行われたとが検討されており動物介護制度と武漢の使用委員会によって承認されました。大学。
1. ソフトウェアのインストール
2. 海綿分析用 3 D データセットを準備します。
3. プロファイルの解析パラメーター
4. 骨梁解析
5. 模擬オブジェクトの定量化
6. 海綿等データ プレゼンテーションの校正: プロファイルの選択した分析領域の骨対策の分布
骨梁の解析プラグインは自動的にセグメント化し、精度と海綿骨を定量化する設計されています。当初、骨の外側の境界が検出され、線引き後に穴充填操作骨外側皮質シェル内の任意の穴が満ちています。その後、外側の皮質骨を除外し、セグメント化された海綿骨を得る侵食操作を実行します。最後に、海綿骨部のセグメント化された地域の対策は明らかにした.
マイクロ CT 画像は本質的にノイズの多い、区切りよく定義済みの任意のパラメーターを使用して正確に骨の外側の境界を識別するために失敗します。満足のいくセグメンテーション パラメーターを選択するためのパラメーターの組み合わせの多くをしようとする退屈な労働集約的です。骨群のパラメーターの共通セットを選択することを容易にも十分なパラメーターの組み合わせを選択する際に支援するために自動的に設定の範囲で 1 つずつパラメーターを変更するためのプラグインをプロファイリング パラメーターを指定するため、サンプル。図 1A良い分割パラメーターをプロファイルに使用設定を示しています。皮質骨のしきい値 (Cort 骨) がしきい値の範囲のパラメーター プロファイル (範囲) と各ステップでのしきい値 (ステップ) の増分量を指定、プロファイルするしきい値のシリーズが生成されます。その後、一連の音と穴の値は、対応するパラメーターを設定することにより同様に生成されます。最後に、海綿骨はすべての可能なパラメーターの組み合わせに対して一度に 1 つのパラメーターを変更することによって分割されます。図 1Bは、別のパラメーターの組み合わせの代表的な結果を示しています。明らかに、いくつかのパラメーターの組み合わせは、骨の外側の境界の輪郭を描くで他より優れている、1 つ以上のパラメーターの組み合わせが良好な領域分割結果を示しています。領域分割結果を視覚的にチェックした後、プロファイリング結果表 (表 1) から満足のいく区切りのパラメーター値を取得できます。
海綿骨の措置を定量化するには、海綿セグメンテーションと分析が実行されます。図 2Aスライスごとで選択した領域の骨が分割し、(図 2B) を抽出し、パラメーターを指定して領域分割結果をチェックできる海綿セグメンテーション用設定ダイアログが表示されます視覚的に。その後、海綿骨は十分なパラメーター (図 3A) を分析します。レポートのオプション、骨のボリューム (BV) の生数量の選択状態に応じてボリューム (テレビ)、グレー値 (強度) の合計を合計し、厚さ測定、2 次元または 3 次元 (表 2) が報告されます。最後に、校正情報はスキャンしたマイクロ CT データセットから抽出され BV、テレビ、BMC、BV/テレビの措置を校正、骨密度計算、選択した分析地域によって層の層に対しての分布をプロファイリングによって続いて位置 (図 4)。
プラグインの品質管理機能として模擬オブジェクトの定量化がサポートされています。既知寸法をシミュレートされた標準的なオブジェクトは、理論値またはマイクロ CT 装置または BoneJ14、無料のオープン ソースのプラグインに付属している商用ソフトウェアなど、ほかのソフトウェアからの対策との比較のためのプラグインにより定量化されています骨画像解析。シミュレーション画像は、ノイズのない高画質とみなされます模擬オブジェクトを分析するためノイズの設定はゼロに設定されます。厚さの異なるオブジェクトはシミュレートされたし、した結果 (図 3表 3)。円、正方形、長方形などのシミュレートされた標準 2D オブジェクトの領域 (テレビまたは BV) の値の正確な厚さが報告される (表 3)。3 D オブジェクトは、キューブ、球体、直方体に正確な厚さ測定を報告、しかし、シリンダーの厚みないシリンダーの両端、近くボクセルの正確な一方、シリンダーの中のスライスのボクセルの厚さは正確予測しました。これは、厚さ測定、基になるアルゴリズム、最大の球の直径または最大のキューブこのボクセルが含まれています、完全に内側の側面の長さによって各オブジェクトのボクセルの厚さを決定する前記の機能、オブジェクト。したがって、オブジェクトの厚さ別の球から成って、両端から中央のスライスの半径距離でシリンダーが正確に測定のみが可能、キューブを正確に測定することができます。
図 1: プロファイリング分析パラメーターの代表の結果。(A) パラメーターの設定ページします。(B) プロファイリング分析パラメーターの代表的な結果。いくつかのパラメーターの組み合わせは、骨の外側の境界を検出するため他のものより優れています。この図の拡大版を表示するのにはここをクリックしてください。
図 2: 海綿セグメンテーション分析代表結果。(A) パラメーターの設定ページします。(B) 代表的なセグメンテーション結果さまざまなレイヤーで海綿骨の。この図の拡大版を表示するのにはここをクリックしてください。
図 3: 梁解析します。(A) パラメーターの設定ページします。(B) シミュレーション 2D オブジェクトの代表的な結果。(C) 3 D ボリューム ビューアーで可視化したシミュレートされた 3 D オブジェクトの代表的な結果。この図の拡大版を表示するのにはここをクリックしてください。
図 4: 選択した分析領域における骨対策の分布.水平軸は、解析領域の開始のスライス層に相対距離を表します。Y 軸の値は、解析領域の調整海綿措置です。この図の拡大版を表示するのにはここをクリックしてください。
スライス | アウトライン境界 | ノイズ低減 φ | 穴充填 φ |
4 | 5200 | 0 | 10 |
147 | 5200 | 2 | 20 |
361 | 7200 | 6 | 16 |
539 | 8000 | 10 | 20 |
表 1: プロファイリング分析パラメーターの代表的な結果。
表 2: 海綿分析の代表的な結果。
オブジェクト | ディメンション | ボリュームb | 表面c | 厚さ |
広場 | 200 X 200 | 40000 | 796 | 200 |
四角形 | 200 X 100 | 20000 | 596 | 100 |
円 | Dia: 200 | 31428 | 796 | 200 |
キューブ | 30 X 30 X 30 | 27000 | 5048 | 30 |
立方体 | 80 X 40 X 30 | 96000 | 13008 | 30 |
球 | Dia: 30 | 14328 | 3944 | 30 |
シリンダーd | Dia:30;H: 100 | 71600 | 12800 | 27.84 |
シリンダーe | Dia:30;H: 100 | 51552 | 9552 | 30 |
a: 結果は、生のボクセルです。Dia: 直径;H: 高さ。 | ||||
b: ボリューム (3 D) または領域 (2 D)。 | ||||
c: 表面 (3 D) や境界 (2 D)。 | ||||
d: スライス 100 にスライス 1 からデータセルは解析に使用されます。 | ||||
e: スライス 15 から 85 のスライスにデータセルは解析に使用されます。 |
表 3: 数量は、シミュレートされたオブジェクトを結果します。
の補足ファイル 1。骨をサンプルします。このファイルをダウンロードするここをクリックしてください。
本研究では、海綿骨部、自動、効率的であるとユーザー ・ フレンドリーを分析するための ImageJ プラグインについて説明します。任意の 2D または 3D オブジェクト領域、ボリューム、および厚さの層によって対策を定量化するプラグインが使えます。現在、骨梁パラメーターごとに単一の測定値のみは、各サンプルの選択した分析領域で測定されたエンティティの特性を表すことはできません完全に標準のマイクロ CT 解析によって報告されます。記載されているプラグイン レポートの各サンプルでは、選択した解析領域、したがってより高度なと敏感な統計的測定のパラメーターの分布情報を完全に保持する各パラメーターの層によって数量アプローチは、このようなデータの解析に適用されます。
演算子バイアス8,9,10より均一を正当化標準、半自動手輪郭セグメンテーション法、異なる手輪郭スライスの間実質的な変動がみ海綿骨領域分割の手法。プラグインは、精度と適切な骨の外側の境界を自動的に検出することができますどこで外側境界アシスト海綿分割を行い、同様にコントラストは最大スキャン背景と骨の外側の境界の間、しきい値、ノイズ、および穴の設定が提供されます。分割パラメーターの決定を容易にする一連のパラメーターの組み合わせは、プロファイリング分析地域から代表的なスライスを使用して 1 つずつ、満足のいく結果が視覚的にチェックされます。その後、海綿骨はプロファイルのパラメーターを分析します。同じ画像の再現性のある結果が報告されます適切な分割パラメーターが設定されていると、プラグインが正常に異なる骨同化剤5と扱われるラット遠位大腿骨から骨を分析に適用されています。オペレーターにとらわれない。
プラグイン現在のみ処理できますサンプルを 1 つずつ。複数のサンプルがイメージに存在する場合、すべてのサンプルは差別することがなく単一のオブジェクトとして定量化しました。したがって、複数のサンプルを同時に単一のサンプル チューブにスキャンする場合、前処理が必要です。定量出力は、ピクセル数または灰色の値は、適切なキャリブレーション情報を使用して手動で調整する必要があります生対策です。
このプロトコルでは最も重要なステップは、海綿の分析のための適切なパラメーターを選択します。一般的に、骨のサンプルのグループは結果に匹敵する同一分割パラメーターを使用して分析されます。このような場合は、注意が必要、セグメンテーション結果すべてのサンプルの調査視覚的に潜在的なエラーの解析を実施してさらに前にことを確認します。
この手法の主な制限は、ノイズリダクションの設定が皮質骨抽出した海綿骨ではなく骨の外側の境界の輪郭を描くためだけに適用されることです。さまざまなフィルタ リング方法は、異なる条件下で骨梁の解析に最適なとして報告された15をされている、1 つのフィルターはすべてイメージ サンプルに適してないです。したがって、海綿骨の定量化する前に組み込み濾過手順はありません。組み込みイメージのろ過ソフトウェアによって出力に対応するには、抽出された海綿骨加工様々 なイメージング ツールを使用して同様に最初にすることができます ImageJ にインポートし、その後プラグインによって量を示されます。また、さまざまなイメージ形式が異なるマイクロ CT ベンダーによって生成されるとそれらのほとんどをインポートまたは ImageJ で直接開くことはできません、したがって ImageJ にさらなる分析のためのさまざまなイメージ形式をインポートするためのモジュールが必要があります。別の制限は、選択範囲の両端近くのスライスで 3 D の厚さは正確なオブジェクトの一部が選択範囲、したがって、任意の利用可能な 3 D 定量パッケージにも当てはまりますスライスのより大きい範囲外にいる場合は定量化する必要があり、選択したスライスの中央部分だけはさらに分析に使用する必要があります。このようなアプローチはすることができます測定量は、このプラグインではなく、測定選択範囲の定量的パラメーターごとに 1 つの値だけが報告されることを示す層によって報告されたときにのみ使用できます。
オブジェクトの厚さの定量化の基になるアルゴリズムの特許が提出されています。著者と共同で他の人は bomomics.com ウェブサイトのコンサルティングと定量化画像解析サービスをリクエストに応じて提供する、無料版プラグインをホスティングを登録します。
この作品はグラント NFSC 81170806 によって部分的に支えられました。著者は、スキャンやラットの大腿骨を分析支援のための武漢大学口腔病学のマイクロ CT コア施設に感謝したいと思います。
Name | Company | Catalog Number | Comments |
ImageJ | NIH | imagej | Any version with a java 1.8 run time |
trabecular analysis plugin | Bomomics | bomomics | free or commercial version |
Micro CT scanner | Scanco | μ-50 | micro CT from any vendor |
Computer System | Lenovo | any brand | |
Windows Operating System | Microsoft | Windows 7 x64 | any 64-bit Windows operating system |
Office Software | Microsoft | Office 2010 | any speadsheet software that has xy chart function |
このJoVE論文のテキスト又は図を再利用するための許可を申請します
許可を申請This article has been published
Video Coming Soon
Copyright © 2023 MyJoVE Corporation. All rights reserved