JP2021500113A - 椎間板モデリング - Google Patents
椎間板モデリング Download PDFInfo
- Publication number
- JP2021500113A JP2021500113A JP2020521914A JP2020521914A JP2021500113A JP 2021500113 A JP2021500113 A JP 2021500113A JP 2020521914 A JP2020521914 A JP 2020521914A JP 2020521914 A JP2020521914 A JP 2020521914A JP 2021500113 A JP2021500113 A JP 2021500113A
- Authority
- JP
- Japan
- Prior art keywords
- vertebral
- dimensional
- neural network
- residual neural
- medical image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4538—Evaluating a particular part of the muscoloskeletal system or a particular medical condition
- A61B5/4566—Evaluating the spine
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
- G06T2207/30012—Spine; Backbone
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Rheumatology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Physical Education & Sports Medicine (AREA)
- Quality & Reliability (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- High Energy & Nuclear Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Prostheses (AREA)
- Image Analysis (AREA)
Abstract
脊椎分析セグメンテーションのための他の方法が開示される。一例では、当該方法は、完全畳み込みネットワークと残差ニューラルネットワークとを結合することを含む。当該方法はまた、残差ニューラルネットワークと共に結合された完全畳み込みネットワークを端から端まで訓練することを含む。当該方法はまた、脊椎分析の少なくとも1つの医用画像を受信することを含む。当該方法はまた、残差ニューラルネットワークと共に完全畳み込みネットワークを少なくとも1つの医用画像に適用し、脊椎分析の少なくとも1つの医用画像から少なくとも1つの椎体をセグメンテーションすることを含む。
Description
[関連出願への相互参照]
この出願は非仮出願であり、2017年10月20日に出願された米国仮出願第62/575,260号に対する優先権を主張し、その全開示を参照により援用する。
この出願は非仮出願であり、2017年10月20日に出願された米国仮出願第62/575,260号に対する優先権を主張し、その全開示を参照により援用する。
[背景技術]
この開示は、椎間板(intervertebral disc)をモデリングする方法に関する。医用画像からの椎骨(vertebrae)の正確なセグメンテーションは、画像支援脊椎手術における重要な準備作業である。椎骨の詳細な形状を記述することは、多数の脊椎疾患の早期診断、手術計画及び経過観察の評価にかなり役立つ可能性がある。近年、磁気共鳴画像法(MRI, magnetic resonance imaging)は、椎骨解析のための価値ある非侵襲的ツールになっている。コンピュータ断層撮影法(CT, computed tomography)とは異なり、MRIは放射線被曝に関連するリスクを引き起こさない。しかし、磁気共鳴画像解析は、異なるスキャンにわたる解像度及び強度範囲における変動のような、セグメンテーションのためのいくつかの課題に直面する。
この開示は、椎間板(intervertebral disc)をモデリングする方法に関する。医用画像からの椎骨(vertebrae)の正確なセグメンテーションは、画像支援脊椎手術における重要な準備作業である。椎骨の詳細な形状を記述することは、多数の脊椎疾患の早期診断、手術計画及び経過観察の評価にかなり役立つ可能性がある。近年、磁気共鳴画像法(MRI, magnetic resonance imaging)は、椎骨解析のための価値ある非侵襲的ツールになっている。コンピュータ断層撮影法(CT, computed tomography)とは異なり、MRIは放射線被曝に関連するリスクを引き起こさない。しかし、磁気共鳴画像解析は、異なるスキャンにわたる解像度及び強度範囲における変動のような、セグメンテーションのためのいくつかの課題に直面する。
上記のニーズは、他のニーズと共に、この開示に記載されるように椎間板をモデリングするための方法の実施形態により対処される。
脊椎分析(脊椎解剖学)セグメンテーションのための方法が開示される。当該方法は、完全畳み込みネットワーク(fully convolutional network)と残差ニューラルネットワーク(residual neural network)とを結合することを含む。当該方法はまた、残差ニューラルネットワークと共に結合された完全畳み込みネットワークを端から端まで訓練することを含む。当該方法はまた、脊椎分析の少なくとも1つの医用画像を受信することを含む。当該方法はまた、残差ニューラルネットワークと共に完全畳み込みネットワークを少なくとも1つの医用画像に適用し、脊椎分析の少なくとも1つの医用画像から少なくとも1つの椎体(vertebral body)をセグメンテーションすることを含む。
脊椎分析セグメンテーションのための他の方法が開示される。当該方法は、完全畳み込みネットワークと残差ニューラルネットワークとを結合することを含む。当該方法はまた、残差ニューラルネットワークと共に結合された完全畳み込みネットワークを端から端まで訓練することを含む。当該方法はまた、脊椎分析の少なくとも1つの磁気共鳴画像を受信することを含む。当該方法はまた、残差ニューラルネットワークと共に完全畳み込みネットワークを受信した少なくとも1つの磁気共鳴画像に適用し、脊椎分析の少なくとも1つの磁気共鳴画像から少なくとも1つの椎体をセグメンテーションすることを含む。当該方法はまた、セグメンテーションされた少なくとも1つの磁気共鳴画像から、(i)上部椎骨に対応する第1の三次元椎骨面及び(ii)下部椎骨に対応する第2の三次元椎骨面を抽出することを含み、下部椎骨は上部椎骨に隣接する。当該方法はまた、第1の三次元椎骨面及び第2の三次元椎骨面を含む面再構成に基づいて、椎間板メッシュをモデリングすることを含む。
本発明の原理の理解を促進する目的で、図面に例示されており以下に記載の明細書に記載されている実施形態を参照する。これによって本発明の範囲に対する限定が意図されるものではないことが更に理解される。本発明は、例示の実施形態に対するいずれかの変更及び修正を含み、本発明が関係する技術分野の当業者に通常思い浮かぶ本発明の原理の更なる適用を含むことが更に理解される。
ここに開示される方法は、磁気共鳴画像の訓練セットから椎体及び椎弓根(pedicle)の外観を学習するタンデム畳み込みニューラルネットワークに基づく自動椎体セグメンテーション方法の改良を提供する。一例では、セグメンテーションが実行されると、椎骨レベルが自動的に識別される。他の例では、セグメンテーションが実行されると、椎骨レベルは、接続成分に基づくアルゴリズムにより画像から抽出される。他の例では、椎間板メッシュは、モデリングされる椎間板に隣接する椎体の終板(endplate)に基づいてモデリングされる。
次に、図面を参照すると、図1は、ここに記載される少なくとも1つの実施形態による椎間板モデリングのための例示的な方法のフロー図である。各図におけるブロックは順番に図示されているが、いくつかの例では、ブロックは、並列に及び/又はここに記載されているものとは異なる順番で実行されてもよい。また、様々なブロックは、所望の実装に基づいて、より少ないブロックに結合されてもよく、更なるブロックに分割されてもよく、及び/又は、除去されてもよい。
ブロック102により示すように、方法100は、完全畳み込みネットワーク(FCN)と残差ニューラルネットワーク(ResNet)とを結合することを含む。一例では、様々な画像スキャナにより取得された画像データセットは、対象物によって異なるデータ分布をもたらす可能性がある。これは、ネットワークを訓練する際に困難を生じる可能性がある。この課題を克服するために、完全畳み込みネットワークを設計する前に、バイアスフィールド(bias field)補正と等方性ボクセル(isotropic voxel)への再サンプリングが画像データセットに対して実行される。例示的な実施形態によれば、画像を平滑化するために三次元異方性拡散フィルタが適用される。さらに、完全畳み込みネットワークは、全ての入力画像の強度分布を互いに類似させるように設計される。
ブロック104により示すように、方法100はまた、残差ニューラルネットワークと共に結合された完全畳み込みネットワークを端から端まで訓練することを含む。一例では、残差ニューラルネットワークは、サイズ512×512の二次元スライスに基づいて、画素毎の予測確率値を示す512×512の出力を生成する。一実施形態では、残差ニューラルネットワークは、図2に示すようなダウンサンプリング及びアップサンプリングプロセスを含む。この実施形態では、残差ニューラルネットワークアーキテクチャのダウンサンプリングプロセスは、3のカーネルサイズ、2つの単純ブロック(simple block)及び3つの基本ブロック(basic block)を有する畳み込み層を含む。これらの層の後に、層間のボトルネック(bottleneck)が続く。アップサンプリングプロセスの後に、1×1の畳み込み分類器(classifier)が続く。
一例では、ネットワークは、勾配降下最適化アルゴリズムを使用して訓練される。例えば、ネットワークは、減衰を0.001に設定したRMSprop勾配降下最適化アルゴリズムを使用して訓練される。この例では、初期学習速度は0.001に設定され、200エポックの後に0.0001まで低下する。一例では、残差ニューラルネットワークのセグメンテーション結果を更に洗練するために、三次元の完全に接続された条件付き確率場(conditional random field)が残差ニューラルネットワーク出力に適用される。それにより、設計された完全畳み込みネットワーク及び残りのニューラルネットワークは、端から端まで訓練される。
ブロック106により示すように、方法100はまた、脊椎分析の少なくとも1つの医用画像を受信することを含む。一例では、医用画像は、外科的処置の前の関連する分析の撮影された磁気共鳴画像である。
ブロック108により示すように、方法100はまた、残差ニューラルネットワークと共に完全畳み込みネットワークを受信した少なくとも1つの医用画像に適用し、脊椎分析の少なくとも1つの医用画像から少なくとも1つの椎体をセグメンテーションすることを含む。一例では、セグメンテーションされた椎体が与えられた場合、椎骨レベルは半自動的にラベル付けされる。例えば、ユーザは、セグメンテーションされた磁気共鳴画像においてL5椎骨のレベルを示す。この例では、三次元接続成分抽出アルゴリズム(three-dimensional connected component extraction algorithm)が、異なる椎骨領域をラベル付けするために適用される。一例では、所定の閾値よりも小さい隔離成分が除去される。この例では、離散マーチングキューブアルゴリズム(discrete marching cube algorithm)が各成分に対して適用され、続いて、周波数領域において適用されるウィンドウサンプリング関数を使用するメッシュ平滑化プロセスが適用される。これは、各ボクセルに対して実行される補間カーネルとして実装される。さらに、どのレベルがユーザにより最下部の椎骨として定義されたかに依存して、残りの椎骨が順番にラベル付けされる。
一実施形態では、三次元椎骨面がセグメンテーションされた椎体から抽出されると、椎間板メッシュは、下部(lower)及び上位(superior)の椎骨の終板に基づいてモデリングされる。一例として、上部(top)椎骨の下位(inferior)終板と、隣接する下部椎骨の上位終板とが与えられた場合、メッシュの各点における法線が計算され、これは、多角形メッシュの法線を生成する。一例では、上部椎骨の下位終板は、所定の閾値未満の法線値を有する点を選択することによってクリッピングされる。この例では、下部椎骨の上位終板は、同じ所定の閾値よりも大きい法線値を有する点を選択することによってクリッピングされる。一例として、これは、隣接する面をクリッピングし、続いて、小さい隔離領域がメッシュのより大きい接続成分をフィルタ除去することにより除去される。1つのシナリオでは、異方性曲率フィルタリング技術を使用してガウス曲率が計算され、これは、2つのプレートから外れ値の曲率値を除去するのに役立つ。このシナリオでは、終板の不規則性を除去するために、15回の反復及び所定の弛緩係数を有するラプラシアンメッシュ平滑化が上部椎骨及び下部椎骨に適用される。さらに、結果の平滑化された終板は、ポアソン面再構成アルゴリズム(Poisson surface reconstruction algorithm)を使用して接続される。
図1のフロー図は、本実施形態の1つの可能な実装の機能及び動作を示す。これに関して、各ブロックは、プロセス内の特定の論理機能又はステップを実装するためにプロセッサにより実行可能な1つ以上の命令を含むモジュール、セグメント又はプログラムコードの一部を表してもよい。プログラムコードは、例えば、ディスク又はハードドライブを含む記憶デバイスのようないずれかの種類のコンピュータ読み取り可能媒体に記憶されてもよい。コンピュータ読み取り可能媒体は、レジスタメモリ、プロセッサキャッシュ又はランダムアクセスメモリ(RAM, Random Access Memory)のような短期間用のデータを記憶する非一時的なコンピュータ読み取り可能媒体、及び/又は、例えば、読み取り専用メモリ(ROM, read only memory)、光若しくは磁気ディスク又はコンパクトディスク読み取り専用メモリ(CD-ROM, compact-disc read only memory)のような永続的な長期ストレージを含んでもよい。コンピュータ読み取り可能媒体は、いずれかの他の揮発性又は不揮発性記憶システムを含むことができてもよく或いは含んでもよい。コンピュータ読み取り可能媒体は、例えば、コンピュータ読取可能記憶媒体、有形の記憶デバイス又は他の製造物と考えられてもよい。
或いは、図1における各ブロックは、プロセス内の特定の論理機能を実行するように配線された回路を表してもよい。図1に示すような例示的な方法は、全体として或いは部分的に、クラウド内の1つ又は複数のコンポーネントにより実行されてもよい。しかし、例示的な方法は、本発明の範囲から逸脱することなく、代わりに他のエンティティ又はエンティティの組み合わせ(すなわち、他のコンピュータデバイス及び/又はコンピュータデバイスの組み合わせ)により実行されてもよいことが理解されるべきである。例えば、図1の方法の機能は、コンピュータデバイス(又は1つ以上のプロセッサのようなコンピュータデバイスのコンポーネント)により完全に実行されてもよく、或いは、コンピュータデバイスの複数のコンポーネントにわたって、複数のコンピュータデバイスにわたって、及び/又はサーバにわたって、分散されてもよい。
図3は、本開示の実施形態による板面のモデリングフレームワークの例示的な説明図である。図3に示すように、画像Aは、例示的な上部椎骨モデルを示す。画像Bは、どのようにファセット(facet)の法線が計算されるかの例を示す。画像Cは、クリッピングされた終板及び計算されたガウス曲率の例示的な図である。画像Dは、最終的にクリッピングされた終板を示す。画像Eは、上部のクリッピングされた終板と下部のクリッピングされた終板との間の板面の再構成を示す。
図4は、ここに記載される少なくとも1つの実施形態による椎間板モデリングのための例示的な方法のフロー図である。各図におけるブロックは順番に図示されているが、いくつかの例では、ブロックは、並列に及び/又はここに記載されているものとは異なる順番で実行されてもよい。また、様々なブロックは、所望の実装に基づいて、より少ないブロックに結合されてもよく、更なるブロックに分割されてもよく、及び/又は、除去されてもよい。
ブロック402により示すように、方法400は、完全畳み込みネットワーク(FCN)と残差ニューラルネットワーク(ResNet)とを結合することを含む。
ブロック404により示すように、方法400はまた、残差ニューラルネットワークと共に結合された完全畳み込みネットワークを端から端まで訓練することを含む。
ブロック406により示すように、方法400はまた、脊椎分析の少なくとも1つの磁気共鳴画像を受信することを含む。
ブロック408により示すように、方法400はまた、脊椎分析の少なくとも1つの磁気共鳴画像を受信し、残差ニューラルネットワークと共に完全畳み込みネットワークを受信した少なくとも1つの時期共鳴画像に適用し、脊椎分析の少なくとも1つの磁気共鳴画像から少なくとも1つの椎体をセグメンテーションすることを含む。
ブロック410により示すように、方法400はまた、セグメンテーションされた少なくとも1つの磁気共鳴画像から、(i)上部椎骨に対応する第1の三次元椎骨面及び(ii)下部椎骨に対応する第2の三次元椎骨面を抽出することを含み、下部椎骨は上部椎骨に隣接する。
ブロック412により示すように、方法400はまた、第1の三次元椎骨面及び第2の三次元椎骨面を含む面再構成に基づいて、椎間板メッシュをモデリングすることを含む。
[例1]
一例では、ここに記載された手法の性能は、2つの公に利用可能なデータセットに対して実行された。第1のデータセットは、下部椎骨のT2-wターボスピンエコー(turbo spin-echo)三次元磁気共鳴画像を有する23個の対象物を含む。画像は、1.5T Siemensスキャナで取得され、2×1.25×1.25mm^3のボクセルサイズに再サンプリングされた。T11からL5までの各椎体について、グラウンドトルース手動セグメンテーションが利用可能であった。異なるスキャナからの画像に対する設計された自動セグメンテーション方法の性能をテストするために、SpineWebからの下部脊椎の第2の公開データセットも使用された。第2のデータセットは、異なる寸法、ボクセル間隔及び強度範囲を有する20個の場合について、T2-w MRI及びCTスキャンを含む。双方のデータセットについて、椎体及び椎弓根の手動アノテーションが実行された。
一例では、ここに記載された手法の性能は、2つの公に利用可能なデータセットに対して実行された。第1のデータセットは、下部椎骨のT2-wターボスピンエコー(turbo spin-echo)三次元磁気共鳴画像を有する23個の対象物を含む。画像は、1.5T Siemensスキャナで取得され、2×1.25×1.25mm^3のボクセルサイズに再サンプリングされた。T11からL5までの各椎体について、グラウンドトルース手動セグメンテーションが利用可能であった。異なるスキャナからの画像に対する設計された自動セグメンテーション方法の性能をテストするために、SpineWebからの下部脊椎の第2の公開データセットも使用された。第2のデータセットは、異なる寸法、ボクセル間隔及び強度範囲を有する20個の場合について、T2-w MRI及びCTスキャンを含む。双方のデータセットについて、椎体及び椎弓根の手動アノテーションが実行された。
データセットのサイズを増加させるために、訓練データセットは、ランダムフリッピング(水平及び垂直)、シアーリング(0.41の最大範囲)、回転(10の最大範囲)、ランダムクロッピング(256×256)及びスプラインワーピングを使用して補強された。脊椎の反りは、粗い3×3の格子に対するランダム変位ベクトルを使用して生成された。変位は、20ピクセルの標準偏差のガウス分布からサンプリングされた。次いで、ピクセル当たりの変位は、バイキュービック補間を使用して計算された。
椎体セグメンテーション評価は、2DCNNと3DCNNと標準ResNetとの間で比較され、データセットを5倍の交差検証に分割した。個別の検証セットを作成するために、データ量の20%が除外された。タンデムディープニューラルネットワークを実装するために、Keras及びTheanoライブラリが使用された。性能検証のために、データセットに対して利用可能な手動セグメンテーションが、グラウンドトルースとして考えられた。さらに、提案した方法からの結果がグラウンドトルースアノテーションにどれだけ近いかについて推定が実行された。自動でセグメンテーションされた構造と対応する手動でアノテーションされた構造との間の重複量を示すDice係数も計算された。さらに、等高線平均距離(CMD, contour mean distance)及びHausdorff距離(HD, Hausdorff distance)が、それぞれグラウンドトルースと自動セグメンテーションとの間の平均及び最大距離として計算された。
最後に、MRI自動セグメンテーション結果をCTグラウンドトルースセグメンテーションと比較するために、提案した出力は、利用可能なCT/MRIを有する20人の患者に対する対応するCT画像空間にレジストレーションされた。いくつかの場合、MRIスキャンはCTに比べてより多くの椎骨レベルを含み、したがって、ResNetの出力から余分なラベルが削除された。
ここに記載された発明の特徴は、目的を達成するための好ましい実施形態に関して記載されているが、本発明の真意又は範囲から逸脱することなく、これらの教示を考慮してバリエーションが達成され得ることは、当業者により認識される。
Claims (12)
- 脊椎分析セグメンテーションのための方法であって、
完全畳み込みネットワークと残差ニューラルネットワークとを結合するステップと、
前記残差ニューラルネットワークと共に前記結合された完全畳み込みネットワークを端から端まで訓練するステップと、
脊椎分析の少なくとも1つの医用画像を受信するステップと、
前記残差ニューラルネットワークと共に前記完全畳み込みネットワークを前記受信した少なくとも1つの医用画像に適用し、前記脊椎分析の前記少なくとも1つの医用画像から少なくとも1つの椎体をセグメンテーションするステップと
を含む方法。 - 前記残差ニューラルネットワークは、ダウンサンプリング及びアップサンプリングプロセスを含む、請求項1に記載の方法。
- 前記残差ニューラルネットワークと共に前記結合された完全畳み込みネットワークは、勾配降下最適化アルゴリズムを含む、請求項1に記載の方法。
- 前記残差ニューラルネットワークの出力に対して三次元の完全に接続された条件付き確率場を適用するステップを更に含む、請求項1に記載の方法。
- ユーザから入力を受信するステップであり、前記入力は、前記セグメンテーションされた少なくとも1つの医用画像に対する椎骨レベルの識別を含む、ステップと、
前記受信した入力に基づいて、三次元接続成分抽出アルゴリズムに従って1つ以上の椎骨レベルを決定するステップと
を更に含む、請求項1に記載の方法。 - ユーザから、1つ以上の決定された椎骨レベルの選択を含む第2の入力を受信するステップと、
前記受信した第2の入力に基づいて、前記選択に対応する前記セグメンテーションされた少なくとも1つの医用画像から1つ以上の三次元椎骨面を抽出するステップと、
前記抽出された1つ以上の三次元椎骨面に基づいて、椎間板メッシュをモデリングするステップと
を更に含む、請求項5に記載の方法。 - 前記セグメンテーションされた少なくとも1つの医用画像から1つ以上の三次元椎骨面を抽出するステップと、
前記抽出された1つ以上の三次元椎骨面に基づいて、椎間板メッシュをモデリングするステップと
を更に含む、請求項1に記載の方法。 - 前記医用画像は、磁気共鳴画像である、請求項1に記載の方法。
- 前記抽出された1つ以上の三次元椎骨面は、第1の面及び第2の面を含み、前記第1の面は、上部椎骨の下位終板に対応し、前記第2の面は、下部椎骨の上位終板に対応し、前記下部椎骨は前記上部椎骨に隣接する、請求項6に記載の方法。
- 前記抽出された1つ以上の三次元椎骨面は、第1の面及び第2の面を含み、前記第1の面は、上部椎骨の下位終板に対応し、前記第2の面は、下部椎骨の上位終板に対応し、前記下部椎骨は前記上部椎骨に隣接する、請求項7に記載の方法。
- 前記第1の面及び前記第2の面に沿って1つ以上の点における法線ベクトルを決定するステップと、
前記法線ベクトルに関連する値に対応する所定の閾値に従って、前記第1の面及び前記第2の面をクリッピングするステップと、
前記クリッピングされた第1の面及び前記クリッピングされた第2の面を含む面再構成に基づいて、前記椎間板メッシュをモデリングするステップと
を更に含む、請求項6に記載の方法。 - 脊椎分析セグメンテーションのための方法であって、
完全畳み込みネットワークと残差ニューラルネットワークとを結合するステップと、
前記残差ニューラルネットワークと共に前記結合された完全畳み込みネットワークを端から端まで訓練するステップと、
脊椎分析の少なくとも1つの磁気共鳴画像を受信するステップと、
前記残差ニューラルネットワークと共に前記完全畳み込みネットワークを前記受信した少なくとも1つの磁気共鳴画像に適用し、前記脊椎分析の前記少なくとも1つの磁気共鳴画像から少なくとも1つの椎体をセグメンテーションするステップと、
前記セグメンテーションされた少なくとも1つの磁気共鳴画像から、(i)上部椎骨に対応する第1の三次元椎骨面及び(ii)下部椎骨に対応する第2の三次元椎骨面を抽出するステップであり、前記下部椎骨は前記上部椎骨に隣接する、ステップと、
前記第1の三次元椎骨面及び前記第2の三次元椎骨面を含む面再構成に基づいて、椎間板メッシュをモデリングするステップと
を含む方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201762575260P | 2017-10-20 | 2017-10-20 | |
US62/575,260 | 2017-10-20 | ||
PCT/US2018/056797 WO2019079778A1 (en) | 2017-10-20 | 2018-10-19 | INTERVERTEBRAL DISC MODELING |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2021500113A true JP2021500113A (ja) | 2021-01-07 |
Family
ID=64427192
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020521914A Pending JP2021500113A (ja) | 2017-10-20 | 2018-10-19 | 椎間板モデリング |
Country Status (5)
Country | Link |
---|---|
US (2) | US11276232B2 (ja) |
EP (2) | EP3872759A1 (ja) |
JP (1) | JP2021500113A (ja) |
AU (2) | AU2018351078B2 (ja) |
WO (1) | WO2019079778A1 (ja) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3872759A1 (en) * | 2017-10-20 | 2021-09-01 | Nuvasive, Inc. | Intervertebral disc modeling |
WO2019241637A1 (en) * | 2018-06-14 | 2019-12-19 | Thomas Jefferson University | A novel, quantitative framework for the diagnostic, prognostic, and therapeutic evaluation of spinal cord diseases |
CN108986891A (zh) * | 2018-07-24 | 2018-12-11 | 北京市商汤科技开发有限公司 | 医疗影像处理方法及装置、电子设备及存储介质 |
CN110246216B (zh) * | 2019-05-28 | 2023-02-28 | 中国科学院深圳先进技术研究院 | 脊柱模型生成方法、脊柱模型生成系统及终端 |
CN110223279B (zh) * | 2019-05-31 | 2021-10-08 | 上海商汤智能科技有限公司 | 一种图像处理方法和装置、电子设备 |
US10990876B1 (en) * | 2019-10-08 | 2021-04-27 | UiPath, Inc. | Detecting user interface elements in robotic process automation using convolutional neural networks |
CN112754458A (zh) * | 2019-11-01 | 2021-05-07 | 上海联影医疗科技股份有限公司 | 磁共振成像方法、系统及存储介质 |
US11157783B2 (en) | 2019-12-02 | 2021-10-26 | UiPath, Inc. | Training optical character detection and recognition models for robotic process automation |
US11980415B2 (en) | 2020-12-11 | 2024-05-14 | Nuvasive, Inc. | Robotic surgery |
CN118058835A (zh) * | 2024-03-22 | 2024-05-24 | 北京长木谷医疗科技股份有限公司 | 一种脊柱手术导航定位系统设计方法及相关设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009207727A (ja) * | 2008-03-05 | 2009-09-17 | Fujifilm Corp | 椎体位置決定装置および方法並びにプログラム |
JP2014158628A (ja) * | 2013-02-20 | 2014-09-04 | Univ Of Tokushima | 画像処理装置、画像処理方法、制御プログラム、および記録媒体 |
JP2016119954A (ja) * | 2014-12-24 | 2016-07-07 | 好民 村山 | 放射線撮影装置 |
US20160300026A1 (en) * | 2015-04-13 | 2016-10-13 | Siemens Medical Solutions Usa, Inc. | Patient Management Based On Anatomic Measurements |
WO2017106645A1 (en) * | 2015-12-18 | 2017-06-22 | The Regents Of The University Of California | Interpretation and quantification of emergency features on head computed tomography |
JP2017158781A (ja) * | 2016-03-09 | 2017-09-14 | 富士フイルム株式会社 | 画像表示制御装置および方法並びにプログラム |
CN108399616A (zh) * | 2018-02-11 | 2018-08-14 | 西安理工大学 | 基于深度残差网络的骨科病病变分类分级方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2083390B1 (de) | 2008-01-24 | 2016-06-01 | Brainlab AG | Verfahren zur Segmentierung eines 3D-Bilddatensatzes, zugehöriges Computerprogrammprodukt und zugehöriges System |
US9875540B2 (en) | 2013-07-11 | 2018-01-23 | Brainlab Ag | Generation of a patient-specific anatomical atlas |
AU2018316801B2 (en) | 2017-08-16 | 2023-12-21 | Mako Surgical Corp. | Ultrasound bone registration with learning-based segmentation and sound speed calibration |
EP3872759A1 (en) * | 2017-10-20 | 2021-09-01 | Nuvasive, Inc. | Intervertebral disc modeling |
-
2018
- 2018-10-19 EP EP21170454.9A patent/EP3872759A1/en not_active Withdrawn
- 2018-10-19 AU AU2018351078A patent/AU2018351078B2/en active Active
- 2018-10-19 EP EP18807461.1A patent/EP3698320B1/en active Active
- 2018-10-19 JP JP2020521914A patent/JP2021500113A/ja active Pending
- 2018-10-19 WO PCT/US2018/056797 patent/WO2019079778A1/en unknown
-
2020
- 2020-04-19 US US16/852,516 patent/US11276232B2/en active Active
-
2021
- 2021-11-04 US US17/519,164 patent/US11798229B2/en active Active
-
2023
- 2023-08-15 AU AU2023216778A patent/AU2023216778A1/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009207727A (ja) * | 2008-03-05 | 2009-09-17 | Fujifilm Corp | 椎体位置決定装置および方法並びにプログラム |
JP2014158628A (ja) * | 2013-02-20 | 2014-09-04 | Univ Of Tokushima | 画像処理装置、画像処理方法、制御プログラム、および記録媒体 |
JP2016119954A (ja) * | 2014-12-24 | 2016-07-07 | 好民 村山 | 放射線撮影装置 |
US20160300026A1 (en) * | 2015-04-13 | 2016-10-13 | Siemens Medical Solutions Usa, Inc. | Patient Management Based On Anatomic Measurements |
WO2017106645A1 (en) * | 2015-12-18 | 2017-06-22 | The Regents Of The University Of California | Interpretation and quantification of emergency features on head computed tomography |
JP2017158781A (ja) * | 2016-03-09 | 2017-09-14 | 富士フイルム株式会社 | 画像表示制御装置および方法並びにプログラム |
CN108399616A (zh) * | 2018-02-11 | 2018-08-14 | 西安理工大学 | 基于深度残差网络的骨科病病变分类分级方法 |
Non-Patent Citations (3)
Title |
---|
AMIR J. ET AL.: "ISSLS PRIZE IN BIOENGINEERING SCIENCE 2017: Automation of reading of radiological features from magn", EUROPEAN SPINE JOURNAL, vol. 26, JPN7022004619, 2017, pages 1374 - 1383, XP036222640, ISSN: 0004888191, DOI: 10.1007/s00586-017-4956-3 * |
FAEZEH F. ET AL.: "A novel objective function based on a generalized kelly criterion for deep learning", 2017 SIGNAL PROCESSING: ALGORITHMS, ARCHITECTURES, ARRANGEMENTS, AND APPLICATIONS (SPA), JPN7022004621, 2017, pages 84 - 89, XP033276612, ISSN: 0004888190, DOI: 10.23919/SPA.2017.8166843 * |
ROBERT K. ET AL.: "Determination of the intervertebral disc space from CT images of the lumbar spine", MEDICAL IMAGING 2014: IMAGE PROCESSING; 90342Q, vol. 9034, JPN7022004620, 2014, pages 1 - 6, ISSN: 0004888192 * |
Also Published As
Publication number | Publication date |
---|---|
EP3698320A1 (en) | 2020-08-26 |
US20200320786A1 (en) | 2020-10-08 |
EP3698320B1 (en) | 2021-07-21 |
AU2023216778A1 (en) | 2023-09-07 |
US11276232B2 (en) | 2022-03-15 |
AU2018351078A1 (en) | 2020-04-30 |
US20220058869A1 (en) | 2022-02-24 |
WO2019079778A1 (en) | 2019-04-25 |
AU2018351078B2 (en) | 2023-07-06 |
EP3872759A1 (en) | 2021-09-01 |
US11798229B2 (en) | 2023-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11798229B2 (en) | Intervertebral disc modeling | |
Horng et al. | Cobb angle measurement of spine from x‐ray images using convolutional neural network | |
Sekuboyina et al. | Attention-driven deep learning for pathological spine segmentation | |
US9898825B2 (en) | Segmentation of magnetic resonance imaging data | |
EP3509013A1 (en) | Identification of a predefined object in a set of images from a medical image scanner during a surgical procedure | |
US20120281904A1 (en) | System and method for automatic recognition and labeling of anatomical structures and vessels in medical imaging scans | |
JP2014505491A (ja) | 医療用画像の非線形分解能の低減 | |
CN110097517B (zh) | 去除图像伪影的方法及装置 | |
FR2862791A1 (fr) | Procede de segmentation d'une structure et systeme d'imagerie associe | |
CN115359257B (zh) | 基于深度学习的脊椎图像分割方法及手术导航定位系统 | |
Matzkin et al. | Cranial implant design via virtual craniectomy with shape priors | |
CA2778599C (en) | Bone imagery segmentation method and apparatus | |
Khandelwal et al. | Spine and individual vertebrae segmentation in computed tomography images using geometric flows and shape priors | |
JP4991748B2 (ja) | 構造のモデルを作成する方法 | |
Koh et al. | Automatic spinal canal detection in lumbar MR images in the sagittal view using dynamic programming | |
Susan et al. | Deep learning inpainting model on digital and medical images-a review. | |
Shakeri et al. | Inter-vertebral disk modelling from pairs of segmented vertebral models using trainable pre-processing networks | |
Liao et al. | Using three-dimensional multigrid-based snake and multiresolution image registration for reconstruction of cranial defect | |
Grossauer | Inpainting of movies using optical flow | |
CN113658119A (zh) | 一种基于vae的人脑损伤检测方法及装置 | |
CN114402395A (zh) | 基于vrds 4d医学影像的脊椎疾病识别方法及相关装置 | |
CN113689454A (zh) | 基于卷积神经网络的3d ct椎体分割算法 | |
Koh et al. | Automated boundary extraction of the spinal canal in mri based on dynamic programming | |
WO2014032041A1 (en) | Method and system for image registration | |
Khandelwal | Spine segmentation in computed tomography images using geometric flows and shape priors |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20211008 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20220928 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20221004 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20230425 |