(57)【特許請求の範囲】
【請求項１】
母材に繊維状のフィラーを含有させた複合材料の３次元画像から前記フィラーの配向情報を抽出する複合材料中の繊維状フィラーの３次元画像処理方法において、
前記複合材料の３次元画像からフィラーを表す画素を含むと推定される候補領域を所定の閾値を参照することにより選出する領域選出工程と、
フィラーの形状モデルの形状と配置を変化させるパラメータをランダムに変化させるモンテカルロ法を用いて、前記領域選出工程によって選出された候補領域に前記形状モデルをフィッティングさせてフィラーの配向情報を抽出する抽出工程と、を備えることを特徴とする複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項２】
前記抽出工程は、前記候補領域から１本のフィラーの配向情報を抽出した後に、前記抽出に関与した領域を除去し、除去後の候補領域に対して前記モンテカルロ法による抽出を繰り返すことを特徴とする請求項１に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項３】
前記３次元画像をボクセルデータとして入力するデータ入力工程を備え、
前記ボクセルデータは、複合材料のＸ線ＣＴによって取得され、前記ボクセルデータの各ボクセルは、Ｘ線の強度値に基づく値をデータ値として有することを特徴とする請求項１または請求項２に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項４】
前記データ入力工程は、前記データ値を各ボクセル間で線形補間した値をデータ値とするボクセルを有する新たなボクセルデータを構成する補間工程を備えることを特徴とする請求項３に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項５】
前記領域選出工程は、前記ボクセルのデータ値とＸ線の強度値に基づく所定の閾値とを比較し、前記閾値よりも大きいデータ値を有するボクセルの集合を候補領域とすることを特徴とする請求項３または請求項４に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項６】
前記候補領域に属するボクセルから、互いに隣接するボクセルで成る群を生成し、前記生成された個々の群をそれぞれ新たな候補領域とすることを特徴とする請求項５に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項７】
前記抽出工程は、前記形状モデルとして仮想円柱を用い、前記候補領域内で前記仮想円柱の形状と配置を表すパラメータを変化させ、その仮想円柱と前記候補領域とのフィッティングの度合いを、前記仮想円柱に含まれるボクセルのデータ値に基づく評価値の積算値によって評価することを特徴とする請求項３乃至請求項６の何れか一項に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項８】
前記抽出工程は、前記仮想円柱によるフィッティングを行った後に、前記形状モデルとして複数の制御点によって定義されるスプライン曲線を中心軸とする仮想パイプを用い、前記候補領域内で前記制御点の座標をパラメータとしてランダムに変化させ、その仮想パイプと前記候補領域とのフィッティングの度合いを評価し、前記仮想円柱による評価よりも前記仮想パイプによる評価が所定割合以上改善する場合に、前記仮想パイプによるフィッティングを採用し、そうでない場合は前記仮想円柱によるフィッティングを採用してフィラーの配向情報を抽出することを特徴とする請求項７に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項９】
前記抽出工程は、前記形状モデルとして複数の制御点によって定義されるスプライン曲線を中心軸とする仮想パイプを用い、前記候補領域内で前記制御点の座標をパラメータとしてランダムに変化させ、その仮想パイプと前記候補領域とのフィッティングの度合いを、前記仮想パイプに含まれるボクセルのデータ値に基づく評価値の積算値によって評価することを特徴とする請求項３乃至請求項６のいずれか一項に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項１０】
前記抽出工程は、前記スプライン曲線の前記複数の制御点について順次１つずつ座標をランダムに変化させ、その１つの制御点を変化させる毎に前記スプライン曲線を決定して前記フィッティングの度合いを評価し、その次の制御点を変化させる前に、前記決定されたスプライン曲線上で前記全ての制御点を等間隔に再配置することを特徴とする請求項８または請求項９に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項１１】
前記抽出工程は、前記形状モデル内に前記ボクセルのサイズ以下の間隔で評価点を設定し、その評価点にボクセルのデータ値に基づく評価値を付与し、前記評価点に付与した評価値の積算値によって前記フィッティングの度合いを評価することを特徴とする請求項７乃至請求項１０の何れか一項に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項１２】
前記抽出工程は、前記形状モデルの中心軸に垂直な面に同心円状に前記評価点を配置することを特徴とする請求項１１に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項１３】
前記評価点に付与する評価値は、前記ボクセルのデータ値を用いて線形補間して得られる値とすることを特徴とする請求項１１または請求項１２に記載の複合材料中繊維状フィラーの３次元画像処理方法。
【請求項１４】
前記抽出工程は、前記評価値から閾値を差し引いた値を新たな評価値として用いることを特徴とする請求項７乃至請求項１３の何れか一項に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項１５】
前記抽出工程は、前記新たな評価値が負の場合に、その値に所定の正数を掛けた値を新たな評価値として用いることを特徴とする請求項１４に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項１６】
前記抽出工程は、前記形状モデルの中心軸に近い位置の評価値ほど大きくなるように重み付けをして前記積算値を算出することを特徴とする請求項７乃至請求項１５の何れか一項に記載の複合材料中の繊維状フィラーの３次元画像処理方法。
【請求項１７】
母材に繊維状のフィラーを含有させた複合材料の３次元画像から前記フィラーの配向情報を抽出する複合材料中の繊維状フィラーの３次元画像処理装置において、
前記複合材料の３次元画像からフィラーを表す画素を含むと推定される候補領域を所定の閾値を参照することにより選出する領域選出手段と、
フィラーの形状モデルの形状と配置を変化させるパラメータをランダムに変化させるモンテカルロ法を用いて、前記領域選出手段によって選出された候補領域に前記形状モデルをフィッティングさせてフィラーの配向情報を抽出する抽出手段と、を備えることを特徴とする複合材料中の繊維状フィラーの３次元画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複合材料中の繊維状フィラーの３次元画像処理方法および３次元画像処理装置に関する。
【背景技術】
【０００２】
近年、樹脂の射出成形による電子部品は、小型化の進む電子機器用に小型精密化と薄肉化が進んでいる。このような電子部品を構成する樹脂射出による成形品は、その剛性を向上させたり線膨張係数を調整したりするために、樹脂母材にガラス繊維等の繊維状フィラーを混入させた複合材料を用いて成形される。成形品は、寸法精度を確保するために、成形後の反り変形を事前に予測して設計すること、および反り変形が発生しにくい条件で成形することが重要である。成形品の反り変形は、成形品の形状、成形品への熱履歴の分布、成形品における線膨張係数やヤング率などの熱的機械的特性の分布、および成形品への熱負荷などの要因が複合して発生する。反り変形は、熱負荷等に対する成形品の応答として発生し、成形品の熱的特性の分布および機械的特性の分布に依存する。これらの特性分布は成形品内の各点における繊維状フィラーの量の分布と繊維状フィラーの長手方向の向きの分布（配向分布）に依存する。従って、成形品中の繊維状フィラーの分布と配向分布とを把握することにより、これらの分布を最適化できるように成形条件や成形金型を最適化して、成形品の反り変形を抑制したり制御したりすることができる。【０００３】
成形品中の繊維状フィラーの配向分布を非破壊で評価する方法として、Ｘ線ＣＴ（コンピュータトモグラフィ）を用いて成形品の３次元画像を求め、その３次元画像を解析して繊維の形態や分布を求めて評価する方法が知られている（例えば、非特許文献１参照）。この評価方法は、互いに複雑につながった繊維同士が形成する３次元的な経路についての迂回度を指標に用いて、繊維配向の傾向を評価するものである。迂回度は、３次元画像から層状の矩形領域を切り出し、矩形領域の繊維の画像の細線化処理を行い、矩形領域の両端に両端点を有する経路について、両端点間の経路のうち最短の経路の経路長と両端点間の直線距離との比を求め、その比によって定義されている。【０００４】
また、半導体パッケージにおける樹脂で封止されたワイヤおよび生物体内の血管や神経系の３次元座標を非破壊的に直接的に抽出する方法が知られている（例えば、特許文献１参照）。この方法は、被検体の所定の方向と交差する複数の断層画像を求め、隣接する２つの画像の各々から被検体内の線状素の断面に対応するものとして求めた２つの構成点を、線状素を表すベクトルの始点および終点とする。断層画像はＸ線ＣＴまたはＭＲＩなどによって得られる。断層画像中の構成点は線状素の断面を表す円または楕円の画像の重心によって決定され、始点と終点の対は、隣接画像間で互いに最短距離となる２点が選択される。線状素のベクトルを、順次、互いに隣接する断層画像間で接続することにより、被検体に含まれる線状素の３次元座標が抽出される。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】中野亮  他著、「Ｘ線ＣＴによる繊維配向観察とシミュレーション」、成形加工、第２０巻、第４号、ｐｐ２３７−２４１（２００８年）
【特許文献】
【０００６】
【特許文献１】特開２０１２−３２２９３号公報
【発明の開示】
【０００７】
しかしながら、上述した非特許文献１に示されるような繊維状フィラーの配向分布の評価方法においては、画像に細線化処理を行うので、画像中のノイズ部分をフィラーとして抽出してしまう可能性がある。これは、細線化の処理によって繊維の太さの情報が失われてしまい、ノイズと繊維とを区別できなくなることによる。また、個々のフィラーを独立に抽出する方法ではないので、繊維配向の傾向という定性的な情報しか得られない。また、上述した特許文献１に示されるような線状素の３次元座標を抽出する方法においては、線状素の方向が被検体の所定の方向として既知であることが前提とされており、フィラーが任意の方向を向いて複雑に配向している場合には適用できない。
【０００８】
本発明は、上記課題を解消するものであって、母材に繊維状のフィラーを含む複合材料の３次元画像から個々のフィラーの配向情報を定量的に抽出できる複合材料中の繊維状フィラーの３次元画像処理方法および３次元画像処理装置を提供することを目的とする。
【０００９】
上記課題を達成するために、本発明の複合材料中の繊維状フィラーの３次元画像処理方法は、母材に繊維状のフィラーを含有させた複合材料の３次元画像から前記フィラーの配向情報を抽出する複合材料中の繊維状フィラーの３次元画像処理方法において、複合材料の３次元画像からフィラーを表す画素を含むと推定される候補領域を所定の閾値を参照することにより選出する領域選出工程と、フィラーの形状モデルの形状と配置を変化させるパラメータをランダムに変化させるモンテカルロ法を用いて、領域選出工程によって選出された候補領域に形状モデルをフィッティングさせてフィラーの配向情報を抽出する抽出工程と、を備えることを特徴とする。
【００１０】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、候補領域から１本のフィラーの配向情報を抽出した後に、抽出に関与した領域を除去し、除去後の候補領域に対してモンテカルロ法による抽出を繰り返してもよい。
【００１１】
この複合材料中の繊維状フィラーの３次元画像処理方法において、３次元画像をボクセルデータとして入力するデータ入力工程を備え、ボクセルデータは、複合材料のＸ線ＣＴによって取得され、ボクセルデータの各ボクセルは、Ｘ線の強度値に基づく値をデータ値として有してもよい。
【００１２】
この複合材料中の繊維状フィラーの３次元画像処理方法において、データ入力工程は、データ値を各ボクセル間で線形補間した値をデータ値とするボクセルを有する新たなボクセルデータを構成する補間工程を備えてもよい。
【００１３】
この複合材料中の繊維状フィラーの３次元画像処理方法において、領域選出工程は、ボクセルのデータ値とＸ線の強度値に基づく所定の閾値とを比較し、閾値よりも大きいデータ値を有するボクセルの集合を候補領域としてもよい。
【００１４】
この複合材料中の繊維状フィラーの３次元画像処理方法において、候補領域に属するボクセルから、互いに隣接するボクセルで成る群を生成し、生成された個々の群をそれぞれ新たな候補領域としてもよい。
【００１５】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、形状モデルとして仮想円柱を用い、候補領域内で仮想円柱の形状と配置を表すパラメータを変化させ、その仮想円柱と候補領域とのフィッティングの度合いを、仮想円柱に含まれるボクセルのデータ値に基づく評価値の積算値によって評価してもよい。
【００１６】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、仮想円柱によるフィッティングを行った後に、形状モデルとして複数の制御点によって定義されるスプライン曲線を中心軸とする仮想パイプを用い、候補領域内で制御点の座標をパラメータとしてランダムに変化させ、その仮想パイプと候補領域とのフィッティングの度合いを評価し、仮想円柱による評価よりも仮想パイプによる評価が所定割合以上改善する場合に、仮想パイプによるフィッティングを採用し、そうでない場合は仮想円柱によるフィッティングを採用してフィラーの配向情報を抽出してもよい。
【００１７】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、形状モデルとして複数の制御点によって定義されるスプライン曲線を中心軸とする仮想パイプを用い、候補領域内で制御点の座標をパラメータとしてランダムに変化させ、その仮想パイプと候補領域とのフィッティングの度合いを、仮想パイプに含まれるボクセルのデータ値に基づく評価値の積算値によって評価してもよい。
【００１８】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、スプライン曲線の複数の制御点について順次１つずつ座標をランダムに変化させ、その１つの制御点を変化させる毎にスプライン曲線を決定してフィッティングの度合いを評価し、その次の制御点を変化させる前に、決定されたスプライン曲線上で全ての制御点を等間隔に再配置してもよい。
【００１９】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、形状モデル内にボクセルのサイズ以下の間隔で評価点を設定し、その評価点にボクセルのデータ値に基づく評価値を付与し、評価点に付与した評価値の積算値によってフィッティングの度合いを評価してもよい。
【００２０】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、形状モデルの中心軸に垂直な面に同心円状に評価点を配置してもよい。
【００２１】
この複合材料中の繊維状フィラーの３次元画像処理方法において、評価点に付与する評価値は、ボクセルのデータ値を用いて線形補間して得られる値としてもよい。
【００２２】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、評価値から閾値を差し引いた値を新たな評価値として用いてもよい。
【００２３】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、新たな評価値が負の場合に、その値に所定の正数を掛けた値を新たな評価値として用いてもよい。
【００２４】
この複合材料中の繊維状フィラーの３次元画像処理方法において、抽出工程は、形状モデルの中心軸に近い位置の評価値ほど大きくなるように重み付けをして積算値を算出してもよい。
【００２５】
本発明の複合材料中の繊維状フィラーの３次元画像処理装置は、母材に繊維状のフィラーを含有させた複合材料の３次元画像からフィラーの配向情報を抽出する３次元画像処理装置において、複合材料の３次元画像からフィラーを表す画素を含むと推定される候補領域を所定の閾値を参照することにより選出する領域選出手段と、フィラーの形状モデルの形状と配置を変化させるパラメータをランダムに変化させるモンテカルロ法を用いて、領域選出手段によって選出された候補領域に形状モデルをフィッティングさせてフィラーの配向情報を抽出する抽出手段と、を備えることを特徴とする。
【００２６】
本発明の複合材料中の繊維状フィラーの３次元画像処理方法によれば、形状と配置を変化させるパラメータを有するフィラーの形状モデルとモンテカルロ法とに基づいて３次元画像から個々のフィラーを抽出するので、フィラーの配向情報を定量的に抽出できる。
【図面の簡単な説明】
【００２７】
【図１】本発明の一実施形態に係る複合材料中の繊維状フィラーの３次元画像処理方法についてのフローチャート。
【図２】同画像処理方法の変形例のフローチャート。
【図３】（ａ）は同画像処理方法の処理対象物の例を示す複合材料の斜視図、（ｂ）は（ａ）の複合材料の一部を示す斜視図。
【図４】図３（ｂ）に示した複合材料のＸ線ＣＴによる３次元画像を示す図。
【図５】（ａ）は複合材料中の繊維状フィラーの３次元画像を閾値処理した例を示す３次元画像を示す図、（ｂ）は（ａ）の画像の白黒反転画像を示す図。
【図６】同画像処理方法の他の変形例のフローチャート。
【図７】（ａ）（ｂ）は同変形例における補間工程を説明するボクセルの平面図。
【図８】同補間工程の他の例を説明するボクセルの立体図。
【図９】同画像処理方法の他の変形例のフローチャート。
【図１０】同変形例における分割工程によって生成された候補領域の例を示す３次元画像を示す図。
【図１１】（ａ）は同画像処理方法で用いられる仮想円柱の斜視図、（ｂ）は同仮想円柱をＸＹＺ座標空間に表現した斜視図、（ｃ）は同仮想円柱を極座標空間に表現した斜視図。
【図１２】同画像処理方法の抽出工程のフローチャート。
【図１３】同抽出工程における評価値の積算を説明するボクセルと形状モデルの断面図。
【図１４】（ａ）（ｂ）（ｃ）はそれぞれ異なる視点から見た候補領域とその領域中の仮想円柱の３次元画像を示す図。
【図１５】同画像処理方法の抽出工程の変形例のフローチャート。
【図１６】同画像処理方法の抽出工程の他の変形例のフローチャート。
【図１７】同画像処理方法の抽出工程の他の変形例のフローチャート。
【図１８】（ａ）は同画像処理方法の抽出工程の他の変形例に関し、仮想円柱の端点の位置を変動させる様子を示す斜視図、（ｂ）は端点を変動させた後の仮想円柱の斜視図、（ｃ）は（ｂ）の仮想円柱に設定する評価点を説明する平面図。
【図１９】同仮想円柱に評価点を設定した状態の３次元画像を示す図。
【図２０】同仮想円柱と評価点とを用いる抽出工程のフローチャート。
【図２１】複合材料中に曲がった繊維状フィラーが含まれる場合の３次元画像を示す図。
【図２２】曲がった繊維状フィラーに対して仮想円柱モデルをフィッティングした様子を模式的に示す斜視図。
【図２３】同画像処理方法の抽出工程の他の変形例のフローチャート。
【図２４】（ａ）は仮想パイプを制御点とスプライン曲線とによって示す斜視図、（ｂ）は同制御点の１つを変動させた場合の同スプライン曲線の斜視図。
【図２５】（ａ）は同仮想パイプのスプライン曲線に沿って評価点設定用の円盤を配置した斜視図、（ｂ）は同円盤に評価点を配置した平面図。
【図２６】仮想パイプに評価点を設定した状態の３次元画像を示す図。
【図２７】仮想パイプを用いたフィッティング途中の３次元画像を示す図。
【図２８】（ａ）（ｂ）（ｃ）はそれぞれ同仮想パイプの各制御点を変動させる互いに異なる手順の例を示す概念図。
【図２９】同画像処理方法の他の変形例を示すフローチャート。
【図３０】（ａ）〜（ｈ）はそれぞれ同仮想パイプを用いて各制御点を変動させてフィッティングを行う手順を説明する概念図。
【図３１】一実施形態に係る複合材料中の繊維状フィラーの３次元画像処理装置のブロック構成図。
【図３２】同画像処理方法および同画像処理装置によって抽出された繊維状フィラーを３次元画像に重ねて表示した３次元画像を示す図。
【図３３】同画像処理方法および同画像処理装置によって抽出された繊維状フィラーの繊維長の度数割合を示す度数分布図。
【発明を実施するための形態】
【００２８】
以下、本発明の一実施形態に係る複合材料中の繊維状フィラーの３次元画像処理方法および３次元画像処理装置について、図面を参照して説明する。図１は、一実施形態に係る複合材料中の繊維状フィラーの３次元画像処理方法（以下、画像処理方法という）を示す。本画像処理方法は、図１に示すように、領域選出工程（＃１）と、抽出工程（＃２）とを備え、これらの工程を実行して、母材に繊維状フィラーを含む複合材料の３次元画像から複合材料中のフィラーの配向情報を抽出する。配向情報は、例えば、フィラーが、どのような場所で、どのような形状と向きで配置されているかという情報である。領域選出工程（＃１）は、複合材料の３次元画像を与えられ、その３次元画像からフィラーを表す画素（立体画素）を含むと推定される候補領域を所定の閾値を参照することにより選出する。抽出工程（＃２）は、決定されたフィラーの形状モデルの形状と配置を規定するパラメータをランダムに変化させるモンテカルロ法を用いて、領域選出工程によって選出された候補領域に形状モデルをフィッティングさせてフィラーの配向情報を抽出する。フィラーの配向情報は、フィッティング結果の形状モデルの形状と配置から抽出される。処理対象の複合材料は、例えば、液晶ポリマーから成る母材に、繊維状のフィラーとして強化用のガラス繊維を含有させたガラス繊維強化樹脂である。この複合材料の場合、本画像処理方法は、液晶ポリマー樹脂中のガラス繊維の配向情報を抽出する。
【００２９】
図２は、上述の画像処理方法の変形例を示す。この画像処理方法は、上述の図１の実施形態において、領域選出工程（＃１）の前工程としてデータ入力工程（＃０）を備え、抽出工程（＃２）が、モデル設定工程（Ｓ１）と、抽出メイン工程（Ｓ２）と、反復処理工程（Ｓ３，Ｓ４）とを備えている。
【００３０】
データ入力工程（＃０）は、複合材料中の繊維状フィラーの３次元画像をボクセルデータとして入力する。ボクセルデータは、複合材料のＸ線ＣＴによって取得され、ボクセルデータの各ボクセルは、Ｘ線の強度値に基づく値をデータ値として有する。例えば、図３（ａ）（ｂ）に示すように、複合材料Ｃは、長い箱状に射出成形されたガラス繊維強化された樹脂成形品１の一部である。このような樹脂成形品１において、矢印ａで示す長手方向に樹脂注入が行われると、ガラス繊維は、大略その矢印ａの方向に沿って配向する。
【００３１】
図４は、図３（ｂ）に示した複合材料ＣのＸ線ＣＴによる３次元画像Ｇ１を示す。３次元画像Ｇ１は、複合材料Ｃに様々な方向からＸ線を透過させ、Ｘ線の吸収量の分布に応じて形成されるプロファイルを収集し、それらのプロファイルから複合材料Ｃの多数の断層像を白黒の濃淡画像として再構成して３次元表示したものである。母材（例えば、液晶ポリマー）と繊維状フィラー（例えば、ガラス繊維）とのＸ線の吸収量の違いにより、繊維状フィラーの３次元画像情報を含む３次元画像Ｇ１が得られる。ボクセルのデータ値は、フィラーの存在率を表す。ボクセルのデータ値は、例えば、繊維状フィラーを表示するために濃淡や白黒の何れを用いるかに応じて、また、データ値を使用する目的に応じて、任意に設定することができる。ボクセルの用語は、本来、３次元画像空間を構成する一般的な立体画素の意味であるが、本画像処理方法においては、ある閾値で選択された画像処理の対象となる立体画素の意味で、例えば、候補領域からボクセルを削除する、等の表現で限定的に用いることがある。
【００３２】
領域選出工程（＃１）は、複合材料Ｃの３次元画像中にフィラーを表す画素を含む領域であると推定される領域を含む候補領域を選出する。領域選出工程（＃１）は、ボクセルのデータ値を候補領域選出用の閾値によって判別し、例えば、その閾値よりも大きいデータ値を有するボクセルの集合を候補領域とする。この場合、データ値はＸ線の強度値であり、その閾値は、Ｘ線の強度値の分布に基づいて設定される。図５（ａ）（ｂ）の３次元画像Ｇ２，Ｇ３は、それぞれこのような閾値処理によって候補領域を選出した例を示す。この領域選出工程（＃１）によって、フィラーを表す画素の領域がバックグラウンド領域から切り出される。
【００３３】
抽出工程（＃２）のモデル設定工程（Ｓ１）は、領域選出工程（＃１）によって設定された候補領域に対し、形状と配置をそれぞれ変化させるパラメータを有してフィラーの形状を再現できるように構成された形状モデルを設定する。形状モデルの形状は、複合材料中に含有されるフィラーの形状に応じて設定される。例えば、円柱、長楕円、角柱、複数の円柱が曲線に沿って連続した曲線チューブなどが形状モデルとされる。いずれも長手方向が定義される配向性を有するものである。
【００３４】
形状モデルの形状のパラメータは、形状モデルの形と大きさを決めるものであり、例えば、円柱の場合は半径と長さ、長楕円の場合は３次元楕円の３軸長、角柱の場合は断面形状と大きさ、および角柱の長さ、などを、所定の範囲内で変化させるように設定する。曲線チューブ（パイプ）の場合は、半径、曲線長、曲線の曲がり具合（例えば、折れ線近似して、各線分長、相互の折れ角度の方向と大きさ等）等を、所定の範囲内で変化させるように、形状のパラメータを設定すればよい。また、形状モデルの配置のパラメータは、３次元空間座標軸が設定された空間における候補領域に対して、形状が定まった形状モデルの位置と姿勢を定めるものであり、例えば、形状モデルの長手方向の両端の位置座標を配置のパラメータとすることができる。配置のパラメータは、形状モデルの配置を一意的に決定できるものであればよく、形状モデルの中心位置座標と方位角などを用いてもよい。
【００３５】
抽出メイン工程（Ｓ２）は、モデル設定工程（Ｓ１）によって設定された形状モデルのパラメータを変化させて候補領域に形状モデルをフィッティングさせ、そのフィッティングした形状モデルの形状と配置の情報をフィラーの配向情報とする。抽出メイン工程（Ｓ２）は、形状モデルのパラメータをランダムに変化させてフィッティングさせる処理をモンテカルロ法に基づいて行う。モンテカルロ法では、フィラーが存在する確率が高いボクセル群に対し、乱数を用いてパラメータを一定範囲内で変動させ、より評価の高くなるパラメータを採用するように、形状モデルのパラメータが最適化される。モンテカルロ法において、パラメータを変化させる毎に候補領域と形状モデルのフィッティングの度合いが所定の方法で求められた全体評価値によって評価される。その全体評価値が改善する方向に更新されなくなった時点で、そのときの形状モデルをフィラーと見倣して、１本のフィラーの抽出処理が収束する。
【００３６】
反復のための処理工程（Ｓ３）は、候補領域から抽出メイン工程（Ｓ２）によって１本のフィラーの配向情報を抽出する処理を繰り返して、複数のフィラーが含まれる候補領域の全域から全てのフィラーの配向情報を抽出するための処理を行う。反復のための処理工程（Ｓ３）は、１本のフィラーが抽出された後に、後続する次回の抽出を効率的に行うために、１つの形状モデルの最終のフィッティングに関与した領域を候補領域から除去する。その領域を除去したあとに、候補領域の残存領域があれば（Ｓ４でＹｅｓ）、モデル設定工程（Ｓ１）からの処理が繰り返され、なければ（Ｓ４でＮｏ）、画像処理が終了する。
【００３７】
本実施形態の画像処理方法によれば、形状と配置を変化させるパラメータを有するフィラーの形状モデルとモンテカルロ法とに基づいて３次元画像から個々のフィラーを抽出するので、フィラーの配向情報を定量的に抽出できる。また、本画像処理方法は、複合材料中の繊維状フィラーの３次元画像が、Ｘ線ＣＴなどによって非破壊で得られることにより、非破壊でフィラーの配向情報を定量的に抽出できる。
【００３８】
図６は、上述した画像処理方法の他の変形例を示す。この変形例は、複合材料中の繊維状フィラーの３次元画像処理方法において、データ入力工程（＃０）が、ボクセルのデータ値を各ボクセル間で線形補間する補間工程（＃０ａ）を備えるものである。その線形補間は、例えば、三重線形補間を用いて行うことができる。
【００３９】
図７（ａ）に示すように、３次元画像データには、データ値を有するボクセルＢｘと、測定条件の影響などにより、本来有しているべきデータ値を有さないボクセルＢ０がある。補間工程（＃０ａ）は、図７（ｂ）に示すように、このようなボクセルＢ０に、その周辺のボクセルＢｘのデータ値に基づくデータ値を付与する。これにより、フィッティングの精度を向上できる。また、補間工程（＃０ａ）は、図８に示すように、各ボクセルＢｘ間でデータ値を線形補間して得た値をデータ値とする新たなボクセルｂｘを有する新たなボクセルデータを構成することができる。すなわち、データ入力工程（＃０）は、ボクセルＢｘを分割して小さなボクセルｂｘを生成し、ボクセルＢｘ間で線形補間した値をボクセルｂｘのデータ値とする。このボクセルの分割補間は、ボクセルのサイズが、フィラーのサイズと同程度か、より大きい場合などに、フィッティングの精度を向上できる。
【００４０】
図９、図１０は、上述した画像処理方法のさらに他の変形例を示す。この変形例は、図９に示すように、上述の領域選出工程（＃１）が、大きな候補領域を、より小さな候補領域に分割する分割工程（＃１ａ）を備えるものである。分割工程（＃１ａ）では、候補領域に属するボクセルから、互いに隣接するボクセルで成る群（グループ）を生成し、個々の群をそれぞれ新たな候補領域とする。ボクセルは、２つのボクセルが互いに隣接する状況として、面間隣接、辺間隣接、および、頂点間隣接がある。そこで、例えば、互いに面間隣接するボクセルの集合であって所定個数以上のボクセルからなる集合を群と認定し、その群を新たな個別の候補領域とする。面間隣接したボクセルであっても、それらの集合が所定個数以上のボクセルを含まない場合、および、面間隣接していない孤立したボクセルは、候補領域から削除する。
【００４１】
もともと閾値より大きいデータ値を有するボクセルの１つの集合として選出されていた候補領域は、この分割工程（＃１ａ）を経て複数の候補領域に細分化される。さらに、フィラーを構成するとは認められない所定個数以下に離散した島状のボクセルは排除される。その所定個数は、フィラーの形状その他の事前の知見、または試験的に抽出メイン工程（Ｓ２）を実行して得られる知見などに基づいて設定される。図１０の３次元画像Ｇ４は、分割工程（＃１ａ）を経て設定された候補領域の例を示す。なお、これらの分割工程の処理に際し、面間隣接だけでなく、辺間隣接、さらには頂点間隣接したボクセルも群に含めるようにしてもよい。
【００４２】
本画像処理方法は、このような分割工程（＃１ａ）を備えることにより、モンテカルロ法による探索の範囲を狭くできるので、効率よく画像処理でき、抽出メイン工程（Ｓ２）の計算負荷を軽減することができる。また、本方法は、隣接するボクセルによって候補領域を設定するので、座標によるグループ分けをする場合における、グループ間をまたがるフィラーが分断されて２本のフィラーとして判定されてしまう、という弊害を回避できる。分割工程（＃１ａ）は、上述の補間工程（＃０ａ）を経たデータに適用できる。
【００４３】
図１１乃至図１４は、上述した画像処理方法のさらに他の変形例を示す。この変形例は、図１１（ａ）（ｂ）（ｃ）に示すように、フィラーの形状を再現する形状モデルとして仮想円柱Ｍを用いるものである。仮想円柱Ｍは、その半径Ｒと長さＬとによって形状が決定され、その中心軸上の両端点Ｐ１，Ｐ２の座標Ｐ１（ｘ１，ｙ１，ｚ１），Ｐ２（ｘ２，ｙ２，ｚ２）によって、ＸＹＺ直交座標空間における空間配置が決定されるので、８つのパラメータを有する。両端点Ｐ１，Ｐ２間の距離が長さＬを与えるので、パラメータを７つにできる。フィラーの直径が既知で、一定の値に固定できる場合は、半径Ｒを固定して、６つのパラメータとすることができる。仮想円柱Ｍの配置は、両端点Ｐ１，Ｐ２の座標による表示に限らず、中心点Ｐｃの座標を用いたり、極座標（ｒ，θ，φ）を用いて仮想円柱の傾きを表示したりすることができる。
【００４４】
このように、仮想円柱のパラメータとして、両端点の３次元座標、中心の３次元座標、方向（方位、高度）、長さ、半径などがあり、これらのうち、独立しているものは一部であり、残りは他の独立変数から導かれる。独立変数と従属変数の関係は相補的なものであり、どれを独立変数として扱うかは、ソフトウェア内で柔軟に変更することができる。例えば、端点の位置を変動する場合は、中心位置、長さ、および方向は端点の位置から導かれる。長さを変動する場合は、端点の位置は中心位置と長さから導かれる。フィラーの直径が既知である場合、円柱のパラメータとして、両端点の３次元座標（３×２＝６自由度）のみを最適化すればよい。
【００４５】
この画像処理方法は、図１２に示すように、抽出工程（＃２）のモデル設定工程（Ｓ１）においてフィラーの形状を再現する形状モデルとして仮想円柱を設定し、抽出メイン工程（Ｓ２）において仮想円柱を用いてフィラーの抽出を行う。抽出メイン工程（Ｓ２）では、候補領域内で仮想円柱Ｍの形状（半径、長さ）と配置（両端の座標等）のパラメータをランダムに変化させ、変化した形状と配置を有する仮想円柱と候補領域とのフィッティングの度合いを評価してフィラーが抽出される。フィッティングの度合いは、図１３に示すように、各ボクセルＢｘの内、仮想円柱Ｍに含まれるボクセルＢｘ（黒丸で示されている）のデータ値を評価値として、その評価値の積算値によって評価される。評価は、例えば、その積算値、仮想円柱Ｍに含まれるボクセルＢｘの個数Ｎ、積算値を個数Ｎで除算して規格化した値などを用いて、総合的に行われる。例えば、規格化した値が一定で個数Ｎが増大傾向にある場合は、フィッティングがまだ収束してなく、仮想円柱Ｍをさらに伸張できる、と判断される。積算値が最大となる仮想円柱Ｍの形状と配置が、フィラーの配向情報として抽出される。総合評価値、例えば、図１４（ａ）（ｂ）（ｃ）の３次元画像Ｇ５は、それぞれ異なる視点から見た候補領域の３次元画像とその領域中でフィッティングされた１つの仮想円柱Ｍを示す。
【００４６】
上述の仮想円柱Ｍを用いるフィラーの抽出処理をより詳細に説明する。フィラーの３次元画像、例えば、上述の３次元画像Ｇ４は、複数の候補領域の群を含み、各群は、それぞれ複数本のフィラーの画像を含んでいると考えられる。そこで、モデル設定工程（Ｓ１）において、それらのうちの１つの群について、複数の仮想円柱Ｍをランダムに生成し、各仮想円柱Ｍについて適合度を評価し、最も評価の高い仮想円柱Ｍを選択し、そのパラメータを初期パラメータとして採用する。次に、抽出メイン工程（Ｓ２）において、その初期パラメータから開始して、所定の変動範囲内で、パラメータをランダムに動かして評価し、評価が改善したら、パラメータを更新する。すなわち、最適化の基本アルゴリズムにはモンテカルロ法を用い、乱数を用いてパラメータを一定範囲内で変動させ、適合度の評価がより高くなるパラメータを採用する。この手順を、選択された仮想円柱について反復し、評価が改善しなくなったら、収束したとしてパラメータを確定し、これにより、１本のフィラーが抽出される。
【００４７】
次に、反復のための処理工程（Ｓ３）において、その抽出されたフィラーに関するボクセル、すなわち仮想円柱Ｍに属するボクセルは、候補領域から削除し、残りの候補領域について、同一群内で再度、初期パラメータを採用する手順から以後の手順を繰り返す。ボクセルが仮想円柱Ｍに属するか否かは、例えば、ボクセルの中心や１頂点などの、ボクセルに予め設定した点が仮想円柱Ｍ内に含まれているか否かで決めればよい。１つの群内でフィラーを抽出できなくなったら、他の群についてフィラーの抽出を行う。全ての群についてフィラーの抽出が終わった時点で、抽出工程（＃２）が終了する。なお、１本のフィラーを抽出した後に候補領域からボクセルを削除する際に、所定の複数個数以下の島状のボクセルの組が発生したら、その島状のボクセルを同時に削除し、その後の処理の計算負荷を減らす。また、群内で、処理対象のボクセルの数が所定の個数以下になったとき、その群内でフィラーを抽出できなくなった、すなわち、群内で抽出が完了したと判断する。
【００４８】
図１５は、上述した画像処理方法のさらに他の変形例を示す。この画像処理方法の抽出工程（＃２）は、上述のモデル設定工程（Ｓ１）の後の抽出メイン工程（Ｓ２）が、評価点設定工程（Ｓ２１）を備えるものである。評価点設定工程（Ｓ２１）は、仮想円柱Ｍ内にボクセルのサイズ以下の間隔で多数の評価点（サンプリング点）を設定し、各評価点にボクセルのデータ値に基づく評価値を付与する。
【００４９】
評価点は、例えば、上述の図１３に示した仮想円柱Ｍ内のボクセルＢｘを上述の図８に示したように細分化して得た小さなボクセルの代表点、および、元の大きなボクセルの代表点とすればよい。各評価点の評価値は、細分化で生成されたボクセルを含む元の大きなボクセルのデータ値である。このような評価点と評価値は、仮想円柱Ｍ内だけで生成される。これにより、ボクセルデータの全体でデータ密度を上げる場合（例えば、図８）に比べて、計算を局所的に済ますことができる利点がある。評価点は、ボクセルの位置とデータ値とで構成されるボクセルデータ空間に固定された点であり、仮想円柱Ｍは、そのパラメータの変動に伴って、ボクセルデータ空間内で位置を変え、伸縮し、取り込んだ評価点によって、評価される。
【００５０】
抽出メイン工程（Ｓ２）は、上述の図１２の抽出メイン工程（Ｓ２）におけるボクセルのデータ値の積算値に代えて、各評価点に付与された評価値の積算値を用いて、フィッティングの度合いを評価する。このようなボクセルの細分化による評価点とその評価値を用いることにより、仮想円柱Ｍ内に含まれるボクセルの体積をより精度良く積算値に取り込んで、フィッティングの度合いを、より精度良く評価することができる。
【００５１】
図１６は、上述した画像処理方法のさらに他の変形例を示す。この画像処理方法の抽出工程（＃２）は、上述の評価点設定工程（Ｓ２１）の後に、評価値補間工程（Ｓ２２）を備えるものである。すなわち、評価値補間工程（Ｓ２２）は、評価点に付与する評価値として、元のボクセルのデータ値を用いて線形補間して得られる値を、評価点の評価値とする。この場合、評価点は、ボクセルの代表点およびボクセルの代表点間に内挿した点であり、内挿による評価点の評価値はボクセルのデータ値を用いた線形補間で与えられる。線形補間の方法として、三重線形補間を用いることができる。例えば、各ボクセル間の中間位置に評価点を設定する場合、その評価点を囲む８つのボクセルのデータ値（評価値）の三重線形補間によって、評価点の評価値を決定すればよい。このような評価値補間工程（Ｓ２２）を備えることにより、フィッティングの度合いを、より精度良く評価することができる。
【００５２】
図１７は、上述した画像処理方法のさらに他の変形例を示す。この画像処理方法の抽出工程（＃２）は、上述の評価値補間工程（Ｓ２２）の後に、評価値減算工程（Ｓ２３）、ペナルティ工程（Ｓ２４）、重み付け工程（Ｓ２５）をさらに備えている。評価値減算工程（Ｓ２３）は、評価値から所定の閾値を差し引いた値を新たに評価値とする。ペナルティ工程（Ｓ２４）は、評価値が負の場合に、その値に所定の正数を掛けた値を新たな評価値とする。例えば、評価値が負の値を取る場合に評価値を１０倍した値を新たな評価値とする。これらにより、仮想円柱がフィラー領域から出ることを抑止でき、フィラーの抽出効率を向上できる。重み付け工程（Ｓ２５）は、仮想円柱の中心軸に近い位置の評価値ほど大きくなるように重み付けをする、例えば、中心軸からの距離の逆数に比例するように重み付けをする工程である。このような重み付けにより、仮想円柱の中心軸が、評価値の高い部分に早く近づくように、処理を加速することができる。
【００５３】
図１８、図１９、図２０は、画像処理方法のさらに他の変形例を示す。この画像処理方法は、上述した方法におけるボクセルデータ空間に固定された評価点を用いることに替えて、形状モデル（仮想円柱Ｍ）の位置、姿勢、形状の変化に合わせて位置が変化する評価点を用いるものである。図１８（ａ）に示すように、パラメータの変動により、長さＬ０の仮想円柱Ｍの端点Ｐ１が位置Ｐｘに移動した場合を考える。フィッティングの評価は、図１８（ｂ）に示す、端点Ｐ１，Ｐ２を有する長さＬの傾いた仮想円柱Ｍについて行われる。端点Ｐ１は、元の位置Ｐ０から位置Ｐｘに移動している。図１８（ｂ）（ｃ）に示すように、評価点ｂは、仮想円柱Ｍ（形状モデル）の中心軸に垂直な円盤Ｂを端点Ｐ１，Ｐ２間に一定の間隔Δ１で設定し、各円盤Ｂ上に同心円状に配置設定される。円盤Ｂ上の評価点ｂは、半径Ｒ方向に一定の間隔Δ２、円周方向に一定の間隔Δ３で配置される。図１９は、仮想円柱Ｍを、高さ、半径、円周方向のそれぞれについて細かく分割し、各分割点に評価点ｂを設定した様子を示す。画像Ｇ６中の仮想円柱Ｍの周辺のグレー部分は、閾値以上のデータ値を有して候補領域を形成しているボクセルを表す。
【００５４】
図２０を参照して、この画像処理方法の抽出工程（＃２）を説明する。モデル設定工程（Ｓ１）では、仮想円柱Ｍの端点Ｐ１，Ｐ２の位置および半径Ｒの初期値が設定される。抽出メイン工程（Ｓ２）は、単一の直線上のフィラーを抽出するため、工程（Ｓ１０１〜Ｓ１１０）を繰り返す。モデル設定工程（Ｓ１）と抽出メイン工程（Ｓ２）とを合わせて、単一直線工程（＃２１）と称することにする。
【００５５】
抽出メイン工程（Ｓ２）の最初に、変数Ｐｉに端点Ｐ１を設定し（Ｓ１０１）、端点Ｐ１の位置を変動させる（Ｓ１０２）。端点Ｐ２と変動後の端点Ｐ１の間に、等間隔に円盤Ｂを配置し（Ｓ１０３）、各円盤Ｂに同心円状に評価点ｂ（サンプリング点）を配置する（Ｓ１０４）。円盤Ｂは、端点Ｐ１，Ｐ２にも配置される。評価点ｂは、仮想円柱Ｍに対して設定されるため、ボクセルデータ空間に対して非固定であり、任意の配置となる。各評価点ｂに、ボクセルのデータ値に基づいて３重線形補間で求めた評価値を付与する（Ｓ１０５）。仮想円柱Ｍ内の全ての評価点ｂの評価値を積算し、その積算値によってフィッティングの度合いが評価される（Ｓ１０６）。
【００５６】
端点Ｐ１の変動に対する評価の結果、モンテカルロ法における処理が収束していなければ（Ｓ１０７でＮｏ）、工程（Ｓ１０２）からの処理が繰り返される。処理が収束ならば（Ｓ１０７でＹｅｓ）、工程（Ｓ１０８）を経由して、変数Ｐｉに端点Ｐ２を設定し（Ｓ１０９）、端点Ｐ２について、端点Ｐ１の場合と同様に工程（Ｓ１０２〜Ｓ１０７）が行われる。両端点Ｐ１，Ｐ２の各々についての収束が確認されると（Ｓ１０８でＹｅｓ）、収束判断のために、工程（Ｓ１０１）からの処理が少なくとも１回強制的に繰り返される。両端点Ｐ１，Ｐ２の変動に対する評価の向上が収束したと判断されると（Ｓ１１０でＹｅｓ）、抽出メイン工程（Ｓ２）が終了する。これにより、単一直線工程（＃２１）が終了する。工程（Ｓ３、Ｓ４）は上記同様である。
【００５７】
上述の円盤Ｂ上の同心円状に配置した評価点ｂを用いることにより、形状モデルすなわち仮想円柱Ｍに含まれる情報を、仮想円柱Ｍの外形境界部分において取りこぼすことなく有効利用できる。上述の工程（Ｓ１０６）におけるフィッティングの度合いの評価には、上述した図１７における評価値減算工程（Ｓ２３）、ペナルティ工程（Ｓ２４）、重み付け工程（Ｓ２５）の手法を用いることができる。
【００５８】
図２１乃至図２８は、画像処理方法のさらに他の変形例を示す。この画像処理方法は、湾曲したフィラーに対応できる形状モデルを用いるものである。図２１の画像Ｇ７に示すように、複合材料中の繊維状フィラーの３次元画像中に湾曲したフィラーＦが発生することがある。これは、複合材料の強度向上のために繊維状のフィラーの長繊維化が図られる結果、フィラーが、より湾曲しやすいことによる。図２２の画像Ｇ８に示すように、このような湾曲したフィラーについて仮想円柱Ｍによるフィッティングを行うと、複数の短い仮想円柱Ｍに分断された結果が得られる。そこで、湾曲したフィラーに対応できる形状モデルとして、複数の制御点によって定義されるスプライン曲線を中心軸とする一定半径の仮想パイプ（チューブ）を用いる。
【００５９】
図２３に示すように、この画像処理方法の抽出工程（＃２）は、モデル設定工程（Ｓ１）と抽出メイン工程（Ｓ２）とを合わせた単一曲線工程（＃２２）と、その後の工程（Ｓ３，Ｓ４）とを備えている。モデル設定工程（Ｓ１）では、端点Ｐ１，Ｐ２間に、ｎ個の制御点｛Ｃｉ：ｉ＝１，・・，ｎ｝，Ｃ１＝Ｐ１，Ｃｎ＝Ｐ２を等間隔に配置し、初期値として直線状のスプライン曲線が設定される。このスプライン曲線に半径Ｒを付加して、スプライン曲線を中心軸とする一定半径Ｒの仮想パイプＭＰが設定される。スプライン曲線は、制御点間を任意の多項式で接続した曲線であり、従って、その曲線は必ず制御点を通る曲線または直線である。
【００６０】
抽出メイン工程（Ｓ２）では、全制御点Ｃｉがスプライン曲線上に等間隔に再配置され（Ｓ２００）、変数ｉ＝１とされ（Ｓ２０１）、端点Ｐ１＝Ｃ１の位置が変動される（Ｓ２０２）。図２４（ａ）に、スプライン曲線Ｓｐと制御点｛Ｃｉ：ｉ＝１，・・，ｎ｝を備えた仮想パイプＭＰが示されている。制御点Ｃ１（端点Ｐ１）の位置が移動される様子が示されている。なお、図２４（ｂ）は、制御点Ｃ１が移動した後、工程（Ｓ２００）によって、スプライン曲線Ｓｐ上で全制御点Ｃｉが等間隔に再配置される様子を示す。
【００６１】
上述の工程（Ｓ２０２）による制御点Ｃ１（一般化して制御点Ｃｉ）の位置変動後に、図２５（ａ）に示すように、スプライン曲線Ｓｐすなわち仮想パイプの中心軸に沿って、その中心軸に垂直に、等間隔に、円盤Ｂが配置される（Ｓ２０３）。各円盤Ｂには、図２５（ｂ）に示すように、同心円状に評価点ｂが配置される（Ｓ２０４）。この円盤Ｂの配置間隔の設定や円盤Ｂ上の評価点の配置は、上述した図１８（ｂ）（ｃ）に示したΔ１，Δ２，Δ３と同様の値を用いて行えばよい。評価点ｂは、仮想パイプＭＰに対して設定されるため、ボクセルデータ空間に対して非固定であり、任意の配置となる。
【００６２】
各評価点ｂには、ボクセルのデータ値に基づいて３重線形補間で求めた評価値が付与される（Ｓ２０５）。仮想パイプＭＰ内の全ての評価点ｂの評価値を積算し、その積算値によってフィッティングの度合いが評価される（Ｓ２０６）。単一の制御点Ｃｉの変動に対する評価の結果、モンテカルロ法における処理が収束していなければ（Ｓ２０７でＮｏ）、工程（Ｓ２０２）からの処理が繰り返される。処理が収束ならば（Ｓ２０７でＹｅｓ）、工程（Ｓ２０８）を経由して、変数ｉについて、ｉ＝ｉ＋１とされ（Ｓ２０９）、次の制御点Ｃｉについて、上記同様に工程（Ｓ２０２〜Ｓ２０７）が行われる。
【００６３】
全制御点Ｃ１〜Ｃｎの各々についての収束が確認されると（Ｓ１０８でＹｅｓ）、収束判断のために、工程（Ｓ２０１）からの処理が少なくとも１回強制的に繰り返される。両端点Ｐ１，Ｐ２の変動に対する評価の向上が収束したと判断されると（Ｓ２１０でＹｅｓ）、抽出メイン工程（Ｓ２）が終了する。これにより、単一曲線工程（＃２２）が終了する。工程（Ｓ３、Ｓ４）は上記同様である。図２６の画像Ｇ９は、仮想パイプＭＰに評価点ｂを設定した状態を示す。図２７の画像Ｇ１０は、仮想パイプＭＰを用いたフィッティング途中の状態を示す。
【００６４】
上記において、制御点Ｃｉの数ｎを可変としてもよい。例えば、最初はｎ＝２、すなわち仮想円柱Ｍとし、制御点Ｃｉを際配置する工程（Ｓ２００〜Ｓ２０８）のタイミングで、ｎを増やせばよい。また、上記において、制御点｛Ｃｉ：ｉ＝１，・・，ｎ｝は、図２８（ａ）に示すように、その並びの順番に、制御点Ｃ１から位置の変動を行い、全体の繰り返しの際も、同じく制御点Ｃ１から位置の変動を行っている。これに対し、図２８（ｂ）に示すように、制御点Ｃ１から位置の変動を行い、全体の繰り返しの際には、逆順で制御点Ｃｎから位置の変動を行うようにしてもよい。また、図２８（ｃ）に示すように、制御点｛Ｃｉ：ｉ＝１，・・，ｎ｝の両端から、交互に、制御点Ｃｉの位置の変動を行うようにしてもよい。制御点Ｃｉの位置の変動の順番は、例えば、変動量の大きいところを先に変動させることで早期に収束させるために、任意に選択することができる。
【００６５】
図２９、図３０は、画像処理方法のさらに他の変形例を示す。この画像処理方法は、図２９に示すように、上述の図２０に示した単一直線工程（＃２１）と、図２３に示した単一曲線工程（＃２２）と、その後の新たな判定工程（＃２３）と、を組み合わせた処理方法である。すなわち、この画像処理方法は、仮想円柱Ｍによる単一のフィラーに対するフィッティングを行い（＃２１）、決定された仮想円柱Ｍの両端点Ｐ１，Ｐ２を端点とする仮想パイプＭＰを用いて単一のフィラーに対するフィッティングを行う（＃２２）。その後、判定工程（＃２３）において、単一直線工程（＃２１）による評価結果と、単一曲線工程（＃２２）による評価結果とを比較する（Ｓ３０）。比較の結果、仮想パイプＭＰによる評価が所定割合、例えば１０％以上、改善している場合に（Ｓ３０でＹｅｓ）、仮想パイプＭＰによるフィッティング結果を採用し（Ｓ３１）、そうでない場合は（Ｓ３０でＮｏ）、仮想円柱Ｍによる結果を採用する（Ｓ３２）。
【００６６】
図３０（ａ）〜（ｈ）を参照して、上述の処理の具体例を説明する。図３０（ａ）に示すように、直線状の形状モデルである仮想円柱Ｍによってフィッティングを行うと４つの仮想円柱Ｍでフィッティングされるような、湾曲した１本のフィラーを想定する。湾曲したフィラーの全長に対して、単一直線工程（＃２１）によって、短い仮想円柱Ｍが抽出されている。図３０（ｂ）に示すように、抽出された仮想円柱Ｍに対して、単一曲線工程（＃２２）により、直線状のスプライン曲線Ｓｐが設定され、端点Ｐ１に位置する制御点Ｃ１の位置が変動される。図３０（ｃ）に示すように、制御点Ｃ１の位置が変動される毎に、制御点Ｃｉの再配置とフィッティングの評価が行われ、さらに、制御点Ｃ１の位置が変動される。図３０（ｄ）に示すように、制御点Ｃ１の位置がフィラーの端部に至ると、制御点Ｃ１の位置の変動が収束する。
【００６７】
その後、図３０（ｅ）〜（ｇ）に示すように、制御点Ｃ２，Ｃ３，・・，Ｃｎと順番に位置の変動と評価が行われ、制御点Ｃｎの位置がフィラーの他端部に至ると、制御点Ｃｎの位置の変動が収束する。この後、再度、制御点Ｃ１から順に、制御点Ｃｎまで位置の変動と評価が行われ、評価の向上の変化具合に基づいて、１本のフィラーに対する抽出の収束が確認される。その結果、図３０（ｈ）に示すように、湾曲した１本のフィラーに対する仮想パイプＭＰが決定される。
【００６８】
図３１は、一実施形態に係る複合材料中の繊維状フィラーの３次元画像処理装置（以下、画像処理装置２という）を示す。画像処理装置２は領域選出部１１、モデル設定部１２、抽出メイン部１３、データ入力部１４、操作部１５、表示部１６、および、これらを制御する制御部１０とを備え、母材に繊維状のフィラーを含有した複合材料の３次元画像からフィラーの配向情報を抽出する。３次元画像は、データ入力部１４から入力される。制御部１０は、コンピュータで構成され、領域選出部１１、モデル設定部１２、および、抽出部メイン１３は、コンピュータ上で動作するプログラムで構成されている。データ入力部１４、操作部１５、および表示部１６は、通常のコンピュータに備えられる、ＵＳＢポート、ＤＶＤプレーヤ、ハードディスク、キーボード、マウス、ポインタ、フラットパネルディスプレイ等の機器からなる。
【００６９】
領域選出部１１は、複合材料の３次元画像中にフィラーを表す領域であると推定される領域を含む候補領域を選出して設定する。モデル設定部１２は、領域選出部１１によって設定された候補領域に対し、形状と配置をそれぞれ変化させるパラメータを有してフィラーの形状を再現する形状モデルを設定する。抽出メイン部１３は、モデル設定部１２によって設定された形状モデルのパラメータをランダムに変化させることによりモンテカルロ法に基づいて候補領域に形状モデルをフィッティングさせる。抽出メイン部１３は、そのフィッティングした形状モデルの形状と配置をフィラーの配向情報として抽出する。制御部１０は、３次元画像中のフィラーの抽出が完了するまで、モデル設定部１２と抽出メイン部１３の動作を反復させる。領域選出部１１、モデル設定部１２、抽出メイン部１３、およびデータ入力部１４は、上述した図２における領域選出工程（＃１）、モデル設定工程（Ｓ１）、抽出メイン工程（Ｓ２）、およびデータ入力工程（＃０）の各処理を実行する。また、画像処理装置２は、図１、図２のフローチャートによって説明した処理に限らず、他のフローチャートによって説明した処理を行うことができる。
【００７０】
（実施例）
図３２、図３３は実施例を示す。図３２は、図３（ｂ）に示した複合材料Ｃについて、本画像処理方法および画像処理装置を用いて抽出した繊維状フィラーの３次元画像Ｇ１１を示す。複合材料Ｃ（供試材料）は、上述の図３（ａ）に示した、外形長さ１７．８×幅１．８３×高さ０．４ｍｍ、開口部長さ１６．５ｍｍ、底面部厚み０．１２ｍｍ、側面厚み０．１６３ｍｍの各寸法を有する箱形状の薄肉射出成形品１の中央部を切り出したものである。複合材料Ｃは、平均繊維長８８μｍのガラス繊維をフィラーとして添加した超耐熱性のＩ型と呼ばれる液晶ポリマーであり、成形は最大型締力１９６ｋＮのプリプラ式射出成形機を用いて行った。
【００７１】
Ｘ線ＣＴによる３次元画像は、マイクロフォーカスＸ線ＣＴを用いて撮像した成形品の断層画像のデータから再構成した。再構成した画像は、上述の図４に３次元画像Ｇ１として示されており、Ｌ字型の成形品形状と、内部のガラス繊維の存在を目視確認できる。３次元画像データの大きさは、３００×３４８×２００ボクセルであり、１ボクセルの１辺の実サイズは約３μｍである。この実施例では、形状モデルとして仮想円柱を用いて抽出を行った。ガラス繊維の直径が１０μｍと既知である。従って、仮想円柱のパラメータとして、両端点の３次元座標（３×２＝６自由度）のみを、モンテカルロ法によって変化させて最適化した。適合度（フィッティングの度合い）の評価は、図１８，図１９に示した評価点ｂ（サンプリング点）を用いる方法で行った。
【００７２】
図３３は、本画像処理方法および画像処理装置を用いて抽出したフィラーの繊維長の分布を示す実施例であり、実測したフィラーの繊維長の分布を比較例として、一緒に示す。比較例の実測値は、供試材料の樹脂部分を焼いてガラス繊維のみを取り出し、その長さを実測して得たものである。実施例と比較例の両者の繊維長分布は、良い一致を示している。また、例えば、度数の最大を与える６０〜８０μｍのビンにおいて、実施例では２２．２％であり、比較例では２７．０％であり、その差は４．８％に納まっている。また、本画像処理方法および画像処理装置を用いる実施例により、供試材料を破壊することなく、非破壊で、フィラーの配向情報を抽出することができた。ここでは、フィラーの配向情報を抽出する際に得られるフィラーの繊維長の分布のみを示したが、複合材料Ｃ（供試材料）におけるフィラーの配向や配向度の分布は、個々のフィラーの配向情報から、直接、定量的に算出することができる。
【００７３】
なお、本発明は、上記構成に限られることなく種々の変形が可能である。例えば、上述した各実施形態および変形例の構成を互いに組み合わせた構成とすることができ、また、各工程の順番を適宜入れ替えた構成とすることができる。例えば、評価点設定工程（Ｓ２１）は、領域選出工程（＃１）の後であって、モデル設定工程（Ｓ１）よりも前に、備える構成としてもよい。この場合、評価点（サンプリング点）は、例えば、ボクセルを細分化した各分割点に設定すればよく、形状モデルに依存することなく３次元画像空間に固定した評価点とすることができる。また、評価値減算工程（Ｓ２３）、ペナルティ工程（Ｓ２４）、および、重み付け工程（Ｓ２５）は、評価点を用いない場合、すなわち、ボクセルのデータ値を評価値とする場合にも適用することができる。
【００７４】
また、３次元画像の分解能が高く、ボクセルのサイズが形状モデル（フィラー）のサイズよりも十分小さい場合には、評価点（サンプリング点）を設定する必要はなく、逆に、ボクセルを統合する工程を備える構成とすることができる。例えば、４つのボクセルをそれらのデータ値の平均値を新たなデータ値とする１つの大きなボクセルとすればよい。３次元画像のデータは、Ｘ線ＣＴによるデータに限らず、任意の３次元画像取得手段によるデータを用いることができる。例えば、ＭＲＩ画像データを用いることができ、また、透明樹脂の場合には光学的ＣＴ画像データを用いることができる。また、それらのデータ形式は、ボクセルに限らず、任意形状の立体画素を用いることができる。
【符号の説明】
【００７５】
１１    領域選出部（領域選出手段）
１２    モデル設定部（モデル設定手段）
１３    抽出メイン部（抽出手段）
２    ３次元画像処理装置
ｂ    評価点（サンプリング点）
Ｂｘ    ボクセル
Ｃ    複合材料
Ｃ１，・・，Ｃｎ，Ｃｉ    制御点
Ｆ    フィラー
Ｇ１〜Ｇ１１    ３次元画像
Ｍ    仮想円柱（形状モデル）
ＭＰ    仮想パイプ（形状モデル）
Ｓｐ    スプライン曲線
【図１】
0005844921.tif 000002
【図２】
0005844921.tif 000003
【図３】
0005844921.tif 000004
【図４】
0005844921.tif 000005
【図５】
0005844921.tif 000006
【図６】
0005844921.tif 000007
【図７】
0005844921.tif 000008
【図８】
0005844921.tif 000009
【図９】
0005844921.tif 000010
【図１０】
0005844921.tif 000011
【図１１】
0005844921.tif 000012
【図１２】
0005844921.tif 000013
【図１３】
0005844921.tif 000014
【図１４】
0005844921.tif 000015
【図１５】
0005844921.tif 000016
【図１６】
0005844921.tif 000017
【図１７】
0005844921.tif 000018
【図１８】
0005844921.tif 000019
【図１９】
0005844921.tif 000020
【図２０】
0005844921.tif 000021
【図２３】
0005844921.tif 000022
【図２４】
0005844921.tif 000023
【図２５】
0005844921.tif 000024
【図２８】
0005844921.tif 000025
【図２９】
0005844921.tif 000026
【図３０】
0005844921.tif 000027
【図３１】
0005844921.tif 000028
【図３３】
0005844921.tif 000029
【図２１】
0005844921.tif 000030
【図２２】
0005844921.tif 000031
【図２６】
0005844921.tif 000032
【図２７】
0005844921.tif 000033
【図３２】
0005844921.tif 000034