CA2564156C - System and method for approximating an editable surface - Google Patents
System and method for approximating an editable surface Download PDFInfo
- Publication number
- CA2564156C CA2564156C CA2564156A CA2564156A CA2564156C CA 2564156 C CA2564156 C CA 2564156C CA 2564156 A CA2564156 A CA 2564156A CA 2564156 A CA2564156 A CA 2564156A CA 2564156 C CA2564156 C CA 2564156C
- Authority
- CA
- Canada
- Prior art keywords
- grid
- value
- point
- points
- new
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 80
- 238000007670 refining Methods 0.000 claims description 16
- 238000009499 grossing Methods 0.000 claims description 11
- 238000003384 imaging method Methods 0.000 claims description 9
- 238000009877 rendering Methods 0.000 claims description 7
- 230000006870 function Effects 0.000 description 15
- 238000002156 mixing Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000000844 transformation Methods 0.000 description 3
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 241000933095 Neotragus moschatus Species 0.000 description 1
- BQCADISMDOOEFD-UHFFFAOYSA-N Silver Chemical compound [Ag] BQCADISMDOOEFD-UHFFFAOYSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000135 prohibitive effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
Abstract
A system and method are disclosed for automatically approximating an editable surface from a 3D data set or a 3D point set, which may be imaged in the form of a NURBS surface.
Description
SYSTEM AND METHOD FOR
APPROXIMATING AN EDITABLE SURFACE
FIELD OF THE INVENTION
The present invention generally relates to systems and methods for automatically approximating an editable surface from a 3D data set or 3D point set, which may be imaged in the form of a NURBS surface.
BACKGROUND OF THE INVENTION
Non-Uniform Rational B-Splines (NURBS) are industry standard tools for the representation and design of geometry. NURBS, as explained by Markus Altmann at www.cs.wpi.edu/-matt/courses/cs563/talks/nurbs.html, are used for a variety of reasons.
They offer one common mathematical form for both standard analytical shapes (e.g., conics) and free-form shapes. NURBS provide the flexibility to design a large variety of shapes and can be evaluated reasonably fast by numerically stable and accurate algorithms. They are invariant under affine as well as perspective transformations and are generalizations of non-rational B-Splines and non-rational and rational Bezier curves and surfaces.
However, one of the drawbacks of NURBS is the need for extra storage to define traditional shapes (e.g., circles). This results from parameters in addition to the control points, but will allow the desired flexibility for defining parametric shapes.
NURBS-shapes are not only defined by control points; weights, associated with each control point, are also necessary. A NURBS curve C(u), for example, which is a vector-valued piecewise rational polynomial function, may be defined as:
APPROXIMATING AN EDITABLE SURFACE
FIELD OF THE INVENTION
The present invention generally relates to systems and methods for automatically approximating an editable surface from a 3D data set or 3D point set, which may be imaged in the form of a NURBS surface.
BACKGROUND OF THE INVENTION
Non-Uniform Rational B-Splines (NURBS) are industry standard tools for the representation and design of geometry. NURBS, as explained by Markus Altmann at www.cs.wpi.edu/-matt/courses/cs563/talks/nurbs.html, are used for a variety of reasons.
They offer one common mathematical form for both standard analytical shapes (e.g., conics) and free-form shapes. NURBS provide the flexibility to design a large variety of shapes and can be evaluated reasonably fast by numerically stable and accurate algorithms. They are invariant under affine as well as perspective transformations and are generalizations of non-rational B-Splines and non-rational and rational Bezier curves and surfaces.
However, one of the drawbacks of NURBS is the need for extra storage to define traditional shapes (e.g., circles). This results from parameters in addition to the control points, but will allow the desired flexibility for defining parametric shapes.
NURBS-shapes are not only defined by control points; weights, associated with each control point, are also necessary. A NURBS curve C(u), for example, which is a vector-valued piecewise rational polynomial function, may be defined as:
(1) C(u) sum(i = 0, n){wi* Pi* Ni,k(u)}
=
sum(i = 0, n){wi* M,k(u)}
where w1 = weights Pi = control points (vector) = normalized B-spline basis functions of degree k.
The B-Splines are defined recursively as:
u 1V , ti+k+i-u (2) Ni,k(u)= ____ i,k ¨1(0+ ______________ ti + k ¨ti ti+k+1¨ti+1 and / 1, if ti u < ti+i Ni,o(u) =<
\ 0, else where ti are the knots forming a knot vector and U = {to, ti, ,t.}.
The Knot Vector The knot vector uniquely determines the B-Splines as demonstrated above relative to equation (2). The relation between the number of knots (m+1), the degree (k) of Ni,ic and the number of control points (n+1) is given by m = n + k + 1.
The sequence of knots in the knot vector U is assumed to be non-decreasing, i.e., ti <= ti+1. Each successive pair of knots represents an interval [ti, t1+1) for the parameter values to calculate a segment of a shape.
For NURBS, the relative parametric intervals (knot spans) need not be the same for all shape segments, i.e., the knot spacing is non-uniform, leading to a non-periodic knot vector of the form:
=
sum(i = 0, n){wi* M,k(u)}
where w1 = weights Pi = control points (vector) = normalized B-spline basis functions of degree k.
The B-Splines are defined recursively as:
u 1V , ti+k+i-u (2) Ni,k(u)= ____ i,k ¨1(0+ ______________ ti + k ¨ti ti+k+1¨ti+1 and / 1, if ti u < ti+i Ni,o(u) =<
\ 0, else where ti are the knots forming a knot vector and U = {to, ti, ,t.}.
The Knot Vector The knot vector uniquely determines the B-Splines as demonstrated above relative to equation (2). The relation between the number of knots (m+1), the degree (k) of Ni,ic and the number of control points (n+1) is given by m = n + k + 1.
The sequence of knots in the knot vector U is assumed to be non-decreasing, i.e., ti <= ti+1. Each successive pair of knots represents an interval [ti, t1+1) for the parameter values to calculate a segment of a shape.
For NURBS, the relative parametric intervals (knot spans) need not be the same for all shape segments, i.e., the knot spacing is non-uniform, leading to a non-periodic knot vector of the form:
(3) U = { a, ... , a, tk+19 ... ,tmk-1, ,b }, where a and b are repeated with multiplicity of k+1. The multiplicity of a knot affects the parametric continuity at this knot. Non-periodic B-Splines, like NURBS, are infinitely and continuously differentiable in the interior of a knot span and k-M-1 times continuously differentiable at a knot, where M is the multiplicity of the knot. In contrast, a periodic knot vector U = {0, 1, , n} is everywhere k-1 times continuously differentiable. Considering the knot vector for NURBS, the end knot points (tk, tn4-1) with multiplicity k+1 coincide with the end control points PO, P.
Since the knot spacing could be non-uniform, the B-Splines are no longer the same for each interval [ti, ti+1) and the degree of the B-Spline may vary.
Considering the whole range of parameter values represented by the knot vector, the different B-Splines build up continuous (overlapping) blending functions Ni,k(u), as defined in equation (2), over this range of parameter values, which are illustrated in Table 1 below.
These blending functions have the following properties:
1. Ni,k(u) >= 0, for all i, k, u;
2. Ni,k(u) = 0, if u not in [ti, ti+k+i), meaning local support of k+1 knot spans, where Ni,k(u) is nonzero;
3. If u in [ti, ti+1), the non-vanishing blending functions are Ni-k,k(u), 4. Sum (j=i-k, i){Ni,k0.0} = sum(i=0, n){Ni,k(u)} = 1, (partition of unity);
and 5. In case of multiple knots, 0/0 is deemed to be zero.
The first and fourth properties, as illustrated in Table 2 below, together, result in a convex hull. The polyline of control points build up to a shape defined by a NURBS
curve. The second and third properties illustrate that k+1 successive control points define a shape segment, and a control point is involved in k+1 neighboring shape segments. Therefore, changing a control point or weight influences just k+1 shape segments, defined over the interval given in equation (2).
IP
tit IP ' .
(\illiri A16144-- -'441112' Curve/Surface Definition The previous definition of a NURBS-curve in equation (1) may be rewritten using rational basis functions:
11),* Nik(u) (4) Rt,k(zt)-= _______________ SUM(j -= 0, n){w,* NJ,k(u)}
into:
(5) C(u) = suni(i = 0, n){13,* 1?,,k(u)} .
A NURBS-surface may be defined in a similar way:
S(u, v) = sum(i = 0, n)sum(j = 0, in) Pij * Ri,k, v) where wij * Ni,k(u) * Nj,l(v) (6) Ri,k, j,l(u,v)=
sum(r = 0, n){sum(s = 0, m){wr,s * Nr,k(u) * Nsi(u)}
The rational basis functions have the same properties as the blending functions.
One point to emphasize, is their invariance under affine and (even) perspective transformations. Therefore, only the control points have to be transformed to get the appropriate transformation of the NURBS shape.
Computational Algorithm NURBS can be evaluated effectively by using homogeneous coordinates. The following steps demonstrate one method to perform the evaluation:
1. Add one dimension to the control points (e.g., P = (x, y) -> P'(x, y, 1)) and multiply them by their corresponding weights, i.e., in 2D: Pi(xi, yi) Pi(Wi *
x, w1 * yi, wi) 2. Calculate NURBS in homogeneous coordinates:
C'(u) = sum(i = 0, n){Pi,(u) * Ni,k(u)}
3. Map "homogeneous" NURBS back to original coordinate system with:
/ ( X1/W, X2/W, , Xn/W ), if W not = 0 map( X1 , X2, ... ,Xn, W) = <
\ ( X1 , X2, , Xn ), if W = 0 = Pi (7) C(u) = map(C '(u)) = sum(i 0, n){wi * *
sum(i = 0, n){wi * Ni,k(0}
For u in [ti, ti+i), the only existing blending functions to consider in evaluation of the curve at u are Ni-k,k(u), u.An effective algorithm for the computation of the non-vanishing blending functions is described in C. deBoor, "A Practical Guide to Splines,"
1978, New York, Springer-Verlag.
The Weights As mentioned above, changing the weight wi of a control point Pi affects only the range [ti, ti+k+i) (in case of a curve). The geometric meaning of the weights is illustrated in Table 3, below.
Defining the points:
B = C(u; = 0);
N = C(u; w= 1); and Bi = C(u; wi not = {0, 1}).
P2 _______________________________________________________ ,r71171111B
N and B, can also be expressed as:
N = (1 - a) * B + a * P, B, = (1 - b)*B+b*P1, where a = R,,k(u; w = 1) b = Ri,k(u) The following identity is obtained from the expression of a and b:
(1 - a)I a : (1 - b)I b = P,N/BN : P,B,/BB, = w1, which is called the cross- or double ratio of the four points Põ B, N, B,.
From these expressions, the effect of shape modification can be derived:
= B, sweeps out on a straight line segment;
= If w, = 0 then Pi has no effect on shape;
= If w, increases, so b and the curve is pulled toward P, and pushed away from PJ, for j not= i;
= If w, decreases, so b and the curve is pushed away from P, and pulled toward PJ, for j not= i ; and = If w, -> infinity then b -> 1 and B, -> Põ if u in [ti, trEk+i) The Problem Various techniques have been attempted for creating a NURBS surface in different fields of art. For example, International Publication No. WO
02/37422, and U.S. Patent No. 6,765,570, propose various techniques for the manual generation of editable (NURBS) surfaces used in analyzing and interpreting seismic events.
Other conventional applications propose creating an editable (NURBS) surface using interpolation techniques well known in the art. Such techniques may be referred to as an "exact-fit" approach to determining the editable surface. An exact-fit approach is more likely to render a less uniform, if not non-uniform, editable surface as compared to a "best-fit" approximation of the editable surface.
Moreover, the exact-fit approach to defining an editable NURBS surface may be impractical, if not cost prohibitive, for large 3D data sets often encountered in the analysis of seismic data.
Other conventional methods for converting unordered points to surfaces are generally described in Shepard, D., "A two dimensional interpolation function for irregular spaced data," 1968, pp. 517-524, Proceedings 23rd ACM National Conference.
One well-known method for interpolating scattered points, defines a surface value based on given points and weight functions. Although the result is relatively rough, the fundamental idea has inspired many other methods.
Another well-known method, generally referred to as the "Thin Plate Spline"
method has been incorporated in medical imaging applications and is capable of precisely guiding the target surface such that it passes through all given points. See Hardy, R. L. Desmarrais, R.N., "Interpolation using surface spline," 1972, pp.
189-197, Journal of Aircraft 9 and Dyn, N., "Interpolation in Scattered Data by Radial Functions,"
1987, pp. 47-61, In Chui, C.K.; Schumaker, L.L.; Ultreres, F.I. (ed) Topics in Multivariate Approximation. However, this technique requires the inversion of large matrices, which is computationally expensive and generally impractical for a large number of input points that may only require a best-fit approximation.
Another conventional method for converting a point cloud to a surface is described in Hoppe, H., et al., "Surface Reconstruction from Unorganized Points, 1992, pp. 71-78, Comput. Graph. 26. However, this method assumes that the input points are evenly distributed over the entire domain, which is also impractical for input points that are densely populated in some areas and non-existent in other areas.
A need therefore, exists for automatically approximating an editable surface from a 3D data set or 3D point set comprising a large volume of unordered and/or unstructured data points, which may be imaged in the form of an editable NURBS surface.
SUMMARY OF THE INVENTION
The present invention meets the above needs and overcomes one or more deficiencies in the prior by providing systems and methods for automatically approximating an editable surface from a 3D data set or a 3D point set that may comprise a large volume of unstructured and/or unordered data points.
In one embodiment, the present invention includes a system for approximating an editable surface from a 3D data set comprising a computer-readable memory medium configured to store a program of instructions capable of being executable to implement:
i) selecting a point set from the 3D data set; ii) determining a best-fit plane for the point set; iii) projecting at least a portion of the point set on to the best-fit plane; iv) determining a boundary for the best-fit plane according to the projected point set; v) projecting a grid on to the best-fit plane within the boundary, the grid containing a plurality of grid points; vi) initializing the grid points; vii) determining a value for at least a portion of the grid points; viii) selecting at least a portion of the grid points with a value; and ix) imaging the editable surface using at least a portion of the selected grid points.
Since the knot spacing could be non-uniform, the B-Splines are no longer the same for each interval [ti, ti+1) and the degree of the B-Spline may vary.
Considering the whole range of parameter values represented by the knot vector, the different B-Splines build up continuous (overlapping) blending functions Ni,k(u), as defined in equation (2), over this range of parameter values, which are illustrated in Table 1 below.
These blending functions have the following properties:
1. Ni,k(u) >= 0, for all i, k, u;
2. Ni,k(u) = 0, if u not in [ti, ti+k+i), meaning local support of k+1 knot spans, where Ni,k(u) is nonzero;
3. If u in [ti, ti+1), the non-vanishing blending functions are Ni-k,k(u), 4. Sum (j=i-k, i){Ni,k0.0} = sum(i=0, n){Ni,k(u)} = 1, (partition of unity);
and 5. In case of multiple knots, 0/0 is deemed to be zero.
The first and fourth properties, as illustrated in Table 2 below, together, result in a convex hull. The polyline of control points build up to a shape defined by a NURBS
curve. The second and third properties illustrate that k+1 successive control points define a shape segment, and a control point is involved in k+1 neighboring shape segments. Therefore, changing a control point or weight influences just k+1 shape segments, defined over the interval given in equation (2).
IP
tit IP ' .
(\illiri A16144-- -'441112' Curve/Surface Definition The previous definition of a NURBS-curve in equation (1) may be rewritten using rational basis functions:
11),* Nik(u) (4) Rt,k(zt)-= _______________ SUM(j -= 0, n){w,* NJ,k(u)}
into:
(5) C(u) = suni(i = 0, n){13,* 1?,,k(u)} .
A NURBS-surface may be defined in a similar way:
S(u, v) = sum(i = 0, n)sum(j = 0, in) Pij * Ri,k, v) where wij * Ni,k(u) * Nj,l(v) (6) Ri,k, j,l(u,v)=
sum(r = 0, n){sum(s = 0, m){wr,s * Nr,k(u) * Nsi(u)}
The rational basis functions have the same properties as the blending functions.
One point to emphasize, is their invariance under affine and (even) perspective transformations. Therefore, only the control points have to be transformed to get the appropriate transformation of the NURBS shape.
Computational Algorithm NURBS can be evaluated effectively by using homogeneous coordinates. The following steps demonstrate one method to perform the evaluation:
1. Add one dimension to the control points (e.g., P = (x, y) -> P'(x, y, 1)) and multiply them by their corresponding weights, i.e., in 2D: Pi(xi, yi) Pi(Wi *
x, w1 * yi, wi) 2. Calculate NURBS in homogeneous coordinates:
C'(u) = sum(i = 0, n){Pi,(u) * Ni,k(u)}
3. Map "homogeneous" NURBS back to original coordinate system with:
/ ( X1/W, X2/W, , Xn/W ), if W not = 0 map( X1 , X2, ... ,Xn, W) = <
\ ( X1 , X2, , Xn ), if W = 0 = Pi (7) C(u) = map(C '(u)) = sum(i 0, n){wi * *
sum(i = 0, n){wi * Ni,k(0}
For u in [ti, ti+i), the only existing blending functions to consider in evaluation of the curve at u are Ni-k,k(u), u.An effective algorithm for the computation of the non-vanishing blending functions is described in C. deBoor, "A Practical Guide to Splines,"
1978, New York, Springer-Verlag.
The Weights As mentioned above, changing the weight wi of a control point Pi affects only the range [ti, ti+k+i) (in case of a curve). The geometric meaning of the weights is illustrated in Table 3, below.
Defining the points:
B = C(u; = 0);
N = C(u; w= 1); and Bi = C(u; wi not = {0, 1}).
P2 _______________________________________________________ ,r71171111B
N and B, can also be expressed as:
N = (1 - a) * B + a * P, B, = (1 - b)*B+b*P1, where a = R,,k(u; w = 1) b = Ri,k(u) The following identity is obtained from the expression of a and b:
(1 - a)I a : (1 - b)I b = P,N/BN : P,B,/BB, = w1, which is called the cross- or double ratio of the four points Põ B, N, B,.
From these expressions, the effect of shape modification can be derived:
= B, sweeps out on a straight line segment;
= If w, = 0 then Pi has no effect on shape;
= If w, increases, so b and the curve is pulled toward P, and pushed away from PJ, for j not= i;
= If w, decreases, so b and the curve is pushed away from P, and pulled toward PJ, for j not= i ; and = If w, -> infinity then b -> 1 and B, -> Põ if u in [ti, trEk+i) The Problem Various techniques have been attempted for creating a NURBS surface in different fields of art. For example, International Publication No. WO
02/37422, and U.S. Patent No. 6,765,570, propose various techniques for the manual generation of editable (NURBS) surfaces used in analyzing and interpreting seismic events.
Other conventional applications propose creating an editable (NURBS) surface using interpolation techniques well known in the art. Such techniques may be referred to as an "exact-fit" approach to determining the editable surface. An exact-fit approach is more likely to render a less uniform, if not non-uniform, editable surface as compared to a "best-fit" approximation of the editable surface.
Moreover, the exact-fit approach to defining an editable NURBS surface may be impractical, if not cost prohibitive, for large 3D data sets often encountered in the analysis of seismic data.
Other conventional methods for converting unordered points to surfaces are generally described in Shepard, D., "A two dimensional interpolation function for irregular spaced data," 1968, pp. 517-524, Proceedings 23rd ACM National Conference.
One well-known method for interpolating scattered points, defines a surface value based on given points and weight functions. Although the result is relatively rough, the fundamental idea has inspired many other methods.
Another well-known method, generally referred to as the "Thin Plate Spline"
method has been incorporated in medical imaging applications and is capable of precisely guiding the target surface such that it passes through all given points. See Hardy, R. L. Desmarrais, R.N., "Interpolation using surface spline," 1972, pp.
189-197, Journal of Aircraft 9 and Dyn, N., "Interpolation in Scattered Data by Radial Functions,"
1987, pp. 47-61, In Chui, C.K.; Schumaker, L.L.; Ultreres, F.I. (ed) Topics in Multivariate Approximation. However, this technique requires the inversion of large matrices, which is computationally expensive and generally impractical for a large number of input points that may only require a best-fit approximation.
Another conventional method for converting a point cloud to a surface is described in Hoppe, H., et al., "Surface Reconstruction from Unorganized Points, 1992, pp. 71-78, Comput. Graph. 26. However, this method assumes that the input points are evenly distributed over the entire domain, which is also impractical for input points that are densely populated in some areas and non-existent in other areas.
A need therefore, exists for automatically approximating an editable surface from a 3D data set or 3D point set comprising a large volume of unordered and/or unstructured data points, which may be imaged in the form of an editable NURBS surface.
SUMMARY OF THE INVENTION
The present invention meets the above needs and overcomes one or more deficiencies in the prior by providing systems and methods for automatically approximating an editable surface from a 3D data set or a 3D point set that may comprise a large volume of unstructured and/or unordered data points.
In one embodiment, the present invention includes a system for approximating an editable surface from a 3D data set comprising a computer-readable memory medium configured to store a program of instructions capable of being executable to implement:
i) selecting a point set from the 3D data set; ii) determining a best-fit plane for the point set; iii) projecting at least a portion of the point set on to the best-fit plane; iv) determining a boundary for the best-fit plane according to the projected point set; v) projecting a grid on to the best-fit plane within the boundary, the grid containing a plurality of grid points; vi) initializing the grid points; vii) determining a value for at least a portion of the grid points; viii) selecting at least a portion of the grid points with a value; and ix) imaging the editable surface using at least a portion of the selected grid points.
In another embodiment, the present invention includes a system for approximating an editable surface from a 3D point set comprising a computer-readable memory medium configured to store a program of instructions capable of being executable to implement: i) determining a best-fit plane for the point set;
ii) projecting the point set on to the best-fit plane; iii) determining a boundary for the best-fit plane according to the point set; iv) projecting a grid on to the best-fit plane within the boundary, the grid containing a plurality of grid points; v) initializing the grid points to zero; vi) determining a value for the grid points; vii) refining the value for at least a portion of the grid points; viii) selecting at least a portion of the grid points based on one of the value and the refined value; and ix) imaging the editable surface using the selected grid points, the editable surface representing a best-fit approximation of the point set to the editable surface.
In another embodiment, the present invention includes a method for approximating an editable surface from a 3D data set comprising the steps of:
i) selecting a point set from 3D data set; ii) determining a best-fit plane for the point set; iii) projecting at least a portion of the point set on to the best-fit plane; iv) determining a boundary for the best-fit plane according to the projected point set; v) projecting a grid on to the best-fit plane within the boundary, the grid containing a plurality of grid points;
vi) initializing the grid points; vii) determining a value for at least a portion of the grid points; viii) selecting at least a portion of the grid points with a value;
and ix) imaging the editable surface using at least a portion of the selected grid points.
In yet another embodiment, the present invention includes a method for approximating an editable surface from a 3D point set comprising the steps of:
i) determining a best-fit plane for the point set; ii) projecting a point set on to the best-fit plane; iii) determining a boundary for the best-fit plane according to the point set; iv) projecting a_ grid on to the best-fit plane within the boundary, the grid containing a plurality of grid points; v) initializing the grid points to zero; vi) determining a value for the grid points; vii) refining the value for at least a portion of the grid points; viii) selecting at least a portion of the grid points based one of the value and the refined value;
and ix) imaging the editable surface using the selected grid points, the editable surface representing a best-fit approximation of the point set to the editable surface.
These and other objects, features and advantages of the present invention will become apparent to those skilled in the art from the following description of the various embodiments and related drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
The present invention is described below with reference to the accompanying drawings in which like elements are referenced with like reference numerals, and in which:
FIG. 1 is a block diagram illustrating one embodiment of a program for implementing the present invention.
FIG. 1A is a schematic diagram generally illustrating the flow of data into the approximation module illustrated in FIG. 1.
FIG. 2 is a flow diagram illustrating one embodiment of a method for implementing the present invention.
FIG. 3 illustrates step 212 in FIG. 2.
FIG. 4 illustrates step 216 and step 218 in FIG. 2.
ii) projecting the point set on to the best-fit plane; iii) determining a boundary for the best-fit plane according to the point set; iv) projecting a grid on to the best-fit plane within the boundary, the grid containing a plurality of grid points; v) initializing the grid points to zero; vi) determining a value for the grid points; vii) refining the value for at least a portion of the grid points; viii) selecting at least a portion of the grid points based on one of the value and the refined value; and ix) imaging the editable surface using the selected grid points, the editable surface representing a best-fit approximation of the point set to the editable surface.
In another embodiment, the present invention includes a method for approximating an editable surface from a 3D data set comprising the steps of:
i) selecting a point set from 3D data set; ii) determining a best-fit plane for the point set; iii) projecting at least a portion of the point set on to the best-fit plane; iv) determining a boundary for the best-fit plane according to the projected point set; v) projecting a grid on to the best-fit plane within the boundary, the grid containing a plurality of grid points;
vi) initializing the grid points; vii) determining a value for at least a portion of the grid points; viii) selecting at least a portion of the grid points with a value;
and ix) imaging the editable surface using at least a portion of the selected grid points.
In yet another embodiment, the present invention includes a method for approximating an editable surface from a 3D point set comprising the steps of:
i) determining a best-fit plane for the point set; ii) projecting a point set on to the best-fit plane; iii) determining a boundary for the best-fit plane according to the point set; iv) projecting a_ grid on to the best-fit plane within the boundary, the grid containing a plurality of grid points; v) initializing the grid points to zero; vi) determining a value for the grid points; vii) refining the value for at least a portion of the grid points; viii) selecting at least a portion of the grid points based one of the value and the refined value;
and ix) imaging the editable surface using the selected grid points, the editable surface representing a best-fit approximation of the point set to the editable surface.
These and other objects, features and advantages of the present invention will become apparent to those skilled in the art from the following description of the various embodiments and related drawings.
BRIEF DESCRIPTION OF THE DRAWINGS
The present invention is described below with reference to the accompanying drawings in which like elements are referenced with like reference numerals, and in which:
FIG. 1 is a block diagram illustrating one embodiment of a program for implementing the present invention.
FIG. 1A is a schematic diagram generally illustrating the flow of data into the approximation module illustrated in FIG. 1.
FIG. 2 is a flow diagram illustrating one embodiment of a method for implementing the present invention.
FIG. 3 illustrates step 212 in FIG. 2.
FIG. 4 illustrates step 216 and step 218 in FIG. 2.
FIG. 5 illustrates step 220 and step 222 in FIG. 2.
FIG. 5A illustrates step 224 in FIG. 2 and represents the area A in FIG. 5.
FIG. 6 illustrates step 228, step 230 and step 232 in FIG. 2.
FIG. 7 illustrates step 236 and step 238 in FIG. 2.
FIG. 8 is an image illustrating the input data points used to approximate the editable surface imaged at step 240 in FIG. 2.
FIG. 9 is another perspective of a portion of the image in FIG. 8.
FIG. 10 is another perspective of a portion of the image in FIG. 9.
DETAILED DESCRIPTION OF THE INVENTION
The subject matter of the present invention is described with specificity, however, the description itself is not intended to limit the scope of the invention.
The claimed subject matter thus, might also be embodied in other ways, to include different steps or combinations of steps similar to the ones described herein, in conjunction with other present or future technologies. Moreover, although the term "step" may be used herein to connote different elements of methods employed, the term should not be interpreted as implying any particular order among or between various steps herein disclosed unless and except when the order of individual steps is explicitly described.
The present invention provides an improved system and method for analyzing 3D
data sets and/or 3D point sets. The invention may be described in the general context of a computer-executable program of instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc., that perform particular tasks or implement particular abstract data types. Moreover, those skilled in the art will appreciate that the invention may be practiced with a variety of computer-system configurations, including hand-held devices, multiprocessor systems, microprocessor-based or programmable-consumer electronics, minicomputers, mainframe computers, and the like. Any number of computer-systems and computer networks are acceptable for use with the present invention. The invention may be practiced in distributed-computing environments where tasks are performed by remote-processing devices that are linked through a communications network. In a distributed-computing environment, program modules may be located in both local and remote computer-storage media including memory storage devices. The computer-useable instructions form an interface to allow a computer to react according to a source of input. The instructions cooperate with other code segments to initiate a variety of tasks in response to data received in conjunction with the source of the received data.
The present invention may therefore, be implemented using hardware, software or a combination thereof, in a computer system or other processing system.
FIG. 1 is a block diagram illustrating one embodiment of a software program 100 for implementing the present invention. At the base of program 100 is an operating system 102.
A
suitable operating system 102 may include, for example, a Windows operating system from Microsoft Corporation, or other operating systems as would be apparent to one of skill in the relevant art.
Menu and windowing software 104 overlays operating system 102. The menu and windowing software 104 are used to provide various menus and windows to facilitate interaction with the user, and to obtain user input and instructions. As would be readily apparent to one of skill in the relevant art, any number of menu and windowing software programs could be used in conjunction with the present invention.
A basic graphics library 106 overlays the menu and windowing software 104.
The basic graphics library 106 is an application programming interface (API) for computer graphics. The functions performed by the basic graphics library 106 may include, for example, geometric and raster primitives, RGBA or color index mode, display list or immediate mode, viewing and modeling transformations, lighting and shading, hidden surface removal, alpha blending (translucency), anti-aliasing, texture mapping, feedback and selection, stencil planes, and accumulation buffer.
A rendering application 108 overlays basic graphics library 106. As will be understood by those skilled in the art, the rendering application 108 may include a suite of tools for 2D/3D seismic data interpretations, including interactive horizon and fault management, 3D visualization, and attribute analysis. For example, Landmark Graphics Corporation's SeisVisionTM platform is a seismic rendering application appropriate for use with the present invention.
Overlaying the other elements of program 100 is an approximation module 110.
The approximation module 110 is configured to interact with 3D data sets representing predetermined objects such as, for example, horizons and faults or 3D point sets comprising arbitrary, and/or unstructured data points. In a manner generally well known in the art, the approximation module 110 interfaces with, and utilizes the functions carried out by, the rendering application 108, the basic graphics library 106, the menu and windowing software 104, and the operating system 102. The approximation module 110 may be written in an object oriented programming language such as, for example, C++ to allow the creation and use of objects and object functionality.
The program 100 illustrated in FIG. 1 may be executed or implemented through the use of a computer system incorporating the program 100 and various hardware components. The hardware components may include, for example, a processor, memory (e.g., random access memory and/or non-volatile memory devices), one or more input devices, one or more display devices, and one or more interface devices. These hardware components may be interconnected according to a variety of configurations.
Non-volatile memory devices may include, for example, devices such as tape drives, disk drives, semiconductor ROM or EEPROM. Input devices may include, for example, devices such as a keyboard, a mouse, a digitizing pad, a track ball, a touch-sensitive pad and/or a light pen. Display devices may include, for example, devices such as monitors, projectors and/or head-mounted displays. Interface devices may be configured to require digital image data from one or more acquisition devices and/or from one or more remote computers or storage devices through a network.
Any variety of acquisition devices may be used depending on the type of object being imaged. The acquisition device(s) may sense various forms of mechanical energy (e.g., acoustic energy, displacement and/or stress/strain) and/or electromagnetic energy (e.g., light energy, radio wave energy, current and/or voltage).
A processor may be configured to reprogram instructions and/or data from RAM
and/or non-volatile memory devices, and to store computational results into RAM and/or non-volatile memory devices. The program instructions direct the processor to operate on 3D data sets and/or 3D point sets based on the methods described herein.
The input data may be provided to the computer system through a variety of mechanisms.
For example, the input data may be acquired into non-volatile memory and/or RAM
using one or more interface devices. As another example, the input data may be supplied to the computer system through a memory medium such as a disk or a tape, which is loaded into/onto one of the non-volatile memory devices. In this case, the input data will have been previously recorded onto the memory medium.
It is noted that the input data may not necessarily be raw sensor data obtained by an acquisition device. For example, the input data may be the result of one or more processing operations using a set of raw sensor data. The processing operation(s) may be performed by the computer system and/or one or more other computers.
The method of the present invention may be realized in one or more software programs or modules, which are stored onto any of a variety of memory media such as CD-ROM, magnetic disk, bubble memory, semiconductor memory (e.g., any of a various types of RAM or ROM). Furthermore, the software program(s) and/or their results may be transmitted over a variety of carrier media such as optical fiber, metallic wire, free space and/or through any of a variety of networks such as the intemet.
In FIG. 1A, for example, a 3D data set 112 may comprise any predetermined object such as, for example, geoanomalies, geobodies, horizons, faults and/or other surfaces or any 3D point set comprising arbitrary and/or unstructured data points. The 3D data set 112 may be converted to a 3D point set 114, if necessary. The entire 3D
point set 114, or a portion thereof, may be selected from the converted 3D
data set 112 before being processed by the approximation module 110 to yield an image of an editable surface.
FIG. 5A illustrates step 224 in FIG. 2 and represents the area A in FIG. 5.
FIG. 6 illustrates step 228, step 230 and step 232 in FIG. 2.
FIG. 7 illustrates step 236 and step 238 in FIG. 2.
FIG. 8 is an image illustrating the input data points used to approximate the editable surface imaged at step 240 in FIG. 2.
FIG. 9 is another perspective of a portion of the image in FIG. 8.
FIG. 10 is another perspective of a portion of the image in FIG. 9.
DETAILED DESCRIPTION OF THE INVENTION
The subject matter of the present invention is described with specificity, however, the description itself is not intended to limit the scope of the invention.
The claimed subject matter thus, might also be embodied in other ways, to include different steps or combinations of steps similar to the ones described herein, in conjunction with other present or future technologies. Moreover, although the term "step" may be used herein to connote different elements of methods employed, the term should not be interpreted as implying any particular order among or between various steps herein disclosed unless and except when the order of individual steps is explicitly described.
The present invention provides an improved system and method for analyzing 3D
data sets and/or 3D point sets. The invention may be described in the general context of a computer-executable program of instructions, such as program modules, being executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc., that perform particular tasks or implement particular abstract data types. Moreover, those skilled in the art will appreciate that the invention may be practiced with a variety of computer-system configurations, including hand-held devices, multiprocessor systems, microprocessor-based or programmable-consumer electronics, minicomputers, mainframe computers, and the like. Any number of computer-systems and computer networks are acceptable for use with the present invention. The invention may be practiced in distributed-computing environments where tasks are performed by remote-processing devices that are linked through a communications network. In a distributed-computing environment, program modules may be located in both local and remote computer-storage media including memory storage devices. The computer-useable instructions form an interface to allow a computer to react according to a source of input. The instructions cooperate with other code segments to initiate a variety of tasks in response to data received in conjunction with the source of the received data.
The present invention may therefore, be implemented using hardware, software or a combination thereof, in a computer system or other processing system.
FIG. 1 is a block diagram illustrating one embodiment of a software program 100 for implementing the present invention. At the base of program 100 is an operating system 102.
A
suitable operating system 102 may include, for example, a Windows operating system from Microsoft Corporation, or other operating systems as would be apparent to one of skill in the relevant art.
Menu and windowing software 104 overlays operating system 102. The menu and windowing software 104 are used to provide various menus and windows to facilitate interaction with the user, and to obtain user input and instructions. As would be readily apparent to one of skill in the relevant art, any number of menu and windowing software programs could be used in conjunction with the present invention.
A basic graphics library 106 overlays the menu and windowing software 104.
The basic graphics library 106 is an application programming interface (API) for computer graphics. The functions performed by the basic graphics library 106 may include, for example, geometric and raster primitives, RGBA or color index mode, display list or immediate mode, viewing and modeling transformations, lighting and shading, hidden surface removal, alpha blending (translucency), anti-aliasing, texture mapping, feedback and selection, stencil planes, and accumulation buffer.
A rendering application 108 overlays basic graphics library 106. As will be understood by those skilled in the art, the rendering application 108 may include a suite of tools for 2D/3D seismic data interpretations, including interactive horizon and fault management, 3D visualization, and attribute analysis. For example, Landmark Graphics Corporation's SeisVisionTM platform is a seismic rendering application appropriate for use with the present invention.
Overlaying the other elements of program 100 is an approximation module 110.
The approximation module 110 is configured to interact with 3D data sets representing predetermined objects such as, for example, horizons and faults or 3D point sets comprising arbitrary, and/or unstructured data points. In a manner generally well known in the art, the approximation module 110 interfaces with, and utilizes the functions carried out by, the rendering application 108, the basic graphics library 106, the menu and windowing software 104, and the operating system 102. The approximation module 110 may be written in an object oriented programming language such as, for example, C++ to allow the creation and use of objects and object functionality.
The program 100 illustrated in FIG. 1 may be executed or implemented through the use of a computer system incorporating the program 100 and various hardware components. The hardware components may include, for example, a processor, memory (e.g., random access memory and/or non-volatile memory devices), one or more input devices, one or more display devices, and one or more interface devices. These hardware components may be interconnected according to a variety of configurations.
Non-volatile memory devices may include, for example, devices such as tape drives, disk drives, semiconductor ROM or EEPROM. Input devices may include, for example, devices such as a keyboard, a mouse, a digitizing pad, a track ball, a touch-sensitive pad and/or a light pen. Display devices may include, for example, devices such as monitors, projectors and/or head-mounted displays. Interface devices may be configured to require digital image data from one or more acquisition devices and/or from one or more remote computers or storage devices through a network.
Any variety of acquisition devices may be used depending on the type of object being imaged. The acquisition device(s) may sense various forms of mechanical energy (e.g., acoustic energy, displacement and/or stress/strain) and/or electromagnetic energy (e.g., light energy, radio wave energy, current and/or voltage).
A processor may be configured to reprogram instructions and/or data from RAM
and/or non-volatile memory devices, and to store computational results into RAM and/or non-volatile memory devices. The program instructions direct the processor to operate on 3D data sets and/or 3D point sets based on the methods described herein.
The input data may be provided to the computer system through a variety of mechanisms.
For example, the input data may be acquired into non-volatile memory and/or RAM
using one or more interface devices. As another example, the input data may be supplied to the computer system through a memory medium such as a disk or a tape, which is loaded into/onto one of the non-volatile memory devices. In this case, the input data will have been previously recorded onto the memory medium.
It is noted that the input data may not necessarily be raw sensor data obtained by an acquisition device. For example, the input data may be the result of one or more processing operations using a set of raw sensor data. The processing operation(s) may be performed by the computer system and/or one or more other computers.
The method of the present invention may be realized in one or more software programs or modules, which are stored onto any of a variety of memory media such as CD-ROM, magnetic disk, bubble memory, semiconductor memory (e.g., any of a various types of RAM or ROM). Furthermore, the software program(s) and/or their results may be transmitted over a variety of carrier media such as optical fiber, metallic wire, free space and/or through any of a variety of networks such as the intemet.
In FIG. 1A, for example, a 3D data set 112 may comprise any predetermined object such as, for example, geoanomalies, geobodies, horizons, faults and/or other surfaces or any 3D point set comprising arbitrary and/or unstructured data points. The 3D data set 112 may be converted to a 3D point set 114, if necessary. The entire 3D
point set 114, or a portion thereof, may be selected from the converted 3D
data set 112 before being processed by the approximation module 110 to yield an image of an editable surface.
Referring now to FIG. 2, a method 200 for approximating an image of an editable surface from a 3D data set or a 3D point set is illustrated.
In step 202, a 3D data set may be used to input predetermined objects or 3D
point sets comprising arbitrary and/or unstructured data points.
In step 204, an object may be selected from the 3D data set such as, for example, geoanomalies, geobodies, horizons, faults and/or other surfaces. Such objects may be predetermined using a number of imaging techniques well known in the art.
In step 206, the selected object is converted to a 3D point set using conversion techniques well known in the art. The converted 3D point set comprises arbitrary and/or unstructured data points corresponding with the initial data points comprising the initial input 3D data set. The 3D point set therefore, may comprise arbitrary and/or unstructured data points forming an arbitrary surface.
In step 208, a 3D point set may be selected from the initial input 3D data set (step 202) or the converted 3D point set (step 206), which may comprise all of the data points from the initial input 3D data set or a portion thereof. Thus, the selected 3D
point set comprises arbitrary and/or unstructured data points corresponding with at least a portion of the initial data points comprising the initial 3D data set.
In step 210, the 3D point set selected in step 208 is analyzed to determine if it contains a dominant (long) axis. If a dominant axis can be determined for the 3D point set, then the next step 212 is unnecessary. For certain predetermined objects selected in step 204, a dominant axis may be readily apparent to one skilled in the art.
In step 212, an ellipsoid 304 is fit to the 3D point set (point cloud) illustrated in FIG. 3 by programming means well known in the art. The 3D point set is represented by data points 302. The ellipsoid 304 includes three cardinal axes, which may be denoted by the long axis (L), the middle axis (M) and the short axis (S). In the case of a perfect circle, however, all three of the cardinal axes will have the same length. The short axis vector of the ellipsoid 304 is compared against the world coordinate axes (x, y, z) to determine the primary direction of a best-fit plane and therefore, which of the following three representations should be used for determining a best-fit plane in step 214:
z =f (x, y) y ¨f (z, x) x = f (y, z) In step 214, the best-fit plane is determined by using a least squares approximation to achieve a best-fit plane equation. For example, let T be the matrix of n data points comprising the 3D point set such that:
r ia13:fl where:
(00,60 (xi) Yif Z1)if z Y) Z dominant (121,15õ ,e) = (zI 7t) if y f(z,x)- Y
dominant 020, )6,4 = ()Plata) if x f(y, z) - X &mina E n]
Let M be the 4x4 matrix formed by T = T' where T' is the transpose of T. The least squares, best-fitting plane, equation is found by inverting the matrix M using LU
decomposition, which yields the coefficients A, B, C, D in the following equation that may be used to determine the best-fit plane for the data points 302 comprising the 3D
point set:
Ax + By + Cy +D = 0 In step 216, at least a portion of the data points 302 comprising the 3D point set are projected onto the best-fit plane 404 illustrated in FIG. 4 by programming means well known in the art. The more data points 302 that are projected onto the best-fit plane 404, the more accurate the results will be. Nevertheless, the number of data points projected onto the best-fit plane 404 may be determined as a matter of preference. The data points 302 projected onto the best-fit plane 404 are illustrated as projected data points 402. The projected data points 402 may be represented in a parametric form z' ¨
f(x',/) where x', y', z' denote the orthogonal coordinate system of the best-fit plane 404 wherein x' and y' are aligned with the best-fit plane 404, and z' is perpendicular to the best-fit plane 404.
In step 218, boundaries for the best-fit plane are determined according to data point extremes, meaning the outermost projected data points 402. In FIG. 4, the best-fit plane 404 is illustrated with boundaries that are evident from the outermost projected data points 402.
In step 220, a grid is projected onto the bounded best-fit plane 404 illustrated in FIG. 5 by programming means well known in the art. The grid projected onto the best-fit plane 404 may be represented by N x M where N and M are predetermined grid parameters. The predetermined grid parameters (N x M) may be the same (N=M) or different (NOM).
In step 202, a 3D data set may be used to input predetermined objects or 3D
point sets comprising arbitrary and/or unstructured data points.
In step 204, an object may be selected from the 3D data set such as, for example, geoanomalies, geobodies, horizons, faults and/or other surfaces. Such objects may be predetermined using a number of imaging techniques well known in the art.
In step 206, the selected object is converted to a 3D point set using conversion techniques well known in the art. The converted 3D point set comprises arbitrary and/or unstructured data points corresponding with the initial data points comprising the initial input 3D data set. The 3D point set therefore, may comprise arbitrary and/or unstructured data points forming an arbitrary surface.
In step 208, a 3D point set may be selected from the initial input 3D data set (step 202) or the converted 3D point set (step 206), which may comprise all of the data points from the initial input 3D data set or a portion thereof. Thus, the selected 3D
point set comprises arbitrary and/or unstructured data points corresponding with at least a portion of the initial data points comprising the initial 3D data set.
In step 210, the 3D point set selected in step 208 is analyzed to determine if it contains a dominant (long) axis. If a dominant axis can be determined for the 3D point set, then the next step 212 is unnecessary. For certain predetermined objects selected in step 204, a dominant axis may be readily apparent to one skilled in the art.
In step 212, an ellipsoid 304 is fit to the 3D point set (point cloud) illustrated in FIG. 3 by programming means well known in the art. The 3D point set is represented by data points 302. The ellipsoid 304 includes three cardinal axes, which may be denoted by the long axis (L), the middle axis (M) and the short axis (S). In the case of a perfect circle, however, all three of the cardinal axes will have the same length. The short axis vector of the ellipsoid 304 is compared against the world coordinate axes (x, y, z) to determine the primary direction of a best-fit plane and therefore, which of the following three representations should be used for determining a best-fit plane in step 214:
z =f (x, y) y ¨f (z, x) x = f (y, z) In step 214, the best-fit plane is determined by using a least squares approximation to achieve a best-fit plane equation. For example, let T be the matrix of n data points comprising the 3D point set such that:
r ia13:fl where:
(00,60 (xi) Yif Z1)if z Y) Z dominant (121,15õ ,e) = (zI 7t) if y f(z,x)- Y
dominant 020, )6,4 = ()Plata) if x f(y, z) - X &mina E n]
Let M be the 4x4 matrix formed by T = T' where T' is the transpose of T. The least squares, best-fitting plane, equation is found by inverting the matrix M using LU
decomposition, which yields the coefficients A, B, C, D in the following equation that may be used to determine the best-fit plane for the data points 302 comprising the 3D
point set:
Ax + By + Cy +D = 0 In step 216, at least a portion of the data points 302 comprising the 3D point set are projected onto the best-fit plane 404 illustrated in FIG. 4 by programming means well known in the art. The more data points 302 that are projected onto the best-fit plane 404, the more accurate the results will be. Nevertheless, the number of data points projected onto the best-fit plane 404 may be determined as a matter of preference. The data points 302 projected onto the best-fit plane 404 are illustrated as projected data points 402. The projected data points 402 may be represented in a parametric form z' ¨
f(x',/) where x', y', z' denote the orthogonal coordinate system of the best-fit plane 404 wherein x' and y' are aligned with the best-fit plane 404, and z' is perpendicular to the best-fit plane 404.
In step 218, boundaries for the best-fit plane are determined according to data point extremes, meaning the outermost projected data points 402. In FIG. 4, the best-fit plane 404 is illustrated with boundaries that are evident from the outermost projected data points 402.
In step 220, a grid is projected onto the bounded best-fit plane 404 illustrated in FIG. 5 by programming means well known in the art. The grid projected onto the best-fit plane 404 may be represented by N x M where N and M are predetermined grid parameters. The predetermined grid parameters (N x M) may be the same (N=M) or different (NOM).
In step 222, each grid point 502, defined by an intersection on the grid illustrated in FIG. 5, is initialized, preferably to zero. In this manner, a grid point value may be determined for each grid point 502 illustrated in FIG. 5.
In step 224, the value for each grid point 502 illustrated in FIG. 5 may be For each projected data point 402 within the predetermined area 5A, a data point I
1.0 di I
where Gi is the grid point value at each location of a grid point 502 based upon the data point values (Ei) for each projected data point 402 lying within the predetermined area 5A and its distance (di) to the grid point 502 being valued. This formula is valid for all distances (di) that do not equal zero. For those distances (di) that equal zero, the grid point value (GO is assigned the data point value (Ei) because it is known exactly. Thus, in FIG. 5A, the grid point value (GO for the grid point 502 is determined using each of the data point values (El, E2, E3, E4, E5) and their respective distance (d1, d2, d3, d4, d5) to the grid point 502 being valued. The same procedure may be applied for determining grid point values for the remaining grid points 502 by shifting the predetermined area 5A
such that each grid point 502 being valued is in the same position within the predetermined area 5A. For grid points 502 that are near or on the borders of the best-fit plane 404, the grid point values are determined in the same manner except that the predetermined area 5A may partially reside outside of the best-fit plane 404 in order to maintain the same grid point position within the predetermined area 5A.
Consequently, there are no projected data points 402 to scan outside of the best-fit plane 404 for consideration in determining the value of grid points 502 lying near or on the border of the best-fit plane 404.
In step 226, the grid point values may be refined in the manner described in reference to steps 228 through 236 once all of the grid point values (GO have been determined.
In step 228, the grid density is increased (e.g., doubled) and new grid point values (Gi') for the new grid points are linearly interpolated between the former grid points 502 using the former grid point values (Gi) and interpolation techniques well known in the art.
In step 224, the value for each grid point 502 illustrated in FIG. 5 may be For each projected data point 402 within the predetermined area 5A, a data point I
1.0 di I
where Gi is the grid point value at each location of a grid point 502 based upon the data point values (Ei) for each projected data point 402 lying within the predetermined area 5A and its distance (di) to the grid point 502 being valued. This formula is valid for all distances (di) that do not equal zero. For those distances (di) that equal zero, the grid point value (GO is assigned the data point value (Ei) because it is known exactly. Thus, in FIG. 5A, the grid point value (GO for the grid point 502 is determined using each of the data point values (El, E2, E3, E4, E5) and their respective distance (d1, d2, d3, d4, d5) to the grid point 502 being valued. The same procedure may be applied for determining grid point values for the remaining grid points 502 by shifting the predetermined area 5A
such that each grid point 502 being valued is in the same position within the predetermined area 5A. For grid points 502 that are near or on the borders of the best-fit plane 404, the grid point values are determined in the same manner except that the predetermined area 5A may partially reside outside of the best-fit plane 404 in order to maintain the same grid point position within the predetermined area 5A.
Consequently, there are no projected data points 402 to scan outside of the best-fit plane 404 for consideration in determining the value of grid points 502 lying near or on the border of the best-fit plane 404.
In step 226, the grid point values may be refined in the manner described in reference to steps 228 through 236 once all of the grid point values (GO have been determined.
In step 228, the grid density is increased (e.g., doubled) and new grid point values (Gi') for the new grid points are linearly interpolated between the former grid points 502 using the former grid point values (Gi) and interpolation techniques well known in the art.
In step 230, multiple new (best-fit) planes are determined using the equations described in reference to step 214 and a predetermined area, which may be the same predetermined area 5A illustrated in FIG. 5. In other words, at each projected data point 402, a new (best-fit) plane is determined (calculated) to fit the projected data points 402 and the new grid points within the surrounding predetermined area. In FIG. 6, for example, x' represents the incremental increase in grid density compared to previous grid spacing. A new (best-fit) plane 602 therefore, passes through the projected data point 402A and new grid point 606A within the surrounding predetermined area. A new (best-fit) plane 604 likewise, passes through projected data point 402B and new grid point 606B within the surrounding predetermined area. This step improves the accuracy of the interpolation performed in step 228, and determines a new (best-fit) plane for each projected data point 402.
In step 232, the interpolated new grid point values (Gi') may be reevaluated by replacing each interpolated new grid point value within the predetermined area with a reevaluated new grid point value based on the weighted average of the value for all new grid points that intersect a new (best-fit) plane that passes through the predetermined area and a projected data point within the predetermined area. In FIG. 6, for example, the reevaluated value for new grid point 612 may be determined by replacing its previously interpolated new grid point value with the weighted average of the value for new grid points 606A, 606B that intersect the new (best-fit) planes 602, 604, respectively, that pass through the predetermined area and contain the projected data points 402A, 402B, respectively.
In step 232, the interpolated new grid point values (Gi') may be reevaluated by replacing each interpolated new grid point value within the predetermined area with a reevaluated new grid point value based on the weighted average of the value for all new grid points that intersect a new (best-fit) plane that passes through the predetermined area and a projected data point within the predetermined area. In FIG. 6, for example, the reevaluated value for new grid point 612 may be determined by replacing its previously interpolated new grid point value with the weighted average of the value for new grid points 606A, 606B that intersect the new (best-fit) planes 602, 604, respectively, that pass through the predetermined area and contain the projected data points 402A, 402B, respectively.
In step 234, a Gaussian kernel may be applied to each of the reevaluated new grid point values in order to smooth the resulting reevaluated new grid point values so that they have a more natural continuity. The Gaussian smoothing of the reevaluated new grid point values may be performed in a manner well known in the art. Various modifications to the Gaussian process of smoothing may be apparent to those skilled in the art and applied as necessary to improve the smoothing process described herein.
In step 236, at least a portion of the reevaluated new grid point values may be further refined by returning to step 222 and repeating some or all of the steps described thus far. Upon reiterating step 224, however, the predetermined area 5A is incrementally reduced by a factor of one such that the surrounding predetermined area to be considered may be represented as (4-1) x (4-1). If, however, the reevaluated new grid point values are not refined, then the process proceeds to step 238. Likewise, if the grid point values are not refined according to step 226, the process may proceed directly to step 238.
In step 238, the grid point values and/or reevaluated new grid point values may be sampled for control points. The sampling process simply selects the last iteration of the process described in reference to step 236 such that the last value determined may be selected as a control point. In FIG. 7, for example, the best-fit plane 700 may include multiple values for grid point 702 that may be represented as Gi, G12 Gin.
Similarly, grid points 704, 706, and 708, for example, may have multiple values.
Accordingly, the last determined value for each grid point may be sampled (selected) in whole, or in part, to produce a set of control points for approximating an editable (NURBS) surface at any arbitrary resolution.
In step 236, at least a portion of the reevaluated new grid point values may be further refined by returning to step 222 and repeating some or all of the steps described thus far. Upon reiterating step 224, however, the predetermined area 5A is incrementally reduced by a factor of one such that the surrounding predetermined area to be considered may be represented as (4-1) x (4-1). If, however, the reevaluated new grid point values are not refined, then the process proceeds to step 238. Likewise, if the grid point values are not refined according to step 226, the process may proceed directly to step 238.
In step 238, the grid point values and/or reevaluated new grid point values may be sampled for control points. The sampling process simply selects the last iteration of the process described in reference to step 236 such that the last value determined may be selected as a control point. In FIG. 7, for example, the best-fit plane 700 may include multiple values for grid point 702 that may be represented as Gi, G12 Gin.
Similarly, grid points 704, 706, and 708, for example, may have multiple values.
Accordingly, the last determined value for each grid point may be sampled (selected) in whole, or in part, to produce a set of control points for approximating an editable (NURBS) surface at any arbitrary resolution.
In step 240, the editable surface may be imaged (rendered) using at least a portion of the selected control points and the rendering techniques generally known for creating a NURBS surface. The more control points that are used to image the editable surface, the more accurate the results will be. In other words, a more accurate best-fit In FIG. 8, for example, the input data points 802 and editable surface 804 are simultaneously imaged to illustrate the spatial relationship between the data points 802 In FIG. 9, the image in FIG. 8 has been edited to illustrate a portion of the data points 802 and the corresponding editable surface 804, which may be used for interpreting seismic events.
In FIG. 10, the data points 802 have been removed to reveal only the editable The present invention therefore, may be applied to the geophysical analysis of seismic data. Moreover, the present invention may be integrated with one or more system components described in reference to International Publication No. WO
In FIG. 10, the data points 802 have been removed to reveal only the editable The present invention therefore, may be applied to the geophysical analysis of seismic data. Moreover, the present invention may be integrated with one or more system components described in reference to International Publication No. WO
and/or U.S. Patent No. 6,765,570 for creating an editable surface from a plurality of seismic data points that may include, for example, x, y, z coordinates and a data value representing a predetermined object such as, for example, a horizon or a fault. The seismic image automatically rendered as a result of the present invention represents an The present invention, however, may also be applied to other types of 3D data sets such, for example, medical data and engineering data. It is therefore, contemplated
Claims (92)
1. A method for approximating an editable surface from a 3D data set comprising the steps of:
selecting a 3D point set from the 3D data set;
determining a best-fit plane for the 3D point set;
projecting at least a portion of the 3D point set onto the best-fit plane;
determining a boundary for the best-fit plane according to the projected 3D
point set;
projecting a grid onto the best-fit plane within the boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points;
determining a value for at least a portion of the grid points, wherein the grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point;
selecting at least a portion of the grid points with a value; and imaging the editable surface using at least a portion of the selected grid points.
selecting a 3D point set from the 3D data set;
determining a best-fit plane for the 3D point set;
projecting at least a portion of the 3D point set onto the best-fit plane;
determining a boundary for the best-fit plane according to the projected 3D
point set;
projecting a grid onto the best-fit plane within the boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points;
determining a value for at least a portion of the grid points, wherein the grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point;
selecting at least a portion of the grid points with a value; and imaging the editable surface using at least a portion of the selected grid points.
2. The method of claim 1, wherein the 3D point set comprises an arbitrary surface.
3. The method of claim 2, further comprising the steps of:
fitting an ellipsoid to the 3D point set; and determining a dominant axis for the 3D point set based on the ellipsoid.
fitting an ellipsoid to the 3D point set; and determining a dominant axis for the 3D point set based on the ellipsoid.
4. The method of claim 3, wherein the dominant axis is used to determine the best-fit plane.
5. The method of claim 1, wherein the 3D point set comprises one of the 3D
data set and a portion of the 3D data set.
data set and a portion of the 3D data set.
6. The method of claim 1, wherein the best-fit plane is determined by a dominant axis for the 3D point set.
7. The method of claim 1, wherein the boundary of the best-fit plane is determined by a plurality of point extremes contained in the projected 3D point set.
8. The method of claim 1, wherein the projected grid is based on predetermined grid parameters.
9. The method of claim 1, wherein the grid points are initialized to zero.
10. The method of claim 1, wherein the editable surface represents a best-fit approximation of the 3D point set to the editable surface.
11. The method of claim 1, wherein the editable surface represents a NURBS
surface.
surface.
12. The method of claim 1, further comprising the steps of:
refining the value for at least a portion of the at least a portion of the grid points;
and selecting at least a portion of the grid points based on one of the value and the refined value.
refining the value for at least a portion of the at least a portion of the grid points;
and selecting at least a portion of the grid points based on one of the value and the refined value.
13. The method of claim 12, further comprising the step of replacing the value with the refined value for the grid points selected with a value and a refined value.
14. The method of claim 1, further comprising the step of converting the 3D
data set to the 3D point set.
data set to the 3D point set.
15. A method for approximating an editable surface from a 3D point set comprising the steps of:
determining a best-fit plane for the 3D point set;
projecting the 3D point set onto the best-fit plane;
determining a boundary for the best-fit plane according to the 3D point set;
projecting a grid onto the best-fit plane within the boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points to zero;
determining a value for the grid points, wherein the grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point;
refining the value for at least a portion of the grid points;
selecting at least a portion of the grid points based on one of the value and the refined value; and imaging the editable surface using the selected grid points, the editable surface representing a best-fit approximation of the 3D point set to the editable surface.
determining a best-fit plane for the 3D point set;
projecting the 3D point set onto the best-fit plane;
determining a boundary for the best-fit plane according to the 3D point set;
projecting a grid onto the best-fit plane within the boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points to zero;
determining a value for the grid points, wherein the grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point;
refining the value for at least a portion of the grid points;
selecting at least a portion of the grid points based on one of the value and the refined value; and imaging the editable surface using the selected grid points, the editable surface representing a best-fit approximation of the 3D point set to the editable surface.
16. The method of claim 15, wherein the best-fit plane is determined by a dominant axis for the 3D point set.
17. The method of claim 15, wherein the boundary of the best-fit plane is determined by a plurality of point extremes contained in the projected 3D point set.
18. The method of claim 15, wherein the projected grid is based on predetermined grid parameters.
19. The method of claim 15, wherein the editable surface represents a NURBS
surface.
surface.
20. The method of claim 15, further comprising the step of replacing the value with the refined value for the grid points selected with a value and a refined value.
21. The method of claim 15, wherein refining the value comprises:
forming a new grid by increasing the grid density, the new grid containing a plurality of new grid points;
determining a value for the new grid points;
fitting a new plane to the new grid points;
reevaluating the value for the new grid points; and smoothing the value for the new grid points.
forming a new grid by increasing the grid density, the new grid containing a plurality of new grid points;
determining a value for the new grid points;
fitting a new plane to the new grid points;
reevaluating the value for the new grid points; and smoothing the value for the new grid points.
22. The method of claim 21, further comprising the steps of:
initializing the new grid points to zero; and determining a new value for the new grid points.
initializing the new grid points to zero; and determining a new value for the new grid points.
23. The method of claim 22, wherein the new grid point value is determined by one or more points in the 3D point set within another predetermined area of the new grid point and a distance from each point in the another predetermined area to the new grid point, wherein the another predetermined area is smaller than the predetermined area.
24. The method of claim 22, further comprising the step of repeating the steps in claim 21.
25. The method of claim 21, wherein the step of reevaluating the value for the new grid points comprises replacing each new grid point value within the predetermined area with a reevaluating new grid point value.
26. The method of claim 21, further comprising the step of refining the new grid point values.
27. The method of claim 26, further comprising the step of repeating one or more of the steps in claim 21.
28. A system for approximating an editable surface from a 3D data set comprising a computer-readable memory medium configured to store a program of instructions, the program instructions being executable to implement:
selecting a 3D point set from the 3D data set;
determining a best-fit plane for the 3D point set;
projecting at least a portion of the 3D point set onto the best-fit plane;
determining a boundary for the best-fit plane according to the projected 3D
point set;
projecting a grid onto the best-fit plane within the boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points;
determining a value for at least a portion of the grid points, wherein the grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point;
selecting at least a portion of the grid points with a value; and imaging the editable surface using at least a portion of the selected grid points.
selecting a 3D point set from the 3D data set;
determining a best-fit plane for the 3D point set;
projecting at least a portion of the 3D point set onto the best-fit plane;
determining a boundary for the best-fit plane according to the projected 3D
point set;
projecting a grid onto the best-fit plane within the boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points;
determining a value for at least a portion of the grid points, wherein the grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point;
selecting at least a portion of the grid points with a value; and imaging the editable surface using at least a portion of the selected grid points.
29. The system of claim 28, wherein the 3D point set comprises an arbitrary surface.
30. The system of claim 29, further comprising the steps of:
fitting an ellipsoid to the 3D point set; and determining a dominant axis for the 3D point set based on the ellipsoid.
fitting an ellipsoid to the 3D point set; and determining a dominant axis for the 3D point set based on the ellipsoid.
31. The system of claim 30, wherein the dominant axis is used to determine the best-fit plane.
32. The system of claim 28, wherein the 3D point set comprises one of the 3D data set and a portion of the 3D data set.
33. The system of claim 28, wherein the best-fit plane is determined by a dominant axis for the 3D point set.
34. The system of claim 28, wherein the boundary of the best-fit plane is determined by a plurality of point extremes contained in the projected 3D point set.
35. The system of claim 28, wherein the projected grid is based on predetermined grid parameters.
36. The system of claim 28, wherein the grid points are initialized to zero.
37. The system of claim 28, wherein the editable surface represents a best-fit approximation of the 3D point set to the editable surface.
38. The system of claim 28, wherein the editable surface represents a NURBS
surface.
surface.
39. The system of claim 28, further comprising the steps of:
refining the value for at least a portion of the at least a portion of the grid points;
and selecting at least a portion of the grid points based on one of the value and the refined value.
refining the value for at least a portion of the at least a portion of the grid points;
and selecting at least a portion of the grid points based on one of the value and the refined value.
40. The system of claim 39, further comprising the step of replacing the value with the refined value for the grid points selected with a value and a refined value.
41. The system of claim 28, further comprising the step of converting the 3D data set to the 3D point set.
42. A system for approximating an editable surface from a 3D point set comprising a computer-readable memory medium configured to store a program of instructions, the program instructions being executable to implement:
determining a best-fit plane for the 3D point set;
projecting the 3D point set onto the best-fit plane;
determining a boundary for the best-fit plane according to the 3D point set;
projecting a grid onto the best-fit plane within the boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points to zero;
determining a value for the grid points, wherein the grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point;
refining the value for at least a portion of the grid points;
selecting at least a portion of the grid points based on one of the value and the refined value; and imaging the editable surface using the selected grid points, the editable surface representing a best-fit approximation of the 3D point set to the editable surface.
determining a best-fit plane for the 3D point set;
projecting the 3D point set onto the best-fit plane;
determining a boundary for the best-fit plane according to the 3D point set;
projecting a grid onto the best-fit plane within the boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points to zero;
determining a value for the grid points, wherein the grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point;
refining the value for at least a portion of the grid points;
selecting at least a portion of the grid points based on one of the value and the refined value; and imaging the editable surface using the selected grid points, the editable surface representing a best-fit approximation of the 3D point set to the editable surface.
43. The system of claim 42, wherein the best-fit plane is determined by a dominant axis for the 3D point set.
44. The system of claim 42, wherein the boundary of the best-fit plane is determined by a plurality of point extremes contained in the projected 3D point set.
45. The system of claim 42, wherein the projected grid is based on predetermined grid parameters.
46. The system of claim 42, wherein the editable surface represents a NURBS
surface.
surface.
47. The system of claim 42, further comprising the step of replacing the value with the refined value for the grid points selected with a value and a refined value.
48. The system of claim 42, wherein refining the value comprises:
forming a new grid by increasing the grid density, the new grid containing a plurality of new grid points;
determining a value for the new grid points;
fitting a new plane to the new grid points;
reevaluating the value for the new grid points; and smoothing the value for the new grid points.
forming a new grid by increasing the grid density, the new grid containing a plurality of new grid points;
determining a value for the new grid points;
fitting a new plane to the new grid points;
reevaluating the value for the new grid points; and smoothing the value for the new grid points.
49. The system of claim 48, further comprising the steps of:
initializing the new grid points to zero; and determining a new value for the new grid points.
initializing the new grid points to zero; and determining a new value for the new grid points.
50. The system of claim 49, wherein the new grid point value is determined by one or more points in the 3D point set within another predetermined area of the new grid point and a distance from each point in the another predetermined area to the new grid point, wherein the another predetermined area is smaller than the predetermined area.
51. The system of claim 49, further comprising the step of repeating the steps in claim 48.
52. The system of claim 48, wherein the step of reevaluating the value for the new grid points comprises replacing each new grid point value within the predetermined area with a reevaluated new grid point value.
53. The system of claim 48, further comprising the step of refining the new grid point values.
54. The system of claim 53, comprising the step of repeating one or more of the steps in claim 48.
55. A computer-implemented method for approximating an editable surface from a projection of a 3D point set onto a best-fit plane for the 3D point set, which comprises:
projecting a grid onto the best-fit plane within a predetermined boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points;
selecting at least a portion of the grid points with a predetermined value, the selected grid points defining control points, wherein the predetermined grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point; and rendering the editable surface using at least a portion of the control points.
projecting a grid onto the best-fit plane within a predetermined boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points;
selecting at least a portion of the grid points with a predetermined value, the selected grid points defining control points, wherein the predetermined grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point; and rendering the editable surface using at least a portion of the control points.
56. The method of claim 55,wherein the 3D point set comprises an arbitrary surface.
57. The method of claim 56, further comprising the steps of:
fitting an ellipsoid to the 3D point set; and determining a dominant axis for the 3D point set based on the ellipsoid.
fitting an ellipsoid to the 3D point set; and determining a dominant axis for the 3D point set based on the ellipsoid.
58. The method of claim 57, wherein the dominant axis is used to determine the best-fit plane.
59. The method of claim 55, wherein the best-fit plane is determined by a dominant axis for the 3D point set.
60. The method of claim 55, wherein the predetermined boundary on the best-fit plane is determined by a plurality of point extremes contained in the projected 3D point set.
61. The method of claim 55; wherein the projected grid is based on predetermined grid parameters.
62. The method of claim 55, wherein the grid points are initialized to zero.
63. The method of claim 55, wherein the editable surface represents a best-fit approximation of the 3D point set to the editable surface.
64. The method of claim 55, wherein the editable surface represents a NURBS
surface.
surface.
65. The method of claim 55, further comprising the steps of:
refining the predetermined value for at least a portion of the grid points with the predetermined value; and selecting at least a portion of the grid points based on one of the predetermined value and the refined value.
refining the predetermined value for at least a portion of the grid points with the predetermined value; and selecting at least a portion of the grid points based on one of the predetermined value and the refined value.
66. The method of claim 65, further comprising the step of replacing the predetermined value with the refined value for the grid points selected with the predetermined value and the refined value.
67. The method of claim 65, wherein refining the predetermined value comprises:
forming a new grid by increasing the grid density, the new grid containing a plurality of new grid points;
determining a value for the new grid points;
fitting a new plane to the new grid points;
reevaluating the value for the new grid points; and smoothing the value for the new grid points.
forming a new grid by increasing the grid density, the new grid containing a plurality of new grid points;
determining a value for the new grid points;
fitting a new plane to the new grid points;
reevaluating the value for the new grid points; and smoothing the value for the new grid points.
68. The method of claim 65, further comprising the steps of:
initializing the new grid points to zero; and determining a new value for the new grid points.
initializing the new grid points to zero; and determining a new value for the new grid points.
69. The method of claim 68, wherein the new grid point value is determined by one or more points in the projected 3D point set within another predetermined area of the new grid point and a distance from each point in the another predetermined area to the new grid point, wherein the another predetermined area is smaller than the predetermined area.
70. The method of claim 68, further comprising the step of repeating forming, determining a value, fitting, reevaluating and smoothing.
71. The method of claim 67, wherein the step of reevaluating the value for the new grid points comprises replacing each new grid point value within the predetermined area with a reevaluated new grid point value.
72. The method of claim 67, further comprising the step of refining the new grid point values.
73. The method of claim 72, further comprising the step of repeating one or more of forming, determining a value, fitting, reevaluating and smoothing.
74. A computer-readable medium having recorded thereon computer executable instructions for execution by a computer for approximating an editable surface from a projection of a 3D point set onto a best-fit plane for the 3D point set, the instructions being executable to implement:
projecting a grid onto the best-fit plane within a predetermined boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points;
selecting at least a portion of the grid points with a predetermined value, the selected grid points defining control points, wherein the predetermined grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point; and rendering the editable surface using at least a portion of the control points.
projecting a grid onto the best-fit plane within a predetermined boundary, the grid containing a plurality of grid points, wherein each of the plurality of grid points is defined by an intersection on the grid;
initializing the grid points;
selecting at least a portion of the grid points with a predetermined value, the selected grid points defining control points, wherein the predetermined grid point value is determined by one or more points in the projected 3D point set within a predetermined area of the grid point and a distance from each point in the predetermined area to the grid point; and rendering the editable surface using at least a portion of the control points.
75. The computer-readable medium of claim 74, wherein the 3D point set comprises an arbitrary surface.
76. The computer-readable medium of claim 75, further comprising instructions executable to implement:
fitting an ellipsoid to the 3D point set; and determining a dominant axis for the 3D point set based on the ellipsoid.
fitting an ellipsoid to the 3D point set; and determining a dominant axis for the 3D point set based on the ellipsoid.
77. The computer-readable medium of claim 76, wherein the dominant axis is used to determine the best-fit plane.
78. The computer-readable medium of claim 74, wherein the best-fit plane is determined by a dominant axis for the 3D point set.
79. The computer-readable medium of claim 74, wherein the predetermined boundary on the best-fit plane is determined by a plurality of point extremes contained in the projected 3D point set.
80. The computer-readable medium of claim 74, wherein the projected grid is based on predetermined grid parameters.
81. The computer-readable medium of claim 74, wherein the grid points are initialized to zero.
82. The computer-readable medium of claim 74, wherein the editable surface represents a best-fit approximation of the 3D point set to the editable surface.
83. The computer-readable medium of claim 74, wherein the editable surface represents a NURBS surface.
84. The computer-readable medium of claim 74, further comprising instructions executable to implement:
refining the predetermined value for at least a portion of the grid points with the predetermined value; and selecting at least a portion of the grid points based on one of the 30 predetermined value and the refined value.
refining the predetermined value for at least a portion of the grid points with the predetermined value; and selecting at least a portion of the grid points based on one of the 30 predetermined value and the refined value.
85. The computer-readable medium of claim 84, further comprising instructions executable to implement replacing the predetermined value with the refined value for the grid points selected with the predetermined value and the refined value.
86. The computer-readable medium of claim 84, wherein refining the predetermined value comprises:
forming a new grid by increasing the grid density, the new grid containing a plurality of new grid points;
determining a value for the new grid points;
fitting a new plane to the new grid points;
reevaluating the value for the new grid points; and smoothing the value for the new grid points.
forming a new grid by increasing the grid density, the new grid containing a plurality of new grid points;
determining a value for the new grid points;
fitting a new plane to the new grid points;
reevaluating the value for the new grid points; and smoothing the value for the new grid points.
87. The computer-readable medium of claim 86, further comprising instructions executable to implement:
initializing the new grid points to zero; and determining a new value for the new grid points.
initializing the new grid points to zero; and determining a new value for the new grid points.
88. The computer-readable medium of claim 87, wherein the new grid point value is determined by one or more points in the projected 3D point set within another predetermined area of the new grid point and a distance from each point in the another predetermined area to the new grid point, wherein the another predetermined area is smaller than the predetermined area.
89. The computer-readable medium of claim 87, further comprising instructions executable to implement repeating the instructions for forming, determining a value, fitting, reevaluating and smoothing.
90. The computer-readable medium of claim 86, wherein the instruction for reevaluating the value for the new grid points comprises replacing each new grid point value within the predetermined area with a reevaluated new grid point value.
91. The computer-readable medium of claim 86, further comprising instructions executable to implement refining the new grid point values.
92. The computer-readable medium of claim 91, further comprising instructions executable to implement repeating one or more of the instructions for forming, determining a value, fitting, reevaluating and smoothing.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US56657404P | 2004-04-29 | 2004-04-29 | |
| US60/566,574 | 2004-04-29 | ||
| PCT/US2005/014937 WO2005108923A2 (en) | 2004-04-29 | 2005-04-29 | System and method for approximating an editable surface |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2564156A1 CA2564156A1 (en) | 2005-11-17 |
| CA2564156C true CA2564156C (en) | 2014-04-08 |
Family
ID=35320836
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA2564156A Active CA2564156C (en) | 2004-04-29 | 2005-04-29 | System and method for approximating an editable surface |
Country Status (10)
| Country | Link |
|---|---|
| US (2) | US7352369B2 (en) |
| EP (1) | EP1756517B1 (en) |
| CN (1) | CN100576250C (en) |
| AU (1) | AU2005241463C1 (en) |
| BR (1) | BRPI0509830A (en) |
| CA (1) | CA2564156C (en) |
| EA (1) | EA012579B1 (en) |
| MX (1) | MXPA06012485A (en) |
| NO (1) | NO339362B1 (en) |
| WO (1) | WO2005108923A2 (en) |
Families Citing this family (24)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101192307B (en) * | 2006-11-17 | 2012-05-23 | 鸿富锦精密工业（深圳）有限公司 | Point cloud triangular topological relations construction method |
| US7904092B2 (en) * | 2007-01-04 | 2011-03-08 | Cisco Technology, Inc. | Locally adjusted radio frequency coverage maps in wireless networks |
| US8284193B2 (en) * | 2007-06-28 | 2012-10-09 | Siemens Product Lifecycle Management Software Inc. | System and method for composition of NURBS surfaces |
| US20090027380A1 (en) * | 2007-07-23 | 2009-01-29 | Vivek Rajan | 3-D visualization |
| CN101377851A (en) * | 2007-08-29 | 2009-03-04 | 鸿富锦精密工业（深圳）有限公司 | System and method for computing a minimum distance from point cloud to point cloud |
| CN101441781B (en) * | 2007-11-23 | 2011-02-02 | 鸿富锦精密工业（深圳）有限公司 | Curved surface overturning method |
| US8253726B1 (en) | 2008-01-09 | 2012-08-28 | Spaceclaim Corporation, Inc. | Systems and methods for modifying three dimensional geometry using an arbitrary cross-section plane |
| FR2929417B1 (en) * | 2008-03-27 | 2010-05-21 | Univ Paris 13 | METHOD FOR DETERMINING A THREE-DIMENSIONAL REPRESENTATION OF AN OBJECT FROM POINTS, COMPUTER PROGRAM AND CORRESPONDING IMAGING SYSTEM |
| US8525871B2 (en) * | 2008-08-08 | 2013-09-03 | Adobe Systems Incorporated | Content-aware wide-angle images |
| US8948476B2 (en) | 2010-12-20 | 2015-02-03 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Determination of cardiac geometry responsive to doppler based imaging of blood flow characteristics |
| US20100168557A1 (en) * | 2008-12-30 | 2010-07-01 | Deno D Curtis | Multi-electrode ablation sensing catheter and system |
| US8900150B2 (en) | 2008-12-30 | 2014-12-02 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Intracardiac imaging system utilizing a multipurpose catheter |
| US9610118B2 (en) * | 2008-12-31 | 2017-04-04 | St. Jude Medical, Atrial Fibrillation Division, Inc. | Method and apparatus for the cancellation of motion artifacts in medical interventional navigation |
| US9186088B2 (en) * | 2011-09-29 | 2015-11-17 | Siemens Aktiengesellschaft | Active catheter reconstruction for interventional magnetic resonance imaging |
| US20130106887A1 (en) * | 2011-10-31 | 2013-05-02 | Christopher Tremblay | Texture generation using a transformation matrix |
| GB2505936A (en) * | 2012-09-17 | 2014-03-19 | Materialise Dental Nv | 3D modelling of scanned body |
| BR112015028562A2 (en) * | 2013-05-15 | 2017-07-25 | Logined Bv | method, system, and one or more computer readable storage media |
| US9412040B2 (en) * | 2013-12-04 | 2016-08-09 | Mitsubishi Electric Research Laboratories, Inc. | Method for extracting planes from 3D point cloud sensor data |
| CN105574935A (en) * | 2014-10-10 | 2016-05-11 | 联想(北京)有限公司 | Information processing method and device |
| WO2017007451A1 (en) * | 2015-07-07 | 2017-01-12 | Landmark Graphics Corporation | Tri-cubic and hybrid interpolation in a 3d texture shader |
| FR3043227A1 (en) * | 2015-11-04 | 2017-05-05 | Services Petroliers Schlumberger | |
| RU2617144C1 (en) * | 2015-12-09 | 2017-04-21 | Акционерное общество "Муромский завод радиоизмерительных приборов" | Method for simulating object movement trajectory |
| US11249208B2 (en) | 2016-10-14 | 2022-02-15 | Schlumberger Technology Corporation | Geologic structural model generation |
| CN110163863B (en) * | 2018-11-06 | 2022-11-04 | 腾讯科技（深圳）有限公司 | Three-dimensional object segmentation method, apparatus, and medium |
Family Cites Families (30)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5321612A (en) * | 1991-02-26 | 1994-06-14 | Swift Energy Company | Method for exploring for hydrocarbons utilizing three dimensional modeling of thermal anomalies |
| US5233568A (en) * | 1991-06-28 | 1993-08-03 | Atlantic Richfield Company | Geopressure analysis system |
| US5988862A (en) * | 1996-04-24 | 1999-11-23 | Cyra Technologies, Inc. | Integrated system for quickly and accurately imaging and modeling three dimensional objects |
| US5903458A (en) * | 1997-06-06 | 1999-05-11 | Ford Global Technologies, Inc. | System and method for forming geometric features using global reparametrization |
| US6037949A (en) * | 1997-08-04 | 2000-03-14 | Pixar Animation Studios | Texture mapping and other uses of scalar fields on subdivision surfaces in computer graphics and animation |
| US6300960B1 (en) * | 1997-08-04 | 2001-10-09 | Pixar Animation Studios | Realistic surface simulation in computer animation |
| US5951475A (en) * | 1997-09-25 | 1999-09-14 | International Business Machines Corporation | Methods and apparatus for registering CT-scan data to multiple fluoroscopic images |
| US6765570B1 (en) * | 1998-07-21 | 2004-07-20 | Magic Earth, Inc. | System and method for analyzing and imaging three-dimensional volume data sets using a three-dimensional sampling probe |
| JP4473358B2 (en) * | 1999-01-21 | 2010-06-02 | 株式会社東芝 | Diagnostic equipment |
| US7333648B2 (en) * | 1999-11-19 | 2008-02-19 | General Electric Company | Feature quantification from multidimensional image data |
| AU2001249765A1 (en) * | 2000-03-30 | 2001-10-15 | Align Technology, Inc. | System and method for separating three-dimensional models |
| US7006085B1 (en) * | 2000-10-30 | 2006-02-28 | Magic Earth, Inc. | System and method for analyzing and imaging three-dimensional volume data sets |
| MXPA03004228A (en) * | 2000-10-30 | 2004-12-02 | Magic Earth Inc | System and method for analyzing and imaging three-dimensional volume data sets. |
| US20020164067A1 (en) * | 2001-05-02 | 2002-11-07 | Synapix | Nearest neighbor edge selection from feature tracking |
| US7023432B2 (en) * | 2001-09-24 | 2006-04-04 | Geomagic, Inc. | Methods, apparatus and computer program products that reconstruct surfaces from data point sets |
| US6931367B2 (en) * | 2001-11-29 | 2005-08-16 | Faurecia Exhaust Systems, Inc. | Optimal rib design method for exhaust components |
| US7248259B2 (en) * | 2001-12-12 | 2007-07-24 | Technoguide As | Three dimensional geological model construction |
| US7062419B2 (en) * | 2001-12-21 | 2006-06-13 | Intel Corporation | Surface light field decomposition using non-negative factorization |
| US20050168460A1 (en) * | 2002-04-04 | 2005-08-04 | Anshuman Razdan | Three-dimensional digital library system |
| US6963671B2 (en) * | 2002-04-17 | 2005-11-08 | Mitsubishi Electric Research Labs, Inc. | Method for determining distances to a surface from a range image |
| AU2003245483A1 (en) * | 2002-06-12 | 2003-12-31 | Spatial Integrated Systems, Inc. | Discrete linear space sampling method and apparatus for generating digital 3d models |
| US7399220B2 (en) * | 2002-08-02 | 2008-07-15 | Kriesel Marshall S | Apparatus and methods for the volumetric and dimensional measurement of livestock |
| JP4185052B2 (en) * | 2002-10-15 | 2008-11-19 | ユニバーシティ オブ サザン カリフォルニア | Enhanced virtual environment |
| US20040075656A1 (en) * | 2002-10-18 | 2004-04-22 | Kimia Benjamin B. | Method and apparatus for multi-dimensional shape representation via shock flows |
| WO2004038654A2 (en) * | 2002-10-18 | 2004-05-06 | Exxonmobil Upstream Research Company | A method for rapid fault interpretation of fault surfaces generated to fit three-dimensional seismic discontinuity data |
| US7129942B2 (en) * | 2002-12-10 | 2006-10-31 | International Business Machines Corporation | System and method for performing domain decomposition for multiresolution surface analysis |
| US20060288756A1 (en) * | 2003-02-21 | 2006-12-28 | De Meurechy Guido D K | Method and apparatus for scanning corrosion and surface defects |
| US7889209B2 (en) * | 2003-12-10 | 2011-02-15 | Sensable Technologies, Inc. | Apparatus and methods for wrapping texture onto the surface of a virtual object |
| US7436988B2 (en) * | 2004-06-03 | 2008-10-14 | Arizona Board Of Regents | 3D face authentication and recognition based on bilateral symmetry analysis |
| US7586489B2 (en) * | 2005-08-01 | 2009-09-08 | Nvidia Corporation | Method of generating surface defined by boundary of three-dimensional point cloud |
-
2005
- 2005-04-29 EP EP05743388.0A patent/EP1756517B1/en active Active
- 2005-04-29 CN CN200580019487A patent/CN100576250C/en not_active Expired - Fee Related
- 2005-04-29 MX MXPA06012485A patent/MXPA06012485A/en active IP Right Grant
- 2005-04-29 EA EA200601977A patent/EA012579B1/en not_active IP Right Cessation
- 2005-04-29 AU AU2005241463A patent/AU2005241463C1/en not_active Ceased
- 2005-04-29 CA CA2564156A patent/CA2564156C/en active Active
- 2005-04-29 BR BRPI0509830-0A patent/BRPI0509830A/en not_active Application Discontinuation
- 2005-04-29 US US11/118,027 patent/US7352369B2/en active Active
- 2005-04-29 WO PCT/US2005/014937 patent/WO2005108923A2/en active Application Filing
-
2006
- 2006-11-23 NO NO20065397A patent/NO339362B1/en unknown
-
2008
- 2008-01-22 US US12/017,959 patent/US7576743B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| US20080129731A1 (en) | 2008-06-05 |
| US7352369B2 (en) | 2008-04-01 |
| EA012579B1 (en) | 2009-10-30 |
| CN100576250C (en) | 2009-12-30 |
| NO339362B1 (en) | 2016-12-05 |
| MXPA06012485A (en) | 2007-06-19 |
| EP1756517A4 (en) | 2011-01-26 |
| EA200601977A2 (en) | 2007-04-27 |
| EP1756517B1 (en) | 2020-05-13 |
| AU2005241463A1 (en) | 2005-11-17 |
| EA200601977A3 (en) | 2008-02-28 |
| NO20065397L (en) | 2007-01-24 |
| WO2005108923A2 (en) | 2005-11-17 |
| BRPI0509830A (en) | 2007-10-16 |
| CN101052997A (en) | 2007-10-10 |
| US7576743B2 (en) | 2009-08-18 |
| US20050246130A1 (en) | 2005-11-03 |
| AU2005241463B2 (en) | 2010-01-28 |
| WO2005108923A3 (en) | 2007-06-14 |
| CA2564156A1 (en) | 2005-11-17 |
| EP1756517A2 (en) | 2007-02-28 |
| AU2005241463C1 (en) | 2010-07-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2564156C (en) | System and method for approximating an editable surface | |
| Haber et al. | A general two‐dimensional, graphical finite element preprocessor utilizing discrete transfinite mappings | |
| RU2282242C2 (en) | Method for evening out a grid of nodes in accordance to characteristic features in digital image | |
| Woodward | Cross-sectional design of B-spline surfaces | |
| JPS63217473A (en) | Geometric model formation by cumulative translational sweep | |
| Liu et al. | Skinning cubic Bézier splines and Catmull-Clark subdivision surfaces | |
| Stahl et al. | Post-processing and visualization techniques for isogeometric analysis results | |
| Zhu et al. | Analytical solutions for sketch-based convolution surface modeling on the GPU | |
| Huhnt | Reconstruction of edges in digital building models | |
| De Paor | Bézier curves and geological design | |
| Hahmann et al. | Volume-preserving FFD for programmable graphics hardware | |
| Engleitner et al. | Lofting with patchwork B-splines | |
| Wu et al. | Representation and discretization of arbitrary surfaces for finite element shell analysis | |
| Sorokina et al. | An explicit quasi-interpolation scheme based on C1 quartic splines on type-1 triangulations | |
| Fang et al. | Reconstruction of smooth parametric surfaces from unorganized data points | |
| Kolb et al. | Fair surface reconstruction using quadratic functionals | |
| Angelidis et al. | Space deformations and their application to shape modeling | |
| Kuriyama et al. | Discrete parameterization for deforming arbitrary meshes | |
| Belyaev et al. | Bump Mapping for Isosurface Volume Rendering | |
| Wu et al. | An automatic three-dimensional geological engineering modeling method based on tri-prism | |
| US20050116951A1 (en) | Using runs of cells to traverse a ray through a volume | |
| Singh et al. | Accurate and fast algorithm for the plotting of contours using eight node quadrilateral meshes | |
| Bobach et al. | Natural neighbor concepts in scattered data interpolation and discrete function approx-imation | |
| Sud | Efficient computation of discrete voronoi diagram and homotopy-preserving simplified medial axis of a three-dimensional polyhedron | |
| Chen | Highlight lines for surface quality control and shape manipulation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request |