Abstract
Purpose
Three-dimensional (3D) printing is now widely available, and its potential applications to surgery are limitless. However, 3D printing is presently performed at only a few large institutions. We developed a 3D printing workflow using an affordable 3D printer and open-source 3D printing software. We tested whether this combination could produce a model that reliably reflects real bone.
Methods
We performed a prospective study with a target sample size of seven. Patients with distal radius fractures were enrolled from October 2021 to February 2022. The 3D-printed models of the fractures were produced using open-source software (3D Slicer [Surgical Planning Laboratory, Harvard Medical School] and Cura [Ultimaker]) and a $600 printer. The anterior-to-posterior (AP) and radial-to-ulnar (RU) widths of the fracture sites were measured on computed tomography (CT) images, in 3D printed models, and in real bones (during surgery). Surgery was simulated using the 3D models; the locations and profiles of implants were compared to those placed during real surgery, which was performed without simulation data.
Results
The fracture AP and RU widths did not differ significantly among the CT, 3D model, and real bone measurements. Interclass correlation coefficients indicated that the measurements were reliable (0.943 [p<0.001] and 0.917 [p<0.001], respectively). When the implant profiles of the simulations and surgical procedures were compared, only the most distal (radial) screws were significantly longer in the simulations (p=0.016). The plate location also differed significantly (p=0.043).
Three-dimensional (3D) printing has made great strides in recent years. Although the models cannot yet replace organs, they are very useful in modelling osseous structures. The affordability and accessibility of 3D printing have increased [1], but it has not yet achieved everyday usage.
The process of 3D printing involves two steps. First, stacked computed tomography (CT) or magnetic resonance images in Digital Imaging and Communications in Medicine (DICOM) format are converted to stereolithography (STL) that describes the 3D surface geometry. Then, a 3D model is printed. Some printers now cost less than $1,000. However, most surgeons are not familiar with the software required for the first step of 3D printing. Most medical researchers use Mimics (Materialise, Leuven, Belgium) [2-6], which shows excellent performance. However, a license costs $10,000 to $20,000 depending on the modules included.
Therefore, we developed a 3D printing workflow using open-source software and an affordable desktop printer. This can be achieved in a typical office setting. As George et al. [7] proposed, we need to establish the accuracy and reproducibility of our workflow. We performed a prospective study on patients with distal radius fractures to investigate whether the printed models resembled actual bones.
Ethics statement: This retrospective study was approved by the Institutional Review Board of Seoul National University College of Medicine and Seoul National University Hospital (No. H-2105-093-1219). Informed consent was written by all enrolled patients.
A distal radius fracture displacement of >2 mm is considered significant [8]. Screw sizes vary in steps of 1 or 2 mm. We thus determined the non-inferiority margin (δ), which was 1. The standard deviation of the anterior-to-posterior (AP) size was the same as in a previous cadaver study [9]. Our hypotheses were as follows:
H0 size at fracture site (3D simulation) minus size at fracture site (real surgery) ≤−|δ|
H1 size at fracture site (3D simulation) minus size at fracture site (real surgery) >−|δ|.
Patients were consecutively enrolled beginning in October 2021. We invited adult patients with distal radius fractures admitted for surgery to participate. Fractures with severe comminutions (AO/OTA classes 23-C2 and C3) were excluded [12]. Seven patients gave written informed consent by February 2022.
An IQon Spectral CT scanner (Philips, Amsterdam, Netherlands) was used. Preoperative wrist CT yields images in DICOM format that can be converted into STL images using 3D Slicer (ver. 4.11; Surgical Planning Laboratory, Harvard Medical School, Boston, MA, USA) [13]. Segmentation converts stacked two-dimensional images (in DICOM format) into 3D surface models (STL format) [14]. The process is as follows: the DICOM files to be converted are selected. Usually, the thinnest axial images are ideal. We obtained 299 images per patient, each with a thickness of 0.8 mm. Then, a segmentation editor is used to select an appropriate threshold for the region of interest; only the osseous portion should undergo volumetric reconstruction. After clicking on the 3D icon, the model appears. If only the distal radius is of interest, the “scissors” and “island” functions can be used to crop away other carpal bones and/or the ulna. Finally, the final 3D model is exported as an STL file (Fig. 1). The file must then be converted into geometric (G)-code to command the printer. We used Cura (ver. 4.11.0; Ultimaker, Utrecht, the Netherlands) to this end. This open-source slicing application converts STL to G-code that instructs the 3D printer. The default settings should be used, i.e. 0.2-mm thickness and 10% occupied volume. Clicking the slice icon yields G-code that is then sent to an Ender-5-Plus 3D printer (Shenzhen Creality 3D Technology, Shenzhen, China) using polylactic acid filament as ink. The product is a printed 3D model of the fractured distal radius. Fig. 2 summarizes these steps.
The AP and radial-to-ulnar (RU) widths (in mm) of the fracture sites were measured on CT images, printed models, and real bones (during surgery) (Fig. 3).
(1) To obtain measurements on CT images, the “ruler” of the picture archiving and communication system was used. The fracture site was selected in an axial cut by browsing the coronal and sagittal sequences. At the appropriate axial image, AP and RU widths between the farthest cortices were measured (Fig. 3A).
(2) After 3D models were produced, the models were separated at the fracture sites. We used a digital caliper to measure the AP and RU widths at the distal end of proximal fragments (fracture site) in the printed models (Fig. 3B).
(3) We used the same caliper (now sterilized) to measure real bones during surgery. Similarly, these measurements were also made by holding the distal end of proximal fragments with the caliper before the reduction (Fig. 3C). A gentle periosteal dissection at the fracture site facilitated the measurement using a caliper, which in turn made the reduction also simple.
We recorded the times required for modeling, printing, and simulation. Modeling refers to the conversion of the DICOM file to G-code. Reduction and fixation of printed distal radii (using plates and screws) were simulated by the junior orthopedic surgeon (author 1). During simulation and surgery, narrow or standard 2.4-mm variable angle locking compression plates and volar distal radius plates (DePuy Synthes, Warsaw, IN, USA) with two shaft holes were used. Each plate was secured with 2.4-mm variable angle locking screws. Plate locations were recorded by measuring three distances (location, sLOC [and rLOC] 1, 2, and 3) (Table 1, Fig. 4A) from the radial margins of the fracture site to the wire holes. If all three distances were similar during simulation and real surgery, we assumed that the plate locations would be similar. The screw lengths (ssc [and rsc] 1 to ssc [and rsc] 8) (Table 2, Fig. 4B) were recorded. The senior orthopedic surgeon (author 2), who was blinded to the simulation data, performed all surgeries. The corresponding plate profiles and screw sizes during the surgery were recorded.
Measurements were compared using the signed-rank test. When no significant difference was apparent, measurement consistency was further evaluated by calculating interclass correlation coefficients (ICCs). All statistical tests were two-tailed and p-values of <0.05 were considered significant.
Modeling required an average of 8 minutes (range, 6-9 minutes). Printing required an average of 149 minutes (range, 99–225 minutes) but was not supervised. The simulation required an average of 19 minutes (range, 12–26 minutes).
The CT, 3D model, and real bone AP and RU widths did not differ significantly (Table 3). The ICC revealed excellent consistency among the measurements.
For cases 1 to 5, narrow plates were used for both simulation and real surgery; while standard plates were employed for cases 6 and 7. Distance 2 was significantly shorter during real surgery than simulation (mean sLOC2, 16.04 mm vs. mean rLOC2, 14.69 mm; p=0.043) (Table 1). In terms of screw sizes, the most distal and radial screw, i.e. screw 1, was significantly longer in the simulation compared with real surgery (mean ssc1, 17.42 mm vs. mean rsc1, 14.29 mm; p=0.016) (Table 2). The sizes of other screws did not differ significantly, but only the two screws for the shaft holes had ICCs indicating excellent consistency (0.927 and 0.842, respectively).
In this study, we introduced an efficient and affordable 3D printing workflow. The only cost is associated with the hardware (about $600). In our prospective study of seven distal radius fractures, the models matched the real bones excellently. CT is generally used to measure bone size [15]. The integrity of CT DICOM files is extremely high. However, a divergence in a given DICOM file may be seen among institutes due to the STL file conversion process [16]. In another study, a single institute using nine different software packages noted no significant divergence [14]. This supports the necessity of establishing the accuracy and reproducibility for each 3D printing workflow [7]. The 3D Slicer free open-source software is supported by international developers [13]. The pioneering Quantitative Imaging Network (QIN) uses 3D Slicer for personalized cancer therapy. The QIN group includes the Brigham and Women’s Hospital, the University of Iowa, and Massachusetts General Hospital. Bagaria and Chaudhary [17] reported their experience with 3D Slicer; subjective feedback from five surgeons indicated that the models accurately reflected the actual anatomies. Our study is the first to show that fractured bone geometry is preserved by a workflow using a 3D Slicer.
3D printing technology is evolving rapidly but is currently used by only a few large centers [18]. Perceived barriers to widespread use include high monetary and time costs. Printers range in price from $100–$10,000. Most devices that cost less than $1,000 can produce bone models. Open-source software, such as 3D Slicer and Cura, show excellent modeling performance, and there are no additional labor costs because clinicians perform all functions. In this study, the average modeling time was 8 minutes. Printing took over 2 hours but was not supervised. One kilogram of ink costs $10, and there are no maintenance costs. Also, the printer is small and we use it to create one or two 3D models daily in our office; even small private practices can use our workflow.
In this study, we focused on introducing and validating our 3D printing workflow rather than an emphasis on preoperative planning. The senior surgeon lacked knowledge of the simulation performed by the junior surgeon. The plate location differed significantly between the simulation and real surgery, but the length of only one screw differed significantly. Preoperative planning is essential; previous studies on die-punch, intercondylar distal humeral, and complex proximal humeral fractures found that preoperative planning (via 3D printing) significantly reduced blood loss, operative time, and the number of intraoperative fluoroscopies [4-6]. Yoshii et al. [19] reported that the 3D preoperative plans (without printing) for the treatment of distal radius fractures exhibited errors of about 2 mm. A recent systematic review on preoperative planning prior to reverse total shoulder arthroplasty also revealed minimal real-world deviations from such plans [20]. The utility of 3D printing when performing bone surgery is thus well-recognized, but accessibility must be improved.
There are limitations to our study. This is a prospective study, but the minimum number of patients according to the power analysis were enrolled. Second, the clinical evidence of simulation using this workflow has not been properly investigated in this study. In the absence of soft tissue, the simulation using the model inherently differs from the actual surgery. The simulation cannot reflect the surgical approach, dissection, and reduction of fracture when there are ligaments and tendons surrounding the target. A more elaborate model reflecting these obstacles should be investigated further. Also, the simulator and the surgeon who operated the actual cases were different. Ideally, the simulator and the surgeon should be the same person, and parameters like operation time, bleeding amount, and radiographic measurement should be compared. However, this kind of investigation belongs to clinical investigation of medical devices. This research will provide the background for our next research which will be the clinical investigation of medical devices for the treatment effect of our 3D printing workflow.
Notes
References
1. Skelley NW, Smith MJ, Ma R, Cook JL. Three-dimensional printing technology in orthopaedics. J Am Acad Orthop Surg. 2019; 27:918–25.
2. Kim HN, Liu XN, Noh KC. Use of a real-size 3D-printed model as a preoperative and intraoperative tool for minimally invasive plating of comminuted midshaft clavicle fractures. J Orthop Surg Res. 2015; 10:91.
3. Kim JW, Lee Y, Seo J, et al. Clinical experience with three-dimensional printing techniques in orthopedic trauma. J Orthop Sci. 2018; 23:383–8.
4. Chen C, Cai L, Zhang C, Wang J, Guo X, Zhou Y. Treatment of die-punch fractures with 3D printing technology. J Invest Surg. 2018; 31:385–92.
5. You W, Liu LJ, Chen HX, et al. Application of 3D printing technology on the treatment of complex proximal humeral fractures (Neer3-part and 4-part) in old people. Orthop Traumatol Surg Res. 2016; 102:897–903.
6. Zheng W, Su J, Cai L, et al. Application of 3D-printing technology in the treatment of humeral intercondylar fractures. Orthop Traumatol Surg Res. 2018; 104:83–8.
7. George E, Liacouras P, Rybicki FJ, Mitsouras D. Measuring and establishing the accuracy and reproducibility of 3D printed medical models. Radiographics. 2017; 37:1424–50.
8. Ng CY, McQueen MM. What are the radiological predictors of functional outcome following fractures of the distal radius? J Bone Joint Surg Br. 2011; 93:145–50.
9. Gasse N, Lepage D, Pem R, et al. Anatomical and radiological study applied to distal radius surgery. Surg Radiol Anat. 2011; 33:485–90.
10. Centre for Clinical Research and Biostatics (CCRB). Sample size estimation [Internet]. Hong Kong: CCRB;c2021 [cited 2021 May 1]. Available from: https://www2.ccrb.cuhk.edu.hk/stat/Means.htm.
11. Faul F, Erdfelder E, Buchner A, Lang AG. Statistical power analyses using G*Power 3.1: tests for correlation and regression analyses. Behav Res Methods. 2009; 41:1149–60.
12. Wæver D, Madsen ML, Rölfing JH, et al. Distal radius fractures are difficult to classify. Injury. 2018; 49 Suppl 1:S29–32.
13. Fedorov A, Beichel R, Kalpathy-Cramer J, et al. 3D Slicer as an image computing platform for the Quantitative Imaging Network. Magn Reson Imaging. 2012; 30:1323–41.
14. Kamio T, Suzuki M, Asaumi R, Kawai T. DICOM segmentation and STL creation for 3D printing: a process and software package comparison for osseous anatomy. 3D Print Med. 2020; 6:17.
15. Rathnayaka K, Momot KI, Noser H, et al. Quantification of the accuracy of MRI generated 3D models of long bones compared to CT generated 3D models. Med Eng Phys. 2012; 34:357–63.
16. Huotilainen E, Jaanimets R, Valášek J, et al. Inaccuracies in additive manufactured medical skull models caused by the DICOM to STL conversion process. J Craniomaxillofac Surg. 2014; 42:e259–65.
17. Bagaria V, Chaudhary K. A paradigm shift in surgical planning and simulation using 3Dgraphy: experience of first 50 surgeries done using 3D-printed biomodels. Injury. 2017; 48:2501–8.
18. Horas K, Hoffmann R, Faulenbach M, Heinz SM, Langheinrich A, Schweigkofler U. Advances in the preoperative planning of revision trauma surgery using 3D printing technology. J Orthop Trauma. 2020; 34:e181–6.
Table 1.
sLOC, simulated location; rLOC, real surgical location (see Fig. 4A. for the numbering); LE, left; R, right; H, holes.
*p<0.05.
Table 2.
ssc, sizes of screws during the simulation; rsc, sizes of screws during the real surgery. The screws are numbered from the most radial (most distal) hole (ssc1 and rsc1). The ssc7 (and rsc7) and ssc8 (and rsc8) screws are for the shaft holes (see Fig. 4B for the numbering).
*p<0.05.