Automatic 4D Reconstruction of Patient-Specific Cardiac Mesh with 1-to-1 Vertex Correspondence from Segmented Contours Lines

We propose an automatic algorithm for the reconstruction of patient-specific cardiac mesh models with 1-to-1 vertex correspondence. In this framework, a series of 3D meshes depicting the endocardial surface of the heart at each time step is constructed, based on a set of border delineated magnetic resonance imaging (MRI) data of the whole cardiac cycle. The key contribution in this work involves a novel reconstruction technique to generate a 4D (i.e., spatial–temporal) model of the heart with 1-to-1 vertex mapping throughout the time frames. The reconstructed 3D model from the first time step is used as a base template model and then deformed to fit the segmented contours from the subsequent time steps. A method to determine a tree-based connectivity relationship is proposed to ensure robust mapping during mesh deformation. The novel feature is the ability to handle intra- and inter-frame 2D topology changes of the contours, which manifests as a series of merging and splitting of contours when the images are viewed either in a spatial or temporal sequence. Our algorithm has been tested on five acquisitions of cardiac MRI and can successfully reconstruct the full 4D heart model in around 30 minutes per subject. The generated 4D heart model conforms very well with the input segmented contours and the mesh element shape is of reasonably good quality. The work is important in the support of downstream computational simulation activities.     

Design of Orthopedic Plates and Its Modification Based on Feature

To quickly construct the orthopedic plates and to conveniently edit it, a novel method for designing the plates is put forward based on feature idea and parameterization. Firstly, attached to the existing or repaired bone model, the region of interest (ROI) is selected as the abutted surface of orthopedic plate, and the ROI is reconstructed to form a CAD surface. Secondly, the CAD surface is to be defined as a surface feature (SF) and then some semantic parameters are configured for it. Lastly, the plate body is constructed through thickening, and some higher parameters are defined for it so as to produce a volumetric feature (VF). In the above process, there exist two main problems: one is parameterization of the abutted surface, and the other is construction of the outer surface. Besides, the mapping relationship has to be built between surface feature parameters and volumetric feature parameters. This method supports the modification of high-level parameters, consequently promoting the quality and efficiency of orthopedic plate design.     

An Improved Approach of Mesh Segmentation to Extract Feature Regions

The objective of this paper is to extract concave and convex feature regions via segmenting surface mesh of a mechanical part whose surface geometry exhibits drastic variations and concave-convex features are equally important when modeling. Referring to the original approach based on the minima rule (MR) in cognitive science, we have created a revised minima rule (RMR) and presented an improved approach based on RMR in the paper. Using the logarithmic function in terms of the minimum curvatures that are normalized by the expectation and the standard deviation on the vertices of the mesh, we determined the solution formulas for the feature vertices according to RMR. Because only a small range of the threshold parameters was selected from in the determined formulas, an iterative process was implemented to realize the automatic selection of thresholds. Finally according to the obtained feature vertices, the feature edges and facets were obtained by growing neighbors. The improved approach overcomes the inherent inadequacies of the original approach for our objective in the paper, realizes full automation without setting parameters, and obtains better results compared with the latest conventional approaches. We demonstrated the feasibility and superiority of our approach by performing certain experimental comparisons.     

Anatomical Feature Segmentation of Femur Point Cloud Based on Medical Semantics

Feature segmentation is an essential phase for geometric modeling and shape processing in anatomical study of human skeleton and clinical digital treatment of orthopedics. Due to various degrees of freedom of bone surface, the existing segmentation algorithms can hardly meet specific medical need. To address this, a novel segmentation methodology for anatomical features of femur model based on medical semantics is put forward. First, anatomical reference objects (ARO) are created to represent typical characteristics of femur anatomy by 3D point fitting in combination with medical priori knowledge. Then, local point clouds between adjacent anatomies are selected according to the AROs to extract boundary feature point (BFP)s. Finally, the complete model of femur is divided into anatomical regions by executing the enhanced watershed algorithm guided with BFPs. Experimental results show that the proposed method has the advantages of automatic segmentation of femoral head, neck and other complex areas, and the segmentation results have better medical semantics. In addition, the slight modification of segmentation results can be achieved by adjusting a few threshold parameter values, which improves the convenience of modification for ordinary users.     

A spatial mechanism with higher pairs for modelling the human knee joint

By generalizing a previous model proposed in the literature, a new spatial kinematic model of the knee joint passive motion is presented. The model is based on an equivalent spatial parallel mechanism which relies upon the assumption that fibers within the anterior cruciate ligament (ACL), the medial collateral ligament (MCL) and the posterior cruciate ligament (PCL) can be considered as isometric during the knee flexion in passive motion (virtually unloaded motion). The articular surfaces of femoral and tibial condyles are modelled as 3-D surfaces of general shapes. In particular, the paper presents the closure equations of the new mechanism both for surfaces represented by means of scalar equations that have the Cartesian coordinates of the points of the surface as variables and for surfaces represented in parametric form. An example of simulation is presented in the case both femoral condyles are modelled as ellipsoidal surfaces and tibial condyles as spherical surfaces. The results of the simulation are compared to those of the previous models and to measurements. The comparison confirms the expectation that a better approximation of the tibiofemoral condyle surfaces leads to a more accurate model of the knee passive motion.     

Surface reconstruction of torsos with and without scoliosis

Visible surface asymmetries such as uneven shoulders, waist and hips, shoulder height differences and a shoulder blade prominence are often the most troublesome features associated with adolescent scoliosis. Treatment considerations are influenced by the severity and changes over time of these asymmetries. Outcomes are judged on how well the asymmetries are improved towards a normal trunk shape. In this paper, a deformable self organizing feature map (SOFM) is used as a geometric surface reconstruction tool to model the torso surface of subjects with and without scoliosis. The proposed parameterization technique provides a means of quantifying the surface asymmetries and assessing the changes due to either natural history or the effects of treatment. For evaluation 10 control subjects without scoliosis and 10 adolescents with scoliosis were scanned and their torsos were reconstructed. This preliminary study demonstrates that in around 5 min a torso scan with 60,000 data points can be transformed into a 2562 nodes mesh using SOFM. The accuracy of the final mesh is around 1.40 mm on average. The high accuracy and speed of this technique, makes it well suitable to be used in a clinical setting to assess surface features of subjects with scoliosis.     

Minimization of the gliding index: criterion for the generation of the surfaces of a knee endoprosthesis

Previous studies on removed failed artificial knees revealed significant degradation of the articular surfaces, including pitting and shredding, as well as burnishing accompanied by score marks and scratches, the latter damage group being related to the gliding motion of the joint. In an attempt to introduce an improved version of an artificial knee joint, we have proposed a general model by which the opposing surfaces of the prosthesis components can be synthesized. The criterion applied was that of minimization of a defined gliding index. The femoral condyles in this model were expressed in terms of torus geometry, and the kinematics of motion fed into the model was that of normal motion in the sagittal plane, including angles as well as the displacement vector in the knee joint. Geometry of the tibial component was obtained from the tangents of the femoral surface, in subsequent positions of motion. The optimal surfaces were those for which the gliding index assumed a minimal value. The solutions obtained for various input motions are presented and discussed.     

From high-resolution CT data to finite element models: development of an integrated modular framework

This work introduces a novel method of automating the process of patient-specific finite element (FE) model development using a mapped mesh technique. The objective is to map a predefined mesh (template) of high quality directly onto a new bony surface (target) definition, thereby yielding a similar mesh with minimal user interaction. To bring the template mesh into correspondence with the target surface, a deformable registration technique based on the FE method has been adopted. The procedure has been made hierarchical allowing several levels of mesh refinement to be used, thus reducing the time required to achieve a solution. Our initial efforts have focused on the phalanx bones of the human hand. Mesh quality metrics, such as element volume and distortion were evaluated. Furthermore, the distance between the target surface and the final mapped mesh were measured. The results have satisfactorily proven the applicability of the proposed method.     

Automated hexahedral meshing of anatomic structures using deformable registration

This work introduces a novel method of automating the process of patient-specific finite element (FE) model development using a mapped mesh technique. The objective is to map a predefined mesh (template) of high quality directly onto a new bony surface (target) definition, thereby yielding a similar mesh with minimal user interaction. To bring the template mesh into correspondence with the target surface, a deformable registration technique based on the FE method has been adopted. The procedure has been made hierarchical allowing several levels of mesh refinement to be used, thus reducing the time required to achieve a solution. Our initial efforts have focused on the phalanx bones of the human hand. Mesh quality metrics, such as element volume and distortion were evaluated. Furthermore, the distance between the target surface and the final mapped mesh were measured. The results have satisfactorily proven the applicability of the proposed method.     

Adaptive generation of multimaterial grids from imaging data for biomedical Lagrangian fluid–structure simulations

Spatial discretization of complex imaging- derived fluid–solid geometries, such as the cardiac environment, is a critical but often overlooked challenge in biomechanical computations. This is particularly true in problems with Lagrangian interfaces, where the fluid and solid phases share a common interface geometrically. For simplicity and better accuracy, it is also highly desirable for the two phases to have a matching surface mesh at the interface between them. We outline a method for solving this problem, and illustrate the approach with a 3D fluid–solid mesh of the mouse heart. An MRI dataset of a perfusion-fixed mouse heart with 50μm isotropic resolution was semi-automatically segmented using a customized multimaterial connected-threshold approach that divided the volume into non-overlapping regions of blood, tissue, and background. Subsequently a multimaterial marching cubes algorithm was applied to the segmented data to produce two detailed, compatible isosurfaces, one for blood and one for tissue. Both isosurfaces were simultaneously smoothed with a multimaterial smoothing algorithm that exactly conserves the volume for each phase. Using these two isosurfaces, we developed and applied novel automated meshing algorithms to generate anisotropic hybrid meshes on arbitrary biological geometries with the number of layers and the desired element anisotropy for each phase as the only input parameters. Since our meshes adapt to the local feature sizes and include boundary layer prisms, they are more efficient and accurate than non-adaptive, isotropic meshes, and the fluid–structure interaction computations will tend to have relative error equilibrated over the whole mesh.     

Merging of intersecting triangulations for finite element modeling

Surface mesh generation over intersecting triangulations is a problem common to many branches of biomechanics. A new strategy for merging intersecting triangulations is described. The basis of the method is that object surfaces are represented as the zero-level iso-surface of the distance-to-surface function defined on a background grid. Thus, the triangulation of intersecting objects reduces to the extraction of an iso-surface from an unstructured grid. In a first step, a regular background mesh is constructed. For each point of the background grid, the closest distance to the surface of each object is computed. Background points are then classified as external or internal by checking the direction of the surface normal at the closest location and assigned a positive or negative distance, respectively. Finally, the zero-level iso-surface is constructed. This is the final triangulation of the intersecting objects. The overall accuracy is enhanced by adaptive refinement of the background grid elements. The resulting surface models are used as support surfaces to generate three-dimensional grids for finite element analysis. The algorithms are demonstrated by merging arterial branches independently reconstructed from contrast-enhanced magnetic resonance images and by adding extra features such as vascular stents. Although the methodology is presented in the context of finite element analysis of blood flow, the algorithms are general and can be applied in other areas as well.     

基于机器学习的蛋白质编码区识别

针对DNA序列编码区的识别问题,本研究提出一个特征向量和逻辑回归的组合模型。首先对DNA序列进行数值处理转化为特征向量,并结合k字符相对频率技术提取特征向量的元素特征,之后利用二分类逻辑回归算法,对编码区和非编码区进行准确区分。选取了HMR195和BG570两个基准数据集进行五折交叉验证,结果表明,平均AUC(Area Under Curve)值分别为0.981 3和0.987 4,明显优于传统的贝叶斯判别法和VOSSDFT等方法。此外,本文提出的特征向量的维度很低,提高了运算效率。因此,本文组合模型能够较为高效准确地识别蛋白质编码区。     

A Time-Series-Based Feature Extraction Approach for Prediction of Protein Structural Class

This paper presents a novel feature vector based on physicochemical property of amino acids for prediction protein structural classes. The proposed method is divided into three different stages. First, a discrete time series representation to protein sequences using physicochemical scale is provided. Later on, a wavelet-based time-series technique is proposed for extracting features from mapped amino acid sequence and a fixed length feature vector for classification is constructed. The proposed feature space summarizes the variance information of ten different biological properties of amino acids. Finally, an optimized support vector machine model is constructed for prediction of each protein structural class. The proposed approach is evaluated using leave-one-out cross-validation tests on two standard datasets. Comparison of our result with existing approaches shows that overall accuracy achieved by our approach is better than exiting methods.     

Prostate segmentation on T2 MRI using Optimal Surface Detection

This paper deals with a T2 MRI prostate segmentation method. We assume to have an initial surface mesh obtained interactively or after a first rough segmentation. The surface of the prostate is then searched within the initial mesh neighborhood using the Optimal Surface Detection algorithm (OSD). This algorithm is based on the construction of a directed graph from the information obtained around the initial mesh. The optimal surface is then obtained by a graph cut. Three different cost functions for the graph have been explored, one based on the local gradient, another on a statistical model of shape and a third on a model of gradient profile. The parameters of this method have been tuned on 33 different T2 MRI volumes.     

An eFTD-VP framework for efficiently generating patient-specific anatomically detailed facial soft tissue FE mesh for craniomaxillofacial surgery simulation

Accurate surgical planning and prediction of craniomaxillofacial surgery outcome requires simulation of soft tissue changes following osteotomy. This can only be achieved by using an anatomically detailed facial soft tissue model. The current state-of-the-art of model generation is not appropriate to clinical applications due to the time-intensive nature of manual segmentation and volumetric mesh generation. The conventional patient-specific finite element (FE) mesh generation methods are to deform a template FE mesh to match the shape of a patient based on registration. However, these methods commonly produce element distortion. Additionally, the mesh density for patients depends on that of the template model. It could not be adjusted to conduct mesh density sensitivity analysis. In this study, we propose a new framework of patient-specific facial soft tissue FE mesh generation. The goal of the developed method is to efficiently generate a high-quality patient-specific hexahedral FE mesh with adjustable mesh density while preserving the accuracy in anatomical structure correspondence. Our FE mesh is generated by eFace template deformation followed by volumetric parametrization. First, the patient-specific anatomically detailed facial soft tissue model (including skin, mucosa, and muscles) is generated by deforming an eFace template model. The adaptation of the eFace template model is achieved by using a hybrid landmark-based morphing and dense surface fitting approach followed by a thin-plate spline interpolation. Then, high-quality hexahedral mesh is constructed by using volumetric parameterization. The user can control the resolution of hexahedron mesh to best reflect clinicians’ need. Our approach was validated using 30 patient models and 4 visible human datasets. The generated patient-specific FE mesh showed high surface matching accuracy, element quality, and internal structure matching accuracy. They can be directly and effectively used for clinical simulation of facial soft tissue change.     

Does a micro-grooved trunnion stem surface finish improve fixation and reduce fretting wear at the taper junction of total hip replacements? A finite element evaluation

The generation of particulate debris at the taper junction of total hip replacements (THRs), can cause failure of the artificial hip. The taper surfaces of femoral heads and trunnions of femoral stems are generally machined to a certain roughness to enhance fixation. However, the effect of the surface roughness of these surfaces on the fixation, wear and consequently clinical outcomes of the design is largely unknown. In this study, we asked whether a micro-grooved trunnion surface finish (1) improves the fixation and (2) reduces the wear rate at the taper junction of THRs. We used 3D finite element (FE) models of THRs to, firstly, investigate the effect of initial fixation of a Cobalt-Chromium femoral head with a smooth taper surface mated with a Titanium (1) micro-grooved and (2) smooth, trunnion surface finishes. Secondly, we used a computational FE wear model to compare the wear evolution between the models, which was then validated against wear measurements of the taper surface of explanted femoral heads. The fixation at the taper junction was found to be better for the smooth couplings. Over a 7 million load cycle analysis in-silico, the linear wear depth and the total material loss was around 3.2 and 1.4 times higher for the femoral heads mated with micro-grooved trunnions. It was therefore concluded that smooth taper and trunnion surfaces will provide better fixation at the taper junction and reduce the volumetric wear rates.