JP2010014541A - System, method and program for recognizing euclid symmetry - Google Patents
System, method and program for recognizing euclid symmetry Download PDFInfo
- Publication number
- JP2010014541A JP2010014541A JP2008174641A JP2008174641A JP2010014541A JP 2010014541 A JP2010014541 A JP 2010014541A JP 2008174641 A JP2008174641 A JP 2008174641A JP 2008174641 A JP2008174641 A JP 2008174641A JP 2010014541 A JP2010014541 A JP 2010014541A
- Authority
- JP
- Japan
- Prior art keywords
- feature region
- pair
- symmetry
- region pair
- region
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Ceased
Links
Images
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Processing Or Creating Images (AREA)
Abstract
Description
本発明は、計測メッシュから対象物のユークリッド対称性を自動認識するユークリッド対称性認識システム、方法及びプログラムに関する。 The present invention relates to a Euclidean symmetry recognition system, method, and program for automatically recognizing Euclidean symmetry of an object from a measurement mesh.
近年、産業用高エネルギX線CT計測器の発達により、機械部品のような対象物品の計測メッシュが入手可能となり、製品開発の様々なプロセスで利用されるようになっている。CADシステム上で機械部品の形状を定義する際には、その機能の観点から、部分ソリッドをコピーし、これを平行移動、回転、平面反射の組み合わせから成る変換により配置し定義することが多い。そのため、計測メッシュ上からこのような対称性を認識することで、対称性を反映させたコンパクトな定義を持つCADモデルが生成でき、製品開発の様々なプロセスでより有効に利用することができる。 In recent years, with the development of industrial high-energy X-ray CT measuring instruments, measurement meshes of target articles such as machine parts have become available and are used in various processes of product development. When defining the shape of a machine part on a CAD system, from the viewpoint of its function, it is often the case that a partial solid is copied and placed and defined by a transformation composed of a combination of translation, rotation, and plane reflection. Therefore, by recognizing such symmetry from the measurement mesh, a CAD model having a compact definition reflecting the symmetry can be generated, and can be used more effectively in various processes of product development.
近年の代表的な計測メッシュの対称性認識手法として、まず、非特許文献1に記載されているボクセル化と投票に基づくPodolakらの手法が挙げられる。このPodolakらの手法は平面反射対称性の認識のみを目的としており、平行移動や回転を含む対称性は認識できない。
As a typical measurement mesh symmetry recognition method in recent years, first, the method of Podolak et al. Based on voxelization and voting described in Non-Patent
これに対してMitraらは、非特許文献2にて、各頂点における主曲率・主方向算出と投票に基づく、平行移動、回転平面反射の任意の組み合わせから成る対称性を認識する手法を提案している。しかしながらこの手法では、主曲率・主方向を頂点毎に計算しているため、計測ノイズに対してロバストに認識することができないといった問題点があり、さらにそれらが一意に定義できない平面領域を多く含むメッシュからは対称性の認識ができないといった問題点がある。
In contrast, Mitra et al. Proposed a method for recognizing symmetry consisting of an arbitrary combination of translation and rotation plane reflection based on the calculation of the principal curvature and direction at each vertex and voting in Non-Patent
これらの従来技術とは別に、Rusinkiewiczらの論文である非特許文献3に記載されているようなICP(Iterative Closest Point)法が知られている。また、溝口氏らの論文である非特許文献4に記載されているようなRegion Growing法(領域拡大法)も知られている。
本発明は、上記従来技術の問題点に鑑み、非特許文献3に記載されているICP法と非特許文献4に記載されているRegion Growing法を組み合わせて得た、平面領域を多く含む機械部品のような対象物品の計測メッシュから、測定ノイズに対してロバスト(robust)に計測メッシュのユークリッド対称性認識が行えるユークリッド対称性認識システム、方法及びプログラムを提供することを目的とする。
In view of the above-described problems of the prior art, the present invention provides a machine part including a large number of planar regions obtained by combining the ICP method described in Non-Patent
尚、本出願において計測メッシュのユークリッド対称性認識とは、平行移動、回転、平面反射の任意の組み合わせからなる変換によって互いに一致するメッシュ上の部分領域ペアとその変換パラメータとを計測メッシュから抽出することを意味する。 In the present application, the Euclidean symmetry recognition of the measurement mesh refers to extracting from the measurement mesh a pair of partial areas on the mesh that coincide with each other by a conversion composed of any combination of translation, rotation, and plane reflection, and its conversion parameters. Means that.
本発明の特徴は、計測メッシュデータと必要な各種のデータを記憶するデータ記憶部と、前記データ記憶部の計測メッシュデータに対してメッシュ各頂点の主曲率を計算するメッシュ主曲率算出部と、前記計測メッシュデータのメッシュ各頂点の主曲率を用いて、平面領域の境界近傍にある頂点集合を特徴領域として抽出する特徴領域抽出部と、前記特徴領域群から、形状が一致する特徴領域ペアをシード特徴領域ペアとして抽出するシード特徴領域ペア抽出部と、ICP法とRegion Growing法により同一の変換パラメータでマッピング可能な頂点ペアを逐次的に各シード特徴領域に追加することで特徴領域ペアをそれぞれ反復的に拡大し、互いにユークリッド対称な領域ペアとそれらの変換パラメータを抽出するユークリッド対称領域ペア抽出部と、前記ユークリッド対称な領域ペアとそれらの変換パラメータを評価し、単一平面反射、2−fold回転若しくは単一平行移動で表現できる領域ペアを見出し、当該領域ペアと共に対称性を表示する対称性表示部とを備えたユークリッド対称性認識システムである。 The features of the present invention are a data storage unit that stores measurement mesh data and various necessary data, a mesh main curvature calculation unit that calculates the main curvature of each vertex of the mesh with respect to the measurement mesh data of the data storage unit, Using the principal curvature of each vertex of the mesh of the measurement mesh data, a feature region extraction unit that extracts a set of vertices in the vicinity of the boundary of the planar region as a feature region, and a feature region pair having the same shape from the feature region group A seed feature region pair extraction unit that extracts as a seed feature region pair, and a vertex pair that can be mapped with the same conversion parameter by the ICP method and the Region Growing method are sequentially added to each seed feature region, respectively. Euclidean symmetric region pair extraction that iteratively expands and extracts Euclidean symmetric region pairs and their transformation parameters And evaluating the Euclidean symmetric region pairs and their transformation parameters, finding region pairs that can be expressed by single plane reflection, 2-fold rotation, or single translation, and displaying symmetry along with the region pairs A Euclidean symmetry recognition system including a display unit.
上記発明のユークリッド対称性認識システムでは、前記シード特徴領域ペア抽出部は、前記各特徴領域ペアの頂点数を比較する頂点数比較部と、前記頂点数の比較にて所定の範囲内にある特徴領域ペアについて、重心・主軸一致によるマッチングとその精度を評価するマッチング評価部と、ICP法によりマッチングし、その精度を評価して所定の範囲内にある特徴領域ペアについてシード特徴領域ペアとして抽出するシード特徴領域ペア決定部とを有するものとすることができる。 In the Euclidean symmetry recognition system of the above invention, the seed feature region pair extraction unit is a feature within a predetermined range in comparison with the vertex number comparison unit that compares the number of vertices of each feature region pair. For the region pair, matching by the center of gravity / principal axis matching and the matching evaluation unit that evaluates the accuracy thereof, matching by the ICP method, the accuracy is evaluated, and a feature region pair within a predetermined range is extracted as a seed feature region pair And a seed feature region pair determination unit.
本発明の別の特徴は、メッシュ各頂点の主曲率を計算するステップと、前記メッシュ各頂点の主曲率を用いて、平面領域の境界近傍にある頂点集合を特徴領域として抽出するステップと、特徴領域群に対して、形状が一致する特徴領域ペアをシード特徴領域ペアとして抽出するステップと、ICP法とRegion Growing法を組み合わせて、同一の変換パラメータでマッピング可能な頂点ペアを逐次的に各シード特徴領域に追加することで領域ペアをそれぞれ反復的に拡大し、互いにユークリッド対称な領域ペアとそれらの変換パラメータを抽出するステップと、前記ステップにて抽出した変換パラメータを評価し、単一平面反射、2−fold回転若しくは単一平行移動の対称性表現できる領域ペアを見出し、当該領域ペアと共に対称性を表示するステップとを有するユークリッド対称性認識方法である。 Another feature of the present invention includes a step of calculating a principal curvature of each vertex of the mesh, a step of extracting a vertex set in the vicinity of the boundary of the planar region as a feature region using the principal curvature of each vertex of the mesh, and a feature For each region group, a feature region pair having the same shape is extracted as a seed feature region pair, and an ICP method and a Region Growing method are combined, and vertex pairs that can be mapped with the same conversion parameters are sequentially added to each seed. Each of the pair of regions is repeatedly expanded by adding to the feature region, and a pair of regions that are Euclidean symmetric and their transformation parameters are extracted, and the transformation parameters extracted in the previous step are evaluated, and single plane reflection is performed. , Find a pair of regions that can express symmetry of 2-fold rotation or single translation, and display symmetry along with the pair of regions The Euclidean symmetry recognition method and a step.
上記発明のユークリッド対称性認識方法では、前記特徴領域群に対して、形状が一致する特徴領域ペアをシード特徴領域ペアとして抽出するステップは、前記各特徴領域ペアの頂点数を比較するステップと、前記頂点数の比較にて所定の範囲内にある特徴領域ペアについて、重心・主軸一致によるマッチングとその精度を評価するステップと、ICP法によりマッチングし、その精度を評価して所定の範囲内にある特徴領域ペアについてシード特徴領域ペアとして抽出するステップとを有するものとすることができる。 In the Euclidean symmetry recognition method of the above invention, the step of extracting the feature region pair having the same shape as the seed feature region pair for the feature region group includes comparing the number of vertices of each feature region pair; For the feature region pair within the predetermined range in the comparison of the number of vertices, matching by the center of gravity / principal axis matching and its accuracy are evaluated, matching is performed by the ICP method, and the accuracy is evaluated to be within the predetermined range. Extracting a feature region pair as a seed feature region pair.
本発明のさらに別の特徴は、メッシュ各頂点の主曲率を計算するステップと、前記メッシュ各頂点の主曲率を用いて、平面領域の境界近傍にある頂点集合を特徴領域として抽出するステップと、特徴領域群に対して、形状が一致する特徴領域ペアをシード特徴領域ペアとして抽出するステップと、ICP法とRegion Growing法を組み合わせて、同一の変換パラメータでマッピング可能な頂点ペアを逐次的に各シード特徴領域に追加することで領域ペアをそれぞれ反復的に拡大し、互いにユークリッド対称な領域ペアとそれらの変換パラメータを抽出するステップと、前記ステップにて抽出した変換パラメータを評価し、単一平面反射、2−fold回転若しくは単一平行移動の対称性表現できる領域ペアを見出し、当該領域ペアと共に対称性を表示するステップとをコンピュータに実行させるユークリッド対称性認識プログラムである。 Still another feature of the present invention is a step of calculating a principal curvature of each vertex of the mesh, a step of extracting a set of vertices in the vicinity of the boundary of the planar region as a feature region using the principal curvature of each vertex of the mesh, For each feature region group, a step of extracting a feature region pair having a matching shape as a seed feature region pair and a combination of the ICP method and the Region Growing method are used to sequentially add vertex pairs that can be mapped with the same transformation parameter. Each of the pair of regions is iteratively expanded by adding to the seed feature region, and a pair of regions that are Euclidean symmetric and their transformation parameters are extracted, and the transformation parameters extracted in the previous step are evaluated to obtain a single plane. A pair of regions that can express symmetry of reflection, 2-fold rotation, or single translation is found, and the symmetry is displayed together with the pair of regions. The Euclidean symmetry recognition program for executing a step in computers.
上記発明のユークリッド対称性認識プログラムでは、前記特徴領域群に対して、形状が一致する特徴領域ペアをシード特徴領域ペアとして抽出するステップは、前記各特徴領域ペアの頂点数を比較するステップと、前記頂点数の比較にて所定の範囲内にある特徴領域ペアについて、重心・主軸一致によるマッチングとその精度を評価するステップと、ICP法によりマッチングし、その精度を評価して所定の範囲内にある特徴領域ペアについてシード特徴領域ペアとして抽出するステップとをコンピュータに実行させるものとすることができる。 In the Euclidean symmetry recognition program of the above invention, for the feature region group, extracting a feature region pair whose shape matches as a seed feature region pair includes comparing the number of vertices of each feature region pair; For the feature region pair within the predetermined range in the comparison of the number of vertices, matching by the center of gravity / principal axis matching and its accuracy are evaluated, matching is performed by the ICP method, and the accuracy is evaluated to be within the predetermined range. The step of extracting a feature region pair as a seed feature region pair may be executed by a computer.
本発明のユークリッド対称性認識システム、方法及びプログラムによれば、対象物品の計測メッシュから平行移動、回転、平面反射の対称性を精度よく認識することができ、しかも平面領域を多く含む計測メッシュからのユークリッド対称性認識ができ、その上に、測定ノイズに対してロバストで、平行移動量に依存しない安定した認識ができ、さらに最大限大きな対称領域ペアを抽出して対称性と共に表現できる。 According to the Euclidean symmetry recognition system, method, and program of the present invention, it is possible to accurately recognize the symmetry of translation, rotation, and plane reflection from the measurement mesh of the target article, and from the measurement mesh including many plane areas. Euclidean symmetry recognition can be performed, and furthermore, stable recognition that is robust to measurement noise and does not depend on the amount of translation can be performed, and a maximum symmetric region pair can be extracted and expressed with symmetry.
以下、本発明の実施の形態を図に基づいて詳説する。 Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
先ず、非特許文献3に記載されているICP法について説明する。計測データペア〈X、Y〉のマッチングを行う手法として広く用いられている。この手法では、対応点の探索と変換パラメータの計算を反復的に行うことにより、点x'i∈X'とその対応点y c(i)∈Yとの点間距離の二乗和を最小化するよう、x iを
のようにx'iに変換するパラメータ
を求めることができる。尚、本明細書中では、ベクトルや行列について、符号に下線を施すことによりベクトルであることを示すことがある。 Can be requested. In the present specification, vectors and matrices may be indicated as vectors by underlining the codes.
領域ペアをマッチングさせるための変換パラメータを求める際に、このICP法を利用することは有効な手段ではあるが、このICP法では、マッチングさせる領域ペアが単一の平面をサンプリングした領域である場合には適切に変換パラメータを計算することが困難である。 It is an effective means to use this ICP method when obtaining conversion parameters for matching region pairs. However, in this ICP method, the region pair to be matched is a region obtained by sampling a single plane. It is difficult to calculate the conversion parameters appropriately.
これを解決したのが本発明のユークリッド対称性認識システム、方法及びプログラムである。図1には、本発明の1つの実施の形態のユークリッド対称性認識システム100を示している。本システム100はコンピュータシステムにより構成されていて、HMIとしてキーボード21、マウスその他のポインティングデバイス22、ディスプレイ23が接続されている。本システム100内にソフトウェアとして組み込まれたユークリッド対称性認識プログラムは、機能毎に区別して表示すると、外部とのデータの送受信、上記HMIとの信号の送受信を行う入出力制御部1、計測メッシュデータその他の演算処理データを記憶するデータ記憶部2、後述するStep1の処理を行うメッシュ主曲率算出部3と特徴領域抽出部4、Step2の処理を行うシード特徴領域ペア抽出部5、Step3の処理を行うICP処理部6と領域拡大処理部7、Step4のユークリッド対称性を認識するユークリッド対称性認識処理部8、そしてこれらの演算処理を総合的に制御する演算制御処理部9から構成されている。
A solution to this problem is the Euclidean symmetry recognition system, method and program of the present invention. FIG. 1 shows an Euclidean
そして、シード特徴領域ペア抽出部5は、各特徴領域ペアの頂点数を比較する頂点数比較部51と、頂点数の比較にて所定の範囲内にある特徴領域ペアについて、重心・主軸一致によるマッチングとその精度を評価するマッチング評価部52と、ICP法によりマッチングし、その精度を評価して所定の範囲内にある特徴領域ペアについてシード特徴領域ペアとして抽出するシード特徴領域ペア決定部53とを有している。
Then, the seed feature region
次に、上記の実施の形態のユークリッド対称性認識システムによるユークリッド対称性認識方法について説明する。図2に示すように、まず、メッシュ各頂点の主曲率を計算し、これを利用して平面領域の境界近傍にある頂点集合を特徴領域(feature region)として抽出する(Step1)。 Next, the Euclidean symmetry recognition method by the Euclidean symmetry recognition system of the above embodiment will be described. As shown in FIG. 2, first, the main curvature of each vertex of the mesh is calculated, and using this, a vertex set near the boundary of the planar region is extracted as a feature region (Step 1).
次に、各特徴領域ペアに対して、後述するStep2−1:頂点数比較、Step2−2:重心・主軸一致によるマッチングとその精度評価、Step2−3:ICP法によるマッチングとその精度評価、の3つの手順を実行し、その形状がおおよそ一致する特徴領域ペアをシード特徴領域ペア(pair of seed feature regions)として抽出する(Step2)。 Next, for each feature region pair, Step 2-1 described below: Comparison of the number of vertices, Step 2-2: Matching by center-of-gravity / main axis matching and its accuracy evaluation, Step 2-3: Matching by ICP method and its accuracy evaluation Three steps are executed, and feature region pairs whose shapes roughly match are extracted as a pair of seed feature regions (Step 2).
次に、ICP法と非特許文献4に記載されているRegion Growing法を組み合わせた手法を用いて、同一の変換パラメータでマッピング可能な頂点ペアを逐次的に各シード特徴領域に追加することで領域ペアをそれぞれ反復的に拡大し、互いにユークリッド対称な領域ペアとそれらの変換パラメータを抽出する(Step3)。
Next, by using a method combining the ICP method and the Region Growing method described in Non-Patent
最後に、Step3で抽出した変換パラメータを評価し、単一平面反射、2−fold回転、単一平行移動で表現できるものは、それらを利用した陽な対称性表現に変換する(Step4)。
Finally, the conversion parameters extracted in
Step1の主曲率計算に基づく特徴領域の抽出処理は、図3のフローチャート、図4の説明図を用いて説明する。まず、各メッシュ頂点viの主曲率(max curvature)を計算し(Step1−1)、求めた主曲率に基づき、シャープ頂点(sharp vertex)と非シャープ頂点(non-sharp vertex)を分類する(Step1−2)。これには、非特許文献4に記載され、また本願発明者らによる先の出願、特願2007−015473に記載された方法を用いる。
The feature region extraction processing based on the calculation of the main curvature in
すなわち、ノイズの多いメッシュモデル上のメッシュ曲率を算出するために、最初に、各頂点x iの周りの頂点群N(i)に対して(1)式で表わされる二次多項式曲面h(u,v)をフィッティングする。
ここで、N(i)は各頂点x iに対して位相的に連結されている頂点群(該頂点x iも含む)であり、(2)式を満たす隣接頂点群である。
ここで、li,avgは頂点群N(i)に含まれるエッジ集合の平均エッジ長さであり、Wは曲面フィットのサイズを特定する変数である。Wは経験的に3.0とするのが好ましい。 Here, l i, avg is the average edge length of the edge set included in the vertex group N (i), and W is a variable that specifies the size of the curved surface fit. W is empirically preferably 3.0.
このフィッティングにて、頂点x iにおける二次多項式曲面h(u,v)の主曲率ki,max,ki,minを算出する。 With this fitting, the principal curvatures k i, max and k i, min of the quadratic polynomial surface h (u, v) at the vertex x i are calculated.
続いて、頂点x iの最大曲率ki,maxが(3)式を満たすなら、その頂点はシャープ頂点であると分類する。
閾値thsharpは対象物品の形状や計測メッシュの解像度に応じて設定されるべきものであるが、10.0に設定すると計測メッシュモデルに対して良い結果が得られる。 The threshold th sharp should be set according to the shape of the target article and the resolution of the measurement mesh, but when set to 10.0, a good result can be obtained for the measurement mesh model.
図4(b)には、これを利用してシャープ頂点、非シャープ頂点(sharp vertical, non-sharp verticle)に分類した結果を示している。この図4(b)では、シャープ頂点群は濃い色にて非シャープ頂点と区別してある。 FIG. 4B shows the result of classification into sharp vertices and non-sharp verticles using this. In FIG. 4B, the sharp vertex group is distinguished from the non-sharp vertex by a dark color.
次に、図4(c)のように、シャープ頂点を要素として含まない、位相的に連結された非シャープ頂点集合を初期領域IRjとして抽出する(Step1−3)。また、抽出された初期領域IRjに対して最小二乗平面をフィットし、その平均誤差eiを計算する。この誤差eiがトレランスτplalavg以下であれば、図4(d)のようにIRjを平面領域PRkとして抽出する(Step1−4)。ここで、lavgはメッシュの平均稜線長である。通常τpla=0.5に設定するのが好ましい。 Next, as shown in FIG. 4C, a topologically connected non-sharp vertex set that does not include a sharp vertex as an element is extracted as an initial region IR j (Step 1-3). Further, the least square plane is fitted to the extracted initial region IR j and the average error e i is calculated. If this error e i is equal to or less than the tolerance τ pla avg , IR j is extracted as a plane region PR k as shown in FIG. 4D (Step 1-4). Here, l avg is the average edge length of the mesh. Usually, it is preferable to set τ pla = 0.5.
次に図4(e)のように、抽出された平面領域PRk内の、位相的に連結された領域境界頂点集合をBRlとして抽出する(Step1−5)。そして最後に、図4(f)のように、各領域境界頂点集合BRlからあるユークリッド距離Wlavg内にあり、平面領域PRk内にはない頂点集合を特徴領域
として抽出する(Step1−6)。通常W=10.0〜20.0に設定するのが好ましい。この例では、W=10.0を用いている。 (Step 1-6). Usually, it is preferable to set W = 10.0 to 20.0. In this example, W = 10.0 is used.
Step2のシード特徴領域ペアの抽出処理を、図5のフローチャートと図6の説明図を用いて説明する。この処理の目的は、ある変換のもとで形状がほぼ一致する特徴領域ペアを抽出することである。そのために、すべての特徴領域ペア〈FRi,FRj〉に対して、以下の手順でその形状類似度を評価し、シード特徴領域ペアを生成する。平面反射対称性をも認識するため、各特徴領域FRiと鏡像の関係にあるFRirefを含めたペア〈FRiref,FRj〉をも処理の対象とするが、これらは〈FRi,FRj〉と同様の手法で扱えるため、説明は特徴領域ペア〈FRi,FRj〉のみに限定する。
The seed feature region pair extraction process in
Step2−1では、頂点数比較を行う。そのために、特徴領域ペア〈FRi,FRj〉の頂点数|FRi|、|FRj|を比較し、これが、1.0-α/100<|FRi|/|FRj|<1.0+α/100を満たすならば、次のStep2−2に進み、そうでなければ処理を終了する(Step2−1A)。通常、α=10.0〜15.0に設定するが、この例ではα=15.0に設定している。 In Step 2-1, the number of vertices is compared. Therefore, the number of vertices | FR i | and | FR j | of the feature region pair <FR i , FR j > are compared, and this is 1.0-α / 100 <| FR i | / | FR j | <1.0 + α If / 100 is satisfied, the process proceeds to the next Step 2-2. Otherwise, the process is terminated (Step 2-1A). Usually, α = 10.0 to 15.0 is set. In this example, α = 15.0 is set.
Step2−2では、重心・主軸一致と平均最近点間距離を評価する。重心b iをb jに、主軸A iをA jに一致させるための変換パラメータ
を計算し、これを利用して特徴領域FRi内の各頂点v i,kを
と変換し、これにより特徴領域FRiを特徴領域FRjにマッチングする。ここで、
であり、特徴領域FRi内の頂点v i,kに対して、
のように変換する。そして、主軸A iとA jを一致させる回転行列R ijを計算し、その後、重心b jの周りに各頂点を次のように回転させることを意味している。
さらに、形状類似度を評価するため、ランダムに選択されたβ(%)の各v"i,kに対して、FRjの中からv"i,kとの点間距離が最小となる頂点v j,c(k)を見つけ、その点間距離eij,k pcaを計算する。(ここで、アンダーラインを付けていない符号v"i,kは頂点ベクトルv"i,kの頂点名である。)通常、β=5.0に設定するのが好ましい。また、c(k)はv"i,kの特徴領域FRj内にある最近点v i,c(k)のインデックスである。最後に、平均点間距離eij pcaを計算し、これがトレランスτpcalavg以下であれば、Step2−3に進み、そうでなければ処理を終了する。通常τpca=3.0に設定するのが好ましい。 Further, in order to evaluate the shape similarity, for each v " i, k of β (%) selected at random, the vertex where the distance between the points v " i, k from FR j is the smallest v Find j, c (k) and calculate the distance e ij, k pca between the points. (Here, the code v ″ i, k without underlining is the vertex name of the vertex vector v ″ i, k .) Normally, it is preferable to set β = 5.0. C (k) is the index of the nearest point v i, c (k) in the feature region FR j of v ″ i, k . Finally, the average point distance e ij pca is calculated and this is the tolerance. If it is equal to or less than τ pca avg , the process proceeds to Step 2-3, otherwise the process is terminated, and it is preferable to set τ pca = 3.0.
Step2−3では、ICP法によるマッチングと精度評価を行う。ここではICP法を用いて特徴領域FR'jをFRjにマッチングし、Step2−2の処理と同様に平均点間距離eij icpを計算する。ここで平均点間距離eij icpがトレランスτicplavg以下であれば、〈FRi,FRj〉をシード特徴領域ペア〈FRi seed,FRj seed〉として抽出する。通常τicp=1.0に設定するのが好ましい。 In Step 2-3, matching by ICP method and accuracy evaluation are performed. Here, the feature region FR ′ j is matched with FR j using the ICP method, and the average point distance e ij icp is calculated in the same manner as in Step 2-2 . Here, if the average inter-point distance e ij icp is less tolerance tau icp l avg, is extracted as <FR i, FR j> a seed feature region pair <FR i seed, FR j seed >. Usually, it is preferable to set τ icp = 1.0.
図6には抽出されたシード特徴領域ペアの例を示している。図4(f)に示す抽出された15の特徴領域群の中から、図6(a)は、平行移動と回転のペア、図6(b)は平行移動と回転と平面反射のペア、図6(c)は平行移動と回転と平面反射のペア、図6(d)は平行移動と回転のペア、図6(e)は平行移動と回転と平面反射のペア、図6(f)は平行移動と回転、さらに平行移動と回転と平面反射のペアを示している。1つの特徴領域に対して複数のペアが抽出され得るので、抽出されるペア数はさらに多数に上る。 FIG. 6 shows an example of the extracted seed feature region pair. Of the extracted 15 feature region groups shown in FIG. 4 (f), FIG. 6 (a) is a translation / rotation pair, FIG. 6 (b) is a translation / rotation / plane reflection pair, 6 (c) is a pair of translation, rotation and plane reflection, FIG. 6 (d) is a pair of translation and rotation, FIG. 6 (e) is a pair of translation, rotation and plane reflection, and FIG. A pair of translation and rotation, and further a pair of translation, rotation and plane reflection is shown. Since a plurality of pairs can be extracted for one feature region, the number of extracted pairs is further increased.
この後、Step3のICP法とRegion Growing法の組み合わせによるユークリッド対称性認識の処理に移行する。このStep3の処理を、図7のフローチャートと図8の説明図を用いて説明する。抽出したすべてのシード特徴領域ペア
を初期ユークリッド対称領域ペア、
The initial Euclidean symmetric region pair,
にまず設定し、以下の手順でそれぞれの領域SRi1 0,SRi2 0を反復的に拡大し、最終的なユークリッド対称領域ペアとその変換パラメータを抽出する。 First, each region SR i1 0 , SR i2 0 is iteratively enlarged by the following procedure, and a final Euclidean symmetric region pair and its transformation parameter are extracted.
そのために、Step3−1において、初期変換パラメータの計算と初期マッチングを行う。この処理は、Step2−2の処理と同様の処理であり、領域SRi1 0,SRi2 0に対して初期変換パラメータ
を計算し、これをもとに頂点を
と更新し、初期マッチングを行う。ここで、この処理の繰り返しを表す変数itrを、itr=0に初期化する。 And update the initial matching. Here, a variable itr representing repetition of this process is initialized to itr = 0.
続く、Step3−2では、変換パラメータの計算とマッチングの処理を行う。そのために、ICP法を用いて
を領域SRi2 0にマッチングする変換パラメータ
を計算し、これをもとに頂点を
と更新し、これよりマッチングを行う。 And update the matching.
次のStep3−3では、領域拡大処理を行う。この領域拡大は図8を用いて説明する。図8(a)において、v"il,k itr+1とv i2,c(k) 0はシード特徴領域内にあり、次の変換にてv"il,k itr+1はv i2,c(k) 0にマッチしている。
次に、図8(b)に示すように、新たにすべての可能なペア
を生成する。各ペアは、
v"il,k itr+1を含むその1近傍隣接頂点集合
v "The set of 1 neighboring neighbors including il, k itr + 1
と、v i2,c(k) 0を含むその1近傍頂点集合
のうちから選択されたものである。この図8の例では、48ペアが生成される。 Is selected from In the example of FIG. 8, 48 pairs are generated.
各頂点間距離ei,sを計算し、その頂点間距離ei,sが閾値τaddlavg以下であれば、図8(c)に示すように、そのv"i1,s itr+1をSR"i1 0に、v i2,c(s) 0をSRi2 0にそれぞれ追加する。そしてこの処理を、条件を満たす頂点が存在しなくなるまで続け、拡大された領域SR"i1 itr+1とSRi2 itr+1を抽出する。通常τadd=1.0に設定するのが好ましい。 If each inter-vertex distance e i, s is calculated and the inter-vertex distance e i, s is less than or equal to the threshold value τ add l avg , as shown in FIG. 8C, the v “ i1, s itr + 1 Is added to SR " i1 0 , and v i2, c (s) 0 is added to SR i2 0 . This process is continued until there are no more vertices satisfying the condition, and the expanded regions SR " i1 itr + 1 and SR i2 itr + 1 are extracted. Usually, it is preferable to set τ add = 1.0.
Step3−4にて終了条件を判断する。もし、|SR"i1 itr+1|>|SR"i1 itr|であれば、itr←itr+1と更新し、Step3−2に戻る。そうでなければ、〈SRi1 itr,SRi2 itr〉をユークリッド対称領域ペアとして抽出し(Step3−5)、さらに、
から最終変換パラメータを
として抽出し、このStep3全体の処理を終了する(Step3−6)。
And the processing of the
続くStep4では、コンパクトな対称性表現の抽出を行う。この手順を図9のフローチャートを用いて説明する。Step3で抽出された数24式の最終変換パラメータを単一平面反射、2−fold回転、単一平行移動といった対称性を表す陽な表現に変換する。
In the
Step4−1,2にて、単一平面反射対称性を判断する。単一平面反射の場合、対称領域ペア〈SRi1 itr,SRi2 itr〉内の各点ペア
の中点はその反射平面上に分布する。そこでこの中点集合に対し最小二乗平面をフィットし、次に領域SRi1 itr内の各点v i1,s itrをこの平面に対して反射させ、反射後の点v i1,s refと対応点v i2,c(s) 0との平均点間距離ei refを計算する。この平均点間距離ei refが閾値threflavg以下であれば、単一平面反射対称性であると判定し、反射平面を抽出する。ここで、全てthref=1.0に設定する。 The midpoints are distributed on the reflection plane. So we fit the least squares plane to this set of midpoints, then reflect each point v i1, s itr in the region SR i1 itr to this plane, and the corresponding points v i1, s ref after reflection v i2, c (s) The average point distance e i ref with 0 is calculated. If this average point distance e i ref is less than or equal to the threshold th ref l avg, it is determined that the single plane reflection symmetry is present, and the reflection plane is extracted. Here, all th ref = 1.0 is set.
同様にStep4−3,4の2−fold回転の場合は、その中点に最小二乗直線をフィットし、その軸周りの回転後の平均点間距離を評価し判定を行い、2−fold回転軸を抽出する。 Similarly, in the case of 2-fold rotation of Steps 4-3 and 4, a least-squares line is fitted to the midpoint, and the average point-to-point distance after rotation around the axis is evaluated and determined, and the 2-fold rotation axis To extract.
Step4−5,6の単一平行移動の場合は、対応点間を結ぶベクトルの平均ベクトルを計算し、そのベクトルに沿った平行移動後の平均点間距離を評価し判定を行い、平行移動ベクトルを抽出する。 In the case of single translation in Steps 4-5 and 6, the average vector of vectors connecting corresponding points is calculated, the average distance between the translated points along the vector is evaluated and determined, and the translation vector is calculated. To extract.
[実施例]
実験結果は図10に示している。図10は、頂点数約12万の、コンロッドのX線CT計測メッシュに対する実験結果を示している。このモデルからは15個の特徴領域が抽出され、平行移動、回転、平面反射の組み合わせから成る20個のユークリッド対称性が認識できた。これらのうち、図10(b)、図10(d)に示す2つの大局的な単一平面反射対称性が認識できた。また、2つの反射平面の法線ベクトル間角度が約89.7degであったことから、本発明方法により、高精度に対称性が認識できることが確認できた。計算時間は、Core2 Duo 2.4GHz CPUで約98秒であった。
[Example]
The experimental results are shown in FIG. FIG. 10 shows experimental results for a connecting rod X-ray CT measurement mesh having about 120,000 vertices. Fifteen feature regions were extracted from this model, and 20 Euclidean symmetries consisting of a combination of translation, rotation, and plane reflection could be recognized. Of these, two global single plane reflection symmetries shown in FIGS. 10B and 10D were recognized. In addition, since the angle between the normal vectors of the two reflection planes was about 89.7 deg, it was confirmed that symmetry could be recognized with high accuracy by the method of the present invention. The calculation time was about 98 seconds with a Core2 Duo 2.4 GHz CPU.
1 入出力制御部
2 データ記憶部
3 メッシュ主曲率算出部
4 特徴領域抽出部
5 シード特徴領域ペア抽出部
6 ICP処理部
7 領域拡大処理部
8 ユークリッド対称性認識処理部
9 演算制御部
51 頂点数比較部
52 マッチング評価部
53 シード特徴領域ペア決定部
100 ユークリッド対称性認識システム
DESCRIPTION OF
Claims (6)
前記データ記憶部の計測メッシュデータに対してメッシュ各頂点の主曲率を計算するメッシュ主曲率算出部と、
前記計測メッシュデータのメッシュ各頂点の主曲率を用いて、平面領域の境界近傍にある頂点集合を特徴領域として抽出する特徴領域抽出部と、
前記特徴領域群から、形状が一致する特徴領域ペアをシード特徴領域ペアとして抽出するシード特徴領域ペア抽出部と、
ICP法とRegion Growing法により同一の変換パラメータでマッピング可能な頂点ペアを逐次的に各シード特徴領域に追加することで特徴領域ペアをそれぞれ反復的に拡大し、互いにユークリッド対称な領域ペアとそれらの変換パラメータを抽出するユークリッド対称領域ペア抽出部と、
前記ユークリッド対称な領域ペアとそれらの変換パラメータを評価し、単一平面反射、2−fold回転若しくは単一平行移動で表現できる領域ペアを見出し、当該領域ペアと共に対称性を表示する対称性表示部とを備えたことを特徴とするユークリッド対称性認識システム。 A data storage unit for storing measurement mesh data and various necessary data;
A mesh principal curvature calculator that calculates a principal curvature of each vertex of the mesh with respect to the measurement mesh data of the data storage unit;
A feature region extraction unit that extracts a set of vertices near the boundary of the planar region as a feature region using the principal curvature of each vertex of the mesh of the measurement mesh data;
From the feature region group, a seed feature region pair extraction unit that extracts a feature region pair having a matching shape as a seed feature region pair;
By adding vertex pairs that can be mapped with the same transformation parameters by the ICP method and Region Growing method to each seed feature region sequentially, each feature region pair is expanded repetitively, and Euclidean symmetric region pairs and their pairs A Euclidean symmetric region pair extractor for extracting transformation parameters;
A symmetry display unit that evaluates the Euclidean symmetric region pairs and their conversion parameters, finds a region pair that can be expressed by single plane reflection, 2-fold rotation, or single translation, and displays symmetry along with the region pair A Euclidean symmetry recognition system characterized by comprising:
前記各特徴領域ペアの頂点数を比較する頂点数比較部と、
前記頂点数の比較にて所定の範囲内にある特徴領域ペアについて、重心・主軸一致によるマッチングとその精度を評価するマッチング評価部と、
ICP法によりマッチングし、その精度を評価して所定の範囲内にある特徴領域ペアについてシード特徴領域ペアとして抽出するシード特徴領域ペア決定部とを有することを特徴とする請求項1に記載のユークリッド対称性認識システム。 The seed feature region pair extraction unit
A vertex number comparison unit for comparing the number of vertices of each feature region pair;
A matching evaluation unit that evaluates the matching by the centroid and principal axis matching and the accuracy of the feature region pair within a predetermined range in the comparison of the number of vertices,
The Euclidean according to claim 1, further comprising: a seed feature region pair determination unit that performs matching by an ICP method, evaluates accuracy thereof, and extracts a feature region pair within a predetermined range as a seed feature region pair. Symmetry recognition system.
前記メッシュ各頂点の主曲率を用いて、平面領域の境界近傍にある頂点集合を特徴領域として抽出するステップと、
特徴領域群に対して、形状が一致する特徴領域ペアをシード特徴領域ペアとして抽出するステップと、
ICP法とRegion Growing法を組み合わせて、同一の変換パラメータでマッピング可能な頂点ペアを逐次的に各シード特徴領域に追加することで領域ペアをそれぞれ反復的に拡大し、互いにユークリッド対称な領域ペアとそれらの変換パラメータを抽出するステップと、
前記ステップにて抽出した変換パラメータを評価し、単一平面反射、2−fold回転若しくは単一平行移動の対称性表現できる領域ペアを見出し、当該領域ペアと共に対称性を表示するステップとを有することを特徴とするユークリッド対称性認識方法。 Calculating the principal curvature of each vertex of the mesh;
Using the principal curvature of each vertex of the mesh, extracting a vertex set in the vicinity of the boundary of the planar region as a feature region;
Extracting a feature region pair having a matching shape as a seed feature region pair for the feature region group; and
By combining the ICP method and the Region Growing method, vertex pairs that can be mapped with the same transformation parameters are sequentially added to each seed feature region, so that each region pair is expanded repeatedly, and Euclidean symmetric region pairs are obtained. Extracting those conversion parameters;
Evaluating the transformation parameters extracted in the above step, finding a region pair that can express symmetry of single plane reflection, 2-fold rotation, or single translation, and displaying symmetry along with the region pair. Euclidean symmetry recognition method characterized by
前記各特徴領域ペアの頂点数を比較するステップと、
前記頂点数の比較にて所定の範囲内にある特徴領域ペアについて、重心・主軸一致によるマッチングとその精度を評価するステップと、
ICP法によりマッチングし、その精度を評価して所定の範囲内にある特徴領域ペアについてシード特徴領域ペアとして抽出するステップとを有することを特徴とする請求項3に記載のユークリッド対称性認識方法。 For the feature region group, extracting a feature region pair having a matching shape as a seed feature region pair,
Comparing the number of vertices of each feature region pair;
For the feature region pair within a predetermined range in the comparison of the number of vertices, the step of evaluating the matching by the center-of-gravity / main axis matching and its accuracy;
The method according to claim 3, further comprising: matching by an ICP method, evaluating accuracy thereof, and extracting a feature region pair within a predetermined range as a seed feature region pair.
前記メッシュ各頂点の主曲率を用いて、平面領域の境界近傍にある頂点集合を特徴領域として抽出するステップと、
特徴領域群に対して、形状が一致する特徴領域ペアをシード特徴領域ペアとして抽出するステップと、
ICP法とRegion Growing法を組み合わせて、同一の変換パラメータでマッピング可能な頂点ペアを逐次的に各シード特徴領域に追加することで領域ペアをそれぞれ反復的に拡大し、互いにユークリッド対称な領域ペアとそれらの変換パラメータを抽出するステップと、
前記ステップにて抽出した変換パラメータを評価し、単一平面反射、2−fold回転若しくは単一平行移動の対称性表現できる領域ペアを見出し、当該領域ペアと共に対称性を表示するステップとをコンピュータに実行させることを特徴とするユークリッド対称性認識プログラム。 Calculating the principal curvature of each vertex of the mesh;
Using the principal curvature of each vertex of the mesh, extracting a vertex set in the vicinity of the boundary of the planar region as a feature region;
Extracting a feature region pair having a matching shape as a seed feature region pair for the feature region group; and
By combining the ICP method and the Region Growing method, vertex pairs that can be mapped with the same transformation parameters are sequentially added to each seed feature region, so that each region pair is expanded repeatedly, and Euclidean symmetric region pairs are obtained. Extracting those conversion parameters;
Evaluating the transformation parameter extracted in the above step, finding a region pair capable of expressing symmetry of single plane reflection, 2-fold rotation, or single translation, and displaying the symmetry together with the region pair to a computer A Euclidean symmetry recognition program characterized by being executed.
前記各特徴領域ペアの頂点数を比較するステップと、
前記頂点数の比較にて所定の範囲内にある特徴領域ペアについて、重心・主軸一致によるマッチングとその精度を評価するステップと、
ICP法によりマッチングし、その精度を評価して所定の範囲内にある特徴領域ペアについてシード特徴領域ペアとして抽出するステップとをコンピュータに実行させることを特徴とする請求項5に記載のユークリッド対称性認識プログラム。 For the feature region group, extracting a feature region pair having a matching shape as a seed feature region pair,
Comparing the number of vertices of each feature region pair;
For the feature region pair within a predetermined range in the comparison of the number of vertices, the step of evaluating the matching by the center of gravity / main axis matching and its accuracy;
6. The Euclidean symmetry according to claim 5, further comprising: causing the computer to perform matching by an ICP method, evaluating accuracy thereof, and extracting a feature region pair within a predetermined range as a seed feature region pair. Recognition program.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008174641A JP2010014541A (en) | 2008-07-03 | 2008-07-03 | System, method and program for recognizing euclid symmetry |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008174641A JP2010014541A (en) | 2008-07-03 | 2008-07-03 | System, method and program for recognizing euclid symmetry |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2010014541A true JP2010014541A (en) | 2010-01-21 |
Family
ID=41700784
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008174641A Ceased JP2010014541A (en) | 2008-07-03 | 2008-07-03 | System, method and program for recognizing euclid symmetry |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2010014541A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011180786A (en) * | 2010-02-26 | 2011-09-15 | Mitsubishi Heavy Ind Ltd | Cad data transmission method, cad data reception method, cad data transmission device, cad data reception device, cad data transmission program, cad data reception program and data structure |
JP2011198330A (en) * | 2010-03-24 | 2011-10-06 | National Institute Of Advanced Industrial Science & Technology | Method and program for collation in three-dimensional registration |
JP2013529819A (en) * | 2010-06-30 | 2013-07-22 | トムソン ライセンシング | Method and apparatus for detecting repetitive structures in a 3D mesh model |
EP4086849A1 (en) | 2021-12-24 | 2022-11-09 | Univerza v Mariboru | Procedure for detecting the reflectional symmetry of geometrical objects |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007241996A (en) * | 2006-02-13 | 2007-09-20 | Hokkaido Univ | Apparatus, method and program for segmentation of analytic curved surface and recording medium |
-
2008
- 2008-07-03 JP JP2008174641A patent/JP2010014541A/en not_active Ceased
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007241996A (en) * | 2006-02-13 | 2007-09-20 | Hokkaido Univ | Apparatus, method and program for segmentation of analytic curved surface and recording medium |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011180786A (en) * | 2010-02-26 | 2011-09-15 | Mitsubishi Heavy Ind Ltd | Cad data transmission method, cad data reception method, cad data transmission device, cad data reception device, cad data transmission program, cad data reception program and data structure |
JP2011198330A (en) * | 2010-03-24 | 2011-10-06 | National Institute Of Advanced Industrial Science & Technology | Method and program for collation in three-dimensional registration |
JP2013529819A (en) * | 2010-06-30 | 2013-07-22 | トムソン ライセンシング | Method and apparatus for detecting repetitive structures in a 3D mesh model |
US10311635B2 (en) | 2010-06-30 | 2019-06-04 | Interdigital Vc Holdings, Inc. | Method and apparatus for detecting repetitive structures in 3D mesh models |
EP4086849A1 (en) | 2021-12-24 | 2022-11-09 | Univerza v Mariboru | Procedure for detecting the reflectional symmetry of geometrical objects |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Maiseli et al. | Recent developments and trends in point set registration methods | |
JP4924165B2 (en) | Shape comparison device and shape comparison method | |
US11017210B2 (en) | Image processing apparatus and method | |
US6044171A (en) | Method and apparatus for pattern recognition and representation using fourier descriptors and iterative transformation-reparametrization | |
EP2589026B1 (en) | Method and apparatus for detecting repetitive structures in 3d mesh models | |
JP6069489B2 (en) | Object recognition apparatus, object recognition method, and program | |
CN109325490B (en) | Terahertz image target identification method based on deep learning and RPCA | |
CN112102294B (en) | Training method and device for generating countermeasure network, and image registration method and device | |
JP2011138388A (en) | Data correction apparatus and method | |
JP5448758B2 (en) | Information processing apparatus, information processing method, information processing system, program, and data structure | |
CN109426828B (en) | Point cloud classification method, device, equipment and storage medium | |
US20030235337A1 (en) | Non-rigid image registration using distance functions | |
CN113095333B (en) | Unsupervised feature point detection method and unsupervised feature point detection device | |
JP2012507793A (en) | Complexity normalization pattern representation, search, and compression | |
KR20160053833A (en) | Improved data comparison method | |
Shao et al. | 3D dynamic facial expression recognition using low-resolution videos | |
JP2010014541A (en) | System, method and program for recognizing euclid symmetry | |
Ma et al. | Kinematic skeleton extraction from 3D articulated models | |
WO2007098579A1 (en) | Method and system for locating landmarks on 3d models | |
KR20170100718A (en) | Apparatus and method for registration of point cloud | |
JP2019194788A (en) | Learning device, recognition device, learning method and computer program | |
Abbas et al. | Automatic classification of facial morphology for medical applications | |
JP6453618B2 (en) | Calculation apparatus, method and program | |
Li et al. | Near-convex decomposition of 2D shape using visibility range | |
US20230143070A1 (en) | Learning device, learning method, and computer-readable medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110629 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20110705 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20121129 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20121204 |
|
A045 | Written measure of dismissal of application [lapsed due to lack of payment] |
Free format text: JAPANESE INTERMEDIATE CODE: A045 Effective date: 20130423 |