Exploring Variability in CT Characterization of Tumors: A Preliminary Phantom Study

PURPOSE: To explore the effects of computed tomography (CT) slice thickness and reconstruction algorithm on quantification of image features to characterize tumors using a chest phantom. MATERIALS AND METHODS: Twenty-two phantom lesions of known sizes (10 and 20 mm), shapes (spherical, elliptical, lobulated, and spiculated), and densities [-630, -10, and +100 Hounsfield Unit (HU)] were inserted into an anthropomorphic thorax phantom and scanned three times with relocations. The raw data were reconstructed using six imaging settings, i.e., a combination of three slice thicknesses of 1.25, 2.5, and 5 mm and two reconstruction kernels of lung and standard. Lesions were segmented and 14 image features representing lesion size, shape, and texture were calculated. Differences in the measured image features due to slice thickness and reconstruction algorithm were compared using linear regression method by adjusting three confounding variables (size, density, and shape). RESULTS: All 14 features were significantly different between 1.25 and 5 mm slice images. The 1.25 and 2.5 mm slice thicknesses were better than 5 mm for volume, density mean, density SD gray-level co-occurrence matrix (GLCM) energy and homogeneity. As for the reconstruction algorithm, there was no significant difference in uni-dimension, volume, shape index 9, and compactness. Lung reconstruction was better for density mean, whereas standard reconstruction was better for density SD. CONCLUSIONS: CT slice thickness and reconstruction algorithm can significantly affect the quantification of image features. Thinner (1.25 and 2.5 mm) and thicker (5 mm) slice images should not be used interchangeably. Sharper and smoother reconstructions significantly affect the density-based features.     

Designing filters for fast-known NcRNA identification

Detecting members of known noncoding RNA (ncRNA) families in genomic DNA is an important part of sequence annotation. However, the most widely used tool for modeling ncRNA families, the covariance model (CM), incurs a high-computational cost when used for genome-wide search. This cost can be reduced by using a filter to exclude sequences that are unlikely to contain the ncRNA of interest, applying the CM only where it is likely to match strongly. Despite recent advances, designing an efficient filter that can detect ncRNA instances lacking strong conservation while excluding most irrelevant sequences remains challenging. In this work, we design three types of filters based on multiple secondary structure profiles (SSPs). An SSP augments a regular profile (i.e., a position weight matrix) with secondary structure information but can still be efficiently scanned against long sequences. Multi-SSPbased filters combine evidence from multiple SSP matches and can achieve high sensitivity and specificity. Our SSP-based filters are extensively tested in BRAliBase III data set, Rfam 9.0, and a published soil metagenomic data set. In addition, we compare the SSPbased filters with several other ncRNA search tools including Infernal (with profile HMMs as filters), ERPIN, and tRNAscan-SE. Our experiments demonstrate that carefully designed SSP filters can achieve significant speedup over unfiltered CM search while maintaining high sensitivity for various ncRNA families. The designed filters and filter-scanning programs are available at our website: www.cse.msu.edu/~yannisun/ssp/.     

CT dosimetry considerations in non typical conditions: The effect of scan field of view and table height selection

PurposeTo experimentally investigate the effect of the scan field of view (SFOV) selection and table height settings on the Computed Tomography Dose Index (CTDI) and the implications concerning patient effective and skin dose.MethodsAir-kerma length product (AKLP) measurements were carried out in a helical CT scanner using a pencil type dosimeter positioned in air and inside the holes of a head and a body phantom, using all available SFOV selections and different table height settings. Furthermore, using radiotherapy verification films placed on the CT table surface, the entrance surface air kerma (ESAK) profiles were derived with different SFOV and table height selections, both with and without a phantom on top of the films.ResultsThe AKLP is strongly dependent on the SFOV selection and the table height settings. Different SFOV selections correspond to the selection of different bowtie filters that shape the X-ray beam intensity, resulting in different ESAK values at the isocenter and at the other points within the scanning plane. With the off-center positioning the calculated CTDI values within the center and the periphery of the phantom change also, as a result of the different intensity and width of the X-ray beam to which are exposed to.ConclusionsThe existing protocols for calculating effective dose are limited to only two patient anatomy-SFOV combinations and cannot account for off-center positioning. Therefore, more work will be required to estimate the effective and skin dose for non-standard SFOV-patient anatomy combinations and off-center patient positioning.     

Statistical learning methods as a preprocessing step for survival analysis: evaluation of concept using lung cancer data

Background

Despite its superb lateral resolution, flat-panel-detector (FPD) based tomosynthesis suffers from low contrast and inter-plane artifacts caused by incomplete cancellation of the projection components stemming from outside the focal plane. The incomplete cancellation of the projection components, mostly due to the limited scan angle in the conventional tomosynthesis scan geometry, often makes the image contrast too low to differentiate the malignant tissues from the background tissues with confidence.

Methods

In this paper, we propose a new method to suppress the inter-plane artifacts in FPD-based tomosynthesis. If 3D whole volume CT images are available before the tomosynthesis scan, the CT image data can be incorporated into the tomosynthesis image reconstruction to suppress the inter-plane artifacts, hence, improving the image contrast. In the proposed technique, the projection components stemming from outside the region-of-interest (ROI) are subtracted from the measured tomosynthesis projection data to suppress the inter-plane artifacts. The projection components stemming from outside the ROI are calculated from the 3D whole volume CT images which usually have lower lateral resolution than the tomosynthesis images. The tomosynthesis images are reconstructed from the subtracted projection data which account for the x-ray attenuation through the ROI. After verifying the proposed method by simulation, we have performed both CT scan and tomosynthesis scan on a phantom and a sacrificed rat using a FPD-based micro-CT.

Results

We have measured contrast-to-noise ratio (CNR) from the tomosynthesis images which is an indicator of the residual inter-plane artifacts on the focal-plane image. In both cases of the simulation and experimental imaging studies of the contrast evaluating phantom, CNRs have been significantly improved by the proposed method. In the rat imaging also, we have observed better visual contrast from the tomosynthesis images reconstructed by the proposed method.

Conclusions

The proposed tomosynthesis technique can improve image contrast with aids of 3D whole volume CT images. Even though local tomosynthesis needs extra 3D CT scanning, it may find clinical applications in special situations in which extra 3D CT scan is already available or allowed.     

Automatic exposure control at single- and dual-heartbeat CTCA on a 320-MDCT volume scanner: Effect of heart rate,exposure phase window setting,and reconstruction algorithm

PurposeTo investigate whether electrocardiogram (ECG)-gated single- and dual-heartbeat computed tomography coronary angiography (CTCA) with automatic exposure control (AEC) yields images with uniform image noise at reduced radiation doses.Materials and methodsUsing an anthropomorphic chest CT phantom we performed prospectively ECG-gated single- and dual-heartbeat CTCA on a second-generation 320-multidetector CT volume scanner. The exposure phase window was set at 75%, 70–80%, 40–80%, and 0–100% and the heart rate at 60 or 80 or corr80 bpm; images were reconstructed with filtered back projection (FBP) or iterative reconstruction (IR, adaptive iterative dose reduction 3D). We applied AEC and set the image noise level to 20 or 25 HU. For each technique we determined the image noise and the radiation dose to the phantom center.ResultsWith half-scan reconstruction at 60 bpm, a 70–80% phase window- and a 20-HU standard deviation (SD) setting, the imagenoise level and -variation along the z axis manifested similar curves with FBP and IR. With half-scan reconstruction, the radiation dose to the phantom center with 70–80% phase window was 18.89 and 12.34 mGy for FBP and 4.61 and 3.10 mGy for IR at an SD setting SD of 20 and 25 HU, respectively. At 80 bpm with two-segment reconstruction the dose was approximately twice that of 60 bpm at both SD settings. However, increasing radiation dose at corr80 bpm was suppressed to 1.39 times compared to 60 bpm.ConclusionAEC at ECG-gated single- and dual-heartbeat CTCA controls the image noise at different radiation dose.     

Assessment of the variation in CT scanner performance (image quality and Hounsfield units) with scan parameters,for image optimisation in radiotherapy treatment planning

PurposeTo define a method and investigate how the adjustment of scan parameters affected the image quality and Hounsfield units (HUs) on a CT scanner used for radiotherapy treatment planning. A lack of similar investigations in the literature may be a contributing factor in the apparent reluctance to optimise radiotherapy CT protocols.MethodA Catphan phantom was used to assess how image quality on a Toshiba Aquilion LB scanner changed with scan parameters. Acquisition and reconstruction field-of-view (FOV), collimation, image slice thickness, effective mAs per rotation and reconstruction algorithm were varied. Changes were assessed for HUs of different materials, high contrast spatial resolution (HCSR), contrast-noise ratio (CNR), HU uniformity, scan direction low contrast and CT dose-index.ResultsCNR and HCSR varied most with reconstruction algorithm, reconstruction FOV and effective mAs. Collimation, but not image slice width, had a significant effect on CT dose-index with narrower collimation giving higher doses. Dose increased with effective mAs. Highest HU differences were seen when changing reconstruction algorithm: 56 HU for densities close to water and 117 HU for bone-like materials. Acquisition FOV affected the HUs but reconstruction FOV and effective mAs did not.ConclusionsAll the scan parameters investigated affected the image quality metrics. Reconstruction algorithm, reconstruction FOV, collimation and effective mAs were most important. Reconstruction algorithm and acquisition FOV had significant effect on HU. The methodology is applicable to radiotherapy CT scanners when investigating image quality optimisation, prior to assessing the impact of scan protocol changes on clinical CT images and treatment plans.     

Potential role of nuclear appearance in pathologic recognition and delimitation of sessile serrated polyps of the colon: a karyometric study

OBJECTIVE: To estimate the extent to which surgical pathologists can rely on abnormal nuclear appearance to recognize sessile serrated adenoma (SSP) and to define its extent. STUDY DESIGN: Digitized images of nuclei of superficial crypt cells from SSPs, banal hyperplastic polyps (BHPs), tubular adenomas (TAs) and normal colonic mucosa (N) in surgical pathology specimens were analyzed for size, shape, area, optical density (summed and average) and 22 Markovian texture characteristics. RESULTS: Statistical classification functions correctly distinguished TA from N nuclear image profiles in 93.3% of cases, SSP from N in 70.0% and BHP from N in 74.1%. SSP nuclear feature vectors were less effectively separated from N than BHP and TA on discriminant analysis of the combined data set, and correct classification was achieved in 79.6% of TAs, 53.5% of SSPs and N and 64.2% of BHPs. CONCLUSION: Karyometry distinguished stained nuclei of TAs, but had limited ability to separate nuclei of SSP from those of N and BHP. This suggests that the pathologist attempting to diagnose or delimit the margins of an SSP will find nuclear appearances less helpful than when examining a TA.     

Development of a novel algorithm for metal artifact reduction in digital tomosynthesis using projection-based dual-energy material decomposition for arthroplasty: A phantom study

In this study, a novel dual-energy (DE) material decomposition reconstruction algorithm (DEMDRA) was developed using projection data with the aim of reducing metal artifacts during digital tomosynthesis (DT) for implants. Using the three-material decomposition method and decomposition projection data specific for each material, a novel DEMDRA was implemented to reduce metal artifacts via weighted hybrid reconstructed images [maximum likelihood expectation maximization (MLEM) and shift-and-add (SAA)]. Pulsed X-ray exposures with rapid switching between low and high tube potential kVp were used for DE-DT imaging, and the images were compared using conventional filtered back projection (FBP), MLEM, the simultaneous algebraic reconstruction technique total variation (SART-TV), virtual monochromatic processing, and metal artifact reduction (MAR)-processing algorithms. The reductions in metal artifacts were compared using an artifact index (AI), Gumbel distribution of the largest variations, and the artifact spread functions (ASFs) for prosthesis phantom. The novel DEMDRA yielded an adequately effective overall performance in terms of the AI, and the resulting images yielded good results independently of the type of metal used in the prosthetic phantom, as well as good noise artifact removal, particularly at greater distances from metal objects. Furthermore, the DEMDRA represented the minimum in the model of largest variations. Regarding the ASF analysis, the novel DEMDRA yielded superior metal artifact reduction when compared with conventional reconstruction algorithms with and without MAR processing. Finally, the DEMDRA was particularly useful for reducing high-frequency artifacts.     

生物组织X光照片的层析合成法的研究——Ⅱ.优化层片特性的模拟

在圆层析合成法中,我们提出了两种X光扫描采样的方法.一种是以零阶贝塞尔函数为窗函数,采用单个的采样园锥.另一种是以圆长球函数为窗函数,用Fourier-Bessel级数展开的系数为权,采用多重采样圆锥.在适当选择窗函数的参数下,给出了采样几何学的设计及对不同窗函数的衰减特性进行了模拟.     

Dependency of prescribed CT dose on table height,patient size,and localizer acquisition for one clinical MDCT

PurposeThe purpose of this study was to quantify the effect that table height, patient size, and localizer acquisition order may have on AEC prescribed dose.Method and materialsThree phantoms were used for this study: the Mercury Phantom, acrylic sheets, and an anthropomorphic phantom. A lateral (LAT) and a posterior-anterior (PA) localizer was acquired for each phantom at different table heights on a MDCT scanner (GE Discovery CT750 HD). AEC scan acquisitions were prescribed for each combination of phantom, localizer orientation, and table height ±4 cm with the center position; the displayed CTDI_vol was recorded. Based on the institutional dose monitoring program, the relationship between change in CTDI_vol and change in table height were studied for LAT and AP localizers for clinical exams.ResultsFor all phantom scans based on the PA localizer, the percent change in ranged between −18% and 42% for table heights 4 cm below and above proper centering; while for the LAT localizer, the percent change in CTDI_vol from ideal were no greater than 12% different for ±4 cm differences in table height. Change in CTDI_vol and change in table height displayed a strong linear relationship for AP localizer exams (P = 0.002), and weak correlation for LAT localizer exams (P = 0.12).ConclusionsSince uncertainty in vertical patient positioning is inherently greater than lateral positioning, the LAT localizer should be utilized to precisely and reproducibly deliver the intended amount of radiation prescribed by CT protocols.     

Ultrasound-driven cardiac MRI

PurposeOne of the challenges of cardiac MR imaging is the compensation of respiratory motion, which causes the heart and the surrounding tissues to move. Commonly-used methods to overcome this effect, breath-holding and MR navigation, present shortcomings in terms of available acquisition time or need to periodically interrupt the acquisition, respectively. In this work, an implementation of respiratory motion compensation that obtains information from abdominal ultrasound and continuously adapts the imaged slice position in real time is presented.MethodsA custom workflow was developed, comprising an MR-compatible ultrasound acquisition system, a feature-motion-tracking system with polynomial predictive capability, and a custom MR sequence that continuously adapts the position of the acquired slice according to the tracked position. The system was evaluated on a moving phantom by comparing sharpness and image blurring between static and moving conditions, and in vivo by tracking the motion of the blood vessels of the liver to estimate the cardiac motion. Cine images of the heart were acquired during free breathing.ResultsIn vitro, the predictive motion correction yielded significantly better results than non-predictive or non-corrected acquisitions (p ≪ 0.01). In vivo, the predictive correction resulted in an image quality very similar to the breath-hold acquisition, whereas the uncorrected images show noticeable blurring artifacts.ConclusionIn this work, the possibility of using ultrasound navigation with tracking for the real-time adaptation of MR imaging slices was demonstrated. The implemented technique enabled efficient imaging of the heart with resolutions that would not be feasible in a single breath-hold.     

Simultaneous achievement of accurate CT number and image quality improvement for myocardial perfusion CT at 320-MDCT volume scanning

PurposeTo investigate differences in image-to-image variations between full- and half-scan reconstruction on myocardial CT perfusion (CTP) study.MethodsUsing a cardiac phantom we performed ECG-gated myocardial CTP on a second-generation 320-multidetector CT volume scanner. The heart rate was set at 60 bpm; once per second for a total of 24 s were performed. CT images were acquired at 80- and 120 kVp and subjected to full- and half-scan reconstruction. On images acquired at the same slice level we then measured image-to-image variations, coefficients of variance (CV), and image noise.ResultsThe image-to-image variations with full- and half-scan reconstruction were 1.3 HU vs. 27.2 HU at 80 kVp (p < 0.001) and 0.70 HU vs. 9.3 HU at 120 kVp (p < 0.001) even though the mean HU value was almost the same for both reconstruction methods. The CV of 80- and 120-kVp images of the left ventricular cavity decreased by 0.16% and 0.17%, respectively, with full-scan reconstruction; with half-scan reconstruction it decreased by 3.34% and 2.30%, respectively. Compared with half-scan reconstruction, the image noise was reduced by 27.2% at 80 kVp and by 28.0% at 120 kVp with full-scan reconstruction.ConclusionMyocardial CTP with full-scan reconstruction substantially decreased image-to-image variations and provided accurate CT attenuation.     

Segmentation improvement through denoising of PET images with 3D-context modelling in wavelet domain

Positron emission tomography (PET) images have been incorporated into the radiotherapy process as a powerful tool to assist in the contouring of lesions, leading to the emergence of a broad spectrum of automatic segmentation schemes for PET images (PET-AS). However, not all proposed PET-AS algorithms take into consideration the previous steps of image preparation. PET image noise has been shown to be one of the most relevant affecting factors in segmentation tasks. This study demonstrates a nonlinear filtering method based on spatially adaptive wavelet shrinkage using three-dimensional context modelling that considers the correlation of each voxel with its neighbours. Using this noise reduction method, excellent edge conservation properties are obtained. To evaluate the influence in the segmentation schemes of this filter, it was compared with a set of Gaussian filters (the most conventional) and with two previously optimised edge-preserving filters. Five segmentation schemes were used (most commonly implemented in commercial software): fixed thresholding, adaptive thresholding, watershed, adaptive region growing and affinity propagation clustering. Segmentation results were evaluated using the Dice similarity coefficient and classification error. A simple metric was also included to improve the characterisation of the filters used for induced blurring evaluation, based on the measurement of the average edge width. The proposed noise reduction procedure improves the results of segmentation throughout the performed settings and was shown to be more stable in low-contrast and high-noise conditions. Thus, the capacity of the segmentation method is reinforced by the denoising plan used.     

Characterization of the imaging settings in screening mammography using a tracking and reporting system: A multi-center and multi-vendor analysis

A tracking and reporting system was developed to monitor radiation dose in X-ray breast imaging. We used our tracking system to characterize and compare the mammographic practices of five breast imaging centers located in the United States and Brazil. Clinical data were acquired using eight mammography systems comprising three modalities: computed radiography (CR), full-field digital mammography (FFDM), and digital breast tomosynthesis (DBT). Our database consists of metadata extracted from 334,234 images. We analyzed distributions and correlations of compressed breast thickness (CBT), compression force, target-filter combinations, X-ray tube voltage, and average glandular dose (AGD). AGD reference curves were calculated based on AGD distributions as a function of CBT. These curves represent an AGD reference for a particular population and system. Differences in AGD and imaging settings were attributed to a combination of factors, such as improvements in technology, imaging protocol, and patient demographics. The tracking system allows the comparison of various imaging settings used in screening mammography, as well as the tracking of patient- and population-specific breast data collected from different populations.     

Spectral bandwidth in automated leukocyte classification

Investigators have repeatedly pointed out the importance of spectral information in the automated classification of white blood cells. In general, monochromatic images recorded through two or three color filters are used to extract this information. Although it has generally been thought that the use of narrow band filters provides "cleaner" color information than is obtainable through wide band filters, the choice has not been fully investigated and the question is far from being settled. The use of wide band filters has the clear practical advantage of increased light levels at the detector, resulting in higher signal-to-noise ratio with less demand on light source design. In order to investigate this issue, a series of 681 leukocytes of the most frequently occurring types were digitized by the use of both narrow (10 nm) and wide (90 nm) band filters. Parameters were extracted independently from both sets of images. These parameters were then used to develop a classifier for each set of images. The choice of features and classifier results indicate that there are no major performance differences between the two types of filters.     

A detection method for streak artifacts and radiological noise in a non-uniform region in a CT image

By using the CT images obtained by subtracting two CT images acquired under the same conditions and slice locations, we have devised a method for detecting streak artifacts in non-uniform regions and only radiological noise components in CT images. A chest phantom was scanned using 16- and 64-multidetector row helical CT scanners with various mAs values at 120 kVp. The upper lung slice image was employed as a target image for evaluating the streak artifacts and radiological noise. One hundred parallel line segments with a length of 80 pixels were placed on the subtracted CT image, and the largest CT value in each CT value profile was employed as a feature variable of the streak artifacts; these feature variables were analyzed with the extreme value theory (Gumbel distribution). To detect only the radiological noise, all CT values contained in the 100 line profile were plotted on normal probability paper and the standard deviation was estimated from the inclination of its fitted line for the CT value plots. The two detection methods devised in this study were able to evaluate the streak artifacts and radiological noise in the CT images with high accuracy.     

Location of radiosensitive organs inside pediatric anthropomorphic phantoms: Data required for dosimetry

PurposeThe aim of this study was to determine the location of radiosensitive organs in the interior of four pediatric anthropomorphic phantoms for dosimetric purposes.MethodsFour pediatric anthropomorphic phantoms representing the average individual as newborn, 1-year-old, 5-year-old and 10-year-old child underwent head, thorax and abdomen CT scans. CT and MRI scans of all children aged 0–16 years performed during a 5-year-period in our hospital were reviewed, and 503 were found to be eligible for normal anatomy. Anterior-posterior and lateral dimensions of twelve of the above children closely matched that of the phantoms' head, thoracic and abdominal region in each four phantoms. The mid-sagittal and mid–coronal planes were drawn on selected matching axial images of patients and phantoms. Multiple points outlining large radiosensitive organs in patient images were identified at each slice level and their orthogonal distances from the mid-sagittal and mid–coronal planes were measured. In small organs, the coordinates of organs' centers were similarly determined. The outlines and centers of all radiosensitive organs were reproduced using the coordinates of each organ on corresponding phantoms’ transverse images.ResultsThe locations of the following radiosensitive organs in the interior of the four phantoms was determined: brain, eye lenses, salivary glands, thyroid, lungs, heart, thymus, esophagus, breasts, adrenals, liver, spleen, kidneys, stomach, gallbladder, small bowel, pancreas, colon, ovaries, bladder, prostate, uterus and rectum.ConclusionsThe production of charts of radiosensitive organs inside pediatric anthropomorphic phantoms was feasible and may provide users reliable data for positioning of dosimeters during direct organ dose measurements.     

Novel glucocorticoid effects on acute inflammation in the CNS

The CNS can mount an inflammatory reaction to excitotoxic insults that contributes to the emerging brain damage. Therefore, anti-inflammatory drugs should be beneficial in neurological insults. In contrast, glucocorticoids (GCs), while known for their anti-inflammatory effects, can exacerbate neurotoxicity in the hippocampus after excitotoxic insults. We investigated the effect of GCs on the inflammatory response after a neurological insult. Intact control (INT; intact stress response GC profile), adrenalectomized/GC-supplemented (ADX; low basal GC profile) and GC-treated (COR; chronically high GC profile) rats were injected with kainic acid into the hippocampal CA3 region. Lesion size was determined 8-72 h later. The inflammatory response was characterized using immunohistochemistry, RNAse protection assay and ELISA. The INT and COR rats developed larger CA3 lesions than ADX rats. We found that GCs surprisingly caused an increase in relative numbers of inflammatory cells (granulocytes, monocytes/macrophages and microglia). Additionally, mRNA and protein (IL-1beta and TNF-alpha) levels of the pro-inflammatory cytokines IL-1alpha, IL-1beta and TNF-alpha were elevated in COR rats compared with INT and ADX rats. These data strongly question the traditional view of GCs being uniformly anti-inflammatory and could further explain how GCs worsen the outcome of neurological insults.     

Accuracy of computer-aided geometric 3D reconstruction based on histological serial microgrinding preparation

Purpose: For our research on computer-optimised and automated cochlear implant surgery, we pursue a model-based approach to overcome the limitations of currently available clinical imaging modalities. A serial cross section preparation procedure has been developed and evaluated concerning accuracy to serve for modelling of a digital anatomic atlas to make delicate soft tissue structures available for pre-operative planning.

Methods: A special grinding tool was developed allowing the setting of a specific amount of abrasion as equidistant slice thickness was considered a crucial step. Additionally, each actual abrasion was accurately measured and used during three-dimensional reconstruction of the serial cross-sectional images obtained via digital photo documentation after each microgrinding step. A well-known reference object was prepared using this procedure and evaluated in terms of accuracy.

Results: Reconstruction of the whole sample was achieved with an error less than 0.4%, and the edge lengths in the direction of abrasion could be reconstructed with an average error of 0.6 ± 0.3 mm; both prove the realisation of equidistant abrasion. Using artificial registration fiducials and a custom-made algorithm for image alignment, parallelism and rectangularity could be preserved with average errors less than 0.4° ± 0.3°.

Conclusion: We present a systematic, practicable and reliable method for the geometrically accurate reconstruction of anatomical structures, which is especially suitable for the middle and inner ear anatomy including soft tissue structures. For the first time, the quality of such a reconstruction process has been quantified and successfully proven for its usability.