特開 2024-72096 画像に含まれる複数の輪郭の分析方法、画像に含まれる複数の輪郭を分析するためのシステム、及び画像に含まれる複数の輪郭を分析するためのプログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公開特許公報(A)

(11)【公開番号】P2024072096

(43)【公開日】2024-05-27

(54)【発明の名称】画像に含まれる複数の輪郭の分析方法、画像に含まれる複数の輪郭を分析するためのシステム、及び画像に含まれる複数の輪郭を分析するためのプログラム

(51)【国際特許分類】

   G06T 7/62 20170101AFI20240520BHJP

【ＦＩ】

G06T7/62

【審査請求】未請求

【請求項の数】19

【出願形態】ＯＬ

(21)【出願番号】P 2022182744

(22)【出願日】2022-11-15

(71)【出願人】

【識別番号】000000033

【氏名又は名称】旭化成株式会社

(74)【代理人】

【識別番号】100079108

【弁理士】

【氏名又は名称】稲葉 良幸

(74)【代理人】

【識別番号】100109346

【弁理士】

【氏名又は名称】大貫 敏史

(74)【代理人】

【識別番号】100117189

【弁理士】

【氏名又は名称】江口 昭彦

(74)【代理人】

【識別番号】100134120

【弁理士】

【氏名又は名称】内藤 和彦

(72)【発明者】

【氏名】根本 ほのか

(72)【発明者】

【氏名】松下 忠史

(72)【発明者】

【氏名】山口 陽平

【テーマコード（参考）】

5L096

【Ｆターム（参考）】

5L096BA13

5L096BA20

5L096FA06

5L096FA65

5L096FA70

5L096GA17

(57)【要約】

【課題】画像に含まれる複数の輪郭の特徴量を評価可能な方法を提供する。
【解決手段】画像の一辺を含む１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出し、１番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、一辺に対して垂直方向における１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より大きい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を狭くしてから、１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より小さい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を広げてから、１番目領域の幅の中央を、１番目領域の中心とすることを含む、画像に含まれる複数の輪郭の分析方法。
【選択図】図２１

【特許請求の範囲】

【請求項1】

（１）複数の輪郭を含む画像を取得することと、
（２）前記画像の一辺を含む１番目領域を特定することと、
（３）前記１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出することと、
（４）前記１番目領域に含まれる複数の輪郭の前記径の変動係数が所定の範囲内にある場合、前記一辺に対して垂直方向における前記１番目領域の幅の中央を、前記１番目領域の中心とし、前記径の変動係数が前記所定の範囲より大きい場合、前記径の変動係数が前記所定の範囲内になるまで前記１番目領域の幅を狭くしてから、前記１番目領域の幅の中央を、前記１番目領域の中心とし、前記径の変動係数が前記所定の範囲より小さい場合、前記径の変動係数が前記所定の範囲内になるまで前記１番目領域の幅を広げてから、前記１番目領域の幅の中央を、前記１番目領域の中心とすることと、
（５）ｎを０を初期値とする整数として、前記一辺に対して垂直方向において、ｎ＋１番目領域の中心から一定の間隔離れたｎ＋２番目領域の中心を設定することと、
（６）前記ｎ＋２番目領域の幅を、前記ｎ＋１番目領域に含まれる複数の輪郭の最大径以上にすることと、
（７）前記ｎ＋２番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出することと、
（８）前記ｎ＋２番目領域に含まれる複数の輪郭の前記径の変動係数が前記所定の範囲内にある場合、前記一辺に対して垂直方向における前記ｎ＋２番目領域の幅を保ち、前記径の変動係数が前記所定の範囲より大きい場合、前記ｎ＋２番目領域の中心を保ったまま前記径の変動係数が前記所定の範囲内になるまで前記ｎ＋２番目領域の幅を狭くし、前記径の変動係数が前記所定の範囲より小さい場合、前記ｎ＋２番目領域の中心を保ったまま前記径の変動係数が前記所定の範囲内になるまで前記ｎ＋２番目領域の幅を広げることと、
（９）ｎに１を加算して（５）から（８）を実施することを、前記ｎ＋２番目領域が前記画像の前記一辺に対向する他辺を超えるまで繰り返すことと、
（１０）前記ｎ＋１番目領域に含まれる前記複数の輪郭の特徴量を取得することと、
を含む、画像に含まれる複数の輪郭の分析方法。

【請求項2】

前記画像が正方形又は矩形である、請求項１に記載の方法。

【請求項3】

前記画像が二値化又は多値化されている、請求項１に記載の方法。

【請求項4】

前記画像が断面画像又は表面画像である、請求項１に記載の方法。

【請求項5】

前記画像が、電子顕微鏡画像、光学顕微鏡画像、プローブ顕微鏡画像、又は放射線画像である、請求項１に記載の方法。

【請求項6】

前記複数の輪郭のそれぞれが、孔、樹脂骨格、海島構造における島、海島構造における海、又はネットワーク構造を示している、請求項１に記載の方法。

【請求項7】

前記ｎ＋２番目領域の幅を、前記ｎ＋１番目領域に含まれる複数の輪郭の最大径の５倍以上にする、請求項１に記載の方法。

【請求項8】

前記特徴量が統計値である、請求項１に記載の方法。

【請求項9】

前記特徴量が、径、平均径、周囲長、及び平均周囲長の少なくともいずれかである、請求項１に記載の方法。

【請求項10】

複数の輪郭を含む画像を取得するための画像取得部と、
前記画像の一辺を含む１番目領域を特定し、前記１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出し、前記１番目領域に含まれる複数の輪郭の前記径の変動係数が所定の範囲内にある場合、前記一辺に対して垂直方向における前記１番目領域の幅の中央を、前記１番目領域の中心とし、前記径の変動係数が前記所定の範囲より大きい場合、前記径の変動係数が前記所定の範囲内になるまで前記１番目領域の幅を狭くしてから、前記１番目領域の幅の中央を、前記１番目領域の中心とし、前記径の変動係数が前記所定の範囲より小さい場合、前記径の変動係数が前記所定の範囲内になるまで前記１番目領域の幅を広げてから、前記１番目領域の幅の中央を、前記１番目領域の中心とする処理を実行するための第１設定部と、
ｎを０を初期値とする整数として、前記一辺に対して垂直方向において、ｎ＋１番目領域の中心から一定の間隔離れたｎ＋２番目領域の中心を設定し、前記ｎ＋２番目領域の幅を、前記ｎ＋１番目領域に含まれる複数の輪郭の最大径以上にし、前記ｎ＋２番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出し、前記ｎ＋２番目領域に含まれる複数の輪郭の前記径の変動係数が前記所定の範囲内にある場合、前記一辺に対して垂直方向における前記ｎ＋２番目領域の幅を保ち、前記径の変動係数が前記所定の範囲より大きい場合、前記ｎ＋２番目領域の中心を保ったまま前記径の変動係数が前記所定の範囲内になるまで前記ｎ＋２番目領域の幅を狭くし、前記径の変動係数が前記所定の範囲より小さい場合、前記ｎ＋２番目領域の中心を保ったまま前記径の変動係数が前記所定の範囲内になるまで前記ｎ＋２番目領域の幅を広げる処理を実行するための第２設定部と、
前記第２設定部に、ｎに１を加算して前記処理を実行することを、前記ｎ＋２番目領域が前記画像の前記一辺に対向する他辺を超えるまで繰り返させる繰り返し指示部と、
前記ｎ＋１番目領域に含まれる前記複数の輪郭の特徴量を取得する特徴量取得部と、
を備える、画像に含まれる複数の輪郭を分析するためのシステム。

【請求項11】

前記画像が正方形又は矩形である、請求項１０に記載のシステム。

【請求項12】

前記画像を二値化又は多値化する多値化部をさらに備える、請求項１０に記載のシステム。

【請求項13】

前記画像が断面画像又は表面画像である、請求項１０に記載のシステム。

【請求項14】

前記画像が、電子顕微鏡画像、光学顕微鏡画像、プローブ顕微鏡画像、又は放射線画像である、請求項１０に記載のシステム。

【請求項15】

前記複数の輪郭のそれぞれが、孔、樹脂骨格、海島構造における島、海島構造における海、又はネットワーク構造を示している、請求項１０に記載のシステム。

【請求項16】

前記第２設定部が、前記ｎ＋２番目領域の幅を、前記ｎ＋１番目領域に含まれる複数の輪郭の最大径の５倍以上にする、請求項１０に記載のシステム。

【請求項17】

前記特徴量が統計値である、請求項１０に記載のシステム。

【請求項18】

前記特徴量が、径、平均径、周囲長、及び平均周囲長の少なくともいずれかである、請求項１０に記載のシステム。

【請求項19】

コンピューターに、
（１）複数の輪郭を含む画像を取得することと、
（２）前記画像の一辺を含む１番目領域を特定することと、
（３）前記１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出することと、
（４）前記１番目領域に含まれる複数の輪郭の前記径の変動係数が所定の範囲内にある場合、前記一辺に対して垂直方向における前記１番目領域の幅の中央を、前記１番目領域の中心とし、前記径の変動係数が前記所定の範囲より大きい場合、前記径の変動係数が前記所定の範囲内になるまで前記１番目領域の幅を狭くしてから、前記１番目領域の幅の中央を、前記１番目領域の中心とし、前記径の変動係数が前記所定の範囲より小さい場合、前記径の変動係数が前記所定の範囲内になるまで前記１番目領域の幅を広げてから、前記１番目領域の幅の中央を、前記１番目領域の中心とすることと、
（５）ｎを０を初期値とする整数として、前記一辺に対して垂直方向において、ｎ＋１番目領域の中心から一定の間隔離れたｎ＋２番目領域の中心を設定することと、
（６）前記ｎ＋２番目領域の幅を、前記ｎ＋１番目領域に含まれる複数の輪郭の最大径以上にすることと、
（７）前記ｎ＋２番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出することと、
（８）前記ｎ＋２番目領域に含まれる複数の輪郭の前記径の変動係数が前記所定の範囲内にある場合、前記一辺に対して垂直方向における前記ｎ＋２番目領域の幅を保ち、前記径の変動係数が前記所定の範囲より大きい場合、前記ｎ＋２番目領域の中心を保ったまま前記径の変動係数が前記所定の範囲内になるまで前記ｎ＋２番目領域の幅を狭くし、前記径の変動係数が前記所定の範囲より小さい場合、前記ｎ＋２番目領域の中心を保ったまま前記径の変動係数が前記所定の範囲内になるまで前記ｎ＋２番目領域の幅を広げることと、
（９）ｎに１を加算して（５）から（８）を実施することを、前記ｎ＋２番目領域が前記画像の前記一辺に対向する他辺を超えるまで繰り返すことと、
（１０）前記ｎ＋１番目領域に含まれる前記複数の輪郭の特徴量を取得することと、
を含む手順を実行させるための、画像に含まれる複数の輪郭を分析するためのプログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、画像分析技術に関し、画像に含まれる複数の輪郭の分析方法、画像に含まれる複数の輪郭を分析するためのシステム、及び画像に含まれる複数の輪郭を分析するためのプログラムに関する。

【背景技術】

【0002】

画像に含まれる文字や輪郭を分析するための様々な方法が提案されている(例えば、特許文献１、２及び非特許文献１参照。)。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開平５－１８９６１０号公報

【特許文献2】特開平９－３０５７５４号公報

【非特許文献】

【0004】

【非特許文献1】T. Hildebrand, et al., “A new method for the model-independent assessment of thickness in three-dimensional images,” Journal of Microscopy, Volume185, Issue 1, January 1997, Pages 67-75

【発明の概要】

【発明が解決しようとする課題】

【0005】

本発明は、画像に含まれる複数の輪郭の特徴量を評価可能な方法、システム、及びプログラムを提供することを課題の一つとする。

【課題を解決するための手段】

【0006】

［１］（１）複数の輪郭を含む画像を取得することと、（２）画像の一辺を含む１番目領域を特定することと、（３）１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出することと、（４）１番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、一辺に対して垂直方向における１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より大きい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を狭くしてから、１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より小さい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を広げてから、１番目領域の幅の中央を、１番目領域の中心とすることと、（５）ｎを０を初期値とする整数として、一辺に対して垂直方向において、ｎ＋１番目領域の中心から一定の間隔離れたｎ＋２番目領域の中心を設定することと、（６）ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径以上にすることと、（７）ｎ＋２番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出することと、（８）ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、一辺に対して垂直方向におけるｎ＋２番目領域の幅を保ち、径の変動係数が所定の範囲より大きい場合、ｎ＋２番目領域の中心を保ったまま径の変動係数が所定の範囲内になるまでｎ＋２番目領域の幅を狭くし、径の変動係数が所定の範囲より小さい場合、ｎ＋２番目領域の中心を保ったまま径の変動係数が所定の範囲内になるまでｎ＋２番目領域の幅を広げることと、（９）ｎに１を加算して（５）から（８）を実施することを、ｎ＋２番目領域が画像の一辺に対向する他辺を超えるまで繰り返すことと、（１０）ｎ＋１番目領域に含まれる複数の輪郭の特徴量を取得することと、を含む、画像に含まれる複数の輪郭の分析方法。

【0007】

［２］画像が正方形又は矩形である、［１］に記載の方法。

【0008】

［３］画像が二値化又は多値化されている、［１］又は［２］に記載の方法。

【0009】

［４］画像が断面画像又は表面画像である、［１］から［３］のいずれかに記載の方法。

【0010】

［５］画像が、電子顕微鏡画像、光学顕微鏡画像、プローブ顕微鏡画像、又は放射線画像である、［１］から［４］のいずれかに記載の方法。

【0011】

［６］複数の輪郭のそれぞれが、孔、樹脂骨格、海島構造における島、海島構造における海、又はネットワーク構造を示している、［１］から［５］のいずれかに記載の方法。

【0012】

［７］ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径の５倍以上にする、［１］から［６］のいずれかに記載の方法。

【0013】

［８］特徴量が統計値である、［１］から［７］のいずれかに記載の方法。

【0014】

［９］特徴量が、径、平均径、周囲長、及び平均周囲長の少なくともいずれかである、［１］から［８］のいずれかに記載の方法。

【0015】

［１０］複数の輪郭を含む画像を取得するための画像取得部と、画像の一辺を含む１番目領域を特定し、１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出し、１番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、一辺に対して垂直方向における１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より大きい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を狭くしてから、１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より小さい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を広げてから、１番目領域の幅の中央を、１番目領域の中心とする処理を実行するための第１設定部と、ｎを０を初期値とする整数として、一辺に対して垂直方向において、ｎ＋１番目領域の中心から一定の間隔離れたｎ＋２番目領域の中心を設定し、ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径以上にし、ｎ＋２番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出し、ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、一辺に対して垂直方向におけるｎ＋２番目領域の幅を保ち、径の変動係数が所定の範囲より大きい場合、ｎ＋２番目領域の中心を保ったまま径の変動係数が所定の範囲内になるまでｎ＋２番目領域の幅を狭くし、径の変動係数が所定の範囲より小さい場合、ｎ＋２番目領域の中心を保ったまま径の変動係数が所定の範囲内になるまでｎ＋２番目領域の幅を広げる処理を実行するための第２設定部と、第２設定部に、ｎに１を加算して処理を実行することを、ｎ＋２番目領域が画像の一辺に対向する他辺を超えるまで繰り返させる繰り返し指示部と、ｎ＋１番目領域に含まれる複数の輪郭の特徴量を取得する特徴量取得部と、を備える、画像に含まれる複数の輪郭を分析するためのシステム。

【0016】

［１１］画像が正方形又は矩形である、［１０］に記載のシステム。

【0017】

［１２］画像を二値化又は多値化する多値化部をさらに備える、［１０］又は［１１］に記載のシステム。

【0018】

［１３］画像が断面画像又は表面画像である、［１０］から［１２］のいずれかに記載のシステム。

【0019】

［１４］画像が、電子顕微鏡画像、光学顕微鏡画像、プローブ顕微鏡画像、又は放射線画像である、［１０］から［１３］のいずれかに記載のシステム。

【0020】

［１５］複数の輪郭のそれぞれが、孔、樹脂骨格、海島構造における島、海島構造における海、又はネットワーク構造を示している、［１０］から［１４］のいずれかに記載のシステム。

【0021】

［１６］第２設定部が、ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径の５倍以上にする、［１０］から［１５］のいずれかに記載のシステム。

【0022】

［１７］特徴量が統計値である、［１０］から［１６］のいずれかに記載のシステム。

【0023】

［１８］特徴量が、径、平均径、周囲長、及び平均周囲長の少なくともいずれかである、［１０］から［１７］のいずれかに記載のシステム。

【0024】

［１９］コンピューターに、（１）複数の輪郭を含む画像を取得することと、（２）画像の一辺を含む１番目領域を特定することと、（３）１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出することと、（４）１番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、一辺に対して垂直方向における１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より大きい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を狭くしてから、１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より小さい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を広げてから、１番目領域の幅の中央を、１番目領域の中心とすることと、（５）ｎを０を初期値とする整数として、一辺に対して垂直方向において、ｎ＋１番目領域の中心から一定の間隔離れたｎ＋２番目領域の中心を設定することと、（６）ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径以上にすることと、（７）ｎ＋２番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出することと、（８）ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、一辺に対して垂直方向におけるｎ＋２番目領域の幅を保ち、径の変動係数が所定の範囲より大きい場合、ｎ＋２番目領域の中心を保ったまま径の変動係数が所定の範囲内になるまでｎ＋２番目領域の幅を狭くし、径の変動係数が所定の範囲より小さい場合、ｎ＋２番目領域の中心を保ったまま径の変動係数が所定の範囲内になるまでｎ＋２番目領域の幅を広げることと、（９）ｎに１を加算して（５）から（８）を実施することを、ｎ＋２番目領域が画像の一辺に対向する他辺を超えるまで繰り返すことと、（１０）ｎ＋１番目領域に含まれる複数の輪郭の特徴量を取得することと、を含む手順を実行させるための、画像に含まれる複数の輪郭を分析するためのプログラム。

【0025】

［２０］画像が正方形又は矩形である、［１９］に記載のプログラム。

【0026】

［２１］画像が二値化又は多値化されている、［１９］又は［２０］に記載のプログラム。

【0027】

［２２］画像が断面画像又は表面画像である、［１９］から［２１］のいずれかに記載のプログラム。

【0028】

［２３］画像が、電子顕微鏡画像、光学顕微鏡画像、プローブ顕微鏡画像、又は放射線画像である、［１９］から［２２］のいずれかに記載のプログラム。

【0029】

［２４］複数の輪郭のそれぞれが、孔、樹脂骨格、海島構造における島、海島構造における海、又はネットワーク構造を示している、［１９］から［２３］のいずれかに記載のプログラム。

【0030】

［２５］ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径の５倍以上にする、［１９］から［２４］のいずれかに記載のプログラム。

【0031】

［２６］特徴量が統計値である、［１９］から［２５］のいずれかに記載のプログラム。

【0032】

［２７］特徴量が、径、平均径、周囲長、及び平均周囲長の少なくともいずれかである、［１９］から［２６］のいずれかに記載のプログラム。

【発明の効果】

【0033】

本発明によれば、画像に含まれる複数の輪郭の特徴量を評価可能な方法、システム、及びプログラムを提供可能である。

【図面の簡単な説明】

【0034】

【図1】実施形態に係る画像に含まれる複数の輪郭を分析するためのシステムを示す模式図である。

【図2】実施形態に係る画像の例を示す図である。

【図3】実施形態に係る画像処理の例を示す図である。

【図4】実施形態に係る画像処理の例を示す図である。

【図5】実施形態に係るＬｏｃａｌ ｔｈｉｃｋｎｅｓｓ法の例を示す図である。

【図6】実施形態に係るＬｏｃａｌ ｔｈｉｃｋｎｅｓｓ法の例を示す図である。

【図7】実施形態に係る輪郭検出法の例を示す図である。

【図8】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図9】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図10】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図11】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図12】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図13】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図14】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図15】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図16】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図17】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図18】実施形態に係る画像に領域を設定する方法の例を示す図である。

【図19】実施形態に係るプロットの例を示す図である。

【図20】実施形態に係るプロットの例を示す図である。

【図21】実施形態に係る画像に含まれる複数の輪郭の分析方法のフローチャートである。

【図22】実施形態に係る画像に含まれる複数の輪郭の分析方法のフローチャートである。

【図23】実施形態に係る画像に含まれる複数の輪郭の分析方法のフローチャートである。

【発明を実施するための形態】

【0035】

以下、本発明を実施するための形態（以下、「本実施形態」という。）について詳細に説明する。なお、本実施形態は、本発明の理解を容易にするためのものであり、本発明を限定して解釈するためのものではない。本発明は、本実施形態に限定されるものではなく、その要旨の範囲内で種々変形して実施することができる。

【0036】

実施形態に係る画像に含まれる複数の輪郭を分析するためのシステムは、図１に示すように、複数の輪郭を含む画像を取得するための画像取得部３０１と、画像の一辺を含む１番目領域を特定し、１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出し、１番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、一辺に対して垂直方向における１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より大きい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を狭くしてから、１番目領域の幅の中央を、１番目領域の中心とし、径の変動係数が所定の範囲より小さい場合、径の変動係数が所定の範囲内になるまで１番目領域の幅を広げてから、１番目領域の幅の中央を、１番目領域の中心とする処理を実行するための第１設定部３０２と、を備える。

【0037】

実施形態に係る画像に含まれる複数の輪郭を分析するためのシステムは、ｎを０を初期値とする整数として、一辺に対して垂直方向において、ｎ＋１番目領域の中心から一定の間隔離れたｎ＋２番目領域の中心を設定し、ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径以上にし、ｎ＋２番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出し、ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、一辺に対して垂直方向におけるｎ＋２番目領域の幅を保ち、径の変動係数が所定の範囲より大きい場合、ｎ＋２番目領域の中心を保ったまま径の変動係数が所定の範囲内になるまでｎ＋２番目領域の幅を狭くし、径の変動係数が所定の範囲より小さい場合、ｎ＋２番目領域の中心を保ったまま径の変動係数が所定の範囲内になるまでｎ＋２番目領域の幅を広げる処理を実行するための第２設定部３０３をさらに備える。

【0038】

実施形態に係る画像に含まれる複数の輪郭を分析するためのシステムは、第２設定部３０３に、ｎに１を加算して処理を実行することを、ｎ＋２番目領域が画像の一辺に対向する他辺を超えるまで繰り返させる繰り返し指示部３０４と、ｎ＋１番目領域に含まれる複数の輪郭の特徴量を取得する特徴量取得部３０５と、をさらに備える。

【0039】

画像取得部３０１、第１設定部３０２、第２設定部３０３、繰り返し指示部３０４、及び特徴量取得部３０５は、例えば、中央演算処理装置（ＣＰＵ）３００に含まれる。ＣＰＵ３００には、キーボード、及びマウス等の入力装置３５１と、ディスプレイ等の出力装置３５２と、情報を記憶するメモリー３５３と、が接続されていてもよい。

【0040】

画像取得部３０１が取得する画像は、例えば、電子顕微鏡画像、光学顕微鏡画像、プローブ顕微鏡画像、及び放射線画像であるが、特に限定されない。電子顕微鏡の例としては、透過型電子顕微鏡（ＴＥＭ）、及び走査型電子顕微鏡（ＳＥＭ）が挙げられる。光学顕微鏡の例としては、位相差顕微鏡、微分干渉顕微鏡、偏光顕微鏡、蛍光顕微鏡、共焦点レーザー顕微鏡、全反射照明蛍光顕微鏡、ラマン顕微鏡、及び非線形光学顕微鏡が挙げられる。プローブ顕微鏡の例としては、走査型トンネル顕微鏡（ＳＴＭ）、及び原子間力顕微鏡（ＡＦＭ）が挙げられる。放射線画像の例としては、レントゲン画像、ＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）画像、放射性同位元素（ＲＩ）画像、及び血管造影画像が挙げられる。また、画像は、磁気共鳴（ＭＲ）画像であってもよい。画像は、例えば、正方形又は矩形である。画像は、断面画像であってもよいし、表面画像であってもよい。画像は、例えば、被写体に対して一定方向に連続して撮影された複数の画像を１枚に合成した画像であってもよい。また、例えば、画像は、ノイズ除去処理された画像であってもよい。

【0041】

画像に含まれる複数の輪郭のそれぞれは、例えば、孔、樹脂骨格、海島構造における島、海島構造における海、又はネットワーク構造を示していてもよい。以下においては、画像の被写体が多孔質膜の断面である例を説明するが、本実施形態はこれに限定されない。

【0042】

実施形態に係るシステムは、画像を二値化又は多値化する多値化部３０６をさらに備えていてもよい。二値化処理においては、各画素において、所定の閾値以上の濃度値を白に変換し、所定の閾値未満の濃度値を黒に変換する。多値化処理においては、各画素において、ｍを自然数として、所定のｍ個の閾値を設定し、濃度値とｍ個の閾値の大小関係に基づいて、濃度値をｍ＋１個の色のいずれかに変換する。二値化はｍが１である場合の多値化の形態である。二値化法は、大津の二値化法であってもよい。多値化部３０６が多値化した画像を第１設定部３０２及び第２設定部３０３が処理してもよいし、第１設定部３０２及び第２設定部３０３が特定したｎ＋１番目領域を、多値化部３０６が多値化してもよい。

【0043】

第１設定部３０２は、図２に示すように、画像２００の一辺２０１を含む１番目領域２１０を特定する。正方形又は矩形の画像２００の外郭は、４つの辺２０１、２０２、２０３、２０４からなり、第１設定部３０２は、４つの辺２０１、２０２、２０３、２０４から一辺を特定する。第１設定部３０２は、入力装置３５１を介したオペレーターからの指示に基づき、画像２００の一辺２０１を特定してもよい。１番目領域２１０は、正方形又は矩形の領域であり、１番目領域２１０の外郭の一辺は、画像２００の外郭の一辺と重複する。１番目領域２１０は、画像２００の一辺２０１に対して垂直方向に、任意の幅を有する。第１設定部３０２は、入力装置３５１を介したオペレーターからの指示に基づき、１番目領域２１０の任意の幅を設定してもよい。

【0044】

図３に示すように、第１設定部３０２が特定した１番目領域２１０を、多値化部３０６が多値化してもよい。第１設定部３０２は、画像の１番目領域２１０又は多値化された画像の１番目領域２１１を解析して、１番目領域２１０に含まれる複数の輪郭の最大径及び径の変動係数を算出する。解析方法の例としては、Ｌｏｃａｌ ｔｈｉｃｋｎｅｓｓ法及び輪郭検出法が挙げられる。径の変動係数Ｖは、例えば、径の標準偏差をＳ、径の平均値をＡとして、下記（１）式で与えられる。
Ｖ＝Ｓ／Ａ （１）

【0045】

Ｌｏｃａｌ ｔｈｉｃｋｎｅｓｓ法においては、図４に示すように、輪郭で囲まれた領域に入る最大の円の大きさを求めて、当該円の径から、輪郭で囲まれた領域の径を求める。なお、図５に示すように、円５０１を入れられる輪郭５００で囲まれた領域は円である必要はない。また、図６に示すように、輪郭５０２で囲まれた領域は、線分で完全に囲まれた閉じた領域である必要はない。Ｌｏｃａｌ ｔｈｉｃｋｎｅｓｓ法を用いる場合、第１設定部３０２は、輪郭で囲まれた複数の領域の径を算出し、複数の径の値のうち最大値を最大径として特定する。また、複数の径の値から径の標準偏差と径の平均値を算出し、例えば上記（１）式に従って、径の変動係数を算出する。

【0046】

輪郭検出法においては、図７に示すように、閉じた領域の輪郭５０５を検出する。例えば、画像が二値化画像であれば、白と黒の境界を検出する。第１設定部３０２は、検出された輪郭で囲まれた閉じた領域の面積Ｃを算出し、例えば下記（２）式に従って、領域の円相当径Ｒを算出する。
Ｒ＝２（Ｃ／π）^１／２ （２）
輪郭検出法を用いる場合、第１設定部３０２は、輪郭で囲まれた複数の領域の円相当径を算出し、複数の円相当径の値のうち最大値を最大径として特定する。また、複数の円相当径の値から円相当径の標準偏差と円相当径の平均値を算出し、例えば上記（１）式に従って、径の変動係数を算出する。

【0047】

第１設定部３０２は、領域の外周Ｌに基づいて、例えば下記（３）式に従って、領域の円形度Ｄを算出してもよい。
Ｄ＝（４πＣ）／Ｌ^２ （３）
また、第１設定部３０２は、領域に外接する四角形の短辺長に対する長辺長の比、あるいは領域に外接する円形の短径に対する長径の比を算出してもよい。これらの値は、１から離れるほど、輪郭が歪んでいることを示している。

【0048】

１番目領域２１０における径の変動係数が大きい場合は、１番目領域２１０の幅が広く、１番目領域２１０における径の分布が大きいことを意味し得る。また、１番目領域２１０における径の変動係数が小さい場合は、１番目領域２１０の幅が狭く、１番目領域２１０における径の分布が小さいことを意味し得る。

【0049】

本実施形態においては、画像を、中心が一定の間隔で離れている複数の領域に抽出し、それぞれの領域に含まれる複数の輪郭の特徴量を取得することを目的としているため、１つの領域に含まれる径の分布が大き過ぎたり、小さ過ぎたりする場合は、領域の位置に対する特徴量の変化を正確に把握することが困難になり得る。また、１つの領域における径の変動係数が小さ過ぎる場合は、１つの領域に含まれる輪郭が少なすぎる可能性や、１つの領域の幅が輪郭で囲まれた領域の径よりも短い可能性がある。

【0050】

適切な変動係数の大きさは、取得しようとする特徴量や画像の被写体の大きさに基づいて変わり得るため、任意である。適切な変動係数の大きさに基づいて設定される変動係数の所定の範囲は、例えば、０．２５以上０．５以下、あるいは０．３０以上０．４５以下であるが、これらに限定されない。変動係数の所定の範囲は、例えば、オペレーターが、図１に示す入力装置３５１を介して、メモリー３５３に記憶させていてもよい。

【0051】

第１設定部３０２は、メモリー３５３から、変動係数の所定の範囲を読み出し、１番目領域に含まれる複数の輪郭の径の変動係数の値が、所定の範囲内にあるか否かを判定する。径の変動係数の値が、所定の範囲内にある場合、１番目領域２１０の幅が適切であったことを意味し得る。そのため、第１設定部３０２は、図８に示すように、１番目領域２１０の幅を変更せずに、１番目領域２１０の幅の中央を、１番目領域２１０の中心２１５に設定する。

【0052】

径の変動係数が所定の範囲の上限より大きい場合、第１設定部３０２は、径の変動係数が所定の範囲内になるまで、１番目領域２１０の幅を狭め、径の変動係数を計算することを繰り返す。第１設定部３０２が、１番目領域２１０の幅を狭める方法は、特に限定されない。例えば、１番目領域２１０の現在の幅に１未満の所定の値をかけてもよいし、１番目領域２１０の現在の幅から所定の値を引いてもよい。あるいは、第１設定部３０２は、例えば、１番目領域に含まれる複数の輪郭の最大径に１より大きいの所定の値をかけて得られる値で、かつ、１番目領域の現在の幅より小さい値に、１番目領域の幅を狭めて、径の変動係数を計算することを繰り返してもよい。また、あるいは、第１設定部３０２は、１番目領域に含まれる複数の輪郭の最大径に所定の値を足して得られる値で、かつ、１番目領域の現在の幅より小さい値に、１番目領域の幅を狭めて、径の変動係数を計算することを繰り返してもよい。第１設定部３０２は、幅を狭める回数に応じて、狭める量を増やしてもよいし、減らしてもよい。第１設定部３０２は、径の変動係数が所定の範囲内になったときの１番目領域２１０の幅の中央を、１番目領域の中心に設定する。

【0053】

径の変動係数が所定の範囲の上限より小さい場合、第１設定部３０２は、径の変動係数が所定の範囲内になるまで、１番目領域２１０の幅を広げ、径の変動係数を計算することを繰り返す。第１設定部３０２が、１番目領域２１０の幅を広げる方法は、特に限定されない。例えば、１番目領域２１０の現在の幅に１より大の所定の値をかけてもよいし、１番目領域２１０の現在の幅に所定の値を足してもよい。あるいは、第１設定部３０２は、例えば、１番目領域に含まれる複数の輪郭の最大径に１より大きいの所定の値をかけて得られる値で、かつ、１番目領域の現在の幅より大きい値に、１番目領域の幅を広げて、径の変動係数を計算することを繰り返してもよい。また、あるいは、第１設定部３０２は、１番目領域に含まれる複数の輪郭の最大径に所定の値を足して得られる値で、かつ、１番目領域の現在の幅より大きい値に、１番目領域の幅を広げて、径の変動係数を計算することを繰り返してもよい。第１設定部３０２は、幅を広げる回数に応じて、広げる量を増やしてもよいし、減らしてもよい。第１設定部３０２は、径の変動係数が所定の範囲内になったときの１番目領域２１０の幅の中央を、１番目領域の中心に設定する。

【0054】

第２設定部３０３は、ｎを０を初期値とする整数として、画像２００の一辺２０１に対して垂直方向において、ｎ＋１番目領域の中心から一定の間隔離れたｎ＋２番目領域の中心を設定する。１番目領域の中心を設定した時点で、ｎは０である。一定の間隔は任意である。一定の間隔は、画像２００におけるピクセル数で表されてもよい。一定の間隔は、画像２００から何個程度の領域に抽出するかに応じて、任意に設定してもよい。一定の間隔は、例えば、オペレーターが、入力装置３５１を介して、メモリー３５３に記憶させていてもよく、第２設定部３０３は、メモリー３５３に保存されている一定の間隔を読み出してもよい。ｎが０である場合、第２設定部３０３は、図９に示すように、１番目領域の中心２１５から一定の間隔Ｇ離れた位置に、２番目領域の中心２２５を設定する。

【0055】

第２設定部３０３は、ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径以上に設定する。なお、ｎ＋２番目領域は、ｎ＋１番目領域と重複してもよいし、ｎ＋１番目領域と重複しなくてもよい。ｎ＋２番目領域の幅は、ｎ＋１番目領域に含まれる複数の輪郭の最大径以上であれば、特に限定されないが、例えば、ｎ＋２番目領域の幅は、ｎ＋１番目領域に含まれる複数の輪郭の最大径の５倍、７倍、１０倍、又は１５倍である。ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径以上に設定に設定することにより、ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数の値が小さくなり過ぎることを抑制可能である。ｎが０である場合、第２設定部３０３は、図１０に示すように、２番目領域２２０の幅を、１番目領域に含まれる複数の輪郭の最大径以上に設定する。

【0056】

第２設定部３０３は、ｎ＋２番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出し、ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数の値が、所定の範囲内にあるか否かを判定する。径の変動係数の値が、所定の範囲内にある場合、ｎ＋２番目領域の幅が適切であったことを意味し得る。そのため、第１設定部３０２は、ｎ＋２番目領域の幅を変更せずに保つ。

【0057】

径の変動係数が所定の範囲の上限より大きい場合、第２設定部３０３は、ｎ＋２番目領域の中心を保ったまま、径の変動係数が所定の範囲内になるまで、ｎ＋２番目領域の幅を狭め、径の変動係数を計算することを繰り返す。第２設定部３０３が、ｎ＋２番目領域の幅を狭める方法は、特に限定されない。例えば、ｎ＋２番目領域の現在の幅に１未満の所定の値をかけてもよいし、ｎ＋２番目領域の現在の幅から所定の値を引いてもよい。あるいは、第２設定部３０３は、例えば、ｎ＋２番目領域の中心を保ったまま、ｎ＋２番目領域に含まれる複数の輪郭の最大径に１より大きいの所定の値をかけて得られる値で、かつ、ｎ＋２番目領域の現在の幅より小さい値に、ｎ＋２番目領域の幅を狭めて、径の変動係数を計算することを繰り返してもよい。また、あるいは、第２設定部３０３は、ｎ＋２番目領域の中心を保ったまま、ｎ＋２番目領域に含まれる複数の輪郭の最大径に所定の値を足して得られる値で、かつ、ｎ＋２番目領域の現在の幅より小さい値に、ｎ＋２番目領域の幅を狭めて、径の変動係数を計算することを繰り返してもよい。第２設定部３０３は、幅を狭める回数に応じて、狭める量を増やしてもよいし、減らしてもよい。第２設定部３０３は、径の変動係数が所定の範囲内になったときのｎ＋２番目領域の幅を、ｎ＋２番目領域の幅と確定する。

【0058】

径の変動係数が所定の範囲の上限より小さい場合、第２設定部３０３は、ｎ＋２番目領域の中心を保ったまま、径の変動係数が所定の範囲内になるまで、ｎ＋２番目領域の幅を広げ、径の変動係数を計算することを繰り返す。第２設定部３０３が、ｎ＋２番目領域の幅を広げる方法は、特に限定されない。例えば、ｎ＋２番目領域の現在の幅に１より大きいの所定の値をかけてもよいし、ｎ＋２番目領域の現在の幅に所定の値を足してもよい。あるいは、第２設定部３０３は、例えば、ｎ＋２番目領域の中心を保ったまま、ｎ＋２番目領域に含まれる複数の輪郭の最大径に１より大きいの所定の値をかけて得られる値で、かつ、ｎ＋２番目領域の現在の幅より大きい値に、ｎ＋２番目領域の幅を広げて、径の変動係数を計算することを繰り返してもよい。また、あるいは、第２設定部３０３は、ｎ＋２番目領域の中心を保ったまま、ｎ＋２番目領域に含まれる複数の輪郭の最大径に所定の値を足して得られる値で、かつ、ｎ＋２番目領域の現在の幅より大きい値に、ｎ＋２番目領域の幅を広げて、径の変動係数を計算することを繰り返してもよい。第２設定部３０３は、幅を狭める回数に応じて、広げる量を増やしてもよいし、減らしてもよい。第２設定部３０３は、径の変動係数が所定の範囲内になったときのｎ＋２番目領域の幅を、ｎ＋２番目領域の幅と確定する。

【0059】

繰り返し指示部３０４は、第２設定部３０３に、ｎに１を加算して処理を実行することを、ｎ＋２番目領域が画像２００の一辺２０１に対向する他辺２０３を超えるまで繰り返させる。これにより、画像２００に複数の領域が設定され、各領域の中心は、隣接する領域の中心から所定の一定の間隔Ｇ離れている。

【0060】

例えば、図１１に示すように、ｎを１として、画像２００の一辺２０１に対して垂直方向において、２番目領域２２０の中心２２５から一定の間隔Ｇ離れた位置に中心２３５を有する３番目領域２３０が設定される。次に、図１２に示すように、ｎを２として、画像２００の一辺２０１に対して垂直方向において、３番目領域２３０の中心２３５から一定の間隔Ｇ離れた位置に中心２４５を有する４番目領域２４０が設定される。次に、図１３に示すように、ｎを３として、画像２００の一辺２０１に対して垂直方向において、４番目領域２４０の中心２４５から一定の間隔Ｇ離れた位置に中心２５５を有する５番目領域２５０が設定される。

【0061】

次に、図１４に示すように、ｎを４として、画像２００の一辺２０１に対して垂直方向において、５番目領域２５０の中心２５５から一定の間隔Ｇ離れた位置に中心２６５を有する６番目領域２６０が設定される。次に、図１５に示すように、ｎを５として、画像２００の一辺２０１に対して垂直方向において、６番目領域２６０の中心２６５から一定の間隔Ｇ離れた位置に中心２７５を有する７番目領域２７０が設定される。次に、図１６に示すように、ｎを６として、画像２００の一辺２０１に対して垂直方向において、７番目領域２７０の中心２７５から一定の間隔Ｇ離れた位置に中心２８５を有する８番目領域２８０が設定される。次に、図１７に示すように、ｎを７として、画像２００の一辺２０１に対して垂直方向において、８番目領域２８０の中心２８５から一定の間隔Ｇ離れた位置に中心２９５を有する９番目領域２９０が設定される。

【0062】

図１８に示すように、ｎを８として、画像２００の一辺２０１に対して垂直方向において、９番目領域２９０の中心２９５から一定の間隔Ｇ離れた位置に中心４０５を有する１０番目領域４００を設定しようとすると、１０番目領域は、画像２００の一辺２０１に対向する他辺２０３を超えるため、繰り返し指示部３０４は、第２設定部３０３に、領域の設定を停止させる。なお、９個の領域を設定するのは単なる例示であり、領域の数はこれに限定されない。また、繰り返し指示部３０４は、ｎ＋２番目領域の中心が画像２００の一辺２０１に対向する他辺２０３を超えたときに、第２設定部３０３に、領域の設定を停止させてもよいし、ｎ＋２番目領域の幅が確定した後に、ｎ＋２番目領域の少なくとも一部が画像２００の一辺２０１に対向する他辺２０３を超えたときに、第２設定部３０３に、領域の設定を停止させてもよい。

【0063】

図１に示す繰り返し指示部３０４は、例えば、ｎ＋２番目領域が画像２００の一辺２０１に対向する他辺２０３を超える直前のｎの最大値ｑをメモリー３５３に保存する。特徴量取得部３０５は、メモリー３５３からｎの最大値ｑを読み出し、ｎ＋１番目領域（ｎ：０からｑ）に含まれる複数の輪郭の特徴量を取得する。輪郭の特徴量は、例えば、輪郭の統計値であり、輪郭の径、輪郭の平均径、輪郭の周囲長、輪郭の平均周囲長、輪郭の円形度、輪郭の真円度、及び輪郭の歪みであってもよい。輪郭の特徴量は、第１設定部３０２及び第２設定部３０３が、ｎ＋１番目領域（ｎ：０からｑ）を設定するときに算出してもよいし、特徴量取得部３０５が算出してもよい。特徴量取得部３０５は、例えば、特徴量取得部３０５は、ｎ＋１番目領域（ｎ：０からｑ）に含まれる輪郭の特徴量のデータを生成してもよい。例えば、１番目領域に含まれる輪郭の特徴量、２番目領域に含まれる輪郭の特徴量、・・・及びｑ番目領域に含まれる輪郭の特徴量を含むデータを生成する。特徴量取得部３０５は、データを表にしてもよい。あるいは、特徴量取得部３０は、データを、図１９に例示するにグラフにプロットしてもよい。特徴量取得部３０５は、輪郭の特徴量を、図２０に示すように、ｎ＋１番目領域（ｎ：０からｑ）の中心の位置に対してプロットしてもよい。これにより、１番目領域からｑ番目領域までの輪郭の特徴量の推移を把握することが可能である。

【0064】

次に、図２１を参照して、実施形態に係る画像に含まれる複数の輪郭の分析方法を説明する。

【0065】

ステップＳ１０１で、画像取得部３０１が、複数の輪郭を含む画像を取得する。ステップＳ１０２で、第１設定部３０２が、画像の一辺を含む１番目領域を特定する。ステップＳ１０３で、第１設定部３０２が、１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出し、１番目領域の中心を設定する。

【0066】

ステップＳ１０３においては、図２２に示すように、ステップＳ２００で、第１設定部３０２が、１番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出する。ステップＳ２０１で、第１設定部３０２が、１番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にあるか否かを判定する。１番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、ステップＳ２０２で、第１設定部３０２が、一辺に対して垂直方向における１番目領域の幅の中央を、１番目領域の中心と設定する。

【0067】

１番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にない場合、ステップＳ２０３で、第１設定部３０２が、径の変動係数が所定の範囲より大きいか否かを判定する。径の変動係数が所定の範囲より大きい場合、ステップＳ２０４で、第１設定部３０２が、１番目領域の幅を狭くする。その後、ステップＳ２００に戻る。径の変動係数が所定の範囲内になるまで、ステップＳ２００、ステップＳ２０１、ステップＳ２０３、及びステップＳ２０４からなるループが繰り返される。

【0068】

ステップＳ２０３で、第１設定部３０２が、径の変動係数が所定の範囲より小さいと判定した場合、ステップＳ２０５で、第１設定部３０２が、１番目領域の幅を広げる。その後、ステップＳ２００に戻る。径の変動係数が所定の範囲内になるまで、ステップＳ２００、ステップＳ２０１、ステップＳ２０３、及びステップＳ２０５からなるループが繰り返される。

【0069】

図２１のステップＳ１０４で、第２設定部３０３が、ｎを０に設定する。ステップＳ１０５で、第２設定部３０３が、ｎ＋１番目領域の中心から一定の間隔離れたｎ＋２番目領域の中心を設定する。ステップＳ１０６で、第２設定部３０３が、ｎ＋２番目領域の幅を、ｎ＋１番目領域に含まれる複数の輪郭の最大径以上に設定する。ステップＳ１０７で、第２設定部３０３が、ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数に基づき、ｎ＋２番目領域の幅を確定する。

【0070】

ステップＳ１０７においては、図２３に示すように、ステップＳ３００で、第２設定部３０３が、ｎ＋２番目領域に含まれる複数の輪郭の最大径及び径の変動係数を算出する。ステップＳ３０１で、第２設定部３０３が、ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にあるか否かを判定する。ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にある場合、ステップＳ３０２で、第２設定部３０３が、現在のｎ＋２番目領域の幅を採用する。

【0071】

ｎ＋２番目領域に含まれる複数の輪郭の径の変動係数が所定の範囲内にない場合、ステップＳ３０３で、第２設定部３０３が、径の変動係数が所定の範囲より大きいか否かを判定する。径の変動係数が所定の範囲より大きい場合、ステップＳ３０４で、第２設定部３０３が、ｎ＋２番目領域の中心を保ったまま、ｎ＋２番目領域の幅を狭くする。その後、ステップＳ３００に戻る。径の変動係数が所定の範囲内になるまで、ステップＳ３００、ステップＳ３０１、ステップＳ３０３、及びステップＳ３０４からなるループが繰り返される。

【0072】

ステップＳ３０３で、第２設定部３０３が、径の変動係数が所定の範囲より小さいと判定した場合、ステップＳ３０５で、第２設定部３０３が、ｎ＋２番目領域の中心を保ったまま、ｎ＋２番目領域の幅を広げる。その後、ステップＳ３００に戻る。径の変動係数が所定の範囲内になるまで、ステップＳ３００、ステップＳ３０１、ステップＳ３０３、及びステップＳ３０５からなるループが繰り返される。

【0073】

図２１のステップＳ１０８で、第２設定部３０３が、ｎ＋２番目領域が、画像の１番目領域が設定された一辺に対向する他辺を超えているか否かを判定する。ｎ＋２番目領域が、画像の１番目領域が接する一辺に対向する他辺を超えていない場合、ステップＳ１１０で、繰り返し指示部３０４が、ｎに１を加算し、第２設定部３０３に、ステップＳ１０５からステップＳ１０８を繰り返させる。ステップＳ１０８で、ｎ＋２番目領域が、画像の１番目領域が接する一辺に対向する他辺を超えている場合、ステップＳ１０９で、特徴量取得部３０５が、ｎ＋１番目領域に含まれる複数の輪郭の特徴量を取得する。

【0074】

実施形態に係る画像に含まれる複数の輪郭の分析システム及び方法によれば、画像を複数の領域に分ける際に、複数の領域のそれぞれの幅を、領域に含まれる輪郭の変動係数が所定の範囲になるように設定するため、画像の一辺に対して垂直方向における、輪郭の特徴量の変化を正確に評価することが可能である。例えば、画像の被写体が中空糸多孔質膜である場合、中空糸多孔膜の断面における内表面側から外表面側への平均孔径の変化を正確に評価することが可能である。

【符号の説明】

【0075】

２００・・・画像、２０１・・・一辺、２０３・・・他辺、２１０、２１１、２２０、２３０、２４０、２５０、２６０、２７０、２８０、２９０、４００・・・領域、２１５、２２５、２３５、２４５、２５５、２６５、２７５、２８５、２９５、４０５・・・中心、３００・・・ＣＰＵ、３０１・・・画像取得部、３０２・・・第１設定部、３０３・・・第２設定部、３０４・・・繰り返し指示部、３０５・・・特徴量取得部、３０６・・・多値化部、３５１・・・入力装置、３５２・・・出力装置、３５３・・・メモリー、５００、５０２、５０５・・・輪郭、５０１・・・円

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】

【図20】

【図21】

【図22】

【図23】