Thoracic radiotherapy treatment planning with cine PET/CT

Purpose: Respiratory motion causes substantial uncertainty in radiotherapy treatment planning. Four-dimensional computed tomography (4D-CT) is a useful tool to image tumor motion during normal respiration. Treatment margins can be reduced by targeting the motion path of the tumor. The expense and complexity of 4D-CT, however, may be cost-prohibitive at some facilities. We developed an image processing technique to produce images from cine CT that contain significant motion information without 4D-CT. The purpose of this work was to compare cine CT and 4D-CT for the purposes of target delineation and dose calculation, and to explore the role of PET in target delineation of lung cancer. Methods: To determine whether cine CT could substitute 4D-CT for small mobile lung tumors, we compared target volumes delineated by a physician on cine CT and 4D-CT for 27 tumors with intrafractional motion greater than 1 cm. We assessed dose calculation by comparing dose distributions calculated on respiratory-averaged cine CT and respiratory-averaged 4D-CT using the gamma index. A threshold-based PET segmentation model of size, motion, and source-to-background was developed from phantom scans and validated with 24 lung tumors. Finally, feasibility of integrating cine CT and PET for contouring was assessed on a small group of larger tumors. Results: Cine CT to 4D-CT target volume ratios were (1.05��0.14) and (0.97��0.13) for high-contrast and low-contrast tumors respectively which was within intraobserver variation. Dose distributions on cine CT produced good agreement (< 2%/1 mm) with 4D-CT for 71 of 73 patients. The segmentation model fit the phantom data with R2 = 0.96 and produced PET target volumes that matched CT better than 6 published methods (-5.15%). Application of the model to more complex tumors produced mixed results and further research is necessary to adequately integrate PET and cine CT for delineation. Conclusions: Cine CT can be used for target delineation of small mobile lesions with minimal differences to 4D-CT. PET, utilizing the segmentation model, can provide additional contrast. Additional research is required to assess the efficacy of complex tumor delineation with cine CT and PET. Respiratory-averaged cine CT can substitute respiratory-averaged 4D-CT for dose calculation with negligible differences.

Differential uptake of (68)Ga-DOTATOC and (68)Ga-DOTATATE in PET/CT of gastroenteropancreatic neuroendocrine tumors

PURPOSE Abundant expression of somatostatin receptors (sst) is a characteristic of neuroendocrine tumors (NET). Thus, radiolabeled somatostatin analogs have emerged as important tools for both in vivo diagnosis and therapy of NET. The two compounds most often used in functional imaging with positron emission tomography (PET) are (68)Ga-DOTATATE and (68)Ga-DOTATOC. Both analogs share a quite similar sst binding profile. However, the in vitro affinity of (68)Ga-DOTATATE in binding the sst subtype 2 (sst2) is approximately tenfold higher than that of (68)Ga-DOTATOC. This difference may affect their efficiency in detection of NET lesions, as sst2 is the predominant receptor subtype on gastroenteropancreatic NET. We thus compared the diagnostic value of PET/CT with both radiolabeled somatostatin analogs ((68)Ga-DOTATATE and (68)Ga-DOTATOC) in the same patients with gastroenteropancreatic NET. PATIENTS AND METHODS Twenty-seven patients with metastatic gastroenteropancreatic NET underwent (68)Ga-DOTATOC and (68)Ga-DOTATATE PET/CT as part of the workup before prospective peptide receptor radionuclide therapy (PRRT). The performance of both imaging methods was analyzed and compared for detection of individual lesions per patient and for eight defined body regions. A region was regarded as positive if at least one lesion was detected in that region. In addition, radiopeptide uptake in terms of the maximal standardized uptake value (SUV(max)) was compared for concordant lesions and renal parenchyma. RESULTS Fifty-one regions were found positive with both (68)Ga-DOTATATE and (68)Ga-DOTATOC. Overall, however, significantly fewer lesions were detected with (68)Ga-DOTATATE in comparison with (68)Ga-DOTATOC (174 versus 179, p < 0.05). Mean (68)Ga-DOTATATE SUV(max) across all lesions was significantly lower compared with (68)Ga-DOTATOC (16.9 �� 6.8 versus 22.1 �� 12.0, p < 0.01). Mean SUV(max) for renal parenchyma was not significantly different between (68)Ga-DOTATATE and (68)Ga-DOTATOC (12.6 �� 2.6 versus 12.6 �� 2.7). CONCLUSIONS (68)Ga-DOTATOC and (68)Ga-DOTATATE possess similar diagnostic accuracy for detection of gastroenteropancreatic NET lesions (with a potential advantage of (68)Ga-DOTATOC) despite their evident difference in affinity for sst2. Quite unexpectedly, maximal uptake of (68)Ga-DOTATOC tended to be higher than its (68)Ga-DOTATATE counterpart. However, tumor uptake shows high inter- and intraindividual variance with unpredictable preference of one radiopeptide. Thus, our data encourage the application of different sst ligands to enable personalized imaging and therapy of gastroenteropancreatic NET with optimal targeting of tumor receptors.

Transconvolution and the virtual positron emission tomograph--a new method for cross calibration in quantitative PET���CT imaging

PURPOSE Positron emission tomography (PET)���computed tomography (CT) measurements on small lesions are impaired by the partial volume effect, which is intrinsically tied to the point spread function of the actual imaging system, including the reconstruction algorithms. The variability resulting from different point spread functions hinders the assessment of quantitative measurements in clinical routine and especially degrades comparability within multicenter trials. To improve quantitative comparability there is a need for methods to match different PET���CT systems through elimination of this systemic variability. Consequently, a new method was developed and tested that transforms the image of an object as produced by one tomograph to another image of the same object as it would have been seen by a different tomograph. The proposed new method, termed Transconvolution, compensates for differing imaging properties of different tomographs and particularly aims at quantitative comparability of PET���CT in the context of multicenter trials. METHODS To solve the problem of image normalization, the theory of Transconvolution was mathematically established together with new methods to handle point spread functions of different PET���CT systems. Knowing the point spread functions of two different imaging systems allows determining a Transconvolution function to convert one image into the other. This function is calculated by convolving one point spread function with the inverse of the other point spread function which, when adhering to certain boundary conditions such as the use of linear acquisition and image reconstruction methods, is a numerically accessible operation. For reliable measurement of such point spread functions characterizing different PET���CT systems, a dedicated solid-state phantom incorporating (68)Ge���(68)Ga filled spheres was developed. To iteratively determine and represent such point spread functions, exponential density functions in combination with a Gaussian distribution were introduced. Furthermore, simulation of a virtual PET system provided a standard imaging system with clearly defined properties to which the real PET systems were to be matched. A Hann window served as the modulation transfer function for the virtual PET. The Hann's apodization properties suppressed high spatial frequencies above a certain critical frequency, thereby fulfilling the above-mentioned boundary conditions. The determined point spread functions were subsequently used by the novel Transconvolution algorithm to match different PET���CT systems onto the virtual PET system. Finally, the theoretically elaborated Transconvolution method was validated transforming phantom images acquired on two different PET systems to nearly identical data sets, as they would be imaged by the virtual PET system. RESULTS The proposed Transconvolution method matched different PET���CT-systems for an improved and reproducible determination of a normalized activity concentration. The highest difference in measured activity concentration between the two different PET systems of 18.2% was found in spheres of 2 ml volume. Transconvolution reduced this difference down to 1.6%. In addition to reestablishing comparability the new method with its parameterization of point spread functions allowed a full characterization of imaging properties of the examined tomographs. CONCLUSIONS By matching different tomographs to a virtual standardized imaging system, Transconvolution opens a new comprehensive method for cross calibration in quantitative PET imaging. The use of a virtual PET system restores comparability between data sets from different PET systems by exerting a common, reproducible, and defined partial volume effect.

Role of A2B adenosine receptor signaling in adenosine-dependent pulmonary inflammation and injury.

Adenosine has been implicated in the pathogenesis of chronic lung diseases such as asthma and chronic obstructive pulmonary disease. In vitro studies suggest that activation of the A2B adenosine receptor (A2BAR) results in proinflammatory and profibrotic effects relevant to the progression of lung diseases; however, in vivo data supporting these observations are lacking. Adenosine deaminase-deficient (ADA-deficient) mice develop pulmonary inflammation and injury that are dependent on increased lung adenosine levels. To investigate the role of the A2BAR in vivo, ADA-deficient mice were treated with the selective A2BAR antagonist CVT-6883, and pulmonary inflammation, fibrosis, and airspace integrity were assessed. Untreated and vehicle-treated ADA-deficient mice developed pulmonary inflammation, fibrosis, and enlargement of alveolar airspaces; conversely, CVT-6883-treated ADA-deficient mice showed less pulmonary inflammation, fibrosis, and alveolar airspace enlargement. A2BAR antagonism significantly reduced elevations in proinflammatory cytokines and chemokines as well as mediators of fibrosis and airway destruction. In addition, treatment with CVT-6883 attenuated pulmonary inflammation and fibrosis in wild-type mice subjected to bleomycin-induced lung injury. These findings suggest that A2BAR signaling influences pathways critical for pulmonary inflammation and injury in vivo. Thus in chronic lung diseases associated with increased adenosine, antagonism of A2BAR-mediated responses may prove to be a beneficial therapy.

IN-VIVO CT DOSIMETRY DURING VIRTUAL COLONOSCOPY

Virtual colonoscopy (VC) is a minimally invasive means for identifying colorectal polyps and colorectal lesions by insufflating a patient���s bowel, applying contrast agent via rectal catheter, and performing multi-detector computed tomography (MDCT) scans. The technique is recommended for colonic health screening by the American Cancer Society but not funded by the Centers for Medicare and Medicaid Services (CMS) partially because of potential risks from radiation exposure. To date, no in���vivo organ dose measurements have been performed for MDCT scans; thus, the accuracy of any current dose estimates is currently unknown. In this study, two TLDs were affixed to the inner lumen of standard rectal catheters used in VC, and in-vivo rectal dose measurements were obtained within 6 VC patients. In order to calculate rectal dose, TLD-100 powder response was characterized at diagnostic doses such that appropriate correction factors could be determined for VC. A third-order polynomial regression with a goodness of fit factor of R2=0.992 was constructed from this data. Rectal dose measurements were acquired with TLDs during simulated VC within a modified anthropomorphic phantom configured to represent three sizes of patients undergoing VC. The measured rectal doses decreased in an exponential manner with increasing phantom effective diameter, with R2=0.993 for the exponential regression model and a maximum percent coefficient of variation (%CoV) of 4.33%. In-vivo measurements yielded rectal doses ranged from that decreased exponentially with increasing patient effective diameter, in a manner that was also favorably predicted by the size specific dose estimate (SSDE) model for all VC patients that were of similar age, body composition, and TLD placement. The measured rectal dose within a younger patient was favorably predicted by the anthropomorphic phantom dose regression model due to similarities in the percentages of highly attenuating material at the respective measurement locations and in the placement of the TLDs. The in-vivo TLD response did not increase in %CoV with decreasing dose, and the largest %CoV was 10.0%.

Evaluation of deformable image registration for improved 4D CT-derived ventilation for image guided radiotherapy

Recent treatment planning studies have demonstrated the use of physiologic images in radiation therapy treatment planning to identify regions for functional avoidance. This image-guided radiotherapy (IGRT) strategy may reduce the injury and/or functional loss following thoracic radiotherapy. 4D computed tomography (CT), developed for radiotherapy treatment planning, is a relatively new imaging technique that allows the acquisition of a time-varying sequence of 3D CT images of the patient's lungs through the respiratory cycle. Guerrero et al. developed a method to calculate ventilation imaging from 4D CT, which is potentially better suited and more broadly available for IGRT than the current standard imaging methods. The key to extracting function information from 4D CT is the construction of a volumetric deformation field that accurately tracks the motion of the patient's lungs during the respiratory cycle. The spatial accuracy of the displacement field directly impacts the ventilation images; higher spatial registration accuracy will result in less ventilation image artifacts and physiologic inaccuracies. Presently, a consistent methodology for spatial accuracy evaluation of the DIR transformation is lacking. Evaluation of the 4D CT-derived ventilation images will be performed to assess correlation with global measurements of lung ventilation, as well as regional correlation of the distribution of ventilation with the current clinical standard SPECT. This requires a novel framework for both the detailed assessment of an image registration algorithm's performance characteristics as well as quality assurance for spatial accuracy assessment in routine application. Finally, we hypothesize that hypo-ventilated regions, identified on 4D CT ventilation images, will correlate with hypo-perfused regions in lung cancer patients who have obstructive lesions. A prospective imaging trial of patients with locally advanced non-small-cell lung cancer will allow this hypothesis to be tested. These advances are intended to contribute to the validation and clinical implementation of CT-based ventilation imaging in prospective clinical trials, in which the impact of this imaging method on patient outcomes may be tested.

The association of genetic groups 1, 2, and 3 of Mycobacterium tuberculosis with pulmonary and extra-pulmonary tuberculosis

Objective. To investigate the association of the three major genetic groups of Mycobacterium tuberculosis with pulmonary and extra-pulmonary tuberculosis in clustered and non-clustered TB cases in the Houston area. ^ Study design. Secondary analysis of an ambi-directional study. ^ Study population. Three hundred fifty-eight confirmed cases of tuberculosis in the Houston that occurred between October 1995 and May 1997, who had been interviewed by the Houston T13 Initiative staff at Baylor College of Medicine, and whose isolates have had their DNA fingerprint and genetic group determined. ^ Exclusions. Individuals whose mycobacterial genotype was unknown, or whose data variables were unavailable. ^ Source of data. Laboratory results, patient interviews, and medical records at clinics and hospitals of the study population. ^ Results. In clustered cases, the majority of both, pulmonary and extra-pulmonary TB cases were caused by genetic group 1. Independent factors were assessed to determine the interactions that may influence the site of infection or increase the risk for one site or another. HIV negative males were protected against extra-pulmonary TB compared to HIV negative females. Individuals ages 1���14 years were at higher risk of having extra-pulmonary TB. Group 3 organisms were found less frequently in the total population in general, especially in extra-pulmonary disease. This supports the evidence in previous studies that this group is the least virulent and genetically distinct from the other two groups. Group 1 was found more frequently among African Americans than other ethnic groups, a trend for future investigations. ^ Among the non-clustered cases, group 2 organisms were the majority of the organisms found in both sites. They were also the majority of organisms found in African Americans, Caucasians, and Hispanics causing the majority of the infections at both sites. However, group 1 organisms were the overwhelming majority found in Asian/Pacific Islander individuals, which may indicate these organisms are either endemic to that area, or that there is an ethnic biological factor involved. This may also be due to a systematic bias, since isolates from individuals from that geographic region lack adequate copies of the insertion sequence IS6110, which leads to their placement in the non-clustered population. ^ The three genetic groups of Mycobacterium tuberculosis were not found equally distributed between sites of infection in both clustered and non-clustered cases. Furthermore, these groups were not distributed in the same patterns among the clustered and non-clustered cases, but rather in distinct patterns. ^

Ultrasound-assisted thrombolysis for acute pulmonary embolism: a systematic review.

Pulmonary embolism remains a common and potentially life-threatening disease. For patients with intermediate- and high-risk pulmonary embolism, catheter-based revascularization therapy has emerged as potential alternative to systemic thrombolysis or surgical embolectomy. Ultrasound-assisted catheter-directed thrombolysis is a contemporary catheter-based technique and is the focus of the present review. Ultrasound-assisted catheter-directed thrombolysis is more effective in reversing right ventricular dysfunction and dilatation in comparison with anticoagulation alone in patients at intermediate risk. However, a direct comparison of ultrasound-assisted thrombolysis with systemic thrombolysis or surgical thrombectomy is not available. Ultrasound-assisted thrombolysis with initial intrapulmonary thrombolytic bolus may also be effective in high-risk patients, but evidence from randomized trials is not available. This review summarizes current data on ultrasound-assisted thrombolysis for acute pulmonary embolism.

Assessment of Transverse Isotropy in Clinical-Level CT Images of Trabecular Bone Using the Gradient Structure Tensor

The aim of this study was to develop a GST-based methodology for accurately measuring the degree of transverse isotropy in trabecular bone. Using femoral sub-regions scanned in high-resolution peripheral QCT (HR-pQCT) and clinical-level-resolution QCT, trabecular orientation was evaluated using the mean intercept length (MIL) and the gradient structure tensor (GST) on the HR-pQCT and QCT data, respectively. The influence of local degree of transverse isotropy (DTI) and bone mineral density (BMD) was incorporated into the investigation. In addition, a power based model was derived, rendering a 1:1 relationship between GST and MIL eigenvalues. A specific DTI threshold (DTI thres) was found for each investigated size of region of interest (ROI), above which the estimate of major trabecular direction of the GST deviated no more than 30�� from the gold standard MIL in 95% of the remaining ROIs (mean error: 16��). An inverse relationship between ROI size and DTI thres was found for discrete ranges of BMD. A novel methodology has been developed, where transversal isotropic measures of trabecular bone can be obtained from clinical QCT images for a given ROI size, DTI thres and power coefficient. Including DTI may improve future clinical QCT finite-element predictions of bone strength and diagnoses of bone disease.

The effects of cardiac output and pulmonary arterial hypertension on volumetric capnography derived-variables during normoxia and hypoxia

The aim of this study was to test the effect of cardiac output (CO) and pulmonary artery hypertension (PHT) on volumetric capnography (VCap) derived-variables. Nine pigs were mechanically ventilated using fixed ventilatory settings. Two steps of PHT were induced by IV infusion of a thromboxane analogue: PHT25 [mean pulmonary arterial pressure (MPAP) of 25��mmHg] and PHT40 (MPAP of 40��mmHg). CO was increased by 50��% from baseline (COup) with an infusion of dobutamine ���5����g��kg(-1)��min(-1) and decreased by 40��% from baseline (COdown) infusing sodium nitroglycerine ���30����g��kg(-1)��min(-1) plus esmolol 500����g��kg(-1)��min(-1). Another state of PHT and COdown was induced by severe hypoxemia (FiO2 0.07). Invasive hemodynamic data and VCap were recorded and compared before and after each step using a mixed random effects model. Compared to baseline, the normalized slope of phase III (SnIII) increased by 32��% in PHT25 and by 22��% in PHT40. SnIII decreased non-significantly by 4��% with COdown. A combination of PHT and COdown associated with severe hypoxemia increased SnIII by 28��% compared to baseline. The elimination of CO2 per breath decreased by 7��% in PHT40 and by 12��% in COdown but increased only slightly with COup. Dead space variables did not change significantly along the protocol. At constant ventilation and body metabolism, pulmonary artery hypertension and decreases in CO had the biggest effects on the SnIII of the volumetric capnogram and on the elimination of CO2.

Multi-phase postmortem CT angiography: recognizing technique-related artefacts and pitfalls

BACKGROUND AND PURPOSE Multi-phase postmortem CT angiography (MPMCTA) is increasingly being recognized as a valuable adjunct medicolegal tool to explore the vascular system. Adequate interpretation, however, requires knowledge about the most common technique-related artefacts. The purpose of this study was to identify and index the possible artefacts related to MPMCTA. MATERIAL AND METHODS An experienced radiologist blinded to all clinical and forensic data retrospectively reviewed 49 MPMCTAs. Each angiographic phase, i.e. arterial, venous and dynamic, was analysed separately to identify phase-specific artefacts based on location and aspect. RESULTS Incomplete contrast filling of the cerebral venous system was the most commonly encountered artefact, followed by contrast agent layering in the lumen of the thoracic aorta. Enhancement or so-called oedematization of the digestive system mucosa was also frequently observed. CONCLUSION All MPMCTA artefacts observed and described here are reproducible and easily identifiable. Knowledge about these artefacts is important to avoid misinterpreting them as pathological findings.

Pulmonary aspergilloma: A rare differential diagnosis to lung cancer after positive FDG PET scan

Early diagnosis and treatment of lung cancer, one of the leading causes of cancer-related death, is important to improve morbidity and mortality. Therefore any suspect solitary pulmonary nodule should prompt the pursuit for a definitive histological diagnosis. We describe the case of a 55-years-old male ex-smoker, who was admitted to our hospital due to recurrent hemoptysis and dry cough. A CT scan showed an irregular nodule of increasing size (28 mm in diameter) in the left lower lobe (LLL). A whole body PET-CT scan (643 MBq F-18 FDG i.v.) was performed and confirmed an avid FDG uptake of the nodule in the LLL, highly suspicious of lung cancer, without any evidence of lymphogenic or hematogenic metastasis. Bronchoscopy was not diagnostic and due to severe adhesions after prior chest trauma and the central location of the nodule, a lobectomy of the LLL was performed. Surprisingly, histology showed a simple aspergilloma located in a circumscribed bronchiectasis with no evidence of malignancy. This is a report of an informative example of an aspergilloma, which presented with symptoms and radiological features of malignant lung cancer.

Reducing CT radiation dose with iterative reconstruction algorithms: the influence of scan and reconstruction parameters on image quality and CTDIvol

OBJECTIVES In this phantom CT study, we investigated whether images reconstructed using filtered back projection (FBP) and iterative reconstruction (IR) with reduced tube voltage and current have equivalent quality. We evaluated the effects of different acquisition and reconstruction parameter settings on image quality and radiation doses. Additionally, patient CT studies were evaluated to confirm our phantom results. METHODS Helical and axial 256 multi-slice computed tomography scans of the phantom (Catphan(��)) were performed with varying tube voltages (80-140kV) and currents (30-200mAs). 198 phantom data sets were reconstructed applying FBP and IR with increasing iterations, and soft and sharp kernels. Further, 25 chest and abdomen CT scans, performed with high and low exposure per patient, were reconstructed with IR and FBP. Two independent observers evaluated image quality and radiation doses of both phantom and patient scans. RESULTS In phantom scans, noise reduction was significantly improved using IR with increasing iterations, independent from tissue, scan-mode, tube-voltage, current, and kernel. IR did not affect high-contrast resolution. Low-contrast resolution was also not negatively affected, but improved in scans with doses <5mGy, although object detectability generally decreased with the lowering of exposure. At comparable image quality levels, CTDIvol was reduced by 26-50% using IR. In patients, applying IR vs. FBP resulted in good to excellent image quality, while tube voltage and current settings could be significantly decreased. CONCLUSIONS Our phantom experiments demonstrate that image quality levels of FBP reconstructions can also be achieved at lower tube voltages and tube currents when applying IR. Our findings could be confirmed in patients revealing the potential of IR to significantly reduce CT radiation doses.

Volumetric analysis of lung nodules in computed tomography (CT): comparison of two different segmentation algorithm softwares and two different reconstruction filters on automated volume calculation.

BACKGROUND A precise detection of volume change allows for better estimating the biological behavior of the lung nodules. Postprocessing tools with automated detection, segmentation, and volumetric analysis of lung nodules may expedite radiological processes and give additional confidence to the radiologists. PURPOSE To compare two different postprocessing software algorithms (LMS Lung, Median Technologies; LungCARE��, Siemens) in CT volumetric measurement and to analyze the effect of soft (B30) and hard reconstruction filter (B70) on automated volume measurement. MATERIAL AND METHODS Between January 2010 and April 2010, 45 patients with a total of 113 pulmonary nodules were included. The CT exam was performed on a 64-row multidetector CT scanner (Somatom Sensation, Siemens, Erlangen, Germany) with the following parameters: collimation, 24x1.2���mm; pitch, 1.15; voltage, 120���kVp; reference tube current-time, 100���mAs. Automated volumetric measurement of each lung nodule was performed with the two different postprocessing algorithms based on two reconstruction filters (B30 and B70). The average relative volume measurement difference (VME%) and the limits of agreement between two methods were used for comparison. RESULTS At soft reconstruction filters the LMS system produced mean nodule volumes that were 34.1% (P���<���0.0001) larger than those by LungCARE�� system. The VME% was 42.2% with a limit of agreement between -53.9% and 138.4%.The volume measurement with soft filters (B30) was significantly larger than with hard filters (B70); 11.2% for LMS and 1.6% for LungCARE��, respectively (both with P���<���0.05). LMS measured greater volumes with both filters, 13.6% for soft and 3.8% for hard filters, respectively (P���<���0.01 and P���>���0.05). CONCLUSION There is a substantial inter-software (LMS/LungCARE��) as well as intra-software variability (B30/B70) in lung nodule volume measurement; therefore, it is mandatory to use the same equipment with the same reconstruction filter for the follow-up of lung nodule volume.

Fusion of fundus images and MRI data of the human eye

Purpose Ophthalmologists are confronted with a set of different image modalities to diagnose eye tumors e.g., fundus photography, CT and MRI. However, these images are often complementary and represent pathologies differently. Some aspects of tumors can only be seen in a particular modality. A fusion of modalities would improve the contextual information for diagnosis. The presented work attempts to register color fundus photography with MRI volumes. This would complement the low resolution 3D information in the MRI with high resolution 2D fundus images. Methods MRI volumes were acquired from 12 infants under the age of 5 with unilateral retinoblastoma. The contrast-enhanced T1-FLAIR sequence was performed with an isotropic resolution of less than 0.5mm. Fundus images were acquired with a RetCam camera. For healthy eyes, two landmarks were used: the optic disk and the fovea. The eyes were detected and extracted from the MRI volume using a 3D adaption of the Fast Radial Symmetry Transform (FRST). The cropped volume was automatically segmented using the Split Bregman algorithm. The optic nerve was enhanced by a Frangi vessel filter. By intersection the nerve with the retina the optic disk was found. The fovea position was estimated by constraining the position with the angle between the optic and the visual axis as well as the distance from the optic disk. The optical axis was detected automatically by fitting a parable on to the lens surface. On the fundus, the optic disk and the fovea were detected by using the method of Budai et al. Finally, the image was projected on to the segmented surface using the lens position as the camera center. In tumor affected eyes, the manually segmented tumors were used instead of the optic disk and macula for the registration. Results In all of the 12 MRI volumes that were tested the 24 eyes were found correctly, including healthy and pathological cases. In healthy eyes the optic nerve head was found in all of the tested eyes with an error of 1.08 +/- 0.37mm. A successful registration can be seen in figure 1. Conclusions The presented method is a step toward automatic fusion of modalities in ophthalmology. The combination enhances the MRI volume with higher resolution from the color fundus on the retina. Tumor treatment planning is improved by avoiding critical structures and disease progression monitoring is made easier.