(57)【要約】【課題】医用画像データから、所定の構造物に対応する木構造を精度よく作成する。【解決手段】医用画像データから、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる所定の構造物を起始部に対応する根ノードＮ0を含む複数のノードＮiと複数のエッジＥによって定義した木構造の候補として抽出し、抽出した木構造の候補の各ノードＮについて、所定の構造物が根ノードから分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて各ノードＮiと根ノードＮ0との関係を評価することにより前記各ノード間の接続されやすさを重み付けするコスト関数を用いて、木構造の定義を再評価して木構造を作成する。【選択図】図３
2012223338.tif 000002
【特許請求の範囲】
【請求項１】
医用画像データから、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる所定の構造物を前記起始部に対応する根ノードを含む複数のノードと複数のエッジによって定義した木構造の候補として抽出する木構造候補抽出手段と、
前記抽出した木構造の候補の各ノードについて、前記所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて前記各ノードと前記根ノードとの関係を評価することにより前記各ノード間の接続されやすさを重み付けするコスト関数を用いて、前記木構造の定義を再評価して木構造を作成する木構造作成手段とを備えたことを特徴とする木構造作成装置。
【請求項２】
前記コスト関数は、前記根ノードと各ノードとを結ぶ線分と、前記各ノードと該ノードの前記離れる方向に隣接するノード間のエッジとのなす角度に基づいて前記角度が小さい程前記各ノードと前記隣接するノードが接続されやすくなるよう重み付けを行うものであることを特徴とする請求項１記載の木構造作成装置。
【請求項３】
前記コスト関数は、前記根ノードと各ノードとを結ぶ線分の長さと、前記各ノードと該ノードの前記離れる方向に隣接するノード間のエッジの長さとの差に基づいて前記長さの差が小さい程前記各ノードと前記隣接するノードが接続されやすくなるよう重み付けを行うものであることを特徴とする請求項１記載の木構造作成装置。
【請求項４】
前記起始部が既知である木構造を表す複数の教師データを機械学習することにより、前記起始部を検出する起始部検出手段をさらに備えたことを特徴とする請求項１から３のいずれか１項記載の木構造作成装置。
【請求項５】
前記所定の構造物は、肺の血管であることを特徴とする請求項１から４いずれか１項記載の木構造作成装置。
【請求項６】
前記所定の構造物は、肝臓の血管であることを特徴とする請求項１から４いずれか１項記載の木構造作成装置。
【請求項７】
医用画像データから、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる所定の構造物を前記起始部に対応する根ノードを含む複数のノードと複数のエッジによって定義した木構造の候補として抽出し、
前記抽出した木構造の候補の各ノードについて、前記所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて前記各ノードと前記根ノードとの関係を評価することにより前記各ノード間の接続されやすさを重み付けするコスト関数を用いて、前記木構造の定義を再評価して木構造を作成することを特徴とする木構造作成方法。
【請求項８】
医用画像データから、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる所定の構造物を前記起始部に対応する根ノードを含む複数のノードと複数のエッジによって定義した木構造の候補として抽出する木構造候補抽出手段と、
前記抽出した木構造の候補の各ノードについて、前記所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて前記各ノードと前記根ノードとの関係を評価することにより前記各ノード間の接続されやすさを重み付けするコスト関数を用いて、前記木構造の定義を再評価して木構造を作成する木構造作成手段とを備えたことを特徴とする木構造作成プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像データから検出した特定の構造物を木構造として作成する木構造作成装置および方法ならびにプログラムに関するものである。
【背景技術】
【０００２】
肝臓・肺などの臓器において患部を切除する手術を行う際には、例えば肝臓の場合、肝臓のＸ線ＣＴ画像から血管や肝臓実質、腫瘍部分を抽出し、抽出された血管の芯線の位置や血管径等に基づいて、抽出された腫瘍がどの血管の支配領域に属するかを特定することによって腫瘍に栄養を供給する血管を特定し、その血管の支配領域を切除すべき部分として適切に決定することが求められる。肝臓の切除手術においては、腫瘍に栄養を送っている門脈の部分と、その門脈の部分によって支配される、がん細胞などの注意すべき物質が送られている可能性のある領域を、肝機能を維持できる範囲で適切に切除するためである。このため、術前にどこを切除部位とすべきか綿密にシミュレートすることが重要であり、このシミュレートのために肺や肝臓を走行する血管の中心経路を精度よく抽出することが必要とされている。
【０００３】
ここで、ＣＴ等で得られた３次元医用画像から気管支等の線状構造を木構造として抽出する画像認識技術として、ヘッセ行列を用いた手法が提案されている。具体的には、まず、３次元医用画像に対して多重解像度変換を行った後、各解像度の画像に対してヘッセ行列の固有値解析を行い、線状構造要素を抽出する。この線状構造要素は、固有値解析によって得られる３つの固有値のうち１つだけが０に近いという特徴を有する。次に、各解像度の画像における解析結果を統合することによって、３次元医用画像中の様々なサイズの線状構造要素（血管）を抽出する。そして、最小全域木アルゴリズム等を用いて、抽出された各線状構造要素を連結することにより、３次元医用画像中の管状構造を表す木構造のデータが得られる。なお、最小全域木アルゴリズムによる線状構造要素の連結の際には、各線状構造要素の位置関係や、上記０に近い固有値に対応する固有ベクトルで表される、各線状構造要素の主軸方向に基づいたコスト関数が用いられる（特許文献１）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−２２０７４２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかし、肝臓の門脈と肝静脈のように、枝分かれ回数が多く、各血管の枝同士が近い位置を絡み合うように走行している構造物では、特許文献１と同様の手法により門脈などの血管を抽出しようとすると、別の血管である肝静脈の枝を門脈の枝として誤って抽出してしまうことがある。
【０００６】
門脈の枝や肺血管の枝は、起始部から分岐を繰り返しながら離れる方向に広がって延びる点で、他の血管とは性質が異なる。本発明においては、この性質を利用して、起始部から分岐を繰り返しながら離れる方向に広がって延びる血管を従来方法よりも効率よく且つ正確に抽出する木構造作成装置および方法ならびにプログラムを提供することを目的とするものである。
【課題を解決するための手段】
【０００７】
本発明における木構造作成装置は、医用画像データから、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる所定の構造物を前記起始部に対応する根ノードを含む複数のノードと複数のエッジによって定義した木構造の候補として抽出する木構造候補抽出手段と、前記抽出した木構造の候補の各ノードについて、前記所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて前記各ノードと前記根ノードとの関係を評価することにより前記各ノード間の接続されやすさを重み付けするコスト関数を用いて、前記木構造の定義を再評価して木構造を作成する木構造作成手段とを備えたことを特徴とするものである。
【０００８】
本発明における木構造作成方法は、医用画像データから、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる所定の構造物を前記起始部に対応する根ノードを含む複数のノードと複数のエッジによって定義した木構造の候補として抽出し、前記抽出した木構造の候補の各ノードについて、前記所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて前記各ノードと前記根ノードとの関係を評価することにより前記各ノード間の接続されやすさを重み付けするコスト関数を用いて、前記木構造の定義を再評価して木構造を作成することを特徴とするものである。
【０００９】
本発明における木構造作成プログラムは、コンピュータを、医用画像データから、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる所定の構造物を前記起始部に対応する根ノードを含む複数のノードと複数のエッジによって定義した木構造の候補として抽出する木構造候補抽出手段と、前記抽出した木構造の候補の各ノードについて、前記所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて前記各ノードと前記根ノードとの関係を評価することにより前記各ノード間の接続されやすさを重み付けするコスト関数を用いて、前記木構造の定義を再評価して木構造を作成する木構造作成手段として機能させることを特徴とするものである。
【００１０】
ここで、所定の構造物とは、ノード（節点）とそれを結ぶエッジ（辺）とにより根つき木構造として形状モデル化可能なオブジェクトであるとともに、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる構造物であれば何でもよい。例えば、複数の所定の構造物は、肺または肝臓の血管であってもよく、さらに、肺における肺動脈、肺静脈、肝臓における門脈、冠動脈、肝静脈などが考えられる。
【００１１】
また、医用画像データとは、たとえばＣＴ、ＭＲ、超音波装置、ＰＥＴ―ＣＴ、ＳＰＥＣＴ、４Ｄ-ＣＴ、ＯＣＴ、Ｘ線撮影装置（ＣＲ、ＤＲ）により撮像された医用画像データであってもよいし、たとえばボリュームデータ等の３次元画像データであってもよい。
【００１２】
また、所定の構造物を木構造の候補として抽出する方法には、所定の構造物を木構造として抽出できるものであれば周知のいかなる方法も適用でき、例えば、本発明の木構造作成装置が、画像データから所定の構造物の画像的特徴を有する領域を検出し、検出された領域中を細線化して得られた細線を分岐点および所定の距離などにより分割することにより抽出してもよい。
【００１３】
「木構造を作成する」とは、それぞれのノードごとに、各ノード間の接続の強さを表す指標値をコストとして算出することにより、各ノード間の接続の強さが強くなるように各ノードを接続できるものであればいかなる木構造作成方法を適用してもよい。例えば、最小全域木法や最短路木法などのスパニング木生成アルゴリズムを適用することが考えられる。
【００１４】
上記本発明の木構造作成装置における前記コスト関数は、所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて前記各ノードと前記根ノードとの関係を評価することにより前記各ノード間の接続されやすさを重み付けするものであれば、いかなるものであってもよい。例えば、前記根ノードと各ノードとを結ぶ線分と、前記各ノードと該ノードの前記離れる方向に隣接するノード間のエッジとのなす角度に基づいて前記角度が小さい程前記各ノードと前記隣接するノードが接続されやすくなるよう重み付けを行うものであることが好ましい。また、上記本発明の木構造作成装置における前記コスト関数は、前記根ノードと各ノードとを結ぶ線分の長さと、前記各ノードと該ノードの前記離れる方向に隣接するノード間のエッジの長さとの差に基づいて前記長さの差が小さい程前記各ノードと前記隣接するノードが接続されやすくなるよう重み付けを行うものであることが好ましい。
【００１５】
また、起始部は、任意の方法により特定したものであってよく、起始部に基づいて周知の方法により起始部に対応する根ノードを特定してよい。例えば、表示された画像上でマウス等の入力装置により起始部を指定してもよく、起始部が既知である所定の構造物を表す複数の教師データを機械学習することにより、前記起始部を検出する起始部検出手段をさらに備え、かかる起始部検出手段により検出してもよい。なお、教師データを機械学習して根ノードを抽出する周知の種々の方法を用いることができ、例えば、アダブースト(Adaboost)法により教師データにおいて既知である起始部の特徴量に基づいて起始部を検出することが考えられる。
【発明の効果】
【００１６】
本発明の木構造作成装置および方法ならびにプログラムによれば、医用画像データから、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる所定の構造物を前記起始部に対応する根ノードを含む複数のノードと複数のエッジによって定義した木構造を抽出し、抽出した木構造の各ノードについて、所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて前記各ノードと前記根ノードとの関係を評価することにより前記各ノード間の接続されやすさを重み付けするため、起始部から分岐を繰り返しながら離れる方向に広がって延びる構造物の幾何学的特徴を利用して木構造の定義を再評価して木構造を作成することにより、効率よく且つ正確に木構造を作成することができる。
【図面の簡単な説明】
【００１７】
【図１】本実施形態における木構造作成装置の機能ブロック図
【図２】構造物として検出された肝臓の血管領域の木構造候補抽出処理を説明する図
【図３】本実施形態における各ノード間のコストの算出方法を説明する図
【図４】抽出された肝臓の門脈と肝静脈を表す図
【図５】第２の実施形態における木構造作成装置の機能ブロック図
【発明を実施するための形態】
【００１８】
以下、図面を参照して本発明の木構造作成装置の実施の形態を詳細に説明する。図１は本発明の第１の実施形態となる木構造作成装置１の概略構成図である。なお、図１のような木構造作成装置１の構成は、補助記憶装置に読み込まれた木構造作成プログラムをコンピュータ上で実行することにより実現される。このとき、この木構造作成プログラムは、ＣＤ−ＲＯＭ等の記憶媒体に記憶され、もしくはインターネット等のネットワークを介して配布され、コンピュータにインストールされる。図１の木構造作成装置１は、肝臓における門脈等の構造物を表す画像データから門脈を表す木構造を作成するものであって、構造物領域検出手段１１、木構造候補抽出手段１２、木構造作成手段１３と表示制御手段１４とを備えている。また、本実施形態における木構造作成プログラムがインストールされたコンピュータは、木構造作成装置１として機能する本体と、ディスプレイからなる表示装置３とマウスやキーボード等の入力装置４とメモリやハードディスク等から構成される記憶手段２を備えている。
【００１９】
構造物領域検出手段１１は、画像データにおいて所定の構造物の一部を構成するか否かを判断することにより候補領域Ｒｃを検出するものである。なお、画像データ７は、データ記憶手段２に記憶されたたとえば撮影装置もしくは放射線検出装置により撮像された２次元の画像もしくは複数の２次元画像から生成された３次元のボリュームデータからなっている。図２は、所定の構造物領域である血管を、血管の候補領域Ｒｃとして抽出し、抽出された血管領域に基づいてグラフ化して木構造の候補を抽出する様子を示す模式図である。
【００２０】
ここで、一例として、所定の構造物が肝臓における門脈であり、ボリュームデータから門脈の候補領域Ｒｃを検出する場合について説明する。
【００２１】
まず、構造物領域検出手段１１は、図２(i)に示すように、ボリュームデータ７を構成するボクセルデータの値に基づいて、門脈の芯線を構成する複数の候補点Ｓｐ（ｐ＝１〜ｎ：ｎは候補点の抽出数）の位置を算出する。なお、ここでは、かかる位置の算出に先立って門脈を表す画素であることが既知であるボクセルデータの値を統計等により取得し、そのボクセルデータの値により門脈らしい画素を候補点Ｓｐとして判断する。
【００２２】
そして、構造物領域検出手段１１は、図２(ii)に示すように、候補点を拡張する。さらに、図２(iii)に示すように、拡張された候補点Ｓｐの門脈を表す画素値を含む所定の範囲の画素値を有する画像データを、門脈領域であると判断し候補領域Ｒｃとして抽出する。
【００２３】
なお、構造物領域検出手段１１は、本実施形態に限られず、候補領域Ｒｃを抽出可能な周知の種々の方法を適用してよい。例えば、候補点周辺のボクセルデータについて門脈らしさを表す特徴量を算出し、算出した特徴量に基づいてそのボクセルデータが門脈領域を表すものであるか否かを判別してもよく、かかる場合、特徴量に基づく判別は、マシンラーニングにより予め取得された評価関数に基づいて行なうことが考えられる。あるいは、構造物領域検出手段１１は特開２０１０−２００９２５号もしくは特許文献１に開示された手法やその他公知の技術により候補領域を検出してもよい。
【００２４】
次に、木構造候補抽出手段１２は、１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる所定の構造物（門脈）を起始部に対応する根ノードを含む複数のノードと複数のエッジによって定義した木構造の候補として抽出する。本実施形態では、まず、木構造候補抽出手段１２は、図２(iii)に示す構造物領域検出手段１１が検出した構造物領域Ｒｃを取得し、図２(iv)に示すように、取得した構造物領域Ｒｃを周知の方法により細線化し、図２(v)に示すように、細線化処理された線を分岐点において分割する。そして、分岐点および端点を複数のノードNi（ｉ＝１〜ｍ：ｍは候補点の抽出数）として定義するとともに、分割された線分を複数のエッジＥとして定義して木構造の候補を抽出する。なお、細線化処理された線を分岐点だけでなく、所定の間隔などの所定の条件で分割してもよい。細線化処理された線のうち、ゆるやかな湾曲形状をなす部分を、湾曲に沿って適宜複数の線分に分割できるようにするためである。
【００２５】
木構造作成手段１３は、抽出した木構造の候補の各ノードNi（ｉ＝１〜ｍ：ｍは候補点の抽出数）について、各ノード間の接続されやすさを重み付けすることにより木構造の定義を再評価して木構造Mを作成する。図３を用いて木構造作成手段１３により各ノードN間のコストを決定する方法の原理を説明する。図３は、上述の木構造候補抽出処理により抽出された門脈をあらわす木構造の候補を表した図である。
【００２６】
ここで、本発明者は、肝臓の門脈などの所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びる幾何学的な特徴を有することに着目した。そして、所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて各ノードと根ノードとの関係を評価することにより各ノード間の接続されやすさを重み付けするコスト関数を用いて、木構造の定義を再評価して好適に木構造を作成可能であることを見いだした。
【００２７】
本実施形態では、図３に示すように、上記幾何学的特徴を、所定の構造物を木構造として表した木構造においては各ノードについて、根ノードNoと各ノードNiとを結ぶ線分と、各ノードNiと該ノードの離れる方向に隣接するノードNi+1間のエッジＥとのなす角度θが小さいという特徴としてさらに捉え、コスト関数を定義した。
【００２８】
さらに詳細には、上記角度θを下記式（１）のように内積によって判断し、この式（１）を用いて角度θが小さい程各ノードNiとノードNi+1が接続されやすくなるよう重み付けを行うことにより、上記幾何学的な特徴を木構造の作成に反映した。
【数１】
2012223338.tif 000003
【００２９】
本実施形態では、上記観点に基づいて、角度θが小さい程各ノードNiとノードNi+1が接続されやすくなるように、式（１）をコスト関数ｆとして定義する。そして、コスト関数ｆにより、Ｎ個のノードについてコストを個々に算出する。なお、本実施形態では、マウス等の入力装置をマニュアル操作することにより指定された起始部の位置を門脈の起始部に対応する根ノードN0として入力する。
【００３０】
そして、木構造作成手段１３は、上記算出された各コストに基づいて各ノードを再接続してmaxΣｆとなるように最小全域木法により木構造を作成する。各ノードのコストを評価するコスト関数に基づいて木構造作成する周知の種々の方法を用いることができ、例えば、周知の最小全域木法や最短路木法などのスパニング木生成アルゴリズムにより、maxΣｆとなるように、最適経路を決定し、木構造Ｍを作成することができる。
【００３１】
図４は、上記処理により作成された木構造Ｍに基づいて抽出された門脈Ｍ１と、同様に別途作成された木構造に基づいて抽出された肝静脈Ｍ２を表す図である。図４に示すように、表示制御手段１４は、ディスプレイ３に、木構造作成手段１３により作成された木構造Ｍに基づいて周知の方法により門脈Ｍ１を表示させる。また、本実施形態では、肺静脈についても門脈同様に木構造作成処理を行い、その結果得られた肝静脈Ｍ２をディスプレイ３に表示させる。
【００３２】
上記第１の実施の形態によれば、肝臓の門脈と肝静脈のように、枝分かれ回数が多く、各血管の枝同士が近い位置を絡み合うように走行している構造物においても、起始部から分岐を繰り返しながら離れる方向に広がって延びる構造物の幾何学的特徴を利用して木構造の定義を再評価して木構造を作成するため、効率よく且つ正確に木構造を作成することができる。
【００３３】
また、上記起始部から分岐を繰り返しながら離れる方向に広がって延びる構造物の幾何学的特徴を根ノードと各ノードを結ぶ線分と、各根ノードと根ノードから分岐を繰り返しながら離れる方向に隣接するノードとの内積により規定することにより、簡易かつ好適に上記特徴を各ノードの接続されやすさに重み付けして木構造を作成することができる。
【００３４】
また、本実施形態の応用例として、式（１）に替えて、所定の構造物を木構造として表した木構造においては各接続されるべき候補点と、その候補点との端点と根ノードとを結ぶ線分との相対的な角度が小さいという特徴を表すいかなる評価方法を用いてもよい。
【００３５】
また、本実施形態に限られず、本発明の木構造作成装置におけるコスト関数は、所定の構造物が１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという幾何学的特徴に基づいて各ノードと根ノードとの関係を評価することにより各ノード間の接続されやすさを重み付けするものであれば、いかなるものであってもよい。例えば、コスト関数は、式（２）に示すように、根ノードと各ノードとを結ぶ線分の長さと、各ノードとそのノードの根ノードから分岐を繰り返しながら離れる方向に隣接するノード間のエッジの長さとの差に基づいて長さの差が小さい程各ノードと前記隣接するノードが接続されやすくなるよう重み付けを行うものであってもよい。
【数２】
2012223338.tif 000004
【００３６】
上記根ノードと各ノードとを結ぶ線分の長さと、ノードとそのノードの根ノードから分岐を繰り返しながら離れる方向に隣接するノード間のエッジの長さとの差に基づいて、長さの差が小さい程各ノードと前記隣接するノードが接続されやすくなるよう重み付けを行う場合には、簡易かつ好適に上記特徴を各ノードの接続されやすさとして重み付けして木構造を作成することができる。
【００３７】
また、上記実施形態に限定されず、コスト関数ｆに、各ノード同士の距離が近いほど接続強度が大きくなるようにさらに重み付けを行ってもよく、さらなる周知の重み付け方法を任意に組み合わせてもよい。
【００３８】
また、第２の実施形態として、図５に示すように、木構造作成装置１は起始部が既知である門脈を表す複数の教師データを機械学習することにより、門脈の起始部を検出する起始部検出手段１５をさらに備え、かかる起始部検出手段１５により検出された起始部を取得し、取得した起始部に対応する根ノードを、周知の方法により特定して上記木構造作成処理を行ってもよい。第２の実施形態では、アダブースト(Adaboost)法により教師データにおいて既知である起始部の特徴量に基づいて起始部を検出し、画像データの座標系における起始部の座標値を取得する。そして、起始部の座標値に基づいて根ノードを特定する。かかる場合には、起始部の入力のためのマニュアル操作を省略することができるため、より効率よく木構造の作成を行うことができる。また、アダブースト(Adaboost)法により起始部検出を行った場合には、所定の構造物から好適に起始部を検出することができる。なお、第２の実施形態に限定されず、本発明に起始部を抽出する周知の種々の方法を適用可能である。
【００３９】
また、ここでは所定の構造物として門脈を基に説明したが、所定の構造物は点とそれを結ぶ候補点とにより木構造として形状モデル化可能なオブジェクトであるとともに１つの起始部から分岐を繰り返しながら離れる方向に広がって延びるという特徴を有するものであれば何でもよい。例えば、所定の構造物は、肺または肝臓の血管であってよく、さらには、肺動脈、肺静脈、肝臓における門脈、肝動脈、肝静脈があげられる。
【符号の説明】
【００４０】
１     木構造作成装置
７     画像データ
１１   構造物領域検出手段
１２   木構造候補抽出手段
１３   木構造作成手段
１４   表示制御手段
１５   起始部検出手段
Ｓｐ   候補点
Ｎ０    根ノード
Ｎ     ノード
Ｅ     エッジ
ｆ     コスト関数
【図１】
2012223338.tif 000005
【図２】
2012223338.tif 000006
【図３】
2012223338.tif 000007
【図５】
2012223338.tif 000008
【図４】
2012223338.tif 000009