Midfacial Advancement Line-A Comparative Evaluation of a New Measurement Method in Orthognathic Surgery.

PURPOSE: In planning intraoral quadrangular Le Fort II osteotomy (IQLFIIO), simulation of hard and soft tissue changes will be important at the infraorbital, Le Fort I, and incisor level. The aim of our study was to evaluate a new method for visualization and quantification. MATERIALS AND METHODS: Three different methods of quantification were compared: the point-to-point (PTP) measurement, which has been viewed as the reference standard; part comparison analysis (PCA); and a new method, the midfacial advancement line (MFAL) measurement. We performed a measurement comparison study using the Bland-Altman method to measure agreement and enrolled patients with midfacial deficiency and Class III malocclusion who had undergone IQLFIIO. The primary predictor variable was the method of measurement. The primary outcome variable was the amount of midfacial advancement. We also investigated the time required, visualization quality, and interobserver agreement. RESULTS: The sample included 12 subjects with a mean age of 21.6 years; 7 patients were male. The PTP and MFAL showed no significant observer dependence. The advancement measured with PTP and MFAL showed no significant differences. However, the advancement measured using MFAL and PCA showed a significant difference. The highest rating of visualization was found for MFAL. The time requirements were similar for all 3 methods. CONCLUSIONS: Our results have shown that the MFAL is a suitable method for visualization and quantification of soft and hard tissue changes at all 3 face levels in 1 image. It could be a valuable tool for virtual planning of midfacial advancement surgery.

Design, Implementation, and Evaluation of a Head and Neck MRI RF Array Integrated with a 511 keV Transmission Source for Attenuation Correction in PET/MR.

The goal of this work is to further improve positron emission tomography (PET) attenuation correction and magnetic resonance (MR) sensitivity for head and neck applications of PET/MR. A dedicated 24-channel receive-only array, fully-integrated with a hydraulic system to move a transmission source helically around the patient and radiofrequency (RF) coil array, is designed, implemented, and evaluated. The device enables the calculation of attenuation coefficients from PET measurements at 511 keV including the RF coil and the particular patient. The RF coil design is PET-optimized by minimizing photon attenuation from coil components and housing. The functionality of the presented device is successfully demonstrated by calculating the attenuation map of a water bottle based on PET transmission measurements; results are in excellent agreement with reference values. It is shown that the device itself has marginal influence on the static magnetic field B0 and the radiofrequency transmit field B1 of the 3T PET/MR system. Furthermore, the developed RF array is shown to outperform a standard commercial 16-channel head and neck coil in terms of signal-to-noise ratio (SNR) and parallel imaging performance. In conclusion, the presented hardware enables accurate calculation of attenuation maps for PET/MR systems while improving the SNR of corresponding MR images in a single device without degrading the B0 and B1 homogeneity of the scanner.

Soft Tissue Changes in Patients Undergoing Intraoral Quadrangular Le Fort II Osteotomy Versus Conventional Le Fort I Osteotomy.

PURPOSE: The aim of this study was to evaluate soft tissue changes after intraoral quadrangular Le Fort II osteotomy (IOQLFII) and correlate those changes to underlying osseous changes. MATERIALS AND METHODS: Twenty-six non-growing patients with midfacial deficiency and Class III malocclusion were analyzed. A study group of 13 patients who underwent IOQLFII was compared with 13 patients who underwent conventional Le Fort I osteotomy (LFI). After fusion of pre- and postoperative computed tomograms, each patient's hard and corresponding soft tissue changes were measured. Measurement points were defined at 3 levels in the IOQLFII group (infraorbital rim [IR], sinus floor [SF], and lateral incisor tip [LI]) and at 2 levels in the LFI group (SF and LI). Linear models were created to test for correlations between hard and soft tissues. RESULTS: The slope (a1 coefficient) between anteroposterior hard and soft tissue changes was found to be highly significant at each measurement point for all groups. In the IOQLFII group, soft tissue advancement was 69% (confidence interval [CI], 62 to 77%) of the hard tissue advancement at the IR, 90% (CI, 84 to 96%) at the SF, and 73% (CI, 64 to 82%) at the LI. In the LFI group, the corresponding percentages were 90% (83 to 97%) at SF and 84% (77 to 90%) at LI. CONCLUSION: IOQLFII results in predictable correction of midfacial deficiency. At the IR, bony advancement always resulted in markedly less soft tissue advancement than at the SF level. These results indicate that the planned infraorbital advancement should not be too conservative because soft tissue changes are smaller in this region.

Equivalent thicknesses of beam hardening filters consisting of aluminium, copper, Al/Cu and Al/Gold combinations and plumbiferous acrylic for 40 to 150 kVp diagnostic spectra.

PURPOSE: Beam hardening filters used to reduce patient doses typically consist of aluminium, copper, or a combination of both. Optically transparent filters containing lead in plumbiferous acrylic became available. One vendor also uses a combination of aluminium and gold. Data is provided to compare filter thicknesses in terms of half-value layer (HVL) for clinically relevant kVp. METHODS: Equivalent filter thicknesses were defined by identical kVp and 1st HVL. Equivalent copper filter thicknesses were calculated for aluminium and typical filters found in radiographic and interventional systems. A verified semi-empirical spectrum calculation programme and National Institute of Standards and Technology (NIST) mass attenuation coefficients were applied. Lead acrylic filters were simulated by a two-component model of acrylic plus lead with mass thicknesses determined by matching 1 HVLs in Al at RQR5 using filter specifications. RESULTS: Coefficients are provided to convert mm Cu to mm Al and vice versa for tube potentials from 40 to 150 kVp. 1 mm Al corresponds to 27.8 ± 1 µm Cu over the entire energy range simulated. Using this simple model as opposed to simulations of all individual filters made from Al/Cu combinations (1 and 2 mm Al, 1 Al + 0.1 and 0.2 Cu, 1.5 Al plus 0.3 and 0.6 Cu, 2 Al plus 0.1 Cu) for the entire energy range results in differences in equivalent Cu thicknesses below 4 µm Cu (3 µm for 50-150 kVp). kVp dependence is larger for filters containing larger Z elements. 1 mm Al plus 10 µm gold used by Shimadzu corresponds to 75-80 µm Cu, depending on kVp; plumbiferous acrylic with nominal filtrations of 1 Al plus 0.1 Cu, and 1 Al plus 0.2 Cu corresponded to 124-132 µm, and 206-232 µm Cu, respectively. CONCLUSIONS: Experimental verification of the equivalence of aluminium and combined aluminium plus copper filters showed excellent agreement between calculated copper equivalent thickness and measurements with copper filters for clinical beams from 40 to 150 kVp.

Optical coherence tomography for planning of follicular unit extraction.

BACKGROUND: As the subcutaneous course of the hair follicle cannot be extrapolated from the angle of hair emergence, the follicular unit extraction (FUE) method bears the risk of follicle transsection potentially affecting the success of the whole intervention. Histological examination can provide this angle information but are too time consuming and invasive to be integrated in the workflow. OBJECTIVE: In this work, the authors present an effective and noninvasive way to determine these subcutaneous follicle angles. METHODS AND MATERIALS: For this purpose, the authors use optical coherence tomography to provide images of the first 2 to 3 mm of the surface of the skin. RESULTS: On 3 probands, they scanned 2 regions. We computed the number of hair follicles in a certain angular region around a mean angle, representing an FUE tool, and found a considerable diversity of the angular deviations. CONCLUSION: As this angular distribution massively influences the transsection rate, the authors suggest such a method to be considered as a future improvement of preoperative assessment.

Patient-specific polyetheretherketone facial implants in a computer-aided planning workflow.

PURPOSE: In the present study, we report an innovative workflow using polyetheretherketone (PEEK) patient-specific implants for esthetic corrections in the facial region through onlay grafting. The planning includes implant design according to virtual osteotomy and generation of a subtraction volume. The implant design was refined by stepwise changing the implant geometry according to soft tissue simulations. MATERIALS AND METHODS: One patient was scanned using computed tomography. PEEK implants were interactively designed and manufactured using rapid prototyping techniques. Positioning intraoperatively was assisted by computer-aided navigation. Two months after surgery, a 3-dimensional surface model of the patient's face was generated using photogrammetry. Finally, the Hausdorff distance calculation was used to quantify the overall error, encompassing the failures in soft tissue simulation and implantation. RESULTS: The implant positioning process during surgery was satisfactory. The simulated soft tissue surface and the photogrammetry scan of the patient showed a high correspondence, especially where the skin covered the implants. The mean total error (Hausdorff distance) was 0.81 ± 1.00 mm (median 0.48, interquartile range 1.11). The spatial deviation remained less than 0.7 mm for the vast majority of points. CONCLUSIONS: The proposed workflow provides a complete computer-aided design, computer-aided manufacturing, and computer-aided surgery chain for implant design, allowing for soft tissue simulation, fabrication of patient-specific implants, and image-guided surgery to position the implants. Much of the surgical complexity resulting from osteotomies of the zygoma, chin, or mandibular angle might be transferred into the planning phase of patient-specific implants.

A spiculated mass target model for clinical image quality control in digital mammography.

OBJECTIVES: Quality assurance of breast imaging has a long history of using test objects to optimize and follow up imaging devices. In particular, the evaluation of new techniques benefits from suitable test objects. The applicability of a phantom consisting of spiculated masses to assess image quality and its dependence on dose in flat field digital mammography (FFDM) and digital breast tomosynthesis systems (DBT) is investigated. METHODS: Two spiculated masses in five different sizes each were created from a database of clinical tumour models. The masses were produced using 3D printing and embedded into a cuboid phantom. Image quality is determined by the number of spicules identified by human observers. RESULTS: The results suggest that the effect of dose on spicule detection is limited especially in cases with smaller objects and probably hidden by the inter-reader variability. Here, an average relative inter-reader variation of the counted number of 31% was found (maximum 83%). The mean relative intra-reader variability was found to be 17%. In DBT, sufficiently good results were obtained only for the largest masses. CONCLUSIONS: It is possible to integrate spiculated masses into a cuboid phantom. It is easy to print and should allow a direct and prompt evaluation of the quality status of the device by counting visible spicules. Human readout presented the major uncertainty in this study, indicating that automated readout may improve the reproducibility and consistency of the results considerably. ADVANCES IN KNOWLEDGE: A cuboid phantom including clinical objects as spiculated lesion models for visual assessing the image quality in FFDM and DBT was developed and is introduced in this work. The evaluation of image quality works best with the two larger masses with 21 spicules.

Realistic 3D printed CT imaging tumor phantoms for validation of image processing algorithms.

Medical imaging phantoms are widely used for validation and verification of imaging systems and algorithms in surgical guidance and radiation oncology procedures. Especially, for the performance evaluation of new algorithms in the field of medical imaging, manufactured phantoms need to replicate specific properties of the human body, e.g., tissue morphology and radiological properties. Additive manufacturing (AM) technology provides an inexpensive opportunity for accurate anatomical replication with customization capabilities. In this study, we proposed a simple and cheap protocol using Fused Deposition Modeling (FDM) technology to manufacture realistic tumor phantoms based on the filament 3D printing technology. Tumor phantoms with both homogenous and heterogeneous radiodensity were fabricated. The radiodensity similarity between the printed tumor models and real tumor data from CT images of lung cancer patients was evaluated. Additionally, it was investigated whether a heterogeneity in the 3D printed tumor phantoms as observed in the tumor patient data had an influence on the validation of image registration algorithms. A radiodensity range between -217 to 226 HUs was achieved for 3D printed phantoms using different filament materials; this range of radiation attenuation is also observed in the human lung tumor tissue. The resulted HU range could serve as a lookup-table for researchers and phantom manufactures to create realistic CT tumor phantoms with the desired range of radiodensities. The 3D printed tumor phantoms also precisely replicated real lung tumor patient data regarding morphology and could also include life-like heterogeneity of the radiodensity inside the tumor models. An influence of the heterogeneity on accuracy and robustness of the image registration algorithms was not found.

Technical note: Realization and uncertainty analysis for an adjustable 3D structured breast phantom in digital breast tomosynthesis.

BACKGROUND: Projection imaging phantoms are often optimized for 2-dimensional image characteristics in homogeneous backgrounds. Therefore, evaluation of image quality in tomosynthesis (DBT) lacks accepted and established phantoms. PURPOSE: We describe a 3D breast phantom with a structured, variable background. The phantom is an adaptable and advanced version of the L1 phantom by Cockmartin et al. Phantom design and its use for quality assurance measurements for DBT devices are described. Four phantoms were compared to assess the objectivity. METHODS: The container size was increased to a diameter of 24 cm and a total height of 53.5 mm. Spiculated masses were replaced by five additional non-spiculated masses for higher granularity in threshold diameter resolution. These patterns are adjustable to the imaging device. The masses were printed in one session with a base layer using two-component 3D printing. New materials compared to the L1 phantom improved the attenuation difference between the lesion models and the background. Four phantoms were built and intra-human observer, inter-human observer and inter-phantom variations were determined. The latter assess the reproducibility of the phantom production. Coefficients of variance (V) were calculated for all three variations. RESULTS: The difference of the attenuation coefficients between the lesion models and the background was 0.20 cm-1 (with W/Al at 32 kV, equivalent to 19-20 keV effective energy) compared to 0.21 cm-1 for 50/50 glandular/adipose breast tissue and cancerous lesions. PMMA equivalent thickness of the phantom was 47.0 mm for the Siemens Mammomat Revelation. For the masses, the V i n t r a $V_{intra}$ for the intra-observer variation was 0.248, the averaged inter-observer variation, V ¯ i n t e r $\overline{V}_{inter}$ was 0.383. V p h a n t o m $V_{phantom}$ for phantom variance was 0.321. For the micro-calcifications, V i n t r a $V_{intra}$ was 0.0429, V ¯ i n t e r = $\overline{V}_{inter}=$ 0.0731 and V p h a n t o m = $V_{phantom}=$ 0.0759. CONCLUSIONS: Position, orientation and shape of the masses are reproducible and attenuation differences appropriate. The phantom presented proved to be a candidate test object for quality control.

Splint sterilization--a potential registration hazard in computer-assisted surgery.

PURPOSE: Registration of preoperative targeting information for the intraoperative situation is a crucial step in computer-assisted surgical interventions. Point-to-point registration using acrylic splints is among the most frequently used procedures. There are, however, no generally accepted recommendations for sterilization of the splint. An appropriate method for the thermolabile splint would be hydrogen peroxide-based plasma sterilization. This study evaluated the potential deformation of the splint undergoing such sterilization. Deformation was quantified using image-processing methods applied to computed tomographic (CT) volumes before and after sterilization. MATERIALS AND METHODS: An acrylic navigation splint was used as the study object. Eight metallic markers placed in the splint were used for registration. Six steel spheres in the mouthpiece were used as targets. Two CT volumes of the splint were acquired before and after 5 sterilization cycles using a hydrogen peroxide sterilizer. Point-to-point registration was applied, and fiducial and target registration errors were computed. Surfaces were extracted from CT scans and Hausdorff distances were derived. Effectiveness of sterilization was determined using Geobacillus stearothermophilus. RESULTS: Fiducial-based registration of CT scans before and after sterilization resulted in a mean fiducial registration error of 0.74 mm; the target registration error in the mouthpiece was 0.15 mm. The Hausdorff distance, describing the maximal deformation of the splint, was 2.51 mm. Ninety percent of point-surface distances were shorter than 0.61 mm, and 95% were shorter than 0.73 mm. No bacterial growth was found after the sterilization process. CONCLUSION: Hydrogen peroxide-based low-temperature plasma sterilization does not deform the splint, which is the base for correct computer-navigated surgery.

X-ray attenuation of bone, soft and adipose tissue in CT from 70 to 140 kV and comparison with 3D printable additive manufacturing materials.

Additive manufacturing and 3D printing are widely used in medical imaging to produce phantoms for image quality optimization, imaging protocol definition, comparison of image quality between different imaging systems, dosimetry, and quality control. Anthropomorphic phantoms mimic tissues and contrasts in real patients with regard to X-ray attenuation, as well as dependence on X-ray spectra. If used with different X-ray energies, or to optimize the spectrum for a certain procedure, the energy dependence of the attenuation must replicate the corresponding energy dependence of the tissues mimicked, or at least be similar. In the latter case the materials' Hounsfield values need to be known exactly to allow to correct contrast and contrast to noise ratios accordingly for different beam energies. Fresh bovine and porcine tissues including soft and adipose tissues, and hard tissues from soft spongious bone to cortical bone were scanned at different energies, and reference values of attenuation in Hounsfield units (HU) determined. Mathematical model equations describing CT number dependence on kV for bones of arbitrary density, and for adipose tissues are derived. These data can be used to select appropriate phantom constituents, compare CT values with arbitrary phantom materials, and calculate correction factors for phantoms consisting of materials with an energy dependence different to the tissues. Using data on a wide number of additive manufacturing and 3D printing materials, CT numbers and their energy dependence were compared to those of the tissues. Two commercially available printing filaments containing calcium carbonate powder imitate bone tissues with high accuracy at all kV values. Average adipose tissue can be duplicated by several off-the-shelf printing polymers. Since suitable printing materials typically exhibit a too high density for the desired attenuation of especially soft tissues, controlled density reduction by underfilling might improve tissue equivalence.

Technical note: A PET/MR coil with an integrated, orbiting 511 keV transmission source for PET/MR imaging validated in an animal study.

BACKGROUND: MR-based methods for attenuation correction (AC) in PET/MRI either neglect attenuation of bone, or use MR-signal derived information about bone, which leads to a bias in quantification of tracer uptake in PET. In a previous study, we presented a PET/MRI specific MR coil with an integrated transmission source (TX) system allowing for direct measurement of attenuation. In phantom measurements, this system successfully reproduced the linear attenuation coefficient of water. PURPOSE: The purpose of this study is to validate the TX system in a clinical setting using animals and to show its applicability compared to standard clinical methods. METHODS: As test subject, a 15-kg piglet was injected with 53 MBq of 18F-NaF. The µ-map obtained with the TX system and the reconstructed activity distribution were compared to four established AC methods: a Dixon sequence, an ultra-short echo time (UTE) sequence, a CT scan, and a 511 keV transmission scan using a Siemens ECAT EXACT HR+ as the reference. The PET/MRI measurements were performed on a Siemens Biograph mMR to obtain the µ-map using the TX system as well as the Dixon and UTE sequence directly followed by the CT and ECAT measurements. RESULTS: The reconstructed activity distribution using the TX system for AC showed similar results compared to the reference (<5% difference in hot regions) and outperformed the MR-based methods as implemented in the PET/MRI system (<10% difference in hot regions). However, the additional hardware of the TX system adds complexity to the acquisition process. CONCLUSION: Our porcine study demonstrates the feasibility of post-injection transmission scans using the developed TX system in a clinical setting. This makes it a useful tool for PET/MRI in cases where transmission information is needed for AC. Potential applications are studies using larger animals where state-of-the-art atlas-based or artificial intelligence AC methods are not available.

Factors for conversion between human and automatic read-outs of CDMAM images.

PURPOSE: According to the European protocol for the quality control of the physical and technical aspects of mammography screening (EPQCM) image quality of digital mammography devices has to be assessed using human evaluation of the CDMAM contrast-detail phantom. This is accomplished by the determination of threshold thicknesses of gold disks with different diameters (0.08-2 mm) and revealed to be very time consuming. Therefore a software solution based on a nonprewhitening matched filter (NPW) model was developed at the University of Nijmegen. Factors for the conversion from automatic to human readouts have been determined by Young et al. [Proc. SPIE 614206, 1-13 (2006) and Proc. SPIE 6913, 69131C1 (2008)] using a huge amount of data of both human and automatic readouts. These factors depend on the observer groups and are purely phenomenological. The authors present an alternative approach to determine the factors by using the Rose observer model. METHODS: Their method uses the Rose theory which gives a relationship between threshold contrast, diameter of the object and number of incident photons. To estimate the conversion factors for the five diameters from 0.2 to 0.5 mm they exposed with five different current-time products which resulted in 25 equations with five unknowns. RESULTS: The theoretical conversion factors (in dependence of the diameters) amounted to be 1.61 ± 0.02 (0.2 mm diameter), 1.67 ± 0.02 (0.25 mm), 1.85 ± 0.02 (0.31 mm), 2.09 ± 0.02 (0.4 mm), and 2.28 ± 0.02 (0.5 mm). The corresponding phenomenological factors found in literature are 1.74 (0.2 mm), 1.78 (0.25 mm), 1.83 (0.31 mm), 1.88 (0.4 mm), and 1.93 (0.5 mm). CONCLUSIONS: They transferred the problem of determining the factors to a well known observer model which has been examined for many years and is also well established. This method reveals to be reproduceable and produces factors comparable to the phenomenological ones.

Validation for 2D/3D registration. II: The comparison of intensity- and gradient-based merit functions using a new gold standard data set.

PURPOSE: A new gold standard data set for validation of 2D/3D registration based on a porcine cadaver head with attached fiducial markers was presented in the first part of this article. The advantage of this new phantom is the large amount of soft tissue, which simulates realistic conditions for registration. This article tests the performance of intensity- and gradient-based algorithms for 2D/3D registration using the new phantom data set. METHODS: Intensity-based methods with four merit functions, namely, cross correlation, rank correlation, correlation ratio, and mutual information (MI), and two gradient-based algorithms, the backprojection gradient-based (BGB) registration method and the reconstruction gradient-based (RGB) registration method, were compared. Four volumes consisting of CBCT with two fields of view, 64 slice multidetector CT, and magnetic resonance-T1 weighted images were registered to a pair of kV x-ray images and a pair of MV images. A standardized evaluation methodology was employed. Targets were evenly spread over the volumes and 250 starting positions of the 3D volumes with initial displacements of up to 25 mm from the gold standard position were calculated. After the registration, the displacement from the gold standard was retrieved and the root mean square (RMS), mean, and standard deviation mean target registration errors (mTREs) over 250 registrations were derived. Additionally, the following merit properties were computed: Accuracy, capture range, number of minima, risk of nonconvergence, and distinctiveness of optimum for better comparison of the robustness of each merit. RESULTS: Among the merit functions used for the intensity-based method, MI reached the best accuracy with an RMS mTRE down to 1.30 mm. Furthermore, it was the only merit function that could accurately register the CT to the kV x rays with the presence of tissue deformation. As for the gradient-based methods, BGB and RGB methods achieved subvoxel accuracy (RMS mTRE down to 0.56 and 0.70 mm, respectively). Overall, gradient-based similarity measures were found to be substantially more accurate than intensity-based methods and could cope with soft tissue deformation and enabled also accurate registrations of the MR-T1 volume to the kV x-ray image. CONCLUSIONS: In this article, the authors demonstrate the usefulness of a new phantom image data set for the evaluation of 2D/3D registration methods, which featured soft tissue deformation. The author's evaluation shows that gradient-based methods are more accurate than intensity-based methods, especially when soft tissue deformation is present. However, the current nonoptimized implementations make them prohibitively slow for practical applications. On the other hand, the speed of the intensity-based method renders these more suitable for clinical use, while the accuracy is still competitive.

3D ultrasound guided navigation system with hybrid image fusion.

A prototype of a navigation system to fuse two image modalities is presented. The standard inter-modality registration is replaced with a tracker-based image registration of calibrated imaging devices. Intra-procedure transrectal US (TRUS) images were merged with pre-procedure magnetic resonance (MR) images for prostate biopsy. The registration between MR and TRUS images was performed by an additional abdominal 3D-US (ab-3D-US), which enables replacing the inter-modal MR/TRUS registration by an intra-modal ab-3D-US/3D-TRUS registration. Calibration procedures were carried out using an optical tracking system (OTS) for the pre-procedure image fusion of the ab-3D-US with the MR. Inter-modal ab-3D-US/MR image fusion was evaluated using a multi-cone phantom for the target registration error (TRE) and a prostate phantom for the Dice score and the Hausdorff distance of lesions . Finally, the pre-procedure ab- 3D-US was registered with the TRUS images and the errors for the transformation from the MR to the TRUS were determined. The TRE of the ab-3D-US/MR image registration was 1.81 mm. The Dice-score and the Hausdorff distance for ab-3D-US and MR were found to be 0.67 and 3.19 mm. The Dice score and the Hausdorff distance for TRUS and MR were 0.67 and 3.18 mm. The hybrid navigation system showed sufficient accuracy for fusion guided biopsy procedures with prostate phantoms. The system might provide intra-procedure fusion for most US-guided biopsy and ablation interventions.

Toward on-the-fly trajectory optimization for C-arm CBCT under strong kinematic constraints.

Cone beam computed tomography (CBCT) has become a vital tool in interventional radiology. Usually, a circular source-detector trajectory is used to acquire a three-dimensional (3D) image. Kinematic constraints due to the patient size or additional medical equipment often cause collisions with the imager while performing a full circular rotation. In a previous study, we developed a framework to design collision-free, patient-specific trajectories for the cases in which circular CBCT is not feasible. Our proposed trajectories included enough information to appropriately reconstruct a particular volume of interest (VOI), but the constraints had to be defined before the intervention. As most collisions are unpredictable, performing an on-the-fly trajectory optimization is desirable. In this study, we propose a search strategy that explores a set of trajectories that cover the whole collision-free area and subsequently performs a search locally in the areas with the highest image quality. Selecting the best trajectories is performed using simulations on a prior diagnostic CT volume which serves as a digital phantom for simulations. In our simulations, the Feature SIMilarity Index (FSIM) is used as the objective function to evaluate the imaging quality provided by different trajectories. We investigated the performance of our methods using three different anatomical targets inside the Alderson-Rando phantom. We used FSIM and Universal Quality Image (UQI) to evaluate the final reconstruction results. Our experiments showed that our proposed trajectories could achieve a comparable image quality in the VOI compared to the standard C-arm circular CBCT. We achieved a relative deviation less than 10% for both FSIM and UQI metrics between the reconstructed images from the optimized trajectories and the standard C-arm CBCT for all three targets. The whole trajectory optimization took approximately three to four minutes.

Evaluation of 3D ultrasound for image guidance.

PURPOSE: In this paper we compared two different 3D ultrasound (US) modes (3D free-hand mode and 3D wobbler mode) to see which is more suitable to perform the 3D-US/3D-US registration for clinical guidance applications. The typical errors with respect to their impact on the final localization error were evaluated step by step. METHODS: Multi-point target and Hand-eye calibration methods were used for 3D US calibration together with a newly designed multi-cone phantom. Pointer based and image based methods were used for 2D US calibration. The calibration target error was computed by using a different multi-cone phantom. An egg-shaped phantom was used as ground truth to compare distortions for both 3D modes along with the measurements of the volume. Finally, we compared 3D ultrasound images acquired by 3D wobbler mode and 3D free-hand mode with respect to their 3D-US/3D-US registration accuracy using both, phantom and patient data. A theoretical step by step error analysis was performed and compared to empirical data. RESULTS: Target registration errors based on the calibration with the 3D Multi-point and 2D pointer/image method have been found to be comparable (∼1mm). They both outperformed the 3D Hand-eye method (error >2mm). Volume measurements with the 3D free-hand mode were closest to the ground truth (around 6% error compared to 9% with the 3D wobbler mode). Additional scans on phantoms showed a 3D-US/3D-US registration error below 1 mm for both, the 3D free-hand mode and the 3D wobbler mode, respectively. Results with patient data showed greater error with the 3D free-hand mode (6.50mm - 13.37mm) than with the 3D wobbler mode (2.99 ± 1.54 mm). All the measured errors were found to be in accordance to their theoretical upper bounds. CONCLUSION: While both 3D volume methods showed comparable results with respect to 3D-US/3D-US registration for phantom images, for patient data registrations the 3D wobbler mode is superior to the 3D free-hand mode. The effect of all error sources could be estimated by theoretical derivations.

Additively Manufactured Patient-Specific Anthropomorphic Thorax Phantom With Realistic Radiation Attenuation Properties.

Conventional medical imaging phantoms are limited by simplified geometry and radiographic skeletal homogeneity, which confines their usability for image quality assessment and radiation dosimetry. These challenges can be addressed by additive manufacturing technology, colloquially called 3D printing, which provides accurate anatomical replication and flexibility in material manipulation. In this study, we used Computed Tomography (CT)-based modified PolyJetTM 3D printing technology to print a hollow thorax phantom simulating skeletal morphology of the patient. To achieve realistic heterogenous skeletal radiation attenuation, we developed a novel radiopaque amalgamate constituting of epoxy, polypropylene and bone meal powder in twelve different ratios. We performed CT analysis for quantification of material radiodensity (in Hounsfield Units, HU) and for identification of specific compositions corresponding to the various skeletal structures in the thorax. We filled the skeletal structures with their respective radiopaque amalgamates. The phantom and isolated 3D printed rib specimens were rescanned by CT for reproducibility tests regarding verification of radiodensity and geometry. Our results showed that structural densities in the range of 42-705HU could be achieved. The radiodensity of the reconstructed phantom was comparable to the three skeletal structures investigated in a real patient thorax CT: ribs, ventral vertebral body and dorsal vertebral body. Reproducibility tests based on physical dimensional comparison between the patient and phantom CT-based segmentation displayed 97% of overlap in the range of 0.00-4.57 mm embracing the anatomical accuracy. Thus, the additively manufactured anthropomorphic thorax phantom opens new vistas for imaging- and radiation-based patient care in precision medicine.

Optimization for customized trajectories in cone beam computed tomography.

PURPOSE: We developed a target-based cone beam computed tomography (CBCT) imaging framework for optimizing an unconstrained three dimensional (3D) source-detector trajectory by incorporating prior image information. Our main aim is to enable a CBCT system to provide topical information about the target using a limited angle noncircular scan orbit with a minimal number of projections. Such a customized trajectory should include enough information to sufficiently reconstruct a particular volume of interest (VOI) under kinematic constraints, which may result from the patient size or additional surgical or radiation therapy-related equipment. METHODS: A patient-specific model from a prior diagnostic computed tomography (CT) volume is used as a digital phantom for CBCT trajectory simulations. Selection of the best projection views is accomplished through maximizing an objective function fed by the imaging quality provided by different x-ray positions on the digital phantom data. The final optimized trajectory includes a limited angular range and a minimal number of projections which can be applied to a C-arm device capable of general source-detector positioning. The performance of the proposed framework is investigated in experiments involving an in-house-built box phantom including spherical targets as well as an Alderson-Rando head phantom. In order to quantify the image quality of the reconstructed image, we use the average full-width-half-maximum (FWHMavg ) for the spherical target and feature similarity index (FSIM), universal quality index (UQI), and contrast-to-noise ratio (CNR) for an anatomical target. RESULTS: Our experiments based on both the box and head phantom showed that optimized trajectories could achieve a comparable image quality in the VOI with respect to the standard C-arm circular CBCT while using approximately one quarter of projections. We achieved a relative deviation <7% for FWHMavg between the reconstructed images from the optimized trajectories and the standard C-arm CBCT for all spherical targets. Furthermore, for the anatomical target, the relative deviation of FSIM, UQI, and CNR between the reconstructed image related to the proposed trajectory and the standard C-arm circular CBCT was found to be 5.06%, 6.89%, and 8.64%, respectively. We also compared our proposed trajectories to circular trajectories with equivalent angular sampling as the optimized trajectories. Our results show that optimized trajectories can outperform simple partial circular trajectories in the VOI in term of image quality. Typically, an angular range between 116° and 152° was used for the optimized trajectories. CONCLUSION: We demonstrated that applying limited angle noncircular trajectories with optimized orientations in 3D space can provide a suitable image quality for particular image targets and has a potential for limited angle and low-dose CBCT-based interventions under strong spatial constraints.

A comparative study on manual and automatic slice-to-volume registration of CT images.

In order to assess the clinical relevance of a slice-to-volume registration algorithm, this technique was compared to manual registration. Reformatted images obtained from a diagnostic CT examination of the lower abdomen were reviewed and manually registered by 41 individuals. The results were refined by the algorithm. Furthermore, a fully automatic registration of the single slices to the whole CT examination, without manual initialization, was also performed. The manual registration error for rotation and translation was found to be 2.7+/-2.8 degrees and 4.0+/-2.5 mm. The automated registration algorithm significantly reduced the registration error to 1.6+/-2.6 degrees and 1.3+/-1.6 mm (p = 0.01). In 3 of 41 (7.3%) registration cases, the automated registration algorithm failed completely. On average, the time required for manual registration was 213+/-197 s; automatic registration took 82+/-15 s. Registration was also performed without any human interaction. The resulting registration error of the algorithm without manual pre-registration was found to be 2.9+/-2.9 degrees and 1.1+/-0.2 mm. Here, a registration took 91+/-6 s, on average. Overall, the automated registration algorithm improved the accuracy of manual registration by 59% in rotation and 325% in translation. The absolute values are well within a clinically relevant range.