Automatic prostate segmentation in cone-beam computed tomography images using rigid registration.

We propose to evaluate automatic three-dimensional gray-value rigid registration (RR) methods for prostate localization on cone-beam computed tomography (CBCT) scans. In total, 103 CBCT scans of 9 prostate patients have been analyzed. Each one was registered to the planning CT scan using different methods: (a) global RR, (b) pelvis bone structure RR, (c) bone RR refined by local soft-tissue RR using the CT clinical target volume (CTV) expanded with a 1, 3, 5, 8, 10, 12, 15 or 20-mm margin. To evaluate results, a radiation oncologist was asked to manually delineate the CTV on the CBCT scans. The Dice coefficients between each automatic CBCT segmentation - derived from the transformation of the manual CT segmentation - and the manual CBCT segmentation were calculated. Global or bone CT/CBCT RR has been shown to yield insufficient results in average. Local RR with an 8-mm margin around the CTV after bone RR was found to be the best candidate for systematically significantly improving prostate localization.

Constraining 3-D electrical resistance tomography with GPR reflection data forimproved aquifer characterization

Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240×40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with then underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (N300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.

Rib fractures at postmortem computed tomography (PMCT) validated against the autopsy

To evaluate the sensitivity of postmortem computed tomography (PMCT) in rib fracture detection validated against autopsy. Fifty-one forensic cases underwent a postmortem CT prior to forensic autopsy. Two image readers (radiologist and forensic pathologist) assessed high resolution CT data sets for rib fractures. Correct recognition rates (CRR), sensitivity and specificity values were calculated over all observations as well as individually for every rib and region. Additionally, for partial rib fractures the sensitivity of autopsy was calculated vice versa. 3876 entries in each study protocol (autopsy, PMCT radiologist and PMCT forensic pathologist) were investigated. A total of 690 fractures (autopsy), 491 (PMCT and radiologist) and 559 (PMCT and forensic pathologist) were detected. The CRR was 0.85. Sensitivity and specificity of PMCT for rib fracture detection were 0.63 (0.58 radiologist, 0.68 forensic pathologist) and 0.97 (both readers 0.97), respectively. Low CRR and sensitivity values were obtained for antero-lateral fractures. Partial rib fractures were better detected by PMCT. PMCT has a rather low sensitivity for rib fracture detection when validated against autopsy and indicates that clinical CT may also demonstrate a reasonable number of false negatives. Partial rib fractures often remain undetected at autopsy.

Immunoscintigraphie avec tomographie d'émission dans les cancers de la thyroïde

Immunoscintigraphy (IS) consists of in vivo body structure imaging using a specific labelled antibody to an antigen concentrated in the structure under study. Technically, the image contrast is better when IS is performed with a computerized emission tomography system. High concentrations of carcinoembryonic antigen (CEA) have been reported in medullary thyroid carcinoma and thyroglobulin (Tg) is a marker for differentiated thyroid carcinoma. We used injections of 131-I labelled monoclonal antibodies to CEA and Tg to detect thyroid tumours. A feasibility trial using anti-CEA antibodies gave very encouraging results. However, only tumours larger than 10 cm3 could be detected. Contradictory results were obtained using anti-Tg antibodies but this data must be considered as preliminary. Various means of improving the method and the concept of accessibility of the antigen to the antibody in vivo are discussed. This study shows IS to be a promising experimental technique. Further studies are required to define its clinical indications before it can be advocated for routine use.

A realistic Computed Tomography simulator for small motion analysis of cerebral aneurysms.

This paper describes a realistic simulator for the Computed Tomography (CT) scan process for motion analysis. In fact, we are currently developing a new framework to find small motion from the CT scan. In order to prove the fidelity of this framework, or potentially any other algorithm, we present in this paper a simulator to simulate the whole CT acquisition process with a priori known parameters. In other words, it is a digital phantom for the motion analysis that can be used to compare the results of any related algorithm with the ground-truth realistic analytical model. Such a simulator can be used by the community to test different algorithms in the biomedical imaging domain. The most important features of this simulator are its different considerations to simulate the best the real acquisition process and its generality.

Validity of the acoustic approximation in full-waveform seismic crosshole tomography

Acoustic waveform inversions are an increasingly popular tool for extracting subsurface information from seismic data. They are computationally much more efficient than elastic inversions. Naturally, an inherent disadvantage is that any elastic effects present in the recorded data are ignored in acoustic inversions. We investigate the extent to which elastic effects influence seismic crosshole data. Our numerical modeling studies reveal that in the presence of high contrast interfaces, at which P-to-S conversions occur, elastic effects can dominate the seismic sections, even for experiments involving pressure sources and pressure receivers. Comparisons of waveform inversion results using a purely acoustic algorithm on synthetic data that is either acoustic or elastic, show that subsurface models comprising small low-to-medium contrast (?30%) structures can be successfully resolved in the acoustic approximation. However, in the presence of extended high-contrast anomalous bodies, P-to-S-conversions may substantially degrade the quality of the tomographic images. In particular, extended low-velocity zones are difficult to image. Likewise, relatively small low-velocity features are unresolved, even when advanced a priori information is included. One option for mitigating elastic effects is data windowing, which suppresses later arriving seismic arrivals, such as shear waves. Our tests of this approach found it to be inappropriate because elastic effects are also included in earlier arriving wavetrains. Furthermore, data windowing removes later arriving P-wave phases that may provide critical constraints on the tomograms. Finally, we investigated the extent to which acoustic inversions of elastic data are useful for time-lapse analyses of high contrast engineered structures, for which accurate reconstruction of the subsurface structure is not as critical as imaging differential changes between sequential experiments. Based on a realistic scenario for monitoring a radioactive waste repository, we demonstrated that acoustic inversions of elastic data yield substantial distortions of the tomograms and also unreliable information on trends in the velocity changes.

Randomized trial of a computerized coronary heart disease risk assessment tool in HIV-infected patients receiving combination antiretroviral therapy.

BACKGROUND: Exposure to combination antiretroviral therapy (cART) can lead to important metabolic changes and increased risk of coronary heart disease (CHD). Computerized clinical decision support systems have been advocated to improve the management of patients at risk for CHD but it is unclear whether such systems reduce patients' risk for CHD. METHODS: We conducted a cluster trial within the Swiss HIV Cohort Study (SHCS) of HIV-infected patients, aged 18 years or older, not pregnant and receiving cART for >3 months. We randomized 165 physicians to either guidelines for CHD risk factor management alone or guidelines plus CHD risk profiles. Risk profiles included the Framingham risk score, CHD drug prescriptions and CHD events based on biannual assessments, and were continuously updated by the SHCS data centre and integrated into patient charts by study nurses. Outcome measures were total cholesterol, systolic and diastolic blood pressure and Framingham risk score. RESULTS: A total of 3,266 patients (80% of those eligible) had a final assessment of the primary outcome at least 12 months after the start of the trial. Mean (95% confidence interval) patient differences where physicians received CHD risk profiles and guidelines, rather than guidelines alone, were total cholesterol -0.02 mmol/l (-0.09-0.06), systolic blood pressure -0.4 mmHg (-1.6-0.8), diastolic blood pressure -0.4 mmHg (-1.5-0.7) and Framingham 10-year risk score -0.2% (-0.5-0.1). CONCLUSIONS: Systemic computerized routine provision of CHD risk profiles in addition to guidelines does not significantly improve risk factors for CHD in patients on cART.

[(18)F]Fluoroethyltyrosine- positron emission tomography-guided radiotherapy for high-grade glioma.

BACKGROUND: To compare morphological gross tumor volumes (GTVs), defined as pre- and postoperative gadolinium enhancement on T1-weighted magnetic resonance imaging to biological tumor volumes (BTVs), defined by the uptake of (18)F fluoroethyltyrosine (FET) for the radiotherapy planning of high-grade glioma, using a dedicated positron emission tomography (PET)-CT scanner equipped with three triangulation lasers for patient positioning. METHODS: Nineteen patients with malignant glioma were included into a prospective protocol using FET PET-CT for radiotherapy planning. To be eligible, patients had to present with residual disease after surgery. Planning was performed using the clinical target volume (CTV = GTV union or logical sum BTV) and planning target volume (PTV = CTV + 20 mm). First, the interrater reliability for BTV delineation was assessed among three observers. Second, the BTV and GTV were quantified and compared. Finally, the geometrical relationships between GTV and BTV were assessed. RESULTS: Interrater agreement for BTV delineation was excellent (intraclass correlation coefficient 0.9). Although, BTVs and GTVs were not significantly different (p = 0.9), CTVs (mean 57.8 +/- 30.4 cm(3)) were significantly larger than BTVs (mean 42.1 +/- 24.4 cm(3); p < 0.01) or GTVs (mean 38.7 +/- 25.7 cm(3); p < 0.01). In 13 (68%) and 6 (32%) of 19 patients, FET uptake extended >or= 10 and 20 mm from the margin of the gadolinium enhancement. CONCLUSION: Using FET, the interrater reliability had excellent agreement for BTV delineation. With FET PET-CT planning, the size and geometrical location of GTVs and BTVs differed in a majority of patients.

Stiffness tomography by atomic force microscopy.

The atomic force microscope is a convenient tool to probe living samples at the nanometric scale. Among its numerous capabilities, the instrument can be operated as a nano-indenter to gather information about the mechanical properties of the sample. In this operating mode, the deformation of the cantilever is displayed as a function of the indentation depth of the tip into the sample. Fitting this curve with different theoretical models permits us to estimate the Young's modulus of the sample at the indentation spot. We describe what to our knowledge is a new technique to process these curves to distinguish structures of different stiffness buried into the bulk of the sample. The working principle of this new imaging technique has been verified by finite element models and successfully applied to living cells.

Update on the non-prewhitening model observer in computed tomography for the assessment of the adaptive statistical and model-based iterative reconstruction algorithms.

The state of the art to describe image quality in medical imaging is to assess the performance of an observer conducting a task of clinical interest. This can be done by using a model observer leading to a figure of merit such as the signal-to-noise ratio (SNR). Using the non-prewhitening (NPW) model observer, we objectively characterised the evolution of its figure of merit in various acquisition conditions. The NPW model observer usually requires the use of the modulation transfer function (MTF) as well as noise power spectra. However, although the computation of the MTF poses no problem when dealing with the traditional filtered back-projection (FBP) algorithm, this is not the case when using iterative reconstruction (IR) algorithms, such as adaptive statistical iterative reconstruction (ASIR) or model-based iterative reconstruction (MBIR). Given that the target transfer function (TTF) had already shown it could accurately express the system resolution even with non-linear algorithms, we decided to tune the NPW model observer, replacing the standard MTF by the TTF. It was estimated using a custom-made phantom containing cylindrical inserts surrounded by water. The contrast differences between the inserts and water were plotted for each acquisition condition. Then, mathematical transformations were performed leading to the TTF. As expected, the first results showed a dependency of the image contrast and noise levels on the TTF for both ASIR and MBIR. Moreover, FBP also proved to be dependent of the contrast and noise when using the lung kernel. Those results were then introduced in the NPW model observer. We observed an enhancement of SNR every time we switched from FBP to ASIR to MBIR. IR algorithms greatly improve image quality, especially in low-dose conditions. Based on our results, the use of MBIR could lead to further dose reduction in several clinical applications.

Computed tomography commissioning programmes: how to obtain a reliable MTF with an automatic approach?

The purpose of this study was to assess the spatial resolution of a computed tomography (CT) scanner with an automatic approach developed for routine quality controls when varying CT parameters. The methods available to assess the modulation transfer functions (MTF) with the automatic approach were Droege's and the bead point source (BPS) methods. These MTFs were compared with presampled ones obtained using Boone's method. The results show that Droege's method is not accurate in the low-frequency range, whereas the BPS method is highly sensitive to image noise. While both methods are well adapted to routine stability controls, it was shown that they are not able to provide absolute measurements. On the other hand, Boone's method, which is robust with respect to aliasing, more resilient to noise and provides absolute measurements, satisfies the commissioning requirements perfectly. Thus, Boone's method combined with a modified Catphan 600 phantom could be a good solution to assess CT spatial resolution in the different CT planes.