Scientific assessment of the quality of OSIRIS images

Context. OSIRIS, the scientific imaging system onboard the ESA Rosetta spacecraft, has been imaging the nucleus of comet 67P/Churyumov-Gerasimenko and its dust and gas environment since March 2014. The images serve different scientific goals, from morphology and composition studies of the nucleus surface, to the motion and trajectories of dust grains, the general structure of the dust coma, the morphology and intensity of jets, gas distribution, mass loss, and dust and gas production rates. Aims. We present the calibration of the raw images taken by OSIRIS and address the accuracy that we can expect in our scientific results based on the accuracy of the calibration steps that we have performed. Methods. We describe the pipeline that has been developed to automatically calibrate the OSIRIS images. Through a series of steps, radiometrically calibrated and distortion corrected images are produced and can be used for scientific studies. Calibration campaigns were run on the ground before launch and throughout the years in flight to determine the parameters that are used to calibrate the images and to verify their evolution with time. We describe how these parameters were determined and we address their accuracy. Results. We provide a guideline to the level of trust that can be put into the various studies performed with OSIRIS images, based on the accuracy of the image calibration.

Personalized X-Ray Reconstruction of the Proximal Femur via Intensity-Based Non-rigid 2D-3D Registration

This paper presents a new approach for reconstructing a patient-specific shape model and internal relative intensity distribution of the proximal femur from a limited number (e.g., 2) of calibrated C-arm images or X-ray radiographs. Our approach uses independent shape and appearance models that are learned from a set of training data to encode the a priori information about the proximal femur. An intensity-based non-rigid 2D-3D registration algorithm is then proposed to deformably fit the learned models to the input images. The fitting is conducted iteratively by minimizing the dissimilarity between the input images and the associated digitally reconstructed radiographs of the learned models together with regularization terms encoding the strain energy of the forward deformation and the smoothness of the inverse deformation. Comprehensive experiments conducted on images of cadaveric femurs and on clinical datasets demonstrate the efficacy of the present approach.

Automatic Scan Planning for Magnetic Resonance Imaging of the Knee Joint

Automatic scan planning for magnetic resonance imaging of the knee aims at defining an oriented bounding box around the knee joint from sparse scout images in order to choose the optimal field of view for the diagnostic images and limit acquisition time. We propose a fast and fully automatic method to perform this task based on the standard clinical scout imaging protocol. The method is based on sequential Chamfer matching of 2D scout feature images with a three-dimensional mean model of femur and tibia. Subsequently, the joint plane separating femur and tibia, which contains both menisci, can be automatically detected using an information-augmented active shape model on the diagnostic images. This can assist the clinicians in quickly defining slices with standardized and reproducible orientation, thus increasing diagnostic accuracy and also comparability of serial examinations. The method has been evaluated on 42 knee MR images. It has the potential to be incorporated into existing systems because it does not change the current acquisition protocol.

Comparison of contrast-enhanced multi-station MR angiography and digital subtraction angiography of the lower extremity arterial disease

PURPOSE: To compare diagnostic accuracy of multi-station, high-spatial resolution contrast-enhanced MR angiography (CE-MRA) of the lower extremities with digital subtraction angiography (DSA) as the reference standard in patients with symptomatic peripheral arterial occlusive disease. MATERIALS AND METHODS: Of 485 consecutive patients undergoing a run-off CE-MRA, 152 patients (86 male, 66 female; mean age, 71.6 years) with suspected peripheral arterial occlusive disease were included into our Institutional Review Board approved study. All patients underwent MRA and DSA of the lower extremities within 30 days. MRA was performed at 1.5 Tesla with a single bolus of 0.1 mmol/kg body weight of gadobutrol administered at a rate of 2.0 mL/s at three stations. Two readers evaluated the MRA images independently for stenosis grade and image quality. Sensitivity and specificity were derived. RESULTS: Sensitivity and specificity ranged from 73% to 93% and 64% to 89% and were highest in the thigh area. Both readers showed comparable results. Evaluation of good and better quality MRAs resulted in a considerable improvement in diagnostic accuracy. CONCLUSION: Contrast-enhanced MRA demonstrates good sensitivity and specificity in the investigation of the vasculature of the lower extremities. While a minor investigator experience dependence remains, it is standardizable and shows good inter-observer agreement. Our results confirm that the administration of Gadobutrol at a standard dose of 0.1 mmol/kg for contrast-enhanced runoff MRA is able to detect hemodynamically relevant stenoses. Use of contrast-enhanced MRA as an alternative to intra-arterial DSA in the evaluation and therapeutic planning of patients with suspected peripheral arterial occlusive disease is well justified. J. Magn. Reson. Imaging 2013;37:1427-1435. © 2012 Wiley Periodicals, Inc.

In Vitro and In Vivo Imaging Characteristics Assessment of Polymeric Coils Compared with Standard Platinum Coils for the Treatment of Intracranial Aneurysms

BACKGROUND AND PURPOSE:Conventional platinum coils cause imaging artifacts that reduce imaging quality and therefore impair imaging interpretation on intraprocedural or noninvasive follow-up imaging. The purpose of this study was to evaluate imaging characteristics and artifact production of polymeric coils compared with standard platinum coils in vitro and in vivo.MATERIALS AND METHODS:Polymeric coils and standard platinum coils were evaluated in vitro with the use of 2 identical silicon aneurysm models coiled with a packing attenuation of 20% each. DSA, flat panel CT, CT, and MR imaging were performed. In vivo evaluation of imaging characteristics of polymeric coils was performed in experimentally created rabbit carotid bifurcation aneurysms. DSA, CT/CTA, and MR imaging were performed after endovascular treatment of the aneurysms. Images were evaluated regarding visibility of individual coils, coil mass, artifact production, and visibility of residual flow within the aneurysm.RESULTS:Overall, in vitro and in vivo imaging showed relevantly reduced artifact production of polymeric coils in all imaging modalities compared with standard platinum coils. Image quality of CT and MR imaging was improved with the use of polymeric coils, which permitted enhanced depiction of individual coil loops and residual aneurysm lumen as well as the peri-aneurysmal area. Remarkably, CT images demonstrated considerably improved image quality with only minor artifacts compared with standard coils. On DSA, polymeric coils showed transparency and allowed visualization of superimposed vessel structures.CONCLUSIONS:This initial experimental study showed improved imaging quality with the use of polymeric coils compared with standard platinum coils in all imaging modalities. This might be advantageous for improved intraprocedural imaging for the detection of complications and posttreatment noninvasive follow-up imaging.

Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and point distribution model

Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and a point distribution model is discussed. We present a 2D/3D reconstruction scheme combining statistical extrapolation and regularized shape deformation with an iterative image-to-model correspondence establishing algorithm, and show its application to reconstruct the surface of proximal femur. The image-to-model correspondence is established using a non-rigid 2D point matching process, which iteratively uses a symmetric injective nearest-neighbor mapping operator and 2D thin-plate splines based deformation to find a fraction of best matched 2D point pairs between features detected from the fluoroscopic images and those extracted from the 3D model. The obtained 2D point pairs are then used to set up a set of 3D point pairs such that we turn a 2D/3D reconstruction problem to a 3D/3D one. We designed and conducted experiments on 11 cadaveric femurs to validate the present reconstruction scheme. An average mean reconstruction error of 1.2 mm was found when two fluoroscopic images were used for each bone. It decreased to 1.0 mm when three fluoroscopic images were used.

A rare anomaly of the anterior communicating artery complex hidden by a large broad-neck aneurysm and disclosed by three-dimensional rotational angiography

Double fenestration of the anterior communicating artery (ACoA) complex associated with an aneurysm is a very rare finding and is usually caused by ACoA duplication and the presence of a median artery of the corpus callosum (MACC). We present a patient in whom double fenestration was not associated with ACoA duplication or even with MACC, representing therefore, a previously unreported anatomic variation. A 43 year old woman experienced sudden headache and the CT scans showed subarachnoid haemorrhage (SAH). On admission, her clinical condition was consistent with Hunt and Hess grade II. Conventional digital subtraction angiography (DSA) was performed and revealed multiple intracranial aneurysms arising from both middle cerebral arteries (MCA) and from the ACoA. Three-dimensional rotational angiography (3D-RA) disclosed a double fenestration of the ACoA complex which was missed by DSA. The patient underwent a classic pterional approach in order to achieve occlusion of both left MCA and ACoA aneurysms by surgical clipping. The post-operative period was uneventful. A rare anatomical variation characterised by a double fenestration not associated with ACoA duplication or MACC is described. The DSA images missed the double fenestration which was disclosed by 3D-RA, indicating the importance of 3D-RA in the diagnosis and surgical planning of intracranial aneurysms.

Aperistaltic effect of hyoscine N-butylbromide versus glucagon on the small bowel assessed by magnetic resonance imaging

The aim of this prospective study was to compare the intraindividual aperistaltic effect of 40 mg hyoscine N-butylbromide (HBB/Buscopan) with that of 1 mg glucagon on small bowel motility by using magnetic resonance imaging (MRI). Ten healthy volunteers underwent two separate 1.5-T MRI studies (HBB/glucagon) after a standardized oral preparation with an aqueous solution of Gd-DOTA and ispaghula (Metamucil). A 2D T1-w GRE sequence was acquired (TR 2.7 ms/TE 1.3 ms, temporal resolution 0.25 s) before and after intravenous (i.v.) drug administration and motility was followed over 1 h. On the resulting images the cross-sectional luminal diameters were assessed and plotted over time. Baseline motility frequency, onset of aperistalsis, duration of arrest, reappearance of motility and return to normal motility were analysed. Significant differences regarding reliability and duration of aperistalsis were observed. In the HBB group aperistalsis lasted a mean of 6.8 +/- 5.3 min compared with 18.3 +/- 7 min after glucagon (p < 0.0001). In 50% of cases HBB did not accomplish aperistalsis, whereas glucagon always succeeded (p = 0.05). There were no significant differences in terms of baseline and end frequencies for the onset of aperistalsis (22.2 +/- 37.5 s HBB/13.4 +/- 9.2 s glucagon, p = 0.1), nor for the return to normal motility. Arrest of small bowel motion is achieved more reliably and lasts significantly longer after i.v. administration of 1 mg glucagon compared with 40 mg HBB.

Scintigraphic and radiographic appearance of the sacroiliac region in horses with gait abnormalities or poor performance

Nuclear bone scintigraphy is commonly used in the diagnosis of sacroiliac disease in horses. The aim of the present retrospective study was to determine if there was an association between radiopharmaceutical uptake pattern and radiographic appearance of the sacroiliac region in horses. Seventy-nine horses undergoing bone scintigraphy with Tc-99 m-HDP and radiography of the pelvis because of lameness or poor performance were studied. Subjective and semiquantitative methods were used to characterize and compare radiopharmaceutical uptake between horses. Ventrodorsal radiographs of the region were obtained and were evaluated. Subjectively, 70 horses (88.6%) had an abnormal uptake pattern. In nine horses, the sacroiliac region was normal (11.4%). There was no association between subjective evaluation of the scintigraphic images and semiquantitative methods. There was a significant association between radiopharmaceutical uptake and conformation (T- or Y-like form) and shape (butterfly-, wing-, leaf-, or horn-like) of the sacrum. The radiopharmaceutical uptake of the tubera sacralia was significantly higher in males (left side P = 0.002, right side P = 0.003). In conclusion, the conformation of the sacrum may play an important role in the scintigraphic appearance and may be the cause of increased radiopharmaceutical uptake.

Magnetic resonance imaging anatomy of the normal equine larynx and pharynx

The purpose of the present study was to describe normal magnetic resonance (MR) imaging anatomy of the equine larynx and pharynx and to present the optimal protocol, sequences, and possible limitations of this examination technique. Using a 0.3 T unit, the laryngeal and pharyngeal regions was imaged in two horses. The protocol consisted of sagittal and transverse T2-weighted (T2w) fast spin echo, transverse T1-weighted (T1w) spin echo, and dorsal high-resolution T1w gradient echo (both pre- and postcontrast enhancement) sequences. Euthanasia was performed at the end of the imaging procedure. Macroscopic anatomy of the cadaver sections were compared with the MR images in transverse, midsagittal, and parasagittal planes. There was good differentiation of anatomic structures, including soft tissues. The laryngeal cartilages, hyoid apparatus, and upper airway muscle groups with their attachments could be clearly identified. However, it was not always possible to delineate individual muscles in each plane. Most useful were both T2w and T1w transverse sequences. Intravenous application of contrast medium was helpful to identify blood vessels. The MR images corresponded with the macroscopic anatomy of cadaver sections.

A complementary radiographic projection of the equine temporomandibular joint

The complexity of the equine skull makes the temporomandibular joint a difficult area to evaluate radiographically. The goal of this study was to determine the optimal angle for a complementary radiographic projection of the equine temporomandibular joint based on a computed tomography (CT) cadaver study. CT was performed on six equine cadaver heads of horses that were euthanized for other reasons than temporomandibular joint disease. After the CT examination, 3D reconstruction of the equine skull was performed to subjectively determine the angle for a complementary radiographic projection of the temporomandibular joint. The angle was measured on the left and right temporomandibular joint of each head. Based on the measurements obtained from the CT images, a radiographic projection of the temporomandibular joint in a rostra-145 degrees ventral-caudodorsal oblique (R45 degrees V-CdDO) direction was developed by placing the X-ray unit 30 degrees laterally, maintaining at the same time the R45 degrees V-CdDO angle (R45 degrees V30 degrees L-CdDLO). This radiographic projection was applied to all cadaver heads and on six live horses. In three of the live horses abnormal findings associated with the temporomandibular joint were detected. We conclude that this new radiographic projection of the temporomandibular joint provides superior visualization of the temporomandibular joint space and the articular surface of the mandibular condyle.

Modeling the heterogeneous ice and gas coma of Comet 103P/Hartley 2

The spectacular images of Comet 103P/Hartley 2 recorded by the Medium Resolution Instrument (MRI) and High Resolution Instrument (HRI) on board of the Extrasolar Planet Observation and Deep Impact Extended Investigation (EPOXI) spacecraft, as the Deep Impact extended mission, revealed that its bi-lobed very active nucleus outgasses volatiles heterogeneously. Indeed, CO2 is the primary driver of activity by dragging out chunks of pure ice out of the nucleus from the sub-solar lobe that appear to be the main source of water in Hartley 2's coma by sublimating slowly as they go away from the nucleus. However, water vapor is released by direct sublimation of the nucleus at the waist without any significant amount of either CO2 or icy grains. The coma structure for a comet with such areas of diverse chemistry differs from the usual models where gases are produced in a homogeneous way from the surface. We use the fully kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J. 685, 659-677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J. 732, 104-120) applied to Comet 103P/Hartley 2 including sublimating icy grains to reproduce the observations made by EPOXI and ground-based measurements. A realistic bi-lobed nucleus with a succession of active areas with different chemistry was included in the model enabling us to study in details the coma of Hartley 2. The different gas production rates from each area were found by fitting the spectra computed using a line-by-line non-LTE radiative transfer model to the HRI observations. The presence of icy grains with long lifetimes, which are pushed anti-sunward by radiation pressure, explains the observed OH asymmetry with enhancement on the night side of the coma.

Spectrophotometric properties of the nucleus of comet 67P/Churyumov-Gerasimenko from the OSIRIS instrument onboard the ROSETTA spacecraft

Context. The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. Aims. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The observations presented here were performed during July and the beginning of August 2014, during the approach phase, when OSIRIS was mapping the surface of the comet with several filters at different phase angles (1.3 degrees-54 degrees). The resolution reached up to 2.1 m/px. Methods. The OSIRIS images were processed with the OSIRIS standard pipeline, then converted into I/F. radiance factors and corrected for the illumination conditions at each pixel using the Lommel-Seeliger disk law. Color cubes of the surface were produced by stacking registered and illumination-corrected images. Furthermore, photometric analysis was performed both on disk-averaged photometry in several filters and on disk-resolved images acquired with the NAC orange filter, centered at 649 ran, using Hapke modeling. Results. The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13 +/- 0.01 in the HG system formalism and an absolute magnitude H-v(1, 1, 0) = 15.74 +/- 0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at similar to 290 rim that is possibly due to SO2 ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11%/( 100 nm) to 16%/(100 nm) in the 1.3 degrees-54 degrees phase angle range. The geometric albedo of the comet is 6.5 +/- 0.2% at 649 nm, with local variations of up to similar to 16% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

Building the View Graph of a Category by Exploiting Image Realism

We propose a weakly supervised method to arrange images of a given category based on the relative pose between the camera and the object in the scene. Relative poses are points on a sphere centered at the object in a given canonical pose, which we call object viewpoints. Our method builds a graph on this sphere by assigning images with similar viewpoint to the same node and by connecting nodes if they are related by a small rotation. The key idea is to exploit a large unlabeled dataset to validate the likelihood of dominant 3D planes of the object geometry. A number of 3D plane hypotheses are evaluated by applying small 3D rotations to each hypothesis and by measuring how well the deformed images match other images in the dataset. Correct hypotheses will result in deformed images that correspond to plausible views of the object, and thus will likely match well other images in the same category. The identified 3D planes are then used to compute affinities between images related by a change of viewpoint. We then use the affinities to build a view graph via a greedy method and the maximum spanning tree.

A comparison between 2D and 3D cephalometry on CBCT scans of human skulls

The purpose of this study was to evaluate whether measurements on conventional cephalometric radiographs are comparable with 3D measurements on 3D models of human skulls, derived from cone beam CT (CBCT) data. A CBCT scan and a conventional cephalometric radiograph were made of 40 dry skulls. Standard cephalometric software was used to identify landmarks on both the 2D images and the 3D models. The same operator identified 17 landmarks on the cephalometric radiographs and on the 3D models. All images and 3D models were traced five times with a time-interval of 1 week and the mean value of repeated measurements was used for further statistical analysis. Distances and angles were calculated. Intra-observer reliability was good for all measurements. The reproducibility of the measurements on the conventional cephalometric radiographs was higher compared with the reproducibility of measurements on the 3D models. For a few measurements a clinically relevant difference between measurements on conventional cephalometric radiographs and 3D models was found. Measurements on conventional cephalometric radiographs can differ significantly from measurements on 3D models of the same skull. The authors recommend that 3D tracings for longitudinal research are not used in cases were there are only 2D records from the past.