JP4748843B2

JP4748843B2 - 臓器の特定領域抽出表示装置及び臓器の特定領域抽出表示方法をコンピュータに実行させるプログラムを記録した記録媒体

Info

Publication number: JP4748843B2
Application number: JP2000334687A
Authority: JP
Inventors: 朋洋永尾; 行矣三浦; 宣夫中尾
Original assignee: Hitachi Medical Corp
Current assignee: Hitachi Healthcare Manufacturing Ltd
Priority date: 2000-01-25
Filing date: 2000-11-01
Publication date: 2011-08-17
Anticipated expiration: 2020-11-01
Also published as: JP2001283191A

Description

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

Claims

医用画像診断装置にて得た被検体の画像を用いて臓器の特定領域を抽出し表示装置に表示する装置において、
前記画像中の目的臓器が撮影された領域を関心領域として抽出する関心領域抽出手段と、
前記関心領域内の領域、かつ当該関心領域に含まれる血管により支配される支配領域を特定する支配領域特定手段であって、前記血管の血管径が大きいほど離れた位置に前記支配領域の境界を設定する、又は複数の血管の間に位置する前記支配領域の境界を前記複数の血管の中間位置に設定する、ことにより前記支配領域を特定する支配領域特定手段と、
前記支配領域において、前記血管上の指定された位置から先の血管が支配する領域を特定領域として抽出する特定領域抽出手段と、
前記特定領域と異なる前記関心領域内の他の領域を非抽出領域とし、該非抽出領域と前記特定領域とを識別可能に前記表示装置へ表示する表示手段と、
を備えたことを特徴とする臓器の特定領域抽出表示装置。
前記特定領域と前記非抽出領域との表示は、非抽出領域を表示しない表示モードと、前記特定領域を非抽出領域と識別して表示する表示モードとを設定し、それらの表示モードを切り替えて表示する、
ことを特徴とする請求項１に記載の臓器の特定領域抽出表示装置。
前記医用画像診断装置にて得た被検体の画像は、三次元画像データであって、前記表示手段は、前記目的臓器の関心領域を基に形成された前記目的臓器の三次元画像内の前記特定領域と前記非抽出領域とに識別可能に表示すると同時に、前記目的臓器の三次元画像に対して設定された任意の切断面における２次元断面像内の前記特定領域と前記非抽出領域とを識別可能に表示する、
ことを特徴とする請求項１又は２に記載の臓器の特定領域抽出表示装置。
医用画像診断装置にて得た被検体の画像を用いて臓器の特定領域を抽出し表示装置に表示する方法をコンピュータに実行させるプログラムにおいて、
前記画像中の目的臓器が撮影された領域を関心領域として抽出するステップと、
前記関心領域内の領域、かつ当該関心領域に含まれる血管により支配される支配領域を特定するステップであって、前記血管の血管径が大きいほど離れた位置に前記支配領域の境界を設定する、又は複数の血管の間に位置する前記支配領域の境界を前記複数の血管の中間位置に設定する、ことにより前記支配領域を特定するステップと、
前記支配領域において、前記血管上の指定された位置から先の血管が支配する領域を特定領域として抽出するステップと、
前記特定領域と異なる前記関心領域内の他の領域を非抽出領域とし、該非抽出領域と前記特定領域とを識別可能に前記表示装置へ表示するステップと、
を含むことを特徴とする臓器の特定領域抽出表示方法をコンピュータに実行させるプログラムを記録した記録媒体。