1. Introduction
The abdominal area is characterized by both bone and soft tissues, and a large amount of free gas, which is: (a) in the lumen of the bowel; (b) in the peritoneal cavity but outside the lumen of the bowel (due to bowel perforation); and (c) in the soft tissues in case of efizoma. The bone, the soft tissues and the free gas or air attenuate the X-rays quite differently, which results in difficulties reproducing the attenuation correctly in the 3D printing process, when anthropomorphic abdomen phantoms are produced. The fused deposition modelling (FDM) technique offers a low-cost method for the production of anthropomorphic radiological phantoms. Commercial FDM materials, such as PLA (Polylactic acid) and ABS (Acrylonitrile Butadiene Styrene) do not suitably represent the X-ray properties of the human bone tissues. For instance, a high-resolution 3D model of pelvic bones was produced for surgical planning [
[1]- Qu Y.
- Li X.
- Yan Z.
- Zhao L.
- Zhang L.
- Liu C.
- et al.
Surgical planning of pelvic tumor using multi-view CNN with relation-context representation learning.
], however, bone structures were replicated by ABS [
[2]- Yin J.
- Li M.Q.
- Dai G.L.
- Zhou H.Z.
- Ma L.
- Zheng Y.X.
3D Printed Multi-material Medical Phantoms for Needle-tissue Interaction Modelling of Heterogeneous Structures.
], a material, which is inappropriate for CT imaging applications. Similarly, PLA was used by Koh et al. [
[3]Koh JC, Jang YK, Seong H, Lee KH, Jun S, Choi JB. Creation of a three-dimensional printed spine model for training in pain procedures. J Int Med Res. 2021;49:3000605211053281 DOI: 10.1177/03000605211053281.
] to reproduce the bone and intervertebral disks. All these models may be successfully used to educate clinicians on surgical techniques and patient preparation. However, for optimisation of reconstruction algorithms [
[4]- Hatamikia S.
- Biguri A.
- Kronreif G.
- Kettenbach J.
- Russ T.
- Furtado H.
- et al.
Optimization for customized trajectories in cone beam computed tomography.
], image segmentation [
[5]- Holmes R.B.
- Negus I.S.
- Wiltshire S.J.
- Thorne G.C.
- Young P.
Creation of an anthropomorphic CT head phantom for verification of image segmentation.
], testing of novel techniques such as dual-energy CT [
[6]Niebuhr NI, Johnen W, Echner G, Runz A, Bach M, Stoll M, et al. The ADAM-pelvis phantom-an anthropomorphic, deformable and multimodal phantom for MRgRT. Phys Med Biol. 2019;64:04NT5 DOI: 10.1088/1361-6560/aafd5f.
], performing accurate radiation dosimetry [
[7]- Giacometti V.
- King R.B.
- McCreery C.
- Buchanan F.
- Jeevanandam P.
- Jain S.
- et al.
3D-printed patient-specific pelvis phantom for dosimetry measurements for prostate stereotactic radiotherapy with dominant intraprostatic lesion boost.
] and QA applications, the bone structures should be represented by 3D printing materials with X-ray properties close to the X-ray properties of the respective bone tissues.
In the last few years, efforts have been put towards developing new materials by the inclusion of add-ins containing a high Z-value element. Reports showed that ABS filaments were doped with BaSO
4, CaCO
3 [
[8]- Hamedani B.A.
- Melvin A.
- Vaheesan K.
- Gadani S.
- Pereira K.
- Hall A.F.
Three-dimensional printing CT-derived objects with controllable radiopacity.
], Bi
2O
3 [
[9]- Darling C.J.
- Curtis C.
- Sciacca B.J.
- Sarkar K.
- Smith D.A.
Fused filament fabrication of complex anatomical phantoms with infill-tunable image contrast.
], and CaTiO
3 [
[10]Price G, Biglin ER, Collins S, Aitkinhead A, Subiel A, Chadwick AL, et al. An open source heterogeneous 3D printed mouse phantom utilising a novel bone representative thermoplastic. Phys Med Biol. 2020;65:10NT02 DOI: 10.1088/1361-6560/ab8078.
] to bone structures, such as pelvis, spine, and cortical bone, with FDM technology. In all these approaches, the bone structures were segmented prior to 3D printing, rather than using clinical CT data to directly print the structures. In addition, all these approaches require knowledge about both materials and necessary equipment to initially produce the new filament. Further, in some cases, the obtained filaments as the one achieved from ABS pellets doped with BaSO
4 resulted in a maximum Hounsfield Unit (HU) of cortical bone in the phantom lower than that of the filament [
[8]- Hamedani B.A.
- Melvin A.
- Vaheesan K.
- Gadani S.
- Pereira K.
- Hall A.F.
Three-dimensional printing CT-derived objects with controllable radiopacity.
]. The inclusion of water in the produced bone did not result in air elimination. Other researchers used PLA doped with gravitational stone [
11- Goodall S.K.
- Rampant P.
- Smith W.
- Waterhouse D.
- Rowshanfarzad P.
- Ebert M.A.
Investigation of the effects of spinal surgical implants on radiotherapy dosimetry: A study of 3D printed phantoms.
,
12- Okkalidis N.
- Marinakis G.
Technical Note: Accurate replication of soft and bone tissues with 3D printing.
], Bi
2O
3 [
[13]- Arnold J.
- Sarkar K.
- Smith D.
3D printed bismuth oxide-polylactic acid composites for radio-mimetic computed tomography spine phantoms.
] and ferromagnetic particles [
14- Savi M.
- Villani D.
- Andrade M.A.B.
- Rodrigues O.
- Potiens M.P.A.
Study on attenuation of 3D printing commercial filaments on standard X-ray beams for dosimetry and tissue equivalence.
,
15- Mille M.M.
- Griffin K.T.
- Maass-Moreno R.
- Lee C.
Fabrication of a pediatric torso phantom with multiple tissues represented using a dual nozzle thermoplastic 3D printer.
]. Amongst the preferred filaments for 3D printing pelvis bone is the StoneFil. Giacometti et al. 2021 [
[7]- Giacometti V.
- King R.B.
- McCreery C.
- Buchanan F.
- Jeevanandam P.
- Jain S.
- et al.
3D-printed patient-specific pelvis phantom for dosimetry measurements for prostate stereotactic radiotherapy with dominant intraprostatic lesion boost.
] represented the cortical bone by 100 % infill solid StoneFil filament, while the trabecular bone corresponded to an infill percentage of 65 %. Okkalidis et al. [
[12]- Okkalidis N.
- Marinakis G.
Technical Note: Accurate replication of soft and bone tissues with 3D printing.
] proposed a method to control the filament extrusion rate during the 3D printing process. The maximum HU values obtained with their approach are an average of 523 ± 259 HU, which is limiting the application of the cortical bone. Iron PLA was used by Mille et al. [
[15]- Mille M.M.
- Griffin K.T.
- Maass-Moreno R.
- Lee C.
Fabrication of a pediatric torso phantom with multiple tissues represented using a dual nozzle thermoplastic 3D printer.
] in printing a torso paediatric phantom. However, all bone structures were printed with Iron PLA with an infill density of 50 %.
The methods for controlling an FDM 3D printer for the replication of the soft and bone tissues presented so far in the literature comprise the use of printing patterns and specific infill densities [
16- Madamesila J.
- McGeachy P.
- Barajas J.E.V.
- Khan R.
Characterizing 3D printing in the fabrication of variable density phantoms for quality assurance of radiotherapy.
,
173D printing methods for radiological anthropomorphic phantoms.
], or by depositing estimated amounts of melted filaments per voxel volume [
12- Okkalidis N.
- Marinakis G.
Technical Note: Accurate replication of soft and bone tissues with 3D printing.
,
18A novel 3D printing method for accurate anatomy replication in patient-specific phantoms.
,
19- Mei K.
- Geagan M.
- Roshkovan L.
- Litt H.I.
- Gang G.J.
- Shapira N.
- et al.
Three-dimensional printing of patient-specific lung phantoms for CT imaging: Emulating lung tissue with accurate attenuation profiles and textures.
]. In all the aforementioned cases, a single material has been used for the fabrication of each organ. A basic limitation of this method is the presence of air gaps in the produced phantom since the more or less melted filament is deposited on each predefined voxel according to the voxeĺs HU. Hence, the properties of the produced phantoms derived from the common filament materials are limited concerning applications in X-ray imaging modalities. Recently, Tino et al. [
[20]- Tino R.
- Yeo A.
- Brandt M.
- Leary M.
- Kron T.
The interlace deposition method of bone equivalent material extrusion 3D printing for imaging in radiotherapy.
] presented a bone-equivalent HU phantom, where two different filaments were combined by adjusting the printed layer ratio of the two materials based on a standard number of layers. However, their method can replicate only 13 different mean HU for each bone. Karin et al. [
[21]- Kairn T.
- Zahrani M.
- Cassim N.
- Livingstone A.G.
- Charles P.H.
- Crowe S.B.
Quasi-simultaneous 3D printing of muscle-, lung- and bone-equivalent media: a proof-of-concept study.
] investigated the use of a dual-head printer to deposit two different filaments, PLA and StoneFil PLA, at several different in-fill densities, to achieve quasi-simultaneous 3D printing of muscle-, lung- and bone-equivalent media. The approach was applied to print one thoracic and one cranial phantom slab. However, the two extruders were used to print different tissue: soft or bone, rather than combining the two filaments in the fabrication of a bone structure.
A study by Ma et al. [
[22]- Ma X.
- Buschmann M.
- Unger E.
- Homolka P.
Classification of X-Ray Attenuation Properties of Additive Manufacturing and 3D Printing Materials Using Computed Tomography From 70 to 140 kVp.
] with 26 FDM materials showed that pure PLA, PET (polyethylene terephthalate) and PETG (polyethylene terephthalate glycol) filaments have higher HUs compared to the HU of the breast tissues [
[23]- Dukov N.
- Bliznakova K.
- Okkalidis N.
- Teneva T.
- Encheva E.
- Bliznakov Z.
Thermoplastic 3D printing technology using a single filament for producing realistic patient-derived breast models.
]. On the other hand, the polypropylene (PP) showed that the adipose tissue can be represented adequately. The latter material, when combined with PLA in appropriate ratios, is suitable to represent all the HU values of the soft tissues in the abdomen region, while failing in the case of bone structures. Finally, the use of PLA filled with 50 % powdered stone filament can be used to replicate the bone tissues.
The overall objective of the interdisciplinary team at the Medical University of Varna is the development of physical anthropomorphic abdominal phantoms for CT applications. The first step in reaching this goal is the development of a methodology for creating physical bone phantoms from appropriate materials, which can be utilized to replicate accurately the X-ray properties of the bone tissue. This study presents a novel approach employing a custom-made dual motor filament extrusion setup for the fabrication of realistic and radiologically equivalent bone tissues. The proposed approach mixes two commercially available filaments during the printing process.
4. Discussion
In this study, we investigated the simultaneously use of two FDM filaments of a different type, i.e. PLA and Stonefil (a filament consisting of 50 % PLA and 50 % powdered stone) during a 3D printing process for the replication of bone tissues to be used in an anthropomorphic abdomen phantom. A custom-made software was used for the direct association of the images with the 3D printing process depositing melted filament voxel-by-voxel according to the corresponding HU from the medical images. The extruder attached to the printer could employ-three different filaments during a 3D printing process, while there was a single nozzle as an exit for the melted filaments (
Fig. 1). In this type of FDM printer, when two materials of different compositions are used simultaneously, the result will be an extrusion of the dominant material only. Therefore, when PLA and Stonefil filaments were used simultaneously, the final extrusion result was only PLA. The reason is the lower viscosity of the melted PLA compared to the melted Stonefil in the heating chamber of the extruder. By increasing the input force of the second material enough to overcome the blocking from the dominant material, it was possible to have an extrusion of both filaments during a 3D printing process.
Three phantoms were produced based on three different scenarios for the evaluation of the new approach. The results from the first scenario (
Fig. 8a) showed that it is possible to produce entities with a homogeneous result instead of using specific infill patterns and infill densities [
8- Hamedani B.A.
- Melvin A.
- Vaheesan K.
- Gadani S.
- Pereira K.
- Hall A.F.
Three-dimensional printing CT-derived objects with controllable radiopacity.
,
16- Madamesila J.
- McGeachy P.
- Barajas J.E.V.
- Khan R.
Characterizing 3D printing in the fabrication of variable density phantoms for quality assurance of radiotherapy.
]. A high degree of correlation was achieved between the various filaments ratios with the corresponding HU resulting in R2 values between 0.9877 up to 0.9998 (
Table 1). This means that the proposed custom-made dual motor filament extrusion setup demonstrated a high performance in keeping constant the chosen ratios during the whole printing process. The CT scan images of the calibration cubes and phantoms (
Fig. 3 and
Fig. 7) showed a high uniformity after combining the two filaments, while the retrieved histograms from the three phantoms (
Fig. 8) revealed the expected uniformity under the chosen parameters shown in the section ‘Phantom’, (a) to (c) scenarios.
It has been shown that there are two ways to produce the same HU for each entity, either by changing the mixing ratio of the two materials or by using a single material and selecting a constant but lower filament extrusion rate per voxel (
Fig. 4). For Phantom 1, the chosen mixing composition was 80 % Stonefil and 20 % PLA with 0.04 mm extrusion rate per voxel resulting in an average value of 700 HU as shown in
Fig. 8a. The same mean value at 700 HU could be achieved by using only Stonefil with a constant filament extrusion rate of approximately 0.0375. The results from Phantom 2 demonstrated that it is also possible to print an entity of a phantom by manipulating the filament extrusion rate per voxel based on the correlated HU of the medical images. The whole hip bone was printed line-by-line (
Fig. 5), controlling the extrusion filament rate voxel-by-voxel [
[18]A novel 3D printing method for accurate anatomy replication in patient-specific phantoms.
] and using the least square line fitted for the mixing composition of 50 % Stonefil and 50 % PLA with a maximum value at 500 HU, approximately (
Fig. 4). Although, the replicated morphology was visually close to an actual hip bone (
Fig. 7b), it can be noticed that the replicated HUs start from 100 HUs approximately, while much lower HUs can be observed in the original patient histogram (
Fig. 5b). This leads to the conclusion that the response of the filament extrusion rate per voxel is low enough to create high concentrations of the Stonefil filament inside the cancellous bone and lower concentrations on the cortical bone making the printed boundaries of the two entities visible after the CT scanning (
Fig. 7b). The first reason is that the specific extruder has a wide heating chamber where the two materials are concentrated before the extrusion creating delays on the deposition of the melted materials per voxels. The second reason is that the replicated entity is thin enough to affect the response of the proposed approach. Finally, the results from Phantom 3 showed that it is possible to replicate bones by 3D printing each entity under different parameters (
Fig. 6c), and still matching the anatomy of a realistic hip bone (
Fig. 5). The cancellous bone (the yellow contour (
Fig. 5a)) was printed with a constant mixing composition of 30 % Stonefil and 70 % PLA with 0.0375 extrusion rate per voxel resulting in an average value of 300 HU. The cortical bone (i.e. the red contour removing the yellow contour (
Fig. 5)) was printed using a concentric/perimetric printing pattern [
[23]- Dukov N.
- Bliznakova K.
- Okkalidis N.
- Teneva T.
- Encheva E.
- Bliznakov Z.
Thermoplastic 3D printing technology using a single filament for producing realistic patient-derived breast models.
] and 100 % Stonefil with 0.04 mm extrusion rate resulting in an average value of 800 HU.
This study demonstrated the possibility to combine two different filaments during the 3D printing process and to produce radiologically anthropomorphic bone tissues. The use of the Stonefil filament was chosen since the produced HU are close to the hip bone of the pelvis. Experiments with metallic filaments were performed as well, since they have been previously investigated and presented in the literature [
[20]- Tino R.
- Yeo A.
- Brandt M.
- Leary M.
- Kron T.
The interlace deposition method of bone equivalent material extrusion 3D printing for imaging in radiotherapy.
]. It has been confirmed that the metallic nature of the filaments can create higher HU values than the Stonefil, however, this may lead to streak artefacts in the X-ray images and may corrupt the calibration process and consequently the final produced phantom. This study showed that it is possible to produce the same HU using two different material compositions. This means that the radiological properties of a bone produced using FDM technology can be replicated in more than one way. We are planning to investigate the use of other materials than PLA such as ABS, HIPS, PETG, PP etc. with Stonefil or other materials of high contrast such as ABS doped with Sulfate Barium. This will provide the ability to manipulate the attenuation coefficients of a phantom voxel-by-voxel but also its mechanical properties.
Until now, it was common to use more than one material for 3D printing different tissues such as in the fabrication of an anthropomorphic breast phantom [
24di Franco F, Mettivier G, Sarno A, Varallo A, Russo P. Manufacturing of physical breast phantoms with 3D printing technology for X-ray breast imaging. IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC) 2019. p. 1-5 DOI: 10.1109/nss/mic42101.2019.9059986.
,
25- Varallo A.
- Sarno A.
- Castriconi R.
- Mazzilli A.
- Loria A.
- del Vecchio A.
- et al.
Fabrication of 3D printed patient-derived anthropomorphic breast phantoms for mammography and digital breast tomosynthesis: Imaging assessment with clinical X-ray spectra.
,
26Bliznakova K. The advent of anthropomorphic three-dimensional breast phantoms for X-ray imaging. Phys Med. 2020;79:145-61 DOI: S1120-1797(20)30299-4 doi: 10.1016/j.ejmp.2020.11.025.
], where they used an FDM printer with two extruders and two printing filaments PET and ABS. A method close to the proposed approach was introduced by Tino et al. [
[20]- Tino R.
- Yeo A.
- Brandt M.
- Leary M.
- Kron T.
The interlace deposition method of bone equivalent material extrusion 3D printing for imaging in radiotherapy.
], where a bone-equivalent HU is produced by two different materials and by adjusting the printed layer ratio to a standard volume. However, this method can replicate only 13 different mean HU for each entity, while the method proposed in this work presented a much higher resolution and the ability to print an entity voxel-by-voxel.
Other approaches for producing bone structures with X-ray properties close to these of human bones are 3D-printed in gypsum, a calcium sulfate mineral (CaSO
4·2H
2O), by using binder jet technology and PolyJet 3D printing technology combined with new curable resins. Of these approaches, the most popular is the use of gypsum to create bone structures. For instance, Javan et al. [
[27]- Javan R.
- Ellenbogen A.L.
- Greek N.
- Haji-Momenian S.
A prototype assembled 3D-printed phantom of the glenohumeral joint for fluoroscopic-guided shoulder arthrography.
] printed osseous structures from gypsum. Further, a pelvis phantom was developed by Niebuhr et al. [
[6]Niebuhr NI, Johnen W, Echner G, Runz A, Bach M, Stoll M, et al. The ADAM-pelvis phantom-an anthropomorphic, deformable and multimodal phantom for MRgRT. Phys Med Biol. 2019;64:04NT5 DOI: 10.1088/1361-6560/aafd5f.
], represented by a PMMA cylindrical phantom and bone structures, which were 3D plastic hollow models covered by strong gypsum in order to produce the needed HU numbers. The cortical bone in a paediatric pelvis phantom developed by Ali et al. [
[28]- Ali A.M.
- Hogg P.
- Johansen S.
- England A.
Construction and validation of a low cost paediatric pelvis phantom.
], and the radiodense of bony lumbosacrum [
[29]- Javan R.
- Bansal M.
- Tangestanipoor A.
A Prototype Hybrid Gypsum-Based 3-Dimensional Printed Training Model for Computed Tomography-Guided Spinal Pain Management.
] were also created by gypsum. However, the approach of utilizing gypsum has limitations concerning the reproducibility of the trabeculae. The other technology, the PolyJet 3D printing, combined with a new ultraviolet curable 3D printing resin with added ZrO
2 nanopowder to print a 3D model of a hand [
[30]- Shannon A.
- O'Connell A.
- O'Sullivan A.
- Byrne M.
- Clifford S.
- O'Sullivan K.J.
- et al.
A Radiopaque Nanoparticle-Based Ink Using PolyJet 3D Printing for Medical Applications.
]. Similarly, Hatamikia et al. [
[31]- Hatamikia S.
- Oberoi G.
- Unger E.
- Kronreif G.
- Kettenbach J.
- Buschmann M.
- et al.
Additively Manufactured Patient-Specific Anthropomorphic Thorax Phantom With Realistic Radiation Attenuation Properties. Front Bioeng.
] developed an epoxy-based amalgamate for printing the various skeletal structures in the thorax. Results showed that structural densities in the range of 42–705HU could be achieved. A similar approach is undertaken by [
[32]- Hazelaar C.
- van Eijnatten M.
- Dahele M.
- Wolff J.
- Forouzanfar T.
- Slotman B.
- et al.
Using 3D printing techniques to create an anthropomorphic thorax phantom for medical imaging purposes.
]. A limitation of these approaches is the cost of the printed anthropomorphic phantoms in comparison to these printed with FDM technology.
The developed method will be exploited in the manufacturing of a CT-compatible abdomen phantom to be used in several major applications such as optimizing CT protocols for colorectal cancer screening, development and tuning of 3D reconstruction techniques [
[33]- Frandon J.
- Akessoul P.
- Hamard A.
- Bezandry E.
- Loffroy R.
- Addala T.
- et al.
Comparison of acquisition and iterative reconstruction parameters in abdominal computed tomography-guided procedures: a phantom study.
], evaluating precisely the dose delivered in megavoltage X-rays [
[34]- Iramina H.
- Nakamura M.
- Mizowaki T.
Actual delivered dose calculation on intra-irradiation cone-beam computed tomography images: a phantom study.
] and checking the accuracy of radiotherapy treatment plans [
[11]- Goodall S.K.
- Rampant P.
- Smith W.
- Waterhouse D.
- Rowshanfarzad P.
- Ebert M.A.
Investigation of the effects of spinal surgical implants on radiotherapy dosimetry: A study of 3D printed phantoms.
]. Furthermore, this method is intended to assess the results of the surgical operation and follow-up of operated by robotic surgery patients, by printing phantoms before and after the operation and thus validating the robotic surgery results. The proposed new methodology opens possibilities for manufacturing of anthropomorphic phantoms, which can be successfully used for surgical planning, training and simulation of robotic needle-positioning and guidance systems for CT-guided puncture [
[35]- Kettenbach J.
- Kara L.
- Toporek G.
- Fuerst M.
- Kronreif G.
A robotic needle-positioning and guidance system for CT-guided puncture: Ex vivo results.
]. The abdomen phantoms, built on CT images would allow the creation of physical models, which accurately replicate the patient-specific anatomy. These physical anthropomorphic models are the perfect choice to initially set up and test the new technologies, such as the navigation system of the surgery robot, before conducting costly cadaveric studies and use of expensive virtual reality simulators.
Article info
Publication history
Published online: September 23, 2022
Accepted:
September 14,
2022
Received in revised form:
September 9,
2022
Received:
January 13,
2022
Copyright
© 2022 Associazione Italiana di Fisica Medica e Sanitaria. Published by Elsevier Ltd.