1
|
Stengl C, Panow K, Arbes E, Muñoz ID, Christensen JB, Neelsen C, Dinkel F, Weidner A, Runz A, Johnen W, Liermann J, Echner G, Vedelago J, Jäkel O. A phantom to simulate organ motion and its effect on dose distribution in carbon ion therapy for pancreatic cancer. Phys Med Biol 2023; 68:245013. [PMID: 37918022 DOI: 10.1088/1361-6560/ad0902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Accepted: 11/02/2023] [Indexed: 11/04/2023]
Abstract
Objective. Carbon ion radiotherapy is a promising radiation technique for malignancies like pancreatic cancer. However, organs' motion imposes challenges for achieving homogeneous dose delivery. In this study, an anthropomorphicPancreasPhantom forIon-beamTherapy (PPIeT) was developed to simulate breathing and gastrointestinal motion during radiotherapy.Approach. The developed phantom contains a pancreas, two kidneys, a duodenum, a spine and a spinal cord. The shell of the organs was 3D printed and filled with agarose-based mixtures. Hounsfield Units (HU) of PPIeTs' organs were measured by CT. The pancreas motion amplitude in cranial-caudal (CC) direction was evaluated from patients' 4D CT data. Motions within the obtained range were simulated and analyzed in PPIeT using MRI. Additionally, GI motion was mimicked by changing the volume of the duodenum and quantified by MRI. A patient-like treatment plan was calculated for carbon ions, and the phantom was irradiated in a static and moving condition. Dose measurements in the organs were performed using an ionization chamber and dosimetric films.Main results. PPIeT presented tissue equivalent HU and reproducible breathing-induced CC displacements of the pancreas between (3.98 ± 0.36) mm and a maximum of (18.19 ± 0.44) mm. The observed maximum change in distance of (14.28 ± 0.12) mm between pancreas and duodenum was consistent with findings in patients. Carbon ion irradiation revealed homogenous coverage of the virtual tumor at the pancreas in static condition with a 1% deviation from the treatment plan. Instead, the dose delivery during motion with the maximum amplitude yielded an underdosage of 21% at the target and an increased uncertainty by two orders of magnitude.Significance. A dedicated phantom was designed and developed for breathing motion assessment of dose deposition during carbon ion radiotherapy. PPIeT is a unique tool for dose verification in the pancreas and its organs at risk during end-to-end tests.
Collapse
Affiliation(s)
- Christina Stengl
- Medical Faculty Heidelberg, Heidelberg University, Im Neuenheimer Feld 672, Heidelberg D-69120, Germany
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| | - Kathrin Panow
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| | - Eric Arbes
- Biomedical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Department for Physics and Astronomy, Heidelberg University, Im Neuenheimer Feld 226, Heidelberg D-69120, Germany
| | - Iván D Muñoz
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
- Department for Physics and Astronomy, Heidelberg University, Im Neuenheimer Feld 226, Heidelberg D-69120, Germany
| | - Jeppe B Christensen
- Department of Radiation Safety and Security, Paul Scherrer Institute (PSI), Forschungsstrasse 111, Villigen PSI 5232, Switzerland
| | - Christian Neelsen
- Department of Radiology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Department of Nuclear Medicine, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Department of Radiology, Charité-Universitätsmedizin Berlin, Charitéplatz 1, Berlin D-10117, Germany
| | - Fabian Dinkel
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| | - Artur Weidner
- Medical Faculty Heidelberg, Heidelberg University, Im Neuenheimer Feld 672, Heidelberg D-69120, Germany
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| | - Armin Runz
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| | - Wibke Johnen
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| | - Jakob Liermann
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
- Department of Radiation Oncology, Heidelberg University Hospital (UKHD), Im Neuenheimer Feld 400, Heidelberg D-69120, Germany
- Heidelberg Ion-Beam Therapy Center (HIT), Department of Radiation Oncology, Heidelberg University Hospital (UKHD), Im Neuenheimer Feld 450, Heidelberg D-69120, Germany
- National Center for Tumor Diseases (NCT), Im Neuenheimer Feld 460, D-69120 Heidelberg, Germany
| | - Gernot Echner
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
| | - José Vedelago
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
- Department of Radiation Oncology, Heidelberg University Hospital (UKHD), Im Neuenheimer Feld 400, Heidelberg D-69120, Germany
| | - Oliver Jäkel
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
- Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology (NCRO), Heidelberg, Germany
- Heidelberg Ion-Beam Therapy Center (HIT), Department of Radiation Oncology, Heidelberg University Hospital (UKHD), Im Neuenheimer Feld 450, Heidelberg D-69120, Germany
| |
Collapse
|
2
|
Fahad HM, Dorsch S, Zaiss M, Karger CP. Multi-parametric optimization of magnetic resonance imaging sequences for magnetic resonance-guided radiotherapy. Phys Imaging Radiat Oncol 2023; 28:100497. [PMID: 37869476 PMCID: PMC10585385 DOI: 10.1016/j.phro.2023.100497] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2023] [Revised: 09/15/2023] [Accepted: 09/29/2023] [Indexed: 10/24/2023] Open
Abstract
Background and purpose Magnetic Resonance Imaging (MRI) is widely used in oncology for tumor staging, treatment response assessment, and radiation therapy (RT) planning. This study proposes a framework for automatic optimization of MRI sequences based on pulse sequence parameter sets (SPS) that are directly applied on the scanner, for application in RT planning. Materials and methods A phantom with seven in-house fabricated contrasts was used for measurements. The proposed framework employed a derivative-free optimization algorithm to repeatedly update and execute a parametrized sequence on the MR scanner to acquire new data. In each iteration, the mean-square error was calculated based on the clinical application. Two clinically relevant optimization goals were pursued: achieving the same signal and therefore contrast as in a target image, and maximizing the signal difference (contrast) between specified tissue types. The framework was evaluated using two optimization methods: a covariance matrix adaptation evolution strategy (CMA-ES) and a genetic algorithm (GA). Results The obtained results demonstrated the potential of the proposed framework for automatic optimization of MRI sequences. Both CMA-ES and GA methods showed promising results in achieving the two optimization goals, however, CMA-ES converged much faster as compared to GA. Conclusions The proposed framework enables for automatic optimization of MRI sequences based on SPS that are directly applied on the scanner and it may be used to enhance the quality of MRI images for dedicated applications in MR-guided RT.
Collapse
Affiliation(s)
- Hafiz Muhammad Fahad
- German Cancer Research Center (DKFZ), Medical Physics in Radiation Oncology, Heidelberg, Germany
- University of Heidelberg, Faculty of Medicine, Heidelberg, Germany
- National Center for Radiation Research in Oncology (NCRO), Heidelberg Institute for Radiation Oncology (HIRO), Heidelberg, Germany
| | - Stefan Dorsch
- German Cancer Research Center (DKFZ), Medical Physics in Radiation Oncology, Heidelberg, Germany
- National Center for Radiation Research in Oncology (NCRO), Heidelberg Institute for Radiation Oncology (HIRO), Heidelberg, Germany
- Department of Radiation Oncology, Heidelberg University Hospital, Heidelberg, Germany
| | - Moritz Zaiss
- Friedrich-Alexander Universität Erlangen-Nürnberg (FAU), Institute of Neuroradiology, University Hospital Erlangen, Erlangen, Germany
- Magnetic Resonance Center, Max- Planck Institute for Biological Cyberrnetics, Tübingen, Germany
| | - Christian P. Karger
- German Cancer Research Center (DKFZ), Medical Physics in Radiation Oncology, Heidelberg, Germany
- National Center for Radiation Research in Oncology (NCRO), Heidelberg Institute for Radiation Oncology (HIRO), Heidelberg, Germany
| |
Collapse
|
3
|
Jin H, Lee SY, An HJ, Choi CH, Chie EK, Wu HG, Park JM, Park S, Kim JI. Development of an anthropomorphic multimodality pelvic phantom for quantitative evaluation of a deep-learning-based synthetic computed tomography generation technique. J Appl Clin Med Phys 2022; 23:e13644. [PMID: 35579090 PMCID: PMC9359037 DOI: 10.1002/acm2.13644] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2021] [Revised: 04/06/2022] [Accepted: 04/28/2022] [Indexed: 11/11/2022] Open
Abstract
PURPOSE The objective of this study was to fabricate an anthropomorphic multimodality pelvic phantom to evaluate a deep-learning-based synthetic computed tomography (CT) algorithm for magnetic resonance (MR)-only radiotherapy. METHODS Polyurethane-based and silicone-based materials with various silicone oil concentrations were scanned using 0.35 T MR and CT scanner to determine the tissue surrogate. Five tissue surrogates were determined by comparing the organ intensity with patient CT and MR images. Patient-specific organ modeling for three-dimensional printing was performed by manually delineating the structures of interest. The phantom was finally fabricated by casting materials for each structure. For the quantitative evaluation, the mean and standard deviations were measured within the regions of interest on the MR, simulation CT (CTsim ), and synthetic CT (CTsyn ) images. Intensity-modulated radiation therapy plans were generated to assess the impact of different electron density assignments on plan quality using CTsim and CTsyn . The dose calculation accuracy was investigated in terms of gamma analysis and dose-volume histogram parameters. RESULTS For the prostate site, the mean MR intensities for the patient and phantom were 78.1 ± 13.8 and 86.5 ± 19.3, respectively. The mean intensity of the synthetic image was 30.9 Hounsfield unit (HU), which was comparable to that of the real CT phantom image. The original and synthetic CT intensities of the fat tissue in the phantom were -105.8 ± 4.9 HU and -107.8 ± 7.8 HU, respectively. For the target volume, the difference in D95% was 0.32 Gy using CTsyn with respect to CTsim values. The V65Gy values for the bladder in the plans using CTsim and CTsyn were 0.31% and 0.15%, respectively. CONCLUSION This work demonstrated that the anthropomorphic phantom was physiologically and geometrically similar to the patient organs and was employed to quantitatively evaluate the deep-learning-based synthetic CT algorithm.
Collapse
Affiliation(s)
- Hyeongmin Jin
- Department of Radiation Oncology, Seoul National University Hospital, Seoul, Republic of Korea
| | - Sung Young Lee
- Department of Radiation Oncology, Seoul National University Hospital, Seoul, Republic of Korea
| | - Hyun Joon An
- Department of Radiation Oncology, Seoul National University Hospital, Seoul, Republic of Korea.,Institute of Radiation Medicine, Seoul National University Medical Research Center, Seoul, Republic of Korea.,Biomedical Research Institute, Seoul National University Hospital, Seoul, Republic of Korea
| | - Chang Heon Choi
- Department of Radiation Oncology, Seoul National University Hospital, Seoul, Republic of Korea.,Institute of Radiation Medicine, Seoul National University Medical Research Center, Seoul, Republic of Korea.,Biomedical Research Institute, Seoul National University Hospital, Seoul, Republic of Korea
| | - Eui Kyu Chie
- Department of Radiation Oncology, Seoul National University Hospital, Seoul, Republic of Korea.,Institute of Radiation Medicine, Seoul National University Medical Research Center, Seoul, Republic of Korea.,Biomedical Research Institute, Seoul National University Hospital, Seoul, Republic of Korea.,Department of Radiation Oncology, Seoul National University College of Medicine, Seoul, Republic of Korea
| | - Hong-Gyun Wu
- Department of Radiation Oncology, Seoul National University Hospital, Seoul, Republic of Korea.,Institute of Radiation Medicine, Seoul National University Medical Research Center, Seoul, Republic of Korea.,Biomedical Research Institute, Seoul National University Hospital, Seoul, Republic of Korea.,Department of Radiation Oncology, Seoul National University College of Medicine, Seoul, Republic of Korea
| | - Jong Min Park
- Department of Radiation Oncology, Seoul National University Hospital, Seoul, Republic of Korea.,Institute of Radiation Medicine, Seoul National University Medical Research Center, Seoul, Republic of Korea.,Biomedical Research Institute, Seoul National University Hospital, Seoul, Republic of Korea.,Department of Radiation Oncology, Seoul National University College of Medicine, Seoul, Republic of Korea.,Robotics Research Laboratory for Extreme Environments, Advanced Institute of Convergence Technology, Suwon, Republic of Korea
| | - Sukwon Park
- Department of Radiation Oncology, Myongji Hospital, Goyang-si, Gyeonggi-do, Republic of Korea
| | - Jung-In Kim
- Department of Radiation Oncology, Seoul National University Hospital, Seoul, Republic of Korea.,Institute of Radiation Medicine, Seoul National University Medical Research Center, Seoul, Republic of Korea.,Biomedical Research Institute, Seoul National University Hospital, Seoul, Republic of Korea
| |
Collapse
|
4
|
Weidner A, Stengl C, Dinkel F, Dorsch S, Murillo C, Seeber S, Gnirs R, Runz A, Echner G, Karger CP, Jäkel O. An abdominal phantom with anthropomorphic organ motion and multimodal imaging contrast for MR-guided radiotherapy. Phys Med Biol 2022; 67. [PMID: 35081516 DOI: 10.1088/1361-6560/ac4ef8] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Accepted: 01/26/2022] [Indexed: 11/12/2022]
Abstract
Purpose. Improvements in image-guided radiotherapy (IGRT) enable accurate and precise treatment of moving tumors in the abdomen while simultaneously sparing healthy tissue. However, the lack of validation tools for newly developed MR-guided radiotherapy hybrid devices such as the MR-Linac is an open issue. This study presents a custom developed abdominal phantom with respiratory organ motion and multimodal imaging contrast to perform end-to-end tests for IGRT treatment planning scenarios.Methods. The abdominal phantom contains deformable and anatomically shaped liver and kidney models made of Ni-DTPA and KCl-doped agarose mixtures that can be reproducibly positioned within the phantom. Organ models are wrapped in foil to avoid ion exchange with the surrounding agarose and to provide stable T1 and T2 relaxation times as well as HU numbers. Breathing motion is realized by a diaphragm connected to an actuator that is hydraulically controlled via a programmable logic controller. With this system, artificial and patient-specific breathing patterns can be carried out. In 1.5 T magnetic resonance imaging (MRI), diaphragm, liver and kidney motion was measured and compared to the breathing motion of a healthy male volunteer for different breathing amplitudes including shallow, normal and deep breathing.Results. The constructed abdominal phantom demonstrated organ-equivalent intensity values in CT as well as in MRI. T1-weighted (T1w) and T2-weighted (T2w) relaxation times for 1.5 T and CT numbers were 552.9 ms, 48.2 ms and 48.8 HU (liver) as well as 950.42 ms, 79 ms and 28.2 HU (kidney), respectively. These values were stable for more than six months. Extracted breathing motion from a healthy volunteer revealed a liver to diaphragm motion ratio (LDMR) of 64.4% and a kidney to diaphragm motion ratio (KDMR) of 30.7%. Well-comparable values were obtained for the phantom (LDMR: 65.5%, KDMR: 27.5%).Conclusions. The abdominal phantom demonstrated anthropomorphic T1 and T2 relaxation times as well as HU numbers and physiological motion pattern in MRI and CT. This allows for wide use in the validation of IGRT including MRgRT.
Collapse
Affiliation(s)
- Artur Weidner
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany.,Faculty of Medicine, University of Heidelberg, Im Neuenheimer Feld 672, Heidelberg D-69120, Germany.,Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology, Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Christina Stengl
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany.,Faculty of Medicine, University of Heidelberg, Im Neuenheimer Feld 672, Heidelberg D-69120, Germany.,Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology, Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Fabian Dinkel
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany.,Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology, Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Stefan Dorsch
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany.,Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology, Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Carlos Murillo
- Division of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Steffen Seeber
- Division of Medical Physics in Radiology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Regula Gnirs
- Division of Radiology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Armin Runz
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany.,Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology, Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Gernot Echner
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany.,Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology, Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Christian P Karger
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany.,Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology, Im Neuenheimer Feld 280, Heidelberg D-69120, Germany
| | - Oliver Jäkel
- Division of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, Heidelberg D-69120, Germany.,Heidelberg Institute for Radiation Oncology (HIRO), National Center for Radiation Research in Oncology, Im Neuenheimer Feld 280, Heidelberg D-69120, Germany.,Heidelberg Ion-Beam Therapy Center (HIT), Im Neuenheimer Feld 450, Heidelberg D-69120, Germany
| |
Collapse
|
5
|
Ma D, Gao R, Li M, Qiu J. Mechanical and medical imaging properties of 3D-printed materials as tissue equivalent materials. J Appl Clin Med Phys 2021; 23:e13495. [PMID: 34878729 PMCID: PMC8833282 DOI: 10.1002/acm2.13495] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Revised: 07/12/2021] [Accepted: 11/18/2021] [Indexed: 12/21/2022] Open
Abstract
Three materials of polylactic acid (PLA), polyamide 12 (PA12), and light curing resin (LCR) were used to construct phantom using 3D printing technology. The mechanical and medical imaging properties of the three materials, such as elastic modulus, density, effective atomic number, X‐ray attenuation coefficient, computed tomography (CT) number, and acoustic properties, were investigated. The results showed that the elastic modulus for PLA was 1.98 × 103 MPa, for PA12 was 848 MPa, for LCR was 1.18×103 MPa, and that of three materials was close to some bones. In the range of 40∼120 kV, the X‐ray attenuation coefficient of three materials decreased with increasing tube voltage. The CT number for PLA, PA12, and LCR was 144, −88, and 312 Hounsfield units at 120 kV tube voltage, respectively. The density and the effective atomic number product (ρ*Zeff) were computed from three materials and decreased in the order of LCR, PLA, and PA12. The acoustic properties of materials were also studied. The speeds of sound of three materials were similar with those of some soft tissues.
Collapse
Affiliation(s)
- Depeng Ma
- Medical Engineering and Technology Center, Shandong First Medical University and Shandong Academy of Medical Sciences, Taian, P. R. China.,Qingdao 3E3D Tech. Co. Ltd., Qingdao, P. R. China
| | - Ronghui Gao
- Health Care Department, Taishan Sanatorium of Shandong, Province, Taian, P. R. China
| | - Minghui Li
- Medical Engineering and Technology Center, Shandong First Medical University and Shandong Academy of Medical Sciences, Taian, P. R. China
| | - Jianfeng Qiu
- Medical Engineering and Technology Center, Shandong First Medical University and Shandong Academy of Medical Sciences, Taian, P. R. China.,Medical Science and Technology Innovation Center, Shandong First Medical University and Shandong Academy of Medical Sciences, Ji'nan, P. R. China.,Qingdao 3E3D Tech. Co. Ltd., Qingdao, P. R. China
| |
Collapse
|