JP6003964B2

JP6003964B2 - 診断装置並びに当該診断装置における画像処理方法及びそのプログラム

Info

Publication number: JP6003964B2
Application number: JP2014227527A
Authority: JP
Inventors: 佐藤　俊次; 俊次佐藤; 中嶋　光康; 光康中嶋; 玲浜田; 田中　勝; 勝田中
Original assignee: Casio Computer Co Ltd
Current assignee: Casio Computer Co Ltd
Priority date: 2014-11-07
Filing date: 2014-11-07
Publication date: 2016-10-05
Anticipated expiration: 2034-11-07
Also published as: EP3023936A1; CA2905061C; EP3023936B1; AU2015230734B2; US9818183B2; US20180040124A1; CA2905061A1; AU2015230734A1; JP2016087272A; US20160133010A1

Description

本発明は、診断装置並びに当該診断装置における画像処理方法及びそのプログラムに関する。

一般に、皮膚病変の診断として視診が必ず行われ、多くの情報を得ることが出来る。しかしながら、肉眼やルーペだけでは、例えば、ホクロとしみの判別さえ難しく、良性腫瘍と悪性腫瘍の鑑別も難しい。このため、従来、ダーモスコープ付きカメラを用いて病変を撮影するダーモスコピー診断が行われている。

ダーモスコープとは、ハロゲンランプ等で病変部を明るく照射し、エコージェルや偏光フィルタなどにより反射光の無い状態に設定し、１０倍程度に拡大して観察する非侵襲性診察器具である。この器具を用いた観察法をダーモスコピー診断と呼んでいる。ダーモスコピー診断については、インターネットＵＲＬ（ｈｔｔｐ：／／ｗｗｗ．ｔｗｍｕ．ａｃ．ｊｐ／ＤＮＨ／ｄｅｐａｒｔｍｅｎｔ／ｄｅｒｍａｔｏｌｏｇｙ／ｄｅｒｍｏｓｃｏｐｙ．ｈｔｍｌ）＜平成２６年９月１日閲覧＞に詳細に記載されている。このダーモスコピー診断によれば、角質による乱反射がなくなることにより、表皮内から真皮浅層までの色素分布が良く見えてくる。

例えば、特許文献１に、上記したダーモスコープで撮影された皮膚画像に対して、色調、テクスチャ、非対称度、円度等の値を用いて診断を行う色素沈着部位の遠隔診断システムの技術が開示されている。特許文献１に開示された技術によれば、ダーモスコープを付したカメラ付き携帯電話を用い、ダーモスコープを通して、メラノーマの心配がある良性色素性母斑などがある皮膚を撮影する。そして、携帯電話のネットワーク接続機能を用いてインターネットに接続し、撮影した皮膚画像を遠隔診断装置に送信して診断を要求するものである。そして、要求に基づき皮膚画像を受信した遠隔診断装置は、メラノーマ診断プログラムを用い、皮膚画像から、それがメラノーマであるか否か、あるいはメラノーマであった場合にどの病期のメラノーマであるかを診断し、その結果を要求者である医師に返信する。

特開２００５−１９２９４４号

ところで、皮膚病については上記したダーモスコープ画像による診断が普及しつつあるが、明瞭な形状変化や模様を得られないことも多く、画像の観察や病変の判断は医師の熟練度に依存しているのが現状である。したがって、ダーモスコープ画像の病変部分を強調する等、画像処理によって、容易、かつ的確に診断可能なツールの出現がのぞまれていた。

本発明は上記した課題を解決するためになされたものであり、医師の診断を容易にするとともに診断精度の向上をはかる、診断装置並びに当該診断装置における画像処理方法及びそのプログラムを提供することを目的とする。

本発明の一態様は、
患部の撮影画像を用いて病変を診断する診断装置であって、
前記撮影画像を記憶する画像記憶部と、
前記画像記憶部に記憶された前記撮影画像を処理する処理部と、を備え、
前記処理部が、
前記撮影画像を輝度成分と色情報成分に分離する分離手段と、
診断の対象の部位を、前記部位のそれらしさを強調するために、前記撮影画像の前記輝度成分又は前記色情報成分により抽出する抽出手段と、を備え、
前記分離手段は、ＲＧＢ色空間の前記撮影画像を、Ｌａｂ色空間に変換し、
前記抽出手段は、Ｌ成分を前記輝度成分とし、ａｂ成分を前記色情報成分とし、前記Ｌａｂ色空間に変換された前記撮影画像から血管領域を尤度Ａとして抽出する
ことを特徴とする診断装置を提供するものである。

本発明によれば、医師の診断を容易にするとともに診断精度の向上をはかる、診断装置並びに当該診断装置における画像処理方法及びそのプログラムを提供することができる。

本発明の実施の形態に係る診断装置の構成を示すブロック図である。本発明の実施の形態に係る診断装置の基本処理動作を示すフローチャートである。図３の血管らしさ抽出処理理動作の詳細を示すフローチャートである。本発明の実施の形態に係る診断装置の表示画面構成の一例を示す図である。

以下、添付図面を参照して、本発明を実施するための形態（以下、「実施形態」と称する）について詳細に説明する。なお、実施形態の説明の全体を通して同じ要素には同じ番号を付している。

（実施形態の構成）
図１は、本実施形態に係る診断装置１００の構成を示すブロック図である。図１に示すように、本実施形態に係る診断装置１００には、ダーモスコープ付き撮影装置１１０が接続されている。ダーモスコープ付き撮影装置１１０は、診断装置１００（処理部１０１）からの指示により患部の撮影を行ない、撮影画像（ダーモスコープ画像）を画像記憶部１０２に格納するとともに表示装置１２０に表示する。また、撮影画像は、処理部１０１により強調処理が施されて画像記憶部１０２に保存されると共に表示装置１２０に表示される。なお、入力装置１３０は、ダーモスコープ画像の撮影開始指示、後述するダーモスコープ画像中の部位選択操作等を行う。

なお、表示装置１２０、例えば、ＬＣＤ（Liquid Crystal Display）モニタであり、入力装置１３０は、例えば、マウス等の入力デバイスである。

処理部１０１は、画像記憶部１０２に記憶されたダーモスコピー画像等、患部の撮影画像を処理するもので、図１に示すように、分離手段１０１ａと、抽出手段１０１ｂと、生成手段１０１ｃと、を含む。

分離手段１０１ａは、撮影画像を輝度成分と色情報成分に分離する手段として機能する。抽出手段１０１ｂは、血管相当部位等、診断の対象の部位を、部位のそれらしさを強調するために、撮影画像の輝度成分又は色情報成分により抽出する手段として機能する。

生成手段１０１ｃは、上記した部位の抽出された結果を背景画像と再合成して生成する手段として機能する。このとき、背景画像は、医師が入力装置１３０を操作することにより選択される、撮影画像、撮影画像のグレー変換画像、撮影画像のコントラスト強調画像、及び撮影画像の輝度成分を大域変動成分と局所変動成分に分離して異なる強調を施した輝度成分強調画像からなる群のうちの少なくとも１つである。

上記した分離手段１０１ａ、抽出手段１０１ｂ、生成手段１０１ｃは、いずれも、処理部１０１が、処理部１０１が有する本実施形態に係るプログラムを逐次読み出し実行することにより、それぞれが持つ上記した機能を実現する。

（実施形態の動作）
以下、図１に示す本実施形態に係る診断装置１００の動作について、図２以降を参照しながら詳細に説明する。

図２に、本実施形態に係る診断装置１００の基本処理動作の流れが示されている。図２によれば、処理部１０１は、まず、ダーモスコープ付き撮影装置１１０で撮影された患部（皮膚病変部位）の撮影画像を取得する（ステップＳ１１）。そして、取得した撮影画像を画像記憶部１０２の所定の領域に格納するとともに、表示装置１２０に表示する（ステップＳ１２）。

図４に表示装置１２０に表示される表示画面イメージの一例が示されている。画面に向かって左に撮影画像表示領域１２１が、右に、例えば、血管の強調画像表示領域１２２が割り当てられ、それぞれの領域に、撮影画像、強調画像が表示される。医師が、図面右下に割り当てられている「撮影開始」ボタン１２３を、入力装置１３０を使用してクリックすることで、ダーモスコープ付き撮影装置１１０による患部の撮影が開始され、処理部１０１が、撮影画像を取得し、表示装置１２０の撮影画像表示領域１２１に表示することができる。

なお、画面左下には、血管領域の表示方法（表示形態）を指定するために用意されたチェックボックス１２４，１２５，１２６が割り付けられており、撮影画像上に赤色で強調して表示する場合、白黒撮影画像上に赤色で強調して表示する場合、血管画像を表示する場合に、医師が入力装置１３０を使用して選択操作（クリック）することで、該当のチェックボックス１２４，１２５，１２６にチェックマークが付される。

説明を図２のフローチャートに戻す。処理部１０１は、取得した撮影画像を表示装置１２０の撮影画像表示領域１２１に表示した後、その撮影画像から血管領域を尤度Ａとして抽出する（ステップＳ１３）。ここで、血管領域は２値の場合、“０”または“１”の値を持ち、多値の場合、０から１の範囲で値を持つ。ここでは、“０”が非血管領域で、“１”が血管領域を表し、この血管領域を尤度Ａとして抽出する処理については後述する。

次に、処理部１０１は、表示装置１２０の表示画面上でのチェックボックス１２４，１２５，１２６の状態から医師により指定される表示方法（表示形態）な関する情報を取得する。医師が、所望のチェックボックス１２４，１２５，１２６のいずれかを選択操作することで、処理部１０１は、その選択された方法（形態）で血管相当部位を強調して表示装置１２０の強調画像表示領域１２２に表示する。

具体的に、チェックボックス１２４が選択されていた場合（ステップＳ１４“撮影画像と合成”）、処理部１０１は、生成手段１０１ｃが、撮影画像ＩＭＧ上に、血管領域を赤色で描画した血管強調画像Ｅを生成する（ステップＳ１５）。ここで、血管強調画像Ｅは、Ｅ＝ＩＭＧ．＊（１−Ａ）＋Ｒｅｄ＊Ａにより生成することができる。ここで、ＩＭＧ．は撮影画像を、Ｒｅｄは赤色を、＊は乗算を、「．＊」は要素毎の乗算を表す。続いて処理部１０１は、生成した血管強調画像Ｅを、表示装置１２０の撮影画像表示領域１２１に表示されている撮影画像とともに強調画像表示領域１２２に並べて表示する（ステップＳ１６）。

一方、チェックボックス１２５が選択されていた場合（ステップＳ１４“グレー画像と合成”）、処理部１０１は、生成手段１０１ｃが、撮影画像をグレー画像に変換した後、血管領域を赤色で描画した血管強調画像Ｅを生成する（ステップＳ１７）。ここで、血管強調画像Ｅは、Ｅ＝Ｇｒｅｙ（ＩＭＧ）．＊（１−Ａ）＋Ｒｅｄ＊Ａにより生成することができる。なお、Ｇｒｅｙ（）は、グレー画像への変換を表している。続いて処理部１０１は、生成した血管強調画像Ｅを、表示装置１２０の撮影画像表示領域１２１に表示されている撮影画像とともに強調画像表示領域１２２に並べて表示する（ステップＳ１６）。

また、チェックボックス１２６が選択されていた場合（ステップＳ１４“血液画像を表示”）、処理部１０１は、生成手段１０１ｃが、血管部分を赤で描画し、非血管部分を黒で描画する血管強調画像を生成する（ステップＳ１８）。ここで、血管強調画像Ｅは、Ｅ＝Ｒｅｄ＊Ａにより生成することができる。なお、血管領域は赤で示しているが、任意であり、赤に限定されない。続いて処理部１０１は、生成した血管強調画像Ｅを、表示装置１２０の撮影画像表示領域１２１に表示されている撮影画像とともに強調画像表示領域１２２に並べて表示する（ステップＳ１６）。

ここで、ステップＳ１３の血管領域を尤度Ａとして抽出する処理について図３のフローチャートを用いて詳述する。図３によれば、処理部１０１は、まず、分離手段１０１ａが、ＲＧＢ色空間の撮影画像を、Ｌａｂ色空間（正確にはＣＩＥ１９７６Ｌ＊ａ＊ｂ色空間）に変換する。Ｌａｂ色空間については、例えば、インターネットＵＲＬ（ｈｔｔｐ：／／Ｊａ．ｗｉｋｉｐｅｄｉａ．ｏｒｇ／ｗｉｋｉ／Ｌａｂ％Ｅ８％８９％Ｂ２％Ｅ７％Ａ９％ＢＡ％Ｅ９％９６％９３）＜平成２６年９月１日閲覧＞にその詳細が記述されている。

続いて、処理部１０１は、抽出手段１０１ｂが、分離手段１０１ａにより取得されたＬａｂ色空間のＬ画像（輝度画像）を利用し、輝度画像Ｌに対して画素毎にヘシアン行列Ｈを取得する（ステップＳ１３２）。ヘシアン行列Ｈは、各画素で、ｘ方向への２階微分、ｙ方向への２階微分、ｘｙ方向への２階部分の３つの要素を持つ。それぞれの要素を、Ｈｘｘ，Ｈｘｙ，Ｈｙｙとすれば、
Ｈｘｘ＝ｓ^２＊（Ｄｘｘ※Ｌ）
Ｈｙｙ＝ｓ^２＊（Ｄｙｙ※Ｌ）
Ｈｘｙ＝ｓ^２＊（Ｄｘｙ※Ｌ）
で示される。

ここで、※は畳み込み（コンポリューション）を示し、Ｄｘｘは、ガウスカーネルをｘ方向に２階微分したもので、
Ｄｘｘ＝１／（２＊ｐｉ＊σ^＊４）＊（（Ｘ．^２）／ｓ^２−１）.＊ｅｘｐ（−（Ｘ．^２＋Ｙ．^２）／（２＊ｓ^２））で与えられる。

ここで、ｓは、検出する血管の太さに応じたスケール値である。なお、複数の太さの血管を抽出する場合には、血管の太さごとに後述する血管らしさＡを求め、画素毎に最大となるＡの値を選択すればよい。ｐｉは円周率、Ｘ及びＹは、カーネル内の位置、「．＊」は、行列の要素毎の乗算、「．^」は要素毎の累乗を表す。

Ｄｙｙは、ガウスカーネルをＹ方向に２階微分したものであり、Ｄｘｘの転置で与えられる。
Ｄｙｙ＝Ｄｘｘ'
ここで、'は転置を表す。

Ｄｘｙは、ガウスカーネルをｘ方向に1階、ｙ方向に1階微分したもので、
Ｄｘｙ＝１／（２＊ｐｉ＊σ^＊６）＊（Ｙ．＊Ｙ）.＊ｅｘｐ（−（Ｘ．^２＋Ｙ．^２）／（２＊ｓ^２））で与えられる。

ヘシアン行列Ｈは、
Ｈ＝［ＨｘｘＨｘｙ；ＨｘｙＨｙｙ］で表される。

次に、抽出手段１０１ｂは、画素毎にヘシアン行列Ｈの固有値λ１、λ２を取得する（ステップＳ１３３)。ヘシアン行列Ｈは、実対象行列であるため、固有値λ１、λ２は必ず実数になり、以下で求めることができる。
λ１＝０．５＊（Ｄｘｘ＋Ｄｙｙ＋ｔｍｐ）
λ２＝０．５＊（Ｄｘｘ；Ｄｙｙ−ｔｍｐ）
ここで、ｔｍｐ＝ｓｑｒｔ（（Ｄｘｘ−Ｄｙｙ）．＊（Ｄｘｘ−Ｄｙｙ）＋４＊Ｄｘｙ．＊Ｄｘｙ）である。

そして、抽出手段１０１ｂは、ステップＳ１３３で取得した各画素での固有値λ１とλ２から以下の式にしたがい血管領域Ａを抽出する（ステップＳ１３５）。血管領域Ａは、Ａ＝１−ｅｘｐ（（Ｋ．^２）／２＊σ^２））で示される。但し、Ｋは、ステップＳ１３４で示されるように、Ｋ＝ｓｑｒｔ（λ１.＊λ１＋λ２.＊λ２）で掲載された値である。ここで、σは調整係数であり、ｓｑｒｔ（）は要素ごとの平方根を表している。

なお、処理部１０１は、血管領域を２値で表す場合、ｔｈ１を閾値として２値化すれば良い。
すなわち、
Ａ＝０ｉｆＡ＜ｔｈ１
Ａ＝１ｅｌｓｅ

（実施形態の効果）
以上説明のように本実施形態に係る診断装置１００によれば、処理部１０１は、分離手段１０１ａが、撮影画像を輝度成分と色情報成分に分離し、抽出手段１０１ｂが、血管相当部位等、診断の対象の部位を、部位のそれらしさを強調するために、撮影画像の輝度成分又は色情報成分により抽出し、抽出された結果を表示装置１２０に表示する。このため、医師は、診断の対象の部位について強調表示された画面を視認することで、容易、かつ的確に診断を行うことができ、その結果、診断精度が向上する。

また、本実施形態に係る診断装置１００によれば、処理部１０１は、生成手段１０１ｃが、部位の抽出された結果を背景画像と再合成して表示装置１２０に表示する。このとき、医師に、背景画像として、（１）撮影画像、（２）撮影画像のグレー変換画像、（３）撮影画像のコントラスト強調画像、及び（４）撮影画像の輝度成分を大域変動成分と局所変動成分に分離して異なる強調を施した輝度成分強調画像からなる群の中から選択できるＵＩ（User Interface）を提供することで、診断目的に応じた表示形態を動的に選択することができ、その結果、医師は、一層容易に、かつ的確に診断を行うことができ、診断精度が一段と向上する。

なお、撮影画像の輝度成分を大域変動成分と局所変動成分に分離して異なる強調を施した輝度成分強調画像が選択されると、分離手段１０１ａが、Ｌａｂ色空間で輝度に相当するＬ画像を用い、大域変動成分と局所変動成分を分離するためにエッジ保存型フィルタ処理を施す。ここで、エッジ保存型フィルタとして、例えば、バイラテラルフィルタを用いる。バイラテラルフィルタについては、例えば、インターネットＵＲＬ（ｈｔｔｐ：／／ｅｎ．ｗｉｋｉｐｅｄｉａ．ｏｒｇ／ｗｉｋｉ／Ｂｉｌａｔｅｒａｌｆｉｌｔｅｒ）＜平成２６年９月１日閲覧＞にその詳細が記載されている。

以上、実施形態を用いて本発明を説明したが、本発明の技術的範囲は上記実施形態に記載の範囲には限定されないことは言うまでもない。上記実施形態及び実施例に、多様な変更又は改良を加えることが可能であることが当業者に明らかである。また、その様な変更又は改良を加えた形態も本発明の技術的範囲に含まれ得ることが、特許請求の範囲の記載から明らかである。

以下に、この出願の願書に最初に添付した特許請求の範囲に記載した発明を付記する。付記に記載した請求項の項番は、この出願の願書に最初に添付した特許請求の範囲の通りである。

〔付記〕
［請求項１］
患部の撮影画像を用いて病変を診断する診断装置であって、
前記撮影画像を記憶する画像記憶部と、
前記画像記憶部に記憶された前記撮影画像を処理する処理部と、を備え、
前記処理部が、
前記撮影画像を輝度成分と色情報成分に分離する分離手段と、
診断の対象の部位を、前記部位のそれらしさを強調するために、前記撮影画像の前記輝度成分又は前記色情報成分により抽出することを特徴とする診断装置。
［請求項２］
前記部位が血管相当部位であることを特徴とする請求項１に記載の診断装置。
［請求項３］
さらに、前記部位の抽出された結果を背景画像と再合成して生成する生成手段を備えることを特徴とする請求項１又は２に記載の診断装置。
［請求項４］
前記背景画像が、前記撮影画像、前記撮影画像のグレー変換画像、前記撮影画像のコントラスト強調画像、及び前記撮影画像の前記輝度成分を大域変動成分と局所変動成分に分離して異なる強調を施した輝度成分強調画像からなる群から選択される少なくとも１つであることを特徴とする請求項３に記載の診断装置。
［請求項５］
患部の撮影画像を用いて病変を診断する診断装置における画像処理方法であって、
前記撮影画像を輝度成分と色情報成分に分離するステップと、
診断の対象の部位を、前記部位のそれらしさを強調するために、前記撮影画像の前記輝度成分又は前記色情報成分により抽出するステップと、を有することを特徴とする画像処理方法。
［請求項６］
患部の撮影画像を用いて病変を診断する診断装置における画像処理のプログラムであって、
コンピュータに、
記憶した前記撮影画像を輝度成分と色情報成分に分離する機能と、
診断の対象の部位を、前記部位のそれらしさを強調するために、前記撮影画像の前記輝度成分又は前記色情報成分により抽出する機能と、を実行させることを特徴とするプログラム。

１００…診断装置、１０１…処理部、１０１ａ…分離手段、１０１ｂ…抽出手段、１０１ｃ…生成手段、１１０…ダーモスコープ付き撮影装置、１２０…表示装置、１２１…撮影画像、１２２…強調画像、１２３…撮影開始ボタン、１２４，１２５，１２６…チェックボックス、１３０…入力装置

Claims

患部の撮影画像を用いて病変を診断する診断装置であって、
前記撮影画像を記憶する画像記憶部と、
前記画像記憶部に記憶された前記撮影画像を処理する処理部と、を備え、
前記処理部が、
前記撮影画像を輝度成分と色情報成分に分離する分離手段と、
診断の対象の部位を、前記部位のそれらしさを強調するために、前記撮影画像の前記輝度成分又は前記色情報成分により抽出する抽出手段と、を備え、
前記分離手段は、ＲＧＢ色空間の前記撮影画像を、Ｌａｂ色空間に変換し、
前記抽出手段は、Ｌ成分を前記輝度成分とし、ａｂ成分を前記色情報成分とし、前記Ｌａｂ色空間に変換された前記撮影画像から血管領域を尤度Ａとして抽出する
ことを特徴とする診断装置。
前記部位が血管相当部位であることを特徴とする請求項１に記載の診断装置。
前記部位が血管相当部位であり、前記部位のそれらしさは、血管らしさであることを特徴とする請求項１に記載の診断装置。
さらに、前記部位の抽出された結果を背景画像と再合成して生成する生成手段を備えることを特徴とする請求項１ないし３いずれかに記載の診断装置。
前記背景画像が、前記撮影画像、前記撮影画像のグレー変換画像、前記撮影画像のコントラスト強調画像、及び前記撮影画像の前記輝度成分を大域変動成分と局所変動成分に分離して異なる強調を施した輝度成分強調画像からなる群から選択される少なくとも１つであることを特徴とする請求項４に記載の診断装置。
前記抽出手段は、前記Ｌａｂ色空間に変換された前記撮影画像から血管領域を尤度Ａとして抽出し、前記尤度Ａは、
Ａ＝１−ｅｘｐ（（Ｋ．^２）／２＊σ^２））
で示されることを特徴とする請求項１に記載の診断装置。
但し、
Ｋ＝ｓｑｒｔ（λ１.＊λ１＋λ２.＊λ２）
である。ここで、σは調整係数であり、ｓｑｒｔ（）は要素ごとの平方根を表す。λ１、λ２は各画素に対応したヘシアン行列（Hessian matrix）の固有値を表す実数である。
前記尤度Ａは血管のそれらしさを示すことを特徴とする請求項１又は６に記載の診断装置。
前記抽出手段は、複数の太さの血管を抽出する場合には、血管の太さに対応した調整係数σにより血管らしさＡを求め、画素毎に最大となるＡの値を選択することを特徴とする請求項６に記載の診断装置。
患部の撮影画像を用いて病変を診断する診断装置における画像処理方法であって、
前記撮影画像を輝度成分と色情報成分に分離するステップと、
診断の対象の部位を、前記部位のそれらしさを強調するために、前記診断装置が、前記撮影画像の前記輝度成分又は前記色情報成分により抽出するステップと、を有し、
前記分離するステップは、ＲＧＢ色空間の前記撮影画像を、Ｌａｂ色空間に変換し、
前記抽出するステップは、Ｌ成分を前記輝度成分とし、ａｂ成分を前記色情報成分とし、前記Ｌａｂ色空間に変換された前記撮影画像から血管領域を尤度Ａとして抽出する
ことを特徴とする診断装置における画像処理方法。
患部の撮影画像を用いて病変を診断する診断装置におけるプログラムであって、
コンピュータに、
記憶した前記撮影画像を輝度成分と色情報成分に分離する機能と、
診断の対象の部位を、前記部位のそれらしさを強調するために、前記撮影画像の前記輝度成分又は前記色情報成分により抽出する機能と、を実行させ、
前記分離する機能は、ＲＧＢ色空間の前記撮影画像を、Ｌａｂ色空間に変換し、
前記抽出する機能は、Ｌ成分を前記輝度成分とし、ａｂ成分を前記色情報成分とし、前記Ｌａｂ色空間に変換された前記撮影画像から血管領域を尤度Ａとして抽出する
ことを特徴とするプログラム。