Surface Reconstruction From Non-parallel Curve Networks

Building surfaces from cross‐section curves has wide applications including bio‐medical modeling. Previous work in this area has mostly focused on connecting simple closed curves on parallel cross‐sections. Here we consider the more general problem where input data may lie on non‐parallel cross‐sections and consist of curve networks that represent the segmentation of the underlying object by different material or tissue types (e.g., skin, muscle, bone, etc.) on each cross‐section. The desired output is a surface network that models both the exterior surface and the internal partitioning of the object. We introduce an algorithm that is capable of handling curve networks of arbitrary shape and topology on cross‐section planes with arbitrary orientations. Our algorithm is simple to implement and is guaranteed to produce a closed surface network that interpolates the curve network on each cross‐section. Our method is demonstrated on both synthetic and bio‐medical examples.     

Curve interpolation in recursively generated B-spline surfaces over arbitrary topology

Recursive subdivision is receiving a great deal of attention in the definition of B-spline surfaces over arbitrary topology. The technique has recently been extended to generate interpolating surfaces with given normal vectors at the interpolated vertices. This paper describes an algorithm to generate recursive subdivision surfaces that interpolate B-spline curves. The control polygon of each curve is defined by a path of vertices of the polyhedral network describing the surface. The method consists of applying a one-step subdivision of the initial network and modifying the topology in the neighborhood of the vertices generated from the control polygons. Subsequent subdivisions of the modified network generate sequences of polygons each of which converges to a curve interpolated by the limit surface. In the case of regular networks, the method can be reduced to a knot insertion process.     

三维形体的形态学运算及应用

本文首先讨论了边界表示的物体的形态学运算问题。文中讨论了任意形状的物体，包括凸物体和非凸物体；引入了负形体的概念；并提出了２Ｄ和３Ｄ连续物体的Ｍｉｎｋｏｗｓｋｉ运算的统一算法。利用这一算法，提出了３Ｄ任意形体间的变形新算法，把边界表示的物体的形态学运算的概念扩展到任意曲线和任意曲面上，并且介绍了３Ｄ物空间中曲线的概念和３Ｄ物空间双参数曲面的概念。文中任意形体的３Ｄ变形算法能很好地解决具有不同拓扑结构的３Ｄ边界表示的物体之间的变形问题。     

Fair Surface Reconstruction Using Quadratic Functionals

An algorithm for surface reconstruction from a polyhedron with arbitrary topology consisting of triangular faces is presented. The first variant of the algorithm constructs a curve network consisting of cubic Bézier curves meeting with tangent plane continuity at the vertices. This curve network is extended to a smooth surface by replacing each of the networks facets with a split patch consisting of three triangular Bézier patches. The remaining degrees of freedom of the curve network and the split patches are determined by minimizing a quadratic functional. This optimization process works either for the curve network and the split patches separately or in one simultaneous step. The second variant of our algorithm is based on the construction of an optimized curve network with higher continuity. Examples demonstrate the quality of the different methods.     

Recovering the position and orientation of free-form objects fromimage contours using 3D distance maps

The accurate matching of 3D anatomical surfaces with sensory data such as 2D X-ray projections is a basic problem in computer and robot assisted surgery, In model-based vision, this problem can be formulated as the estimation of the spatial pose (position and orientation) of a 3D smooth object from 2D video images. The authors present a new method for determining the rigid body transformation that describes this match. The authors' method performs a least squares minimization of the energy necessary to bring the set of the camera-contour projection lines tangent to the surface. To correctly deal with projection lines that penetrate the surface, the authors consider the minimum signed distance to the surface along each line (i.e., distances inside the object are negative). To quickly and accurately compute distances to the surface, the authors introduce a precomputed distance map represented using an octree spline whose resolution increases near the surface. This octree structure allows the authors to quickly find the minimum distance along each line using best-first search. Experimental results for 3D surface to 2D projection matching are presented for both simulated and real data. The combination of the authors' problem formulation in 3D, their computation of line to surface distances with the octree-spline distance map, and their simple minimization technique based on the Levenberg-Marquardt algorithm results in a method that solves the 3D/2D matching problem for arbitrary smooth shapes accurately and quickly     

三维水下监视传感器网络的拓扑生成算法

当前大多数传感器网络研究假设传感器节点在二维平面部署,然而许多水下传感器网络应用要求节点在三维空间分布。针对三维水下监视应用,提出一种新颖的分布式传感器网络拓扑生成算法（ETG）。基于格理论,ETG算法将节点的移动控制与活动设备的调度相结合。初始时刻监视设备密集部署在二维海平面上,ETG算法根据局部信息选择活动设备,并通过控制其传感器节点在垂直方向的移动,形成一个三维水下网络。仿真实验表明,ETG算法能够以较小的平均节点移动距离覆盖较大的空间,从而有效减少网络建立阶段的能量开销。     

基于最佳圆弧样条逼近的快速距离曲面计算

距离曲面是一种常用的隐式曲面,它在几何造型和计算机动画中具有重要的应用价值,但以往往在对距离曲面进行多边形化时速较慢,为了提高点到曲线最近距离计算的效率,提出了一种基于最佳圆弧样条逼近的快速线骨架距离曲面计算方法,该算法对于一条任意的二维NURBS曲线,在用户给定的误差范围内,先用最少量的圆弧样条来逼近给定的曲线,从而把点到NURBS曲线最近距离的计算问题转化为点到圆弧样条最近距离的计算问题,由于在对曲面进行多边形化时,需要大量的点到曲线最近距离的计算,而该处可以将点到圆弧样条最近距离很少的计算量来解析求得,故该算法效率很高,该实验表明,算法简单实用,具有很大的应用价值。     

G interpolation of free form curve networks by biquintic Gregory patches

The problem of interpolating a free form curve network with irregular topology is investigated in order to create a curvature continuous surface. The spanning curve segments are parametric quintic polynomials, the interpolating surface elements are biquintic Gregory patches. A necessary compatibility condition is formulated and proved which need to be satisfied at each node of the curve network. Constraints are derived for assuring G² continuity between biquintic Gregory patches, which share a common side or a common corner point. The above conditions still leave certain geometric freedom for defining the entire G² surface, so following some analysis a particular construction is presented, by which after computing the principle curvatures at each node the free parameters are locally set for each interpolating Gregory patch.     

基于连接点的3D多角弧匹配

为了满足实际应用的要求，快速地进行多角弧匹配是必要的，因此，提出了一种快速的3D多角弧表示和匹配的方法。该方法是在分析多角弧几何形状基础上，引入连接点的概念，同时，通过在连接点处建立局部直角坐标系来得到每一连接点处的球面坐标。并用连接点的球面坐标集表示多角弧，因这一表示在旋转和平移变换下是不变的，所以可取该球面坐标集作为多角弧匹配的特征集，用该特征集可保持多角弧的几何属性和拓扑结构。这样3D多角弧匹配就降为1D数值串匹配，从而使使匹配变得简单快速。其测量函数为对应连接点间的均方差。实验结果表明，该匹配算法效果良好，并且对于数值污染具有健壮性。     

Decentralized control of traffic networks

An arbitrary interconnected network of intersecting streets with one-way traffic flow is considered in which some or all intersections are assumed to be operating at maximum capacity. In this case, it is desired to find a robust decentralized controller at each traffic intersection, so that for all perturbed external traffic flow sources/sinks, the resulting perturbed queues at each intersection are asymptotically "balanced". The only local knowledge required of the controllers at each intersection is that re the queue lengths of the intersection. It is shown that the necessary and sufficient conditions for a solution to the problem to exist depend solely on the topology of the network and the percentage of vehicles turning at each intersection, and are identical to the conditions required for a solution to exist for the centralized controller case. In particular, it is shown that if turns always occur at an intersection, then a solution to the problem always exists. A characterization of controllers which solve the problem are then given. In particular, it is shown that the local controllers which solve the problem are identical to each other, which implies that if the network is expanded (assuming a solution exists) then one can use the same local controllers to solve the problem. It is also shown that the local controllers can, in principle, give "perfect control" to the problem. Some examples are included to illustrate the results.     

Intersecting a freeform surface with a general swept surface

We present efficient and robust algorithms for intersecting a rational parametric freeform surface with a general swept surface. A swept surface is given as a one-parameter family of cross-sectional curves. By computing the intersection between a freeform surface and each cross-sectional curve in the family, we can solve the intersection problem. We propose two approaches, which are closely related to each other. The first approach detects certain critical points on the intersection curve, and then connects them in a correct topology. The second approach converts the intersection problem to that of finding the zero-set of polynomial equations in the parameter space. We first present these algorithms for the special case of intersecting a freeform surface with a ruled surface or a ringed surface. We then consider the intersection with a general swept surface, where each cross-sectional curve may be defined as a rational parametric curve or as an implicit algebraic curve.     

Simultaneous aligning and smoothing of surface triangulations

In this work we develop a procedure to deform a given surface triangulation to obtain its alignment with interior curves. These curves are defined by splines in a parametric space and, subsequently, mapped to the surface triangulation. We have restricted our study to orthogonal mapping, so we require the curves to be included in a patch of the surface that can be orthogonally projected onto a plane (our parametric space). For example, the curves can represent interfaces between different materials or boundary conditions, internal boundaries or feature lines. Another setting in which this procedure can be used is the adaption of a reference mesh to changing curves in the course of an evolutionary process. Specifically, we propose a new method that moves the nodes of the mesh, maintaining its topology, in order to achieve two objectives simultaneously: the piecewise approximation of the curves by edges of the surface triangulation and the optimization of the resulting mesh. We will designate this procedure as projecting/smoothing method and it is based on the smoothing technique that we have introduced for surface triangulations in previous works. The mesh quality improvement is obtained by an iterative process where each free node is moved to a new position that minimizes a certain objective function. The minimization process is done on the parametric plane attending to the surface piece-wise approximation and to an algebraic quality measure (mean ratio) of the set of triangles that are connected to the free node. So, the 3-D local projecting/smoothing problem is reduced to a 2-D optimization problem. Several applications of this method are presented.     

Polynomial surfaces interpolating arbitrary triangulations

Triangular Bezier patches are an important tool for defining smooth surfaces over arbitrary triangular meshes. The previously introduced 4-split method interpolates the vertices of a 2-manifold triangle mesh by a set of tangent plane continuous triangular Bezier patches of degree five. The resulting surface has an explicit closed form representation and is defined locally. In this paper, we introduce a new method for visually smooth interpolation of arbitrary triangle meshes based on a regular 4-split of the domain triangles. Ensuring tangent plane continuity of the surface is not enough for producing an overall fair shape. Interpolation of irregular control-polygons, be that in 1D or in 2D, often yields unwanted undulations. Note that this undulation problem is not particular to parametric interpolation, but also occurs with interpolatory subdivision surfaces. Our new method avoids unwanted undulations by relaxing the constraint of the first derivatives at the input mesh vertices: The tangent directions of the boundary curves at the mesh vertices are now completely free. Irregular triangulations can be handled much better in the sense that unwanted undulations due to flat triangles in the mesh are now avoided.     

Texture mapping on surfaces of arbitrary topology using norm preserving-based optimization

A simple and yet highly efficient, high-quality texture mapping method for surfaces of arbitrary topology is presented. The new method projects the given surface from the 3D object space into the 2D texture space to identify the 2D texture structure that will be used to texture the surface. The object space to texture space projection is optimized to ensure minimum distortion of the texture mapping process. The optimization is achieved through a commonly used norm preserving minimization process on edges of the surface. The main difference here is, by using an initial value approach, the optimization problem can be set up as a quadratic programming problem and, consequently, solved by a linear least squares method. Three methods to choose a good initial value are presented. Test cases show that the new method works well on surfaces of arbitrary topology, with the exception of surfaces with exceptionally abnormal curvature distribution. Other advantages of the new method include uniformity and seamlessness of the texture mapping process. The new method is suitable for applications that do not require precise texture mapping results but demand highly efficient mapping process such as computer animation or video games.     

Interactive 3D visualization for wireless sensor networks

Wireless sensor networks open up a new realm of ubiquitous computing applications based on distributed large-scale data collection by embedded sensor nodes that are wirelessly connected and seamlessly integrated within the environment. 3D visualization of sensory data is a challenging issue, however, due to the large number of sensors used in typical deployments, continuous data streams, and constantly varying network topology. This paper describes a practical approach for interactive 3D visualization of wireless sensor network data. A regular 3D grid is reconstructed using scattered sensor data points and used to generate view-dependent 2D slices that are consequently rendered with commodity graphics hardware leading to smooth visualization over space and time. Furthermore, the use of efficient space partitioning data structures and the independent processing of sensor data points facilitates interactive rendering for large number of sensors while accommodating constantly changing network topology. The practical value of the proposed method is demonstrated and results obtained for visualizing time-varying temperature distributions in an urban area are presented.     

G2 quasi-developable Bezier surface interpolation of two space curves

Surface development is used in many manufacturing planning operations, e.g., for garments, ships and automobiles. However, most freeform surfaces used in design are not developable, and therefore the developed patterns are not isometric to the original design surface. In some domains, the CAD model is created by interpolating two given space curves. In this paper, we propose a method to obtain a G² quasi-developable Bezier surface interpolating two arbitrary space curves. The given curves are first split into a number of piecewise Bezier curves and elemental Bezier patches each of which passes through four splitting points are constructed. All neighboring elemental patches are G² connected and they are assembled optimally in terms of the degree of developability (the integral Gaussian curvature). Experiments show that the final composite Bezier surface is superior to a lofted one which is defined regardless of the final surface developability.     

Complete Tensor Field Topology on 2D Triangulated Manifolds embedded in 3D

This paper is concerned with the extraction of the surface topology of tensor fields on 2D triangulated manifolds embedded in 3D. In scientific visualization topology is a meaningful instrument to get a hold on the structure of a given dataset. Due to the discontinuity of tensor fields on a piecewise planar domain, standard topology extraction methods result in an incomplete topological skeleton. In particular with regard to the high computational costs of the extraction this is not satisfactory. This paper provides a method for topology extraction of tensor fields that leads to complete results. The core idea is to include the locations of discontinuity into the topological analysis. For this purpose the model of continuous transition bridges is introduced, which allows to capture the entire topology on the discontinuous field. The proposed method is applied to piecewise linear three‐dimensional tensor fields defined on the vertices of the triangulation and for piecewise constant two or three‐dimensional tensor fields given per triangle, e.g. rate of strain tensors of piecewise linear flow fields.     

3D anatomical shape atlas construction using mesh quality preserved deformable models

3D anatomical shape atlas construction has been extensively studied in medical image analysis research, owing to its importance in model-based image segmentation, longitudinal studies and populational statistical analysis, etc. Among multiple steps of 3D shape atlas construction, establishing anatomical correspondences across subjects, i.e., surface registration, is probably the most critical but challenging one. Adaptive focus deformable model (AFDM) [1] was proposed to tackle this problem by exploiting cross-scale geometry characteristics of 3D anatomy surfaces. Although the effectiveness of AFDM has been proved in various studies, its performance is highly dependent on the quality of 3D surface meshes, which often degrades along with the iterations of deformable surface registration (the process of correspondence matching). In this paper, we propose a new framework for 3D anatomical shape atlas construction. Our method aims to robustly establish correspondences across different subjects and simultaneously generate high-quality surface meshes without removing shape details. Mathematically, a new energy term is embedded into the original energy function of AFDM to preserve surface mesh qualities during deformable surface matching. More specifically, we employ the Laplacian representation to encode shape details and smoothness constraints. An expectation–maximization style algorithm is designed to optimize multiple energy terms alternatively until convergence. We demonstrate the performance of our method via a set of diverse applications, including a population of sparse cardiac MRI slices with 2D labels, 3D high resolution CT cardiac images and rodent brain MRIs with multiple structures. The constructed shape atlases exhibit good mesh qualities and preserve fine shape details. The constructed shape atlases can further benefit other research topics such as segmentation and statistical analysis.     

实平面奇异代数曲线的全局B样条逼近

提出了一种用k次B样条曲线全局逼近实平面k次代数曲线的算法,每个连通部分用一条B样条曲线逼近.它适合于任意亏格的不可约的实平面代数曲线(包括含奇异点的曲线).这种逼近建立在所提出的代数曲线胀开采样的基础上,这种胀开采样算法从本质上解决了奇异点周围采样难的问题.实验结果表明,该方法的逼近精度高于已有算法.     

基于三维表面模型的任意切割算法*

通过对三维表面模型进行切割,人们可以方便地观察物体的内部结构。针对三维表面模型,提出了一种任意切割算法。任意切割过程首先是交互生成切割曲线,生成切割面;然后通过切割面对表面模型进行切割,在切割面上生成交线序列;再由交线序列与切割面边界生成封闭的边界轮廓,确定各边界轮廓间的包含关系;最后对边界轮廓包围的截面区域进行Delaunay三角剖分并着色,得到完整的剖面。实验结果证明了该算法的有效性和可行性。