Point cloud denoising is a fundamental yet challenging task in computer graphics. Existing solutions typically rely on supervised training on synthesized noise. However, real-world noise often exhibits greater complexity, causing learning-based methods trained on synthetic noise to struggle when encountering unseen noise-a phenomenon we refer to as noise misalignment. To address this challenge, we propose LaPDA (Latent-space Point cloud Denoising with Adaptivity), a neural network explicitly designed to mitigate noise misalignment and enhance denoising robustness. LaPDA consists of two key stages. First, we adaptively model noise in the latent space, aligning unseen noise distributions with the known training distributions or adjusting them toward distributions with lower noise scales. Training objectives at this stage are formulated based on controlled synthetic noise with varying intensity levels. Second, we introduce a gradual noise removal module that optimizes the spatial distribution of the adaptively adjusted noisy points. Extensive experiments conducted on both synthetic and scanned datasets demonstrate that LaPDA achieves enhanced accuracy and robustness compared to state-of-the-art methods.

In magnetic resonance imaging (MRI), anisotropic volumes with low through-plane resolution are typically acquired. Recently, diffusion models have shown strong performance in anisotropic MRI super-resolution (SR). However, iterative sam-pling limits the clinical applicability of diffusion models. To address this problem, we propose a One-Step Denoising Diffusion Probabilistic Model (OS-DDPM) for anisotropic MRI SR, which can generate isotropic High-Resolution (HR) data in a single sampling step. Firstly, we construct a Dual-Domain One-Step Generator (DDOS-Generator) comprising a student network for initial reconstruction in the spatial domain and a wavelet denoising module for noise suppression and detail refinement in the wavelet domain. Secondly, variational score distillation is applied to distill the generative prior from the pre-trained multi-step DDPM (teacher) to OS-DDPM, which can produce high-quality one-step SR data comparable to multi-step SR data. Finally, diffusion-based noise-aware contrastive learning is designed to bridge the distribution mismatch between one-step SR data and HR data. Extensive experiments on two public datasets demonstrate that the proposed method offers a fast inference speed and outperforms existing anisotropic SR methods and its teacher diffusion model in most metrics.

Probabilistic Partial Least Squares (PPLS) extends Partial Least Squares (PLS) by incorporating a probabilistic framework and identifiability constraints. PPLS models the relationship between two datasets as the sum of a joint component comprising correlated latent vectors and a noise component of isotropic normal random vectors. However, estimating PPLS model parameters involves a challenging nonconvex constrained optimization problem. We propose the Scalar Likelihood Method (SLM) to estimate the parameters of the PPLS model using a scalar likelihood function of the observed variables, derived via a novel rank n update technique. This technique avoids introducing additional constraints that would compromise the model’s dimension reduction properties. In addition, we derive a closed-form solution for the noise distribution in the observed data, significantly reducing parameter coupling in the objective function, thereby greatly improving parameter estimation quality. Through simulation studies, as well as association analysis and prediction tasks, we demonstrate the effectiveness and efficiency of SLM, highlighting its potential in practical applications.

Intracranial aneurysms (IAs) pose a significant risk due to their potential to rupture, leading to severe clinical outcomes. Accurate quantification of aneurysm morphology and assessment of rupture risk are crucial for timely intervention and treatment. In this study, we present an approach for the geometric quantification of IAs and clinical rupture risk evaluation based on the relationship between these measurements. Our approach includes a comprehensive set of geometric characteristics, such as height, width, neck width, arterial diameter, area, and volume. We applied these measurements to the public IntrA dataset, and to our knowledge, this is the first geometric analysis conducted on this dataset. By integrating these morphological characteristics with expert assessments, we developed a predictive model for IA rupture risk. Our findings reveal a strong correlation between aneurysm depth and rupture risk, with even stronger associations observed for higher-order dimension metrics like surface area and volume. This highlights the critical role of 3D automatic quantification in evaluating rupture risk. This research provides a foundation for further geometric analysis of IAs and offers potential advancements in automated diagnostics and precision medicine.

Saliency is an important characteristic in 3D analysis, and saliency detection for 3D meshes has been extensively studied in visual computing. However, insufficient feature description poses a significant challenge for 3D mesh saliency maps, which consist with human visual perception that is independent of the surroundings. We propose a multi-level fusion saliency detection method that leverages multiple eigenfunctions of the Laplacian-Beltrami operator(LBO) through a non-linear suppression operator, effectively integrating information across various levels. This method is based on reconstruction error and sparse matrix within the feature space constructed through dictionary learning (DL) at each level. The effectiveness of this method has been validated on the SchellingData, as demonstrated through both visualization and quantization. The method can be widely applied to down-stream tasks such as mesh simplification, object detection and so on.

Accurate segmentation of cerebral aneurysms from 3D vessel meshes is an essential yet challenging task, complicated by diverse imaging modalities that produce multi-resolution representations. Existing methods often struggle to handle meshes of varying granularity and orientations while maintaining segmentation accuracy. In this paper, we propose an end-to-end multi-scale graph-based segmentation framework that incorporates SE(3)-equivariance and an uncertainty-aware loss function. Our approach constructs a multi-scale graph representation on the 3D mesh and leverages a self-attention mechanism over graph edges to achieve adaptive neighborhood awareness, enabling the network to effectively handle multiple mesh resolutions simultaneously. By introducing an SE(3)-equivariant backbone, rotational variations in aneurysm orientation are naturally accommodated, ensuring relevant and effective feature learning. Furthermore, we develop an uncertainty-aware loss that adaptively emphasizes ambiguous regions, improving segmentation quality and confidence. Experimental results on the public datasets IntrA demonstrate that our method outperforms existing techniques, offering improved accuracy, consistency and stability across different mesh resolutions for 3D aneurysm segmentation. Code is available on https://github.com/Dolphin4mi/se3meshseg.

Tubular-like system shape analysis is quite difficult in geometry and topology, while it is widely used in plants and organs analysis in practice. However, traditional discrete representations such as voxels and point clouds often require substantial storage and may lead to the loss of fine-grained geometric and topological details. To address these challenges, we propose SE(3)-BBSCformerGCN, a novel framework for learning shape-aware representations from continuous tubular topological manifolds with equivariance to rotations and translations. Our approach leverages Ball B-Spline Curve (BBSC) to define tubular manifolds and its functional space. We provide a formal mathematical definition and analysis of the resulting manifolds and the BBSC functional space, and incorporate an equivariant mapping that preserves geometric and topological stability. Compared to the point cloud and voxel based representations, our manifold-based formulation significantly reduces data complexity while preserving geometric attributes together with topological features. We validate our method on the branch classification task for Circle of Willis (CoW) on the TopCoW 2024 dataset and the clinical dataset. Our method consistently outperforms voxel and point cloud based baselines in terms of classification performance, generalization ability, convergence speed, and robustness to overfitting.

ABSTRACT Accurate 3D craniofacial landmark detection is critical for applications in medicine and computer animation, yet remains challenging due to the complex geometry of craniofacial structures. In this work, we propose a coarse‐to‐fine framework for anatomical landmark localization on 3D craniofacial models. First, we introduce a Diffused Two‐Stream Network (DTS‐Net) for heatmap regression, which effectively captures both local and global geometric features by integrating pointwise scalar flow, tangent space vector flow, and spectral features in the Laplace‐Beltrami space. This design enables robust representation of complex anatomical structures. Second, we propose a heat kernel‐based energy optimization method to extract landmark coordinates from the predicted heatmaps. This approach exhibits strong performance across various geometric regions, including boundaries, flat surfaces, and high‐curvature areas, ensuring accurate and consistent localization. Our method achieves state‐of‐the‐art results on both a 3D cranial dataset and the BU‐3DFE facial dataset.

This article presents Diffusing Winding Gradients (DWG) for reconstructing watertight surfaces from unoriented point clouds. Our method exploits the alignment between the gradients of the screened generalized winding number (GWN) field—a robust variant of the standard GWN field—and globally consistent normals to orient points. Starting with an unoriented point cloud, DWG initially assigns a random normal to each point. It computes the corresponding screened GWN field and extracts a level set whose iso-value is the average of GWN values across all input points. The gradients of this level set are then utilized to update the point normals. This cycle of recomputing the screened GWN field and updating point normals is repeated until the screened GWN level sets stabilize and their gradients cease to change. Unlike conventional methods, DWG does not rely on solving linear systems or optimizing objective functions, which simplifies its implementation and enhances its suitability for efficient parallel execution. Experimental results demonstrate that DWG significantly outperforms existing methods in terms of runtime performance. For large-scale models with 10 to 20 million points, our CUDA implementation on an NVIDIA GTX 4090 GPU achieves speeds 30 to 120 times faster than iPSR, the leading sequential method, tested on a high-end PC with an Intel i9 CPU. Furthermore, by employing a screened variant of GWN, DWG demonstrates enhanced robustness against noise and outliers and proves effective for models with thin structures and real-world inputs with overlapping and misaligned scans. For source code and additional results, visit our project webpage: https://dwgtech.github.io/ .

A fast and fully automatic design of 3D cranial implants is highly desired in cranioplasty, and is key to the treatment of skull trauma. We have defined the repair of skull defects as a 3D shape completion task by proposing a two-stage diffusion model based on the representation of 3D shapes using signed distance function (SDF). Specifically, we design a diffusion model conditioned on partial shapes, we compress the 3D shape into a compact latent representation using the encoder in the vector quantized variational autoencoder (VQ-VAE) and learn the diffusion model based on this compressed discrete representation. Encoding the latent space with the autoencoder can achieve high-quality 3D cranial shape completion. In order to accurately capture local and fine-grained shape details, the training data is geometrically encoded from a compactly learned code-book. The two-stage diffusion generator with a coarse-to-fine approach possesses precise and expressive structural modeling capabilities to ensure the supplementation of detailed geometric information. Experimental results verified sufficient expressiveness of our model with generating high-fidelity results with fine-grained local details, outperforming the state-of-the-art methods.

Automated, accurate and robust segmentation of CT scans on abdominal organs remains a formidable challenge due to the blurred boundaries and subtle features of adjacent structures of abdominal organs. MedSAM, the fully fine-tuned model based on the Segment Anything Model (SAM) for medical image, fails to perform well relying solely on box prompt. Besides, utilizing the contextual information in the 3D volume data is tricky. To address above issues, we develop a MedSAM-based model for abdominal multi-organ CT segmentation. A feature selection network is introduced to extract entropy and curvature information. Together with the image embedding from the image encoder of MedSAM, the extracted information is processed during the inference. Low-Rank Adaptation (LoRA) is also applied to the image encoder and mask decoder of MedSAM for efficient fine-tuning on the FLARE2022 dataset. During the inference, mask prompts from segmentation results of last adjacent structure together with the expanded box prompt are used to improve the performance in continued context, thereby achieving an automated workflow. The experimental results suggest that the proposed method demonstrates a state-of-the-art performance and a potential clinical promise.

Unsupervised learning-based action recognition methods have shown immense potential by leveraging vast amounts of unlabeled data, yielding competitive performance in action recognition tasks. To capture semantic information in action, many efforts have focused on multi-stream fusion techniques, which makes the models heavier and less flexible to train. However, existing single-stream methods struggle to provide rich semantic information through data augmentation. To address these challenges, we propose a novel joint-level semantic augmentation method based on a constructed motion manifold, which achieves significant performance gains with virtually no additional training costs. Our approach introduces random perturbations on the motion manifold constructed based on evolutionary metrics. We propose several semantic augmentation strategies, including directional semantic perturbations, magnitude semantic perturbations, and initial pose perturbations. To pay greater attention to samples with significant semantic differences, we introduce the Semantic Adaptive Weighted Loss (SAWL) based on motion manifold distance. SAWL encourages the model to pay more attention to these samples, leading to the learning of an embedding space with semantic invariance. Extensive experiments were conducted on three large-scale datasets, i.e., NTU-60, NTU-120, and PKU-MMD II. The results demonstrate that our single-stream approach, empowered by joint-level semantic augmentation, achieves state-of-the-art (SOTA) performance among single-stream methods and competitive performance among multi-stream methods.

3D human motion prediction is a challenging task due to the highly non-linear nature of movements. Existing deep learning-based action prediction methods emphasize the design of sophisticated network architectures to achieve state-of-the-art performance on datasets. However, in real-life scenarios, changes in skeletal size lead to a shift in data distribution, which presents challenges for accurate motion prediction. Additionally, the presence of stretching artifacts in predicted bone sequences significantly impacts the quality of the predictions. To address these issues, we propose a framework that combines geometric encoding with neural networks to achieve size-invariant and high-quality motion prediction without stretching artifacts. We consider the constraint of bone length and construct a motion space using Riemannian manifold theory, which remains unaffected by changes in skeletal size and can fully represent human motion. Furthermore, we propose the trajectory transport square-root velocity function to encode motion sequences into a flattened space. This transformation simplifies the distance calculation and linearizes the optimization problem in non-flattened space. Experiments on the Human 3.6M and CMU MoCap datasets demonstrated that the proposed method has achieved competitive performance without any stretching artifacts and exhibits robustness to changes in skeletal size.

This paper presents a geometric model for achieving automatic 3D quantitative calculation of the femoroacetabular impingement (FAI) index on computed tomography (CT) images. The result is then used in subsequent work to reduce errors due to perspective differences in existing clinical measurements. This is the first automated quantitative method for a comprehensive FAI diagnostic index that does not rely on datasets. First, the geometric description and equation expression of the hip point cloud were obtained from CT images, fitting the key points and outlines required for the diagnostic index according to such geometric properties as Gaussian curvature. Then, an objective quantitative expression of the FAI index was calculated based on the 3D morphological definition. Next, we statistically analyzed 37 clinical data cases to verify our method’s effectiveness. Especially for the cases of hip dysplasia and acetabular boundary, there were unclear, strong correlations between manual and automatic measures (r = 0.88, ICC = 0.84).

Automatic detecting intracranial aneurysms (IAs) poses significant challenges due to their diversity, varying locations, and complex classifications by size, shape, and phenotype. Current shape-based IAs detection methods, while promising, often neglect the topological connectivity of IAs vertices and the variable traits of the aneurysm's neck, leading to fragmented detections. To address these issues, we present a mesh convolutional neural network based on gauge equivariant convolution to leverage the topological and geometric features of 3D mesh models. Our network comprises four key components: anisotropic message passing (AMP) on mesh surfaces, gauge equivariant convolution (GEC), vector-aware feature reconstruction (VFR), and a pooling-free convolutional architecture. AMP ensures accurate detection of IAs from surrounding vessels by utilizing topological connectivity and anisotropic relationships between mesh vertices. GEC offers rotational equivariance for consistently learning geometric features, improving feature learning stability and efficiency. VFR preserves the geometric and directional integrity of the vector features, enriching the representational capacity of the network. The pooling-free convolutional architecture captures local and global geometric nuances of 3D meshes, achieving precise IAs detection and producing sharper IAs boundaries. Tests on the IntrA dataset show our method outperforms the current best by 1.83% and 1.02% in mIoU and mDSC, respectively.

Traditional sign language teaching methods face challenges such as limited feedback and diverse learning scenarios. Although 2D resources lack reality sences, classroom teaching is constrained by a scarcity of teacher and methods based on VR and AR have relatively primitive interaction feedback mechanisms. This study proposes an innovative teaching model that uses real-time monocular vision and mixed reality technology. First, we introduce an improved hand-posture reconstruction method to achieve sign language semantic retention and real-time feedback. Second, a ternary system evaluation algorithm is proposed for a comprehensive assessment, maintaining good consistency with experts in sign language. Furthermore, we use mixed reality technology to construct a scenario-based 3D sign language classroom and explore the user experience of scenario teaching. Overall, this paper presents a novel teaching method that provides an immersive learning experience, advanced posture reconstruction, and precise feedback, achieving positive feedback on user experience and learning effectiveness.

Abstract Due to the widespread applications of curves on n‐dimensional spheres, fitting curves on n‐dimensional spheres has received increasing attention in recent years. However, due to the non‐Euclidean nature of spheres, curve fitting methods on n‐dimensional spheres often struggle to balance fitting accuracy and curve fairness. In this paper, we propose a new fitting framework, GLS‐PIA, for parameterized point sets on n‐dimensional spheres to address the challenge. Meanwhile, we provide the proof of the method. Firstly, we propose a progressive iterative approximation method based on geodesic least squares which can directly optimize the geodesic least squares loss on the n‐sphere, improving the accuracy of the fitting. Additionally, we use an error allocation method based on contribution coefficients to ensure the fairness of the fitting curve. Secondly, we propose an adaptive knot placement method based on geodesic difference to estimate a more reasonable distribution of control points in the parameter domain, placing more control points in areas with greater detail. This enables B‐spline curves to capture more details with a limited number of control points. Experimental results demonstrate that our framework achieves outstanding performance, especially in handling imbalanced data points. (In this paper, “sphere” refers to n‐sphere (n≥ 2) unless otherwise specified.)

Abstract Due to the widespread applications of 2‐dimensional spherical designs, there has been an increasing requirement of modeling on the 𝕊2 manifold in recent years. Due to the non‐Euclidean nature of the sphere, it has some challenges to find a method to represent 2D regions on 𝕊2 manifold. In this paper, a skeleton‐based representation method of regions on 𝕊2, disk B‐spline(DBSC) on 𝕊2 is proposed. Firstly, we give the definition and basic algorithms of DBSC on 𝕊2. Then we provide the calculation method of DBSC on 𝕊2, which includes calculating the boundary points, internal points and their corresponding derivatives. Based on that, we give some modeling methods of DBSC on 𝕊2, including approximation, deformation. In the end, some stunning application examples of DBSC on 𝕊2 are shown. This work lays a theoretical foundation for further applications of DBSC on 𝕊2.

Magnetic resonance imaging (MRI) scans often suffer from noise and low-resolution (LR), which affect the diagnosis and treatment results obtained for patients. LR images and noise come together with MRI, and the existing methods solve image super-resolution (SR) reconstruction and denoising tasks in a step-by-step manner, which influences the overall real distribution of the MRI data. In this paper, we present a simultaneous SR and denoising algorithm based on a stochastic normalizing flow (SNF), named the MR image SR and denoising model based on an SNF (SRDSNF). SRDSNF adds the encoded information of the input image as the conditional information to each reverse step of the stochastic normalizing flow, which realizes a consistent description of the spatial distribution between the reconstruction result and the input image. We introduce rangenull space decomposition and subsequence sampling strategies to enhance the consistency of the input and output data and increase the generation speed of the model. Simultaneous SR and denoising tasks experiment is carried out using the BrainWeb and NFBS datasets. The experimental results show that good SR and denoising results are obtained with fewer sampling steps, these results are consistent with the ground truths, and the structural similarity and peak signal-to-noise ratio of the results are also higher than those of the comparison methods. The proposed method demonstrates potential clinical promise.

We present an effective spectral matching method based on a shape association graph for finding region correspondences between two cel animation keyframes. We formulate the correspondence problem as an adapted quadratic assignment problem, which comprehensively considers both the intrinsic geometric and topology of regions to find the globally optimal correspondence. To simultaneously represent the geometric and topological similarities between regions, we propose a shape association graph (SAG), whose node attributes indicate the geometric distance between regions, and whose edge attributes indicate the topological distance between combined region pairs. We convert topological distance to geometric distance between geometric objects with topological features of the pairs, and introduce Kendall shape space to calculate the intrinsic geometric distance. By utilizing the spectral properties of the affinity matrix induced by the SAG, our approach can efficiently extract globally optimal region correspondences, even if shapes have inconsistent topology and severe deformation. It is also robust to shapes undergoing similarity transformations, and compatible with parallel computing techniques.

Image denoising is an important topic in the digital image processing field. This paper theoretically studies the validity of the classical non-local mean filter (NLM) for removing Gaussian noise from a novel statistic perspective. By regarding the restored image as an estimator of the clear image from the statistical view, we gradually analyse the unbiasedness and effectiveness of the restored value obtained by the NLM filter. Then, we propose an improved NLM algorithm called the clustering-based NLM filter (CNLM) that derived from the conditions obtained through the theoretical analysis. The proposed filter attempts to restore an ideal value using the approximately constant intensities obtained by the image clustering process. Here, we adopt a mixed probability model on a prefiltered image to generate an estimator of the ideal clustered components. The experimental results show that our algorithm obtains considerable improvement in peak signal-to-noise ratio (PSNR) values and visual results when removing Gaussian noise. On the other hand, the considerable practical performance of our filter shows that our method is theoretically acceptable as it can effectively estimates ideal images.

This paper presents an example-guided facial reconstruction method for creating a 3D refined human face model from sparse landmarks, with a given example dataset. It is challenging to rapidly and accurately generate a high-precision 3D face model from raw scan data with a simple setup and processing. The main roadblock is that existing 3D face databases are far from adequate to describe the full variability of faces. To address these problems, we analyse the characteristics of examples from a relatively small sample set for more reliable reconstruction knowledge, as well as a new landmark marking method to simplify the description of human face shape. Principal component analysis is used to extract the feature patterns from samples and simplify data representation. Then the 3D face model and the landmarks are correlated via a mapping matrix. An effective mapping algorithm is devised to learn the transformation relation from landmarks to 3D face shapes. Compared with existing methods, the proposed method can generate a high-precision 3D face model from sparse landmarks more accurately. The application and extensive experimental evaluations on the Chinese craniofacial database and FaceWarehouse database show that our method can achieve high accuracy, effectiveness and robustness in 3D face reconstruction.

Epidermal growth factor receptor (EGFR) plays an essential role in tumor cell proliferation, angiogenesis and apoptosis inhibition; it is a crucial factor leading to cancer occurrence. For example, EGFR tyrosine kinase inhibitors in treating lung cancer patients have an excellent therapeutic effect. Targeted therapy based on EGFR gene mutation is one of the mainstream lung cancer treatment methods. Recent studies have shown that pulmonary nodules' characteristics are associated with the mutant status of EGFR, which provides the possibility of using CT images of patients with pulmonary nodules to predict the mutant status of EGFR. This study used the deep learning algorithm to establish the EGFR mutation type prediction model based on CT image recognition. The data sets used for model training and testing included 121 labeled CT images from hospital patients with pulmonary nodules. The research results showed that the model could be used for the non-invasive EGFR mutation type based on CT images.

We present an automated treadmill featuring physiological-data-driven feedback-PhysioTreadmill, which allows its user to easily control running settings based on their physical condition and can also motivate them through real-time physiological computing. We developed a robust exercise intensity self-adaptive adjustment algorithm using physiological data processing to adjust the user's physiological state accurately. We also designed exergames with physiological-data-driven visual/audio feedback in PhysioTreadmill. With PhysioTreadmill, we can ideally increase exercise duration, and enhance exercise performance and safety. Two user studies involving 42 participants showed that PhysioTreadmill is user-friendly and can effectively extend users' training duration.

Abstract Curve extension is a useful tool in the computer‐aided design (CAD) community. A given B‐spline curve usually needs to be extended by another curve to reach one or more target points. In this work, we aim to enlarge the representation domain of the extending part to achieve the optimal extending result in the global solution space. Inspired by this, we have made three contributions. First, we use piecewise polynomial, that is, a nonuniform B‐spline, instead of one polynomial segment to extend the original curve. Compared to one polynomial segment, curves consisting of piecewise polynomial have stronger modeling ability and therefore expand the solution space of the problem. For extension to multiple target points, we are the first to directly give extension results based on all target points rather than extending to every target step by step. Third, we exploit the matrix representation of B‐splines to obtain an explicit solution for this extension problem. The detailed formula derivations and experimental results are provided to show the validity and effectiveness of our method.

Abstract In this article, we propose a novel 3D face modeling method which constructs a new 3D face model from a low‐dimensional feature space consisted of a large set of blend shapes based on the discrete shape space theory. The details of original face features are completely retained during the modeling process and a large number of new natural faces are constructed by several face samples. The optimization process of our method is independently decoupled for different facial attributes (identity, expression, and head pose), which improves the application flexibility and reduces the probability of it falling into a local optimal situation. The new facial data with new attributes are constructed based on the geodesic path search in discrete shape space with sufficient freedom and accuracy. In experiments and applications based on public databases (Helen, LFW, and CUFS), the modeling results show our method can provide high‐quality 3D face model, with enough freedom for face expression editing and natural facial expression animation from a small facial sample set.

The cover image is based on the Original Article Dynamic disk B-spline curves by Zhongke Wu et al., https://doi.org/10.1002/cav.1955.

Abstract A disk B‐spline curve (DBSC) is an extension of a B‐spline curve and is used to represent a two‐dimensional (2D) region. DBSC is a useful 2D geometric representation and is widely applied in the 2D art design area, such as computer calligraphy, 2D computer animation, and nonphotorealistic rendering. To enhance the flexibility of DBSC, in this article, we propose a physics‐based generalization of DBSC–dynamic DBSC (D‐DBSC), which extends the traditional DBSC in the time domain. We give the mathematical expression of D‐DBSC and prove its several mathematical properties. We derive the motion equations of D‐DBSC based on Lagrangian mechanics and investigate the motion equations when it is under linear geometric constraints. Last, a D‐DBSC physical simulation system based on finite difference method is presented.

We propose a novel framework for 3D facial similarity measurement and its application in facial organization. The construction of the framework is based on Kendall shape space theory. Kendall shape space is a quotient space that is constructed by shape features. In Kendall shape space, the shape features can be measured and is robust to similarity transformations. In our framework, a 3D face is represented by the facial feature landmarks model (FFLM), which can be regarded as the facial shape features. We utilize the geodesic in Kendall shape space to represent the FFLM similarity measurement, which can be regarded as the 3D facial similarity measurement. The FFLM similarity measurement is robust to facial expressions, head poses, and partial facial data. In our experiments, we compute the distance between different FFLMs in two public facial databases: FRGC2.0 and BosphorusDB. On average, we achieve a rank-one facial recognition rate of 98%. Based on the similarity results, we propose a method to construct the facial organization. The facial organization is a hierarchical structure that is achieved from the facial clustering by FFLM similarity measurement. Based on the facial organization, the performance of face searching in a large facial database can be improved obviously (about 400% improvement in experiments).

Segmenting the cerebral vessels precisely from the time-of-flight magnetic resonance angiography (TOF-MRA) images is important for the diagnosis and therapy of the cerebrovascular diseases. Since the complex structures of cerebral vessels, the current cerebrovascular segmentation algorithms based on statistical model have less accuracy for stenotic vessels and are quite time-consuming. In this paper, we propose a novel automatic cerebrovascular segmentation algorithm based on focused Multi-Gaussians (FMG) model and weighted 3D Markov Random Field. As far as our knowledge, this is the first time to adopt multi-Gaussians distributions as vascular model with the purpose of modeling the vascular tissue more accurately. Furthermore, the fitting range is narrowed to local region related to vessels in order to make the model focus on the vascular tissue and simplify the finite mixture model. To incorporate precise local character of images to the model, we design a new weighted 3D MRF by a weighted neighborhood system (W-NBS). Finally, the particle swarm optimization (PSO) algorithm of parameter estimation has been implemented parallelly based on GPUs and the execution speed was improved by about 70 times. The experimental results show that the algorithm can produce detailed segmentation results especially for stenotic vessels.

Flexible indoor scene synthesis is a popular topic in computer graphics and virtual reality research due to its wide-ranging applications in home design, games and automated robotics training. We propose a novel approach to automatic and flexible indoor scene synthesis using an energy-based method. We regard indoor scene synthesis as a multiple-object optimization problem with furniture location and orientation according to the user's intention, as a constraint on the energy of the optimization problem. Based on the relationship of objects, the embedded aesthetic criterion, the design criterion for proper placement and human movement in a scene, we design five energy functions, the overlap constraint, pairwise constraint, wall constraint, aisle constraint, angle constraint and penalty item, are proposed. We use a multi-object particle swarm intelligence optimization method with a Markov chain Monte Carlo algorithm to solve this optimization problem and obtain a Pareto-optimal solution. 3D gestures are used as the medium of interaction between the user and the system. Our method significantly enhances the existing weighted energy optimization method by allowing a joint optimization of various energy functions. The experiments confirm that all the energy functions can converge at the same time and that the proposed method obtains results superior to those of the weighted methods. The proposed method is general which can be used to obtain layouts for various kind of rooms with different furniture.

The 3D skull is a well preserved bone under the effect of fire, humidity, temperature changes, and it is a important biological characteristic in the fields of archaeology, forensic science and anthropology. In particular, measuring the 3D skull similarity is a challenging and meaningful task. 3D skulls are geometric models with multiple holes and complex topologies. It is difficult to correctly calculate the similarity between 3D skulls because the general 3D shape similarity measurement is sensitive to boundaries. In this paper, we provide an effective pipeline for measuring the 3D skull similarity by calculating the cosine distance between the harmonic wave kernel signature (HWKS) values of 3D skulls. Based on the wave kernel signature, the HWKS is a shape descriptor which is involved the Laplace-Beltrami operator that can effectively extract geometrical and topological information from the 3D skulls. And the HWKS simultaneously describes the local and global properties of a skull compared to the wave kernel signature. In addition, our method is more flexible, and can be generalized to other 3D shapes. Several experiments show our method achieves good results and can correctly calculate the similarity between 3D skulls.

We propose a skeleton tree based method for classifying and retrieving non-rigid shapes. Firstly, based on the extracted skeletons of a non-rigid shape and geodesic distance computation, the center point in skeleton is defined and detected. Then, a skeleton tree is constructed based on the connection between the center point and other discrete points in skeleton. After that, a correspondence between the skeleton tree and the area distribution of the non-rigid shape is established. The skeleton tree features are achieved. The advantages of our method can be summarized as follows: (1) Scale-Invariant; (2) Low computational complexity; (3) Automatic topology repair. The experimental results show that our method is more accurate than existing methods.

Abstract Team formation with realistic crowd simulation behavior is a challenge in computer graphics, multiagent control, and social simulation. In this study, we propose a framework of crowd formation via hierarchical planning, which includes cooperative‐task, coordinated‐behavior, and action‐control planning. In cooperative‐task planning, we improve the grid potential field to achieve global path planning for a team. In coordinated‐behavior planning, we propose a time–space table to arrange behavior scheduling for a movement. In action‐control planning, we combine the gaze‐movement angle model and fuzzy logic control to achieve agent action. Our method has several advantages. (1) The hierarchical architecture is guaranteed to match the human decision process from high to low intelligence. (2) The agent plans his behavior only with the local information of his neighbor; the global intelligence of the group emerges from these local interactions. (3) The time–space table fully utilizes three‐dimensional information. Our method is verified using crowds of various densities, from sparse to dense, employing quantitative performance measures. The approach is independent of the simulation model and can be extended to other crowd simulation tasks.

Facial expression editing in face sketch is an important and challenging problem in computer vision community as facial animation and modeling. For criminal investigation and portrait drawing, automatic expression editing tools for face sketch improve work efficiency obviously and reduce professional requirements for users. In this paper, we propose a novel method for facial expression editing in face sketch using shape space theory. The new facial expressions in the sketch images can be regenerated automatically. The method includes two components: 1) face sketch modeling; 2) expression editing. The face sketch modeling constructs 3D face sketch data from 3D facial database to match the 2D face sketch. Using facial landmarks, the "shape" of the face sketch is represented in shape space. The shape space is a manifold space which removes the rigid transform group. In shape space, the accurate 3D face sketch model is obtained which is consistent to the original 2D face sketch. For expression editing, we change the parameters of 3D face sketch model in the shape space to obtain new expressions. The expression transfer in 3D face sketch model can be mapped into the 2D face sketch. The advantages of our method are: full-automatic in modeling process; no requirements of drawing skills to user and friendly interaction; robustness to head poses and different scales. In experiments, we use the 3D facial database, FaceWareHouse, to construct the 3D face sketch model and use face sketch images from database: CUHK Face sketch Database (CUFS) to show the performance of expression editing. Experimental results demonstrate that our method can effectively edit facial expressions in face sketch with high consistency and fidelity.

Scattered data fitting is always a challenging problem in the fields of geometric modeling and computer aided design. As the skeleton based three-dimensional solid model representation, the Ball B-Spline Curve is suitable to fit the tubular scattered data points. We study the problem of fitting the scattered data points with Ball B-spline curves (BBSCs) and propose the corresponding fitting algorithm based on the Particle Swarm Optimization (PSO) algorithm. In this process, we face three critical and difficult sub problems: (1) parameterization of the data points, (2) determination of the knot vector and (3) calculation of the control radii. All of them are multidimensional and nonlinear, especially the calculation of the parametric values. The parallelism of the PSO algorithm provides a high optimization, which is more suitable for solving nonlinear, nondifferentiable and multi-modal optimization problems. So we use it to solve the scattered data fitting problem. The PSO is applied in three steps to solve them. Firstly, we determine the parametric values of the data points with PSO. Then we compute the knot vector based on the parametric values of the data points. At last, we get the radius function. The experiments on the shell surface, the crescent surface and the real-world models verify the accuracy and flexibility of the method. The research can be widely used in the computer aided design, animation and model analysis.

The tsunami of video services brings huge pressure on communication systems, especially on some bandwidth-limited scenarios like wireless networks. Represented by H.264, current video compression approaches can save the transmission bandwidth by reducing the intra- and inter-frame redundancy. However, the huge amount of video content redundancy (VCR) in the transmission is still existing. In this paper, we propose a novel Content-Slimming System (CSS) framework based on the convergence of Computing, Communication and Cache to avoid the transmission of VCRs. The main idea of CSS is to detect VCRs, generate VCR models and clip them from the original frames by Computing, then transmit the necessary video content and semantic difference description by Communication, finally reconstruct the full video based on the received video content, semantic difference description and VCR models stored by Cache. Moreover, we also investigate the Video Monitoring application based on the CSS framework (CSS-VM) in which the background of monitoring frames is modeled, detected and sheared to reduce the transmission bandwidth consumption. The simulation results show that the bandwidth consumption of CSS-VM can be reduced at least by half compared to H.264, while the video quality and visual experience of CSS-VM are even better.

The analysis of cerebrovascular shape is important for the diagnose and pathologic identification. But as the limitation of the segmentation algorithm, the complete cerebrovascular volume data are difficult to obtain. So the triangle mesh of the vessel model generated for the medical images may appear many gaps. In the paper, we present a extension algorithm for Ball B-Spline curve with G2 continuity to repair the cerebrovascular structure from time-of-flight (TOF) magnetic resonance angiography (MRA) data. Ball B-Spline curve has its distinct advantages in representing a 3D tube like organs. A ball Bezier segment is used to construct the extending part and G2-continuity is applied to describe the smoothness at the joints. Fairness of the extending ball Bezier curve segment is achieved by minimizing energy objective functions for the center curve and the radius function separately. New control balls are computed by unclamping algorithm to represent the whole extended ball B-Spline curve. The experimental results demonstrate the effectiveness of our algorithm. The final results show that the proposed method provides good blending result, especially for those blood vessels of small size.

With the development of computer techniques and computer graphics, computer aided design (CAD) technology has been applied widely in the textile printing industry. There are many good methods for simulating weft knitted fabric, however, no work on knitted fabric had been performed using ball B-spline curves. In this paper, we have proposed a novel algorithm based on ball B-spline curves and Hooke's Law for modelling weft knitted fabric in 3D. Ball B-spline curves are based on B-spline curves that use data balls instead of data points. By changing the position and the radius of balls of the ball B-spline curves, we can change the thickness and position of the yarn. The loop is the basic unit of knitted fabric, it requires 7 data points to be represented by a ball B-spline curve, and the model is notably fast. At the same time, the model maintains both its accuracy and the geometrical and topological features of the fabric. With VC++ and Open GL, we added texture and light to simulate different patterns of weft knitted fabric in 3D, including stockinette knit, garter knit, 2-2 rib knit, wing knit, heart knit and stripe knit. Because of the nested structure of knitted fabric, each data point of each loop has changes in a different way when a force is applied to the knitted fabric. The deformation and force on the knitted fabric is modelled using Hooke's Law, which results in a more realistic and faster simulation. We use the parameterization method, and the parameters, such as the radius of ball, the height of a loop, the number of courses, the number of wales and the magnitude of the force, are easy to change. This method is a new means of naturally building and simulating weft knitted fabric in a CAD system.

In the field of computer graphics and human simulation, hair simulation is one of the most challenging physics and rendering problem. This paper presents a curly hair modeling utilizing the Nelder-Mead method and being based on the ball B-Spline Curves (BBSCs). The ball B-Spline Curves (BBSCs) is used to reconstruct hair model and the Nelder-Mead method is for calculating the equilibrium shape of each hair. The main advantages of our method used in this experiment are: (1) We construct hair model for each individual hair. (2) It can provide a wide range of curly hair types with the BBSCs for its good properties that it is flexible for modifying, editing and deforming. (3) The method can obtain efficient curly hair type because of the fast calculation and easy implementation of Nelder-Mead method. The hypothesis in this paper is supported by several credible evidences and the model can be widely used in the similar systems (e.g. Fabrics, green fields).

We present an automatic statistical intensity-based approach to extract the 3D cerebrovascular structure from time-of flight (TOF) magnetic resonance angiography (MRA) data. We use the finite mixture model (FMM) to fit the intensity histogram of the brain image sequence, where the cerebral vascular structure is modeled by a Gaussian distribution function and the other low intensity tissues are modeled by Gaussian and Rayleigh distribution functions. To estimate the parameters of the FMM, we propose an improved particle swarm optimization (PSO) algorithm, which has a disturbing term in speeding updating the formula of PSO to ensure its convergence. We also use the ring shape topology of the particles neighborhood to improve the performance of the algorithm. Computational results on 34 test data show that the proposed method provides accurate segmentation, especially for those blood vessels of small sizes.

Touch less interaction has received considerable attention in recent years with benefit of removing the burden of physical contact. To achieve mid-air interaction, several strategies are available. However, since most of these techniques directly map the 2D WIMP GUI to 3D user interface, they lead unnatural result. In this paper, interaction gestures and tools for exploring volume dataset are designed to perform the similar tasks in the real world. We mainly employ the idea of focus + context based on GPU volume ray casting by trapezoid-shaped transfer function when designing interaction tools. User studies are conducted to demonstrate the usability and intuitiveness of our method. The experimental results show a significant advantage in completion time after a short period of training.

Curve blending is an essential task in geometric modeling, while a ball B-spline curve (BBSC) has its advantages in representing freeform tubular objects. This paper proposes a blending algorithm for ball B-Spline curve with G2 continuity, which is used to describe the smoothness of the joint point. An original BBSC is extended smoothly to join another one, such that no additional blending curve is created and the two original curves are not changed. The shape of the extended curve is then determined by minimizing strain energy. The corresponding scalar function of the control balls is determined through applying G2-continuity conditions to the scalar function. In order to ensure the radii of the control balls are positive, we make a decision about the range of the G2-continuity parameter and then determine it by minimizing the strain energy in the affected area. The experiment results demonstrate our method for blending BBSC is effective. Moreover, some G2 blending results of the BBSC in simulating the tubular objects are given.

Curve extension is a useful function in CAD systems. Disk B-Spline curve has its distinct advantages in representing a 2D region. This paper presents an algorithm for extending the disk B-Spline curve. A disk Bezier segment is used to construct the extending part and G2-continuity can be used to describe the smoothness of the joint disk. Fairness of the extending disk Bezier segment is achieved by minimizing an energy objective function. New control disks are computed by unclamping algorithm to represent the whole extended disk B-Spline curve. The experimental results demonstrate the effectiveness of our method.

Curve extension is a useful function in shape modeling for cyber worlds, while a ball B-spline curve (BBSC) has its advantages in representing freeform tubular objects. In this paper, an extension algorithm for ball B-Spline curve with G2-continuity is investigated. We apply the extending method of B-Spline curves to the center curve of BBSC through generalizing a minimal strain energy method from 2D to 3D. And the initial value of G2-continuity parameter was selected by minimizing the approximate energy function which is a problem with O(1) time complexity. The corresponding scalar function of control balls is determined through applying G2-continuity condition to scalar function. In order to ensure the radii of the control balls are positive, we make a decision about the range of the G2-continuity parameter and then determine it by minimizing the strain energy in the affected area. Some experiments for comparing our method with other methods are given. The results show our method for extending BBSC is effective.