(57)【特許請求の範囲】
【請求項１】
医用画像診断装置にて得た被検体の目的臓器が撮像された画像を読み込む手段と、
前記画像中に含まれる前記被検体の血管が撮影された領域を抽出するための条件と、前記画像中に含まれる前記被検体の腫瘍が撮影された領域を抽出するための条件と、の設定を受け付ける手段と、
前記画像中に含まれる前記被検体の血管が撮影された領域を抽出するための条件に従って前記画像中の第一領域を抽出するとともに、前記画像中に含まれる前記被検体の腫瘍が撮影された領域を抽出するための条件に従って、前記画像中の第二領域を抽出する手段と、
前記抽出された第一領域に連なる領域を支配領域として特定する支配領域特定手段と、
前記第一領域と第二領域とを表示する表示手段と、
前記表示された第一領域上の任意の位置の指定を受け付ける指定手段と、
前記任意の位置から末梢部までの支配領域を、前記目的臓器の他の領域と識別可能に表示する手段と、
を備えたことを特徴とする臓器の特定領域抽出表示装置。
【請求項２】
前記支配領域特定手段は、前記抽出された第一領域に相当する血管からの距離、又は前記抽出された第一領域に相当する血管の血管径に基づいて前記支配領域を特定する、
ことを特徴とする請求項１に記載の臓器の特定領域抽出表示装置。
【請求項３】
医用画像診断装置にて得た被検体の目的臓器が撮像された画像を用いて臓器の特定領域を抽出し表示する臓器の特定領域抽出表示装置における表示方法であって、
前記特定領域抽出表示装置が、前記被検体の画像を読み込むステップと、
前記特定領域抽出表示装置が、前記画像中に含まれる前記被検体の血管が撮影された領域を抽出するための条件と、前記画像中に含まれる前記被検体の腫瘍が撮影された領域を抽出するための条件と、の設定を受け付けるステップと、
前記特定領域抽出表示装置が、前記被検体の血管が撮影された領域を抽出するための条件に従って、前記画像中から第一領域を抽出するとともに、前記画像中に含まれる前記被検体の腫瘍が撮影された領域を抽出するための条件に従って、前記画像中から第二領域を抽出するステップと、
前記特定領域抽出表示装置が、前記第一領域に連なる領域を支配領域として特定するステップと、
前記特定領域抽出表示装置が、前記第一領域と前記第二領域とを表示するステップと、
前記特定領域抽出表示装置が、前記表示された第一領域上の任意の位置の指定を受け付けるステップと、
前記特定領域抽出表示装置が、前記任意の位置から末梢部までの支配領域を、前記目的臓器の他の領域と識別可能に表示するステップと、
を含むことを特徴とする臓器の特定領域抽出表示装置の表示方法。
【発明の詳細な説明】
【０００１】
【発明が属する技術分野】
本発明は、生体内の臓器の特定領域抽出表示方法及び装置に関するもので、特に、医師が臨床において肝臓のような臓器の診断や治療を行う際に医師を支援するためのシミュレーション方法及び装置に関するものである。
【０００２】
【従来の技術】
近年、Ｘ線撮影装置やＸ線ＣＴ装置やＭＲＩ装置などの医用画像診断装置で得られる画像を診断のみならず治療に用いることが盛んに行われるようになっている。治療には被検体にカテーテルを挿入して患部を切除するカテーテル術と、従来どおり切開手術により患部を切除する外科手術とがある。このうち、外科手術では、手術前に前記医用画像診断装置により患部の画像を得て、その切除する部分を決めておくことが通常行われる。
【０００３】
上記切除する部分を決めるための表示画像は三次元画像を用いている。これは人体の形態により近いので直感的に切除部分が決められるからである。一方、従来の画像上での臓器の抽出に関しては、目的とする臓器を他の臓器と分離する方法が研究されており、単一の臓器内の特定領域を抽出する方法としては、例えば臓器の三次元画像中に領域特定用の幾何学的平面又は曲面を医師等が解剖学的知識に基づいて設定して特定領域を設定するという方法が行われていた。
【０００４】
【発明が解決しようとする課題】
上記のような単一臓器内の特定領域を抽出する方法はあるが、実際の臨床の場において、シミュレーションのように臓器を幾何学的な平面や曲面で切除することは行われておらず、そのようなシミュレーションでは実際の臨床の場ではあまり役立つものとは言えないものであった。
【０００５】
また、肝臓のような臓器の切除領域は、門脈に略一致して流れる動脈だけでなく、静脈など様々な血管の走行や腫瘍の位置を考慮して設定されなければならないが、上記シミュレーションでは配慮されていなかった。
【０００６】
また、手術において、上記肝臓のような臓器の腫瘍を栄養する血管のみを最小限に切除して、臓器の機能を温存したいという要望があった。
【０００７】
本発明の目的は、複数種類の血管等が複雑に入り込んだ臓器の切除非抽出領域と識別可能に表示することができる方法と装置を提供することにある。
【０００８】
【課題を解決するための手段】
上記目的は、医用画像診断装置にて得た被検体の画像を用いて臓器の特定領域を抽出し表示装置に表示する方法において、前記画像中の目的臓器を構成する複数の組織の情報からそれぞれの組織に対応する条件を設定するステップと、前記条件を満たす組織情報に連なる領域を切除候補領域として抽出するステップと、該抽出された切除候補領域を前記目的臓器の非抽出領域と識別可能に表示するステップとを含むことを特徴とする臓器の特定領域抽出表示方法によって達成される。
ここでいう複数の組織とは、例えば目的臓器を肝臓とした場合、門脈、静脈、腫瘍等を選択的に設定できるものである。
【０００９】
また、前記表示ステップは、血管もしくは血管の中心線を、まわりの血管もしくはまわりの血管の中心線と異なる色若しくはは点滅表示などの異なる形態で表示することで達成される。
【００１０】
また、前記表示ステップは、前記医用画像診断装置より得られた画像を３次元的に積み上げ、積み上げ３次元画像を作成し、しきい値や各画素値に任意の不透明度を与えて三次元的に可視化を行う画像再構成方法により作成される画像についてその臓器の抽出領域を識別可能に表示することで達成される。
【００１１】
また、医用画像診断装置にて得た被検体の画像を用いて臓器の特定領域を抽出しその特定領域を表示する装置において、前記画像中の目的臓器を構成する複数の組織の情報からそれぞれの組織に対応する条件を設定する手段と、前記条件を満たす組織情報に連なる領域を切除候補領域として抽出する手段と、該抽出された切除候補領域を前記目的臓器の非抽出領域と識別可能に表示する手段とを備えたことを特徴とする臓器の特定領域抽出表示装置によって達成される。
【００１２】
また、前記抽出領域と非抽出領域との表示は、非抽出領域を表示しない表示モードか、抽出領域を非抽出領域と識別して表示するモードを設定し、その設定により表示モードを切り替えて表示することによって、非抽出領域と表示したとき際には全体との位置関係が、抽出領域のみを表示したときには抽出領域のみを際立たせて表示することができるので、診断能を向上させることができる。
【００１３】
【発明の実施の形態】
以下、添付図面に従って本発明に係る臓器の特定領域抽出表示方法及び装置の好ましい実施の形態について説明する。
【００１４】
実施の形態として、肝臓の造影撮影を行ったＸ線ＣＴ画像を処理対象画像として用い門脈の走行情報を利用して肝臓の領域特定を行う手順及びそれを用いて肝臓の切除部分を表示する方法を以下に説明する。
【００１５】
図１（ａ）に示すように、Ｘ線ＣＴ装置やＭＲＩ装置等の三次元計測の可能な画像診断装置で取得した複数の断層像１１を積み上げて図１（ｂ）に示すような積み上げ三次元画像１２とし、処理対象を三次元化する。積み上げ三次元画像１２は肝臓の組織と門脈を含み、ここには図示しない二次元の投影面に陰影付けして投影処理された擬似三次元画像として例えばモニタへ表示される。
【００１６】
図２は本発明の臓器の特定領域抽出表示方法を示すフローチャートである。本発明の領域特定処理８０は、関心領域抽出処理８１、距離値変換処理８２、細線化処理８３又は表面画素検出処理８４、支配領域特定処理８５、指定抽出処理８６及び切除処理８７から成る。
【００１７】
読み込まれた積み上げ三次元画像データに対し関心領域抽出処理８１が行われ、図３に示すように関心領域（対象臓器：肝臓実質３１、門脈３２）が抽出される。ここでは、簡単に説明するために、門脈３２だけを抽出したが、腫瘍や静脈も同様に抽出可能である。図４に示すように、抽出された門脈３２に対して距離値変換処理８２が行われ、更に細線化処理８３又は表面画素検出処理８４が行われる。距離値変換処理８２と細線化処理８３又は表面画素検出処理８４との処理は、矢印が示すように細線化処理を先に行いそれに続いて距離値変換処理を、また表面画素検出処理を先に行いそれに続いて距離値変換処理を行っても良い。処理結果を利用して、肝臓実質３１のうち門脈３２が支配する支配領域特定処理８５を行う。この支配領域特定処理８５においても、腫瘍や静脈も同様に領域特定可能である。つまり、支配領域特定処理８５を行う際に、門脈だけでなく、静脈や腫瘍についても行い、その腫瘍の支配領域に対応する門脈を求めれば、腫瘍に栄養する門脈を選択して切除することができるので、腫瘍の位置を考慮したシミュレーションとなり、腫瘍を栄養する血管のみを最小限に切除して、臓器の機能を温存することに配慮できるのである。
【００１８】
次に、門脈３２のうちの切除領域を特定する門脈枝を特定し抽出する関心領域部特定処理８６を行い、前記指定抽出処理８６により定義される切除領域を特定する切除処理８７を行う。
【００１９】
図５は本発明の臓器の特定領域抽出表示方法の実施に即した詳細な手順の一例を示すフローチャートである。まず積み上げ三次元画像データを読み込み（ステップ２１）、図３に示すように読み込まれた三次元画像データから関心領域（門脈３２、肝臓実質３１）の抽出処理を行う（ステップ２２）。この抽出処理には、積み上げ三次元画像データに対して、画像処理の分野においては公知の閾値による二値化を利用したセグメンテーションや領域拡張法を利用する。
【００２０】
抽出された門脈３２に対して距離値変換処理を行う（ステップ２３）。距離値変換については門脈を構成する各抽出画素に対する背景画素からの最短距離を求める方法を用いる。この公知技術として、「画像理解のためのディジタル画像処理（II）：鳥脇純一郎著（昭晃堂）：３．５距離変換とスケルトン」を挙げる。この公知技術を三次元方向に拡張して利用することで、三次元的な距離値変換を行う。ステップ２３の距離値変換処理とステップ２４の芯線抽出結果とを組み合わせ、芯線を構成する画素の位置における距離値を利用することで、抽出データ（門脈）の血管径を定義することができる。
【００２１】
更に抽出された門脈に対して細線化処理を行う（ステップ２４）。芯線抽出については抽出データに対して細線化処理を行い、細線化結果を芯線として利用する。細線化の手法は一般的に利用されている二次元の細線化手法、例えばＨｉｌｄｉｔｃｈのアルゴリズムを三次元に拡張する手法を用いたり、三次元的な薄面化処理を拡張した細線化手法を用いても良い。ここでも公知技術として「画像理解のためのディジタル画像処理（II）：鳥脇純一郎著（昭晃堂）：アルゴリズム３．５」を挙げることができる。
【００２２】
また、ステップ２３の距離値変換結果を利用し、背景画素との距離が、８近傍の場合は１、１８近傍の場合は１と√２、２６近傍の場合は１と√２と√３、すなわち背景画素と隣り合っている（接している）画素を選択することで、ステップ２４において表面抽出処理を行うこともできる。
【００２３】
抽出した門脈に対してこのような処理を行った後に、ステップ２３、２４の処理結果を利用して肝臓実質３１のうち門脈３２が支配する支配領域特定処理を行う（ステップ２５）。支配領域の特定方法としては、門脈３２の表面画素を利用する方法、門脈３２の芯線の位置を利用する方法、門脈３２の芯線の位置と距離情報（血管径）を合わせて利用する方法等が考えられる。血管径を利用する方法では単純に血管径を利用する方法、血管径を利用して芯線を構成する画素の地点における血管断面積を利用する方法、さらに門脈３２を構成する画素数（血管体積）を利用する方法などが考えられる。ここでは門脈３２の芯線の位置と血管径を合わせて利用する方法を説明する。また、門脈に限らず腫瘍や静脈も同様に支配領域特定処理を行うことが可能である。
【００２４】
抽出された肝臓実質データ集合をＬ＝[Ｌijk]、抽出された門脈３２の芯線データ集合をＰ＝[Ｐijk]、Ｐijkにおける距離値変換値（径）をＲijkと定義すると、肝臓実質構成画素Ｌijkを支配する芯線画素Ｐijkは式１を満たすものとして定義できる。
Ｐijk＝ｍｉｎ（p,q,r）[(Ｌijk-Ｐpqr)2／（α×Ｒpqr）]…(1)
ここでαは係数である。
【００２５】
すなわち、肝臓実質構成画素Ｌijkと芯線画素Ｐpqrとの三次元的距離を元にした値を門脈径Ｒpqrに比例した値で割った相対値が最も小さくなる芯線画素Ｐpqrが、肝臓実質構成画素Ｌijkを支配する芯線画素Ｐijkとして定義されるとになる。この場合、門脈枝５１、５２、５３の芯線毎に設定される肝臓実質の支配領域の境界５４、５５、５６は、図６に示すように血管径が大きい門脈枝５２程境界５５、５６が離れた位置に設定されることになる。なお、上式において、距離値変換値（上記式における分母）を径から近似として計算できる血管断面の縁の長さ（円周）、もしくは断面積を与え、条件式とすることも可能である。
【００２６】
一方、門脈枝の表面情報、もしくは芯線の位置情報のみを利用して同様な支配領域を設定する場合、門脈の表面データ集合をＳ＝[Ｓijk]とすると、肝臓実質構成画素Ｌijkを支配する門脈表面画素Ｓijkは以下の式２を満たすものとして定義される。
Ｓijk＝ｍｉｎ（p,q,r）[(Ｌijk-Ｐpqr)2]…(2)
【００２７】
すなわちこの場合、門脈枝６１、６２、６３毎に設定される肝臓実質の支配領域の境界６４、６５、６６は図７に示すように門脈枝の径によらず、境界６４、６５、６６は各門脈枝との中間位置に設定されることになる。このようにして求められた肝臓実質を支配する画素情報Ｐijk（またはＳijk）を各肝臓実質画素Ｌijk毎に持たせておく。
【００２８】
次に、上記設定した門脈枝のうちの支配領域特定処理によって求めた腫瘍の支配領域を表示する（ステップ２６）。
【００２９】
次に、切除領域を特定する抽出門脈枝、例えばグリソン鞘と呼ばれる部分についての特定処理を行う（ステップ２７）。切除領域を決定するために利用される門脈枝を設定する。門脈枝の指定には前述の領域拡張法を利用して、指定した位置から末梢部までの門脈枝を設定する方法や、３次元的にクリッピング等の処理を利用して門脈枝を切り出す方法等を利用する。
【００３０】
さらに、この処理を全ての門脈枝に対して行い、抽出門脈枝全体をグループ化し、それぞれの門脈枝グループを構成する総画素数を求める。この画素数を式１における分母に利用することで、門脈枝毎の体積による支配領域設定処理を行うこともできる。この場合、ステップ２６の処理をステップ２４の前に行う必要がある。
【００３１】
切除領域データを特定する（ステップ２８）。最後に切り出した門脈枝における芯線画素を検索し、この芯線画素に支配さているものとして定義される肝臓実質画素を検索し、削除することで、所望の切除領域を設定することができる。すなわち、切り出した芯線画素データ集合をＣ＝[Ｃijk]とすると、集合Ｌにおいて集合Ｃの情報を含むデータを検索し、削除する（計算結果として表示しない）または区別できる画素に置き換えて表示する処理を行う。以上の処理により図８に示されるような切除領域７１を設定することができる。
【００３２】
図８は、切除領域（関心領域）７１を支配する血管７２、７３を前処理結果から求めた場合を示す。この血管７２、７３は関心領域７１のみを対象に求められており、必然的に領域７４は関心領域７３を最低限囲むための最小の大きさとなる。すなわち、図８に示す方法は関心領域７１から血管７２、７３にアプローチする方向である。図８では関心領域７１に関与する血管は２本であるが、実際は関心領域７１の近傍には３次元的に複数血管が存在するので血管の特定が困難であることが多いが、図８に示す２次元的なアプローチでは容易に特定することが可能となる。
【００３３】
このようにして得られたデータを表示用データとして合成し、表示データ処理を行う（ステップ２９）。ステップ２８の画像合成については一般的な画像再構成方法である平行投影法を利用したサーフェースレンダリング法やボリュームレンダリング法等を利用する。
【００３４】
処理結果を表示する（ステップ３０）。この時、抽出された特定領域、即ち切除領域とその他の非抽出領域とを画像の濃度値や色相を異ならせて、画像観察者が特定領域をその他の非抽出領域と識別可能に表示すると良い。あるいは切除領域を点滅表示させることで特定領域とその他の非抽出領域と識別可能に表示すると良い。これで一連の処理は完了する。
【００３５】
図９には、本発明のシステムが実現可能であるハードウエア例の構成図を示す。このシステムは、ＣＰＵ９２、主メモリ９０、磁気ディスク９１、表示メモリ９３、ＣＲＴ９４、コントローラ９５、マウス９６、及び共通バス９７から成る。磁気ディスク９１には、各断層像が格納されており、主メモリ９０の投影表示ソフトウェア（図５）に従ってＣＰＵ９２が所定の処理を行う。この処理では、マウス９６やコントローラ９５に付加されているキーボードを利用して入出力処理や処理操作が行われる。積み上げ三次元画像は表示メモリ９３を介してＣＲＴ９４に表示され、オペレータの操作を利用して図５の処理がなされ、閾値条件にあった画像が得られる。また、表示内容は磁気ディスク９１に格納され、再表示に利用される。
【００３６】
以上本発明の実施の形態を説明したが、本発明の手法はＸ線ＣＴ装置だけでなく、磁気共鳴イメージング装置や超音波診断装置などの他の画像診断装置により取得した三次元画像に対しても用いることができる。また、対象臓器としては上記実施の形態中で説明した肝臓の他に人体の多くの部位について適用可能である。
【００３７】
本発明の実施の形態によれば、切除領域と非切除領域とを識別可能に表示することによって、例えば肝臓切除シミュレーション等を行う際に、より臨床に近い形での手術計画、切除シミュレーション、切除率の計算や３次元的な可視化が可能となる。
【００３８】
【発明の効果】
本発明によれば、複数種類の血管等が複雑に入り込んだ臓器の切除非抽出領域と識別可能に表示することができる方法と装置を提供することができる。
【図面の簡単な説明】
【図１】断層像とデータの関係を示す図。
【図２】本発明の処理方法の概要を説明するフローチャート図。
【図３】断層像データからの関心領域抽出を示す図。
【図４】抽出データに対する各種処理を示す図。
【図５】本アルゴリズムの基本処理フローの一例を説明する図。
【図６】本発明における抽出血管の距離値変換処理と細線化結果を組み合わせて利用して領域設定を行った際の領域境界の位置関係を示す図。
【図７】本発明における抽出血管の細線化結果、もしくは表面情報のみを利用して領域設定を行った際の領域境界の位置関係を示す図。
【図８】本発明による切除領域設定結果を示す図。
【図９】本発明を実施可能なハードウエア構成例を示す図。
【符号の説明】
１１…断層像、１２…積み上げ三次元像、３１，３２…抽出データ、５１，５２，５３…抽出血管、５４，５５，５６…支配領域境界、６１，６２，６３…抽出血管、６４，６５，６６…支配領域境界、７１…切除領域
【図１】
0004688361.tif 000002
【図２】
0004688361.tif 000003
【図３】
0004688361.tif 000004
【図４】
0004688361.tif 000005
【図５】
0004688361.tif 000006
【図６】
0004688361.tif 000007
【図７】
0004688361.tif 000008
【図８】
0004688361.tif 000009
【図９】
0004688361.tif 000010