(57)【要約】【課題】ＣＭＣの３次元画像から繊維束の配向を容易に解析し得る画像解析装置、画像解析方法及びプログラムを提案する。【解決手段】Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像からＸ糸及びＹ糸の繊維束の配向を解析する画像解析装置において、３次元画像を２値化する２値化部と、２値化された画像からＸ糸及びＹ糸が垂直に立体交差する重なり領域を抽出する重なり領域抽出部と、重なり領域に含まれる各ボクセルの重なり方向を平均化し、平均化した方向をリファレンス方向として決定するリファレンス方向決定部と、リファレンス方向に対して垂直なリファレンス平面上で方向距離法を適用し、２値化された画像からＺ糸を除去するＺ糸除去部と、Ｚ糸が除去された画像に対してリファレンス平面上で方向距離法を再度適用し、適用時に計算される方向距離に基づいて、Ｘ糸及びＹ糸の繊維束の配向を推定する繊維束配向推定部とを備えることを特徴とする。
2016052489.tif 000002
【特許請求の範囲】
【請求項１】
Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像から前記Ｘ糸及び前記Ｙ糸の繊維束の配向を解析する画像解析装置において、
前記３次元画像を２値化する２値化部と、
前記２値化された画像から前記Ｘ糸及び前記Ｙ糸が垂直に立体交差する重なり領域を抽出する重なり領域抽出部と、
前記重なり領域に含まれる各ボクセルの重なり方向を平均化し、平均化した方向をリファレンス方向として決定するリファレンス方向決定部と、
前記リファレンス方向に対して垂直なリファレンス平面上で方向距離法を適用し、前記２値化された画像から前記Ｚ糸を除去するＺ糸除去部と、
前記Ｚ糸が除去された画像に対して前記リファレンス平面上で方向距離法を再度適用し、適用時に計算される方向距離に基づいて、前記Ｘ糸及び前記Ｙ糸の繊維束の配向を推定する繊維束配向推定部と
を備えることを特徴とする画像解析装置。
【請求項２】
前記２値化部は、
前記３次元画像を所定の閾値に基づいて２値化し、
前記２値化後の画像に対してクロージング処理を実行する
ことを特徴とする請求項１に記載の画像解析装置。
【請求項３】
前記重なり領域抽出部は、
前記クロージング処理後の画像に対してオープニング処理を実行し、
前記オープニング処理後の画像に対して膨張処理を実行し、
前記膨張処理後の画像と、前記クロージング処理後の画像との積集合を計算して、前記重なり領域を抽出する
ことを特徴とする請求項２に記載の画像解析装置。
【請求項４】
前記リファレンス方向決定部は、
前記重なり領域に対して方向距離法を適用し、適用時に計算される各ボクセルの方向距離に基づいて、前記重なり領域に含まれる各ボクセルの重なり方向を推定し、
前記重なり領域に対して収縮処理を実行して中心領域を抽出し、
前記中心領域に含まれる各ボクセルの前記重なり方向を平均化し、平均化した方向を前記リファレンス方向として決定する
ことを特徴とする請求項１〜３の何れか一項に記載の画像解析装置。
【請求項５】
前記Ｚ糸除去部は、
前記リファレンス平面上で方向距離法を適用し、適用時に計算される各ボクセルの方向距離に基づいて、方向テンソルの２番目に大きい固有値に対する固有ベクトルを計算し、
前記固有ベクトルに基づいて、前記２値化された画像を前記Ｘ糸又は前記Ｙ糸の領域と、前記Ｚ糸の領域とに分離し、
前記Ｚ糸の領域に含まれるボクセルを除去することにより前記２値化された画像から前記Ｚ糸を除去する
ことを特徴とする請求項１〜４の何れか一項に記載の画像解析装置。
【請求項６】
前記Ｚ糸除去部は、
注目ボクセルの方向距離を計算する際、前記注目ボクセルの固有ベクトルが示す方向及び前記固有ベクトルが示す方向とは逆方向にボクセルを進行し、背景に到達した場合又は次に進行するボクセルの固有ベクトルが示す方向とのなす角が所定の閾値よりも大きい場合に進行を停止し、最後に到達したボクセル間の距離を条件付き方向距離として計算し、
前記条件付き方向距離に基づいて、前記２値化された画像を前記Ｘ糸又は前記Ｙ糸の領域と、前記Ｚ糸の領域とに分離する
ことを特徴とする請求項５に記載の画像解析装置。
【請求項７】
前記Ｚ糸除去部は、
前記２値化された画像から前記Ｚ糸を除去した後、膨張処理を実行し、
前記膨張処理後の画像と、前記膨張処理前の画像との積集合を計算する
ことを特徴とする請求項１〜６の何れか一項に記載の画像解析装置。
【請求項８】
前記繊維束配向推定部は、
前記Ｚ糸が除去された画像に対して前記リファレンス平面上で方向距離法を再度適用し、前記Ｘ糸及び前記Ｙ糸の繊維束の配向を推定し、
前記繊維束の配向を推定したボクセルのうち、中心部近傍のボクセルを抽出し、
前記中心部近傍のボクセルについてクラスタリング処理を実行し、
所属するボクセルが最小のクラスタを削除してノイズを除去する
ことを特徴とする請求項１〜７の何れか一項に記載の画像解析装置。
【請求項９】
前記Ｘ糸及び前記Ｙ糸の繊維束の配向を示す画像を表示する表示部を備える
ことを特徴とする請求項１〜８の何れか一項に記載の画像解析装置。
【請求項１０】
Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像から前記Ｘ糸及び前記Ｙ糸の繊維束の配向を解析する画像解析方法において、
コンピュータが、
前記３次元画像を２値化する第１のステップと、
前記２値化された画像から前記Ｘ糸及び前記Ｙ糸が垂直に立体交差する重なり領域を抽出する第２のステップと、
前記重なり領域に含まれる各ボクセルの重なり方向を平均化し、平均化した方向をリファレンス方向として決定する第３のステップと、
前記リファレンス方向に対して垂直なリファレンス平面上で方向距離法を適用し、前記２値化された画像から前記Ｚ糸を除去する第４のステップと、
前記Ｚ糸が除去された画像に対して前記リファレンス平面上で方向距離法を再度適用し、適用時に計算される方向距離に基づいて、前記Ｘ糸及び前記Ｙ糸の繊維束の配向を推定する第５のステップと
を備えることを特徴とする画像解析方法。
【請求項１１】
Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像から前記Ｘ糸及び前記Ｙ糸の繊維束の配向を解析するプログラムにおいて、
コンピュータに、
前記３次元画像を２値化する第１のステップと、
前記２値化された画像から前記Ｘ糸及び前記Ｙ糸が垂直に立体交差する重なり領域を抽出する第２のステップと、
前記重なり領域に含まれる各ボクセルの重なり方向を平均化し、平均化した方向をリファレンス方向として決定する第３のステップと、
前記リファレンス方向に対して垂直なリファレンス平面上で方向距離法を適用し、前記２値化された画像から前記Ｚ糸を除去する第４のステップと、
前記Ｚ糸が除去された画像に対して前記リファレンス平面上で方向距離法を再度適用し、適用時に計算される方向距離に基づいて、前記Ｘ糸及び前記Ｙ糸の繊維束の配向を推定する第５のステップと
を実行させるためのプログラム。
【請求項１２】
Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像を２値化する２値化部と、
前記２値化された画像から前記Ｘ糸及び前記Ｙ糸が立体交差する重なり領域を抽出する重なり領域抽出部と、
前記抽出された重なり領域の形態を解析する重なり領域形態解析部と
を備えることを特徴とする画像解析装置。
【請求項１３】
前記重なり領域形態解析部は、更に、前記重なり領域の形態の解析の結果に基づいて、前記Ｘ糸又は前記Ｙ糸の繊維束の配向異常を検出する
ことを特徴とする請求項１２に記載の画像解析装置。
【請求項１４】
前記重なり領域形態解析部は、前記重なり領域の体積を算出する
ことを特徴とする請求項１２又は１３に記載の画像解析装置。
【請求項１５】
前記重なり領域形態解析部は、前記重なり領域が伸長する方向を算出する
ことを特徴とする請求項１２乃至１４のいずれか一項に記載の画像解析装置。
【請求項１６】
前記重なり領域形態解析部は、複数の前記重なり領域の重心位置を算出する
ことを特徴とする請求項１２乃至１５のいずれか一項に記載の画像解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像解析装置、画像解析方法及びプログラムに関し、特に繊維強化複合材料に含まれる繊維束の配向を解析する画像解析装置、画像解析方法及びプログラムに適用して好適なものである。
【背景技術】
【０００２】
近年、繊維強化複合材料の一種であるセラミック基複合材料（CMC : Ceramic Matrix Composites）の開発が進められている。ＣＭＣは、セラミックの繊維を母材（マトリクス）で強化した複合材料であり、軽量かつ耐熱性に優れるという特徴を有する。この特徴を利用して、例えばＣＭＣを航空機用エンジン部材に用いることが検討されており、現在実用化が進められている。ＣＭＣを航空機用エンジン部材として用いることで、大幅な燃費向上を期待することができる。【０００３】
ＣＭＣの一般的な形成プロセスは、次の通りである。まず数百本程度のセラミックの繊維を束ねて繊維束を作製し、この繊維束を編み込んで織物を作製する。繊維束の編み込み手法には、例えば３次元織り又は平織りがある。３次元織りは、ＸＹＺ方向の３方向から繊維束を編み込んで織物を作製する手法であり、平織りは、ＸＹ方向の２方向から繊維束を編み込んで織物を作製する手法である。【０００４】
織物を作製した後はＣＶＩ（Chemical Vapor Infiltration）及びＰＩＰ（Polymer Impregnation and Pyrolysis）によりマトリクスを形成し、最後に機械加工及び表面コーティング等を行うことにより、ＣＭＣが形成される。ここで、このとき形成されたＣＭＣにおける繊維束の配向がＣＭＣの強度に大きく影響する。【０００５】
すなわち繊維束が本来直線であるべき箇所で蛇行していたり、本来配置されるべき基準軸から全体的にずれて偏向していたり、或いは途中で断裂していたりする場合、ＣＭＣの強度が低下する。これに対し繊維束が蛇行、偏向又は断裂せずに一定の方向に向きをそろえて適切に配列している場合、強度が高く耐熱性に優れる。よって形成されたＣＭＣの強度が十分であるか否かを確認するためには、繊維束の配向を評価することが重要となる。【０００６】
特許文献１には、樹脂成形品のスライス画像を２値化して２値画像を取得し、この２値画像をフーリエ変換してパワースペクトル画像を取得し、このパワースペクトル画像により描かれる楕円と直交する楕円の主軸方向を樹脂成形品に含まれる充填剤（繊維）の配向方向とする配向解析方法が開示されている。【０００７】
また非特許文献１には、繊維束を編み込んだ織物をＸ線ＣＴ装置により撮像してＸ線ＣＴ画像を取得し、このＸ線ＣＴ画像に対して特殊なフィルタ関数を用いて計算することにより、繊維束を構成する繊維の一本一本の配向を解析する技術が開示されている。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２０１２‐２５４７号公報
【非特許文献】
【０００９】
【非特許文献１】T. Shinohara, J. Takayama, S. Ohyama, and A. Kobayashi, “Extractionof Yarn Positional Information from a Three‐dimensional CT Image of Textile Fabricusing Yarn Tracingwith a Filament Model for Structure Analysis”,Textile Research Journal , Vol. 80, No. 7, pp.623-630 (2010)
【発明の概要】
【発明が解決しようとする課題】
【００１０】
しかし特許文献１に記載の技術では、スライス画像に含まれる充填剤（繊維）の配向について、一方向しか解析結果として得ることができない。従って例えば３次元織りや平織りのように繊維束が複数の方向に配列している場合、それぞれの繊維束の配向を解析結果として得ることはできない。また繊維束が蛇行、偏向又は断裂せずに一定の方向に向きをそろえて適切に配列しているか否かを解析することはできない。【００１１】
また非特許文献１に記載の技術では、繊維束を構成する繊維の一本一本が識別できる程度の高解像度のＸ線ＣＴ画像を取得する必要がある。この場合、Ｘ線ＣＴ画像を得るための撮像時間が長くなり、製品検査に利用することができないため実用的ではない。また断面が円形である繊維に対して有効な技術であり、断面が扁平である繊維束の配向を解析する技術としてそのまま利用することはできない。更にはＸ線ＣＴ画像内の各繊維の開始点を入力する必要があり操作が煩雑になるという問題がある。【００１２】
本発明の開示は以上の点を考慮してなされたもので、ＣＭＣの３次元画像から繊維束の配向を容易に解析し得る画像解析装置、画像解析方法及びプログラムを提案するものである。
【課題を解決するための手段】
【００１３】
かかる課題を解決するために、本発明の開示においては、Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像からＸ糸及びＹ糸の繊維束の配向を解析する画像解析装置において、３次元画像を２値化する２値化部と、２値化された画像からＸ糸及びＹ糸が垂直に立体交差する重なり領域を抽出する重なり領域抽出部と、重なり領域に含まれる各ボクセルの重なり方向を平均化し、平均化した方向をリファレンス方向として決定するリファレンス方向決定部と、リファレンス方向に対して垂直なリファレンス平面上で方向距離法を適用し、２値化された画像からＺ糸を除去するＺ糸除去部と、Ｚ糸が除去された画像に対してリファレンス平面上で方向距離法を再度適用し、適用時に計算される方向距離に基づいて、Ｘ糸及びＹ糸の繊維束の配向を推定する繊維束配向推定部とを備えることを特徴とする。【００１４】
またかかる課題を解決するために、本発明の開示においては、Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像からＸ糸及びＹ糸の繊維束の配向を解析する画像解析方法において、コンピュータが、３次元画像を２値化する第１のステップと、２値化された画像からＸ糸及びＹ糸が垂直に立体交差する重なり領域を抽出する第２のステップと、重なり領域に含まれる各ボクセルの重なり方向を平均化し、平均化した方向をリファレンス方向として決定する第３のステップと、リファレンス方向に対して垂直なリファレンス平面上で方向距離法を適用し、２値化された画像からＺ糸を除去する第４のステップと、Ｚ糸が除去された画像に対してリファレンス平面上で方向距離法を再度適用し、適用時に計算される方向距離に基づいて、Ｘ糸及びＹ糸の繊維束の配向を推定する第５のステップとを備えることを特徴とする。【００１５】
またかかる課題を解決するために、本発明の開示においては、Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像からＸ糸及びＹ糸の繊維束の配向を解析するプログラムにおいて、コンピュータに、３次元画像を２値化する第１のステップと、２値化された画像からＸ糸及びＹ糸が垂直に立体交差する重なり領域を抽出する第２のステップと、重なり領域に含まれる各ボクセルの重なり方向を平均化し、平均化した方向をリファレンス方向として決定する第３のステップと、リファレンス方向に対して垂直なリファレンス平面上で方向距離法を適用し、２値化された画像からＺ糸を除去する第４のステップと、Ｚ糸が除去された画像に対してリファレンス平面上で方向距離法を再度適用し、適用時に計算される方向距離に基づいて、Ｘ糸及びＹ糸の繊維束の配向を推定する第５のステップとを実行させることを特徴とする。【００１６】
またかかる課題を解決するために、本開示の画像解析装置は、Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像を２値化する２値化部と、前記２値化された画像から前記Ｘ糸及び前記Ｙ糸が立体交差する重なり領域を抽出する重なり領域抽出部と、前記抽出された重なり領域の形態を解析する重なり領域形態解析部とを備えることを特徴とする。【００１７】
またかかる課題を解決するために、本開示の画像解析方法は、Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像を２値化するステップと、前記２値化された画像から前記Ｘ糸及び前記Ｙ糸が立体交差する重なり領域を抽出するステップと、前記抽出された重なり領域の形態を解析するステップと、を備えることを特徴とする。【００１８】
またかかる課題を解決するために、本開示のプログラムは、コンピュータに、Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像を２値化するステップと、前記２値化された画像から前記Ｘ糸及び前記Ｙ糸が立体交差する重なり領域を抽出するステップと、前記抽出された重なり領域の形態を解析するステップとを実行させるためのプログラムである。
【発明の効果】
【００１９】
本発明の開示によれば、ＣＭＣの３次元画像から繊維束の配向を容易に解析することができる。
【図面の簡単な説明】
【００２０】
【図１】画像解析装置の全体構成図である。
【図２】３次元織りの織物の概念図である。
【図３】画像解析処理の全体フローチャートである。
【図４】２値化処理の詳細フローチャートである。
【図５】２値化処理により作成される処理画像である。
【図６】重なり領域抽出処理の詳細フローチャートである。
【図７】重なり領域抽出処理により作成される処理画像である。
【図８】リファレンス方向決定処理の詳細フローチャートである。
【図９】通常の方向距離法を適用して方向距離を推定する処理の概要図である。
【図１０】リファレンス方向決定処理により作成される処理画像である。
【図１１】Ｚ糸除去処理の詳細フローチャートである。
【図１２】通常の方向距離法を適用した場合の概念図である。
【図１３】リファレンス付き方向距離法の概念図である。
【図１４】リファレンス付き方向距離法を適用した場合の概念図である。
【図１５】条件付き方向距離を推定する処理の概要図である。
【図１６】Ｘ糸又はＹ糸の領域とＺ糸の領域とを示す概念図である。
【図１７】繊維束配向推定処理の詳細フローチャートである。
【図１８】入力画像が模擬画像である場合の解析結果である。
【図１９】入力画像が高解像度のＸ線ＣＴ画像である場合の解析結果である。
【図２０】入力画像が低解像度のＸ線ＣＴ画像である場合の解析結果である。
【図２１】入力画像が欠陥模擬画像である場合の解析結果である。
【図２２】入力画像が曲面模擬画像である場合の解析結果である。
【図２３】画像解析装置の全体構成図である。
【図２４】画像解析処理の全体フローチャートである。
【図２５】ピラーの体積の算出の例について示す図である。
【図２６】ピラーが伸長する方向の算出の例について示す図である。
【図２７】ピラー折れ曲がっている場合の画像処理について示す図である。
【図２８】ピラーの一部が不規則に並んでいる例について示す図である。
【図２９】Ｘ糸及びＹ糸の繊維束の理想的な配置についてＸＹ平面で示す図である。
【図３０】図２９の繊維束の配置におけるＸＺ平面（又はＹＺ平面）における重なり領域抽出画像を模式的に表す図である。
【図３１】図２９の繊維束の配置におけるＸＹ平面における重なり領域抽出画像を模式的に表す図である。
【図３２】ＸＹ平面において、繊維束の一部が延伸すべき方向とは別のＸＹ平面内の方向に伸びている状態について示す図である。
【図３３】図３２の繊維束の状態における、ＹＺ平面の重なり領域抽出画像を模式的に示す図である。
【図３４】図３２の繊維束の状態における、ＸＹ平面の重なり領域抽出画像を模式的に示す図である。
【図３５】ＸＹ平面において、繊維束の一部が延伸すべき方向とは別のＺ方向に伸びている状態について示す図である。
【図３６】図３５の繊維束の状態における、ＹＺ平面の重なり領域抽出画像を模式的に示す図である。
【図３７】ＸＹ平面において、繊維束の一部が延伸すべき方向とは別の、ＸＹ平面内の方向成分及びＺ方向成分を共に有する方向に伸びている状態について示す図である。
【図３８】図３７の繊維束の状態における、ＹＺ平面の重なり領域抽出画像を模式的に示す図である。
【図３９】ＸＹ平面において、それぞれ繊維束が全体として傾いた方向に伸びている状態について示す図である。
【図４０】ＸＹ平面において、それぞれ繊維束が全体として傾いた方向に伸びている状態について示す図である。
【図４１】図３９又は４０に示される繊維束の状態における、ＸＹ平面の重なり領域抽出画像を模式的に示す図である。
【図４２】ＸＹ平面において、繊維束の一部（例えばＸ糸）がない状態について示す図である。
【図４３】ＸＹ平面において、繊維束の一部の太さが他の繊維束より太く又は細くなっている場合について示す図である。
【図４４】図４２又は図４３の繊維束の状態における、ＹＺ平面の重なり領域抽出画像を模式的に示す図である。
【図４５】ＸＹ平面において、繊維束の一部（例えばＸ糸）が折り返されている状態について示す図である。
【図４６】図４５の繊維束の状態における、ＸＹ平面の重なり領域抽出画像を模式的に示す図である。
【図４７】画像解析装置の全体構成図である。
【図４８】画像解析処理の全体フローチャートである。
【発明を実施するための形態】
【００２１】
以下図面について、本発明の一実施の形態を詳述する。
【００２２】
（１）画像解析装置１の全体構成
図１は、本実施の形態における画像解析装置１の全体構成を示す。画像解析装置１は、ＣＰＵ（Central Processing Unit）１１、入力部１２、記憶部１３、表示部１４及びメモリ１５を備えて構成されるコンピュータである。
【００２３】
ＣＰＵ１１は、メモリ１５に格納されている各種プログラムと協働して、画像解析装置１の動作を統括的に制御するプロセッサである。入力部１２は、ユーザからの入力を受け付けるインタフェースであり、例えばキーボードやマウスである。本実施の形態における入力部１２は、ＣＭＣ（Ceramic Matrix Composites）を構成する織物のＸ線ＣＴ画像Ｇ１０を入力するインタフェースでもある。
【００２４】
ここでＣＭＣとは、数百本程度のセラミックの繊維を束ねて繊維束を作製し、この繊維束を編み込んで織物を作製した後、繊維表面をカーボン等で被覆し、マトリクスを形成するためのＣＶＩ（Chemical Vapor Infiltration）プロセス及びＰＩＰ（Polymer Impregnation and Pyrolysis）プロセス等を経て形成される繊維強化複合材料である。
【００２５】
織物を作製する際の編み込み手法には、３次元織り又は平織りと呼ばれるものがある。３次元織りは、ＸＹＺ方向の３方向から繊維束を編み込んで織物を作製する手法であり、平織りは、ＸＹ方向の２方向から繊維束を編み込んで織物を作製する手法である。
【００２６】
図２は、３次元織りの織物の概念図を示す。本実施の形態においては、特に複雑な３次元織りにより作製された織物又はこの織物から形成されるＣＭＣをＸ線ＣＴ装置で撮像し、得られたＸ線ＣＴ画像Ｇ１０から繊維束の配向を自動解析しようとするものである。
【００２７】
図２に示すように３次元織りの織物は、等間隔に配列したＸ糸の繊維束からなる繊維層と、等間隔に配列したＹ糸の繊維束からなる繊維層とを交互に積層し、積層した多数の繊維層がＣＭＣの形成過程で剥離しないようにＺ糸の繊維束により固定して形成される。
【００２８】
この織物により形成されるＣＭＣは、通常Ｘ糸方向又はＹ糸方向に伸縮することを想定して設計される。よってＸ糸及びＹ糸に略垂直に編み込まれるＺ糸は、ＣＭＣの強度に直接的に影響しない。一方でＺ糸の存在は、Ｘ糸及びＹ糸の配向を解析する際、精度不良を及ぼす場合がある。
【００２９】
そこで本実施の形態においては、３次元織りの織物のＸ線ＣＴ画像Ｇ１０からＺ糸を除去し、Ｘ糸及びＹ糸の繊維束の配向を精度良く解析しようとするものである。
【００３０】
なお一般に配向とは、一定の方向に向きをそろえて配列すること又は配列した状態を意味する用語であり、本実施の形態においても同様の意味で用いる。蛇行、偏向又は断裂して配列した状態であっても、一定の方向に向きをそろえて配列した状態であれば配向と呼ぶものとする。
【００３１】
図１に戻り、記憶部１３は、入力部１２から入力されたＸ線ＣＴ画像Ｇ１０及びこのＸ線ＣＴ画像Ｇ１０に各種画像処理を施した処理画像を記憶する記憶媒体である。表示部１４は、Ｘ線ＣＴ画像Ｇ１０画像及び処理画像を表示するＬＣＤ（Liquid Crystal Display）等の表示装置である。例えば表示部１４は、Ｘ線ＣＴ画像Ｇ１０から繊維束の配向を自動解析し、繊維束配向推定画像Ｇ１００を表示画面上に表示する。
【００３２】
メモリ１５は、ＣＰＵ１１と協働して画像解析処理を実行するための各種プログラムを格納する記憶媒体である。各種プログラムには、２値化部１５１、重なり領域抽出部１５２、リファレンス方向決定部１５３、Ｚ糸除去部１５４及び繊維束配向推定部１５５がある。これら各種プログラムにより実行される画像解析処理（図３）については後述する。
【００３３】
（２）画像解析処理のフローチャート
図３は、本実施の形態における画像解析処理Ｐ１の全体フローチャートを示す。この画像解析処理は、ユーザからの実行指示を入力部１２が受け付けたことを契機として、ＣＰＵ１１とメモリ１５に格納されている各種プログラムとの協働により実行される。以下説明の便宜上、処理主体を各種プログラムとして説明する。
【００３４】
まず２値化部１５１は、入力部１２を介してＸ線ＣＴ画像Ｇ１０を入力すると（ＳＰ１）、入力したＸ線ＣＴ画像Ｇ１０を所定の閾値に基づいて２値化し、前景にＸ糸、Ｙ糸及びＺ糸の各繊維束が示される２値画像を作成する（ＳＰ２）。
【００３５】
次いで重なり領域抽出部１５２は、Ｘ糸とＹ糸とが垂直に立体交差する重なり領域を抽出し（ＳＰ３）、リファレンス方向決定部１５３は、抽出された重なり領域の重なり方向をリファレンス方向として決定する（ＳＰ４）。
【００３６】
ここで重なり領域を抽出する理由は、抽出した重なり領域に対して通常の方向距離法を適用し、重なり方向を推定するためである。通常の方向距離法については後述する。
【００３７】
また重なり方向をリファレンス方向として決定する理由は、このリファレンス方向に垂直な平面上で２次元の通常の方向距離法を適用し、繊維束の配向を推定するためである。
【００３８】
リファレンス方向に垂直な平面をリファレンス平面と呼び、リファレンス平面上で２次元の通常の方向距離法を適用する手法をリファレンス付き方向距離法と呼ぶ。
【００３９】
リファレンス平面上にはＸ糸又はＹ糸が存在するため、リファレンス付き方向距離法を適用することにより、Ｘ糸及びＹ糸の繊維束の配向を精度良く推定することができる。リファレンス付き方向距離法については後述する。
【００４０】
次いでＺ糸除去部１５４は、２値画像に対してリファレンス付き方向距離法を適用し、Ｘ糸、Ｙ糸及びＺ糸の繊維束の配向を推定する。そしてＺ糸除去部１５４は、推定した配向に基づいて２値画像に含まれるＺ糸を除去する（ＳＰ５）。
【００４１】
次いで繊維束配向推定部１５５は、Ｚ糸を除去した２値画像に対して再びリファレンス付き方向距離法を適用し、Ｘ糸及びＹ糸の繊維束の配向を推定する（ＳＰ６）。
【００４２】
そして繊維束配向推定部１５５は、繊維束配向推定画像Ｇ１００を作成し、表示部１４により表示して（ＳＰ７）、この画像解析処理を終了する。
【００４３】
（３）各処理の詳細
以下図４〜図１７を参照して、また計算式（式１〜式３）を用いて、図３で説明した各処理（ＳＰ２〜ＳＰ６）の詳細について説明する。なおＸ線ＣＴ画像Ｇ１０の入力処理（ＳＰ１）及び繊維束配向推定画像Ｇ１００の表示処理（ＳＰ７）については、一般的な手法によるものであるため説明を省略する。
【００４４】
（３−１）２値化処理
図４は、２値化処理の詳細なフローチャートを示す。２値化部１５１は、まず入力部１２を介して入力したＸ線ＣＴ画像Ｇ１０を所定の閾値に基づいて２値化し、２値画像を作成する（ＳＰ２１）。作成された２値画像は、繊維束内の繊維の密度にムラがあることを原因とした黒抜けが生じる場合がある。すなわち本来背景であるはずのボクセルが前景ボクセルとされた個所が存在する。
【００４５】
この本来背景であるはずの前景ボクセルを背景ボクセルに戻すため、２値化部１５１はモルフォロジー処理のうちのクロージング処理を実行する（ＳＰ２２）。クロージング処理を実行することにより、２値画像に生じた黒抜けを除去することができる。２値化部１５１は、黒抜けを除去してクリーニングした２値画像を作成して、本処理を終了する。
【００４６】
図５は、２値化処理により作成される処理画像を示す。具体的にはクロージング処理前の２値画像Ｇ２１と、クロージング処理後の２値画像Ｇ２２とを示す。
【００４７】
図５に示すように２値画像Ｇ２１には本来背景であるはずのボクセルが前景ボクセルとされた個所がある。クロージング処理を実行することで、黒抜けを除去したクロージング処理後の２値画像Ｇ２２を得ることができる。
【００４８】
（３−２）重なり領域抽出処理
図６は、重なり領域抽出処理の詳細なフローチャートを示す。重なり領域抽出部１５２は、Ｘ糸とＹ糸とが立体交差する重なり領域を抽出するため、モルフォロジー処理のうちのオープニング処理を実行する（ＳＰ３１）。
【００４９】
オープニング処理を実行することにより、大まかな重なり領域を抽出することができる。オープニング処理後の２値画像の形状は変化しており、本来重なり領域であるはずの前景ボクセルが重なり領域ではない背景ボクセルとされた個所が存在する。
【００５０】
この本来重なり領域である背景ボクセルを前景ボクセルに戻すため、重なり領域抽出部１５２は、モルフォロジー処理のうちの膨張処理を実行する（ＳＰ３２）。
【００５１】
次いで重なり領域抽出部１５２は、膨張処理後の２値画像とオープニング処理前の２値画像との積集合を計算して、正確な重なり領域を抽出する（ＳＰ３３）。重なり領域抽出部１５２は、重なり領域を抽出した重なり領域抽出画像を作成して、本処理を終了する。
【００５２】
図７は、重なり領域抽出処理により作成される処理画像を示す。具体的にはオープニング処理前の２値画像Ｇ３、オープニング処理後の２値画像Ｇ３１、膨張処理後の２値画像Ｇ３２及び重なり領域抽出画像Ｇ３３を示す。
【００５３】
図７に示すようにオープニング処理前の２値画像Ｇ３と、膨張処理後の２値画像Ｇ３２との積集合を計算することで、正確な重なり領域が抽出された重なり領域抽出画像Ｇ３３を得ることができる。
【００５４】
（３−３）リファレンス方向決定処理
図８は、リファレンス方向決定処理の詳細なフローチャートを示す。リファレンス方向決定部１５３は、重なり領域抽出画像Ｇ３３（図７）に対して通常の方向距離法を適用し、重なり領域における各ボクセルについて重なり方向を推定する（ＳＰ４１）。
【００５５】
図９は、通常の方向距離法を適用して注目ボクセルの方向距離を推定する処理の概要を示す。２値化された入力画像Ｇ４１１が入力された場合、注目ボクセルについて、背景までの方向距離を計算する。注目ボクセルに近接するボクセルへの方向は、２次元の場合には８方向、３次元の場合には２６方向となる。
【００５６】
ここでは２次元の場合について説明する。向きが逆方向のベクトルを除いたものをｖｓｉとすると、４方向に画像上を進行して背景に到達した時点で進行を停止する。−ｖｓｉについても同様に４方向に画像上を進行して背景に到達した時点で進行を停止する。
【００５７】
このときｖｓｉ方向に進んだ距離をｄ（ｖｓｉ）、−ｖｓｉ方向に進んだ距離をｄ（−ｖｓｉ）とすると、±ｖｓｉ方向の方向距離ｄ（±ｖｓｉ）は、下記式１により示される。
【００５８】
【数１】
2016052489.tif 000003
【００５９】
また方向ベクトルＰｉ及び方向テンソルＭをそれぞれ下記式２及び３により定義する。方向テンソルＭを固有値分解した場合に最大の固有値に対する固有ベクトルが繊維束の配向を示す。また２番目に大きい固有値に対する固有ベクトルが繊維束の幅方向を示す。３次元の場合には最小の固有値に対する固有ベクトルが厚さ方向を示す。
【００６０】
【数２】
2016052489.tif 000004
【００６１】
【数３】
2016052489.tif 000005
【００６２】
図８に戻り、ステップＳＰ４１で推定した重なり方向は、方向距離法の性質上、重なり領域と背景との境界部分のボクセルについては正しく計算されない。この正しく計算されない境界部分のボクセルを取り除くため、リファレンス方向決定部１５３は、モルフォロジー処理のうちの収縮処理を実行し、重なり領域から中心領域を抽出する（ＳＰ４２）。
【００６３】
次いでリファレンス方向決定部１５３は、重なり領域のそれぞれが独立した連結成分を構成していることから、抽出した各中心領域に対して６近傍のラベリング処理を実行して領域分割し、重なり領域を分離する（ＳＰ４３）。
【００６４】
次いでリファレンス方向決定部１５３は、ステップＳＰ４１で推定した重なり方向のうち、中心領域の各ボクセルにおける重なり方向を平均化し、平均化して得られた方向をリファレンス方向として決定して（ＳＰ４４）、本処理を終了する。
【００６５】
図１０は、リファレンス方向決定処理により作成される処理画像を示す。具体的には重なり方向推定画像Ｇ４１、中心領域抽出画像Ｇ４２、ラベリング画像Ｇ４３及びリファレンス方向画像Ｇ４４を示す。
【００６６】
図１０に示すように重なり領域抽出画像Ｇ３３（図７）に対し、通常の方向距離法を適用することで重なり方向推定画像Ｇ４１を得ることができる。また重なり方向推定画像Ｇ４１に対し、収縮処理を実行することで中心領域抽出画像Ｇ４２を得ることができる。
【００６７】
また中心領域抽出画像Ｇ４２に対し、ラベリング処理を実行することでラベリング画像Ｇ４３を得ることができる。また中心領域における重なり方向を平均化することでリファレンス方向画像Ｇ４４を得ることができる。
【００６８】
（３−４）Ｚ糸除去処理
図１１は、Ｚ糸除去処理の詳細なフローチャートを示す。Ｚ糸除去部１５４は、ステップＳＰ２で得られた２値画像Ｇ２２（図５）に対してリファレンス付き方向距離法を適用し、方向距離を計算する（ＳＰ５１）。
【００６９】
リファレンス付き方向距離法とは、方向距離を計算する平面をリファレンス方向に垂直な平面（リファレンス平面）に回転させて、このリファレンス平面上で２次元の通常の方向距離法により方向距離を計算する手法である。
【００７０】
図１２は、通常の方向距離法を適用して方向距離を計算する場合の概念図を示す。Ｘ糸とＹ糸とが立体交差する重なり領域において通常の方向距離法を適用して方向距離を計算すると、正確な方向距離を計算することができないという問題がある。
【００７１】
具体的には繊維束の配向に近い方向は、Ｘ糸及びＹ糸の繊維束の厚みが薄いためにすぐに背景に到達して距離が短くなる一方で、リファレンス方向に近い方向は、重なり領域の断面が太いためなかなか背景に到達せずに距離が長くなる。
【００７２】
従って繊維束の配向に近い方向の距離よりも、リファレンス方向に近い方向の距離の方が長くなる。この結果、リファレンス方向に近い方向の方向距離が重なり領域における方向距離として計算されるという問題がある。そこで本実施の形態においてはリファレンス付き方向距離法を適用して方向距離を計算するようにしている。
【００７３】
図１３は、リファレンス付き方向距離法の概念図を示す。リファレンス付き方向距離法とは、方向距離の計算対象となる平面を回転させて、回転させた後の平面上で２次元の方向距離法を適用して方向距離を計算する手法である。
【００７４】
具体的には方向±ｖｓ０をリファレンス方向に一致させるように回転させ、±ｖｓ０を回転させた分だけ、±ｖｓ０に垂直な他の方向±ｖｓ１及び±ｖｓ２も回転させる。回転させた後の方向をそれぞれ±ｖｓ０´、±ｖｓ１´及び±ｖｓ２´とすると、Ｘ糸又はＹ糸は、±ｖｓ１´及び±ｖｓ２´により定義される平面（リファレンス平面）上に存在する。
【００７５】
このリファレンス平面上で２次元の方向距離法を適用して方向距離を計算することにより、重なり領域の影響を受けずに繊維束の配向を推定することができる。
【００７６】
図１４は、リファレンス付き方向距離法を適用してリファレンス平面上のＸ糸又はＹ糸の方向距離を計算する場合の概念図を示す。重なり領域内に位置する注目ボクセルについて方向距離を計算する場合、リファレンス付き方向距離法を適用して方向距離を計算する。
【００７７】
具体的には上記の通り、方向±ｖｓ０をリファレンス方向に一致させるように回転させ、±ｖｓ０と垂直な他の方向±ｖｓ１、±ｖｓ２、±ｖｓ３及び±ｖｓ４も同程度回転させる。
【００７８】
回転させた後の方向をそれぞれ±ｖｓ０´±ｖｓ１´、±ｖｓ２´、±ｖｓ３´及び±ｖｓ４´とすると、±ｖｓ１´、±ｖｓ２´、±ｖｓ３´及び±ｖｓ４´により定義されるリファレンス平面上で２次元の方向距離法を適用して方向距離を計算する。
【００７９】
リファレンス付き方向距離法を適用して方向距離を計算する場合についても、上記式１〜３に示す関係が成立する。すなわち方向テンソルＭを固有値分解した場合に最大の固有値に対する固有ベクトルが繊維束の配向を示す。また２番目に大きい固有値に対する固有ベクトルが繊維束の幅方向を示す。３次元の場合には最小の固有値に対する固有ベクトルが厚さ方向を示す。
【００８０】
図１１に戻り、Ｚ糸除去部１５４は、リファレンス付き方向距離法の計算の際に得られる２番目に大きい固有値に対する固有ベクトルを用いて、条件付き方向距離を計算する（ＳＰ５２）。そしてＺ糸除去部１５４は、条件付き方向距離に基づいて、Ｘ糸又はＹ糸の領域と、Ｚ糸の領域とを分離する（ＳＰ５３）。
【００８１】
条件付き方向距離とは、注目ボクセルから画像上を進行して背景に到達した時点又は次に進行する位置にあるボクセルの固有ベクトルが示す方向と、進行方向とのなす角が所定の閾値よりも大きい時点で、進行を停止するとの条件下で計算される方向距離である。
【００８２】
図１５は、条件付き方向距離を推定する処理の概要を示す。位置ｘにある注目ボクセルを起点として、注目ボクセルの固有ベクトルｖ２（ｘ）が示す方向に画像上を進行し、背景に到達した時点で進行を停止する。最後に到達したボクセルから注目ボクセルまでの距離をｘｅ１とする。
【００８３】
一方で逆方向のベクトル−ｖ２（ｘ）の方向についても画像上を進行し、次に進行する位置ｘ´にあるボクセルの固有ベクトルｖ２（ｘ´）が示す方向と、−ｖ２（ｘ）が示す方向とのなす角が所定の閾値よりも大きい時点で進行を停止する。最後に到達したボクセルから注目ボクセルまでの距離をｘｅ２とする。そして距離ｘｅ１及びｘｅ２の絶対値の合計を条件付き方向距離ｄｃ（ｘ）とする。
【００８４】
図１６は、条件付き方向距離に基づいてＸ糸又はＹ糸の領域とＺ糸の領域とを分離した概念図を示す。Ｘ糸又はＹ糸の幅Ｗｘｙは、Ｚ糸の断面の対角線距離Ｗｚよりも大きいことから、Ｘ糸又はＹ糸の領域におけるボクセルｐの条件付き方向距離ｄｃ（ｐ）は、Ｚ糸の領域におけるボクセルｑの条件付き方向距離ｄｃ（ｑ）よりも大きい。
【００８５】
例えばＺ糸の断面の対角線距離Ｗｚよりも大きく、かつ、Ｘ糸又はＹ糸の幅Ｗｘｙよりも小さい値を閾値として、この閾値よりも大きい条件付き方向距離のボクセルをＸ糸又はＹ糸の領域に含め、閾値よりも小さい条件付き方向距離のボクセルをＺ糸の領域に含めることができる。
【００８６】
図１１に戻り、Ｚ糸除去部１５４は、Ｚ糸の領域に含まれるボクセルを削除することにより、２値画像Ｇ２２（図５）からＺ糸を除去する（ＳＰ５４）。Ｚ糸を除去するとＸ糸又はＹ糸の領域の一部も除去される場合がある。
【００８７】
そこでＺ糸除去部１５４は、Ｚ糸除去後の画像に対してモルフォロジー処理のうちの膨張処理を実行し（ＳＰ５５）、膨張処理後の画像と膨張処理前の画像との積集合を計算してＸ糸及びＹ糸を抽出し（ＳＰ５６）、本処理を終了する。
【００８８】
（３−５）繊維束配向推定処理
図１７は、繊維束配向推定処理の詳細なフローチャートを示す。繊維束配向推定部１５５は、Ｘ糸及びＹ糸を抽出した画像に対して再度リファレンス付き方向距離法を適用し、Ｘ糸及びＹ糸の繊維方向を推定する（ＳＰ６１）。
【００８９】
次いで繊維束配向推定部１５５は、繊維方向を推定したボクセルのうち、比較的精度が良いとされる中心部近傍のボクセルだけを残して他のボクセルを削除し、中心部近傍のボクセルを抽出する（ＳＰ６２）。
【００９０】
次いで繊維束配向推定部１５５は、中心部近傍のボクセルについて、同程度の方向のボクセルを繋いでいき、繋がったボクセル群を同一クラスタに所属させ、繊維束の途中でベクトルが急激に変化しているボクセルは異なるクラスタに所属させるクラスタリング処理を実行する。
【００９１】
繊維束の配向は数ボクセル単位の短い距離で急激に変化することはない。よって繊維束の途中でベクトルが急激に変化しているボクセルはノイズであると判断することができる。そこで繊維束配向推定部１５５は、所属するボクセルが少ないクラスタを削除してノイズを除去して（ＳＰ６３）、本処理を終了する。
【００９２】
（４）解析結果
以下図１８〜図２１を参照して、各種入力画像に対して上記説明してきた本実施の形態における画像解析処理を実行した場合に得られる処理画像について説明する。
【００９３】
図１８は、入力画像が模擬画像Ｇ１１である場合に画像解析処理を実行して得られる処理画像を示す。模擬画像Ｇ１１は、Ｘ糸及びＹ糸の繊維束の長さ、幅、厚さ、Ｘ糸の繊維層及びＹ糸の繊維層の積層数、繊維束同士の間隔並びにＺ糸の断面形状等を設定して作成される３次元画像である。すなわち各種情報が既知の画像である。
【００９４】
この模擬画像Ｇ１１に対して本実施の形態における画像解析処理を実行した場合、３次元画像Ｇ１１０、Ｘ糸断面画像Ｇ１１１及びＹ糸断面画像Ｇ１１２を得ることができる。３次元画像Ｇ１１０は、Ｘ糸及びＹ糸の繊維束の配向を示す３次元画像である。
【００９５】
またＸ糸断面画像Ｇ１１１は、Ｘ糸の繊維束の配向を示す２次元画像であり、Ｙ糸断面画像Ｇ１１２は、Ｙ糸の繊維束の配向を示す２次元画像である。
【００９６】
これらの処理画像Ｇ１１０〜Ｇ１１２を得る際に計算された繊維束の方向と、模擬画像Ｇ１１の作成の際に設定された繊維束の方向とのなす角を誤差として計算した結果、誤差の最大値は８９．９度、平均値は４．９度であった。
【００９７】
なお一般に方向距離法自体の平均誤差は、４．３度であることが知られている。また全ボクセルのうち、９４．７％のボクセルの誤差が６度以内であった。以上の結果から本実施の形態における画像解析処理の有効性を確認することができる。
【００９８】
図１９は、入力画像がＸ線ＣＴ画像Ｇ１０である場合に画像解析処理を実行して得られる処理画像を示す。Ｘ線ＣＴ画像Ｇ１０は、実際の織物をＸ線ＣＴ装置で撮影して得られる３次元画像である。またＸ線ＣＴ画像Ｇ１０は、図２０に示すＸ線ＣＴ画像Ｇ１２と比較して高解像度である。
【００９９】
このＸ線ＣＴ画像Ｇ１０に対して本実施の形態における画像解析処理を実行した場合、３次元画像Ｇ１００（繊維束配向推定画像Ｇ１００）、Ｘ糸断面画像Ｇ１０１及びＹ糸断面画像Ｇ１０２を得ることができる。３次元画像Ｇ１００は、Ｘ糸及びＹ糸の繊維束の配向を示す３次元画像である。
【０１００】
またＸ糸断面画像Ｇ１０１は、Ｘ糸の繊維束の配向を示す２次元画像であり、Ｙ糸断面画像Ｇ１０２は、Ｙ糸の繊維束の配向を示す２次元画像である。これらの処理画像Ｇ１００〜Ｇ１０２を参照することにより、Ｘ糸及びＹ糸の繊維束の配向を容易に確認することができる。
【０１０１】
また計算時間は２４３．８秒程度であった。従来同程度のボクセル数の画像を解析しようとすると、２０倍程度の時間が必要であったことから、本実施の形態における画像解析処理を実行することで、計算時間の短縮化を図ることができる。
【０１０２】
図２０は、入力画像がＸ線ＣＴ画像Ｇ１２である場合に画像解析処理を実行して得られる処理画像を示す。Ｘ線ＣＴ画像Ｇ１２は、図１９に示すＸ線ＣＴ画像Ｇ１０と比較して低解像度である点で異なる。
【０１０３】
このＸ線ＣＴ画像Ｇ１２に対して本実施の形態における画像解析処理を実行した場合、３次元画像Ｇ１２０、Ｘ糸断面画像Ｇ１２１及びＹ糸断面画像Ｇ１２２を得ることができる。これらの処理画像Ｇ１２０〜Ｇ１２２を参照すると、一部の繊維束が欠落しているものの、おおよその繊維束の配向を容易に確認することができる。
【０１０４】
図２１は、入力画像が欠陥模擬画像Ｇ１３である場合に画像解析処理を実行して得られる処理画像を示す。欠陥模擬画像Ｇ１３は、図１８に示す模擬画像Ｇ１１と比較してＸ糸の繊維束が途中で１０度折れ曲がっている欠陥を有する点で異なる。
【０１０５】
この欠陥模擬画像Ｇ１３に対して本実施の形態における画像解析処理を実行した場合、３次元画像Ｇ１３０及びＸ糸断面画像Ｇ１３１を得ることができる。これらの処理画像Ｇ１３０及びＧ１３１を参照すると、多少の誤差はあるものの、折れ曲がり地点からＸ糸の繊維束が折れ曲がっていることを容易に確認することができる。
【０１０６】
図２２は、入力画像が曲面模擬画像Ｇ１４である場合に画像解析処理を実行して得られる処理画像を示す。曲面模擬画像Ｇ１４は、図１８に示す模擬画像Ｇ１１と比較して円弧上に変形した曲面部分を有する点で異なる。
【０１０７】
この曲面模擬画像Ｇ１４に対して本実施の形態における画像解析処理を実行した場合、３次元画像Ｇ１４０及びＸ糸断面画像Ｇ１４１を得ることができる。これらの処理画像Ｇ１４０及びＧ１４１を参照すると、繊維束が円弧上に配向していることを容易に確認することができる。
【０１０８】
（５）本実施の形態による効果
以上のように本実施の形態における画像解析装置、画像解析方法及びプログラムによれば、３次元織りにより作製された織物のＸ線ＣＴ画像に対してリファレンス付き方向距離法を適用してＺ糸を除去し、Ｚ糸を除去した画像に対して再度リファレンス付き方向距離法を適用してＸ糸及びＹ糸の繊維束の配向を推定するようにしたので、Ｚ糸の影響を除去して繊維束の配向を精度良く短時間で推定することができる。また曲面形状を有する織物のＸ線ＣＴ画像についても繊維束の配向を推定することができる。よって実際の製品検査に利用することができる。
【０１０９】
次に、本発明の一実施の形態である画像解析装置２について図２３〜図４６を参照して説明する。
【０１１０】
（１）画像解析装置２の全体構成
この実施形態における画像解析装置２は、図２３に示されるように、図１の画像解析装置１と同様のＣＰＵ１１、入力部１２、記憶部１３、表示部１４及びメモリ１５を備え、メモリ１５は２値化部１５１及び重なり領域抽出部１５２を有している。これらの構成は、図１において同一の名称及び符号を有する構成と同様であり、同様の処理をおこなうものであるため、重複する説明を省略する。ここで、本実施形態の画像解析装置２では、メモリ１５内に更に、重なり領域形態解析部２５３を有している。重なり領域形態解析部２５３は、重なり領域抽出部１５２の処理により得られた重なり領域抽出画像Ｇ３３における「重なり領域」の形態を解析する。ここで、重なり領域の形態には、一つの重なり領域の形態だけでなく、複数の点在する重なり領域を組み合わせた形態も含むものとする。また、図面において、画像解析装置２は、表示部１４を有しているが、表示は行わずに解析を行うこととした表示部１４を有さない構成であってもよい。
【０１１１】
（２）画像解析処理Ｐ２のフローチャート
図２４には、画像解析装置２で行われる画像解析処理Ｐ２のフローチャートが示されている。このフローチャートに示されるように、画像解析処理Ｐ２では、まず、Ｘ線ＣＴ画像の入力（ＳＰ１）、２値化処理（ＳＰ２）、及び重なり領域抽出処理（ＳＰ３）が行われる。これらの処理は、上述の図３に示されたＳＰ１〜ＳＰ３の処理と同様であるため、説明を省略する。次に、重なり領域形態解析処理（ＳＰ１４）において、重なり領域形態解析部２５３により、重なり領域抽出処理により得られた重なり領域抽出画像Ｇ３３（図７参照）を用いて重なり領域２１の形態解析の処理を行う。ここで、本実施形態においては、以降、重なり領域抽出画像Ｇ３３における個々の重なり領域２１を「ピラー２１」と呼ぶこととする。
【０１１２】
（３）重なり領域形態解析処理の例
重なり領域形態解析部２５３は、例えば、３次元画像のピラー２１に含まれるボクセル数を数える等により各ピラー２１の体積を算出することとしてもよい。図２５においては、２次元画像において黒点でしめされるピラー２１内のピクセルを例示的に示しているが、３次元画像においても同様にボクセル数を数えることができる。この場合には、例えば、複数のピラー２１の体積の平均を参照値として、特定のピラー２１の体積と比較することにより、特定のピラー２１の体積が適当であるかどうかを判定することとしてもよい。ここで体積の参照値には、複数のピラー２１の平均だけでなく、設計上の基準値その他の値を用いることができる。また、ピラー２１の体積を算出することとしたが、２次元画像を用いて解析する場合には面積を算出することとしてもよい。このような体積を算出する処理により、周囲のピラー２１の体積と比較して異なる体積である等、繊維束の配向異常の検出に有用な情報を提供することができる。また提供された情報に基づいて、例えば、立体交差すべきＸ糸又はＹ糸の数の不足の検出等の繊維束の配向異常を検出することができる。
【０１１３】
また、重なり領域形態解析部２５３は、ピラー２１が伸長する方向を算出することとしてもよい。ピラー２１が伸長する方向は、例えば図２６に示されるように、モルフォロジー画像処理の収縮処理を用いて伸長方向を顕在化し（Ｓ２１の斜線部）、ベクトル化してもよいし（Ｓ２２）、その他の画像処理を用いて伸長する方向を求めることができる。また、例えば、重なり領域形態解析部２５３は、複数のピラー２１の伸長する方向の平均を参照値として、特定のピラー２１の伸長する方向と比較することにより、ピラー２１の伸長する方向が適当であるかどうかを判定することができる。ここで伸長する方向の参照値には、複数のピラー２１の伸長する方向の平均だけでなく、設計上の基準値等その他の値を用いることができる。また、ピラー２１の伸長する方向とは別に、ピラー２１の中立軸の形状を求めることとしてもよい。中立軸は、例えば、図２７の斜線部に示されるように、モルフォロジー画像処理の収縮処理を用いて抽出することができる。この場合には、ピラー２１がある一定方向に伸長することなく、折れ曲がっている場合等についても検出することができ、更に参照形状と比較することにより、適当な形状かどうか判定することができる。このようなピラー２１の伸長する方向や中立軸を算出する処理により、例えば、積層にずれが生じている部分等その他の繊維束の配向の異常の検出に有用な情報を提供することができる。また提供された情報に基づいて繊維束の配向異常を検出することができる。
【０１１４】
また、重なり領域形態解析部２５３は、複数のピラー２１の重心位置を算出することとしてもよい。隣り合うピラー２１の重心位置を結ぶ線は、理想的には、Ｘ糸及びＹ糸の伸びる方向となるため、ＸＹ平面において、これらが滑らかな線上に並ばない等を検出し、不規則な重心位置を求めることにより、重心位置が適当であるかどうか判定することができる。図２８は、斜線で示されるピラー２１の重心Ｇが、周囲のピラー２１の重心Ｇと比較して滑らかに並んでいない様子を示している。また、ある領域を定め、それに含まれる重心Ｇの数を数えること等により、領域内に重心Ｇが適当な数だけ存在するかどうかを求めることにより、不規則な配置のピラー２１が存在するかどうかを判定することとしてもよい。このような重心位置を算出する処理により、繊維束の配向の異常の検出に有用な情報を提供することができる。また提供された情報に基づいて繊維束の配向異常を検出することができる。
【０１１５】
上述の重なり領域形態解析処理の例では、繊維束の配向異常の検出の処理を含む場合についても記載したが、重なり領域形態解析部２５３の処理には、このような検出を含まずに、ピラー２１の体積を算出、ピラー２１の伸長する方向を算出、又はピラー２１の重心位置を算出等の重なり領域形態解析の処理のみを行うこととし、例えば、表示部１４に出力する、又は他の装置に算出されたデータを送信する等により、処理された情報を提供することのみであってもよい。このような処理により、繊維束の配向異常の検出に有用な情報を提供することができる。
【０１１６】
（４）繊維束の配向異常の具体例
以下、図２９〜４６を用いてピラー２１の形態解析により検出される繊維束の配向異常の具体例について説明する。繊維束の配向異常の検出の具体例においては、重なり領域形態解析部２５３により行われる重なり領域形態解析処理において、繊維束の配向異常の検出を含む場合についても記載するが、重なり領域形態解析部２５３の処理には、このような判定を含まずに、ピラー２１の形態解析により得られた繊維束の配向についての情報を提供することのみでもよい。なお、図３０，３３，３６，３８，４０，４１，４３，４４及び４６は、ＦＡＡ（Federal Aviation Administration）発行のレポート”DOT/FAA/AR-06/10” を参考に記載している。
【０１１７】
図２９は、Ｘ糸及びＹ糸の繊維束の理想的な配置の例について、ＸＹ平面で示す図であり、図３０及び３１は、それぞれ、図２９のような理想的な繊維束の配置における、ＸＺ平面（又はＹＺ平面）及びＸＹ平面における重なり領域抽出画像を模式的に表す図である。図３０に示されるように、Ｘ糸及びＹ糸の繊維束の理想的な配置においては、ピラー２１はＸ糸及びＹ糸の繊維束の積層方向（Ｚ方向）に途切れなく伸び、図３１に示されるように、ピラー２１（又はピラー２１の重心Ｇ）は規則性を有して並んでいる。
【０１１８】
図３２（−Ｙ方向）は、ＸＹ平面において、繊維束の一部が延伸すべき方向とは別のＸＹ平面内の方向に伸びている状態について示す図である。図３３及び３４は、それぞれ図３２の繊維束の状態における、ＹＺ平面及びＸＹ平面の重なり領域抽出画像を模式的に示す図である。これらの図に示されるように、図３３の重なり領域抽出画像においては、−Ｙ方向に伸びる繊維束（斜線部）の影響により、それを有するピラー２１１の重心位置Ｇが−Ｙ方向に移動している。図３４の重なり領域抽出画像においては、ピラー２１１の重心位置Ｇが、近傍のピラー２１の重心位置Ｇの規則性から外れている。したがって、重なり領域形態解析部２５３は、例えば、重なり領域抽出画像におけるそれぞれのピラー２１の重心位置Ｇを解析し、繊維束の一部が延伸すべき方向とは別のＸＹ平面内の方向に伸びている領域を検出することができる。
【０１１９】
図３５は、ＸＹ平面において、繊維束の一部が延伸すべき方向とは別のＺ方向に伸びている状態について示す図である。図３６は、図３５の繊維束の状態における、ＹＺ平面の重なり領域抽出画像を模式的に示す図である。この図に示されるように、図３６の重なり領域抽出画像においては、Ｚ方向に伸びる繊維束（斜線部）の影響により、ピラー２１２が複数に分割されている。したがって、重なり領域形態解析部２５３は、重なり領域抽出画像を用いて、繊維束の一部が延伸すべき方向とは別のＺ方向に伸びている領域を検出することができる。
【０１２０】
図３７は、ＸＹ平面において、繊維束の一部が延伸すべき方向とは別の、ＸＹ平面内の方向成分及びＺ方向成分を共に有する方向に伸びている状態について示す図である。図３８は、図３７の繊維束の状態における、ＹＺ平面の重なり領域抽出画像を模式的に示す図である。この図に示されるように、図３８の重なり領域抽出画像においては、−Ｚ方向及び−Ｙ方向に伸びる繊維束（斜線部）の影響により、ピラー２１３が、複数に分割されている。また、分割されたピラー２１３の一方の重心位置Ｇは、＋Ｙ方向に移動している。したがって、重なり領域形態解析部２５３は、重なり領域抽出画像を用いて、繊維束の一部が延伸すべき方向とは別のＸＹ平面内の方向成分及びＺ方向成分を共に有する方向に伸びている領域を検出することができる。
【０１２１】
図３９及び４０は、ＸＹ平面において、それぞれ繊維束が全体として傾いた方向に伸びている状態について示す図である。図４１は、図３９又は４０に示される繊維束の状態における、ＸＹ平面の重なり領域抽出画像を模式的に示す図である。この図に示されるように、図４１の重なり領域抽出画像においては、規則的に並ぶピラー２１の並びの一部の列であるピラー２１４が、ずれて配置されている。したがって、重なり領域形態解析部２５３は、重なり領域抽出画像を用いて、例えば、ピラー２１４の重心位置Ｇの並びが、所定の規則性から外れること等を検出することにより、繊維束が全体として傾いた方向に伸びている領域を検出することができる。
【０１２２】
図４２は、ＸＹ平面において、繊維束の一部（例えばＸ糸）がない状態について示す図である。図４３は、ＸＹ平面において、繊維束の一部の太さが他の繊維束より太く又は細くなっている場合について示す図である。図４４は、図４２又は図４３の繊維束の状態における、ＹＺ平面の重なり領域抽出画像を模式的に示す図である。図４４の重なり領域抽出画像においては、一部のピラー２１５が、複数に分割され、分割された複数のピラー２１５の体積を合わせても、通常のピラー２１の体積よりも小さくなっている。したがって、重なり領域形態解析部２５３は、重なり領域抽出画像を用いて、繊維束の一部がない場合、及び繊維束の一部の太さが他の繊維束より太く又は細くなっている領域について検出することができる。
【０１２３】
図４５は、ＸＹ平面において、繊維束の一部（例えばＸ糸）が折り返されている状態について示す図である。図４６は、図４５の繊維束の状態における、ＸＹ平面の重なり領域抽出画像を模式的に示す図である。図４６の重なり領域抽出画像においては、繊維束の一部の折り返しにより、折り返しの周囲のピラー２１６が押しのけられ、その周りのピラー２１の規則性のある重心位置Ｇと比較して、ピラー２１６の重心位置Ｇに乱れが生じている。したがって、重なり領域形態解析部２５３は、重なり領域抽出画像を用いて、これらの重心位置Ｇの乱れを解析することにより、繊維束の一部が折り返されている領域について検出することができる。
【０１２４】
（５）本実施の形態の効果
以上のように本実施の形態の画像解析装置２では、２値化部１５１が、Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像を２値化し、重なり領域抽出部１５２が、２値化された画像からＸ糸及びＹ糸が立体交差する重なり領域を抽出し、重なり領域形態解析部２５３が、抽出された重なり領域の形態を解析する。これにより処理時間の短い画像処理の組み合わせのみで解析を行うため、より簡易に短時間でＸ糸及びＹ糸の繊維束の配向についての情報を提供することができる。
【０１２５】
また、重なり領域形態解析部２５３が、更に繊維束の配向異常を検出する場合には、より簡易に短時間で繊維束の配向異常の箇所を見つけることができる。また、繊維束の配向異常について詳細な解析が必要な場合であっても、検出された配向異常を含む領域についてのみ別途詳細な解析を行うこととすればよいため、全体として解析時間の短縮を図ることができる。
【０１２６】
また、本実施の形態における重なり領域の形態を解析は、繊維束全体の形状が平面である場合、曲面である場合に関わらず、用いることができる。重なり領域は、繊維束全体がどのような形状であっても抽出することができるため、重なり領域の形態の解析処理を適用することができる。また、ピラー２１は、繊維束全体の形状に沿って規則的に配置されるため、繊維束全体の形状が平面である場合、曲面である場合に関わらず、繊維束の配向の異常等についても、重なり領域の形態を解析することにより検出することができる。
【０１２７】
また、重なり領域形態解析部２５３は、重なり領域の体積を算出してもよい。この場合に、重なり領域の「体積の参照値」と比較することができ、更にこの「体積の参照値」は、複数の重なり領域の体積の平均値とすることができる。
【０１２８】
また、重なり領域形態解析部２５３は、重なり領域が伸長する方向を算出してもよい。この場合に、重なり領域の「向きの参照値」と比較することができ、この「向きの参照値」は、複数の前記重なり領域の方向の平均値とすることができる。
【０１２９】
また、重なり領域形態解析部２５３は、複数の重なり領域の重心位置を算出してもよい。この場合に、重心位置の配置が不規則である領域を検出することができる。例えば重心位置が並ぶ参照線を求め、その参照線からどの程度離れているかを算出、又はある領域において含まれる重心位置の数を算出する等により、重心位置の配置が不規則である領域を検出することとしてもよい。
【０１３０】
また、重なり領域形態解析部２５３は、重なり領域の中立軸を算出し、参照形状と比較してもよい。この場合に、中立軸の算出には、モルフォロジーの収縮処理を用いることができる。
【０１３１】
また、本実施の形態の画像解析方法は、Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像を２値化するステップと、前記２値化された画像から前記Ｘ糸及び前記Ｙ糸が立体交差する重なり領域を抽出するステップと、前記抽出された重なり領域の形態を解析するステップと、を備えることを特徴とする画像解析方法である。
【０１３２】
また、本実施の形態のプログラムは、コンピュータに、Ｘ糸、Ｙ糸及びＺ糸の繊維束からなる織物の３次元画像を２値化するステップと、前記２値化された画像から前記Ｘ糸及び前記Ｙ糸が立体交差する重なり領域を抽出するステップと、前記抽出された重なり領域の形態を解析するステップとを実行させるためのプログラムである。
【０１３３】
（６）他の実施の形態
上述の各実施の形態においては、図１に示される画像解析装置１による画像解析処理（図３）、及び図２３に示される画像解析装置２による画像解析処理Ｐ２（図２４）を別々に使用する場合について述べたが、これらが関連して動作するものとすることができる。図４７には、これらの処理を組み合わせて行う実施の形態の一つである画像解析装置３の構成が示されている。画像解析装置３は、上述の画像解析装置１又は２と同様に、ＣＰＵ１１、入力部１２、記憶部１３、表示部１４及びメモリ１５を備えており、メモリ１５は、２値化部１５１、重なり領域抽出部１５２、リファレンス方向決定部１５３、Ｚ糸除去部１５４、繊維束配向推定部１５５及び重なり領域形態解析部２５３を有している。各構成は、上述した対応する構成の処理と同様の処理を行う。
【０１３４】
図４８には、この画像解析装置３の画像解析処理Ｐ３のフローチャートが示されている。画像解析処理Ｐ３は、まず、図２４の画像解析処理Ｐ２と同様に、Ｘ線ＣＴ画像の入力（ＳＰ１）、２値化処理（ＳＰ２）、重なり領域抽出処理（ＳＰ３）及び重なり領域形態解析処理（ＳＰ１４）を行う。ここで重なり領域形態解析処理では、重なり領域の形態を解析だけでなく、繊維束の配向異常の検出まで行う。次に、検出された繊維束の配向異常を含む領域について、図３の画像解析処理Ｐ１における対応する処理と同様に、リファレンス方向決定処理（ＳＰ４）、Ｚ糸除去処理（ＳＰ５）、繊維束配向推定処理（ＳＰ６）を行う。
【０１３５】
このように、処理時間の短い画像解析処理Ｐ２を、例えば検査対象となるＸ線ＣＴ画像に対して適用し、画像解析処理Ｐ２で検出された、繊維束の配向異常を含む領域について、比較的処理時間が長く、詳細な解析を行う画像解析処理Ｐ１を適用することにより、より効率的にかつ詳細に繊維束の配向の異常について解析することができる。
【符号の説明】
【０１３６】
１      画像解析装置
１１    ＣＰＵ
１２    入力部
１３    記憶部
１４    表示部
１５    メモリ
１５１  ２値化部
１５２  重なり領域抽出部
１５３  リファレンス方向決定部
１５４  Ｚ糸除去部
１５５  繊維束配向推定部
２      画像解析装置
２１    ピラー
２５３  重なり領域形態解析部
３      画像解析装置
【図１】
2016052489.tif 000006
【図２】
2016052489.tif 000007
【図３】
2016052489.tif 000008
【図４】
2016052489.tif 000009
【図５】
2016052489.tif 000010
【図６】
2016052489.tif 000011
【図７】
2016052489.tif 000012
【図８】
2016052489.tif 000013
【図９】
2016052489.tif 000014
【図１０】
2016052489.tif 000015
【図１１】
2016052489.tif 000016
【図１２】
2016052489.tif 000017
【図１３】
2016052489.tif 000018
【図１４】
2016052489.tif 000019
【図１５】
2016052489.tif 000020
【図１６】
2016052489.tif 000021
【図１７】
2016052489.tif 000022
【図１８】
2016052489.tif 000023
【図１９】
2016052489.tif 000024
【図２０】
2016052489.tif 000025
【図２１】
2016052489.tif 000026
【図２２】
2016052489.tif 000027
【図２３】
2016052489.tif 000028
【図２４】
2016052489.tif 000029
【図２５】
2016052489.tif 000030
【図２６】
2016052489.tif 000031
【図２７】
2016052489.tif 000032
【図２８】
2016052489.tif 000033
【図２９】
2016052489.tif 000034
【図３０】
2016052489.tif 000035
【図３１】
2016052489.tif 000036
【図３２】
2016052489.tif 000037
【図３３】
2016052489.tif 000038
【図３４】
2016052489.tif 000039
【図３５】
2016052489.tif 000040
【図３６】
2016052489.tif 000041
【図３７】
2016052489.tif 000042
【図３８】
2016052489.tif 000043
【図３９】
2016052489.tif 000044
【図４０】
2016052489.tif 000045
【図４１】
2016052489.tif 000046
【図４２】
2016052489.tif 000047
【図４３】
2016052489.tif 000048
【図４４】
2016052489.tif 000049
【図４５】
2016052489.tif 000050
【図４６】
2016052489.tif 000051
【図４７】
2016052489.tif 000052
【図４８】
2016052489.tif 000053
【国際調査報告】
2016052489.tif 000054
2016052489.tif 000055
2016052489.tif 000056
2016052489.tif 000057