PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 11048.
Published online 2017 September 8. doi:  10.1038/s41598-017-08579-6
PMCID: PMC5591222

Freeform fabrication of tissue-simulating phantom for potential use of surgical planning in conjoined twins separation surgery

Abstract

Preoperative assessment of tissue anatomy and accurate surgical planning is crucial in conjoined twin separation surgery. We developed a new method that combines three-dimensional (3D) printing, assembling, and casting to produce anatomic models of high fidelity for surgical planning. The related anatomic features of the conjoined twins were captured by computed tomography (CT), classified as five organ groups, and reconstructed as five computer models. Among these organ groups, the skeleton was produced by fused deposition modeling (FDM) using acrylonitrile-butadiene-styrene. For the other four organ groups, shell molds were prepared by FDM and cast with silica gel to simulate soft tissues, with contrast enhancement pigments added to simulate different CT and visual contrasts. The produced models were assembled, positioned firmly within a 3D printed shell mold simulating the skin boundary, and cast with transparent silica gel. The produced phantom was subject to further CT scan in comparison with that of the patient data for fidelity evaluation. Further data analysis showed that the produced model reassembled the geometric features of the original CT data with an overall mean deviation of less than 2 mm, indicating the clinical potential to use this method for surgical planning in conjoined twin separation surgery.

Introduction

Conjoined twinning is a rare congenital malformation15. Although most of the conjoined twins succumb in utero or are stillborn at birth, approximately 1: 30,000–1: 200,000 are born alive. These conjoined twins often display life-threatening conditions and severe disabilities, posing significant challenges to both healthcare providers and patients6. Surgical separation is possible in ~1 out of 650–900 conjoined twins born alive1, 2. Success of separation surgery is affected by many factors including the clinicians experience and an accurate understanding of the anatomy of the affected organs7, 8. In order to reduce surgical uncertainty and increase the success rate, it is important to acquire detailed anatomic information and carry out accurate surgical planning prior to separation surgery.

Clinically, various medical imaging techniques and simple medical models are used to support surgical planning before surgical separation of conjoined twins. Typically, anatomic details of conjoined twins are obtained by conventional medical imaging modalities810. For example, ultrasonography is used to assess the brain11; fluoroscopy is used to obtain the anatomic features of anal fistula10; and high-resolution computed tomography (CT) and magnetic resonance imaging (MRI) are used to reveal the geometric and positional details of the internal organs12, 13. These imaging tools have been used to support successful surgical planning2, 7, 14, 15.

In addition to these imaging tools, conjoined twin models made of foam core and silicone rubber have also been used to guide separation surgery16. However, in most cases, the phantoms fabricated using traditional methods are restricted by inability to simulate the multidimensionality of the body, with no appropriate fabrication method for reproducing three-dimensional (3D) shape.

More recently, 3D printing can be used to produce the models of conjoined twins for surgical planning. 3D printing represents a family of additive manufacturing processes that are capable of producing complex geometric structures through precise layer-by-layer deposition of materials17, 18. Owing to the flexibility for producing freeform geometric features, 3D printing has been widely used in many medical applications, such as medical education, dental customization, skull defect repair and cardiac bypass surgical planning1924.

In this paper, we present a new method that combines 3D printing, assembling, and casting to produce anatomic models of high fidelity for surgical planning of conjoined twin separation. The related anatomic features of the conjoined twins were captured by CT and classified into the following five groups: skeleton group, spinal nerve group, colon group, kidney-bladder group, and other tissue group. The geometric and positional features of these groups were reconstructed from a 3D computer model 25. The casting and the assembling processes were designed based on the reconstructed computer models of these organ groups. With acrylonitrile-butadiene-styrene (ABS), silica gel, contrast agents and plasticine, individual internal phantoms were fabricated using a freeform method that combines 3D printing and casting. The produced phantoms were then assembled, positioned firmly within a 3D printed shell mold system simulating the skin boundary, and casted with transparent silica gel. The produced conjoined twin phantom was further scanned by CT for fidelity evaluation. Subsequent data analysis showed that the produced phantom had an overall mean deviation of less than 2 mm in comparison with the original models, closely reassembling the anatomic features of the conjoined twins. These results support the technical potential of simulating structural, optical and CT properties of relevance in conjoined twins.

Results

A conjoined twin model produced by casting individually fabricated phantoms

Fabrication of the conjoined twin model was based on casting the phantoms of individual organ groups (Fig. 1a), assembling these phantoms in an assembly (Fig. 1b), and placing the assembly in a shell mold for further casting (Fig. 1c). To fabricate the phantom of each organ group, a shell mold was first produced by the FDM process using the CT data. In order to simulate the visual and CT contrasts of the actual tissue components, different contrast agents were added in silica gel. To produce organ groups 2 and 6, nanoscale tantalum at concentrations of 0.024 g/ml and 0.04 g/ml was added in silica gel respectively to simulate different CT contrasts. To produce organ groups 3 and 5, iomeprol at concentrations of 0.02 ml/ml and 0.04 ml/ml was added respectively. To produce organ group 4, green plasticine was used. To produce organ group 1, the hard skeleton was directly printed by the FDM process using acrylonitrile butadiene styrene (ABS). For the casting of the soft tissue phantoms, the inner surface of the shell molds were costed with a thin layer of release agent to facilitate demolding. The two parts of the shell mold were sealed, filled with silica gel, and placed in a sonication chamber for 5 hrs to prevent precipitation of contrast agents. After 72 hrs, the fully cured phantoms were demolded and assembled. To ensure the positioning accuracy of the assembly, multiple positioning holes and columns were designed on the phantoms (Fig. 1b–g). To produce the final conjoined twin model, a shell mold as shown in Fig. 1f was printed in advance. The thickness of the shell mold was 2 mm to ensure structural strength and stability. As the organ group assembly was placed inside the shell mold, the positioning holes on the assembly and on the shell molds were aligned in order to ensure the positioning accuracy of the organ groups. Transparent silica gel was mixed with the curing agent at a ratio of 10:1 (V/V) and poured into the shell mold to produce the green part. After the green part was cured for 120 hrs at room temperature, it was demolded. The redundant fixtures, the supporting surfaces (both top and bottom), the remaining pieces of the shell mold, and the protective layer were subsequently removed to obtain the final conjoined twin model.

Figure 1
Illustration of the fabrication process for a conjoined twin model: (a) Fabrication of a soft tissue phantom (using the kidney-bladder group as an example). From left to right: computer aided design (CAD) model of the kidney-bladder group, shell mold ...

Simulating visual and CT contrasts of individual organs in a high fidelity fabricated conjoined twin model

The produced conjoined twin model enable visual delineation of different organ groups, and resemble CT contrasts of biologic tissue, as shown in Fig. 2a,b. The model was subsequently examined by a 64-slice CT (Light Speed VCT, GE Healthcare) at the First Affiliated Hospital of Anhui Medical University. Figure 2b shows a cross-sectional CT image of the conjoined twin model. According to the figure, the averaged CT intensities for individual organ groups are −256, 710, 350, −600, 590, 1050, and 85 HU, respectively. Through both visualization and CT scan, the internal organ groups can be clearly distinguished, indicating the clinical utility of using such a model to guide the surgical planning.

Figure 2
Visual, CT contrasts and fidelity maps of the produced conjoined twin model: (a) Isometric view of the conjoined twin model where the internal organs can be clearly differentiated. (b) The CT scan image of the conjoined twin model where different organ ...

To ensure accurate guidance of the surgical procedure, the produced conjoined twin model should resemble the patient anatomy with fidelity. We have defined a quantitative method for fidelity evaluation based on an iterative closest point (ICP) algorithm. Figure 2d shows the fidelity maps of individual organ group phantoms in comparison with their original CT data, in which the deviation chromatography on each fidelity map can help to analyze error at the exact location. Considering that both the fidelity of individual organ group phantoms and the relative positions of these phantoms in the assembly contribute to the overall fidelity of the conjoined twin model, we have also calculated the fidelity of individual phantoms within the assembly. Table 1 lists the statistics of fidelity maps (Fig. 2d) of individual phantoms before and after assembly. No fidelity was calculated for the spinal nerve group as it was hand-made. As shown from the Table 1, the mean fidelity error of any individual phantom was less than 0.5 mm before assembly and less than 1.5 mm after assembly.

Table 1
Errors of the individual phantoms before and after assembling.

Discussion

Separation surgery for conjoined twins is always challenging because the sharing of organs in each of them is unique8. Regardless of the low survival rate of conjoined twins, many of those born alive will require separation surgery in order to live a normal life or just to survive. The feasibility of separation surgery after birth needs to be assessed because the distribution of organs is entirely different in each case26. Various imaging techniques provide abundant anatomic information conducive to the assessment of conjoined twins, and the anatomical information obtained from the imaging is of great help in the planning the operative pathway8, 14. However, these imaging tools are not interactive and cannot reflect complex 3D information of the tissues. Abdominal ultrasound, plain radiography, and 1D fluoroscopic examinations may provide minimal information. By contrast, advanced computed tomography (CT) and Magnetic Resonance Imaging (MRI) scanning can provide detailed 3D information, but resolution limitation and virtualization cannot provide maximal information10. This paper reports a freeform fabrication method for fabricating tissue-simulating phantoms by combining 3D printing and casting. The phantoms simulate not only the structural features of conjoined twins but also CT contrasts between different organs. Moreover, the phantom also allows clinicians to visually distinguish detailed cross information about the organs. Most notably, an ICP algorithm-based fidelity calculation algorithm for 3D shape has been developed. The algorithm offers an understanding of the phantom fabrication quality in terms of shape and position, and represents a possible direction in 3D shape similarity measurement for other applications. The comprehensive value of the conjoined twin model depends on the fabrication quality and assembling fidelity of the individual parts. Compared with the replicable fabrication of phantoms with the same contrast difference, the fidelity of the conjoined twin model is affected by the independent fabrication processes for producing individual organ group phantoms with distinct contrast differences. The primary factors contributing to the fidelity of the final conjoined twin model include the FDM fabrication of the protective layer, the FDM fabrication of the shell molds, and the assembly of individual organ group phantoms. The secondary factors contributing to the overall fidelity include the reconstruction of the solid model, the accuracy of CT scan, and the geometric deformation of silica gel during the curing process. In order to facilitate easy assembly of individual phantoms with high fidelity, the design of the positioning structure should consider the following requirements: (1) the positioning columns must remain perpendicular to the contact surface of the phantoms; (2) three positioning columns are required for each phantom; (3) the length of the positioning columns should be sufficient in order to support secure fitting with the matching holes in the neighboring phantoms; and (4) the center of each positioning hole should coincide with that of the corresponding positioning column. To evaluate the fidelity error caused by the modeling algorithm, we used different people to generate the CAD model of the phantoms and the resultant maximum arithmetic average error in repetitive model reconstruction was 0.35 mm with a deviation of less than 0.4 mm. In addition, the pixel/slice thickness of the CT was 0.825/0.5 mm, and the diameter of the 3D printing nozzle of the FDM printer was 0.25 mm. Further systematic analysis is needed to understand the error sources and their contributions to the fidelity of the final model.

Conclusion

We have demonstrated the feasibility of producing a high fidelity conjoined twin model with both visual and CT contrasts by combining 3D printing and casting. In the tissue classification procedure, the related organs were divided into five organ groups, with the CAD models reconstructed based on the CT scanning data. Among the five organ groups, the skeleton group was produced by the FDM printer using ABS. The rest four groups of soft tissue were produced by casting silica gel in FDM-produced shell molds. The produced organ group phantoms were assembled, positioned firmly within a shell mold, and cast with transparent silica gel to produce the conjoined twin model. In order to introduce visual and CT contrasts for different organ groups in the model, we have explored a variety of available contrast agents and added them in silica gel for different organ groups. We have also carried out the fidelity analysis for the produced organ group phantoms in comparison with the original CT data. Our study showed that the produced model reassembled the geometric features of the patient with an overall mean deviation of less than 2 mm, indicating the clinical potential to use this method for surgical planning in conjoined twin separation surgery. Consequently, the produced model could be applied in acquiring organ distribution information of the conjoined twin. Whereas the preliminary investigations conducted in this study focused on macro structures in this biological system, phantoms with more complex microstructures can be investigated in future studies.

Materials and Methods

Methodology for phantom design and fabrication

In the case of conjoined twin separation surgery, it is likely that improved understanding of organ distribution leads to higher chances of success. Therefore, it is of great significance to fabricate multifunctional phantoms to help surgeons better understand anatomic information.

The conjoined twin model was fabricated following a cascade of processes as illustrated in Fig. 3. Considering that the routinely used CT scan has a positional accuracy as high as 0.3 mm, the CT scan data were used as the basis for assessing the organ distribution of the conjoined parts. Taking into account the differences in surgical requirements and the internal organs of different interests involved in surgical procedures, organs should be simulated at different accuracy and spatial resolution levels. In this specific surgical case, we divided the internal organs into the following five groups: skeleton, spinal nerve, colon, kidney-bladder, and other tissue. Considering the possibility for future surgical simulation, these organ groups were further divided into hard and soft types based on their mechanical characteristics. The hard phantoms were directly printed by the FDM process, while the soft phantoms were produced by casting silica gel in a 3D printed shell mold based on the CT data. The individual organ group phantoms were then assembled and casted in a shell mold to produce the final conjoined twin model.

Figure 3
Schematic illustration of the processes involved in fabrication of a conjoined twin model.

Materials

The contrast enhancement pigments were added in transparent silica gel to simulate both visual and CT contrasts of individual organ groups. For visual differentiation of different organ groups, dyes of different colors were added. To simulate the contrast differences between individual organ groups, positive and negative CT contrast agents at different concentration levels were added. The CT contrast is defined as the X-ray attenuation difference in tissue as Hounsfield units (HU) in CT scanning slices. It works with the visual contrast to provide quantitative boundary information for tissue localization and differentiation in a surgical procedure. On the other hand, the conjoined twin model with organ groups of different contrasts enables subsequently quantitative evaluation of the structural fidelity in comparison with the original CT data. To facilitate quantitative guidance for separation surgery, the phantoms should simulate geometric features of individual organ groups with strong visual contrasts and large differences in the HU values27, 28. To facilitate easy visualization of the internal organ structures, the organ groups should be embedded in a transparent phantom material.

Transparent silica gel, contrast agents with different colors, ABS and plasticine were chosen as the construction materials for the organ group phantoms. The transparent silica gel (Nanjing Hua Cheng Chemical Co., Ltd, Nanjing, China) is transparent, non-toxic, chemically stable, and with mechanical properties tunable by changing the mixing ratio between the curing agent and the silica gel, making it an ideal base material for constructing tissue-simulating phantoms29, 30. To simulate both visual and CT contrasts of the individual organ groups, contrast agents with different colors and different Hu levels were added in the silica gel27, 28. Figure 4a shows the silica gel samples cast in small beakers after mixing with nanoscale aluminum, nanoscale tantalum, nanoscale silver, and nanoscale titanium respectively at concentrations of 0.008, 0.016, 0.024, 0.032, and 0.04 g/ml. Similarly, silica gel samples with iomeprol concentrations of 0.02, 0.04, 0.06, 0.08, and 0.10 ml/ml were also produced. After casting, the samples were placed in a sonication chamber for curing in order to prevent precipitation of contrast enhancing pigments. The produced phantoms were subject to a 64-slice CT scan (Light Speed VCT, GE Healthcare) for Hu contrast evaluation, as plotted in Fig. 4b. The casting experiments were triplicated for statistic analysis. According to the figure, we speculated that the CT contrast of the phantom was linearly correlated with the concentration level of the contrast agent within a certain range. In order to achieve simultaneous visual and CT contrasts, we simulated the colon organ group by mixing silica gel with 0.024 and 0.04 g/ml nanoscale tantalum (Shanghai Shui Tian Material Technology Co., Ltd., China); simulating the kidney-bladder organ group by mixing silica gel with 0.02 and 0.04 ml/ml iomeprol (Bracco Imaging Italia SRL, Italia)31, 32; simulating the solid skeleton group using light yellow ABS with HU CT value of −333; and simulating the spinal nerve group using a green plasticine with a HU CT value of −75025.

Figure 4
Silica gel phantoms with different types and concentrations of contrast enhancement pigments to simulate different visual and CT contrasts: (a) Silica gel mixed with a variety of CT contrast agents at different concentration levels; and (b) Quantitative ...

Design and 3D printing of the shell molds for organ groups

The models (Fig. 5a) of the organ group were reconstructed based on the conjoined twins computed tomography (CT) data with the help of the Mimics software package (Mimics 16; Materialise Dental)33. They were subsequently exported into editable initial graphics exchange specification format files23, 34. Unlike the hard phantoms were printed directly, soft phantoms were casted with the aid of the molds designed using Unigraphics NX package (UG 9.0, Unigraphic Solutions Inc., USA) and Geomagic software package (GeoMagic Studio 2012; GeoMagic USA). Just as Fig. 5b shows, by employing operations such as shell construction etc., molds for inter soft tissue phantom casting were designed respectively. Each of them was divided into 2–4 parts for later demolding35. However, the successful fabrication of individual phantoms is not enough, the accuracy of the conjoined twin model also depends on the assembly precision, an appropriate positioning project ensuring assembly precision were thereby expected. Mathematical theory demonstrates that three points in different planes position a rigid object in three dimensional (3D) space. Figure 5d–f show that three positioning columns or holes were designed for each inter individual phantom with the base of hard skeletons. Particularly, the positioning structures were designed to exactly reach the surface of the connected organ, and groove-like structures (Fig. 5f) were helpful in assembling work. Moreover, 2 mm-thick molds (Fig. 5h) for tissue phantom casting were finally designed, and positioning holes (as shown in Fig. 5g) on top and bottom surfaces were coaxial with that on skeleton phantoms, providing precise position of inter individual phantoms inside tissue phantom. All the molds were printed by FDM process. However, limited by the accuracy of the FDM printer, internal surfaces of molds were grainy, which would result in the surface defect of casted phantoms, hence a 0.21 mm-thick smooth protective silica gel layer was fabricated on internal surface of each mold. As a result, all the molds were printed 1.01 times of the original size with a 3D printer (Dimension 1200es, Stratasys, US) to balance the shrinking caused by the protective layer.

Figure 5
Illustration of the reconstruction of organ computer aided design (CAD) models, design of positioning project for phantom assembling, and fabrication process of molds for soft phantom casting: (a) Reconstructed organ CAD models of the conjoined twins. ...

Fidelity calculation based on iterative closest point (ICP) algorithm

To facilitate accurate simulation of tissue anatomy by the 3D printed model, it is critically important to develop a quantitative method for fidelity assessment. Based on the Merriam-Webster Dictionary, fidelity refers to the degree to which something matches or copies something else14. For many years, fidelity research has been a topic for better understanding of the human visual perception system36, 37. Shape fidelity-based retrieval of 3D data has been a focus of research in various disciplines, such as computer vision34, mechanical engineering38, artifact searching39, molecular biology35, and chemistry33. The shape fidelity-based retrieval method has been used in numerous applications for model similarity calculation. However, none of these applications are similar to that of the conjoined twin model. We propose a method for fidelity assessment that involves similarity analysis for both the shape and the position of the produced models.

Theoretically, the term fidelity as described in this paper is defined as the averaged error between the reference vertices in the original model and the duplicated model after optimal alignment. Practically, it is difficult to obtain optimal alignment between two models, not only because the corresponding models are unordered and without a unified world coordinate, but also because both models in this research were in STL format with numerous vertices. To calculate the fidelity parameters and the corresponding color map fo an organ group phantom, we developed a procedure as shown by the flowchart in Fig. 6a. This procedure was further validated by a Geomagic software package (GeoMagic Studio 2012; GeoMagic USA)40.

Figure 6
Definition and verification of the fidelity assessment method based on an ICP algorithm: (a) Flowchart displaying the processes for fidelity calculation. (b) Standard cube model. (ce) Random deformations that applied random movement, random rotation, ...

Before fidelity is calculated, the models were first aligned in a unified coordinate system using an ICP algorithm41, 42. As shown in the flowchart (Fig. 6a), the ICP alignment is implemented by calculating the best suitable rotation matrix and translation matrix. Takes the k-th alignment for example. τ is the tolerance. Pk = {pki} and Zk = {zki} (i = 1, 2  N) are stochastic point sets of corresponding models, where N are the number of the stochastic point, and point pki corresponds to point zki. Rkj and Tkj(j = 1, 2) are rotation matrix and translation matrix respectively. Fk(Rkj, Tkj) (j = 1, 2) is standard least-squares distance between point sets.

According to the best fitting principle, there should be an ideal rotation matrix Rkj and translation matrix Tkj, so that the set RkjPk + Tkj aligns best with Zk 42, and the least-squares distance function (1) is minimized.

Fk(Rkj,Tkj)=i=1N|pkiRkj+Tkjzki|2
1

Calculation of ideal rotation matrix Rkj and translation matrix Tkj begins from the initialization of Rk1 and Tk1. Lets Cp in (2) be the center of the point set Pk, and Cz in (3) to be the center of the point set ZK.

Cp=i=1Npki
2

Cz=i=1Nzki
3

As for the nest rotation matrix Rk2 and translation matrix Tk2, Tk2 = Cz  Rk2Cp.

Lets P˜ki=pkiCp, Z˜ki=zkiCz. To compute the next rotation Rk2, the problem is equivalently simplified into maximizing the one-variable objective function (4).

FK(R)=i=1NP˜kiR,Z˜ki
4

Let H=i=1NP˜ki(Zki)T, which is a 3 × 3 matrix. According to both orthogonal property of rotation matrix and singular value decomposition principle41, 43, there should be 3 × 3 orthogonal matrices, for example U and V, and 3 × 3 diagonal matrix Λ with non-negative elements exist such that H = UΛV. According to the singular value decomposition (SVD) algorithm in the least-squares method41, to maximize the function FK(R), the next rotation matrices are given by (5).

Rk2={VUTdet(VUT)=1V(100010001)UTdet(VUT)=1
5

The rotation matrices Rk1, Rk1 and the translation matrices Tk1, Tk2 that satisfy both function (6) and function (7) are calculated and saved.

Fk(Rk2,Tk2)<Fk(Rk1,Tk1)
6

Fk(Rk1,Tk1)Fk(Rk2,Tk2)<τ
7

Then, the new stochastic point sets Pk+1 and Zk+1 are chosen, and the corresponding least distance Fk+1(Rk+12, Tk+12) is calculated. The calculation iterates until both function (8) and function (9) are satisfied. Finally, the errors are the calculated and the fidelity maps are generated.

Fk+1(Rk+11,Tk+11)<Fk(Rk1,Tk1)
8

Fk(Rk1,Tk1)Fk+1(Rk+11,Tk+11)<τ
9

Models of various complexities and their random deformations are used to validate the feasibility of the aforementioned fidelity calculation algorithm. To he shape deformed models were obtained by the free-form deformation (FFD) technique that deforms an object enclosed within a cube by moving the control points on its edges44. Thedeformation of the model is defined by a matrix random (M, R, V), where three independent variables M, R, and V represent translation, rotation, and distortion, respectively. The default value of M (or R) is 0 and it changes to 1 as the deformation involves random translation (or rotation). The value of V (ranging from 0 to 1) is defined as proportion of deformed control points in the FFD model. For validation of the fidelity calculation algorithm in solid models with moderate deformation, the random motion of the control points in the FFD model is restricted within 10% of the longest axis. As an example, fidelity analysis is carried out in a standard cube model (Fig. (Fig.6b),6b), and the resultant fidelity maps after translation, rotation, and distortion operations are shown in Figs. Figs.6c,6c, ,6d,6d, and and6e,6e, respectively. According to the fidelity maps, neither translation nor rotation affect the fidelity of the solid model, while distortion induces the distributed fideltiy changes. Table Table22 lists the fidelity analysis results for 3D models of various geometrical feastures and deformations. The results in Table 2 indicate that neither random translation nor rotation induce any fidelty issue in a solid model, while distortion induces fidelity errors proportional to the positional deviation of the control points. It concludes that the ICP -based fidelity calculation algorithm is an effective method to assess the similarity between solid models in the case of moderate deformation.

Table 2
Fidelity analysis for various geometric features and deformations.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (Grant Nos 11002139 and 81327803) and the Fundamental Research Funds for the Central Universities. The authors would like to thank Dr. Pengfei Shao and Dr. Ligang Liu of the University of Science and Technology of China, the laboratory students Qiumin Jia, Xiang Lv, Xiangqian Du and Yameng Li for their support of this work.

Author Contributions

Author Contributions

R.X.X., E.B.D., B.C., B.A., G.B. and S.W.S. conceived and designed research. H.L.W., L.Y. and Z.H.Z. performed experiment. B.L. and W.D.L. contributed to the materials and CT scanning. H.L.W., S.W.S. and X.M.Z. contributed to data processing and analyses. S.W.S., Y.X. and E.B.D. contributed to the development of fidelity calculation algorithm. S.W.S., G.B., E.B.D., and R.X.X. wrote the main manuscript text and reviewed the manuscript. All authors discussed the results and commented the manuscript.

Notes

Competing Interests

The authors declare that they have no competing interests.

Footnotes

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Erbao Dong, nc.ude.ctsu@gnodbe.

Ronald X. Xu, moc.liamtoh@dlanor.ux.

References

1. Chatterjee SK, Debmaulik T, Talukder B, Sen B, Sarangi B. Separation of pygopagus twins. Pediatric surgery international. 1986;1:74–78.
2. Inamdar SA, Goel SS, Subhedar VS. A thoracophagus conjoined twins with myelomeningocele: an unusual case. Int J Reprod Contracept Obstet Gynecol. 2012;1:61–63. doi: 10.5455/2320-1770.ijrcog001312. [Cross Ref]
3. Stuart GM, Black AE, Howard RF. The anaesthetic management of conjoined twins. Seminars in Pediatric Surgery. 2015;24:224–228. doi: 10.1053/j.sempedsurg.2015.06.005. [PubMed] [Cross Ref]
4. Spencer R. Theoretical and analytical embryology of conjoined twins: part i: embryogenesis. Clinical Anatomy. 2000;13:36–53. doi: 10.1002/(SICI)1098-2353(2000)13:1<36::AID-CA5>3.0.CO;2-3. [PubMed] [Cross Ref]
5. Hall JG. Twinning. The Lancet. 2003;362:735–743. doi: 10.1016/S0140-6736(03)14237-7. [PubMed] [Cross Ref]
6. Spitz L. Ethics in the management of conjoined twins. Seminars in Pediatric Surgery. 2015;24:263–264. doi: 10.1053/j.sempedsurg.2015.06.008. [PubMed] [Cross Ref]
7. Awasthi R, Iyengar R, Rege S, Jain N. Surgical management of pygopagus conjoined twins with spinal bifida. European Spine Journal. 2015;24:560–563. doi: 10.1007/s00586-014-3707-y. [PubMed] [Cross Ref]
8. Tannuri ACA, Batatinha JAP, Velhote MCP, Tannuri U. Conjoined twins: twenty years’ experience at a reference center in brazil. Clinics. 2013;68:371–377. doi: 10.6061/clinics/2013(03)OA14. [PMC free article] [PubMed] [Cross Ref]
9. Kamata M, et al. Case report anesthetic management of pygopagus conjoined twins: techniques to evaluate cross-circulation. Int J Clin Exp Med. 2016;9:6902–6910.
10. Watson SG, McHugh K. Conjoined twins: Radiological experience. Seminars in Pediatric Surgery. 2015;24:212–216. doi: 10.1053/j.sempedsurg.2015.07.002. [PubMed] [Cross Ref]
11. Petersen S, et al. 3d assisted prenatal sonographic diagnosis of dicephalic conjoined twins and subsequent planned vaginal delivery. Case Reports in Women’s Health. 2014;1:1–3.
12. Saito K, et al. Construction of a computed tomographic phantom for a japanese male adult and dose calculation system. Radiation and environmental biophysics. 2001;40:69–76. doi: 10.1007/s004110000082. [PubMed] [Cross Ref]
13. Castadot P, Geets X, Lee JA, Christian N, Grégoire V. Assessment by a deformable registration method of the volumetric and positional changes of target volumes and organs at risk in pharyngo-laryngeal tumors treated with concomitant chemo-radiation. Radiotherapy and Oncology. 2010;95:209–217. doi: 10.1016/j.radonc.2010.03.007. [PubMed] [Cross Ref]
14. Rode H, et al. Four decades of conjoined twins at red cross children’s hospital-lessons learned. South African Medical Journal. 2006;96:931–940. [PubMed]
15. Jr O, et al. Surgical experience with thirteen conjoined twins. Annals of surgery. 1988;208:299. doi: 10.1097/00000658-198809000-00007. [PubMed] [Cross Ref]
16. Rhodes JL, Yacoe M. Preoperative planning for the separation of omphalopagus conjoined twins the role of a multicomponent medical model. Journal of Craniofacial Surgery. 2013;24:175–177. doi: 10.1097/SCS.0b013e3182646b85. [PubMed] [Cross Ref]
17. Sachs E, Cima M, Williams P, Brancazio D, Cornie J. Three dimensional printing: rapid tooling and prototypes directly from a cad model. Journal of Engineering for Industry. 1992;114:481–488. doi: 10.1115/1.2900701. [Cross Ref]
18. Lipson, H. & Kurman, M. Fabricated: The new world of 3D printing (John Wiley & Sons, 2013).
19. Butscher A, Bohner M, Hofmann S, Gauckler L, Müller R. Structural and material approaches to bone tissue engineering in powder-based three-dimensional printing. Acta Biomaterialia. 2011;7:907–920. doi: 10.1016/j.actbio.2010.09.039. [PubMed] [Cross Ref]
20. Wu, W., DeConinck, A. & Lewis, J. A. Omnidirectional printing of 3d microvascular networks. Advanced Materials23 (2011). [PubMed]
21. Gao Q, He Y, Fu J-Z, Liu A, Ma L. Coaxial nozzle-assisted 3d bioprinting with built-in microchannels for nutrients delivery. Biomaterials. 2015;61:203–215. doi: 10.1016/j.biomaterials.2015.05.031. [PubMed] [Cross Ref]
22. Murphy SV, Atala A. 3d bioprinting of tissues and organs. Nature biotechnology. 2014;32:773–785. doi: 10.1038/nbt.2958. [PubMed] [Cross Ref]
23. Klammert U, et al. 3d powder printed calcium phosphate implants for reconstruction of cranial and maxillofacial defects. Journal of Cranio-Maxillofacial Surgery. 2010;38:565–570. doi: 10.1016/j.jcms.2010.01.009. [PubMed] [Cross Ref]
24. Wang, H. et al. Three-dimensional virtual surgery models for percutaneous coronary intervention (pci) optimization strategies. Scientific reports5, 10945–10945 (2014). [PMC free article] [PubMed]
25. Shen S, et al. Freeform fabrication of tissue-simulating phantoms by combining three-dimensional printing and casting. Proc. SPIE. 2016;9700:970009–6. doi: 10.1117/12.2212305. [Cross Ref]
26. Spitz L, Kiely E. Experience in the management of conjoined twins. British journal of Surgery. 2002;89:1188–1192. doi: 10.1046/j.1365-2168.2002.02193.x. [PubMed] [Cross Ref]
27. Schneider U, Pedroni E, Lomax A. The calibration of ct hounsfield units for radiotherapy treatment planning. Physics in medicine and biology. 1996;41:111. doi: 10.1088/0031-9155/41/1/009. [PubMed] [Cross Ref]
28. Hebb AO, Poliakov AV. Imaging of deep brain stimulation leads using extended hounsfield unit ct. Stereotactic and functional neurosurgery. 2009;87:155–160. doi: 10.1159/000209296. [PubMed] [Cross Ref]
29. MacFarlane M, Rosen J, Hannaford B, Pellegrini C, Sinanan M. Force-feedback grasper helps restore sense of touch in minimally invasive surgery. Journal of Gastrointestinal Surgery. 1999;3:278–285. doi: 10.1016/S1091-255X(99)80069-9. [PubMed] [Cross Ref]
30. Azar FS, Metaxas DN, Miller RT, Schnall MD. Methods for predicting mechanical deformations in the breast during clinical breast biopsy. Methods. 2000;4:8–2000.
31. Norrish K, Hutton JT. An accurate x-ray spectrographic method for the analysis of a wide range of geological samples. Geochimica et cosmochimica acta. 1969;33:431–453. doi: 10.1016/0016-7037(69)90126-4. [Cross Ref]
32. Galperin A, et al. Radiopaque iodinated polymeric nanoparticles for x-ray imaging applications. Biomaterials. 2007;28:4461–4468. doi: 10.1016/j.biomaterials.2007.06.032. [PubMed] [Cross Ref]
33. Lee K-S, Kim S-H. Non-uniform deformation of an stl model satisfying error criteria. Computer-Aided Design. 2010;42:238–247. doi: 10.1016/j.cad.2009.11.006. [Cross Ref]
34. Huang S-H, Zhang L-C, Han M. An effective error-tolerance slicing algorithm for stl files. The International Journal of Advanced Manufacturing Technology. 2002;20:363–367. doi: 10.1007/s001700200164. [Cross Ref]
35. Nee A, Fu M, Fuh J, Lee K, Zhang Y. Automatic determination of 3-d parting lines and surfaces in plastic injection mould design. CIRP Annals-Manufacturing Technology. 1998;47:95–98. doi: 10.1016/S0007-8506(07)62793-0. [Cross Ref]
36. Yantis, S. Visual perception: Essential readings (Psychology Press, 2001).
37. Fechner, G. et al. Elements of psychophysics (Holt, Rinehart & Winston, 1965).
38. Bassoli E, Gatto A, Iuliano L, Grazia Violante M. 3d printing technique applied to rapid casting. Rapid Prototyping Journal. 2007;13:148–155. doi: 10.1108/13552540710750898. [Cross Ref]
39. Lauriks, W. Low density foams. Low density cellular plastics: Physical basis of behaviour319 (2012).
40. Brassard, G. & Salvail, L. Lecture notes in computer science. In Advances in Cryptology Eurocrypt93, vol. 765, 410–23 (1994).
41. Arun, K. S., Huang, T. S. & Blostein, S. D. Least-squares fitting of two 3-d point sets. IEEE Transactions on pattern analysis and machine intelligence 698–700 (1987). [PubMed]
42. Ying S, Peng J, Du S, Qiao H. A scale stretch method based on icp for 3d data registration. IEEE Transactions on Automation Science and Engineering. 2009;6:559–565. doi: 10.1109/TASE.2009.2021337. [Cross Ref]
43. Yan P, Bowyer KW. A fast algorithm for icp-based 3d shape biometrics. Computer Vision and Image Understanding. 2007;107:195–202. doi: 10.1016/j.cviu.2006.11.001. [Cross Ref]
44. Sederberg TW, Parry SR. Free-form deformation of solid geometric models. ACM SIGGRAPH computer graphics. 1986;20:151–160. doi: 10.1145/15886.15903. [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group