(57)【要約】【課題】形状を予測するためのシステムおよび方法を提供する。【解決手段】該システムは、データベースと、トレーニング解析モジュールと、主体解析モジュールと、予測モジュールとを含むことができる。該データベースは、第１および第２パラメータ（例えば骨および軟骨）によってそれぞれ特徴付けられた２つのトレーニングモデル群、および、第１パラメータによって特徴付けられた主体モデル（例えば骨モデル）を記憶可能である。これらモデル間の関係を、トレーニング解析モジュールおよび主体解析モジュールを用いて決定することができる。予測モジュールは、これらの関係に基づき、第２パラメータによって特徴付けられる予測形状（例えば、骨モデルに対応する予測の軟骨モデル）を生成することができる。２つのトレーニングモデル群は、二次元モデル群と三次元モデル群とであってよい。【選択図】図４Ａ
2019069158.tif 000002
【特許請求の範囲】
【請求項１】
形状を予測するためのシステムであって、
第１パラメータによって特徴付けられた第１のトレーニングモデル群と、第２パラメータによって特徴付けられた第２のトレーニングモデル群と、前記第１パラメータによって特徴付けられた主体モデルと、を記憶するように構成された少なくとも１つのデータベースと、
前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との関係を決定するように構成されたトレーニング解析モジュールと、
前記主体モデルと前記第１のトレーニングモデル群との関係を決定するように構成された主体解析モジュールと、
前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係、および前記主体モデルと前記第１のトレーニングモデル群との前記関係に基づいて予測形状を生成するように構成された予測モジュールと、を備え、
前記第１のトレーニングモデル群内の各モデルが前記第２のトレーニングモデル群内のモデルにそれぞれ対応しており、
前記第１のトレーニングモデル群は、二次元イメージまたは２次元モデルを含み、
前記第２のトレーニングモデル群は、三次元モデルを含むシステム。
【請求項２】
前記トレーニング解析モジュールは、
前記第１のトレーニングモデル群についての平均形状を決定し、
前記第１のトレーニングモデル群の分解行列を決定し、
前記第２のトレーニングモデル群についての平均形状を決定し、
前記第２のトレーニングモデル群の分解行列を決定するように構成されている、請求項１に記載のシステム。
【請求項３】
前記主体解析モジュールは、近似ベクトルを決定するように構成されている、請求項１または２に記載のシステム。
【請求項４】
前記予測モジュールは、主体モデルを修正することによって予測形状を生成するように構成されている、請求項１〜３の何れか１項に記載のシステム。
【請求項５】
前記予測モジュールは、近似ベクトルに少なくとも部分的に基づいて、予測形状を生成するように構成されている、請求項１〜４の何れか１項に記載のシステム。
【請求項６】
形状を予測するコンピュータ実行方法であって、
第１パラメータによって特徴付けられた第１のトレーニングモデル群と、第２パラメータによって特徴付けられた第２のトレーニングモデル群との関係を決定するステップと、
前記第１パラメータによって特徴付けられた主体モデルと、前記第１のトレーニングモデル群との関係を決定するステップと、
前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係、および前記主体モデルと前記第１のトレーニングモデル群との前記関係に基づき、前記第２パラメータによって特徴付けられる予測形状を生成するステップと、を含み、
前記第１のトレーニングモデル群内の各モデルが前記第２のトレーニングモデル群内のモデルにそれぞれ対応しており、
前記第１のトレーニングモデル群は、二次元イメージまたは２次元モデルを含み、
前記第２のトレーニングモデル群は、三次元モデルを含むコンピュータ実行方法。
【請求項７】
前記第１のトレーニングモデル群は解剖学上のモデルを含み、前記第１パラメータは骨である、請求項６に記載のコンピュータ実行方法。
【請求項８】
前記第２のトレーニングモデル群は解剖学上のモデルを含み、前記第２パラメータは軟骨である、請求項６または７に記載のコンピュータ実行方法。
【請求項９】
前記第１のトレーニングモデル群および前記第２のトレーニングモデル群はＭＲＩ画像に由来するものであり、前記主体モデルはＸ線画像に由来するものである、請求項６〜８の何れか１項に記載のコンピュータ実行方法。
【請求項１０】
前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係を決定するステップは、
前記第１のトレーニングモデル群についての平均形状を決定するステップと、
前記第１のトレーニングモデル群の分解行列を決定するステップと、
前記第２のトレーニングモデル群についての平均形状を決定するステップと、
前記第２のトレーニングモデル群の分解行列を決定するステップと、を含む請求項６〜９の何れか１項に記載のコンピュータ実行方法。
【請求項１１】
前記主体モデルと前記第１のトレーニングモデル群との関係を決定するステップは、近似ベクトルを決定するステップを含む、請求項６〜１０の何れか１項に記載のコンピュータ実行方法。
【請求項１２】
予測形状を生成するステップは、前記主体モデルを修正するステップを含む、請求項６〜１１の何れか１項に記載のコンピュータ実行方法。
【請求項１３】
予測形状を生成するステップは、前記近似ベクトルに少なくとも部分的に基づいて行われる、請求項１１に記載のコンピュータ実行方法。
【請求項１４】
予測形状のうち精度レベルの閾値を超える精度を有する少なくとも１つの部位を同定するステップをさらに含む、請求項６〜１３の何れか１項に記載のコンピュータ実行方法。
【請求項１５】
プロセッサによって実行されるときにコンピュータに形状予測方法を実行させるコンピュータ実行可能な命令を記憶した、コンピュータ読み取り可能な非一時的記録媒体であって、
前記形状予測方法は、
第１パラメータによって特徴付けられた第１のトレーニングモデル群と、第２パラメータによって特徴付けられた第２のトレーニングモデル群との関係を決定するステップと、
前記第１パラメータによって特徴付けられた主体モデルと、前記第１のトレーニングモデル群との関係を決定するステップと、
前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係、および前記主体モデルと前記第１のトレーニングモデル群との前記関係に基づき、前記第２パラメータによって特徴付けられる予測形状を生成するステップと、を含み、
前記第１のトレーニングモデル群内の各モデルが前記第２のトレーニングモデル群内のモデルにそれぞれ対応している、
前記第１のトレーニングモデル群は、二次元イメージまたは２次元モデルを含み、
前記第２のトレーニングモデル群は、三次元モデルを含む記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は、形状モデルから形状を予測することに関する。より具体的には、本願は、統計学的な形状モデルから、患者特定化形状または症例特定化形状を予測することに関する。
【背景技術】
【０００２】
統計学的なモデル形状は、時には画像の自動判読のために用いられる。所定種類の画像（画像の種類としては、例えば顔画像、手の画像などがある）内構造における形状上、および空間的関係上の「法則的な」変化パターンをトレーニング群から確立するのが、モデルについて用いられる基本的な考えである。また、法則的な変化パターンを効率的にパラメータ化するために、統計学的な解析が用いられ、形状に対する簡潔な表現が実現されている。統計学的な解析は、解析される画像内構造の形状が所望の対象種類のもっともらしい例であるかどうかを決定するために用いられる形状制約も提供している。
【０００３】
統計学的な形状モデルは、一般的に２Ｄまたは３Ｄ形状であるトレーニング形状の画像を分割することによって構築され得る。また、統計学的な形状モデルは、医療分野に利用することができる。この場合、モデルの発展に用いられる画像は、例えばＸ線、ＣＴスキャンまたはＭＲＩによる画像であってもよい。さらに、統計学的な形状モデルは、患者の人体組織の形状解析、または、患者の解剖学的構造のモデル作成、例えば患者特定化外科ガイドの設定など、様々な医療用途に利用することができる。
【発明の概要】
【０００４】
本発明の実施形態によれば、添付の請求項に記載されたように形状を予測するための、システム、方法、および、命令を運ぶコンピュータ読み取り可能な記録媒体を提供する。
【０００５】
より具体的には、本発明の第１態様では、形状を予測するシステムを提供する。該システムは、第１パラメータによって特徴付けられた第１のトレーニングモデル群と、第２パラメータによって特徴付けられた第２のトレーニングモデル群と、前記第１パラメータによって特徴付けられた主体モデルとを記憶するように構成された少なくとも１つのデータベースを含み、前記第１のトレーニングモデル群内の各モデルが前記第２のトレーニングモデル群内のモデルにそれぞれ対応している。また、該システムは、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との関係を決定するように構成されたトレーニング解析モジュールをさらに含む。また、該システムは、前記主体モデルと前記第１のトレーニングモデル群との関係を決定するように構成された主体解析モジュールと、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係、および前記主体モデルと前記第１のトレーニングモデル群との前記関係に基づいて予測形状を生成するように構成された予測モジュールと、をさらに含む。また、該システムにおいて、前記第１のトレーニングモデル群は、二次元イメージまたは２次元モデルを含み、前記第２のトレーニングモデル群は、三次元モデルを含む。
【０００６】
第１態様の一実施形態において、前記トレーニング解析モジュールは、前記第１のトレーニングモデル群についての平均形状を決定し、前記第１のトレーニングモデル群の分解行列を決定し、前記第２のトレーニングモデル群についての平均形状を決定し、前記第２のトレーニングモデル群の分解行列を決定するように構成されている。また、第１態様の別の実施形態において、前記主体解析モジュールは、近似ベクトルを決定するように構成されている。前記第１態様の別の実施形態において、前記予測モジュールは、主体モデルを修正することによって予測形状を生成するように構成されている。第１態様の更なる別の実施形態において、前記予測モジュールは、近似ベクトルに少なくとも部分的に基づいて、予測形状を生成するように構成されている。
【０００７】
本発明の第２態様では、形状を予測するコンピュータ実行方法を提供する。該方法は、第１パラメータによって特徴付けられた第１のトレーニングモデル群と、第２パラメータによって特徴付けられた第２のトレーニングモデル群との関係を決定するステップを含み、前記第１のトレーニングモデル群内の各モデルが前記第２のトレーニングモデル群内のモデルにそれぞれ対応している。また、該方法は、前記第１パラメータによって特徴付けられた主体モデルと、前記第１のトレーニングモデル群との関係を決定するステップをさらに含む。また、該方法は、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係、および前記主体モデルと前記第１のトレーニングモデル群との前記関係に基づき、前記第２パラメータによって特徴付けられる予測形状を生成するステップをさらに含む。また、該方法において、前記第１のトレーニングモデル群は、二次元イメージまたは２次元モデルを含み、前記第２のトレーニングモデル群は、三次元モデルを含む。
【０００８】
第２態様の別の実施形態において、前記第１のトレーニングモデル群は解剖学上のモデルを含み、前記第１パラメータは骨である。第２態様の別の実施形態において、前記第２のトレーニングモデル群は解剖学上のモデルを含み、前記第２パラメータは軟骨である。第２態様の更なる別の実施形態において、前記第１のトレーニングモデル群および前記第２のトレーニングモデル群がＭＲＩ画像に由来するものであり、前記主体モデルがＸ線画像に由来するものである。第２態様の更なる別の実施形態において、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との関係を決定するステップは、前記第１のトレーニングモデル群についての平均形状を決定するステップと、前記第１のトレーニングモデル群の分解行列を決定するステップと、前記第２のトレーニングモデル群についての平均形状を決定するステップと、前記第２のトレーニングモデル群の分解行列を決定するステップとを含む。第２態様の更なる別の実施形態において、前記主体モデルと前記第１のトレーニングモデル群との関係を決定するステップは、近似ベクトルを決定するステップを含む。第２態様の更なる別の実施形態において、予測形状を生成するステップは、前記主体モデルを修正するステップを含む。第２態様の更なる別の実施形態において、予測形状を生成するステップは、前記近似ベクトルに少なくとも部分的に基づいて行われる。第２態様の更なる別の実施形態において、該方法は、予測形状のうち精度レベルの閾値を超える精度を有する少なくとも１つの部位を同定するステップをさらに含んでもよい。
【０００９】
本発明の第３態様では、プロセッサによって実行されるときにコンピュータに形状予測方法を実行させるコンピュータ実行可能な命令を記憶した、コンピュータ読み取り可能な非一時的記録媒体を提供する。該方法は、第１パラメータによって特徴付けられた第１のトレーニングモデル群と、第２パラメータによって特徴付けられた第２のトレーニングモデル群との関係を決定するステップを含んでもよく、前記第１のトレーニングモデル群内の各モデルが、前記第２のトレーニングモデル群内のモデルにそれぞれ対応している。また、該方法は、前記第１パラメータによって特徴付けられた主体モデルと、前記第１のトレーニングモデル群との関係を決定するステップをさらに含んでもよい。また、該方法は、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係、および前記主体モデルと前記第１のトレーニングモデル群との前記関係に基づき、前記第２パラメータによって特徴付けられる予測形状を生成するステップをさらに含んでもよい。また、該システムにおいて、前記第１のトレーニングモデル群は、二次元イメージまたは２次元モデルを含み、前記第２のトレーニングモデル群は、三次元モデルを含む。
【図面の簡単な説明】
【００１０】
【図１】本明細書に記載の各実施形態に適するコンピューティング環境の一例のブロック図である。
【図２】１つ以上の実施形態に用いてもよいコンピューティングシステムの高レベルシステム図である。
【図３】１つ以上の実施形態に用いることが可能な第１および第２のトレーニング画像群の例を例示する図である。
【図４Ａ】１つ以上の実施形態に係る形状予測方法を例示するフローチャートである。
【図４Ｂ】１つ以上の実施形態に係る形状予測方法のサブステップを例示するフローチャートである。
【図５Ａ】一実施形態において、骨モデルの表面から突出するベクトルとして表現された対応する予測軟骨形状を有する該骨モデルの一例を示す図である。
【図５Ｂ】図５Ａのベクトルで示された予測の軟骨形状を示す図である。
【発明を実施するための形態】
【００１１】
本願に係る実施形態は、主体モデルから形状を予測するシステムおよび方法に関する。本発明の様々な実施形態において、２つのトレーニングモデル群間の関係の決定と、前記主体モデルと前記第１のトレーニングモデル群との関係の決定と、これらの関係に基づく予測形状の生成とによって、形状を予測することができる。前記第１のトレーニングモデル群および前記主体モデル群は、第１パラメータによって特徴付けられてもよく、前記第２のトレーニングモデル群および前記予測形状は、第２のおよび対応する（若しくは潜在的に対応する）パラメータによって特徴付けられてもよい。例示的な医療上の応用として、前記第１および第２パラメータは、それぞれ骨および軟骨であってもよく、前記主体モデルは、特定関節の骨輪郭のモデルであってもよい。本発明の実施形態によれば、同関節における軟骨輪郭の予測を容易にすることができる。したがって、医療分野において、本発明の実施形態は、特定患者の骨輪郭しか映っていないＸ線画像から該患者の軟骨輪郭（それ以外ではＭＲＩ画像からしか見えなかろう）の形状を予測することに用いることができる。このような状況における形状予測により、ＭＲＩの費用を節約できると共に、患者特定化の解剖学的モデルまたは患者特定化外科ガイドを発展させることが可能となる。また、医療上の別の応用として、第１および第２のトレーニングモデル群は、心拍の第１相における心臓モデル、および心拍の第２相における心臓モデルを含んでもよい。このような状況における形状予測により、患者の心臓の相が１つしか映っていない主体画像を用いて心拍パターンを予測することが可能となる。
【００１２】
本明細書に記載のシステムおよび方法は、１つ以上のコンピュータシステムによって実行されてもよい。図１には、様々な実施形態の実行に適したコンピュータシステム１００の一例が示されている。該コンピュータステム１００は、一般的に、本明細書に記載の１つ以上の実施形態に係る何らかのプロセスおよび命令を実行するように構成されたコンピュータハードウェアという形式を採用してもよい。該コンピュータハードウェアは、単独のコンピュータであってもよく、協働するように構成された複数のコンピュータであってもよい。コンピュータシステム１００はプロセッサ１０２を含む。該プロセッサは、一般的に、形状予測に関する何らかのタスクを完遂させるためのコンピュータ命令を実行するように構成されている。該プロセッサ１０２は、標準規格のパソコンのプロセッサ、例えばＩｎｔｅｌ社、Ａｄｖａｎｃｅｄ  Ｍｉｃｒｏ  Ｄｅｖｉｃｅｓ社、Ａｐｐｌｅ社、ＡＲＭ社やＭｏｔｏｒｏｌａ社によって設計および／または発売されたプロセッサ等であってもよい。また、該プロセッサ１０２は、画像の処理および解析に特化したプロセッサであってもよい。システム１００は、メモリ１０４を含んでもよい。該メモリ１０４は、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｍｏｒｙ）のような形式の揮発性メモリ１０４Ａを含んでもよい。揮発性メモリ１０４Ａは、実行可能なソフトウェア型モジュールをメモリ内にロードし、該ソフトウェア型モジュールが本分野の周知方式でプロセッサ１０２によって実行されるという構成となっていてもよい。また、該ソフトウェア型モジュールは、不揮発性メモリ１０４Ｂ内に記憶されていてもよい。該不揮発性メモリ１０４Ｂは、ハードディスクドライブ、フラッシュメモリ、固体ハードドライブ、または別タイプの不揮発性メモリという形式を採用してもよい。また、該不揮発性メモリ１０４Ｂは、実行不可能なデータ、例えばデータベースファイル等の記憶に用いてもよい。
【００１３】
また、コンピュータシステム１００は、ネットワークインターフェース１０６を含んでもよい。該ネットワークインターフェースは、システム１００をネットワーク（例えばインターネットなど）にアクセスさせるように構成されたネットワークインターフェースカードおよびそれに対応するソフトウェア型ドライバ、並びに／或いは、ファームウェア、という形式を採用してもよい。ネットワークインターフェースカード１０６は、異なるタイプの様々なネットワークにアクセスするように構成されてもよい。例えば、ネットワークインターフェースカード１０６が、公的にはアクセス不能な私的ネットワークにアクセスするように構成されていてもよい。また、ネットワークインターフェースカード１０６は、例えばＥＶＤＯ、ＷｉＭａｘやＬＴＥネットワークのような無線データ伝送技術等を用いた無線ネットワークにアクセスするように構成されていてもよい。なお、図２に１つのネットワークインターフェースカード１０６が示されているが、異なるタイプのネットワークにアクセスするための複数のネットワークインターフェースカード１０６が存在してもよい。さらに、ネットワークインターフェースカード１０６は、異なるタイプの複数のネットワークへの接続が可能となるように構成されていてもよい。
【００１４】
また、コンピュータシステム１００にオペレーティングシステム１０８が含まれていてもよい。該オペレーティングシステム１０８は、コンピュータソフトウェア形式のアプリケーションをプロセッサ１０２に実行させることのできるプラットフォームを備えるように設計された、例えばＬｉｎｕｘ（登録商標）、Ｗｉｎｄｏｗｓ（登録商標）、ＣｈｒｏｍｅＯＳまたはＭａｃＯＳのような一般周知のオペレーティングシステムであってもよい。あるいは、該オペレーティングシステム１０８は、画像の処理解析アプリケーションのために特別に設計された特定用途のオペレーティングシステムであってもよい。さらに、該オペレーティングシステム１０８は、Ａｎｄｒｏｉｄ（登録商標）、Ｗｉｎｄｏｗｓ  ＭｏｂｉｌｅやＩＯＳのようなモバイル用オペレーティングシステムであってもよい。
【００１５】
オペレーティングシステム１０８上においてウェブサーバ・ソフトウェア１３０が動作していてもよい。該ウェブサーバ・ソフトウェア１１０は、例えばアパッチ（Ａｐａｃｈｅ）、インターネット情報サーバ（Ｉｎｔｅｒｎｅｔ  Ｉｎｆｏｒｍａｔｉｏｎ  Ｓｅｒｖｅｒ）、または他のウェブサーバ・ソフトウェアなど、現在市販の標準ウェブサーバ製品であってもよい。あるいは、該ウェブサーバは、オペレーティングシステム１０８の一部であってもよく、インターネットや他のＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）またはＷＡＮ（Ｗｉｄｅ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）などのネットワークを介してソフトウェアを閲覧するための画像関連または解析関連ウェブページを配信するように特別に構成された専門のＨＴＴＰサーバであってもよい。該ウェブサーバ・ソフトウェア１３０は、オペレーティングシステム１０８に備えられているオペレーティン・プラットフォーム上においてプロセッサ１０２によってアクセスされて実行されるように、メモリ１０４内に記憶されていてもよい。
【００１６】
コンピュータシステム１００は、アプリケーションサーバ１１２をさらに含んでもよい。アプリケーションサーバ１１２は、コンピュータシステム１００上に元々動作可能なアプリケーション、またはウェブサーバ１１０を介してアクセスされるようなアプリケーション、またはこの両方に該当するネットワークベース型アプリケーションを提供するように構成されたコンピュータハードウェアおよび／またはソフトウェアを含んでもよい。また、該アプリケーションサーバは、後に詳述する形状予測システムの配布、利用およびメンテナンスが可能となるように構成されてもよい。
【００１７】
図２は、一実施形態に係る、形状を予測するためのコンピュータベース型システム２００を例示している。該システム２００は、図１に示された前述のコンピュータ１００のような１つ以上のコンピュータから成り得る。また、システム２００はデータベース２０２を含むことができる。該データベース２０２は、システム２００内の様々な機能、例えばキャプチャ、記憶、および／または、ソフトウェアアプリケーションや他のデバイスからデータベース内の画像データ若しくはモデルデータへのアクセス提供が実行されるように構成されてもよい。データベース２０２は、関係データベース、オブジェクト指向データベース、オブジェクト関係データベース、または意図された機能の実行に適する他のいかなるタイプのデータベースであってもよい。また、データベース２０２は、未処理の２Ｄおよび／または３Ｄ画像データ、例えば、写真、ＣＴスキャン、ＭＲＩ、Ｘ線などに由来する医学的画像データを記憶するように構成されてもよい。代わりにまたは加えて、データベース２０２は、２Ｄおよび／または３Ｄモデル、例えば、骨の形状、軟骨の形状および／または軟組織の形状に関するトレーニングモデルおよび／もしくは患者特定化モデルを記憶するように構成されてもよい。また、データベース２０２は、データベース内に記憶されている画像またはモデルデータに関係する出所特定化または形状特定化データ（匿名データか、保護されたデータ）を含んでもよい。さらに、データベース２０２は、画像またはモデルデータに関する形状特定化データ、例えば、画像またはモデルに関連付けられた識別点および／もしくは湾曲部に関する情報を含んでもよい。該形状特定化データは、例えば変換技術者によって手動で示されてもよく、例えばコンピュータ補助設計（ＣＡＤ）ツールや前処理ツールによって自動で示されてもよく、原来の形状における切断箇所または切断部位（例えば骨の交点および骨の切断面）を単純に示してもよい。一部の実施形態において、画像またはモデルのデータは、トレーニング画像群またはトレーニングモデル群を１つ以上含み、かつ、形状の予測に用いられる主体画像または主体モデル（例えば患者特定化モデル）を１つ以上含んでもよい。トレーニング形状群は、画像またはモデルを幾つでも含んでよい。例えば、トレーニング群には、約１０個、５０個、１００個、５００個または１０００個の形状が含まれてもよく、これらの数よりも少ない数若しくは多い数、またはこれらの数のうち何れかで定義される範囲内の数の形状が含まれてもよい。
【００１８】
システム２００は、図２に示されるように、データベース２０２内の２Ｄおよび／または３Ｄのトレーニング画像若しくはトレーニングモデル（任意に加えて、関連付けられた形状特定化データ）を解析すると共に、その関連の様々な機能を実行するように構成されたトレーニング解析モジュール２０４をさらに含んでもよい。トレーニング解析モジュール２０４は、例えば、ＣＡＤ、前処理、パラメータ化、網目化（ｍｅｓｈｉｎｇ）（例えば三角網）および／または他のモデル化機能を実行するように構成されてもよい。また、トレーニング解析モジュール２０４は、例えば主成分解析機能を含む様々な解析機能を実行するように構成されてもよい。また、トレーニング解析モジュール２０４は、自オペレーティングシステム内にインストールされた２Ｄおよび／または３Ｄ画像処理ソフトウェアおよび／若しくはモデル化ソフトウェアを有するコンピューティング装置、という形式を採用することができる。トレーニング解析モジュール２０２は、主にまたは専らソフトウェアから構成されてもよく、ハードウェアとソフトウェアとの組み合わせから構成されてもよく、更に別の実施形態として、ＡＳＩＣや他のタイプのマイクロプロセッサのような専用ハードウェアであってもよい。一部の実施形態では、何らかのモデル化機能性および／または解析機能性を一つのソフトウェアアプリケーションによって提供すると共に、他のモデル化機能性および／または解析機能性を１つ以上の別個のコンピュータアプリケーションによって提供してもよい。また、代替的な構成として、これらすべての機能性を単一のコンピュータプログラムによって提供してもよい。
【００１９】
また、システム２００は、図２に示されるように、データベース２０２内の２Ｄおよび／または３Ｄの主体画像若しくは主体モデル（任意に加え、関連付けられた形状特定化データ）を解析すると共に、その関連の様々な機能を実行するように構成された主体解析モジュール２０６をさらに含んでもよい。主体解析モジュール２０６は、例えば、ＣＡＤ、前処理、パラメータ化、網目化（例えば三角網）および／または他のモデル化機能を実行するように構成されてもよい。また、主体解析モジュール２０６は、例えば主成分解析機能を含む様々な解析機能を実行するように構成されてもよい。また、主体解析モジュール２０６は、自オペレーティングシステム内にインストールされた２Ｄおよび／または３Ｄ画像処理ソフトウェアおよび／若しくはモデル化ソフトウェアを有するコンピューティング装置、という形式を採用することができる。主体解析モジュール２０６は、主にまたは専らソフトウェアから構成されてもよく、ハードウェアとソフトウェアとの組み合わせから構成されてもよい。一部の実施形態では、何らかのモデル化機能性および／または解析機能性を一つのソフトウェアアプリケーションによって提供すると共に、他のモデル化機能性および／または解析機能性を１つ以上の別個のコンピュータアプリケーションによって提供してもよい。また、代替的な構成として、これらすべての機能性を単一のコンピュータプログラムによって提供してもよい。なお、図２には別個のモジュールとして例示されているが、一部の実施形態では、トレーニング解析モジュール２０４および主体解析モジュール２０６は、共に単一のコンピュータプログラムによって実現されることができる。
【００２０】
さらに、システム２００は、図２に示されるように予測モジュール２０８をさらに含んでもよい。予測モジュール２０８は、主体画像または主体モデルに関する情報（該情報は主体解析モジュール２０６から受信されてもよい）、トレーニング画像またはトレーニングモデルに関する情報（該情報はトレーニング解析モジュール２０４から受信されてもよい）、および、該主体モデルと該トレーニングモデルとの関係に関する情報（該情報はトレーニング解析モジュール２０４および／または主体解析モジュール２０６から受信されもよい）を基に、形状を予測するように構成されてもよい。主体解析モジュール２０６は、例えば、ＣＡＤ、前処理、パラメータ化、網目化（例えば三角網）および／または他のモデル化機能を実行するように構成されてもよい。予測モジュール２０８は、自オペレーティングシステム内にインストールされた２Ｄおよび／または３Ｄ画像の処理ソフトウェアおよび／若しくはモデル化ソフトウェアを有するコンピューティング装置、という形式を採用することができる。また、予測モジュール２０８は、主にまたは専らソフトウェアから構成されてもよく、ハードウェアとソフトウェアとの組み合わせから構成されてもよい。一部の実施形態では、何らかのモデル化機能性および／または解析機能性を一つのソフトウェアアプリケーションによって提供すると共に、他のモデル化機能性および／または解析機能性を１つ以上の別個のコンピュータアプリケーションによって提供してもよい。また、代替的な構成として、これらすべての機能性を単一のコンピュータプログラムによって提供してもよい。なお、図２には別個のモジュールとして例示されているが、一部の実施形態では、トレーニング解析モジュール２０４、主体解析モジュール２０６および予測モジュール２０８は、共に単一のコンピュータプログラムによって実現されることができる。
【００２１】
一部の実施形態において、トレーニング解析モジュール２０４、主体解析モジュール２０６、および／または予測モジュール２０８は、数値コンピューティングシステム、例えば米マサチューセッツ州のネイティック市のＭａｔｈＷｏｒｋｓ社から入手できるＭＡＴＬＡＢ（登録商標）環境を備えていてもよい。Ｃ＋＋、Ｐｙｔｈｏｎ、Ｍａｐｌｅ、ＦｒｅｅＭａｔ、ＩＧＯＲ  Ｐｒｏ、ｊＢＥＡＭ等を含む他の数値解析ソフトウェアパッケージを活用してもよいが、これらに限定されない。
【００２２】
図２を参照して説明したシステム２００は、単に本明細書に開示のシステムおよび方法を実行する数多くの好適な環境の１つであるに過ぎないことが理解されるべきである。システム２００は、コンピュータネットワークと接続されていてもよく、いかなるコンピュータネットワークからも隔離され独立していてもよいことが、当業者に容易に理解される。また、本明細書では、図２に記載の様々なコンポーネントを別個のものとして説明したが、これらのコンポーネントは、機能的に統合された少数のモジュールであってもよく、あるいは、分離されている追加モジュールであってもよいことが当業者に理解される。
【００２３】
図３は、１つ以上の実施形態に用いることが可能な第１および第２のトレーニング画像群の例を説明する図である。第１のトレーニング画像群におけるトレーニング画像３０２（ａ）、３０４（ａ）、３０６（ａ）は、３名の異なる患者の膝関節近傍の大腿骨および脛骨の形状を表した骨モデルであり、軟骨無しに骨のみによって特徴付らけれている。第２のトレーニング画像群におけるトレーニング画像３０２（ｂ）、３０４（ｂ）、３０６（ｂ）は、それぞれ、同３名の患者の大腿骨および脛骨の形状を表した軟骨モデルであり、軟骨によって特徴付けられている。図３のように、大腿骨と脛骨との間の隙間近傍の膝関節について、軟骨モデルのほうが、対応する骨モデルより僅かに大きい傾向がある。画像３０２（ａ）、３０２（ｂ）、３０４（ａ）、３０４（ｂ）、３０６（ａ）３０６（ｂ）は、例えば市販で入手できるＣＡＤ（コンピュータ補助設計）および／または他の画像処理ソフトウェアを用い、患者の例えば医療スキャン（例えばＭＲＩ）から構築することができる。
【００２４】
本発明の実施形態において、各々のトレーニング画像のモデルを形成するために第１および第２のトレーニング画像群を分割してもよい。この処理には、例えば、各々の画像のモデルを形成するために各画像群内の画像について分割する工程と、第１画像群および第２画像群からの対であるモデル間の対応を探す工程と、（例えばデカルト座標、極座標、球座標を用いて）各モデルをベクトルで定義する工程とが含まれ得る。そして、トレーニングモデルについてのベクトルは、行列として結合することができる。なお、主体画像または画像は、１つ以上のベクトルで定義してもよい。
【００２５】
図４Ａを参照すると、１つ以上の実施形態に係る、形状を予測するためのプロセス４００が例示されている。プロセス４００は、第１のモデル群と第２のモデル群との関係を決定するブロック４０２から開始する。また、プロセス４００において、主体モデルと第１のモデル群との関係の決定を含むブロック４０４に移行することができる。一部の実施形態では、ブロック４０４において実行される処理（すなわち、主体モデルと第１のモデル群との関係の決定）は、例えば近似アルゴリズムを用いた、近似ベクトルの決定を含んでもよい。近似ベクトルは、例えば、第１のトレーニング群について主体モデルを定義するための適切な係数ベクトルであってもよい。そして、プロセス４００において、第１のモデル群と第２のモデル群との関係、および主体モデルと第１のモデル群との関係に基づいて予測形状を生成するブロック４０６に移行することができる。
【００２６】
本発明の実施形態において、プロセス４００における１つ以上のステップおよび／またはサブステップを、主成分解析法を用いて実行してもよい。例えば、図４Ｂにさらに例示されるように、ブロック４０２において実行される処理（すなわち、第１のモデル群と第２のモデル群との関係の決定）は、第１のトレーニングモデル群についての平均分割（ｓｅｇｍｅｎｔａｔｉｏｎ）または平均形状を決定するブロック４１０を含んでもよい。また、プロセス４０２は、第１のトレーニングモデル群についての分解行列を決定するブロック４１２を含んでもよい。そして、プロセス４０２は、第２のトレーニングモデル群についての平均分割または平均形状を決定するブロック４１４を含んでもよい。また、プロセス４０２は、第２のトレーニングモデル群についての分解行列を決定するブロック４１６を含んでもよい。一部の実施形態では、主成分解析の結果（例えば、ステップ４１０、４１２、４１４および４１６の結果）を利用してブロック４０６（すなわち予測形状の生成）を実行してもよい。本発明の実施形態において、以上に言及された関係（例えば、第１のモデル群と第２のモデル群との関係）を一度決定しておき、それを異なる主体モデルと共に繰返して用いることにより、各々の主体モデルの予測形状を生成してもよい。
【００２７】
以下は、一部の実施形態において予測形状の生成に利用可能な方程式の例である。
【００２８】
方程式１：ＳＳ＝Ａ＊ｃ
方程式２：ＳＳ＝μ＋Ｕ１＊ｅ
方程式３：ＳＰ＝Ｂ＊ｄ
方程式４：ＳＰ＝η＋Ｕ２＊ｅ
式中、
ＳＳは、主体モデルを定義するベクトルである。
【００２９】
Ａは、第１のトレーニングモデル群を定義する行列である。
【００３０】
ｃは、主体モデルと第１のトレーニングモデル群との関係を表す近似ベクトル（例えば係数ベクトル）である。
【００３１】
μは、第１のトレーニングモデル群の加重平均分割（例えば、第１群についての全ベクトルの合計をこれらベクトルの数で割って得られた、均一な加重平均分割）である。
【００３２】
Ｕ１は、第１のトレーニングモデル群についての分解行列である。
【００３３】
ｅは、主体モデルと第１のトレーニングモデル群との関係を表す、代替の近似ベクトル（例えば係数ベクトル）である。
【００３４】
ＳＰは、予測モデルまたは予測形状を定義するベクトルである。
【００３５】
Ｂは、第２のトレーニングモデル群を定義する行列である。
【００３６】
ｄは、予測形状モデルと第２のトレーニングモデル群との関係を表す近似ベクトル（例えば係数ベクトル）である。
【００３７】
ηは、第２のトレーニングモデル群の加重平均分割（例えば、第２群についての全ベクトルの合計をこれらベクトルの数量で割って得られた、均一な加重平均分割）である。
【００３８】
Ｕ２は、第２のトレーニングモデル群についての分解行列である。
【００３９】
一実施形態において、近似アルゴリズムを用い、主体モデルＳＭに関する情報、および第１のトレーニングモデル群に関する情報に基づいて近似ベクトルｅを得る（例えば方程式２を参考）ことができる。そして、第２のトレーニングモデル群についての平均形状ηについて予測形状ＳＰの定義を可能にする関係Ｕ２を、主体モデルと第１のトレーニングモデル群との関係ｅを用いて決定する（例えば方程式４を参考）ことができる。
【００４０】
また、一実施形態において、以下のようにｃを決定することができる。
【００４１】
方程式５：Ａ＊ｃ＝μ＋Ｕ１＊ｅ
方程式６：ｃ＝Ａ−１＊（μ＋Ｕ１＊ｅ）
式中、Ａ−１は、行列Ａの仮の逆行列である。予測形状ＳＰは、単純に方程式３中のベクトルｄを方程式６によるベクトルｃと置き換えた下記方程式を用いて、決定することができる。
【００４２】
方程式７：ＳＰ＝Ｂ＊ｃ＝Ｂ＊Ａ−１＊（μ＋Ｕ１＊ｅ）＝Ｂ＊Ａ−１＊μ＋Ｂ＊Ｕ１＊ｅ
それから、ηおよびＵ２を下記のように定義することができる。
【００４３】
方程式８：η＝Ｂ＊Ａ−１＊μ
方程式９：Ｕ２＝Ｂ＊Ｕ１
そして、予測形状ＳＰは下記のようになる。
【００４４】
方程式１０：ＳＰ＝η＋Ｕ２＊ｅ
この実施形態では、ηは、第２のトレーニングモデル群の算術平均の代わりに、要素Ａ−１＊μで決定された重みでの第２のトレーニングモデル群の加重平均であってもよい。
【００４５】
また、別の実施形態において、ηは第２のトレーニングモデル群の実数算術平均であってもよい。この実施形態では、例えば下記のように行列Ａ内の各列からベクトルμを減算する特異値分解を用い、行列Ａとして行列[Ａ−μ]を定義することができる。
【００４６】
方程式１１：［Ａ−μ］＝Ｕ１＊Σ１＊Ｖ１Ｔ
上記式は、トレーニングモデル群についての分解行列Ｕ１の一演算方法を実現する方程式１および方程式２と整合するものである。また、方程式２は、下記のように書き直すことができる。
【００４７】
方程式１２：ＳＳ＝μ＋［Ａ−μ］＊ｆ
式中のｆは、下記を満たす。
【００４８】
方程式１３：ｆ＝Ｖ１＊Σ１−１＊ｅ
また、方程式２を下記何れの関係のように書き直すこともできる。
【００４９】
方程式１４：ＳＳ＝μ＋［Ａ−μ］＊ｆ
方程式１５：ＳＳ＝μ＋（Ｕ１＊Σ１＊Ｖ１Ｔ）＊ｆ
方程式１６：ＳＳ＝μ＋（Ｕ１＊Σ１＊Ｖ１Ｔ）＊（Ｖ１＊Σ１−１＊ｅ）
そして、下記のようにベクトルｆを用いてＳＰを予測することができる。
【００５０】
方程式１７：ＳＰ＝η＋［Ｂ−η］＊ｆ
したがって、Ｕ２を下記のように定義すると、
方程式１８：Ｕ２＝［Ｂ−η］＊Ｖ１＊Σ１−１
方程式４を下記何れの関係のように書き直すこともできる。
【００５１】
方程式１９：ＳＰ＝η＋［Ｂ−η］＊ｆ
方程式２０：ＳＰ＝η＋［Ｂ−η］＊Ｖ１＊Σ１−１＊ｅ
図５Ａは、一実施形態において、骨モデルの表面から突出するベクトルとして表現された対応する予測軟骨形状を有する主体モデルの一例を示す図である。骨モデルは、画像処理中の膝およびＸ線画像の３Ｄポジショニングを容易にするためにＸ線撮像時に用いられた目印と共に、Ｘ線画像から生成することができる。一部の実施形態では、骨モデルを生成するために、Ｘ線に基づく輪郭（例えば手動きで）を統計学的な形状モデルに近似させることができる。図５Ｂは、図５Ａのベクトルで示された予測の軟骨形状を示す図である。
【００５２】
一部の実施形態では、予測形状を用いて外科的計画（例えば、穿孔（ｄｒｉｌｌｉｎｇ）または組織切除（ｄｅｂｒｉｄｅｍｅｎｔ）のための領域を計画する）を行うことができる。予測された軟骨形状を用いて、１つ以上の患者特定化外科ガイド（例えば、外科的処置中に患者の人体組織上へ設置することと、外科的手順のガイドに用いることとができる構造）を設定することもできる。また、一部の実施形態では、予測形状について、高精度および低精度の部位（例えば、典型的な高信頼性の精度の部位、および、ガイド設定への使用を回避すべき部位または他の信頼できない部位）を示す精度マップを生成することができる。
【００５３】
〔まとめ〕
より具体的には、本発明の第１態様では、形状を予測するシステムを提供する。該システムは、第１パラメータによって特徴付けられた第１のトレーニングモデル群と、第２パラメータによって特徴付けられた第２のトレーニングモデル群と、前記第１パラメータによって特徴付けられた主体モデルとを記憶するように構成された少なくとも１つのデータベースを含み、前記第１のトレーニングモデル群内の各モデルが前記第２のトレーニングモデル群内のモデルにそれぞれ対応している。また、該システムは、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との関係を決定するように構成されたトレーニング解析モジュールをさらに含む。また、該システムは、前記主体モデルと前記第１のトレーニングモデル群との関係を決定するように構成された主体解析モジュールと、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係、および前記主体モデルと前記第１のトレーニングモデル群との前記関係に基づいて予測形状を生成するように構成された予測モジュールと、をさらに含む。第１態様の一実施形態において、前記トレーニング解析モジュールは、前記第１のトレーニングモデル群についての平均形状を決定し、前記第１のトレーニングモデル群の分解行列を決定し、前記第２のトレーニングモデル群についての平均形状を決定し、前記第２のトレーニングモデル群の分解行列を決定するように構成されている。また、第１態様の別の実施形態において、前記主体解析モジュールは、近似ベクトルを決定するように構成されている。前記第１態様の別の実施形態において、前記予測モジュールは、主体モデルを修正することによって予測形状を生成するように構成されている。第１態様の更なる別の実施形態において、前記予測モジュールは、近似ベクトルに少なくとも部分的に基づいて、予測形状を生成するように構成されている。
【００５４】
本発明の第２態様では、形状を予測するコンピュータ実行方法を提供する。該方法は、第１パラメータによって特徴付けられた第１のトレーニングモデル群と、第２パラメータによって特徴付けられた第２のトレーニングモデル群との関係を決定するステップを含み、前記第１のトレーニングモデル群内の各モデルが前記第２のトレーニングモデル群内のモデルにそれぞれ対応している。また、該方法は、前記第１パラメータによって特徴付けられた主体モデルと、前記第１のトレーニングモデル群との関係を決定するステップをさらに含む。また、該方法は、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係、および前記主体モデルと前記第１のトレーニングモデル群との前記関係に基づき、前記第２パラメータによって特徴付けられる予測形状を生成するステップをさらに含む。第２態様の別の実施形態において、前記第１のトレーニングモデル群は解剖学上のモデルを含み、前記第１パラメータは骨である。第２態様の別の実施形態において、前記第２のトレーニングモデル群は解剖学上のモデルを含み、前記第２パラメータは軟骨である。第２態様の更なる別の実施形態において、前記第１のトレーニングモデル群および前記第２のトレーニングモデル群がＭＲＩ画像に由来するものであり、前記主体モデルがＸ線画像に由来するものである。第２態様の更なる別の実施形態において、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との関係を決定するステップは、前記第１のトレーニングモデル群についての平均形状を決定するステップと、前記第１のトレーニングモデル群の分解行列を決定するステップと、前記第２のトレーニングモデル群についての平均形状を決定するステップと、前記第２のトレーニングモデル群の分解行列を決定するステップとを含む。第２態様の更なる別の実施形態において、前記主体モデルと前記第１のトレーニングモデル群との関係を決定するステップは、近似ベクトルを決定するステップを含む。第２態様の更なる別の実施形態において、予測形状を生成するステップは、前記主体モデルを修正するステップを含む。第２態様の更なる別の実施形態において、予測形状を生成するステップは、前記近似ベクトルに少なくとも部分的に基づいて行われる。第２態様の更なる別の実施形態において、該方法は、予測形状のうち精度レベルの閾値を超える精度を有する少なくとも１つの部位を同定するステップをさらに含んでもよい。
【００５５】
本発明の第３態様では、プロセッサによって実行されるときにコンピュータに形状予測方法を実行させるコンピュータ実行可能な命令を記憶した、コンピュータ読み取り可能な非一時的記録媒体を提供する。該方法は、第１パラメータによって特徴付けられた第１のトレーニングモデル群と、第２パラメータによって特徴付けられた第２のトレーニングモデル群との関係を決定するステップを含んでもよく、前記第１のトレーニングモデル群内の各モデルが、前記第２のトレーニングモデル群内のモデルにそれぞれ対応している。また、該方法は、前記第１パラメータによって特徴付けられた主体モデルと、前記第１のトレーニングモデル群との関係を決定するステップをさらに含んでもよい。また、該方法は、前記第１のトレーニングモデル群と前記第２のトレーニングモデル群との前記関係、および前記主体モデルと前記第１のトレーニングモデル群との前記関係に基づき、前記第２パラメータによって特徴付けられる予測形状を生成するステップをさらに含んでもよい。
【００５６】
何れか１つの実施形態について説明された特徴は、単独に用いられても、上記説明された他の特徴と組み合わせて用いられてもよく、他のいかなる実施形態の１つ以上の特徴と組み合わせて用いられても、他のいかなる実施形態のいかなる組合せと組み合わせて用いられてもよい、と理解されるべきである。さらに、以上に説明しなかった均等的構成および変更は、本発明の範囲から逸脱しない範囲で応用してもよい。本発明の範囲は、添付の請求項によって定義される。
【図１】
2019069158.tif 000003
【図２】
2019069158.tif 000004
【図３】
2019069158.tif 000005
【図４Ａ】
2019069158.tif 000006
【図４Ｂ】
2019069158.tif 000007
【図５Ａ】
2019069158.tif 000008
【図５Ｂ】
2019069158.tif 000009