特許 7162793 超音波拓本技術に基づく脊椎画像生成システム及び脊柱手術用のナビゲーション・位置確認システム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2022-10-21

(45)【発行日】2022-10-31

(54)【発明の名称】超音波拓本技術に基づく脊椎画像生成システム及び脊柱手術用のナビゲーション・位置確認システム

(51)【国際特許分類】

   A61B 8/14 20060101AFI20221024BHJP

   A61B 5/055 20060101ALI20221024BHJP

   A61B 6/00 20060101ALI20221024BHJP

   A61B 6/03 20060101ALI20221024BHJP

   A61B 34/20 20160101ALI20221024BHJP

【ＦＩ】

A61B8/14

A61B5/055 390

A61B6/00 370

A61B6/03 377

A61B34/20

【請求項の数】 17

(21)【出願番号】P 2020504215

(86)(22)【出願日】2018-07-25

(65)【公表番号】

(43)【公表日】2020-09-03

(86)【国際出願番号】 CN2018096999

(87)【国際公開番号】W WO2019020048

(87)【国際公開日】2019-01-31

【審査請求日】2020-03-11

【審判番号】

【審判請求日】2021-03-10

(31)【優先権主張番号】201710630811.3

(32)【優先日】2017-07-28

(33)【優先権主張国・地域又は機関】CN

【早期審査対象出願】

(73)【特許権者】

【識別番号】505072650

【氏名又は名称】浙江大学

【氏名又は名称原語表記】ＺＨＥＪＩＡＮＧ ＵＮＩＶＥＲＳＩＴＹ

(74)【代理人】

【識別番号】100146374

【弁理士】

【氏名又は名称】有馬 百子

(72)【発明者】

【氏名】劉 天健

(72)【発明者】

【氏名】朱 永堅

(72)【発明者】

【氏名】陳 高

【合議体】

【審判長】三崎 仁

【審判官】福島 浩司

【審判官】樋口 宗彦

(56)【参考文献】

【文献】特表２０１６－５４０６０４（ＪＰ，Ａ）

【文献】特開２０００－１０２５３８（ＪＰ，Ａ）

【文献】国際公開第２０１６／１７８５７９（ＷＯ，Ａ１）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

IPC A61B 8/00 - 8/15

(57)【特許請求の範囲】

【請求項1】

超音波拓本技術に基づく脊椎画像生成システムであって、二次元的な脊椎表面構造の超音波画像に基づいて超音波拓本を生成し、予め取得しておいた三次元のデジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーが得られ、
収集ユニットと処理ユニットとを含み、
前記収集ユニットは、位置確認タグを備える超音波プローブによって、患者の手術状況でのリアルタイムな手術中体位に対応する前記二次元的な脊椎表面構造の超音波画像、を取得するものであり、
前記処理ユニットは、前記収集ユニットで取得された前記患者の手術状況でのリアルタイムな手術中体位に対応する前記二次元的な脊椎表面構造の超音波画像、を重ね合わせた重ね合わせ超音波画像を得た後、前記重ね合わせ超音波画像について三次元再構成処理を行って三次元の超音波拓本を形成し、前記超音波拓本に含まれる脊椎表面のリアルタイムな動的空間位置情報と、前記デジタル医用画像に含まれる脊椎表面の空間位置情報とを三次元レベルで輪郭マッチングを行うことで、前記患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィー、を得るものであり、
前記超音波画像は、脊椎領域の深部筋組織と脊椎表面とを含むすべてのエコー情報から構成され、前記脊椎表面は横突起と棘突起とを含み、
前記超音波拓本は、前記位置確認タグを備える超音波プローブの与える患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像であり、
前記脊椎表面トポグラフィーは、前記超音波拓本が有する前記リアルタイムな動的空間位置情報を前記デジタル医用画像に付与してなる立体表面画像であり、前記立体表面画像は、脊椎表面の空間情報に基づいて形成される見掛け上の三次元画像であり、
前記超音波拓本が有する前記リアルタイムな動的空間位置情報は、脊椎表面と筋骨格境界とともに構成される空間情報であり、この空間情報は、超音波走査による脊椎表面、深部筋及び脂肪組織の超音波画像情報と空間位置情報を含み、リアルタイムで動的に表示され得るものであることを特徴とする超音波拓本技術に基づく脊椎画像生成システム。

【請求項2】

前記超音波画像は、その輪郭エッジと輪郭内部の情報を含むことを特徴とする、請求項１に記載の超音波拓本技術に基づく脊椎画像生成システム。

【請求項3】

前記脊椎表面は、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分的のいずれか又は任意の複数種をさらに含むことを特徴とする、請求項１に記載の超音波拓本技術に基づく脊椎画像生成システム。

【請求項4】

前記超音波画像の取得方法としては、患者の空間位置確認情報を有する超音波走査機器により患者体表を繰り返し走査し、筋骨格境界が識別できたら、手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を抽出し、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とし、前記手術中の二次元空間位置パラメータは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータであることを特徴とする、請求項１～３のいずれか１つに記載の超音波拓本技術に基づく脊椎画像生成システム。

【請求項5】

前記輪郭マッチングの実施方法としては、前記超音波拓本に含まれる脊椎表面のすべての骨エコー情報をデジタル医用画像の輪郭と一点ずつマッチングさせることで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーが得られることを特徴とする、請求項１～請求項４のいずれかに記載の超音波拓本技術に基づく脊椎画像生成システム。

【請求項6】

前記収集ユニットは超音波画像走査モジュールと超音波画像情報抽出モジュールとを含み、
前記超音波画像走査モジュールは、患者体表を繰り返し走査し、筋骨格境界が識別できたら、生データを取得するものであり、前記超音波画像情報抽出モジュールは、前記生データについて抽出を行い、手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を取得し、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とするものであり、前記手術中の二次元空間位置パラメータは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータであることを特徴とする、請求項１～５のいずれか１つに記載の超音波拓本技術に基づく脊椎画像生成システム。

【請求項7】

前記処理ユニットは、超音波画像最適化モジュール、超音波画像重ね合わせモジュール、超音波拓本生成モジュール及び画像輪郭マッチングモジュールを含み、
前記超音波画像最適化モジュールは、前記収集ユニットで得られた患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を最適化させることで、干渉ノイズが除去され最適化した超音波画像を得るものであり、
前記超音波画像重ね合わせモジュールは、前記超音波画像最適化モジュールで得られた最適化超音波画像を重ね合わせ、超音波骨エコー情報における強エコーと弱エコーの間の相違を一層大きくすることで、画像効果の強調した重ね合わせ超音波画像を得るものであり、
前記超音波拓本生成モジュールは、前記超音波画像重ね合わせモジュールで得られた重ね合わせ超音波画像について処理することで、超音波拓本を形成するものであり、
前記画像輪郭マッチングモジュールは、前記超音波拓本生成モジュールで得られた超音波拓本をデジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを形成するものであることを特徴とする、請求項１～６のいずれか１つに記載の超音波拓本技術に基づく脊椎画像生成システム。

【請求項8】

前記超音波画像最適化モジュールは超音波フィルタ強調技術により前記収集ユニットで取得された超音波画像を最適化させることで、干渉ノイズが除去され最適化した超音波画像を得ることを特徴とする、請求項７に記載の超音波拓本技術に基づく脊椎画像生成システム。

【請求項9】

超音波拓本技術に基づく脊椎画像生成方法であって、超音波画像取得、超音波拓本生成及び脊椎表面トポグラフィー生成を含み、
１）超音波画像取得では、位置確認タグを備える超音波プローブによって、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を取得し、前記超音波画像は脊椎領域の深部筋組織と脊椎表面とを含むすべてのエコー情報から構成され、前記脊椎表面は横突起と棘突起とを含み、
２）超音波拓本生成では、ステップ１）で取得された前記患者の手術状況でのリアルタイムな手術中体位に対応する前記二次元的な脊椎表面構造の超音波画像、を重ね合わせ重ね合わせ超音波画像を得た後、前記重ね合わせ超音波画像について三次元再構成処理を行って）三次元の超音波拓本を生成し、前記超音波拓本は、前記位置確認タグを備える超音波プローブの与える患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像であり、
３）脊椎表面トポグラフィー生成では、ステップ２）で得られた前記超音波拓本に含まれる脊椎表面のリアルタイムな動的空間位置情報と、予め取得しておいた三次元のデジタル医用画像に含まれる脊椎表面の空間位置情報とを三次元レベルで輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得、前記脊椎表面トポグラフィーは、前記超音波拓本が有する前記リアルタイムな動的空間位置情報をデジタル医用画像に付与してなる立体表面画像であり、前記立体表面画像は、脊椎表面の空間情報に基づいて形成される見掛け上の三次元画像であり、前記超音波拓本が有する前記リアルタイムな動的空間位置情報は、脊椎表面と筋骨格境界とともに構成される空間情報であり、この空間情報は、超音波走査による脊椎表面、深部筋及び脂肪組織の超音波画像情報と空間位置情報を含み、リアルタイムで動的に表示され得るものであることを特徴とする超音波拓本技術に基づく脊椎画像生成方法。

【請求項10】

前記超音波画像は、その輪郭エッジと輪郭内部の情報を含むことを特徴とする、請求項９に記載の超音波拓本技術に基づく脊椎画像生成方法。

【請求項11】

前記脊椎表面は、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分的のいずれか又は任意の複数種をさらに含むことを特徴とする、請求項９に記載の超音波拓本技術に基づく脊椎画像生成方法。

【請求項12】

前記超音波画像の取得方法としては、患者の空間位置確認情報を有する超音波走査機器により患者体表を繰り返し走査し、筋骨格境界が識別できたら、手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を抽出し、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とし、前記手術中の二次元空間位置パラメータは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータであることを特徴とする、請求項９～１１のいずれか１つに記載の超音波拓本技術に基づく脊椎画像生成方法。

【請求項13】

前記脊椎表面トポグラフィーの生成方法としては、前記超音波拓本に含まれる脊椎表面のすべての骨エコー情報をデジタル医用画像の輪郭と一点ずつマッチングさせることで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーが得られることを特徴とする、請求項９～１２のいずれか１つに記載の超音波拓本技術に基づく脊椎画像生成方法。

【請求項14】

超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムであって、ナビゲーションモジュールと請求項１～８のいずれか１つに記載の超音波拓本技術に基づく脊椎画像生成システムとを含み、前記患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィー、を得て、この脊椎表面トポグラフィーに基づいて手術中ナビゲーションをリアルタイムで行うことができ、
前記超音波拓本技術に基づく脊椎画像生成システムでは、前記患者の手術状況でのリアルタイムな手術中体位に対応する前記二次元的な脊椎表面構造の超音波画像、について処理して超音波拓本を形成し、前記デジタル医用画像と前記輪郭マッチングを行うことで、前記患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィー、を形成し、前記超音波拓本は、前記患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像、であり、
前記ナビゲーションモジュールは、前記患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィー、に含まれる空間位置情報に基づいてトポグラフィー座標系を確立し、位置確認機能を有する手術器械のソリッド幾何モデルを表面モデルへ変換させ、前記脊椎画像生成システムで得られた前記患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィー、に示し、また、患者の手術中体位の前記空間位置確認情報から患者座標系を確立し、前記トポグラフィー座標系と前記患者座標系における対応点間の関係を特定してから統一座標系にし、また、変換された前記手術器械の前記表面モデルを前記統一座標系に組み込んで、リアルタイムな手術中ナビゲーションが実現され、手術操作を行うように手術者を指導することを特徴とする超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システム。

【請求項15】

前記デジタル医用画像は、ＣＴボリュームレンダリング再構成画像や磁気共鳴ＭＲ、コンピュータＸ線撮影ＣＲ、デジタルＸ線撮影ＤＲによる画像から選ばれることを特徴とする、請求項１４に記載の超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システム。

【請求項16】

前記位置確認機能を有する手術器械は、磁気的な空間位置確認タグを有する手術器械であることを特徴とする、請求項１４に記載の超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システム。

【請求項17】

脊柱手術操作システムであって、請求項１～８のいずれか１つに記載の超音波拓本技術に基づく脊椎画像生成システム、及び／又は、請求項１４～１６のいずれか１つに記載の超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムを含むことを特徴とする脊柱手術操作システム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は医学の分野に関し、具体的には、脊椎画像生成システム、脊椎画像生成方法、並びに、脊柱手術用のナビゲーション・位置確認システム、脊柱手術用のナビゲーション・位置確認方法に関する。

【背景技術】

【0002】

現在、脊柱変性関連疾患の発症率や患者数は日増しに増えており、人々の健康や労働能力にひどく影響を与えている。脊柱手術は椎間板突出や腫瘍等の脊柱疾患を治療する主要手段であり、その中では、脊柱内視鏡に代表される低侵襲手術手技は、科学技術の進歩による手術器械の進化でますます確立されてきているだけでなく、創傷が小さく、患者の回復速度が高いという優位性から、より多くの分野において認められるようになっている。しかしながら、これらの低侵襲技術は、病変部位の位置確認精度に対する要求が高く、また、人体脊柱の特殊性、つまり骨性構造で構成され、解剖箇所が狭くて長く、重要な血管や神経が豊富である等の複雑性から、低侵襲手技を習得する敷居が高く、主としては、手術の難易度及びリスクが大きく、手術者の育成にはかなりの時間を要し、短期間に普及することは難しい。ナビゲーション・位置確認技術によれば、低侵襲手術手技の敷居を下げることができ、手術操作の簡便性の向上や手術の安全性の保障に役立つ。

【0003】

超音波で人体の骨構造を走査することは以下の特性がある。１）超音波が骨構造を透過できないため、超音波画像には骨表面による線形の強エコーしか表示できず、得られる情報が少なく、臨床応用の範囲が制限されており、前記線形の強エコーは医用超音波学の概念で、超音波画像に対する特徴的描写を言い、線形のハイライト画像として現れている。２）超音波で骨構造を走査する場合、回折によるノイズが発生しやすく、骨構造の超音波画像の表示にばらつきが生じるようになり、前記ノイズは、医用超音波学の概念で、超音波画像に現れるぼけやばりのある不鮮明な画像を言い、骨構造のノイズは、ぼけやばりのある不鮮明な線形又は片状のハイライト画像として現れている。現在の臨床に応用されていた超音波走査技術（超音波走査機器及び超音波処理技術を含む）は上記問題を解決できず、そのため、骨構造関連疾患の診療においての応用はかなり制限されていた。

【0004】

上記した超音波で骨構造を走査することの特性は、超音波で脊椎を走査する場合にも同様に存在しており、それ以外、人体中の他の骨構造、例えば上腕骨や大腿骨等とは異なり、脊椎構造を超音波で走査する場合には以下のような特殊性がある。例えば図６（１）に示されるように、人体脊柱は狭くて長い骨性構造の領域で、体の支柱となるものであり、頚椎セグメント、胸椎セグメント、腰椎セグメント、仙椎セグメント及び尾椎セグメントの五つの部分に分けられ、３３個の椎骨（頚椎骨７個、胸椎骨１２個、腰椎骨５個、仙骨と尾骨の合計９個）が靱帯や関節及び椎間板で接続されてなる。脊柱内部には、上から下へ縦列の脊柱管が１本形成されており、脊柱管には脊髄がある。脊椎は、横突起、棘突起、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分で構成される。超音波で脊椎を走査することの特殊性は以下のとおりである。１）超音波走査による脊椎骨情報が少なく、脊柱は狭くて長い骨性構造で、複数個の椎骨で構成されるものから、一つの超音波走査画像にすべての椎骨セグメントを表すことができず、情報が少ないながら骨表面エコーに関する情報のみとなる。２）他の骨構造に比べて、超音波走査による脊椎骨情報にはより大きくて多くのばらつきがあり、例えば図６（２）、図６（３）に示されるように、脊椎は、数多くの骨性突起、例えば棘突起、横突起、関節突起等を有する一方、超音波画像において骨エコーが強く、超音波画像は骨構造に特有の明るくてノイズのたくさんある不鮮明な画像として現れることが多く、骨エコーに起因した回折とノイズにより、脊椎の解剖構造（横突起、棘突起、関節突起、椎板及び脊椎を構成する他の骨格部分を含む）の超音波画像での表示に大きなばらつきが発生する。よく発生するばらつきとしては、前記解剖構造中のある点が回折の原因で超音波画像において一つの面として現れることや、前記解剖構造中のある面と超音波画像中の面との間の大きさや割合が正確ではないことや、前記解剖構造中の鮮明な隣接・接続構造、例えば関節突起関節や横突起と椎体との接続根元等が超音波画像においてぼけた面状構造として現れ、接続箇所が不鮮明になることがある。上記のばらつきにより、一般的な超音波走査方式では脊椎の骨構造に対する超音波画像の鮮明度が低いという問題を解決できず、超音波で骨構造を走査する場合の多くの問題の他、超音波で脊椎骨を走査する場合の情報にはより多くて大きな問題がある。

【0005】

あるセグメントの脊椎、例えば図６（２）に示される腰椎については、一個の腰椎が有する解剖構造は、横突起、棘突起、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分を含み、異なる解剖構造は超音波画像での表示が著しく相違している。つまり、１）棘突起は、椎体における人体の背側へ伸びる骨構造であり、その構造が狭くて長く、骨が薄く、棘突起先端という識別されやすい特徴があり、前記骨が薄いこととは、棘突起の人体の横断面に平行な方向では、左側の骨表面から右側の骨表面までの厚さが薄いことを言い、前記棘突起先端は棘突起の解剖学的構造であり、棘突起の脊柱方向から離れた根元は尖状構造をなし、棘突起の末端となることで、棘突起が脊椎における最も特徴的な解剖構造として超音波画像において完全に表示されることは非常に容易になる。２）横突起、関節突起などの構造は、超音波画像での表示が棘突起よりも悪いものであり、これは、その両者は棘突起のように背側へ伸びて椎体と大きな夾角が出来るものではなく、椎体との夾角が小さく、棘突起先端のような識別されやすい特徴がない一方、体位や走査角度に制限されることで、両者と隣り合う骨構造例えば椎体から鮮明に見分けることができない。また、上記の超音波で骨構造を走査する場合のばらつきの原因で、超音波走査による両者の骨エコー情報は完全ではなく、十分に鮮明かつ正確なものにならないためであり、前記した体位や走査角度に制限されることとは、横突起や関節突起の骨エコー情報を取得する場合、脊椎の他の構造や脊椎に隣り合う器官を走査してしまうことは不可避的であるから、脊椎の他の構造や脊椎に隣り合う器官の画像により横突起や関節突起の骨エコー情報に干渉を与えて、横突起や関節突起の骨エコー情報は完全ではなく、十分に鮮明かつ正確なものにならないことを言う。３）椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分については、上記した超音波で骨構造を走査する場合の数多くの欠陥及びこれらの構造間の特殊な解剖学的位置の関係の原因で、これらの構造まで超音波が走査できず、あるいは、超音波画像においてこれらの構造のすべての骨エコー情報を表示することは不可能である。これらの構造の一部の情報を走査できたとしても、情報の正確度が不十分でばらつきが発生している問題がある。具体的には、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分の解剖構造では、ある点が回折の原因で超音波画像において一つの面として現れることや、前記解剖構造中のある面と超音波画像中の対応する面との間の大きさや割合が正確ではないことや、前記解剖構造中の接続部位例えば横突起と椎体との接続根元（図６（２）の１５を参照）が超音波画像においてぼけた面状構造として現れ、接続箇所が不鮮明になり、正確に識別できなくなることがある。

【0006】

現在、臨床での脊柱低侵襲手術は以下のような手順となる。つまり、１）手術前に患者に対してＣＴ走査を行い、２）手術のための位置確認を行い、即ち、現在の手術状況での患者の手術中体位を例えばＸ線透視技術により把握し、３）手術のための位置確認が完了した後、低侵襲手術や開放手術を含む手術操作を行う。そのうち、一回当たりのＸ線透視では一つの静的画像しか生成できず、手術のための位置確認によっても、ある時刻での患者の手術中体位しか把握できないが、正常人の脊柱は一定の活動度があり、各部位の活動度も異なり、頚椎セグメントや腰椎セグメントは活動度が大きく、胸椎セグメントは活動度が極めて小さく、仙椎セグメントはほとんど活動度がなく、そして、脊柱の天然の活動度及び正常の呼吸運動により、脊柱の位置や湾曲度は手術操作中に変位するようになり、一つの静的画像だけでは、リアルタイムな患者体位を正確に反映することはできない。医者による手術には、患者の手術状況でのリアルタイムな手術中体位を把握することが要求されており、現在では、医者は経験的に手術を行わなければならず、医者の経験及び手術時の状態等の主観的要因に対する手術の依存度が高く、信頼性が高くない。

【0007】

以上から分かるように、超音波で骨構造を走査することの特性についてまとめると、１）超音波で骨構造を走査する場合に問題があることや、２）他の骨構造に比べて、脊椎を超音波で走査する難易度が非常に高いことや、３）臨床向けの従来の超音波走査技術（超音波走査機器と超音波処理技術を含む）では走査情報が完全ではなく、鮮明度や正確度が高くはないという欠陥があり、脊柱手術の臨床応用に適応できず、画像問題が脊柱低侵襲手術に基づく治療の発展の最大の障害となることがある。

【0008】

現在、臨床的な脊柱手術のための位置確認の最も代表的なものは、Ｃアーム型透視・位置確認技術である。この技術の操作手順は図１に示されるようにまとめることができ、即ち、「穿刺－照射－評価－調整－照射による検証」である。具体的なフローとしては、手術前に患者の手術部位に穿刺針を刺入して、一回目のＸ線照射を行った後、穿刺針の患者体内での刺入深さや角度が分かるようになり、そして調整してから二回目のＸ線照射を行い、再度穿刺針の患者体内での刺入深さや角度が分かるようになり、刺入深さや角度が予定部位になると、位置確認が完了し、手術操作を行う。そうでなければ、穿刺位置確認が完了するまで、「評価－調整－照射による検証」の工程を繰り返す。Ｃアーム型透視・位置確認技術により、正確な位置確認が確保されているが、以下のような問題点がある。

【0009】

１）手術者は経験的に判断して、患者の体表に何度も繰り返して穿刺する必要があるから、患者に大きな創傷を与えてしまう。位置確認中では、患者の手術部位に穿刺針を刺入する必要があり、穿刺針が脊椎表面まで達しなければならず、望ましい穿刺深さや角度にするには、一般に何度も調整を行いながら患者体内で繰り返して穿刺する必要があり、手術感覚が悪い。

【0010】

２）手術者は、患者の手術部位に穿刺針を刺入する度にＸ線照射を行い、穿刺針の患者体内での刺入深さや角度を把握してから調整を行い、その後再度Ｘ線照射を行い、穿刺針の患者体内での刺入深さや角度が再度分かるようになり、そして、刺入深さや角度が予定部位になると、位置確認が完了する。位置確認中にはＸ線透視が必要なものであり、それによる損傷が非常に大きく、患者と手術者に対して潜在的な放射線被害を与えてしまう。

【0011】

３）この位置確認技術は手術者の経験に対する依存度が非常に大きなものである。低侵襲手技では、病変部位の位置確認精度が高く要求されており、手術者の経験は操作中の穿刺や照射の回数と正相関の関係にある。現在の臨床業務では、一般に、十年程度の手術経験を持つ手術者でも手術中において穿刺針を６－１０回刺入する必要であり、経験の浅い手術者であれば、より回数の多い穿刺と位置確認を繰り返して行う必要になり、効率が低いだけでなく、患者と手術者に対するより多くのＸ線放射や、患者に対するより大きな穿刺損傷を与えてしまい、この技術の普及と応用に不利になっていた。

【0012】

４）Ｃアーム型透視・位置確認技術による静的画像によれば、後続の手術に対するナビゲーションを行うことができない。Ｃアーム型透視・位置確認技術による画像は、患者の現在の手術状況での手術中体位で、静的Ｘ線画像であり、手術中の脊柱の動きによる変化を正確に反映できず、後続の手術操作をリアルタイムで指導できない。そのため、Ｃアーム型透視・位置確認技術は、手術前の位置確認を行い、手術のアプローチポイントを決定するためだけに適用されており、決定されたものに基づく手術が開始されると、患者の手術中のどの時点の手術状況における手術中体位も提示することができなくなり、手術者は完全に経験に基づき手術操作を行うことになる。

【0013】

Ｃアーム型透視・位置確認技術の優位性や問題点に鑑みて、近年、複数種類の画像から融合画像を形成し（劉彦斌等、「超音波ボリュームナビゲーション技術により経皮的後外側アプローチでの完全内視鏡下の腰椎低侵襲手術のための椎間孔を介する穿刺を誘導する応用と検討」、修士学位論文、第二軍医大学、２０１５年５月）、さらに位置確認を指導する方法が提案されており、その代表的なものは超音波ボリュームナビゲーション技術である。この技術の操作手順は図２に示されるようにまとめることができる。つまり、超音波走査により取得された超音波画像を手術前に収集されたＣＴ又はＭＲＩ三次元データと画像融合を行う。画像融合が完了した後、電磁追従システム（電磁トラッカーから超音波プローブまでの距離＜８０ｃｍ、手術を通してトラッカーの位置が変わらない）を用いて、超音波専門家の協力で脊柱手術者はリアルタイム画像から病変部位を見つける。具体的なフローとしては、手術前にＣＴ走査を行って三次元再構成を行い、三次元ＣＴデータを超音波計測器に導入し、三点法で画像レジストレーションを行う。同一画面にＣＴ画像と超音波画像を表示し、患者体表における４つの標識点から任意の３点（３点は同一直線上にない）を選択し、３点に関連するＣＴ画像及び超音波画像のそれぞれを順に精確に探してロックすることで、外標記点のレジストレーションが完了し、外標記点は、皮膚を切開又は穿刺して金属ピン等の器材により設置される必要のあるものである。外標記点のレジストレーションが完了した後、内標記点のレジストレーションを行い、ＣＴ図から任意の椎体の棘突起の頂点を探してロックし、そして超音波画像から対応点をリアルタイムで探してロックし、レジストレーション誤差を計測する。誤差＞５ｍｍであれば、他の椎体の棘突起の頂点を別の内標記点として選択し補正する。内外標記点の二重のレジストレーションにより、画像レジストレーションの誤差は５ｍｍ以内に縮小され、画像融合が完了した。画像融合の後、電磁トラッカーと患者との相対位置が変わらないようにし、消毒済みの隔離シースでプローブを包んで腹腔鏡を用いてリアルタイムで穿刺を誘導する。ＣＴ画像における病変セグメントの椎間孔の最下縁を刺入点としてロックし、超音波画像によりこの点をリアルタイムで表示しながら精確に探してロックし、超音波計測器によるナビゲーションにより穿刺針を病変セグメントの平面で穿刺させるようにリアルタイムで誘導する。穿刺前の位置確認及び穿刺中の穿刺方向の確認は、Ｃアーム型透視装置の協力により行うことができる。穿刺中、患者に下肢の放散痛が生じると、Ｃアーム型透視装置により透視を行う必要があり、穿刺方向がずれると、誘導穿刺針が確実に所定の位置に達するまで、標記点に基づいて再度画像レジストレーション及び融合を行う。穿刺が成功した後、Ｃアーム型透視装置により透視を行い、さらに穿刺針の位置を確認する。

【0014】

超音波ボリュームナビゲーション技術によれば、手術者はレジストレーション画像に基づいて穿刺と手術操作を行うことが可能になり、Ｃアーム型透視・位置確認技術に比べて、Ｘ線による放射線被害が大幅に低減されるが、以下のような問題点がある。

【0015】

１）患者体表に何度も繰り返して穿刺する必要があり、患者に対して大きな創傷を与えてしまい、操作が複雑で、手順が冗長且つ面倒なものとなり、時間がかかってしまう。ＣＴ画像と超音波画像とのレジストレーション融合では３つの標記点を設置する必要があり、また、皮膚を切開又は穿刺し、針やネジ又は金属釘等の器材により外標記点を設置する必要がある。内標記点は、直接に画像を見ながら画像から選択したものである。外標記点と内標記点のレジストレーションを行い補正することで、画像融合が完了した。融合した画像により穿刺針を誘導して標的部位まで達するように穿刺させ、穿刺方向がずれると、穿刺針が確実に所定の位置に達するまで、外標記点と内標記点に基づいて再度画像レジストレーション及び融合を行う。

【0016】

２）手術を通して超音波専門家の協力は必要なものである。超音波画像における解剖部位の識別、レジストレーション融合に必要な超音波画像及びＣＴ画像の選択、並びに、画像融合後の病変位置の確認のすべての手順において、手術者と超音波専門家が緊密に協力しなければならないから、手術者は迅速且つ独自に手術をすることは不可能である。

【0017】

３）この位置確認技術は、脊柱手術者の経験に対する依存度が大きなものである。この位置確認技術によるナビゲーション画像の精度が高くないが、低侵襲手技は病変部位の位置確認精度に対する要求が高いため、手術者は一定の操作経験を持っていなければ、正確な操作を確実に行えず、そのため、この技術は経験に対する依存度がＣアーム型透視・位置確認技術よりも小さいが、手術者に対する依存から脱却できなかった。

【0018】

４）ナビゲーション画像の精度が高くない。超音波ボリュームナビゲーション技術による融合画像について、ナビゲーション・位置確認の時、その画像の精度は主に三つのことによって決まる。第一は、点レジストレーションモードは精度が低く、超音波画像とＣＴ又はＭＲＩ画像には多大な情報が含まれており、幾つかの特徴点のみを選択してレジストレーションを行うと、精度を保証できず、手術に悪影響を与えてしまうことである。第二は、以上の幾つかの特徴点の選択は超音波専門家に依存し、主観的要因に影響されており、画像の精度に悪影響を与えてしまうことである。第三は、超音波画像の精度が低く、脊柱は狭くて長い骨性構造の領域であるが、超音波画像中の骨エコーが強く、超音波画像は骨構造に特有の明るくてノイズのたくさんある不鮮明な画像として現れることが多く、骨エコーに起因した回折とノイズにより、解剖構造の超音波画像での表示にばらつきが発生することである。よく発生しているばらつきとしては、解剖構造中のある点が超音波画像において一つの面として現れることや、解剖構造中の面と超音波画像中の面との間の大きさや割合が正確ではないことや、解剖構造中の鮮明な隣接・接続構造、例えば関節突起関節や横突起と椎体との接続根元等が超音波画像においてぼけた面状構造として現れ、接続箇所が不鮮明になることがある。上記のばらつきのため、一般的な超音波走査方式では、骨構造の超音波画像の鮮明度が低いという問題を解決できず、ＣＴ又はＭＲＩ画像の鮮明度よりもはるかに低いものとなり、超音波画像から選択した標記点とＣＴ又はＭＲＩ画像から選択した標記点は完全に正確に一対一に対応できず、ナビゲーション画像の精度が高くなくなった。

【0019】

５）超音波ボリュームナビゲーションによる位置確認技術は操作が面倒で、時間と労力がかかる。脊柱内視鏡を用いて椎板後部アプローチ髄核摘出術により腰椎椎間板ヘルニアを治療することを例とすれば、穿刺を開始してから標的部位に到達するまでは、Ｃアーム型透視装置により透視を行う必要があり、穿刺方向がずれると、誘導穿刺針が確実に所定の位置に達するまで、標記点に基づいて再度画像レジストレーション及び融合を行う必要がある。穿刺が成功した後も、Ｃアーム型透視装置により透視を行い、さらに穿刺針の位置を確認する必要がある。上記の位置確認と穿刺が完了した後、脊柱内視鏡を入れて深くまで徐々に降下させ、手術者は、臨床経験に基づき、棘突起、黄色靱帯、椎板、神経根、椎間板等の解剖構造を識別し、突出した髄核又は線維輪等に到達し摘出を行って関連手術操作が完了するまで、手術者は臨床経験に基づき手術器械を用いて操作を行うことしかできない。

【0020】

上記をまとめると、従来の脊柱ナビゲーション・位置確認技術は、主にＣアーム型透視・位置確認技術により精度の高い穿刺位置確認が実現されているが、複数回のＸ線照射が行われることで放射線被害が大きく、複数回の穿刺が行われることで穿刺創傷が大きく、位置確認技術の手術者の経験に対する依存度が大きく、手術時間が長くて時間と労力がかかるという問題点がある。そして得られた静的画像は後続の手術ナビゲーションの助けになることができず、また、超音波ボリュームナビゲーション技術によれば、Ｘ線による放射線被害は大幅に低減されるが、同様に、複数回の穿刺が行われることで穿刺創傷が大きく、超音波専門家の協力が必要となり、ナビゲーション画像の精度が低いという欠点がある。そのため、この二つの技術に共通の欠点としては、操作が複雑で、手術者の経験に対する依存度が大きいこと、患者体表に何度も繰り返して穿刺し、損傷が大きいこと、手術者と患者へのＸ線放射による放射線被害が大きいこと、位置確認操作が面倒で、時間と労力がかかること、体表での穿刺位置確認しか実現できず、後続の関連手術操作に対するナビゲーションが実現できない（ナビゲーション画像の精度が高くなく、ナビゲーションシステムによる支援がない）ことがある。

【0021】

このため、臨床においては、正確かつ迅速なもので放射がなく、創傷のないリアルタイムなナビゲーション・位置確認手段により脊柱手術を指導することは望ましく、それによれば、脊柱手術の精確な治療操作が実現可能になり、そのためには、少なくとも必要な部位（横突起と棘突起）に関する正確な空間位置情報を取得する必要があるが、従来技術では棘突起の正確な空間位置情報のみが得られ、横突起の正確な空間位置情報を取得できないことで、脊柱手術の精確な治療操作の正確度が高くなくなっていた。さらには、脊柱手術において特定部位の正確な空間位置情報を取得する必要もあり、前記特定部位とは、関節突起、椎板、椎板間隙、椎間孔のいずれか又は任意の複数種を言う。要するに、従来技術では、脊柱手術に必要な関連部位（必要部位と特定部位を含む）の正確な空間位置情報を完全かつ正確に取得できなかった。そのため、脊柱手術の精確な治療操作を実現するには、手術に必要とする必要部位の情報を正確でなければならず、前記必要部位の情報が正確である場合のみにおいて脊柱手術の精確な治療操作の正確性を確保でき、また、手術に必要な特定部位の空間位置情報を正確でなければならず、前記特定部位の情報が多くかつ正確であればあるほど、脊柱手術の精確な治療操作の正確性が向上することができる。従来技術では、手術に必要な特定部位の十分かつ正確な空間位置情報を取得して、ナビゲーション技術により脊柱手術の精確な治療操作を指導することはまだできていなかった。

【発明の概要】

【発明が解決しようとする課題】

【0022】

従来技術の問題点に対して、本発明は超音波拓本技術に基づく脊椎画像生成システムと方法を提供し、つまり、進歩性について述べると、取得された患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像について処理し、超音波拓本が形成され、前記超音波拓本は患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化した三次元的な超音波骨格画像であり、前記超音波拓本をデジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーがさらに得られ、それにより、従来技術による脊柱手術に必要な関連部位（必要部位と特定部位を含む）の正確な空間位置情報を取得できない問題、及び、従来技術によるリアルタイム性が遅れリアルタイムで動的に表示できない問題は徹底的に解決されるようになり、ナビゲーション技術により脊柱手術の精確な治療操作を指導することは実現されている。本発明は、超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムと方法をさらに提供し、つまり、前記した超音波拓本技術に基づく脊椎画像生成システムで得られた患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーによれば、リアルタイムな手術中ナビゲーションが実現可能になり、手術者は経験に依存することなく、超音波専門家の協力を必要とせず、患者の体に複数回穿刺することによる複数回の創傷がなくなり、Ｘ線放射による患者と手術者への損傷を最大限低減しながら、一回の穿刺により標的部位まで正確かつ迅速に到達することを可能にし、ナビゲーション技術により脊柱手術の精確な治療操作を指導することは実現されている。

【0023】

本発明による脊柱手術用のナビゲーション・位置確認システムによれば、標的部位まで到達するように穿刺した後、後続の関連手術操作を終えるように手術者を指導して、後続の手術のためのナビゲーションを継続して行うこともでき、操作が簡単で、時間も労力も節約されている。脊柱内視鏡を用いて椎板後部アプローチ髄核摘出術により腰椎椎間板ヘルニアを治療することを例とすれば、標的部位まで到達するように穿刺した後、手術ナビゲーション・位置確認システムによる指導で脊柱内視鏡を深くまで徐々に降下させ、棘突起、黄色靱帯、椎板、神経根、椎間板等の解剖構造の識別を行い、突出した髄核又は線維輪等に到達し摘出を行って関連手術操作を終え、それにより、従来技術によるＸ線と穿刺に起因した損傷が大きく、操作経験に対する依存度が高く、正確にリアルタイムでナビゲーション画像を表示し難い等の問題点を克服した。

【課題を解決するための手段】

【0024】

超音波拓本技術に基づく脊椎画像生成システムであって、二次元的な脊椎表面構造の超音波画像に基づいて超音波拓本を生成し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーが得られ、
収集ユニットと処理ユニットとを含み、
前記収集ユニットは、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を取得するものであり、前記処理ユニットは、収集ユニットで取得された患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像について処理して超音波拓本を形成し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得るものであり、前記超音波画像は脊椎領域の深部筋組織と脊椎表面とを含むすべてのエコー情報から構成され、前記脊椎表面は横突起と棘突起とを含み、前記超音波拓本は患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像である。

【0025】

本発明による超音波拓本技術に基づく脊椎画像生成システムに基づいて、超音波拓本技術に基づく脊椎画像生成方法をさらに提供する。前記の脊椎画像生成方法は、超音波画像取得、超音波拓本生成及び脊椎表面トポグラフィー生成を含み、
１）超音波画像取得では、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を取得し、前記超音波画像は脊椎領域の深部筋組織と脊椎表面とを含むすべてのエコー情報から構成され、前記脊椎表面は横突起と棘突起とを含み、
２）超音波拓本生成では、ステップ１）で取得された患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像について処理して超音波拓本を生成し、前記超音波拓本は患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像であり、
３）脊椎表面トポグラフィー生成では、ステップ２）で得られた前記超音波拓本をデジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得る。

【0026】

本発明による超音波拓本技術に基づく脊椎画像生成システムに基づいて、超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムをさらに提供する。前記脊柱手術用のナビゲーション・位置確認システムは、ナビゲーションモジュールと超音波拓本技術に基づく脊椎画像生成システムとを含み、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得て、この脊椎表面トポグラフィーに基づいてリアルタイムな手術中ナビゲーションを行うことができ、
前記超音波拓本技術に基づく脊椎画像生成システムでは、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像について処理して超音波拓本を形成し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを形成し、前記超音波拓本は、患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像である。好ましくは、前記デジタル医用画像は、ＣＴボリュームレンダリング再構成画像や磁気共鳴ＭＲ、コンピュータＸ線撮影ＣＲ、デジタルＸ線撮影ＤＲによる画像である。

【0027】

本発明による超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムに基づいて、超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認方法をさらに提供する。前記方法はステップ１）～ステップ３）を含み、
ステップ１）では、手術器械のソリッド幾何モデルを確立し表面モデルへ変換させてから、脊椎画像生成システムで得られた患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに示すことで、手術器械に対するリアルタイムな追従が実現され、前記ソリッド幾何モデルの表面モデルへの変換方法は三次元ＣＡＤモデリングが好ましく、
ステップ２）では、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに含まれる空間位置情報に基づいてトポグラフィー座標系を確立し、患者の手術中体位の空間位置情報から患者座標系を確立し、トポグラフィー座標系と患者座標系における対応点間の関係を特定してから統一座標系にし、前記したトポグラフィー座標系と患者座標系における対応点間の関係を特定する方法は反復最近接点アルゴリズムが好ましく、
ステップ３）では、ステップ１）で得られた手術器械の表面モデルをステップ２）で形成された統一座標系に表示し、手術器械と患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーとが統一座標系にあるようになり、リアルタイムな手術中ナビゲーションが実現される。

【0028】

脊柱手術操作システムをさらに提供し、前記脊柱手術操作システムは、前記いずれかの超音波拓本技術に基づく脊椎画像生成システム及び／又は前記いずれかの超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムを含む。この操作システムは脊柱部位の開放・低侵襲手術に適用され、椎間板の突出部分の摘出、腫瘍切除、神経剥離、血腫除去等のみに限られない様々な手術操作を行うことができる。

【発明の効果】

【0029】

従来技術に比べて、本発明による超音波拓本技術に基づく脊椎画像生成システムと方法によれば、従来技術による脊柱手術に必要な関連部位（必要部位と特定部位を含む）の正確な空間位置情報を取得できない問題、及び、従来技術によるリアルタイム性が遅れリアルタイムで動的に表示できない問題は解決され、ナビゲーション技術により脊柱手術の精確な治療操作を指導することは実現されている。本発明による超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムと方法によれば、リアルタイムな手術中ナビゲーションが実現可能になり、手術者は経験に依存することなく、超音波専門家の協力を必要とせず、患者の体に複数回穿刺することによる複数回の創傷がなくなり、Ｘ線放射による患者と手術者への損傷を最大限低減しながら、一回の穿刺により標的部位まで正確かつ迅速に到達することを可能にし、ナビゲーション技術により脊柱手術の精確な治療操作を指導することは実現されている。

【0030】

この脊柱手術用のナビゲーション・位置確認システムによれば、標的部位まで到達するように穿刺した後、後続の関連手術操作を終えるように手術者を指導して、後続の手術のためのナビゲーションを継続して行うことができ、操作が簡単で、時間も労力も節約されている。脊柱内視鏡を用いて椎板後部アプローチ髄核摘出術により腰椎椎間板ヘルニアを治療することを例とすれば、標的部位まで到達するように穿刺した後、手術ナビゲーション・位置確認システムによる指導で脊柱内視鏡を深くまで徐々に降下させ、棘突起、黄色靱帯、椎板、神経根、椎間板等の解剖構造の識別を行い、突出した髄核又は線維輪等に到達し摘出を行って関連手術操作を終える。具体的にまとめると、以下の効果がある。

【0031】

１）Ｘ線による患者と手術者への放射線被害を最大限低減している。

【0032】

２）患者体表に１回穿刺すればよく、穿刺による患者への創傷を最大限低減している。

【0033】

３）手術者の経験に対する依存度が非常に小さくなり、脊椎画像生成システムによれば、患者体位の変動と共にリアルタイムで更新される個人化の脊椎表面トポグラフィーが提供可能になり、及び、この脊椎表面トポグラフィーに基づいて手術中ナビゲーションをリアルタイムで行う脊柱手術用のナビゲーション・位置確認システムによれば、手術者は経験に依存することなく、超音波専門家の協力を必要とせず、脊柱手術の普及と応用に役立っている。

【0034】

４）患者体位の変動と共にリアルタイムで更新される個人化の脊椎表面トポグラフィーは、後続の手術ナビゲーションの大きな助けになる。

【0035】

５）本発明による脊椎画像生成システム及び脊柱手術用のナビゲーション・位置確認システムによれば、操作が簡単になり、時間も労力も節約されている。

【0036】

６）ナビゲーション画像は精度が高く、患者のリアルタイムな手術中体位に合致している。

【0037】

７）本発明による脊椎画像生成システム及び脊柱手術用のナビゲーション・位置確認システムによれば、位置確認技術におけるリアルタイム性が遅れ動的に表示できない問題がなくなる。Ｃアーム型透視装置はある程度のリアルタイム性があるが、遅れることが発生し、動的に表示できないという問題がある。

【0038】

８）体表での穿刺位置確認の後でも、さらに後続の関連手術操作を指導することが可能である。脊柱内視鏡を用いて椎板後部アプローチ髄核摘出術により腰椎椎間板ヘルニアを治療することを例とすれば、標的部位まで到達するように穿刺した後（位置確認・穿刺工程）、手術者は内視鏡を深くまで徐々に降下させ、臨床経験に基づき、棘突起、黄色靱帯、椎板、神経根、椎間板等の解剖構造を識別できたら、脊柱手術用のナビゲーション・位置確認システムによる指導で手術者は正確な角度や距離で手術器械を操作し、突出した髄核又は線維輪等に到達し摘出を行って関連手術操作を終えることができる。

【0039】

従来技術に比べて、本発明の優位性としては、第一は、外標記点の設置を不要とし、創傷を大幅に低減していることであり、第二は、超音波画像の精度が向上し、さらに患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーというナビゲーション画像の精度が向上していることであり、第三は、臨床経験に対する依存度が低下し、ナビゲーション画像によれば、手術の難易度が大幅に低下し、手術者の操作経験や超音波専門家に対する依存度が低下していることである。第四は、従来技術の問題点を克服し、本発明のシステムによれば、正確な位置確認と穿刺が実現され得るばかりか、脊柱手術の後続操作に対するリアルタイムなナビゲーションも実現可能になることである。

【0040】

本発明による脊椎画像生成システム及び脊柱手術用のナビゲーション・位置確認システムによれば、脊柱内視鏡手術に代表される脊柱低侵襲手技は急速に発展し短期間で普及していくようになり、脊柱関連疾患に悩む民衆の福祉が図られている。本発明は脊柱内視鏡に代表される脊柱低侵襲手技に適用されるが、これに限らず、超音波拓本を取得し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得られるものであれば、脊柱関連疾患についての低侵襲・開放手術にも同様に適用される。

【図面の簡単な説明】

【0041】

【図1】Ｃアーム型透視・位置確認技術の操作フロー模式図である。

【図2】超音波ボリュームナビゲーション技術の操作フロー模式図である。

【図3】本発明にかかる超音波拓本技術に基づく脊椎画像生成システムの構成模式図である。

【図4】本発明にかかる超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムの操作フロー模式図である。

【図5】本発明にかかる超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認方法のフロー図である。

【図6】人体脊柱及び脊椎解剖の模式図であり、（１）は脊柱の全面視での模式図であり、Ａは前面視、Ｂは後面視、Ｃは側面視であり、（２）は腰椎解剖の模式図であり、Ａは右側面視、Ｂは上面視であり、（３）は頚椎の上面視である。

【図7】体積レンダリング中のシアー・ワープ法の基本原理の模式図である。

【発明を実施するための形態】

【0042】

以下、図面と実施例に合わせて本発明についてさらに説明する。これらの実施例は本発明を説明するためのものに過ぎず、本発明の範囲の制限となるものではないことは理解されたい。また、本発明において示唆された内容のもとに、当業者は本発明について様々な変形や変更を行うことができ、これらの等価のものは同様に本願中の特許請求の範囲に含まれていることは理解されたい。

【0043】

超音波拓本技術に基づく脊椎画像生成システムを提供し、前記システムは、二次元的な脊椎表面構造の超音波画像に基づいて超音波拓本を生成し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーが得られ、
収集ユニットと処理ユニットとを含み、
前記収集ユニットは、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を取得するものであり、前記処理ユニットは、収集ユニットで取得された患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像について処理して超音波拓本を形成し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得るものであり、前記超音波画像は脊椎領域の深部筋組織と脊椎表面とを含むすべてのエコー情報から構成され、前記脊椎表面は横突起と棘突起とを含み、前記超音波拓本は、患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像である。

【0044】

好ましくは、前記超音波画像はその輪郭エッジと輪郭内部の情報を含む。好ましくは、前記脊椎表面は、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分的のいずれか又は任意の複数種をさらに含む。

【0045】

好ましくは、前記超音波画像の取得方法としては、患者の空間位置確認情報を有する超音波走査機器により患者体表を繰り返し走査し、筋骨格境界が識別できたら、手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を抽出し、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とし、前記手術中の二次元空間位置パラメータとは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータを言う。さらに好ましくは、前記患者の空間位置確認情報を有する超音波走査機器は位置確認タグを有する超音波プローブである。

【0046】

好ましくは、前記超音波拓本の取得方法としては、前記患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を最適化させ重ね合わせてから、さらに処理して超音波拓本を生成する。

【0047】

好ましくは、前記輪郭マッチングの実施方法としては、前記超音波拓本に含まれる脊椎表面のすべての骨エコー情報をデジタル医用画像の輪郭と一点ずつマッチングさせることで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得る。

【0048】

好ましくは、前記脊椎表面トポグラフィーは、超音波拓本が有するリアルタイムな動的空間情報をデジタル医用画像に付与してなる立体表面画像である。好ましくは、前記超音波拓本が有するリアルタイムな動的空間情報は、脊椎表面と筋骨格境界とともに構成される空間情報であり、この空間情報は超音波走査による脊椎表面、深部筋及び脂肪組織の超音波画像情報と空間位置情報を含み、リアルタイムで動的に表示され得るものである。好ましくは、前記立体的な表面様子は脊椎表面の空間情報に基づいて形成される見掛け上の三次元画像であり、前記見掛け上の三次元画像中のグルーヴは、前記すべての超音波画像情報に含まれる脊椎表面のテクスチャである。

【0049】

好ましくは、前記デジタル医用画像は、ＣＴボリュームレンダリング再構成画像や磁気共鳴ＭＲ、コンピュータＸ線撮影ＣＲ又はデジタルＸ線撮影ＤＲによる画像であってもよい。

【0050】

さらには、前記収集ユニットは超音波画像走査モジュールと超音波画像情報抽出モジュールとを含む。

【0051】

前記超音波画像走査モジュールは、患者体表を繰り返し走査し、筋骨格境界が識別できたら、生データを取得するものであり、前記超音波画像情報抽出モジュールは、前記生データについて抽出を行い、手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を取得し、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とするものであり、前記手術中の二次元空間位置パラメータは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータである。前記生データは、患者の空間位置確認情報を有する超音波走査機器から発する無線周波信号であり、前記走査は各種の異なる角度や方向で行われてもよく、前記方向は上下、前後、左右、斜方向を含み、前記角度は患者体表に垂直又は任意の傾斜角を含み、走査順とは無関係である。

【0052】

好ましくは、前記患者の空間位置確認情報を有する超音波走査機器は、位置確認タグを有する超音波プローブであり、前記生データは、超音波プローブから発する無線周波信号である。

【0053】

さらには、前記処理ユニットは、超音波画像最適化モジュール、超音波画像重ね合わせモジュール、超音波拓本生成モジュール及び画像輪郭マッチングモジュールを含む。

【0054】

前記超音波画像最適化モジュールは、前記収集ユニットで得られた患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を最適化させることで、干渉ノイズが除去され最適化した超音波画像を得るものである。好ましくは、前記超音波画像最適化モジュールは、超音波フィルタ強調技術により前記収集ユニットで取得された超音波画像を最適化させることで、干渉ノイズが除去され最適化した超音波画像を得るものであり、前記超音波フィルタ強調技術は適応中央値フィルタと数学的モルフォロジーフィルタとを組み合わせたものが好ましい。

【0055】

前記超音波画像重ね合わせモジュールは、前記超音波画像最適化モジュールで得られた最適化超音波画像を重ね合わせ、超音波骨エコー情報における強エコーと弱エコーの間の相違を一層大きくすることで、画像効果の強調した重ね合わせ超音波画像を得るものである。好ましくは、前記超音波画像重ね合わせに単一モダリティレジストレーション法を採用する。

【0056】

前記超音波拓本生成モジュールは、前記超音波画像重ね合わせモジュールで得られた重ね合わせ超音波画像について処理することで、超音波拓本を形成するものである。好ましくは、前記処理の方法は三次元再構成であり、さらに好ましくは、前記三次元再構成の方法は体積レンダリング法であり、体積レンダリング加速技術であることが好ましい。

【0057】

前記画像輪郭マッチングモジュールは、前記超音波拓本生成モジュールで得られた超音波拓本をデジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを形成するものである。前記輪郭マッチングの方法は、画素階調アルゴリズムのうちの相互情報アルゴリズムに基づくマルチモダリティレジストレーション法が好ましく、前記マルチモダリティレジストレーション法は、画像の空間変換、階調補間、類似度評価、検索最適化を含む。好ましくは、前記デジタル医用画像は、ＣＴボリュームレンダリング再構成画像や磁気共鳴ＭＲ、コンピュータＸ線撮影ＣＲ、デジタルＸ線撮影ＤＲによる画像から選ばれる。

【0058】

本発明による超音波拓本技術に基づく脊椎画像生成システムに基づいて、超音波拓本技術に基づく脊椎画像生成方法をさらに提供し、前記の脊椎画像生成方法は、超音波画像取得、超音波拓本生成及び脊椎表面トポグラフィー生成を含み、
１）超音波画像取得では、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を取得し、前記超音波画像は脊椎領域の深部筋組織と脊椎表面とを含むすべてのエコー情報から構成され、前記脊椎表面は横突起と棘突起とを含み、
２）超音波拓本生成では、ステップ１）で取得された患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像について処理して超音波拓本を生成し、前記超音波拓本は、患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像であり、
３）脊椎表面トポグラフィー生成では、ステップ２）で得られた前記超音波拓本をデジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得る。

【0059】

好ましくは、前記超音波画像は、その輪郭エッジと輪郭内部の情報を含む。好ましくは、前記脊椎表面は、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分的のいずれか又は任意の複数種をさらに含む。

【0060】

好ましくは、前記超音波画像の取得方法としては、患者の空間位置確認情報を有する超音波走査機器により患者体表を繰り返し走査し、筋骨格境界が識別できたら、手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を抽出し、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とし、前記手術中の二次元空間位置パラメータは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータである。さらに好ましくは、前記患者の空間位置確認情報を有する超音波走査機器は位置確認タグを有する超音波プローブである。

【0061】

好ましくは、前記超音波拓本の生成方法としては、前記患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を最適化させ重ね合わせてから、さらに処理して超音波拓本を生成する。

【0062】

好ましくは、前記脊椎表面トポグラフィーの生成方法としては、前記超音波拓本に含まれる脊椎表面のすべての骨エコー情報をデジタル医用画像の輪郭と一点ずつマッチングさせることで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得る。

【0063】

好ましくは、前記脊椎表面トポグラフィーは、超音波拓本が有するリアルタイムな動的空間情報をデジタル医用画像に付与してなる立体的な表面様子である。好ましくは、前記超音波拓本が有するリアルタイムな動的空間情報は、脊椎表面と筋骨格境界とともに構成される空間情報であり、この空間情報は、超音波走査による脊椎表面、深部筋及び脂肪組織の超音波画像情報と空間位置情報を含み、リアルタイムで動的に表示され得るものである。好ましくは、前記立体的な表面様子は脊椎表面の空間情報に基づいて形成される見掛け上の三次元画像であり、前記見掛け上の三次元画像中のグルーヴは、前記すべての超音波画像情報に含まれる脊椎表面のテクスチャである。

【0064】

好ましくは、前記デジタル医用画像は、ＣＴボリュームレンダリング再構成画像や磁気共鳴ＭＲ、コンピュータＸ線撮影ＣＲ又はデジタルＸ線撮影ＤＲによる画像であってもよい。

【0065】

さらには、前記ステップ１）は超音波画像走査と超音波画像情報抽出とを含む。

【0066】

超音波画像走査では、患者の空間位置確認情報を有する超音波走査機器により患者体表を繰り返し走査し、筋骨格境界が識別できたら、生データを取得し、前記生データは、患者の空間位置確認情報を有する超音波走査機器から発する無線周波信号であり、前記走査は各種の異なる角度や方向で行われてもよく、前記方向は上下、前後、左右、斜方向を含み、前記角度は患者体表に垂直又は任意の傾斜角を含み、走査順とは無関係である。好ましくは、位置確認タグを有する超音波プローブで走査して生データを取得し、前記生データは、超音波プローブから発する無線周波信号である。

【0067】

超音波画像情報抽出では、超音波画像走査モジュールで取得された生データから抽出を行い、手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を取得し、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とし、前記手術中の二次元空間位置パラメータとは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータを言う。

【0068】

さらには、前記ステップ２）は超音波画像最適化、超音波画像重ね合わせ及び超音波拓本生成を含む。

【0069】

超音波画像最適化では、ステップ１）で得られた患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を最適化させることで、干渉ノイズが除去され最適化した超音波画像を得る。好ましくは、超音波画像最適化に超音波フィルタ強調技術を採用し、干渉ノイズが除去され最適化した超音波画像を得る。さらに好ましくは、前記超音波フィルタ強調技術は中央値フィルタと数学的モルフォロジーフィルタとを組み合わせたものであってもよい。

【0070】

超音波画像重ね合わせでは、超音波画像最適化により得られた最適化超音波画像を重ね合わせ、超音波骨エコー情報における強エコーと弱エコーの間の相違を一層大きくすることで、画像効果の強調した重ね合わせ超音波画像を得る。好ましくは、前記超音波画像重ね合わせに単一モダリティレジストレーション法を採用する。

【0071】

超音波拓本生成では、超音波画像重ね合わせを行った重ね合わせ超音波画像について処理して超音波拓本を形成し、好ましくは、処理の方法は三次元再構成である。さらに好ましくは、前記三次元再構成の方法は体積レンダリング法であり、体積レンダリング加速技術であることが好ましい。

【0072】

さらには、前記ステップ３）に記載の輪郭マッチング方法は、画素階調アルゴリズムのうちの相互情報アルゴリズムに基づくマルチモダリティレジストレーション法が好ましく、前記マルチモダリティレジストレーション法は画像の空間変換、階調補間、類似度評価、検索最適化を含む。

【0073】

本発明による超音波拓本技術に基づく脊椎画像生成システムに基づいて、超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムをさらに提供し、前記脊柱手術用のナビゲーション・位置確認システムは、ナビゲーションモジュールと超音波拓本技術に基づく脊椎画像生成システムとを含み、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得て、この脊椎表面トポグラフィーに基づいてリアルタイムな手術中ナビゲーションを行うことができる。

【0074】

前記超音波拓本技術に基づく脊椎画像生成システムでは、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像について処理して超音波拓本を形成し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを形成し、前記超音波拓本は、患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像である。好ましくは、前記デジタル医用画像はＣＴボリュームレンダリング再構成画像や磁気共鳴ＭＲ、コンピュータＸ線撮影ＣＲ、デジタルＸ線撮影ＤＲによる画像である。

【0075】

前記ナビゲーションモジュールは、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに含まれる空間位置情報に基づいてトポグラフィー座標系を確立し、位置確認機能を有する手術器械のソリッド幾何モデルを表面モデルへ変換させ、そして脊椎画像生成システムで得られた患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに示し、また、患者の手術中体位の空間位置情報から患者座標系を確立し、トポグラフィー座標系と患者座標系における対応点間の関係を特定してから統一座標系にし、また、変換された手術器械の表面モデルを統一座標系に組み込んで、リアルタイムな手術中ナビゲーションが実現され、手術操作を行うように手術者を指導する。

【0076】

さらには、前記位置確認機能を有する手術器械は、磁気的な空間位置確認タグを有する手術器械である。

【0077】

本発明による超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムに基づいて、超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認方法をさらに提供する。前記方法は、ステップ１）～ステップ３）を含み、
ステップ１）では、手術器械のソリッド幾何モデルを確立し表面モデルへ変換させてから、脊椎画像生成システムで得られた患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに示すことで、手術器械に対するリアルタイムな追従が実現され、前記ソリッド幾何モデルの表面モデルへの変換方法は三次元ＣＡＤモデリングが好ましく、
ステップ２）では、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに含まれる空間位置情報に基づいてトポグラフィー座標系を確立し、患者の手術中体位の空間位置情報から患者座標系を確立し、トポグラフィー座標系と患者座標系における対応点間の関係を特定してから統一座標系にし、前記したトポグラフィー座標系と患者座標系における対応点間の関係を特定する方法は反復最近接点アルゴリズムが好ましく、
ステップ３）では、ステップ１）で得られた手術器械の表面モデルをステップ２）で形成された統一座標系に表示し、手術器械と患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーとが統一座標系にあるようになり、リアルタイムな手術中ナビゲーションが実現される。

【0078】

脊柱手術操作システムをさらに提供し、前記脊柱手術操作システムは、前記いずれかの超音波拓本技術に基づく脊椎画像生成システム及び／又は前記いずれかの超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムを含む。この操作システムは脊柱部位の開放・低侵襲手術に適用され、椎間板の突出部分の摘出、腫瘍切除、神経剥離、血腫除去等のみに限られない様々な手術操作を行うことができる。

【0079】

例えば図３に示されるように、本発明にかかる超音波拓本技術に基づく脊椎画像生成システムでは、二次元的な脊椎表面構造の超音波画像に基づいて超音波拓本を生成し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーが得られ、前記システムは収集ユニットと処理ユニットとを含む。

【0080】

前記収集ユニットは、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を取得するものであり、前記処理ユニットは、収集ユニットで取得された患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像について処理して超音波拓本を形成し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得るものである。前記超音波拓本は、患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像である。前記超音波画像は、脊椎領域の深部筋組織と脊椎表面とを含むすべてのエコー情報から構成される。前記超音波画像は、その輪郭エッジと輪郭内部の情報を含む。前記脊椎表面は横突起と棘突起とを含み、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格のいずれか又は任意の複数種をさらに含んでもよい。前記脊椎表面トポグラフィーは、超音波拓本が有するリアルタイムな動的空間情報をデジタル医用画像に付与してなる立体的な表面様子である。前記超音波拓本が有するリアルタイムな動的空間情報とは、脊椎表面と筋骨格境界とともに構成される空間情報を言い、この空間情報は、超音波走査による脊椎表面、深部筋及び脂肪組織の超音波画像情報と空間位置情報を含み、リアルタイムで動的に表示され得るものである。前記立体的な表面様子は、脊椎表面の空間情報に基づいて形成される見掛け上の三次元画像であり、前記見掛け上の三次元画像中のグルーヴとは、前記すべての超音波画像情報に含まれる脊椎表面のテクスチャを言う。

【0081】

前記デジタル医用画像とは、数値で表すことのできる画像情報を言い、コンピュータにより画像の記憶、再構成、計測、識別、処理を行う。前記デジタル医用画像は、ＣＴボリュームレンダリング再構成画像や磁気共鳴ＭＲ、コンピュータＸ線撮影ＣＲ又はデジタルＸ線撮影ＤＲによる画像であってもよい。

【0082】

前記手術中体位は手術中の患者がなす体位であり、主に仰臥位、腹臥位又は側臥位の一つが採用されるが、手術の需要を満たすための特殊体位もある。前記脊椎領域の深部筋組織は臨床解剖学的な概念であり、具体的には、人体背部の脊柱起立筋及びその深部構造を言う。前記脊椎表面とは、人体脊椎の体表に近い側の骨表面を言う。前記の棘突起、横突起、関節突起、椎板、椎板間隙、椎間孔は医学的概念であり、具体的には、人体脊柱を構成する各々の椎体における解剖学的意味のある骨性構造を言う。背景技術に述べたように、従来の超音波走査機器や処理技術では、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分のすべての超音波画像情報を得ることあるいは正確に得ることが難しかったが、本発明のシステムによれば、脊柱手術に必要な関連部位の正確な空間位置情報を取得することが可能にある。

【0083】

前記収集ユニットは、超音波画像走査モジュールと超音波画像情報抽出モジュールとを含み、前記超音波画像走査モジュールは、走査を行って生データを取得するものであり、前記超音波画像情報抽出モジュールは、前記超音波画像走査モジュールによる走査で取得された生データに含まれる手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を抽出し、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とするものであり、前記手術中の二次元空間位置パラメータとは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータを言う。前記生データとは、患者の空間位置確認情報を有する超音波走査機器から発する無線周波信号を言い、前記走査は各種の異なる角度や方向で行われてもよく、前記方向は上下、前後、左右、斜方向を含み、前記角度は患者体表に垂直又は任意の傾斜角を含み、走査順とは無関係である。

【0084】

超音波画像の取得方法としては、患者の空間位置確認情報を有する超音波走査機器により患者体表を繰り返し走査し、筋骨格境界が識別できたら、手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を抽出し、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とするようにしてもよく、前記手術中の二次元空間位置パラメータとは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータを言う。

【0085】

前記患者の空間位置確認情報を有する超音波走査機器は位置確認タグを有する超音波プローブであり、前記生データとは、超音波プローブから発する無線周波信号を言う。

【0086】

前記処理ユニットは、超音波画像最適化モジュール、超音波画像重ね合わせモジュール、超音波拓本生成モジュール及び画像輪郭マッチングモジュールを含む。前記超音波画像最適化モジュールは、前記収集ユニットで得られた患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を最適化させることで、干渉ノイズが除去され最適化した超音波画像を得るものである。前記超音波画像重ね合わせモジュールは、前記超音波画像最適化モジュールで得られた最適化超音波画像を重ね合わせ、超音波骨エコー情報における強エコーと弱エコーの間の相違を一層大きくすることで、画像効果の強調した重ね合わせ超音波画像を得るものである。前記超音波拓本生成モジュールは、前記超音波画像重ね合わせモジュールで得られた重ね合わせ超音波画像について処理することで、超音波拓本を形成するものである。

【0087】

超音波拓本の取得方法としては、前記患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を最適化させ重ね合わせてから、さらに処理して超音波拓本を生成するようにしてもよい。

【0088】

前記超音波画像最適化モジュールは、超音波フィルタ強調技術により前記収集ユニットで取得された超音波画像を最適化させることで、干渉ノイズが除去され最適化した超音波画像が得られ、前記フィルタ強調は、信号中の特定周波帯に属する周波数をろ過して除去する操作であり、干渉を抑制し防止する重要な措置となる。

【0089】

前記超音波フィルタ強調技術は、適応中央値フィルタと数学的モルフォロジーフィルタとを組み合わせた方法であってもよく、具体的には、まず適応中央値フィルタにより画像前処理を行うことで、スペックルノイズを抑制しながら必要な細部を保ち、そして数学的モルフォロジー方法により画像について二次フィルタとコントラスト強調を行うことで、一層ノイズをろ過し除去しながら境界の細部を保つようにする。

【0090】

前記中央値フィルタは非線形フィルタであり、画像中のインパルスノイズやごま塩ノイズを除去可能で、画像エッジのぼけが生じないように良好に保護でき、二値画像や階調画像の処理に利用され得る。その基本原理は、デジタル画像中のある点の値を、この点の一つの近傍にある各点の値の中央値に置き換えることである。その基本的な技術思想としては、ウインドウ内のすべての画素階調をソートし、ソート結果の中央値を元のウインドウの中央における画素の階調とすることである。中央値フィルタは、周りの画素値が大きく相違している値を周りの画素値に近いものである中央値に変更することで、特異なノイズ点を除去することをその主要機能としている。前記適応中央値フィルタによれば、超音波画像中のスペックルノイズを抑制しながら境界の細部を良好に保護できるが、計算量が大きくて、並列性を持つアルゴリズムではないから、速度が遅くなり、リアルタイム性の要求を満たすことが困難であり、その応用は超音波画像の後処理に限られている。適応中央値フィルタは、一定領域内の画像の統計的特性を考慮して中央値フィルタを改良したものであり、そのろ波手法としては、ウインドウを使用し、ろ波中に画像の局所の統計的特性に基づいて自動的にウインドウの大きさや形状を変え、ウインドウ中央の画素がノイズであると判断した場合、その値を中央値に置き換え、そうでない場合、現在の画素値を変更しないようにする。ウインドウ内の各画素の重みは調節可能なものであり、ウインドウの選択は分割や統合等の方法により行われてもよく、形状が制限されておらず、そして、適応フィルタの設計の技術思想は以下のとおりである。

【0091】

画素毎の局所領域での平均値μと推定分散σ^２は以下の式となる。

【数1】

【0092】

σ^２をノイズ分散とし、すべての局所領域の推定分散の平均値に置き換えれば、適応フィルタの推定式は以下の式となる。

【数2】

式中、ηは画像中のＭ＊Ｎ画素からなる矩形の局所領域であり、ＭとＮのそれぞれはこの近傍画素行列の行数と列数であり、ｎ_１、ｎ_２は画素座標であり、ａ（ｎ_１、ｎ_２）は元画像中の（ｎ_１、ｎ_２）画素の階調値であり、ｂ（ｎ_１、ｎ_２）は画像の適応フィルタ後の（ｎ_１、ｎ_２）画素の階調値である。

【0093】

前記数学的モルフォロジーフィルタは幾何構造の領域フィルタを行うものであり、即ち、一つのノイズブロックでは一つの点が判定条件を満たしていれば、領域全体のノイズをすべて除去できるものである。このため、この方法によれば、モルフォロジーフィルタの品質を保証でき、ノイズがよく除去可能で、画像の細部を良好に保つことができるが、満足できるフィルタ効果を達成するには、適切な構造要素を選択する必要があり、要求の高いリアルタイム性を満たし得る。数学的モルフォロジーフィルタは新規な非線形フィルタであり、幾何学を基礎として画像の幾何構造を検討することに着目するものである。その基本思想は、画像を集合と見なして、予め定義された構造要素により画像を検知し、構造要素を良好に画像内部に入れ得るかどうかを観測するとともに、構造要素を入れる方法が正確であるかどうかを検証する。図中における構造要素を入れるのに適する位置をマークすれば、画像に関する構造情報が得られ、そして画像の構造情報について解析することで、次に選択すべき構造要素を判定できる。そうすれば、より効率よくろ波によって興味のある情報をさらに取得できる。このように最適な構造要素を見つけ出して現在最も望ましい情報を抽出し、興味のない情報を抑えることで、すべてのノイズをろ過し除去して有用な情報が残る。前記構造要素は一定の寸法を有する背景画像であり、大きさや形状、階調、色度等の情報を持ち、その寸法が検討される目標よりもはるかに小さく、画像において変換されてフィルタウインドウに類似するものとなり、構造要素が適切に選択されているか否かは入力画像についての処理結果に直接影響しており、形態変換は形態フィルタを構築する基礎であり、基本的な形態変換は膨張、浸食、開演算及び閉演算を含み、それぞれ以下のように定義されている。

【数3】

【0094】

上式では、ｆは入力画像、ｂは構造要素を示し、前記構造要素とは、情報を収集する「プローブ」を言い、画像を検知するものであり、画像において構造要素を移動させていくことで、構造要素と画像対象との相互作用により画像の各部分間の関係を検討し、さらに画像処理に必要な幾何パラメータを判定することができる。Ｄ_ｆとＤ_ｂはそれぞれｆとｂの定義域である。実際の応用では、開演算は入力画像を平滑化させるものであり、入力画像のエッジのばりや孤立したスペックル、即ち、ろ過信号中の正パルスノイズを除去でき、閉演算は入力画像をろ過するものであり、入力画像のクラックやホールを埋め、つまり信号中の負パルスのノイズをろ過することができる。

【0095】

上記した適応中央値フィルタと数学的モルフォロジーフィルタとを組み合わせた方法によれば、本モジュールでは、収集ユニットで得られた患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を入力し、超音波フィルタ強調により最適化させ、最終的に、干渉ノイズが除去され最適化した超音波画像を出力するようになる。

【0096】

前記超音波画像重ね合わせモジュールは、前記超音波画像最適化モジュールで得られた最適化超音波画像を重ね合わせ、超音波骨エコー情報における強エコーと弱エコーの間の相違を一層大きくすることで、画像効果を強調するものである。前記超音波画像重ね合わせに単一モダリティレジストレーション法を採用してもよい。画像レジストレーションは、一方の医学画像について一種又は一連の空間変換を行うことで、他方の医学画像における対応点と空間的に一致させることであり、レジストレーションの結果、両方の画像におけるすべての解剖点、あるいは、少なくともすべての診断意味のある点及び手術上の興味のある点がいずれもマッチングするようになるべきであり、前記の空間的に一致させることは、人体における同一解剖点が二つのマッチング画像において同様な空間位置にあるようにすることを言う。単一モダリティレジストレーションでは、同一物体は画像取得条件の相違又は物体自身に発生する空間位置変化により、結像に相違が発生することがある。この時、レジストレーション待ち画像間の画素の階調分布はいつも一定の線形関係にあるから、相関係数法によりレジストレーションを行うと良く、前記相関係数（ＣＣ）は二つの変数の線形相関の近さを記述する量であり、以下のように定義されている。

【数4】

【0097】

ただし、

【数5】

はｆ（ｘ_ｎ）、ｇ（ｘ_ｎ）の平均値であり、医学画像の単一モダリティレジストレーションでは、同一結像方式（超音波）で同一組織を結像させれば、相関性を計算する時に誤りが生じることがないため、相関係数法は医学画像の単一モダリティレジストレーションに適している。

【0098】

上記処理により、本モジュールでは、超音波画像最適化モジュールで得られた干渉ノイズが除去され最適化した超音波画像を入力し重ね合わせていくことに従って、重ね合わせた画像において超音波骨エコー情報の強エコーと弱エコーの間の相違が大きくなり、それにより、出力された重ね合わせ超音波画像は効果が強調され、相違が明らかになり、干渉ノイズが除去されている。前記線形の強エコーは医用超音波学概念で、超音波画像を特徴的に記述するものであり、線形のハイライト画像として現れる。

【0099】

前記超音波拓本生成モジュールは、前記超音波画像重ね合わせモジュールで得られた重ね合わせ超音波画像について処理することで、超音波拓本を形成するものである。前記超音波拓本は、超音波画像重ね合わせモジュールで得られた超音波画像について三次元再構成処理を行って得られたものであり、患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像である。超音波画像の三次元再構成は、一定のデータ処理方法により、超音波プローブ走査による二次元画像を三次元画像に変換するようにしてもよく、一連の的二次元画像から三次元画像モデルを再構成し、定性・定量分析を行う技術であり、二次元画像から三次元構造情報を抽出することは実現されている。三次元再構成のレンダリングにおいて、データ記述方法には、主に表面レンダリング法、体積レンダリング法及びハイブリッドレンダリング法の三つの種類がある。前記超音波拓本生成モジュールは、体積レンダリング法及び体積レンダリング加速技術により、脊椎エコーの重ね合わせによる画像について三次元再構成を行うことで、二次元的な脊椎表面構造の超音波画像の超音波拓本への変換を実現することができ、正確性と低コスト化の要求が達成される。体積レンダリング法は、光線が三次元体のデータ場を通過する場合のボクセルとの相互作用を直接検討し、皮膚や骨格、筋等の物質情報を一括して表示し、等値面の輪郭を抽出する必要のない方法である。体積レンダリングは、ボリュームデータの表面情報や内部情報のいずれをも表示でき、主に投影、隠線除去、レンダリング及び統合等を含む。本実施例では、体積レンダリング法中のシアー・ワープ法を採用してもよく、それは、二次元離散データ場の射影変換を二次元データ場せん断変換と二次元画像変形の二つの段階に分けて行うものであり、それにより、二次元空間のリサンプリングが二次元平面のリサンプリングに変換され、計算量が大幅に低減する。シアー・ワープ法によれば、最終画像の品質を著しく低下させることなく、二次元データ場の体積レンダリングをグラフィックスワークステーションにおいてほぼリアルタイムに行うことが可能になる。もちろん、体積レンダリング法はハードウェア要求が高く、多大なメモリを必要とし、計算量が大きくなるという短所もあり、これについては、体積レンダリング加速技術により解決可能である。体積レンダリング加速技術は、画像空間に基づく体積レンダリングアルゴリズム加速、物体空間再構成に基づくアルゴリズム加速、及び、ハードウェア加速を含む。前記画像空間に基づく体積レンダリングアルゴリズム加速では、第一、画像空間の相関性により射線的の数を低減する。隣り合う画素同士は相関性があり、類似した色を有するため、画像平面におけるすべての画素から光線を透過させる必要がなく、間隔的に光線を投影させればよい。第二、対象空間の相関性から不必要なサンプリング点数を低減する。ボリュームデータには数多くの空ボクセルつまり透明ボクセルが存在しており、これらの空ボクセルは体積レンダリングによる最終画像になんら関与せず、そして、ボリュームデータの記憶構造を工夫することで、空ボクセルをスキップし、サンプリング点数を低減することができる。前記物体空間再構成に基づくアルゴリズム加速については、現在、成熟のアルゴリズムとしてボクセル投影法やサブ領域投影法等があり、また、特殊なボリュームデータ構造、例えばソート済みボリュームデータ、ランレングス符号化等を確立することによって、ボリュームデータに対するトラバースを早くすることもできる。前記ハードウェア加速については、超音波拓本の結像範囲は全脊柱セグメントのみにされており、中規模のボリュームデータ場となるため、総合的に考慮すれば、ハードウェア加速を導入したテクスチャレンダリングアルゴリズムが好適であり、このアルゴリズムによれば、他の加速アルゴリズムよりも迅速な体積レンダリングが出来ている。将来、ボリュームデータ場の急速な拡大が必要になる場合、並列計算関連アルゴリズムを導入して、さらなる加速を実現することも可能である。

【0100】

以下、シアー・ワープ法を例（図７）としてその基本原理を説明する。シアー・ワープ法の主な思想は、三次元離散データ場の直接射影を三次元データ場のせん断変換と二次元画像の変形に分けることである。例えば図７に示されるように、まず、視線方向（射影方向）に応じて一つの中間平面を決定し、この平面は物体空間のある軸方向に垂直であり、物体空間をせん断変換させてから垂直に中間平面に射影させることで、中間画像Ａ’Ｂ’が得られる。垂直射影となるから、その計算量は必ず斜射影よりもかなり低くなる。最後に、得られた中間画像Ａ’Ｂ’について変形変換を行うことで、最終画像ＡＢが得られる。元の射影変換のビュー変換Ｍ_ｖｉｅｗを二つの行列、つまりズレ行列Ｍ_{ｓｈｅａｒ}と変形Ｍ_ｗａｒｐ行列に分けている。従って、Ｍ_ｖｉｅｗ＝Ｍ_{ｓｈｅａｒ}×Ｍ_ｗａｒｐ。さらに、平行射影については、一組の水平線で三次元離散データ場の一つの横断面を示し、観測方向を座標系の一つの軸と平行にするとともに、三次元データ場の各データ平面と垂直にすると、せん断変換では、各データ平面はある距離を平行移動するようになり、この距離は各データ平面が所在する座標値に比例するだけでなく、初期観測方向に対する定義にもよるものであり、また、透視射影については、三次元離散データ場がせん断物体空間に変換される場合、データ平面毎は平行移動変換ばかりかスケール変換をも行う必要になり、このため、せん断変換したデータ場の射影による画像は最終画像ではなく、単に中間画像になり、さらなる画像変形は必要である。

【0101】

上記した三次元離散データ場のせん断変換及び二次元画像変形は、以下のような観測変換行列に統合されてもよい。

【数6】

【0102】

式中、Ｐは座標変換行列であり、それにより、座標系のＺ軸と三次元離散データ場の主要観測方向は重なり合うようになり、Ｓは三次元データ場をせん断物体空間に変換させるもので、Ｍ_ｗａｒｐはせん断物体空間を画像空間に変換させるものである。平行射影については、せん断はＺ軸に垂直な方向に行われるため、以下の式が出来る。

【数7】

【0103】

そのうち、Ｓ_ｘ、Ｓ_ｙはＭ_ｖｉｅｗのパラメータから求まる。透視射影については、以下の式が出来る。

【数8】

【0104】

上式から分かるように、ｚ＝ｚ_０のデータ平面を物体空間からせん断物体空間に変換させた場合、このデータ平面は（Ｚ_０Ｓ’_ｘ，Ｚ_０Ｓ’_ｙ）の距離を平行移動し、係数 １／（１＋Ｚ_０Ｓ’_ｗ）に応じてスケール変換を行うべきである。上式における最後の項のＭ_ｗａｒｐはせん断物体空間を画像空間に変換させる行列である。

【数9】

【0105】

上記アルゴリズムによれば、前記超音波画像重ね合わせモジュールで得られた超音波画像について処理し、つまり、Ｚ軸方向を三次元離散データ場の主要観測方向とするように座標系の軸と重なり合わせ、それと平行な部分について平行射影を行い、脊椎表面のすべての骨エコー情報の平行移動のみが行われ、透視射影については、上記平行移動が完了した後、さらにスケール変換を行う。前記変換は、せん断物体空間で、脊椎表面のすべての骨エコー情報をせん断してからせん断物体空間の中間画像平面に投影させることで、中間画像が形成され、そして中間画像平面での画像を二次元画像変形Ｍ_ｗａｒｐにより画像空間に変換させて、最終画像が得られることである。

【0106】

上記によれば、重ね合わせ超音波画像は、体積レンダリング法及び加速技術により超音波拓本になり、患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像として次の処理に用いられる。

【0107】

前記画像輪郭マッチングモジュールは、前記超音波拓本生成モジュールで得られた超音波拓本をデジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを形成するものである。前記超音波拓本が有するリアルタイムな動的空間情報とは、脊椎表面と筋骨格境界とともに構成される空間情報を言い、この空間情報は、脊椎表面と深部筋組織の超音波画像情報と空間位置情報を含み、リアルタイムで動的に表示され得るものである。前記脊椎表面トポグラフィーは、超音波拓本が有するリアルタイムな動的空間情報をデジタル医用画像に付与してなる立体的な表面様子であり、前記立体的な表面様子は、脊椎表面の空間情報に基づいて形成される見掛け上の三次元画像であり、前記見掛け上の三次元画像中のグルーヴとは、前記すべての超音波画像情報に含まれる脊椎表面のテクスチャを言う。前記輪郭マッチングとは、前記超音波拓本に含まれる脊椎表面のすべての骨エコー情報をデジタル医用画像の輪郭と一点ずつマッチングさせることを言う。以下、ＣＴボリュームレンダリング再構成画像を例として輪郭マッチング工程について説明する。超音波拓本の輪郭に含まれるリアルタイムな動的空間位置情報とＣＴボリュームレンダリング再構成画像の輪郭に含まれる空間位置情報とについて、対応点一つ毎にマッチングさせることで、画像輪郭のマッチングが実現され、さらに三次元レベルでのマッチングが実現される。

【0108】

前記輪郭マッチングの実施方法は、前記超音波拓本に含まれる脊椎表面のすべての骨エコー情報をデジタル医用画像の輪郭と一点ずつマッチングさせることで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得るようにしてもよい。前記輪郭マッチングの方法は、画素階調アルゴリズムのうちの相互情報アルゴリズムに基づくマルチモダリティレジストレーション法が好ましく、前記マルチモダリティレジストレーション法は画像の空間変換、階調補間、類似度評価、検索最適化を含む。本実施例で採用される輪郭マッチングは、画素階調アルゴリズムのうちの相互情報アルゴリズムに基づくマルチモダリティレジストレーション法が好ましく、前記マルチモダリティレジストレーション法は、脊椎画像の空間変換、階調補間、類似度評価、検索最適化を含む。前記相互情報アルゴリズムは、レジストレーション待ち画像の独立エントロピー及び結合エントロピーを計算することで得られた相互情報量に基づいて、画像レジストレーションの効果を判断する。つまり、二つの画像の重なり合い度合いが最高になる場合、結合エントロピーが最も小さく、相互情報量が最も大きくなり、最良のレジストレーション効果が達成され、また、相互情報アルゴリズムは画像の階調の統計情報を利用しているから、ノイズによる干渉が少なく、良好なロバスト性を有する。前記ロバスト性とは、一定の構造や大きさのパラメータによる摂動でも、制御システムは何らかの性能を維持する特性を言う。そして、超音波結像とＣＴボリュームレンダリング結像は、その結像方式と原理が異なるが、画像における対応位置での階調分布は類似しており、異なるモダリティの画像のそれぞれの密度変化が線形関係にあるかどうかは相互情報を解くことに関係なく、また、脊柱の各々の椎体は剛体に近いものと見なすことができるため、脊椎画像マッチングは剛体レジストレーションに属し、剛体は歪みが生じにくい特性があり、相互情報アルゴリズムによる剛体レジストレーションの効果がより良好になる。

【0109】

画像マッチング工程は、主に画像の空間変換、階調補間、類似度評価、検索最適化を含み、この工程では、空間変換策、階調補間アルゴリズム、類似性測定方法及び最適化アルゴリズムを選択する必要がある。前記画像の空間変換とは、一定の規則で画像ｆ（ｘ，ｙ）中の画素点を新しい位置に変換させることで、新しい画像ｆ（ｘ’，ｙ’、）が得られることを言い、即ち、

【数10】

であり、画像の空間変換方式の自由度によっては、剛体変換、アフィン変換及び非線形変換の三つに分けられている。好ましくは、画像の特徴に応じて異なる方式を選択し、それらを組み合わせて全体の画像変換を実現する。脊柱椎体の骨表面については、剛体変換とアフィン変換の組合せが採用され、剛体変換は画像の平行移動と回転を含み、画像中の物体の位置と方向のみが変更され、すべての物体の長さや角度、面積及び体積はいずれも変わらない。このような変換は一つの３×３の行列で表すことができ、つまり、変換行列Ｔにより画像中の元の座標（ｘ，ｙ）の点を新しい座標に変換させ、以下の表現式となる。

【数11】

【0110】

ただし、θは画像の回転角度であり、ｄｙ、ｄｘは画像のｚ軸方向とｙ軸方向での平行移動量である。剛体変換では、画像中の物体の相対位置は変化していない。前記アフィン変換は画像の真直度と平行度を保持する空間変換である。前記真直度とは、変換後の画像中の直線が湾曲せず、円弧の弧度も変わらないことを言う。前記平行度とは、画像中の物体間の平行関係が変わらないが、物体の大きさや物体間の角度が変わり得ることを言う。剛体変換に対して、前記アフィン変換には変換自由度つまり各座標方向におけるズーム係数が加えられ、各方向のズーム比が一致する場合、一様アフィン変換と称され、そうでない場合、非一様アフィン変換又はせん断変換と称されるようになる。一様アフィン変換はレンズ結像システムに用いられることが多く、このシステムによる画像の大きさは撮影中の物体からカメラまでの距離と密接に関係しており、せん断変換は機器に起因した画像の歪みを校正するために用いられることが多い。アフィン変換は以下の式で表される。

【数12】

【0111】

マルチモダリティ画像レジストレーションでは、アフィン変換は、結像方法の違いによる画像中の物体の大きさの不一致を直すことや、操作上の人的ミスによる画像の歪みを直すこと、及び、結像角度の違いによる画像表現の違いをなくすことに用いられる。

【0112】

前記非線形変換について、人体軟組織の画像では、異なる時刻や角度で撮影された組織は捻られたり歪みが発生したりことが多いため、レジストレーションにて非線形変換がよく採用されている。超音波走査において、プローブと接触する人体組織は小さな局所歪みが発生するため、弾性モデルに適用される。前記弾性モデルとは、浮動画像からテンプレート画像への変換を弾性材料の伸縮過程と見なすことができ、外部からの材料への外力や材料内部からの反発による内力はこの過程に亘って発生し、外力により材料に歪みが生じる。また、内力は外力の受け入れに応じた材料からのフィードバックであり、内力と外力がバランスを取ると、レジストレーションが終了する。その歪みは、以下のＮａｖｉｅｒ偏微分方程式で表すことができる。

【数13】

【0113】

ただし、μ（ｘ，ｙ）は物体のｘ方向とｙ方向での歪みを示し、λとψは物体の弾性性質を記述するＬａｍｅ弾性定数であり、物体に作用する二つの方向での外力を示す。弾性モデルでは、弾力は歪みの大きさに比例して変化し、局所歪みが大きくなるほど、弾力が大きくなる。

【0114】

前記階調補間アルゴリズムについて、マッチング工程が完了し、浮動画像について空間変換を行った後、整数座標にあった画素は新しい画像において非整数座標になる可能性があり、その問題を解決し、歪んだ画像をより良好に記述するために、階調補間アルゴリズムを採用する。レジストレーション工程において、浮動画像について一定の空間変換を行うことで、可能な限りテンプレート画像に類似させる必要があり、変換後、一般的には、新しい画像中の整数座標点について逆変換を行うことで、この点の元画像での対応座標が得られ、得られた点は補間点となる。階調補間アルゴリズムでは、補間点及びその周りの点の階調値の計算により新しい画像における座標の階調値を得ることができるため、画像の元の情報を最大限保持する。相互情報からレジストレーションを行うアルゴリズムでは、画像の階調分布情報を十分に保持するために、ＰＶ補間アルゴリズムを使用することが多い。前記ＰＶ補間アルゴリズムは、補間により新しい画像の階調分布を直接得るものであり、空間変換前の浮動画像中の画素毎は画像の階調分布にある程度関与しているため、ＰＶ補間は、すべての補間点及びその周りの近傍画素の階調値を考慮したものとなる。一般に、階調ヒストグラムの計算において、階調値ｎ_０が出現すると、統計する度に、ヒストグラムにおけるガンマに対応する数値に１が加算されるようになる。双線形補間と同様に、新しい画像中の画素点（ｘ，ｙ）毎は元画像中の一つの補間点（ｘ’，ｙ’）に対応しており、補間点周りの四つの整数点の階調値をｎ_１１、ｎ_１２、ｎ_２１、ｎ_２２とし、補間点からこれらの点までの距離に基づいて各点の重みｗ_ｉｊを求める。そして、新しい画像の階調ヒストグラムｈ（ａ）において、その分、階調値ｎ_ｉｊの出現回数にｗ_ｉｊを加算する。

【数14】

【0115】

テンプレート画像中の階調をｍとすると、補間点のジョイント階調ヒストグラムＨ（ａ，ｂ）による計算は以下の式となる。

【数15】

【0116】

すべての補間点について計算し終えると、新しい階調ヒストグラムが得られ、それにより、浮動画像の変換後のマージン階調分布とジョイント階調分布が得られる。ＰＶ補間アルゴリズムでは、直接補間する代わりに、画像の階調統計情報を直接変えることで新しい画像を生成するため、相互情報のような階調分布に非常に敏感なレジストレーションアルゴリズムに適している。

【0117】

前記類似度評価とは、医学画像レジストレーション、特にマルチモダリティ医学画像のレジストレーションでは、二つの画像間の類似度を評価する標準を定義する必要があることを言い、二つの画像の類似度が最大になると、二つの画像の最良のレジストレーション効果が達成されるとする。好ましくは、上記画像の空間変換と階調補間アルゴリズムが完了した後、画素階調アルゴリズムのうちの相互情報アルゴリズムに基づいて類似度の計測を行う。階調値のもとに、二つのレジストレーション待ち画像ＡとＢについて、その対応点の階調値特徴同士を階調ペアとし、相互情報の計算は、画像Ａに画像Ｂの情報がどれだけ含まれているかを計算することになる。二つの画像の最良のマッチングが達成されると、それらの階調分布の相関性が最高になる。

【0118】

従来の相互情報アルゴリズムは画像レジストレーションにおいて大きな欠陥がある。つまり、画像に存在する空間情報が無視されているため、相互情報から二つの画像の類似度を評価する場合に、二つの画像間に共通である空間関係が無視されることで、レジストレーションの誤りが生じる場合がある。従来のアルゴリズムを改良しても、膨大な計算量により計算時間が長くなるため、空間情報を相互情報計算に導入する方法、即ち、条件相互情報を提案した。前記条件相互情報は、レジストレーション待ち画像を異なる小領域に区画し、二つの画像における対応する小領域のすべてのレジストレーション効果のそれぞれにより、画像全体のレジストレーション効果が達成されることであり、演算量が増加することなく、画像中の空間情報を考慮した。条件相互情報ＣＭＩは、ジョイントヒストグラムを拡張し、画像中の小領域の空間分布と階調分布とを結合したものである。

【数16】

【0119】

ただし、ｐ（ｓ）は、レジストレーションの時に画像中の領域ｓを選択してレジストレーションを行う確率であり、レジストレーション待ち画像を同じ大きさの領域に区画すれば、この確率は同じになる。一つの画像では、階調値ｉの画素点は画像の異なる領域に分布している可能性がある。条件相互情報の計算において、異なる領域に分布している階調値ｉの画素点による相互情報量全体への寄与はその点の空間位置に関係している。条件相互情報は、領域ｓにおける画像Ａに含まれる画像Ｂの情報量を示しており、すべての領域の情報量の加重平均値を、画像類似度を評価する最終の基準とする。

【0120】

上記により、本モジュールでは、超音波拓本生成モジュールで得られた超音波拓本が入力されてデジタル医用画像と輪郭マッチングが行われ、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーが出力される。この画像は、手術者による患者の手術中体位での手術に関連する解剖構造の位置関係の観測・分析や、穿刺位置確認への指導及び後続の手術操作を支援するものになる。

【0121】

例えば図４に示されるように、本発明による超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認システムは、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを得て、この脊椎表面トポグラフィーに基づいてリアルタイムな手術中ナビゲーションを行うことができ、前記脊柱手術用のナビゲーション・位置確認システムは、ナビゲーションモジュールと超音波拓本技術に基づく脊椎画像生成システムとを含む。

【0122】

前記超音波拓本技術に基づく脊椎画像生成システムでは、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像について処理して超音波拓本を形成し、デジタル医用画像と輪郭マッチングを行うことで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーを形成し、前記超音波拓本は、患者の空間位置確認情報を有しながら患者の位置変動と共にリアルタイムで動的に更新される個人化の三次元的な超音波骨格画像である。前記デジタル医用画像は、ＣＴボリュームレンダリング再構成画像や磁気共鳴ＭＲ、コンピュータＸ線撮影ＣＲ、デジタルＸ線撮影ＤＲによる画像から選ばれてもよい。

【0123】

前記ナビゲーションモジュールは、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに含まれる空間位置情報に基づいてトポグラフィー座標系を確立し、位置確認機能を有する手術器械のソリッド幾何モデルを確立し表面モデルへ変換させ、そして脊椎画像生成システムで得られた患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに示し、また、患者の手術中体位の空間位置情報から患者座標系を確立し、トポグラフィー座標系と患者座標系における対応点間の関係を特定してから統一座標系にし、また、変換された手術器械の表面モデルを統一座標系に組み込んで、リアルタイムな手術中ナビゲーションが実現され、手術操作を行うように手術者を指導するようになる。前記位置確認機能を有する手術器械は、磁気的な空間位置確認タグを有する手術器械である。

【0124】

例えば図５に示されるように、本発明による超音波拓本技術に基づく脊柱手術用のナビゲーション・位置確認方法は以下のステップがある。

【0125】

ステップ１では、手術者は手術適応症に合致する患者について手術の用意をする。前記手術適応症は、脊柱関連疾患、例えば脊柱腫瘍や椎間板突出を含み、診療規則に規定される標準的なものであり、非手術治療では治癒できず、手術により治療を促進できる。

【0126】

ステップ２では、手術治療を予定している患者に対して手術前ＣＴ走査を行い、それによる画像に対してボリュームレンダリング再構成を行うことで、この患者のデジタル医用画像が得られる。ここでＣＴ走査を行う時、患者の体位を手術にて採用する予定の体位と一致させるべきである。例えば、腹臥位で手術を行う予定の患者については、手術前ＣＴ走査の時にも腹臥位で走査を行うようにする。

【0127】

ステップ３では、患者のＣＴ走査とボリュームレンダリング再構成画像に基づいて、手術チームは手術前分析を行うことで、手術標的部位が特定される。前記手術標的部位とは、手術にて到達しようとする脊椎の骨性の解剖学的位置を言い、例えば、脊柱内視鏡を用いて椎板後部アプローチ髄核摘出術により腰椎椎間板ヘルニアを治療する場合、手術標的部位は目標腰椎の椎板構造となる。

【0128】

ステップ４では、通常の手術準備及び麻酔を行った後、手術者は手術目標領域に対して超音波走査を行う。つまり、患者の空間位置確認情報を有する超音波走査機器により患者体表を繰り返し走査し、筋骨格境界が識別できたら、手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を抽出して、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とする。前記手術中の二次元空間位置パラメータとは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータを言う。前記患者の空間位置確認情報を有する超音波走査機器は、位置確認タグを有する超音波プローブであってもよい。前記手術目標領域とは、手術が行われる関連領域を言い、手術が行われる具体的な部位及びその隣接構造を含む。例えば、脊柱内視鏡を用いて腰椎椎間板ヘルニア手術を行う場合、目標領域は、病変髄核の所在する腰椎椎体、隣り合う上下の椎体、及び、対応する腰椎の近傍領域となる。前記超音波画像は、脊椎領域の深部筋組織と脊椎表面とを含むすべてのエコー情報から構成され、前記超音波画像は、前記超音波画像輪郭と超音波画像エッジの情報を含む。前記脊椎表面は、横突起と棘突起とを含み、関節突起、椎板、椎板間隙、椎間孔及び脊椎を構成する他の骨格部分的のいずれか又は任意の複数種をさらに含んでもよい。

【0129】

ステップ４は、超音波画像走査と超音波画像情報抽出とを含んでもよい。

【0130】

超音波画像走査では、患者の空間位置確認情報を有する超音波走査機器により患者体表を繰り返し走査し、筋骨格境界が識別できたら、生データを取得する。前記生データとは、患者の空間位置確認情報を有する超音波走査機器から発する無線周波信号を言う。前記走査は各種の異なる角度や方向で行われてもよく、前記方向は上下、前後、左右、斜方向を含み、前記角度は患者体表に垂直又は任意の傾斜角を含み、走査順とは無関係である。患者の空間位置確認情報を有する超音波走査機器は位置確認タグを有する超音波プローブであってもよく、走査により生データが取得され、前記生データとは、超音波プローブから発する無線周波信号を言う。

【0131】

超音波画像情報抽出では、超音波画像走査により取得された生データから手術中の二次元空間位置パラメータを含む脊椎表面のすべての骨エコー情報を抽出して、患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像とする。前記手術中の二次元空間位置パラメータとは、患者の手術状況での脊椎のリアルタイムな二次元空間位置パラメータを言う。

【0132】

ステップ５では、前記患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を最適化させ重ね合わせてから、さらに処理して超音波拓本を生成する。

【0133】

ステップ５は、超音波画像最適化、超音波画像重ね合わせ及び超音波拓本生成を含んでもよい。

【0134】

超音波画像最適化では、ステップ４で得られた患者の手術状況でのリアルタイムな手術中体位に対応する二次元的な脊椎表面構造の超音波画像を最適化させることで、干渉ノイズが除去され最適化した超音波画像を得る。好ましくは、超音波画像最適化に超音波フィルタ強調技術を採用し、干渉ノイズが除去され最適化した超音波画像が得られる。さらに好ましくは、前記超音波フィルタ強調技術は、中央値フィルタと数学的モルフォロジーフィルタとを組み合わせたものであってもよい。

【0135】

超音波画像重ね合わせでは、前記超音波画像最適化を行った超音波画像を重ね合わせ、超音波骨エコー情報における強エコーと弱エコーの間の相違を一層大きくすることで、画像を強調し、前記超音波画像重ね合わせに単一モダリティレジストレーション法を採用してもよい。

【0136】

超音波拓本生成では、前記超音波画像最適化・重ね合わせを行った超音波画像について三次元再構成処理を行うことで、超音波拓本を生成し、前記三次元再構成の方法は体積レンダリング法及び体積レンダリング加速技術であってもよい。

【0137】

ステップ６では、ステップ５で得られた前記超音波拓本をデジタル医用画像と輪郭マッチングを行う。つまり、超音波拓本に含まれる脊椎表面のすべての骨エコー情報をデジタル医用画像の輪郭と一点ずつマッチングさせることで、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーが得られる。前記輪郭マッチングの方法は、画素階調アルゴリズムのうちの相互情報アルゴリズムに基づくマルチモダリティレジストレーション法が好ましく、前記マルチモダリティレジストレーション法は、画像の空間変換、階調補間、類似度評価、検索最適化を含む。前記輪郭マッチングとは、前記超音波拓本に含まれる脊椎表面のすべての骨エコー情報をデジタル医用画像の輪郭と一点ずつマッチングさせることを言う。前記脊椎表面トポグラフィーは、超音波拓本が有するリアルタイムな動的空間情報をデジタル医用画像に付与してなる立体表面画像である。前記超音波拓本が有するリアルタイムな動的空間情報とは、脊椎表面と筋骨格境界とともに構成される空間情報を言い、この空間情報は、脊椎表面と深部筋の超音波画像情報と空間位置情報を含み、リアルタイムで動的に表示され得るものであり、前記立体的な表面様子は脊椎表面の空間情報に基づいて形成される見掛け上の三次元画像であり、前記見掛け上の三次元画像中のグルーヴとは、前記すべての超音波画像情報に含まれる脊椎表面のテクスチャを言う。前記リアルタイムな動的空間情報は、超音波走査による脊椎表面、深部筋及び脂肪組織の超音波画像情報と空間位置情報を含み、位置確認タグを有する超音波プローブの移動と共に動的にリアルタイムで表示され得るものである。前記デジタル医用画像はＣＴボリュームレンダリング再構成画像や磁気共鳴ＭＲ、コンピュータＸ線撮影ＣＲ、デジタルＸ線撮影ＤＲによる画像である。

【0138】

ステップ７では、手術器械のソリッド幾何モデルを確立し表面モデルへ変換させてから、脊椎画像生成システムで得られた患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに示すことで、手術器械に対するリアルタイムな追従が実現される。前記ソリッド幾何モデルの表面モデルへの変換方法は三次元ＣＡＤモデリングが好ましい。

【0139】

ステップ８では、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーに含まれる空間位置情報に基づいてトポグラフィー座標系を確立し、患者の手術中体位の空間位置情報から患者座標系を確立し、トポグラフィー座標系と患者座標系における対応点間の関係を特定してから統一座標系にする。前記したトポグラフィー座標系と患者座標系における対応点間の関係を特定する方法は反復最近接点アルゴリズムが好ましい。

【0140】

ステップ９では、変換された手術器械の表面モデルを統一座標系に表示し、手術器械と患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーとが統一座標系にあるようにし、リアルタイムな手術中ナビゲーションが実現される。

【0141】

ステップ１０では、ナビゲーション画像に基づいて、手術者は統一座標系において座標系が表示する手術器械と患者の脊柱解剖構造との位置関係に応じて後続操作を行う。このステップでは、手術者は、患者の手術中体位に合致するリアルタイムで更新される個人化の脊椎表面トポグラフィーによる誘導で、位置確認タグを有する手術器械を用いて両手のみ又はロボットによる補助で位置確認と穿刺を行った後、継続してナビゲーション画像により、画像に示されている手術器械と脊柱の骨性構造との位置関係及び手術者から見えている解剖構造の両者による確認で、後続の手術操作を終える。脊柱内視鏡を用いて椎板後部アプローチ髄核摘出術により腰椎椎間板ヘルニアを治療することを例とすれば、その後続の手術操作としては、手術者は、臨床経験に基づき、棘突起、黄色靱帯、椎板、神経根、椎間板等の解剖構造を識別し、突出した髄核又は線維輪等に到達し摘出を行って関連手術操作が完了するまで、手術者は脊柱手術用のナビゲーション・位置確認システムによる指導で手術器械を用いて正確な角度と距離で操作を行う。

【符号の説明】

【0142】

１ 頚椎、２ 胸椎、３ 腰椎、４ 仙骨、５ 尾骨、６ 上関節突起、７ 棘突起、８ 下関節突起、９ 横突起、１０ 椎体、１１ 乳様突起、１２ 椎孔、１３ 横突孔、１４ 椎板、１５ 横突起と椎体との接続根元

【図1】

【図2】

【図3】

【図4】

【図5】

【図6(1)】

【図6(2)】

【図6(3)】

【図7】