3D volumetric intensity reconsturction from 2D x-ray images using partial least squares regression

Reconstruction of shape and intensity from 2D x-ray images has drawn more and more attentions. Previously introduced work suffers from the long computing time due to its iterative optimization characteristics and the requirement of generating digitally reconstructed radiographs within each iteration. In this paper, we propose a novel method which uses a patient-specific 3D surface model reconstructed from 2D x-ray images as a surrogate to get a patient-specific volumetric intensity reconstruction via partial least squares regression. No DRR generation is needed. The method was validated on 20 cadaveric proximal femurs by performing a leave-one-out study. Qualitative and quantitative results demonstrated the efficacy of the present method. Compared to the existing work, the present method has the advantage of much shorter computing time and can be applied to both DXA images as well as conventional x-ray images, which may hold the potentials to be applied to clinical routine task such as total hip arthroplasty (THA).

Investigating the sensitivity of hurricane intensity and trajectory to sea surface temperatures using the regional model WRF

The influence of sea surface temperature (SST) anomalies on the hurricane characteristics are investigated in a set of sensitivity experiments employing the Weather Research and Forecasting (WRF) model. The idealised experiments are performed for the case of Hurricane Katrina in 2005. The first set of sensitivity experiments with basin-wide changes of the SST magnitude shows that the intensity goes along with changes in the SST, i.e., an increase in SST leads to an intensification of Katrina. Additionally, the trajectory is shifted to the west (east), with increasing (decreasing) SSTs. The main reason is a strengthening of the background flow. The second set of experiments investigates the influence of Loop Current eddies idealised by localised SST anomalies. The intensity of Hurricane Katrina is enhanced with increasing SSTs close to the core of a tropical cyclone. Negative nearby SST anomalies reduce the intensity. The trajectory only changes if positive SST anomalies are located west or north of the hurricane centre. In this case the hurricane is attracted by the SST anomaly which causes an additional moisture source and increased vertical winds.

Thoracic target volume delineation using various maximum-intensity projection computed tomography image sets for stereotactic body radiation therapy

The motion of lung tumors during respiration makes the accurate delivery of radiation therapy to the thorax difficult because it increases the uncertainty of target position. The adoption of four-dimensional computed tomography (4D-CT) has allowed us to determine how a tumor moves with respiration for each individual patient. Using information acquired during a 4D-CT scan, we can define the target, visualize motion, and calculate dose during the planning phase of the radiotherapy process. One image data set that can be created from the 4D-CT acquisition is the maximum-intensity projection (MIP). The MIP can be used as a starting point to define the volume that encompasses the motion envelope of the moving gross target volume (GTV). Because of the close relationship that exists between the MIP and the final target volume, we investigated four MIP data sets created with different methodologies (3 using various 4D-CT sorting implementations, and one using all available cine CT images) to compare target delineation. It has been observed that changing the 4D-CT sorting method will lead to the selection of a different collection of images; however, the clinical implications of changing the constituent images on the resultant MIP data set are not clear. There has not been a comprehensive study that compares target delineation based on different 4D-CT sorting methodologies in a patient population. We selected a collection of patients who had previously undergone thoracic 4D-CT scans at our institution, and who had lung tumors that moved at least 1 cm. We then generated the four MIP data sets and automatically contoured the target volumes. In doing so, we identified cases in which the MIP generated from a 4D-CT sorting process under-represented the motion envelope of the target volume by more than 10% than when measured on the MIP generated from all of the cine CT images. The 4D-CT methods suffered from duplicate image selection and might not choose maximum extent images. Based on our results, we suggest utilization of a MIP generated from the full cine CT data set to ensure a representative inclusive tumor extent, and to avoid geometric miss.

Is a 3-mm intrafractional margin sufficient for daily image-guided intensity-modulated radiation therapy of prostate cancer?

PURPOSE: To determine whether a 3-mm isotropic target margin adequately covers the prostate and seminal vesicles (SVs) during administration of an intensity-modulated radiation therapy (IMRT) treatment fraction, assuming that daily image-guided setup is performed just before each fraction. MATERIALS AND METHODS: In-room computed tomographic (CT) scans were acquired immediately before and after a daily treatment fraction in 46 patients with prostate cancer. An eight-field IMRT plan was designed using the pre-fraction CT with a 3-mm margin and subsequently recalculated on the post-fraction CT. For convenience of comparison, dose plans were scaled to full course of treatment (75.6 Gy). Dose coverage was assessed on the post-treatment CT image set. RESULTS: During one treatment fraction (21.4+/-5.5 min), there were reductions in the volumes of the prostate and SVs receiving the prescribed dose (median reduction 0.1% and 1.0%, respectively, p<0.001) and in the minimum dose to 0.1 cm(3) of their volumes (median reduction 0.5 and 1.5 Gy, p<0.001). Of the 46 patients, three patients' prostates and eight patients' SVs did not maintain dose coverage above 70 Gy. Rectal filling correlated with decreased percentage-volume of SV receiving 75.6, 70, and 60 Gy (p<0.02). CONCLUSIONS: The 3-mm intrafractional margin was adequate for prostate dose coverage. However, a significant subset of patients lost SV dose coverage. The rectal volume change significantly affected SV dose coverage. For advanced-stage prostate cancers, we recommend to use larger margins or improve organ immobilization (such as with a rectal balloon) to ensure SV coverage.

Ultraviolet radiation intensity predicts the relative distribution of dermatomyositis and anti-Mi-2 autoantibodies in women.

OBJECTIVE: Because studies suggest that ultraviolet (UV) radiation modulates the myositis phenotype and Mi-2 autoantigen expression, we conducted a retrospective investigation to determine whether UV radiation may influence the relative prevalence of dermatomyositis and anti-Mi-2 autoantibodies in the US. METHODS: We assessed the relationship between surface UV radiation intensity in the state of residence at the time of onset with the relative prevalence of dermatomyositis and myositis autoantibodies in 380 patients with myositis from referral centers in the US. Myositis autoantibodies were detected by validated immunoprecipitation assays. Surface UV radiation intensity was estimated from UV Index data collected by the US National Weather Service. RESULTS: UV radiation intensity was associated with the relative proportion of patients with dermatomyositis (odds ratio [OR] 2.3, 95% confidence interval [95% CI] 0.9-5.8) and with the proportion of patients expressing anti-Mi-2 autoantibodies (OR 6.0, 95% CI 1.1-34.1). Modeling of these data showed that these associations were confined to women (OR 3.8, 95% CI 1.3-11.0 and OR 17.3, 95% CI 1.8-162.4, respectively) and suggests that sex influences the effects of UV radiation on autoimmune disorders. Significant associations were not observed in men, nor were UV radiation levels related to the presence of antisynthetase or anti-signal recognition particle autoantibodies. CONCLUSION: This first study of the distribution of myositis phenotypes and UV radiation exposure in the US showed that UV radiation may modulate the clinical and immunologic expression of autoimmune disease in women. Further investigation of the mechanisms by which these effects are produced may provide insights into pathogenesis and suggest therapeutic or preventative strategies.

The Galactic Cosmic Ray Intensity over the Past 106–109 Years as Recorded by Cosmogenic Nuclides in Meteorites and Terrestrial Samples

Concentrations of stable and radioactive nuclides produced by cosmic ray particles in meteorites allow us to track the long term average of the primary flux of galactic cosmic rays (GCR). During the past ∼10 Ma, the average GCR flux remained constant over timescales of hundreds of thousands to millions of years, and, if corrected for known variations in solar modulation, also during the past several years to hundreds of years. Because the cosmic ray concentrations in meteorites represent integral signals, it is difficult to assess the limits of uncertainty of this statement, but they are larger than the often quoted analytical and model uncertainties of some 30%. Time series of concentrations of the radionuclide 10Be in terrestrial samples strengthen the conclusions drawn from meteorite studies, indicating that the GCR intensity on a ∼0.5 million year scale has remained constant within some ±10% during the past ∼10 million years. The very long-lived radioactive nuclide 40K allows to assess the GCR flux over about the past one billion years. The flux over the past few million years has been the same as the longer-term average in the past 0.5–1 billion years within a factor of ∼1.5. However, newer data do not confirm a long-held belief that the flux in the past few million years has been higher by some 30–50% than the very long term average. Neither does our analysis confirm a hypothesis that the iron meteorite data indicate a ∼150 million year periodicity in the cosmic ray flux, possibly related to variations in the long-term terrestrial climate.

Assessment of Collimator Jaw Optimization in Reducing Normal Tissue Irradiation with Intensity Modulated Radiation Therapy

Purpose: To evaluate normal tissue dose reduction in step-and-shoot intensity-modulated radiation therapy (IMRT) on the Varian 2100 platform by tracking the multileaf collimator (MLC) apertures with the accelerator jaws. Methods: Clinical radiation treatment plans for 10 thoracic, 3 pediatric and 3 head and neck patients were converted to plans with the jaws tracking each segment’s MLC apertures. Each segment was then renormalized to account for the change in collimator scatter to obtain target coverage within 1% of that in the original plan. The new plans were compared to the original plans in a commercial radiation treatment planning system (TPS). Reduction in normal tissue dose was evaluated in the new plan by using the parameters V5, V10, and V20 in the cumulative dose-volume histogram for the following structures: total lung minus GTV (gross target volume), heart, esophagus, spinal cord, liver, parotids, and brainstem. In order to validate the accuracy of our beam model, MLC transmission measurements were made and compared to those predicted by the TPS. Results: The greatest change between the original plan and new plan occurred at lower dose levels. The reduction in V20 was never more than 6.3% and was typically less than 1% for all patients. The reduction in V5 was 16.7% maximum and was typically less than 3% for all patients. The variation in normal tissue dose reduction was not predictable, and we found no clear parameters that indicated which patients would benefit most from jaw tracking. Our TPS model of MLC transmission agreed with measurements with absolute transmission differences of less than 0.1 % and thus uncertainties in the model did not contribute significantly to the uncertainty in the dose determination. Conclusion: The amount of dose reduction achieved by collimating the jaws around each MLC aperture in step-and-shoot IMRT does not appear to be clinically significant.

EVALUATION OF INTENSITY MODULATED RADIATION THERAPY (IMRT) DELIVERY ERROR DUE TO IMRT TREATMENT PLAN COMPLEXITY AND IMPROPERLY MATCHED DOSIMETRY DATA

Intensity modulated radiation therapy (IMRT) is a technique that delivers a highly conformal dose distribution to a target volume while attempting to maximally spare the surrounding normal tissues. IMRT is a common treatment modality used for treating head and neck (H&N) cancers, and the presence of many critical structures in this region requires accurate treatment delivery. The Radiological Physics Center (RPC) acts as both a remote and on-site quality assurance agency that credentials institutions participating in clinical trials. To date, about 30% of all IMRT participants have failed the RPC’s remote audit using the IMRT H&N phantom. The purpose of this project is to evaluate possible causes of H&N IMRT delivery errors observed by the RPC, specifically IMRT treatment plan complexity and the use of improper dosimetry data from machines that were thought to be matched but in reality were not. Eight H&N IMRT plans with a range of complexity defined by total MU (1460-3466), number of segments (54-225), and modulation complexity scores (MCS) (0.181-0.609) were created in Pinnacle v.8m. These plans were delivered to the RPC’s H&N phantom on a single Varian Clinac. One of the IMRT plans (1851 MU, 88 segments, and MCS=0.469) was equivalent to the median H&N plan from 130 previous RPC H&N phantom irradiations. This average IMRT plan was also delivered on four matched Varian Clinac machines and the dose distribution calculated using a different 6MV beam model. Radiochromic film and TLD within the phantom were used to analyze the dose profiles and absolute doses, respectively. The measured and calculated were compared to evaluate the dosimetric accuracy. All deliveries met the RPC acceptance criteria of ±7% absolute dose difference and 4 mm distance-to-agreement (DTA). Additionally, gamma index analysis was performed for all deliveries using a ±7%/4mm and ±5%/3mm criteria. Increasing the treatment plan complexity by varying the MU, number of segments, or varying the MCS resulted in no clear trend toward an increase in dosimetric error determined by the absolute dose difference, DTA, or gamma index. Varying the delivery machines as well as the beam model (use of a Clinac 6EX 6MV beam model vs. Clinac 21EX 6MV model), also did not show any clear trend towards an increased dosimetric error using the same criteria indicated above.

Patient classification system: An integrated method for measuring nursing intensity and optimizing resource allocation

A patient classification system was developed integrating a patient acuity instrument with a computerized nursing distribution method based on a linear programming model. The system was designed for real-time measurement of patient acuity (workload) and allocation of nursing personnel to optimize the utilization of resources.^ The acuity instrument was a prototype tool with eight categories of patients defined by patient severity and nursing intensity parameters. From this tool, the demand for nursing care was defined in patient points with one point equal to one hour of RN time. Validity and reliability of the instrument was determined as follows: (1) Content validity by a panel of expert nurses; (2) predictive validity through a paired t-test analysis of preshift and postshift categorization of patients; (3) initial reliability by a one month pilot of the instrument in a practice setting; and (4) interrater reliability by the Kappa statistic.^ The nursing distribution system was a linear programming model using a branch and bound technique for obtaining integer solutions. The objective function was to minimize the total number of nursing personnel used by optimally assigning the staff to meet the acuity needs of the units. A penalty weight was used as a coefficient of the objective function variables to define priorities for allocation of staff.^ The demand constraints were requirements to meet the total acuity points needed for each unit and to have a minimum number of RNs on each unit. Supply constraints were: (1) total availability of each type of staff and the value of that staff member (value was determined relative to that type of staff's ability to perform the job function of an RN (i.e., value for eight hours RN = 8 points, LVN = 6 points); (2) number of personnel available for floating between units.^ The capability of the model to assign staff quantitatively and qualitatively equal to the manual method was established by a thirty day comparison. Sensitivity testing demonstrated appropriate adjustment of the optimal solution to changes in penalty coefficients in the objective function and to acuity totals in the demand constraints.^ Further investigation of the model documented: correct adjustment of assignments in response to staff value changes; and cost minimization by an addition of a dollar coefficient to the objective function. ^

An Ice-Core Proxy for Antarctic Circumpolar Zonal Wind Intensity

Using US National Centers for Environmental Prediction/US National Center for Atmospheric Research re-analysis data, we investigate the relationships between crustal ion (nssCa(2+)) concentrations from three West Antarctic ice cores, namely, Siple Dome (SD), ITASE00-1 (IT001) and ITASE01-5 (IT015), and primary components of the climate system, namely, air pressure/geopotential height, zonal (u) and meridional (v) wind strength. Linear correlation analyses between nssCa(2+) concentrations and both air-pressure and wind fields for the period of overlap between records indicate that the SD nssCa(2+) variation is positively correlated with spring circumpolar zonal wind, while IT001 nssCa(2+) has a positive correlation with circumpolar zonal wind throughout the year (r > 0.3, p < 0.01). Intensified Southern Westerlies circulation is conducive to transport of more crustal aerosols to both sites. Further correlation analyses between nssCa(2+) concentrations from SD and IT001 and atmospheric circulation suggest that the high inland plateau (represented by core IT001) is largely influenced by transport from the upper troposphere. IT015 nssCa(2+) is negatively correlated with westerly wind in October and November, suggesting that stronger westerly circulation may weaken the transport of crustal species to IT015. Correlations of nssCa(2+) from the three ice cores with the Antarctic Oscillation index are consistent with results developed from the wind-field investigation. In addition, calibration between nssCa(2+) concentration and the multivariate El Nino-Southern Oscillation (ENSO) index shows that crustal species transport to IT001 is enhanced during strong ENSO events.

Resistance of a lowland rain forest to increasing drought intensity in Sabah, Borneo

Occasional strong droughts are an important feature of the climatic environment of tropical rain forest in much of Borneo. This paper compares the response of a lowland dipterocarp forest at Danum, Sabah, in a period of low (LDI) and a period of high (HDI) drought intensity (1986-96, 9.98 y;1996-99, 2.62 y). Mean annual drought intensity was two-fold higher in the HDI than LDI period (1997 v. 976 mm), and each period had one moderately strong main drought (viz. 1992, 1998). Mortality of `all' trees greater than or equal to 10 cm gbh (girth at breast height) and stem growth rates of `small' trees 10less than or equal to50 cm gbh were measured in sixteen 0.16-ha subplots (half on ridge, half on lower slope sites) within two 4-ha plots. These 10-50-cm trees were composed largely of true understorey species. A new procedure was developed to correct for the effect of differences in length of census interval when comparing tree mortality rates. Mortality rates of small trees declined slightly but not significantly between the LDI and HDI periods (1.53 to 1.48% y(-1)): mortality of all trees showed a similar pattern. Relative growth rates declined significantly by 23% from LDI to HDI periods (11.1 to 8.6 mm m(-1) y(-1)): for absolute growth rates the decrease was 28% (2.45 to 1.77 mm y(-1)). Neither mortality nor growth rates were significantly influenced by topography. For small trees, across subplots, absolute growth rate was positively correlated in the LDI period, but negatively correlated in the HDI period, with mortality rate. There was no consistent pattern in the responses among the 19 most abundant species (n greater than or equal to 50 trees) which included a proposed drought-tolerant guild. In terms of tree survival, the forest at Danum was resistant to increasing drought intensity, but showed decreased stem growth attributable to increasing water stress.

Effects of Thoratec pulsatile ventricular assist device (PVAD) timing on the abdominal aortic wave intensity pattern

Arterial waves are seen as possible independent mediators of cardiovascular risks, and the wave intensity analysis (WIA) has therefore been proposed as a method for patient selection for ventricular assist device (VAD) implantation. Interpreting measured wave intensity (WI) is challenging and complexity is increased by the implantation of a VAD. The waves generated by the VAD interact with the waves generated by the native heart, and this interaction varies with changing VAD settings. Eight sheep were implanted with a pulsatile VAD (PVAD) through ventriculo-aortic cannulation. The start of PVAD ejection was synchronized to the native R-wave and delayed between 0 % - 90 % of the cardiac cycle in 10 % steps or phase shifts (PS). Pressure and velocity signals were registered, using a combined Doppler and pressure wire positioned in the abdominal aorta, and used to calculate the WI. Depending on the PS, different wave interference phenomena occurred. Maximum unloading of the left ventricle (LV) coincided with constructive interference and maximum blood flow pulsatility, and maximum loading of the LV coincided with destructive interference and minimum blood flow pulsatility. We believe, that non-invasive WIA could potentially be used clinically to assess the mechanical load of the LV, and to monitor the peripheral hemodynamics such as blood flow pulsatility and risk of intestinal bleeding.

Effect of high-intensity statin therapy on atherosclerosis in non-infarct-related coronary arteries (IBIS-4): a serial intravascular ultrasonography study

AIM The effect of long-term high-intensity statin therapy on coronary atherosclerosis among patients with acute ST-segment elevation myocardial infarction (STEMI) is unknown. The aim of this study was to quantify the impact of high-intensity statin therapy on plaque burden, composition, and phenotype in non-infarct-related arteries of STEMI patients undergoing primary percutaneous coronary intervention (PCI). METHODS AND RESULTS Between September 2009 and January 2011, 103 STEMI patients underwent intravascular ultrasonography (IVUS) and radiofrequency ultrasonography (RF-IVUS) of the two non-infarct-related epicardial coronary arteries (non-IRA) after successful primary PCI. Patients were treated with high-intensity rosuvastatin (40 mg/day) throughout 13 months and serial intracoronary imaging with the analysis of matched segments was available for 82 patients with 146 non-IRA. The primary IVUS end-point was the change in per cent atheroma volume (PAV). After 13 months, low-density lipoprotein cholesterol (LDL-C) had decreased from a median of 3.29 to 1.89 mmol/L (P < 0.001), and high-density lipoprotein cholesterol (HDL-C) levels had increased from 1.10 to 1.20 mmol/L (P < 0.001). PAV of the non-IRA decreased by -0.9% (95% CI: -1.56 to -0.25, P = 0.007). Patients with regression in at least one non-IRA were more common (74%) than those without (26%). Per cent necrotic core remained unchanged (-0.05%, 95% CI: -1.05 to 0.96%, P = 0.93) as did the number of RF-IVUS defined thin cap fibroatheromas (124 vs. 116, P = 0.15). CONCLUSION High-intensity rosuvastatin therapy over 13 months is associated with regression of coronary atherosclerosis in non-infarct-related arteries without changes in RF-IVUS defined necrotic core or plaque phenotype among STEMI patients.

Peak Last Glacial weathering intensity on the North American continent recorded by the authigenic Hf isotope composition of North Atlantic deep-sea sediments

We have retrieved radiogenic hafnium (Hf) isotope compositions (ɛHf) from authigenic Fe–Mn oxyhydroxides of deep northwest Atlantic sediments deposited over the past 26 ka to investigate the oceanic evidence of changes in dissolved weathering inputs from NE America during the last deglaciation. The extraction of seawater-derived Hf isotopic compositions from Fe–Mn oxyhydroxides is not a standard procedure. Comparisons between the Al/Hf ratios and Hf isotopic compositions of the chemically extracted authigenic phase on the one hand, and those of the corresponding detrital fractions on the other, provide evidence that the composition of past seawater has been reliably obtained for most sampled depths with our leaching procedures. This is endorsed most strongly by data for a sediment core from 4250 m water depth at the deeper Blake Ridge, for which consistent replicates were produced throughout. The Hf isotopic composition of the most recent sample in this core also closely matches that of nearby present day central North Atlantic seawater. Comparison with previously published seawater Nd and Pb isotope compositions obtained on the same cores shows that both Hf and Pb were released incongruently during incipient chemical weathering, but responded differently to the deglacial retreat of the Laurentide Ice Sheet. Hafnium was released more congruently during peak glacial conditions of the Last Glacial Maximum (LGM) and changed to typical incongruent interglacial ɛHf signatures either during or shortly after the LGM. This indicates that some zircon-derived Hf was released to seawater during the LGM. Conversely, there is no clear evidence for an increase in the influence of weathering of Lu-rich mineral phases during deglaciation, possibly since relatively unradiogenic Hf contributions from feldspar weathering were superimposed. While the authigenic Pb isotope signal in the same marine sediment samples traced peak chemical weathering rates on continental North America during the transition to the Holocene a similar incongruent excursion is notably absent in the Hf isotope record. The early change towards more radiogenic ɛHf in relation to the LGM may provide direct evidence for the transition from a cold-based to a warm-based Laurentide Ice Sheet on the Atlantic sector of North America.