An optimized leaf-setting algorithm for beam intensity modulation using dynamic multileaf collimators

A leaf-setting algorithm is developed for generating arbitrary beam intensity profiles in discrete levels using dynamic multileaf collimators (DMLCs). The algorithm starts with the algebraic expression for the area under the beam profile. It is shown that the coefficients in this expression can be transformed into the specifications for the leaf-setting sequence. It is proven that the algorithm optimizes beam delivery time and total monitor units for the DMLC leaf setting for intensity modulated radiotherapy (IMRT). The algorithm is demonstrated to be applicable to both the 'step-and-shoot' and 'dynamic' type of beam delivery. The graphical interpretation and numerical implementation scheme of the algorithm is illustrated using a simplified example.     

Efficient scatter modelling for incorporation in maximum likelihood reconstruction

Definition of a simplified model of scatter which can be incorporated in maximum likelihood reconstruction for single-photon emission tomography (SPET) continues to be appealing; however, implementation must be efficient for it to be clinically applicable. In this paper an efficient algorithm for scatter estimation is described in which the spatial scatter distribution is implemented as a spatially invariant convolution for points of constant depth in tissue. The scatter estimate is weighted by a space-dependent build-up factor based on the measured attenuation in tissue. Monte Carlo simulation of a realistic thorax phantom was used to validate this approach. Further efficiency was introduced by estimating scatter once after a small number of iterations using the ordered subsets expectation maximisation (OSEM) reconstruction algorithm. The scatter estimate was incorporated as a constant term in subsequent iterations rather than modifying the scatter estimate each iteration. Monte Carlo simulation was used to demonstrate that the scatter estimate does not change significantly provided at least two iterations OSEM reconstruction, subset size 8, is used. Complete scatter-corrected reconstruction of 64 projections of 40?128 pixels was achieved in 38 min using a Sun Sparc20 computer.     

A method for implementing dynamic photon beam intensity modulation using independent jaws and a multileaf collimator

A mathematical model is derived for digitally controlled linear accelerators to deliver a desired photon intensity distribution by combining collimator motion and machine dose rate variations. It shows that, at any instant, the quotient of the machine dose rate and the speed of collimator motion is proportional to the gradient of the desired in-air photon fluence distribution. The model is applicable for both independently controlled collimator jaws and multileaf collimators and can be implemented by controlling different parameters to accommodate linear accelerators from different manufactures. For independent jaws, each pair of jaws creates photon fluence variations along the direction of the jaw movement. For multileaf collimators, where each leaf is independently controlled, any two-dimensional (2D) photon fluence distribution can be delivered. The model has been implemented for wedged isodose distributions using independent jaws, and 2D intensity modulation using a multileaf collimator. One-dimensional (1D) wedged isodose distributions are created by moving an independent jaw at constant speed while varying machine dose rate. 2D intensity modulation has been implemented using a 'dynamic stepping' scheme, which controls the leaf progression during irradiation at constant machine dose rate. With this automated delivery scheme, the beam delivery time for dynamic intensity modulation, which depends on the complexity of the desired intensity distribution, approaches that of conventional beam modifiers. This paper shows the derivation of the model, its application, and our delivery scheme. Examples of 1D dynamic wedges and 2D intensity modulations will be given to illustrate the versatility of the model, the simplicity of its application, and the efficiency of beam delivery. These features make this approach practical for delivering conformal therapy treatments.     

用静电场高斯定理求电场强度和磁场环路定律求磁场强度的讨论

本文要讨论两个问题：1运用电场高斯定律求场强的条件和高斯面选取的方法;2．运用磁场环路定律求磁场H的条件和环路选取的方法,并通过举例说明它们的具体分析办法。     

An equivalent square field formula for determining head scatter factors of rectangular fields

A simple formula is derived for the calculation of an equivalent square field that gives the same head scatter factor as a given rectangular field. This formula is based strictly on the configuration of a medical linear accelerator treatment head. The geometric parameters used are the distances between the target and the top of each field-defining aperture. The formula accounts for both the effect of field elongation and the collimator exchange effect. This method predicts the output to within 1% accuracy for both open and wedged fields and does not require any new measured data other than the field size dependence of head scatter for a range of square field sizes. Interestingly, the formula we derived has the same format as the formula that was empirically obtained by Vadash and Bj?rngard [Med. Phys. 20, 733-734 (1993)].     

Measurements of phantom scatter factors for small field sizes in high energy x rays

An alternate method to measure the phantom scatter factor in small fields is provided for high energy photon beams. The measurement technique is based on the density scaling theorem described by O'Connor [Phys. Med. Biol. 1, 352-369 (1957)]. The phantom scatter factor (Sp) is measured in balsa and cedar wood to give effective field sizes in 3 x 3 to 0.5 x 0.5 cm2 water. The extrapolated zero area phantom scatter factor from the average data of balsa and cedar is 0.45. This indicates that the variation of output is largely due to a sharp decrease in the phantom scatter factor for small field sizes.     

Electrodermal activity and schizophrenia: The problem of stimulus intensity modulation.

Recent reviews on electrodermal activity and schizophrenia have cited both hyperreactivity and hyporeactivity in groups of schizophrenics. One problem concerning the interpretation of hyporeactivity is that of stimulus-intensity modulation. This problem is defined, and it is suggested that researchers interested in basic neurophysiological activity of schizophrenics employ a nonstressful experimental procedure. (PsycINFO Database Record (c) 2010 APA, all rights reserved)     

带电圆盘导体电场的电势和场强

导出了带电薄圆盘导体空间电场的电势和场强表达式，并作出了相应的分布曲线．     

通道式感应加热五流中间包流场的水力学模拟

 通道式感应加热能有效补偿连铸过程中间包钢水的温降,且具加热效率高、设备简单等优点,是近年来得到快速推广应用的中间包新技术。配置感应加热器后的中间包结构与常规中间包差异较大,其加热通道的设计对中间包流场具有明显影响。针对实际应用的某5流直通道感应加热中间包边部流和中间流温差大、钢水浇铸过热度偏高的问题,通过水模拟试验对该中间包流场进行优化,提出了一种新型的分口通道方案。该方案可将中间包的死区比例由原型的29.50%降低到20.33%,钢水平均停留时间较原型延长了40 s。工业试验表明,分口式新型通道设计中间包各流浇铸温度一致性得到改善,其中边部流和中间流的平均温差较原型降低了3.6 ℃,实现了感应加热中间包应用效果的进一步提升。     

Analgesic effect of vaginal stimulation in rats: modulation by graded stimulus intensity and hormones

Regulation of nitrogen fixation by ammonium and glutamate was examined in Rhizobium sp. 32H1 growing in defined liquid media. Whereas nitrogenase synthesis in Klebsiella pneunoniae is normally completely repressed during growth on NH4+, nitrogenase activity was detected in cultures of Rhizobium sp. grown with excess NH4+. However, an "ammonium effect" on activity was invariably observed in cultures grown on NH4+ as sole nitrogen source; the nitrogenase activity was, depending on conditions, 14 to 36% of that of comparable glutamate-grown cultures. Glutamate inhibited utilization of exogenous NH4+ and, in one of two procedures described, glutamate partially alleviated the ammonium effect on nitrogenase activity. NH4+, apparently produced from N2, was excreted into the culture medium when growth was initiated on glutamate, but not when NH4+ was thesole source of fixed nitrogen for growth. These findings are discussed in relation to nitrogen fixation by Rhizobium bacteroids.     

Perfusion and thermal field during hyperthermia. Experimental measurements and modelling in recurrent breast cancer

Recurrences of malignant tumours in the chest wall are proposed as a valuable model of tissue mainly perfused by small size vessels (the so-called 'phase III' vessels). Invasive thermal measurements have been performed on two patients affected by cutaneous metastasis of malignant tumours during hyperthermic sessions. Thermal probes were inserted into catheters implanted into the tissue at different depths. In one of the catheters a probe connected with laser-Doppler equipment was inserted to assess blood perfusion in the tumour periphery. The perfusion was monitored throughout the sessions, and a noticeable temporal variability was observed. The effect of the perfusion on the thermal map in the tissue was evaluated locally and the 'effective conductivity' of the perfused tissue was estimated by means of the numerical integration of the 'bio-heat' equation. The tumour temperature, at the site where the perfusion probe is located, can be predicted by the numerical model provided two free parameters, alpha and beta, are evaluated with a fitting procedure. Alpha is related to the effective conductivity and beta to the SAR term of the bio-heat equation. The model aimed at estimating the 'effective conductivity' K(eff) of the perfused tissue, and average values of K(eff) of 0.27 +/- 0.03 W m(-1) degrees C(-1) in Patient 1 and of 0.665 +/- 0.005 W m(-1) degrees C(-1) in Patient 2 were obtained throughout the treatment. However, when the average temperature in a larger tumour volume is to be predicted but only a single, 'local' measurement of the perfusion is available and is assumed to be representative for the whole region, the model results are far less satisfactory. This is probably due to the fact that changes of blood perfusion throughout hyperthermic sessions occur to different extents within the tumour volume, and the differences in perfusion cannot be ignored. The above result suggests that, in addition to the 'temperature map', also a 'perfusion map' within the heated volume should be monitored routinely throughout hyperthermic sessions.     

Relations between scatter factor, quality index and attenuation for x-ray beams

In a human in-vitro fertilization (IVF) programme, the effect of co-culture of embryos with human fibroblasts was evaluated with respect to pregnancy rate and embryo development. Patients were included in the study after giving informed written consent. The IVF treatments were randomly assigned by stratification of both age (<36 versus > or =36 years) and previous IVF attempts (yes versus no). After fertilization was established, the zygotes were transferred to a 4-well dish with or without fibroblasts and cultured for 2 days. On the third day after ovum pick-up (OPU), cell number and quality [5 (good) to 1 (poor)] of the embryos were scored and a maximum of three embryos was transferred. Supernumerary embryos of good quality were cryopreserved. The design of this study was a group sequential trial with the objective of detecting differences between pregnancy rates following IVF with conventional incubation or incubation in co-culture with fibroblasts. This design included one evaluation at half-way data collection. In the study, 148 patients had an OPU, of whom 77 were allocated to the co-culture group. There was no statistically significant difference in pregnancy rate, cell number and embryo quality between the two groups. The ongoing pregnancy rate per embryo transfer was 27% in co-culture and 30% in the conventional culture group. The implantation rates per transferred embryo were 17 and 18% respectively. Using a multivariate logistic regression model for the probability of ongoing pregnancies, the odds ratio of co-culture, adjusted for age and previous IVF attempts, was not statistically significant. In conclusion, co-culture with human fibroblasts does not contribute to an improvement of embryo quality nor to a higher pregnancy rate after IVF in an unselected group of patients.     

Secondary shielding barriers for diagnostic x-ray facilities: scatter and leakage revisited

Methods are presented for determining the thickness requirements for barriers against scatter and leakage radiation generated in a diagnostic x-ray facility. Equations are developed that express the shielded doses due to scatter and leakage with explicit dependence on the operating potential used in clinical settings. The 1972 scatter experiment of Trout and Kelley is revisited, with suggested values for the scatter fraction that are somewhat different from those used in Report Number 49 of the NCRP. The dose from leakage radiation due to operation at clinical kVp settings is seen to be orders of magnitude below that predicted by the model in NCRP49. The unshielded secondary radiation dose in air at unit distance and workload is presented as a function of operating potential and for clinical workload distributions. The net transmission of secondary radiations through common shielding materials is presented. Compared to that for the primary beam at the same potential, the transmission of secondary radiation is more penetrating due to the inclusion of leakage radiation that has been hardened by the tube housing. Sample shielding calculations illustrate the utility of this scheme.     

An adaptive control algorithm for optimization of intensity modulated radiotherapy considering uncertainties in beam profiles, patient set-up and internal organ motion

A new general beam optimization algorithm for inverse treatment planning is presented. It utilizes a new formulation of the probability to achieve complication-free tumour control. The new formulation explicitly describes the dependence of the treatment outcome on the incident fluence distribution, the patient geometry, the radiobiological properties of the patient and the fractionation schedule. In order to account for both measured and non-measured positioning uncertainties, the algorithm is based on a combination of dynamic and stochastic optimization techniques. Because of the difficulty in measuring all aspects of the intra- and interfractional variations in the patient geometry, such as internal organ displacements and deformations, these uncertainties are primarily accounted for in the treatment planning process by intensity modulation using stochastic optimization. The information about the deviations from the nominal fluence profiles and the nominal position of the patient relative to the beam that is obtained by portal imaging during treatment delivery, is used in a feedback loop to automatically adjust the profiles and the location of the patient for all subsequent treatments. Based on the treatment delivered in previous fractions, the algorithm furnishes optimal corrections for the remaining dose delivery both with regard to the fluence profile and its position relative to the patient. By dynamically refining the beam configuration from fraction to fraction, the algorithm generates an optimal sequence of treatments that very effectively reduces the influence of systematic and random set-up uncertainties to minimize and almost eliminate their overall effect on the treatment. Computer simulations have shown that the present algorithm leads to a significant increase in the probability of uncomplicated tumour control compared with the simple classical approach of adding fixed set-up margins to the internal target volume.