Dense active appearance models using a bounded diameter minimum spanning tree

We present a method for producing dense Active Appearance Models (AAMs), suitable for video-realistic synthesis. To this end we estimate a joint alignment of all training images using a set of pairwise registrations and ensure that these pairwise registrations are only calculated between similar images. This is achieved by defining a graph on the image set whose edge weights correspond to registration errors and computing a bounded diameter minimum spanning tree (BDMST). Dense optical flow is used to compute pairwise registration and we introduce a flow refinement method to align small scale texture. Once registration between training images has been established we propose a method to add vertices to the AAM in a way that minimises error between the observed flow fields and a flow field interpolated between the AAM mesh points. We demonstrate a significant improvement in model compactness using the proposed method and show it dealing with cases that are problematic for current state-of-the-art approaches.

Computer-aided design of PV/wind hybrid system

A complete set of match calculation methods for optimum sizing of PV/wind hybrid system is presented. In this method, the more accurate and practical mathematic models for characterizing PV module, wind generator and battery are adopted; combining with hourly measured meteorologic data and load data, the performance of a PV/wind hybrid system is determined on a hourly basis; by fixing the capacity of wind generators, the whole year's LPSP (loss of power supply probability) values of PV/wind hybrid systems with different capacity of PV array and battery bank are calculated, then the trade-off curve between battery bank and PV array capacity is drawn for the given LPSP value; the optimum configuration which can meet the energy demand with the minimum cost can be found by drawing a tangent to the trade-off curve with the slope representing the relationship between cost of PV module and that of the battery. According to this match calculation method, a set of match calculation programs for optimum sizing of PV/wind hybrid systems have been developed. Applying these match calculation programs to an assumed PV/wind hybrid system to be installed at Waglan island of Hong Kong, the optimum configuration and its hourly, daily, monthly and yearly performances are given. (C) 2003 Elsevier Science Ltd. All rights reserved.

Calculation of the hourly and daily radiation incident on three step tracking planes

Mathematical formulas for estimating the hourly and daily radiation incident on planes of azimuth three step tracking and hour angle three step tracking have been derived in this paper. Based on the hourly solar radiation data of an average day in each month at Er-Lian-Hao-Te city, the hourly and monthly radiation received by planes of these two kinds of tracking have been calculated. The results show that in this district, one axis azimuth three step tracking and hour angle three step tracking could, respectively, obtain 66.5% and 63.3% higher radiation than that on the horizontal surface all year. Moreover, a two axis azimuth three step tracking plane could receive 72% more radiation than the horizontal surface. (C) 2002 Elsevier Science Ltd. All rights reserved.

Resonant frequencies and quality factors for optical equilateral triangle resonators calculated by FDTD technique and the Pade approximation

The mode wavelength and quality factor (Q-factor) for resonant modes in optical equilateral triangle resonators (ETR's) are calculated by the finite-difference time-domain (FDTD) technique and the Pade approximation, For an ETR with the side length of 3 mu m and the refractive index of 3.2, we get the mode wavelength interval of about 70 nm and the Q-factor of the fundamental mode over 10(3), The results show that the ETR is suitable to realize single-mode operation, and that the radiation loss in the corner regions of ETR is rather low, In addition, the numerical results of the mode wavelength agree very well with our analytical formula.

CALCULATED ELECTRONIC AND MAGNETIC-STRUCTURES OF THE NEW TERNARY RARE-EARTH IRON NITRIDE ND2FE17N3

The electronic and magnetic structures of Nd2Fe17 and Nd2Fe17N3 have been calculated using the first-principle, spin-polarized orthogonalized linear combination of atomic orbitals method. Comparative studies of the two materials reveal important effects of the nitrogen atoms (at 9e site) on the electronic and magnetic structures. Results are presented for the total density of states, site-projected partial density of states and the spin magnetic moments on four nonequivalent Fe sites. The highest magnetic moments are found to be located on the 6c site for Nd2Fe17 and on the 9d site for Nd2Fe17N3, in agreement with the neutron and Mossbauer experiments. The variation trends of the magnetic moments on different Fe sites are discussed in terms of the separation between Fe and N atoms. Compared with Nd2Fe17, an increase in the exchange splitting of the Fe d band is found in Nd2Fe17N3, which accounts for its higher Curie temperature as observed in experiments. The calculated results show that the nitrogen atoms are charge acceptors in these compounds.

黄土高原(陕西部分)雨水资源化潜力初步分析

根据黄土高原（ 陕西部分） 的１４５ 个雨量观测站的３０ａ 逐月降水资料，分析了陕西雨水资源的时空分布特征，并在次基础上，分别对西北地区的作物地土壤拦蓄潜力，陕北、渭北、关中３ 地区居民工矿与交通用地的雨水汇流潜力进行了测算。初步结果为：两北地区作物地年总无效蒸发耗水达３３ ×１０８ ｍ３ ，若采取集流保墒措施，年可减少蒸发损失６－４ ×１０８ ｍ３ ；３ 地区居民工矿与交通用地年可汇流６－２７ ×１０８ ｍ ３ ，若一半用于粮食生产，并采取节水灌溉，初步推算黄土高原地区可增产粮食约２８－８ ×１０８ ｋｇ。

A field study on the conversion ratio of phytoplankton biomass carbon to chlorophyll-a in Jiaozhou Bay, China

A one-year field study was conducted to determine the conversion ratio of phytoplankton biomass carbon (Phyto-C) to chlorophyll-a (Chl-a) in Jiaozhou Bay, China. We measured suspended particulate organic carbon (POC) and phytoplankton Chl-a samples collected in surface water monthly from March 2005 to February 2006. The temporal and spatial variations of Chl-a and POC concentrations were observed in the bay. Based on the field measurements, a linear regression model II was used to generate the conversion ratio of Phyto-C to Chl-a. In most cases, a good linear correlation was found between the observed POC and Chl-a concentrations, and the calculated conversion ratios ranged from 26 to 250 with a mean value of 56 A mu g A mu g(-1). The conversion ratio in the fall was higher than that in the winter and spring months, and had the lowest values in the summer. The ratios also exhibited spatial variations, generally with low values in the near shore regions and relatively high values in offshore waters. Our study suggests that temperature was likely to be the main factor influencing the observed seasonal variations of conversion ratios while nutrient supply and light penetration played important roles in controlling the spatial variations.

Seasonal variation in metabolism of cultured Pacific oyster, Crassostrea gigas, in Sanggou Bay, China

Rates of respiration and excretion of the Pacific oyster, Crassostrea gigas, were measured seasonally from June 2002 to July 2003 under ambient conditions of food, water temperature, pH, and salinity in Sanggou Bay, an important mariculture coast in north China. The aim of this study is to obtain fundamental data for further establishing an energy budget model and assessing the carrying capacity for cultivation of C. gigas in north China. Oysters were collected monthly or bimonthly from the integrated culture areas of bivalve and kelp in the bay. Oxygen consumption and ammonium and phosphorus excretion rates were measured, and ratios of O/N and NIP were calculated. One-way ANOVA was applied to determine differences among these parameters that act as a function of seasonal variation. All the physiological parameters yielded highly significant variations with season (P<0.01) The rate of respiration varied seasonally, with the highest oxygen consumption rate in July and the lowest rate in January, ranging from 0.07 to 2.13 mg O-2 h(-1) g(-1) dry tissue weight (DW). Maximum and minimum ammonium excretion rates were recorded in August and January, respectively, ranging from 0.51 to 5.40 mu mol NH4-N h(-1) g(-1) DW. Rates of phosphorus excretion varied from 0.11 (in January) to 0.64 (in July) mu mol PO4-P h(-1) g(-1) DW. The O/N and N/P ratios changed from 9.2 (in January) to 59.8 (in July) and from 4.6 (in January) to 10.9 (in August), respectively. For each season, the allometric relationship between the physiological response (e.g., rate of oxygen consumption, ammonium and phosphorus excretion) and DW of the animal was estimated using the formula: Y=a x DWb. (C) 2005 Elsevier B.V. All rights reserved.

Seasonal changes in weight and body composition of yak grazing on alpine-meadow grassland in the Qinghai-Tibetan plateau of China

Forty-five male yaks (born April 2001) were studied to determine how seasonal changes on the Qinghai-Tibetan plateau affected BW and body composition. Thirty yaks were weighed monthly from birth to 26 mo of age to determine seasonal changes in BW. The remaining 15 yaks were allocated randomly to five groups (three yaks per group), designated for slaughter at 13, 15, 18, 22, and 25 mo to measure seasonal effects on body chemical composition. All yaks were grazed on the alpine-meadow grassland of the plateau without any supplementation. All BW and body composition data were calculated on an individual basis. Body weight and body composition data were both compared across seven growth periods spanning 2 yr and defined by season. From April (birth) to December 2001 of the first growing season, yak BW increased (P < 0.01); however, during the subsequent cold season (December 2001 to May 2002), BW decreased (P < 0.01). The second growing season ran from May 2002 (13 mo of age) to October 2002 (18 mo of age), and the second live weight-loss season ran from October 2002 until May 2003. The weight loss experienced by yaks during the first weight loss season was 25.64% of the total weight gain in the first growing season. The weight loss experienced by yaks during the second weight-loss season was 29.73% of the total weight gain in the second growing season. Energy retention in the second growing season was 291.07 MJ, 50.8% of which was consumed during the subsequent cold season. Energy accumulation in the summer (from May to July) and fall (from July to October) of the second growing season did not differ (5.01 and 6.30 MJ/kg of EBW gain, respectively; P = 0.63). The energy mobilized during the second winter (from October 2002 to February 2003) was 16.49 MJ/kg of EBW, and in the second spring (from February to May 2003), it was 9.06 MJ/kg of EBW. These data suggest that the decrease in grazing yak BW during the first cold season is much less than during the second cold season, and that the energy content per unit of BW mobilized is greater (P = 0.02) in winter than in spring. Results from this study demonstrate highly efficient compensatory growth in grazing yaks following the first weight loss period during the first cold season. This benefit could be exploited by herders to improve yak production. Yaks may have developed a type of self-protection mechanism to overcome the long cold seasons in the Qinghai-Tibetan plateau.

投资时间序列的混沌分维研究

以某市固定资产投资同比增长数据为根据 ,利用混沌动力学中的相空间重构技术、幅值谱分析、庞加莱截面等方法进行分析 ,并计算其关联分维 ,得出其饱和嵌入维数为 5、关联分维约为 1 7,初步确定某市固定资产投资增长过程具有混沌特性 ,有 5个因素影响着某市固定资产投资的同比增长幅度 ,其中 2个因素起决定作用·该法可用于复杂的经济系统的分析研究·

电离层参量M(3000)F2和hmF2的经验正交函数方法建模研究

The ionospheric parameter M(3000)F2 (the so-called transmission factor or the propagation factor) is important not only in practical applications such as frequency planning for radio-communication but also in ionospheric modeling. This parameter is strongly anti-correlated with the ionospheric F2-layer peak height hmF2，a parameter often used as a key anchor point in some widely used empirical models of the ionospheric electron density profile (e.g., in IRI and NeQuick models). Since hmF2 is not easy to obtain from measurements and M(3000)F2 can be routinely scaled from ionograms recorded by ionosonde/digisonde stations distributed globally and its data has been accumulated for a long history, usually the value of hmF2 is calculated from M(3000)F2 using the empirical formula connecting them. In practice, CCIR M(3000)F2 model is widely used to obtain M(3000)F2 value. However, recently some authors found that the CCIR M(3000)F2 model has remarkable discrepancies with the measured M(3000)F2, especially in low-latitude and equatorial regions. For this reason, the International Reference Ionosphere (IRI) research community proposes to improve or update the currently used CCIR M(3000)F2 model. Any efforts toward the improvement and updating of the current M(3000)F2 model or newly development of a global hmF2 model are encouraged. In this dissertation, an effort is made to construct the empirical models of M(3000)F2 and hmF2 based on the empirical orthogonal function (EOF) analysis combined with regression analysis method. The main results are as follows: 1. A single station model is constructed using monthly median hourly values of M(3000)F2 data observed at Wuhan Ionospheric Observatory during the years of 1957–1991 and compared with the IRI model. The result shows that EOF method is possible to use only a few orders of EOF components to represent most of the variance of the original data set. It is a powerful method for ionospheric modeling. 2. Using the values of M(3000)F2 observed by ionosondes distributed globally, data at grids uniformly distributed globally were obtained by using the Kriging interpolation method. Then the gridded data were decomposed into EOF components using two different coordinates: (1) geographical longitude and latitude; (2) modified dip (Modip) and local time. Based on the EOF decompositions of the gridded data under these two coordinates systems, two types of the global M(3000)F2 model are constructed. Statistical analysis showed that the two types of the constructed M(3000)F2 model have better agreement with the observational M(3000)F2 than the M(3000)F2 model currently used by IRI. The constructed models can represent the global variations of M(3000)F2 better. 3. The hmF2 data used to construct the hmF2 model were converted from the observed M(3000)F2 based on the empirical formula connecting them. We also constructed two types of the global hmF2 model using the similar method of modeling M(3000)F2. Statistical analysis showed that the prediction of our models is more accurate than the model of IRI. This demonstrated that using EOF analysis method to construct global model of hmF2 directly is feasible. The results in this thesis indicate that the modeling technique based on EOF expansion combined with regression analysis is very promising when used to construct the global models of M(3000)F2 and hmF2. It is worthwhile to investigate further and has the potential to be used to the global modeling of other ionospheric parameters.

双聚焦偏移速度估计与可分表示叠前弹性反演

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.

Development of process control methodology for tracking the quality and safety of pain, agitation and sedation management in critical care units

Objective: To develop sedation, pain, and agitation quality measures using process control methodology and evaluate their properties in clinical practice. Design: A Sedation Quality Assessment Tool was developed and validated to capture data for 12-hour periods of nursing care. Domains included pain/discomfort and sedation-agitation behaviors; sedative, analgesic, and neuromuscular blocking drug administration; ventilation status; and conditions potentially justifying deep sedation. Predefined sedation-related adverse events were recorded daily. Using an iterative process, algorithms were developed to describe the proportion of care periods with poor limb relaxation, poor ventilator synchronization, unnecessary deep sedation, agitation, and an overall optimum sedation metric. Proportion charts described processes over time (2 monthly intervals) for each ICU. The numbers of patients treated between sedation-related adverse events were described with G charts. Automated algorithms generated charts for 12 months of sequential data. Mean values for each process were calculated, and variation within and between ICUs explored qualitatively. Setting: Eight Scottish ICUs over a 12-month period. Patients: Mechanically ventilated patients. Interventions: None. Measurements and Main Results: The Sedation Quality Assessment Tool agitation-sedation domains correlated with the Richmond Sedation Agitation Scale score (Spearman [rho] = 0.75) and were reliable in clinician-clinician (weighted kappa; [kappa] = 0.66) and clinician-researcher ([kappa] = 0.82) comparisons. The limb movement domain had fair correlation with Behavioral Pain Scale ([rho] = 0.24) and was reliable in clinician-clinician ([kappa] = 0.58) and clinician-researcher ([kappa] = 0.45) comparisons. Ventilator synchronization correlated with Behavioral Pain Scale ([rho] = 0.54), and reliability in clinician-clinician ([kappa] = 0.29) and clinician-researcher ([kappa] = 0.42) comparisons was fair-moderate. Eight hundred twenty-five patients were enrolled (range, 59-235 across ICUs), providing 12,385 care periods for evaluation (range 655-3,481 across ICUs). The mean proportion of care periods with each quality metric varied between ICUs: excessive sedation 12-38%; agitation 4-17%; poor relaxation 13-21%; poor ventilator synchronization 8-17%; and overall optimum sedation 45-70%. Mean adverse event intervals ranged from 1.5 to 10.3 patients treated. The quality measures appeared relatively stable during the observation period. Conclusions: Process control methodology can be used to simultaneously monitor multiple aspects of pain-sedation-agitation management within ICUs. Variation within and between ICUs could be used as triggers to explore practice variation, improve quality, and monitor this over time

Patient Retention in Antiretroviral Therapy Programs in Sub-Saharan Africa: A Systematic Review

Paper published in PLoS Medicine in 2007.

Long-term dynamics of death rates of emphysema, asthma, and pneumonia and improving air quality.

BACKGROUND: The respiratory tract is a major target of exposure to air pollutants, and respiratory diseases are associated with both short- and long-term exposures. We hypothesized that improved air quality in North Carolina was associated with reduced rates of death from respiratory diseases in local populations. MATERIALS AND METHODS: We analyzed the trends of emphysema, asthma, and pneumonia mortality and changes of the levels of ozone, sulfur dioxide (SO2), nitrogen dioxide (NO2), carbon monoxide (CO), and particulate matters (PM2.5 and PM10) using monthly data measurements from air-monitoring stations in North Carolina in 1993-2010. The log-linear model was used to evaluate associations between air-pollutant levels and age-adjusted death rates (per 100,000 of population) calculated for 5-year age-groups and for standard 2000 North Carolina population. The studied associations were adjusted by age group-specific smoking prevalence and seasonal fluctuations of disease-specific respiratory deaths. RESULTS: Decline in emphysema deaths was associated with decreasing levels of SO2 and CO in the air, decline in asthma deaths-with lower SO2, CO, and PM10 levels, and decline in pneumonia deaths-with lower levels of SO2. Sensitivity analyses were performed to study potential effects of the change from International Classification of Diseases (ICD)-9 to ICD-10 codes, the effects of air pollutants on mortality during summer and winter, the impact of approach when only the underlying causes of deaths were used, and when mortality and air-quality data were analyzed on the county level. In each case, the results of sensitivity analyses demonstrated stability. The importance of analysis of pneumonia as an underlying cause of death was also highlighted. CONCLUSION: Significant associations were observed between decreasing death rates of emphysema, asthma, and pneumonia and decreases in levels of ambient air pollutants in North Carolina.