Research of 3D Irregular Fragment Reassembly Technique in Reverse Engineering
Total Page:16
File Type:pdf, Size:1020Kb
International Journal of Research in Engineering and Science (IJRES) ISSN (Online): 2320-9364, ISSN (Print): 2320-9356 www.ijres.org Volume 4 Issue 5 ǁ May. 2016 ǁ PP. 06-10 Research of 3D Irregular Fragment Reassembly Technique in Reverse Engineering NIE Bo-lin, ZHANG Xu, CHE Xuan-lin (School of Mechanical Engineering, Shanghai University of Engineering Science, Shanghai, 201620) Abstract:- in order to achieve some pieces of important research value of restructuring, this paper aiming at the shortcomings of the current existing algorithm proposed a new spatial 3 d irregular fragments stitching method. First Uses the method to establish point cloud based on kd - tree search space the topology relationship, for k neighborhood search quickly, further realize the boundary extraction of point cloud fragments; Secondly using Nurbs curve interpolation method to deal with boundary feature points, complete boundary fit at the same time the curve interpolation points of curvature and torsion; Finally based on the calculation of curvature and torsion, according to the longest common subsequence complete curve matching, and then complete the point cloud fragments of stitching. Through examples show that the proposed algorithm good robustness, high matching precision. Keywords:- Reverse Engineering ;Curve fitting ;Curve matching ;The Longest Common uence (LCS) ; fragment reassembly I. INTRODUCTION So-called 3 d fragments split is some irregular space debris joining together into a complete model of the initial. In our country, ceramic relics buried in the ground, because of geographical or man-made reasons, the unearthed will be destroyed into large and small pieces. As a result of these pieces of cultural relics still has the very high archaeological value, therefore how to use modern technology to restore these cultural relics has important research significance. Currently, the main space debris in the fracture restoration contour curve matching, as part of the curve matching problem. In recent years, many Chinese and foreign scholars on curve matching this sort of question to do a lot of research., the earliest Besl[1] and proposed recently to 3 d space curve matching point iteration (iterative closest point, ICP) method, but as a result of this method requires an empty a subset of the curve is another empty asked curve, but in pieces joining together, the two can match the intersection, the space curve can be not as the empty set, rather than an outline is another subset of contour, so this method is difficult to directly used for contour curve matching of fragments merging. Ucoluk G[2] the longest common subsequence algorithm is proposed, such as sampling points on the curve curvature and torsion as the characteristic, according to the similar feature points of similar matrix, using the dynamic programming algorithm to find the longest matching sequence, as the optimal matching, the algorithm deficiency is high time complexity, the ideal smooth curve matching effect is good, does not apply to any free curve shape. Wolfson[3] to curve resampling in the first place, such as computing the curvature of the sample point, the curvature to characters value, and then according to the matching hash algorithm. Steven E[4] according to the curvature of the contour curve feature points such as value, take the Pearson coefficient calculation according to the similarity of two curves. Shin H[5] through statistics such as mural contour curve geometric features include information such as length, area, four crown, compare their similarity. Oxholm G[6] computing the curvature of the contour curve of sampling points, such as burning rate and the color of character string, to take the longest common substring realize contour matching. Ding Xianfeng[7] on existing two-dimensional curve matching algorithm to carry on the induction and summary. Pan Rongjiang[8] etc is put forward based on the longest common subsequence (LCS) fragments of the contour curve matching algorithm is used to match, according to the characteristics of point contour curve segment, and find out the LCS between two piecewise curve, the curve of the spell of overlapping tests results, it is concluded that the optimal matching. Lv Ke etc[9] will outline curve segment, and put forward the measure curve segment similarity hash vector, using contour line segment matching algorithm based on Fourier transform, by comparing the two hash of the profile of the vector for the analysis of curve segment is similar. Jian wang[10], etc according to the curvature of the sampling points get two curves on match point, through the alignment of the matching point Frenet frame curve matching. Shu-cheng zhou[11] presents a using wavelet multiscale description of contour curve matching algorithm. Zhu Yanjuan etc[12] are given based on the sampling point of curvature and torsion of the contour curve matching algorithm, according to total curvature to extract www.ijres.org 6 | Page Research Of 3d Irregular Fragment Reassembly Technique In Reverse Engineering the feature points, and the contour curve segment, according to the feature point on the piecewise curve curvature and torsion judgment piecewise curves is similar, and then according to the method of vector for further verification, the curve of the high similarity of 3 d transform at the same time, they are aligned, so as to realize the splice pieces. II. FRAGMENT REASSEMBLY 2.1 boundary extraction and boundary fitting 2.1.1 boundary extraction This article adopts the method of tree search based on kd - a point cloud space the topology relationship, k neighborhood search quickly. In the mining point P and the search for k neighboring points Mj (j = 0,..., k - 1) as a local reference data, using the least squares fitting out the tangent plane. Assume that F (x, y, z) = a1 + a2 + a3 y z + a4 x = 0, for the micro tangent plane of k + 1 point, written in matrix form is: x0 y0 z0 1 a1 x y z 1 a 1 1 1 2 = 0 (1) a3 xk yk zk 1 a4 T ,a=[a1,a2,a3,a4] 。So the Aa = 0.Singular value A= decomposition of ATA: ∆ 0 퐴 = 푈 푉퐻 (2) 0 0 Type in the U and V as the unitary matrix, ∆= 푑푖푎푔[ 휆1, 휆2, ⋯ , 휆푟 ] and the singular value of matrix A T is 휆푖 (i=1,2,3,…,r), which are characteristic values for A A. r is the number of singular values, and that number is eigenvalue of the matrix. ATA eigenvectors corresponding to the minimum eigenvalue is the least squares solution of Aa = 0.The micro tangent plane normal vector is n = (a1, a2, a3). The sampling points and their neighborhood K corresponding to the projection plane microdissection to ′ ′ obtain the projected point 푃푖 and 푀푖 (j = 0,1,2, ⋯, k -1), the point set scattered into general distribution. Calculated projection coordinates, assuming space unorganized points coordinates Ni (xi, yi, zi) (i = 0,1, ⋯, n), ′ ′ ′ which coordinates the micro-projection plane is cut (푥푖 , 푦푖 , 푧푖 ), so there are: ′ ′ ′ 푎1푥푖 + 푎2푦푖 + 푎3푧푖 + 푎4 = 0 (3) ′ ′ ′ 푥푖 − 푥푖 푦푖 − 푦푖 푧푖 − 푧푖 = = = 푘 (4) 푎1 푎2 푎3 In the equation, 푎 푥 + 푎 푦 + 푎 푧 푘 = 1 푖 2 푖 3 푖 2 2 2 푎1 + 푎2 + 푎3 ′ ′ ′ Simultaneous equations (1) and (2) can be obtained projection point coordinates (푥푖 , 푦푖 , 푧푖 ). Construction xP'x plane coordinate system, the micro-cut plane to P'i origin to P'i to the vector direction from the farthest point P'i posed for the x-axis direction, perpendicular to the x-axis direction the y direction, the projection on the tangent plane of the micro-sample point, if its K neighbors projection points evenly distributed around the sample point, then the sampling point non-boundary characteristic points; if K neighborhood projection points around the sampling point unevenly distributed, it is considered that the sampling point for the [13] boundary characteristic point . In P'i as a starting point for its K points in the neighborhood on the micro- projection tangent plane as the end point defined vectors P'iM'j, the three-dimensional coordinates of the projection point conversion to cut the micro plane, and the plane parameters microdissection of. Projection points x, y coordinates are in P'iM'j vector corresponding to the projected length Mjx and Mjy of x-axis and y- axis, thereby obtaining projection point coordinates M'j the parameterized (Mjx, Mjy). Since Scattered clouds in the form of discrete spatial distribution of vector projected point is generally composed of random distribution, in order to facilitate the calculation, the first projection vector is normalized: www.ijres.org 7 | Page Research Of 3d Irregular Fragment Reassembly Technique In Reverse Engineering 푃′푀′ 푀′ , 푀′ ′ " 푖 푗 푗푋 푗푌 푃푖 푀푗 = ′ " = = 푋푖 , 푌푗 (5) 푃푖 푀푗 ′ 2 ′ 2 푀푗푋 + 푀푗푌 ′ " In the equation, j =0,1,2,⋯,k-1. and j≠i. For vector 푃푖 푀푗 , giving it a force 퐹푖푗 , the direction of the force vector in ′ " the same direction, the size of the force is equal to | 푃푖 푀푗 |, then by calculating the K points in the neighborhood of sampling points for force 퐹푖푗 in component A and B on X and Y axes to identify the boundaries of the feature point, component calculated by the following equation: 퐾−1 퐾−1 퐾−1 ′ " 퐹푖푗 = 푃푖 푀푗 = 푋푗 (6) 푗 =0 푗 =0 푗 =0 푗 ≠푖 푗 ≠푖 푗 ≠푖 푋 푋 퐾−1 퐾−1 퐾−1 ′ " 퐹푖푗 = 푃푖 푀푗 = 푌푗 (7) 푗 =0 푗 =0 푗 =0 푗 ≠푖 푗 ≠푖 푗 ≠푖 푌 푌 First, set a threshold σ, when 퐾−1 푗 =0 퐹푖푗 푗 ≠푖 푋 > 휎 퐾 or 퐾−1 푗 =0 퐹푖푗 푗 ≠푖 푌 > 휎 퐾 it can be determined that the sample point as a boundary point.