Advanced Mathematics and Numerical Modeling of IoT (Internet of Things)
View this Special IssueResearch Article  Open Access
Dongmin Seo, Hanmin Jung, WonKyung Sung, Dukyun Nam, "Development of the Korean Spine Database and Automatic Surface Mesh Intersection Algorithm for Constructing eSpine Simulator", Journal of Applied Mathematics, vol. 2014, Article ID 471756, 11 pages, 2014. https://doi.org/10.1155/2014/471756
Development of the Korean Spine Database and Automatic Surface Mesh Intersection Algorithm for Constructing eSpine Simulator
Abstract
By 2026, Korea is expected to surpass the UN’s definition of an aged society and reach the level of a superaged society. With an aging population come increased disorders involving the spine. To prevent unnecessary spinal surgery and support scientific diagnosis of spinal disease and systematic prediction of treatment outcomes, we have been developing eSpine, which is a computer simulation model of the human spine. In this paper, we present the Korean spine database and automatic surface mesh intersection algorithm to construct eSpine. To date, the Korean spine database has collected spine data from 77 cadavers and 298 patients. The spine data consists of 2D images from CT, MRI, or Xray, 3D shapes, geometry data, and property data. The volume and quality of the Korean spine database are now the world’s highest ones. In addition, our triangular surface mesh intersection algorithm automatically remeshes the spineimplant intersection model to make it valid for finite element analysis (FEA). This makes it possible to run the FEA using the spineimplant mesh model without any manual effort. Our database and surface mesh intersection algorithm will offer great value and utility in the diagnosis, treatment, and rehabilitation of patients suffering from spinal diseases.
1. Introduction
eHealth technologies supporting high efficiency and lowcost medical service based on ITBT convergence technology have grown in importance because medical expense is increasing and many people are asking for customerdriven medical service in today’s aging society. In particular, chronic diseases such as degenerative spinal diseases, high blood pressure, diabetes, and cancer cannot be treated adequately and require steady selfcare. Therefore, eHealth technologies are being widely used to treat chronic diseases widely. Nowadays, next generation chronic disease management technologies are being extensively researched based on mobile service, cloud computing service, social network service, big data analysis service, genome sequencing service, and computer simulation service [1–3]. By 2026, Korea is expected to surpass the UN’s definition of an aged society and reach the level of a “superaged society.” As a result, degenerative spinal diseases and related surgical procedures will increase exponentially. As of 2007, medical expenses incurred due to spinal surgery in Korea totaled 178.6 billion won/year, and the treatment duration reached 1.82 million days/year as shown in Figure 1. The resulting medical burden and economic loss are increasing at a rapid rate. Spinal diseases make everyday life of people impossible and impede economic activities, resulting in a compromised quality of life. Between 2002 and 2004, spinal surgery increased at a particularly high rate among the older demographic, by 68.2% among ages 60–69 and by 94.6% among those aged 70 and older. Among the leading causes of hospitalization for ages 65 and older, spinal diseases ranked number 2 with over 65,000 instances [4].
To prevent unnecessary spinal surgery resulting from overtreatment, systematic prediction of treatment outcomes is necessary, including scientific diagnosis, scientific effect analysis, and analysis of spinal rehabilitation exercises. Computer simulations have been utilized in biomechanical research for the past three decades. Today, advances in computer hardware and software are bringing continually increasing simulation accuracy. We have been developing eSpine which is composed of a spine database (DB) and a simulation environment for running a computer simulation model of human spines created by mathematically calculating images, geometries, and properties of human spines and will allow virtual testing without using a real human spine. For this, the high quality spine data and a surface mesh intersection algorithm are essential for realizing eSpine (The proposed intersection algorithm is applicable to a triangular surface mesh. In this paper, the surface mesh means the triangular surface mesh). In virtual surgery with an implant, surgeons or biomechanics researchers select the desired model from the spine DB and try to merge the spine model with an implant model as a virtual case. At this point, the two models of the spine and implant are automatically merged by means of the surface mesh merging functionality. With the merged model, they run the computer simulation and analyze the simulation result. To build the spine database (DB) for eSpine, we have produced and collected many images, geometry, and property data of spines from the Korean cadavers and patients with normal spine or degenerative spinal diseases [5, 6]. For automatic mesh merging functionality, we implement the surface mesh intersection algorithm used in the procedure of remeshing the spineimplant intersection model for finite element analysis (FEA). This makes it possible to run the FEA using the spineimplant mesh model without any manual effort.
In this paper, we present our spine database and surface mesh intersection algorithm in detail. The rest of the paper is organized as follows. Section 2 explores eSpine. Section 3 explains the Korean spine database obtained from cadavers and patients with degenerative spinal diseases. Section 4 describes our surface mesh intersection algorithm on the Korean spine data. Section 5 shows our experimental results. Finally, Section 6 presents the conclusion.
2. eSpine
eSpine is a computerrun simulation model created by mathematically modeling collected human spinal image data, which allows virtual testing without using a real spine. Figure 2 shows a comparison with a vehicle’s navigation system. For optimal and safe driving, a navigation system collects map information, models a map, and predicts a route. Similarly, for optimal treatment, eSpine collects spine images, models 3D spine models, and predicts virtual testing and results. The expected effect of eSpine is as follows:(i)acquisition of reliable, economical, advanced ITbased medical support technologies that can be used in the diagnosis and treatment of degenerative spinal diseases;(ii)strengthened market competitiveness for Korea’s medical equipment industry through the utilization of eSpine;(iii)reduction of medical expenses and improving the quality of life during old age by making available reliable, affordable ITbased medical technologies.
3. Korean Spine Database
3.1. Spine Sample Selection
There are distinct differences, such as facial features, skin color, and hair, between the human races. Differences in the organs of the body are also evident. Generally, spine data obtained from hospitals are only CT, MRI, or Xray, which focus on particular parts of the spine with degenerative diseases. The data are only available for diagnosing spinal diseases. Therefore, to construct eSpine and support the Korean spine research, we constructed the Korean spine data for degenerative spinal diseases. Figure 3 shows the process of the Korean spine database construction. The metadata and schema of our Korean spine database are accepted as a standard from Telecommunications Technology Association and are discussed in detail in [7]. To date, we have collected spine data from 77 cadavers and 298 patients. Among the leading causes of hospitalization for ages 65 and older, spinal diseases ranked number 2 next to cataracts. Therefore, most of the data were obtained from aged cadavers and patients over 50 but some were gathered from younger patients because even young people occasionally suffer from degenerative diseases of their vertebras.
3.2. Spine 2D Image and 3D Shape Data
Our spine 2D images consist of various images from Xray, CT, MRI, and BMD because different types of images are useful for diagnosing different types of diseases. These images are stored in the Digital Imaging and Communication in Medicine (DICOM) file format, which is a standard format for storing medical data. Some images like those from CT often comprise a series of images produced with small intervals and these images need to be stored and managed as a group in our database for efficient search and management. To accomplish this, we add series numbers to the end of each image file name while sharing the prefix of each file name and manage the file name prefix with start and end of series number as additional metadata. From a series of CT images, we can make a 3D shape model of the spine body by piling up image series in order and filling up small triangles to correct intervals between the images. This process was done with Maya, 3D animation software. Therefore, users can analyze the correlation between crosssectional CT images and 3D shape models.
3.3. Spine Geometry and Property Data
Spine geometry data are lengths and angles of key elements for representing shapes and features of the spine. We selected 481 spine geometry data items in [8–10] and measured those of our cadavers. The selected spine geometry data items are accepted as a standard from Telecommunications Technology Association and are discussed in detail in [11]. We did not measure the geometry data of a patient because we did not harvest the spine from a patient. The selected spine geometry data are utilized to analyze characteristics of the Korean spine. Figure 4 shows the selected elements in the spine geometry data related to cervical vertebras.
(a) Cervical vertebra 2
(b) Cervical vertebras 3–7
Figure 5 shows the selected elements in the spine geometry data related to thoracic vertebras.
Figure 6 shows the selected elements in the spine geometry data related to lumbar vertebras.
To obtain spine property data, we distribute a spine to cervical vertebras (C3–C7), thoracic vertebras (T1–T6, T7–T12), and lumbar vertebras (L1S1) and measured flexionextension, lateral bending, twist, and disc compression between bones. Table 1 shows the test conditions and Figure 7 shows the testing devices consisting of spine simulator of MTS and Liberty of Polhemus for measuring spine property data.

3.4. Degenerative Diseases on Spine
The human spine is divided into three parts: cervical, thoracic, and lumbar vertebras, starting from head. As humans age, they often suffer from various degenerative diseases of the cervical and lumbar vertebras or intervertebral discs, while they rarely suffer from diseases of the thoracic vertebras. Therefore, we decided to target the following degenerative diseases frequently occurring at the cervical or lumbar vertebras of older people and further decided to classify these diseases into four or five grades (i.e., grade 0 to grade 3 or grade 4) through consulting with several specialists in the human spine. The larger the grade number is, the more severe the disease is. Grade 0 indicates normal and grades 1, 2, 3, and 4 correspond to mild, moderate, severe, and very severe grade of disease, respectively. Below is a brief description of the selected degenerative diseases on the spine.(i)Osteophyte. This occurs at both cervical and lumbar vertebras and can be further classified into anterior osteophyte and posterior osteophyte. This can be generally diagnosed by Xray and CT.(ii)Ligament calcification (or ossification). There are three kinds of ligaments around cervical vertebras: anterior and posterior longitudinal ligaments and ligamentum nuchae. Thus, this disease can be further divided into three subdiseases: ossification of the anterior longitudinal ligament (OALL), ossification of the posterior longitudinal ligament (OPLL), and ossification of ligamentum nuchae (OLN). This can generally be diagnosed by Xray and CT.(iii)Endplate sclerosis. This occurs at both cervical and lumbar endplates and can generally be diagnosed by Xray, CT, and MRI.(iv)Disc height reduction. This occurs at both cervical and lumbar intervertebral discs and can generally be diagnosed by Xray and CT.(v)Disc herniation. This occurs at both cervical and lumbar intervertebral discs and can generally be diagnosed by Xray, CT, and MRI. This is usually classified into five grades: normal, bulging, protrusion, extrusion, and sequestration.(vi)Disc degeneration. This occurs at both cervical and lumbar intervertebral discs and can generally be diagnosed by Xray, CT, and MRI. This is also classified into five grades.(vii)Facet joint degeneration. This occurs at both cervical and lumbar facet joints and can generally be diagnosed by Xray and CT.(viii)Spondylolisthesis. This occurs at both cervical and lumbar vertebras, but most commonly at lumbar vertebras, and can generally be diagnosed by Xray.(ix)Osteoporosis. This occurs at any vertebras and can be generally diagnosed by Xray and BMD (bone mineral density).
3.5. Statistics about the Korean Spine Data
To date, we have collected the Korean spine data from 77 cadavers and 298 patients with normal spine or degenerative spinal diseases. The detailed statistics on the collected data is shown in Table 2.

4. Automatic Surface Mesh Intersection Algorithm
To support automatic mesh merging functionality, we implement the mesh intersection algorithm. We applied a tracking algorithm in order to rapidly and accurately explore an intersection. The tracking algorithm finds the intersection regions along the intersection curve from a valid tracking point [12, 13]. After finding an initial intersection point, the algorithm starts at the initial intersection point and creates an intersection curve along the direction to new intersection points. While this happens, the algorithm also searches for intersection regions. We use the existing data structure in [12] and add one factor which can distinguish between spines and implant meshes.
4.1. Finding Intersection Points
4.1.1. Searching for Intersection Regions
A plane equation is derived from the outer product of the points of mesh elements. An intersection point is calculated using the topological relation between two intersection planes. After that, we derive the angle from the inner product of three points of a triangle element. If the sum of the derived angles is 2, the intersection point is inside the triangle element. The intersection of planes on threedimensional space has various cases. When the tracking algorithm searches for an intersection, we determine the intersection from an interrelation of a line and a side and there are cases of intersection between intersected triangles in Figure 8.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
4.1.2. Tolerance
The intersection cannot be mathematically determined because the numerical calculation of a computer does not work on consecutive space. We have to define a tolerance to calculate an intersection. Suppose that the gap has the difference as in Figure 9. If , two triangles are not intersected. Otherwise, they are intersected. When two meshes intersect, the shape of intersected meshes changes. This is because the intersection points are derived from the tolerance by moving existing shapes. The shape of an implant in a spineimplant insertion model should not change using the tolerance. Algorithm 1 defines an intersection generation procedure with the tolerance. In case of an intersection point between edge and face , the point’s position depends on which one is an implant mesh between edge and face . If is an element of an implant mesh, the intersection point is the start or end point of an edge . On the other hand, if is an element of an implant mesh, the intersection point is generated by intersecting and the extension of . New mesh points are located on the edge or face of implant mesh because our unique constraint is that implant mesh does not change as much as possible. The tolerance is automatically set up according to the size of a model or by user’s configuration.

4.2. Tracing Algorithm
In the model for the volume finite element, an intersection curve is always closed, so tracking an intersection can start from an initial intersection point to the intersection progress direction. After finding an initial intersection point, the intersection points continuously search along the intersection progressing direction [12–14]. An intersection curve is continuously generated by inputting the continuously generated intersection points to an intersection curve . Using the intersection curve data, we retrieve intersection regions and generate new mesh of the intersection region.
4.3. Meshing
To generate a spineimplant finite element intersection model, the generated mesh on an intersection region should satisfy three conditions.(i)The shape of the spine model can be changed but the shape of the implant cannot be changed. This is because this spineimplant structure analysis simulates the process of inserting a strong durable implant into a relatively weak durable spine. If we allow the shape change of the implant, nonintended stress concentration can occur so the implant shape does not allow the change.(ii)We do not generate more meshes than the user needs. The size and the shape of a mesh work are important factors. A small size mesh is densely formed in a complicated shape, but if the number of meshes is blindly large, we waste more analysis resources than necessary.(iii)A mesh should have the right quality for the finite element analysis. The intersection point, which can only be created by an intersection search, has an unbalanced gap. A mesh with the intersection points is not enough for finite element analysis, so we need to add fixed points along the intersection curve to create a mesh for finite element analysis.
When generating the intersected mesh model, we have to satisfy the above three requirements. First of all, the points of the intersection curves should be reorganized. Algorithm 2 reorganizes the points of the intersection curve. Figure 10(b) shows the points to define the shape of an implant among the points of an intersection curve. is defined as the minimum gap of all other points. The unequally distributed points on the crosscurve are reorganized equally based on the standard.

(a)
(b)
(c)
(d)
The intersection points created on the points or the lines of an implant element are presented as follows:
The algorithm automatically creates a new element network to keep the shape of the original model by using Delaunay triangulation in Figure 10(d). The element created by Delaunay triangulation may be not suitable for finite element analysis. Therefore, the algorithm reorganizes the triangulation element in a spine model through a remesh process. In addition, we remesh the intersection element with neighboring elements to prevent a sharp form. The tracking algorithm searches for the continuous two points and to satisfy this condition. In Figure 11, we generate the points with a uniformed gap. The algorithm stores the points on the new generated curve and constructs the line , between newly generated points. The new organized line is sequentially stored as a hard edge, . The hard edge is fixed and a base form on mesh or remesh processes.
(a)
(b)
4.4. Inserting Implant
After processing triangulation and remesh, we remove the existing elements of intersection to insert the reorganized elements. The algorithm intersects two finite element models by creating the inserted part inside a target element model. We search for the intersection part with the outer product of surface elements along the intersection curve. In Figure 12, the algorithm removes the region where an implant is inserted in the spine model after searching for the part of the implant inserted in the spine model. It finishes an automatic intersection processing of an implant and a spine by inserting the implant according to the direction to the spine model.
(a)
(b)
5. Performance Evaluation
Figure 13 shows an accurate validation of motion properties with the extension and flexion on L5S1 on a spine. The line with the symbol is the values of motion properties obtained by [15]. The solid line is obtained by [16]. The line with the squared symbol is obtained by our database. The motion properties from [15] and our database are constructed through experiments. The motion properties from [16] are constructed through computational models. In conclusion, we can confirm that our spine data is accurate because the shape of each line is similar.
A mesh is closely related to an analysis result and the intersection processing in the intersection regions influences the existing shape. Therefore, in most cases, a mesh is manually created to guarantee the mesh quality. We show that the proposed algorithm creates an appropriate mesh for a structural mesh. The algorithm creates a mesh by automatic intersecting of spine and implant models. As the spine model, we use a threelayered spine model consisting of a vertebral arch, an outer vertebral body, and an inner vertebral body. To evaluate the usefulness of the automatic intersection algorithm, we prepare three automatic intersection models with different sizes of implants, as mesh sizes 1.0, 1.5, and 2.0 in Figure 14.
(a)
(b)
(c)
In Figure 15, we set up the analysis conditions: 10,000 load to the direction on top of the spine with fixed 6 degrees of freedom (DOF). We performed the analysis using an ABAQUS solver.
Figure 16 shows the common patterns even though there is little derivation according to the element size. We cannot find the abnormal results due to intersection regions. Therefore, we expect that the proposed automatic intersection algorithm can be used without human intervention.
6. Conclusion
We collected various Korean spine data items from 77 cadavers and 298 patients with normal spine or degenerative spinal diseases to provide a wealth of information related to spine to medical students, physicians, and biomedical engineers. We also propose the automatic surface mesh intersection algorithm for spine and implant models. Our algorithm automatically remeshes the spineimplant intersection model to make it valid for finite element analysis (FEA). A spineimplant intersection model is manually created so far. The automatic intersection procedure using the proposed intersection algorithm reduces the manual labor time for spineimplant model. Therefore, this makes it possible to run the FEA using the spineimplant mesh model without any manual effort. We show the validation of intersection mesh quality in the simulation. In the near future, we plan to define criteria to check the quality of a spineimplant intersection mesh and perform the numerical analysis based on the defined criteria.
We will offer our spine data and surface mesh intersection algorithm to many researchers and doctors to foster spine research development. In conclusion, our spine data and surface mesh intersection algorithm will be used to realize reliable, economical, and advanced ITbased medical support technologies that can be used in the diagnosis and treatment of degenerative spinal diseases. Furthermore, our technical skills will be used to vitalize the related research fields and industries by providing the developed human spinal information database, model, and virtual simulator to relevant researchers.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by 2014 National Agenda Project (NAP) funded by Korea Research Council of Fundamental Science & Technology (NAP092).
References
 N. Howard and E. Cambria, “Intention awareness: improving upon situation awareness in humancentric environments,” HumanCentric Computing and Information Sciences, vol. 3, no. 9, pp. 2–17, 2013. View at: Google Scholar
 P. Gargiulo, T. Helgason, P. Ingvarsson, W. Mayr, H. Kern, and U. Carraro, “Medical image analysis and 3d modeling to quantify changes and functional restoration in denervated muscle undergoing electrical stimulation treatment,” HumanCentric Computing and Information Sciences, vol. 2, no. 10, pp. 2–11, 2012. View at: Publisher Site  Google Scholar
 M. Brahami, B. Atmani, and N. Matta, “Dynamic knowledge mapping guided by data mining: application on Healthcare,” Journal of Information Processing Systems, vol. 9, no. 1, pp. 1–30, 2013. View at: Publisher Site  Google Scholar
 National Health Insurance, “2009 Major surgery statistics,” Press Release, vol. 12, no. 7, pp. 1–15, 2010. View at: Google Scholar
 “eSpine project,” http://www.espine.kr/. View at: Google Scholar
 D. M. Seo, S. W. Lee, S. B. Lee, S. H. Lee, H. M. Jung, and W. K. Sung, “Implementation of spine database system on Korean patient with degenerative spinal diseases,” in Proceedings of International Smart Media Application (SMS '12), pp. 1–3, 2012. View at: Google Scholar
 S. B. Lee and S. H. Lee, “Metadata for managing the human body data,” Telecommunications Technology Association, TTAStandard (TTAK.KO10.0723), 2013. View at: Google Scholar
 M. M. Panjabi, J. Duranceau, V. Goel, T. Oxland, and K. Takata, “Cervical human vertebrae: quantitative threedimensional anatomy of the middle and lower regions,” Spine, vol. 16, no. 8, pp. 861–869, 1991. View at: Publisher Site  Google Scholar
 M. M. Panjabi, K. Takata, V. Goel et al., “Thoracic human vertebrae: quantitative threedimensional anatomy,” Spine, vol. 16, no. 8, pp. 888–901, 1991. View at: Publisher Site  Google Scholar
 M. M. Panjabi, V. Goel, T. Oxland et al., “Human lumbar vertebrae: quantitative threedimensional anatomy,” Spine, vol. 17, no. 3, pp. 299–306, 1992. View at: Publisher Site  Google Scholar
 S. B. Lee and S. H. Lee, “The measuring elements for representing the human spine shape,” Telecommunications Technology Association, TTAStandard (TTAK.KO10.0724), 2013. View at: Google Scholar
 A. H. Elsheikh and M. Elsheikh, “A reliable triangular mesh intersection algorithm and its application in geological modelling,” Engineering with Computers, vol. 30, no. 1, pp. 143–157, 2014. View at: Publisher Site  Google Scholar
 S. H. Lo and W. X. Wang, “A fast robust algorithm for the intersection of triangulated surfaces,” Engineering with Computers, vol. 20, no. 1, pp. 11–21, 2004. View at: Publisher Site  Google Scholar
 S. H. Lo, “Automatic mesh generation over intersecting surfaces,” Numerical Methods in Engineering, vol. 38, no. 6, pp. 943–954, 1995. View at: Publisher Site  Google Scholar
 M. M. Panjabi, T. R. Oxland, I. Yamamoto, and J. J. Crisco, “Mechanical behavior of the human lumbar and lumbosacral spine as shown by threedimensional loaddisplacement curves,” Journal of Bone and Joint Surgery, vol. 76, no. 3, pp. 413–424, 1994. View at: Google Scholar
 V. Moramarco, A. P. del Palomar, C. Pappalettere, and M. Doblaré, “An accurate validation of a computational model of a human lumbosacral segment,” Journal of Biomechanics, vol. 43, no. 2, pp. 334–342, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Dongmin Seo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.