(57)【要約】【課題】対象物を構成する部品のセグメンテーションを効率よく行うことができる技術を提供する。【解決手段】ＣＴ画像分割部１は、ＣＴ値を用いてＣＴ画像を複数のセルに分割する。階層構造生成部２は、複数のセルを融合して新しいセルを生成し、融合段階に応じた階層構造を生成する。階層選択受付部３は、ユーザからの階層の選択を受け付ける。セグメンテーション受付部４は、選択された階層におけるいずれかのセルに対するセグメンテーションの指令を受け付ける。【選択図】図１
2018040790.tif 000002
【特許請求の範囲】
【請求項１】
画像構成単位に割り当てられたＣＴ値を有するＣＴ画像におけるセグメンテーションを行うための部品形状抽出システムであって、
ＣＴ画像分割部と、階層構造生成部と、階層選択受付部と、セグメンテーション受付部とを備えており、
前記ＣＴ画像分割部は、前記ＣＴ値を用いて前記ＣＴ画像を複数のセルに分割する構成となっており、
階層構造生成部は、前記複数のセルを融合して新しいセルを生成し、融合段階に応じた階層構造を生成する構成となっており、
前記階層選択受付部は、ユーザからの階層の選択を受け付ける構成となっており、
前記セグメンテーション受付部は、選択された階層におけるいずれかの前記セルに対するセグメンテーションの指令を受け付ける構成となっている
部品形状抽出システム。
【請求項２】
前記ＣＴ画像はＣＴボリュームであり、前記画像構成単位はボクセルである
請求項１に記載の部品形状抽出システム。
【請求項３】
前記ＣＴ画像分割部は、前記ＣＴ値の勾配のノルムをモース関数値とするモース複体の生成により、前記ＣＴ画像の分割を行う構成となっている
請求項１又は２に記載の部品形状抽出システム。
【請求項４】
前記モース関数値として、近傍の画像構成単位におけるモース関数値の偏差値が用いられている
請求項３に記載の部品形状抽出システム。
【請求項５】
前記モース関数値は、前記セグメンテーションにおける境界の方向に沿って平滑化されている
請求項３又は４に記載の部品形状抽出システム。
【請求項６】
画像構成単位に割り当てられたＣＴ値を有するＣＴ画像におけるセグメンテーションを行うための部品形状抽出方法であって、
前記ＣＴ値を用いて前記ＣＴ画像を複数のセルに分割するステップと、
前記複数のセルを融合して新しいセルを生成し、融合段階に応じた階層構造を生成するステップと、
ユーザからの階層の選択を受け付けるステップと、
選択された前記階層におけるいずれかの前記セルに対するセグメンテーションの指令を受け付けるステップと
を備える部品形状抽出方法。
【請求項７】
請求項６記載の部品形状抽出方法をコンピュータにより実行するためのコンピュータプログラム。
【請求項８】
画像構成単位に割り当てられたＣＴ値を有するＣＴ画像におけるセグメンテーションを行うためのユーザインタフェース装置であって、
階層選択受付部と、セグメンテーション受付部とを備えており、
前記ＣＴ画像は、複数のセルに分割されており、かつ、前記ＣＴ画像は、前記複数のセルの融合段階に応じた階層構造を備えており、
前記階層選択受付部は、ユーザからの階層の選択を受け付ける構成となっており、
前記セグメンテーション受付部は、選択された階層におけるいずれかの前記セルに対するセグメンテーションの指令を受け付ける構成となっている
ユーザインタフェース装置。
【請求項９】
さらに表示部を備えており、
前記表示部は、前記セグメンテーションの結果として得られた部品ごとに識別された画像を表示する構成となっている
請求項８に記載のユーザインタフェース装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、実物体についての放射線投影像に基づいて、実物体を構成する部品の形状を抽出する技術に関するものである。
【背景技術】
【０００２】
Ｘ線ＣＴ（Computed Tomography）は、被写体（対象物）に多方向からＸ線を照射し、得られたＸ線強度分布（Ｘ線投影像）に基づいて、被写体内部のＣＴ画像を再構成することができる技術である。ここで、この明細書において投影像とは、特に言及しない限り、いわゆる透過像を含む意味である。このようなＸ線ＣＴは、従来から、医療用途において広く用いられているが、近年では、産業用にも用いられるようになった。Ｘ線ＣＴを用いると、外部からは見えない形状（例えば被写体内部の空孔の形状）を特定することができるので、例えば、複雑な形状を持つ製品の寸法検査やリバースエンジニアリングにおいて有用である（下記特許文献１参照）。
【０００３】
ところで、Ｘ線ＣＴにおけるＣＴ画像の再構成においては、例えば３次元ＣＴを前提にすると、空間を埋めるボクセルにＣＴ値を割り当てることで、再構成を行っている。２次元ＣＴの場合は、ボクセルの代わりにピクセルにＣＴ値が割り当てられる。
【０００４】
ここで、被写体として、複数部品の組み合わせにより構成されるアセンブリ品（例えば自動車用エンジン）を用い、アセンブリ品のＣＴ画像を用いて、部品ごとの領域に分割すること（いわゆるセグメンテーション）を考える。セグメンテーションを自動的に行う手法として、ＣＴ値の等値面を用いる方法が考えられる。しかしながら、この方法では、適切にセグメンテーションを行うための閾値の設定が難しく、特に複雑な形状の場合は、非常に困難となる。また、ＣＴ画像には、様々な要因によるアーチファクトが存在しており、この場合にも、自動的なセグメンテーションは困難となる。
【０００５】
そこで、従来は、ボクセル単位でのセグメンテーション指定をオペレータが行う手法が主流となっている。しかしながら、この作業は非常に煩雑であり、長期間を要することが通常である。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】国際公開ＷＯ２０１３／００５４５５公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明は、前記の状況に鑑みてなされたものである。本発明は、対象物を構成する部品のセグメンテーションを効率よく行うことができる技術を提供することを目的としている。
【課題を解決するための手段】
【０００８】
前記した課題を解決する手段は、以下の項目のように記載できる。
【０００９】
（項目１）
画像構成単位に割り当てられたＣＴ値を有するＣＴ画像におけるセグメンテーションを行うための部品形状抽出システムであって、
ＣＴ画像分割部と、階層構造生成部と、階層選択受付部と、セグメンテーション受付部とを備えており、
前記ＣＴ画像分割部は、前記ＣＴ値を用いて前記ＣＴ画像を複数のセルに分割する構成となっており、
階層構造生成部は、前記複数のセルを融合して新しいセルを生成し、融合段階に応じた階層構造を生成する構成となっており、
前記階層選択受付部は、ユーザからの階層の選択を受け付ける構成となっており、
前記セグメンテーション受付部は、選択された階層におけるいずれかの前記セルに対するセグメンテーションの指令を受け付ける構成となっている
部品形状抽出システム。
【００１０】
（項目２）
前記ＣＴ画像はＣＴボリュームであり、前記画像構成単位はボクセルである
項目１に記載の部品形状抽出システム。
【００１１】
（項目３）
前記ＣＴ画像分割部は、前記ＣＴ値の勾配のノルムをモース関数値とするモース複体の生成により、前記ＣＴ画像の分割を行う構成となっている
項目１又は２に記載の部品形状抽出システム。
【００１２】
（項目４）
前記モース関数値として、近傍の画像構成単位におけるモース関数値の偏差値が用いられている
項目３に記載の部品形状抽出システム。
【００１３】
（項目５）
前記モース関数値は、前記セグメンテーションにおける境界の方向に沿って平滑化されている
項目３又は４に記載の部品形状抽出システム。
【００１４】
（項目６）
画像構成単位に割り当てられたＣＴ値を有するＣＴ画像におけるセグメンテーションを行うための部品形状抽出方法であって、
前記ＣＴ値を用いて前記ＣＴ画像を複数のセルに分割するステップと、
前記複数のセルを融合して新しいセルを生成し、融合段階に応じた階層構造を生成するステップと、
ユーザからの階層の選択を受け付けるステップと、
選択された前記階層におけるいずれかの前記セルに対するセグメンテーションの指令を受け付けるステップと
を備える部品形状抽出方法。
【００１５】
（項目７）
項目６記載の部品形状抽出方法をコンピュータにより実行するためのコンピュータプログラム。
【００１６】
（項目８）
画像構成単位に割り当てられたＣＴ値を有するＣＴ画像におけるセグメンテーションを行うためのユーザインタフェース装置であって、
階層選択受付部と、セグメンテーション受付部とを備えており、
前記ＣＴ画像は、複数のセルに分割されており、かつ、前記ＣＴ画像は、前記複数のセルの融合段階に応じた階層構造を備えており、
前記階層選択受付部は、ユーザからの階層の選択を受け付ける構成となっており、
前記セグメンテーション受付部は、選択された階層におけるいずれかの前記セルに対するセグメンテーションの指令を受け付ける構成となっている
ユーザインタフェース装置。
【００１７】
（項目９）
さらに表示部を備えており、
前記表示部は、前記セグメンテーションの結果として得られた部品ごとに識別された画像を表示する構成となっている
項目８に記載のユーザインタフェース装置。
【発明の効果】
【００１８】
本発明によれば、対象物を構成する部品のセグメンテーションを効率よく行うことが可能になる。
【図面の簡単な説明】
【００１９】
【図１】本発明の一実施形態における部品形状抽出システムの概略的な構成を示すブロック図である。
【図２】本発明の一実施形態における部品形状抽出方法の概略的な手順を示すフローチャートである。
【図３】本発明の一実施形態におけるモース複体の概念を説明するための説明図であって、グレイスケールの濃淡で示される曲面は、ｘｙ平面上に定義されたモース関数の値をｚ方向にプロットした結果を表す。この図において、四角が極小点、丸が極大点、×印が鞍点を表す。
【図４】ＣＴボリュームの断面に対して、ＣＴ値の勾配のノルムをモース関数としてセル分割した結果（左上）とその拡大図（左下）、モース関数値を曲面表示した上にセル境界を重ねて表示したもの（右）である。
【図５】パーシステンスを用いたセルの統合の様子を説明するための説明図である。
【図６】パーシステンスを用いたセルの統合の様子を説明するための説明図である。
【図７】セル分割の例を模式的に説明するための説明図である。
【図８】対応するセルを統合する過程およびそのエンコードを示す説明図である。
【図９】セル数の指定による簡略化レベルの変更手法を説明するための説明図である。
【図１０】セル分割の詳細度変更により部品を決定する例を説明するための説明図である。
【図１１】実ＣＴボリュームの断面（上）とその拡大（下）におけるセル分割簡略化の様子を示す説明図である。
【図１２】セグメンテーション用GUIにおける画面表示を説明するための説明図である。
【図１３】セグメンテーションを操作する手順を説明するための説明図である。
【図１４】スライダを用いてセル分割の詳細度を調整する場合の画面表示例を示す説明図である。
【図１５】解析的手法によるＣＴ画像の断面（左側）と、逐次近似法によるＣＴ画像の断面（右側）を示す説明図である。
【図１６】ＣＴアーチファクトにより適切な詳細レベルが存在しないＣＴボリュームの例を示す説明図であって、左側の図は、ＣＴボリュームの一断面、中央の図は、左側の断面図におけるＣＴ値の勾配ノルムのグレイスケール表示、右側の図は、計測物のＣＴボリュームのボリュームレンダリングを示す。
【図１７】コスト付きパーシステンスによるセルの統合を説明するための説明図である。左側の図は、異素材（空気とアルミ）のセルｉ，ｊの統合を示し、右側の図は、同素材（アルミとアルミ）のセルｉ，ｋの統合を示す。
【図１８】左側の図は、ＣＴボリュームの一断面を示し、左から２番目の図は、そのＣＴ値のノルムのグレイスケール画像を示し、右側の二つの図は、その拡大図を示す。
【図１９】モース関数の正規化の様子を説明するための説明図であって、左側の図は、正規化前のモース関数値のグレイスケール画像を示し、右側の図は、正規化後のモース関数値のグレイスケール画像を示す。なお、これらの図においては、値の変化が見やすいよう、コントラストが個別に調整されている。
【図２０】モース関数の平滑化の効果を示す説明図であって、左側の図は、正規化したモース関数値のグレイスケール表示を示し、右側の図は、正規化および平滑化したモース関数値のグレイスケール表示を示す。
【発明を実施するための形態】
【００２０】
以下、添付図面を参照しながら、本発明の実施形態に係る部品形状抽出システム（以下単に「システム」と称することがある）について説明する。本実施形態のシステムは、画像構成単位としてのボクセルに割り当てられたＣＴ値を有するＣＴ画像におけるセグメンテーションを行うためのものである。
【００２１】
（本実施形態の構成）
本実施形態のシステム（図１参照）は、ＣＴ画像分割部１と、階層構造生成部２と、階層選択受付部３と、セグメンテーション受付部４とを基本的な構成要素として備えている。さらに、このシステムは、ＣＴ画像取得部５と、表面メッシュ生成部６と、記憶部７と、表示部８とを追加的に備えている。
【００２２】
ＣＴ画像取得部５は、ＣＴ画像を取得するための機能要素である。ＣＴ画像とは、二次元ＣＴ画像と三次元ＣＴ画像（いわゆるＣＴボリューム）とを含む概念であるが、この実施形態の説明においては、ＣＴボリュームを前提として説明する。したがって、本実施形態のＣＴ画像取得部５は、空間を埋めるボクセルにＣＴ値が割り当てられることで構成されるＣＴ画像（ＣＴボリューム）を生成するようになっている。このようなＣＴ画像取得部５としては、従来から使用されている各種の装置を用いることができるので、これ以上詳しい説明は省略する。
【００２３】
ＣＴ画像分割部１は、ＣＴ画像におけるＣＴ値を用いてＣＴ画像を複数のセルに分割する構成となっている。本実施形態のＣＴ画像分割部１は、ＣＴ値の勾配のノルムをモース関数値とするモース複体の生成により、ＣＴ画像の分割を行う構成となっている。
【００２４】
階層構造生成部２は、複数のセルを融合して新しいセルを生成し、融合段階に応じた階層構造を生成する構成となっている。
【００２５】
階層選択受付部３は、ユーザからの階層の選択を受け付ける構成となっている。
【００２６】
セグメンテーション受付部４は、選択された階層におけるいずれかのセルに対するセグメンテーションの指令を受け付ける構成となっている。
【００２７】
表面メッシュ生成部６は、セグメンテーションごとの表面メッシュを生成する構成となっている。
【００２８】
記憶部７は、生成された階層構造を記憶し、ユーザからの指定に対応して、特定階層のセル集合の情報を階層選択受付部３、セグメンテーション受付部４、及び表示部８のいずれか又は全てに送ることができるようになっている。
【００２９】
表示部８は、セグメンテーションの結果として得られた部品ごとに識別された画像を表示する構成となっている（後述の図１２参照）。
【００３０】
前記した各構成要素の詳細は、以下における部品形状抽出方法の説明において記載する。
【００３１】
（本実施形態における動作）
以下、図２をさらに参照して、本実施形態における部品形状抽出方法を説明する。
【００３２】
（図２のステップＳＡ−１）
まず、ＣＴ画像取得部５により、ＣＴボリュームを取得する。この取得の手法は従来と同様でよいので、詳しい説明は省略する。ＣＴボリュームにおいては、ボクセルにＣＴ値が割り当てられているので、任意断面におけるＣＴ断面像を生成して、ユーザに呈示することができる。
【００３３】
（図２のステップＳＡ−２）
次に、ＣＴ画像分割部１により、ＣＴ画像（本例ではＣＴボリューム）を分割して、複数のセル（複数セルを総称してセル集合と呼ぶことがある）を生成する。ここで、セルとは、複数ボクセルの集合である。以下、セルの生成手法の詳細を説明する。
【００３４】
（モース複体）
曲面上で定義されたスカラー関数ｆを考える。これは曲面上の各点に対し実数を与えるということである。図３に、２次元平面（ｘ、ｙ）を定義域とするスカラー関数ｆの例を示す。グレイスケールで表された面は、関数値を平面に垂直な方向ｚにプロットした結果である。関数値の１階微分が全て０になる点を臨界点といい、臨界点は、極大点、極小点、及び鞍点のいずれかである。図３では、四角が極小点、丸が極大点、×印が鞍点を表す。
【００３５】
臨界点における関数ｆのHesse行列の行列式が０のとき、その臨界点を退化した臨界点と呼び、それ以外の臨界点を非退化な臨界点と呼ぶ。
【００３６】
関数ｆの臨界点がすべて非退化なとき、その関数をモース関数と呼ぶ（ただし後述の通り、本実施形態ではこの概念を拡張している）。
【００３７】
この関数値のプロット曲面の形状を調べる。曲面上の各点から、最急降下方向に進む経路を辿って極小点に到達するまで進む。同一の極小点に至る最急降下経路に属する曲面上の点の集合のことを、その極小点の不安定多様体という。曲面上の各極小点に対する不安定多様体は互いに重なりがなく、また不安定多様体の境界上の点を除き、曲面上のいずれの点もただ一つの不安定多様体に属す。したがって不安定多様体による曲面の分割は複体となる。図３においては、曲線が不安定多様体の境界を表しており、すなわちこれが曲面の分割となっている。
【００３８】
（モース複体によるセル分割）
本実施形態では、システムでは境界をエッジ点の集合とする。より詳しくは、本実施形態におけるエッジ点とはＣＴ値Ｉの勾配方向において勾配ベクトルのノルムが最大となる点であり、以下で定義される。
【００３９】
2018040790.tif 000003
【００４０】
したがって本システムではＣＴ値の勾配のノルム
2018040790.tif 000004
をモース関数とする。簡単のため、以下、不安定多様体をセル、不安定多様体の複体をモース複体あるいはセル分割と呼ぶことにする。セルの境界すなわち多様体の境界が部品の境界候補の役割を果たす。図４にセル分割の結果を示す。ここでは、モース関数値を曲面表示したものにセル境界を重ねて表示した。図４の左下及び右側の図における破線は、部品の境界であるセル境界を示す。
【００４１】
また本実施形態では、一般の数学的な定義を拡張して、モース関数の臨界点の退化を許す。ＣＴボリュームにおいて、退化した臨界点はＣＴ値が変化しない領域中の点に該当するが、こういった領域は同一素材からなる部品中に一般的に現れるものである。そのため、退化した臨界点からなる連結領域は、単一の臨界点と同一視し、同じセルに属するものとして扱う。
【００４２】
（図２のステップＳＡ−３）
ついで、生成されたセル集合の融合により、階層構造生成部２が階層構造を生成する。この手法を以下に詳しく説明する。
【００４３】
（セルの統合による階層構造の構築）
図４のセル分割には小さなセルが多数あり、1つの部品を構成する全セルを選択するには非常に手間がかかることが予想される。そこで同一部品を表現する微細なセルは統合し、部品抽出の際のセル選択の手数を減らすことを考える。
【００４４】
本発明者らの知見によれば、同一部品に属する隣接する２セルが共有する境界の尾根（境界におけるモース関数値）は、他の境界よりも相対的に低い。したがって、共有する境界の尾根が低い隣接する２セルを順次統合していくことで、よりセル数の少ないモース複体を構成することができ、セル数の異なる（つまり詳細度の異なる）一連のモース複体すなわち階層的モース複体を得ることができる。
【００４５】
図５及び図６は、この統合の様子を１次元の模式図で示したものである。２つの隣接するセルｃｉ，ｃｊ（ここでｉ，ｊは任意の自然数）の間の尾根の高さを表す指標として、パーシステンスという概念を導入し、パーシステンスＰを以下のように定義する。
【００４６】
2018040790.tif 000005
【００４７】
ここでＳｉ，ｊは、セルｃｉ，ｃｊ間の鞍点、ｍｉ，ｍｊは、ｃｉ，ｃｊのモース関数
2018040790.tif 000006
の極小点を表す。ある尾根のパーシステンスは、尾根を形成する複数の谷底から測った尾根の高さのうち最小のものである。
【００４８】
統合の準備として、まず全ての隣接セルの組み合わせにおいてパーシステンスを計算しておく。統合はパーシステンスの小さい順、すなわち尾根が低い順に行うこととする。部品境界は関数値プロット曲面において高い尾根で表現されていることが期待されるため、パーシステンスが小さい順に統合を行うことは、部品境界らしくない境界候補から順に消していくことを意味する。
【００４９】
セルの統合と更新を、図５の１次元の複体の例を用いて説明する。図５における縦線（破線）は、セル境界を模式的に示す。図５では、最小のパーシステンスをもつ尾根は、図中に付記したように、真ん中の尾根である。図中における左右のセルを統合するとは、この尾根を切り崩して両側のセルを統合することである。すると、左にあったセルの左の尾根と右にあったセルの右の尾根とで囲まれる大きなセルが新しく生じる（図６における「統合されたセル」を参照）。なお、図６において統合により変更されたモース関数値を二点鎖線で示している。
【００５０】
新しいセルの左の尾根と右の尾根は谷底からの高さが変わった可能性があるため、パーシステンスを再計算する必要がある。この図では、右の尾根におけるパーシステンスは変化しないが、左の尾根では更新が起こる。左の尾根は、右にあった浅い谷が非常に深い谷に取って代わった今、左の谷底から測った高さを新しいパーシステンスとすることになる。更新後、パーシステンスを再度ソートし、その中で最小のものを選んで両側のセルを統合する。この作業を、最終的にセル数が１つだけの分割が得られるまで繰り返す。
【００５１】
図７に、あるセル分割の例を示し、図８に、その統合の様子の模式図を示す。図８におけるノードは、初期分割のセルおよび統合過程におけるセルを、辺は各段階において存在するセルを表す。階層構造の最下段が初期分割（最も細かい分割）であり、上に行くにつれてセルが一つひとつ統合され消えてゆく。結果として階層構造の統合グラフは下にいくほど広がる木構造となる。統合の過程は、統合の順位と統合先のセルのIDの組でコード化して保存・再現可能であり、木構造の下に表で示してある。事前にセル分割の階層構造を計算しコード化しておくことで、次の段階でGUIを用いてセグメンテーションを行う際に、リアルタイムでインタラクティブなセグメンテーションが可能になる。セグメンテーションにより得られた階層構造の情報は、例えば記憶部７に格納しておくことができる。
【００５２】
（図２のステップＳＡ−４）
ついで、ユーザは、階層選択受付部３を介して、階層構造を選択する。この階層構造の選択手法の詳細を、以下に説明する。
【００５３】
ユーザは、例えばセルの数を用いて詳細度を指定することができる。図９は、図８におけるセル分割において、事前に計算しておいたコード（例えば記憶部７に格納されていたコード）を用いてセル数を変更するときの様子である。説明のため、現在の詳細度のレベルを表す一点鎖線をグラフに重ねて描いてある。この一点鎖線と交差する辺が、そのとき分割に残っているセルを表す。この図の例では、現在の詳細度はセル数４である。コードを用いて、初期状態からセル数４のセル分割を得る手法を説明する。図９下部のコードの「順番」を見ると、統合順序１のセルはセル０であることがわかる。したがってセル０が最初に統合され消える。コードの「統合先」のセル０に対応する箇所を見ると、統合先はセル１であることがわかる。すなわちセル０はセル１に統合され、初期状態においてセル０に属するボクセルは統合後のセル分割ではセル１に属するボクセルとして扱われる。次は統合順序２のセルであるセル５が統合先であるセル４に統合される。次いで統合順序３のセル１（セル０が統合済みであることに注意）が統合先のセル２に統合される。ところで、初期状態（最もセル数が多い分割）のセル数（７）―統合回数（３）＝残りセル数（４）が成り立つ。３回セルを統合した時点で、残りセル数（４）＝指定セル数（４）となったので統合がストップする。
【００５４】
上記の統合例を再度図９を用いて模式的に説明する。一点鎖線と交差した辺を見ると、セル2,3,4,6が生き残っていることがわかる。このとき、全セル数（7）−指定セル数（4）＝3以下の統合順位をもつセル0，1，5が統合され簡略化されている。セル0，1，5の統合先は、コードを用いて各セルの統合先を辿って行きついた統合されないセルとなる。具体的にはセル０、１、５は、セル０→セル１→セル２、セル１→セル２、セル５→セル４のように統合操作を受ける。初期状態におけるセル０、１、５に属するボクセルは、統合終了後はセル２または４に属することになる。図９の模式図によれば、セル２、４に対応する辺はいずれも一点鎖線と交差しており、統合終了後にも生き残ったセルであることがわかる。
【００５５】
図１０は、セル数を変更して、アルミでできた連結した一部品または一部品の一部を１つのセル（ラベル２）で表している状態を示している。なお、本実施形態において、部品とは、同じ材質で構成された一体の部品の一部である場合を含む。図１０における水平方向の一点鎖線は、アルミの部品に単一のラベル２を与えることができる詳細レベルを示している。分割レベルをより詳細にする（つまり水平方向の一点鎖線の位置を下げて、より多くのセルを表示する）と、アルミはラベル２のセルだけでなく、他のラベル（例えばラベル０やラベル１）を持つ、より多くのセルで表現されることになる。このとき、同一部品の抽出にはこれらすべてを選択する必要があり、作業量が増えてしまうことがわかる。
【００５６】
図１１には実際のＣＴ画像の断面におけるセル分割の統合の様子を示す。また、図１２の左側に示すように、選択された詳細度のセルは、そのセルを含むＣＴ断面像に重畳されて、表示部８によりユーザに呈示される（後述）。図１３のステップ１は、図２のステップＳＡ−４に対応する。
【００５７】
図２のステップＳＡ−４では、所望の境界に線が現れかつなるべくセル数が少なくなるように詳細度を調節する。本実施形態では、詳細度を示すスライダを表示部８において表示することにより、ユーザによる詳細度の指定操作を容易としている（図１４の右側の図参照）。あるいは、ショートカットキーにより詳細度指定を可能にすることもできる。なお、本実施形態においては、一つの融合ステップごとに一つの詳細度レベルが割り当てられるので、詳細度レベルの数は、セルの分割数にほぼ等しい。しかし、分割数が多い場合における詳細度レベルの相違は、セグメンテーションにあまり寄与しないことが多いので、詳細度（分割数）を対数表示することにより、セグメンテーション操作の効率化を図ることができる。
【００５８】
（図２のステップＳＡ−５）
階層を選択した後、ユーザは、セグメンテーション受付部４及び表示部８を用いて、セルごとのセグメンテーションを行う。
【００５９】
このセグメンテーションの手順を以下に詳しく説明する。図１２は、表示部８によってユーザに呈示される基本画面表示である。この表示は、ＣＴボリュームの一断面とそのセル分割を表示するスライスビュー画面（図中左）と、３次元表示画面（図中右）の２画面を有する。３次元表示画面ではＣＴボリュームのボリュームレンダリングまたは境界面（セグメンテーション指定の結果としての境界面）のいずれかを表示することができるようになっている。
【００６０】
このステップＳＡ−５では、未選択のセルの位置を３次元表示画面で確認しながら、そのセルを含む断面をスライスビュー画面に表示し、目的のセルをクリックして選択する。選択が終了したら、選択したセルに、部品に応じたラベルをユーザ操作（例えば、プルダウンメニューによるラベル番号の選択）により付ける。これらの指定は、セグメンテーション受付部４により、システム側に受け付けられる。
【００６１】
また、本実施形態では、作業効率が向上するように、ラベルごとにユーザが色を指定できるようになっている。この指定も、セグメンテーション受付部４により、システム側に受け付けられる。図１３では、色の違いをグレイスケールにおける濃度の違いにより示しており、図中右下と右上はそれぞれ、色分けによりラベル付けされた２部品が表示されたスライスビュー画面と３次元表示画面の様子である。なお、本実施形態では、選択されたセルは、図１２右側の３次元画面表示において、対応する色（好ましくは画面左側に表示されたスライスビュー画面と同じ色）が表示されるようになっている。なお、図１３のステップ２が図２のステップＳＡ−５に対応する。
【００６２】
（効率化のための工夫）
ユーザの作業量と部品抽出にかかる時間を減らすため、このステップＳＡ−５では、指定したセルのＣＴ値から所定範囲内のＣＴ値をもつセルを一括して選択できる（つまりラベル付けできる）機能を実装することが好ましい。このようにすると、例えば、樹脂と金属からなるアセンブリ品のＣＴボリュームに対して、本機能を用いて、金属部分だけ一括して抽出し、ラベル付けを行うことができる。
【００６３】
また３次元表示画面に境界面を表示し、その一点をクリックすると、クリックした点を含むセルをスライスビュー画面に表示する機能を実装することも好ましい。このようにすると、境界面を参照しつつセル指定を行うことが可能になり、ユーザの作業効率を向上させることができる。
【００６４】
（図２のステップＳＡ−６）
セグメンテーション指定されるべきセルが残っている場合は、ユーザは、ステップＳＡ−４に戻って、再び階層選択を行い、以降の動作を繰り返す。もちろん、同じ階層におけるセグメンテーション指定を継続したい場合は、同じ階層におけるセグメンテーション指定を繰り返すことができる。セグメンテーション指定が完了した場合は、ステップＳＡ−７に進む。
【００６５】
（図２のステップＳＡ−７及びＳＡ−８）
つぎに、表面メッシュ生成部６により、セグメンテーションごとに表面メッシュを生成する。この手順を本明細書では「ポリゴン化」と称することがある。以下において詳しく説明する。
【００６６】
前記したセグメンテーションの手順（ステップＳＡ−５）では、目的とする部品に該当するボクセルを抽出した。ステップＳＡ−７では、抽出されたセグメンテーション（部品に対応）の表面をポリゴン化する。ここで、ＣＴボリュームにおける表面メッシュ生成のために従来からよく用いられるのは、ＣＴ値の等値面であるが、これはＣＴ値のゆらぎの影響を受けやすく、精度が低いという問題がある。本実施形態では、ステップＳＡ−２において部品の境界を求めたのと同様に、「ＣＴ値Ｉの勾配ベクトルｇのノルムの勾配方向の極大点集合」を部品表面と定義する。極大点集合は以下のように書き表せる。
【００６７】
2018040790.tif 000007
2018040790.tif 000008
【００６８】
上記の極大点集合をボクセル上で線形補間により抽出し、つないで表面メッシュを生成することができる。得られた表面メッシュはボクセルの粗さよりも高い精度を有しており、サブボクセル精度のメッシュと言える。したがって、本実施形態では、高精度の表面形状を生成することができ、例えば部品の寸法測定を高精度で行うことができる。
【００６９】
得られた表面メッシュは表示部８を介してユーザに呈示される（図１２右側参照）。
【００７０】
（変形例１）
以下、前記した実施形態のオプションとしての変形例１を説明する。
【００７１】
工業製品のＣＴボリュームでは、アーチファクトがひどいために、適切な境界をもったセル分割が得られない場合がある。そのようなとき、サイノグラムが得られるのであれば、逐次近似法によるＣＴ再構成を行い、得られたＣＴボリュームを入力として上述のセル分割手法（ステップＳＡ−２）を適用することが好ましい。
【００７２】
ＣＴ画像を積分変換により投影像に変換することをラドン変換と呼ぶ。ＣＴ再構成法は、その計算方法により解析的手法と逐次近似法の二つに分けられる。解析的手法はラドン変換の逆変換を数式によりモデル化してＣＴ画像を計算する手法、逐次近似法はX線投影をモデル化して得た線形システムを最適化により解く手法である。
【００７３】
解析的手法は、計算時間が短いため多くの場面で使われているが、ＣＴアーチファクト等に弱いという欠点がある。それに対し逐次近似法は、アーチファクトの原因となる物理現象を最適化条件に含めることができるため、散乱線等の影響の少ないＣＴボリュームを得ることができる。ただし計算時間が解析的手法の数倍から10倍程度かかるため、必要に応じてこの２者を使い分ける必要がある。
【００７４】
解析的手法のなかでもFBP法と呼ばれる手法は、良好なサイノグラムが入力として与えられた場合比較的アーチファクトの少ないＣＴ画像を生成することができる。計測物が重金属などを含んでおりメタルアーチファクトが生じる場合は、逐次近似法を使うことが好ましい。
【００７５】
図１５は、ラジコンカーのＣＴボリュームを解析的手法と逐次近似法でそれぞれ計算し一断面を表示したものである。解析的手法による結果（図１５左側）では、金属であるモーターやシャフトの周辺にアーチファクトが生じ色むらが起こっているが、逐次近似法による結果（図１５右側）では各部品が均一な値（色）で表現されていることが観察できる。
【００７６】
（変形例２：セルの階層構造の改善）
（上記実施形態の課題と改善点）
上記のアルゴリズムは、ＣＴボリュームを用いてアセンブリの部品ごとにセグメンテーションを行うことを可能にするが、アーチファクトを起こしているＣＴボリュームに対しては適切なセル分割を行うことは難しい。このような場合、オペレータは多数のセル選択を強いられる恐れがある。
【００７７】
図１６は、鉄（中央の円筒部分）とアルミニウム（鉄の周辺部分）からなる製品のＣＴボリュームを示す。この図では、重金属である鉄によりＣＴアーチファクトを起こしている。ＣＴ値の勾配ノルムのグレイスケール表示（中央）を見ると、境界面の一部が消えかけていることが分かる。これはＣＴアーチファクトの影響により物質間の境界面においてＣＴ値が連続的に変化してしまい、勾配ノルムが十分に大きくならないためである。
【００７８】
このような場合は、セルの階層構造を生成しても、セグメンテーションに最適な階層（分割数）が存在しないことがある。詳細度を上げて必要な境界面を表示すればセル数は必要以上に増え、詳細度を下げればセル数の数は減る代わりに、必要な境界面が消えてしまう。
【００７９】
これは、モース関数にＣＴ値の微分値を用いているためである。ＣＴ値の微分値はＣＴボリュームに含まれるノイズに強く影響されるため、ＣＴアーチファクトを含むＣＴボリュームに対して上記のような問題が起こるのである。
【００８０】
これらの問題を解決するために、この変形例２では、各部品間の境界が互いにほぼ等しい強さで出現するようなセル分割を得ることを目指し、パーシステンスの定義の改良を行う。なお、後述の変形例３では、モース関数の修正により前記の問題に対応する。
【００８１】
（パーシステンスの定義の改良）
（ＣＴアーチファクトがセル境界に及ぼす影響）
例えば、メタルアーチファクトが生じると、単一素材内であってもＣＴ値がゆらぎ、部品境界面においてＣＴ値が連続的に変化する。結果として、同一部品内に複数のセルが存在する一方で、部品間の境界面が一部消滅するという現象を生じる。これは、メタルアーチファクトの影響により、前記した実施形態で定義したモース関数が正しい境界において大きなピークを持たないためである。その結果として、適切な統合順序をもつ階層的セル分割の構築に失敗することがある。
【００８２】
（コストつきパーシステンス）
すでに説明したパーシステンスの定義は位相特性を考慮したものであったが、ＣＴアーチファクトによりＣＴ値がゆらぐと、適切な統合順序が得られない場合がある。そこで、この変形例２では、隣接する二つのセルにおけるＣＴ値の値域によりパーシステンスの値をスケーリングし、アーチファクトの影響を軽減することを提案する。
【００８３】
隣接するセルｃｉ，ｃｊに対し以下のコスト付きパーシステンスＰ'を定義する。
【００８４】
2018040790.tif 000009
【００８５】
ここで
2018040790.tif 000010
はそれぞれセルｃｉ，ｃｊ中のＣＴ値の最大値と最小値を表し、
2018040790.tif 000011
は前記した実施形態で定義したパーシステンス
2018040790.tif 000012
である。
【００８６】
図１７に、コスト付パーシステンスを用いたセル統合を模式的に示す。図中において白い弧は鉄の円筒の一部、黒色のセルは空気、その他の灰色部分はアルミを表す。図１７左は異素材を表す２セル（セルｉ，セルｊ）の統合の様子を、図１７右は同素材に対応する２セル（セルｉ，セルｋ）の統合の様子を示している。
【００８７】
異素材に対応する２セルでは、
2018040790.tif 000013
（ＣＴ値がより大きい素材に対応するセル中のＣＴ値の最大値であると予想される）と
2018040790.tif 000014
（ＣＴ値がより小さい素材に対応するセル中のＣＴ値の最小値であると予想される）の差は大きく、Ｐ'の値は大きくなり、統合の順序が遅くなることが期待される。
【００８８】
それに対し、同素材に対応する２セルでは
2018040790.tif 000015
と
2018040790.tif 000016
は近い値となるためＰ'の値は小さくなり、早い段階で統合されると期待される。
【００８９】
コスト付きパーシステンスＰ'によるセル分割を行った場合には、局所的なＣＴ値の値域によるスケーリングを行うために、より少ないセル数でも、同程度の正確さで部品の形状を表すことができると考えられる。
【００９０】
（変形例３：モース関数の改良）
（モース関数値のゆらぎによる問題）
前記した実施形態において、セルの統合を進めて分割の詳細度を下げると、ある種の（例えば鉄の）部品の境界面だけが残り、他の種の（例えばアルミニウムと空気とを分ける）境界面が全て消えてしまうことが生じうる。このような、部品間の境界のいわば「強さ」の偏りがある状態において、すべての部品の境界を表せる詳細度を選択すると、セルの数が過剰になってしまうという問題を生じる。
【００９１】
こういった問題が起こる理由としては、モース関数として採用したＣＴ値の微分値のノルム
2018040790.tif 000017
が同一素材内で不均一であることが挙げられる。図１８に、ＣＴ値の勾配ノルムのグレイスケール画像およびその拡大図を示す。部品境界が白くはっきりしている箇所がある一方、非常に薄く殆ど視認できない程度の箇所も存在する。
【００９２】
そこでこのような部品間の境界の強さにむらが生じる問題を軽減するために、以下の２点について勾配ノルム場の改良を行う。
・モース関数の正規化。境界を形成する２材質のＣＴ値の違いに起因する境界の「強さ」の偏りを軽減しセル分割に適したモース関数を生成する。
・モース関数の平滑化。正規化された勾配ノルムから、エッジだと推測される部分を補強する。
【００９３】
通常の２次元画像における境界の抽出において勾配ノルムの正規化と平滑化が有効なことは先行研究（Pravin Bhat, C. Lawrence, Zitnick Michael Cohen, and Brian Curless. Gradientshop: A gradient-domain optimization framework for image and video filtering. ACM Transactions on Graphics, 2010.）で示されており、本変形例ではこれをＣＴ値のノルム場に適用する。
【００９４】
モース関数の正規化では、注目する勾配ノルムの周囲に対する大きさを考慮することで部品間の境界の強さの偏りをなくし、いずこにおいても同程度の強さの境界を得ることを狙う。モース関数の平滑化では、正規化された勾配ノルムの値を、境界が繋がっていると推測される方向にそって拡散する。これにより必要な境界のみ抽出することが可能になる。
【００９５】
（モース関数の正規化）
任意のボクセルｉについて、
2018040790.tif 000018
個の近傍ボクセル集合をＳｉとする。ボクセルｉにおけるモース関数値
2018040790.tif 000019
をその偏差値
2018040790.tif 000020
で置き換えることにより正規化し、値
2018040790.tif 000021
を新たなモース関数値とする。
2018040790.tif 000022
2018040790.tif 000023
2018040790.tif 000024
【００９６】
これにより強すぎる境界を弱め、逆に微かなエッジは補強することができる。なお、前記の例は、モース関数値として、近傍の画像構成単位におけるモース関数値の偏差値を用いる場合の具体例に対応する。
【００９７】
図１９は図１６に示したＣＴボリュームに対し、正規化なしのモース関数値ｆｉと正規化したモース関数値
2018040790.tif 000025
を計算した結果の比較である。正規化なしのモース関数では、鉄でできた中心円筒の境界が消えかけている一方、正規化したモース関数ではすべてのエッジが比較的くっきりした白い線で表示されている。なお、当然のことながら、モース関数値の正規化は、必須の手順ではなく、用途に応じて適宜省略することができる。
【００９８】
ところで図１９を観察すると、正規化モース関数
2018040790.tif 000026
では、正規化なしモース関数ｆｉの結果にはなかったところにエッジが生じており、また勾配ノルムが平均化されたために全体のコントラストが低くグレーがかった画像となっている。これらもまた適切なモース複体を計算する障害となることが予想されるため、正規化したモース関数
2018040790.tif 000027
をさらに平滑化することで解決をはかる。
【００９９】
（モース関数の平滑化）
モース関数の正規化のおかげで境界の強さを均一化することには成功した一方で、本来無い場所に境界が生じたり、全体的にコントラストが低い画像となってしまうという問題が生じた。正規化したモース関数の値を境界方向のみに沿って平滑化することで、このような問題を解決する。
【０１００】
現在注目しているボクセルに対し、境界上のボクセルの値を大きなコストで足し込み、境界面上にないボクセルは小さなコストで足し込むことにする。注目するボクセルをｉとし、その正規化モース関数値を
2018040790.tif 000028
とする。このボクセルの最近傍の6つのボクセルの集合をＡｉとし、注目ボクセルｉからＡｉ中のいずれかのボクセルｊへ向かう単位ベクトルをｅｉｊとする。また、
2018040790.tif 000029
はボクセルｉにおけるＣＴ値の勾配ベクトルを表す。ボクセルｉの平滑化したモース関数値
2018040790.tif 000030
を以下で与える。
【０１０１】
2018040790.tif 000031
【０１０２】
重みＷｉｊはボクセルｉに足し込む際のボクセルｊのコストであり、
2018040790.tif 000032
に対するｅｉｊの正射影の逆数を表す。
【０１０３】
ボクセルｊが境界面付近にあるとき、
2018040790.tif 000033
とｅｉｊは直角に近い方向を向き、正射影は小さくなり重みＷｉｊは大きい値となる。逆にボクセルｊが境界面から離れているとき、
2018040790.tif 000034
とｅｉｊは同じ方向を向き、正射影は大きくなるためＷｉｊの値は小さくなる。なお、前記の例は、モース関数値を、セグメンテーションにおける境界の方向に沿って平滑化する場合の具体例に対応する。
【０１０４】
図２０は上記の平滑化の結果を示している。正規化したモース関数では不要な境界が生じ全体的にコントラストが低くなっていたが、平滑化後は真の境界は太く補強され、それ以外の領域は黒色が濃くなりその面積も広がっている。空気に該当する箇所に生じた不要なエッジも補強されているように見えるが、これらは線分状の形状をしており、環状となって領域の境界をなしてはいない。そのため、本実施形態の手法におけるセルの計算法においてはセル境界として機能しない。したがってこの正規化・平滑化により、部品間の境界を適切に表すセル分割が得られると期待される。
【０１０５】
前記した実施形態における階層選択受付部３と、セグメンテーション受付部４と、表示部８とは、本発明のユーザインタフェース装置における具体例に対応する。
【０１０６】
なお、本発明は、前記した実施の形態に限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々の変更を加え得るものである。
【０１０７】
例えば、前記した各構成要素は、機能ブロックとして存在していればよく、独立したハードウエアとして存在しなくても良い。また、実装方法としては、ハードウエアを用いてもコンピュータソフトウエアを用いても良い。さらに、本発明における一つの機能要素が複数の機能要素の集合によって実現されても良く、本発明における複数の機能要素が一つの機能要素により実現されても良い。
【０１０８】
また、本発明を構成する各機能要素は、離散して存在しても良い。離散して存在する場合には、例えばネットワークを介して必要なデータを受け渡すことができる。各部の内部における各機能も、同様に、離散して存在することが可能である。例えば、グリッドコンピューティングやクラウドコンピューティングを用いて、本実施形態における各機能要素あるいはその一部分を実現することも可能である。
【０１０９】
また、前記した例では、画像構成単位としてボクセルを用いたが、ピクセルを用いて２次元ＣＴ画像を得る構成であってもよい。ボクセルの生成方法にも特に制約はない。画像構成単位としては、八分木構造として分割生成されたオクタントであってもよい。さらに、多面体を要素とするメッシュ構造（例えば四面体メッシュや六面体メッシュ）を画像構成単位として用いることもできる。
【符号の説明】
【０１１０】
１  画像分割部
２  階層構造生成部
３  階層選択受付部
４  セグメンテーション受付部
５  画像取得部
６  表面メッシュ生成部
７  記憶部
８  表示部
【図１】
2018040790.tif 000035
【図２】
2018040790.tif 000036
【図３】
2018040790.tif 000037
【図４】
2018040790.tif 000038
【図５】
2018040790.tif 000039
【図６】
2018040790.tif 000040
【図７】
2018040790.tif 000041
【図８】
2018040790.tif 000042
【図９】
2018040790.tif 000043
【図１０】
2018040790.tif 000044
【図１１】
2018040790.tif 000045
【図１２】
2018040790.tif 000046
【図１３】
2018040790.tif 000047
【図１４】
2018040790.tif 000048
【図１５】
2018040790.tif 000049
【図１６】
2018040790.tif 000050
【図１７】
2018040790.tif 000051
【図１８】
2018040790.tif 000052
【図１９】
2018040790.tif 000053
【図２０】
2018040790.tif 000054