Cosine Of Angle Research Articles

Evaluation of Radiometric and Atmospheric Correction Algorithms for Aboveground Forest Biomass Estimation Using Landsat 5 TM Data

Solar radiation is affected by absorption and emission phenomena during its downward trajectory from the Sun to the Earth’s surface and during the upward trajectory detected by satellite sensors. This leads to distortion of the ground radiometric properties (reflectance) recorded by satellite images, used in this study to estimate aboveground forest biomass (AGB). Atmospherically-corrected remote sensing data can be used to estimate AGB on a global scale and with moderate effort. The objective of this study was to evaluate four atmospheric correction algorithms (for surface reflectance), ATCOR2 (Atmospheric Correction for Flat Terrain), COST (Cosine of the Sun Zenith Angle), FLAASH (Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes) and 6S (Second Simulation of Satellite Signal in the Solar), and one radiometric correction algorithm (for reflectance at the sensor) ToA (Apparent Reflectance at the Top of Atmosphere) to estimate AGB in temperate forest in the northeast of the state of Durango, Mexico. The AGB was estimated from Landsat 5 TM imagery and ancillary information from a digital elevation model (DEM) using the non-parametric multivariate adaptive regression splines (MARS) technique. Field reference data for the model training were collected by systematic sampling of 99 permanent forest growth and soil research sites (SPIFyS) established during the winter of 2011. The following predictor variables were identified in the MARS model: Band 7, Band 5, slope (β), Wetness Index (WI), NDVI and MSAVI2. After cross-validation, 6S was found to be the optimal model for estimating AGB (R2 = 0.71 and RMSE = 33.5 Mg·ha−1; 37.61% of the average stand biomass). We conclude that atmospheric and radiometric correction of satellite images can be used along with non-parametric techniques to estimate AGB with acceptable accuracy.

Pressure-Transient Analysis of a Finite-Conductivity Inclined Fracture Connected to a Slanted Wellbore

Summary The principal focus of this work is on pressure-transient analysis of a finite-conductivity inclined fracture connected to a slanted wellbore, on the basis of a semianalytical model. Detailed analysis of unsteady-state pressure behavior of a fully penetrating inclined fracture in an infinite-slab reservoir was provided. The study has shown that a finite-conductivity inclined fracture may exhibit five flow regimes: bilinear flow, formation linear flow, early radial flow, compound linear flow, and pseudoradial flow. The characteristics of bilinear flow and formation linear flow are predominantly determined by fracture conductivity. In the case of a low formation-thickness/fracture-half-length ratio and small inclination angle, both early radial flow and compound linear flow may be absent. Analytical solutions for transient responses during different flow regimes are similar to that for a fully penetrating vertical fracture and can be correlated with the cosine of fracture-inclination angle with consideration of permeability anisotropy. Effect of inclination angle and reservoir-permeability anisotropy on transient responses is strong, which extends to pseudoradial-flow period. Formation thickness mainly influences the middle to late flow periods. In addition, the pseudoskin factor is also investigated in detail.

A new approach to the real-time assessment of the clear-sky direct normal irradiance

In a context of sustainable development, interest for concentrating solar power and concentrating photovoltaic technologies is growing rapidly. One of the most challenging topics is to improve solar resource assessment and forecasting in order to optimize power plant operation. Since clear sky defines the nominal operating conditions of the plants, improving their management requires the use in real-time of clear-sky direct normal irradiance (DNI) models. Typically, accuracy is best achieved by considering water vapor and aerosol concentrations in the atmosphere separately. However, measuring such physical quantities is not easy and requires a weather station close to the considered site. When these data are not available, the attenuating effects can be modeled by atmospheric turbidity factors which can be obtained from DNI under clear-sky conditions. So, the main purpose of the present paper is to propose an efficient approach to assess the clear-sky DNI in real time. This approach combines an existing empirical model, proposed by Ineichen and Perez, with a new methodology for the computation of atmospheric turbidity. It takes advantage of the fact that changes in atmospheric turbidity are relatively small throughout the day in comparison to changes in DNI, even when the sky is free of clouds. In the present study, we considered data from two experimental sites (Golden, in the USA, and Perpignan, in France) and used a wavelet-based multi-resolution analysis as a clear-sky DNI detection tool. In addition, we compared the proposed approach with several combinations of empirical models and ways of computing atmospheric turbidity. The first model is a polynomial of the cosine of the solar zenith angle, whereas the two other models use atmospheric turbidity as an additional input. Regarding its calculation, monthly and daily mean values have been considered. Moreover, we defined a procedure in order to evaluate the accuracy of all the considered approaches. This procedure allows changes in DNI caused by clouds to be simulated using a noisy signal applied to clear-sky periods. In both sites, our approach to the real-time assessment of the clear-sky DNI outperforms the other approaches. In the worst case, the mean absolute error is reduced by 8 W m−2 in comparison to the approaches based on monthly mean values of atmospheric turbidity, and reduced by about 30 W m−2 when taking the polynomial-based model as a reference.

Phytochemical screening and chemical variability in volatile oils of aerial parts of Morinda morindoides

Morinda morindoides is an important Liberian traditional medicine for the treatment of malaria, fever, worms etc. The plant was subjected to integrated approaches including phytochemical screening and gas chromatography mass spectrometry (GC–MS) analyses. Phytochemical investigation of the powdered plant revealed the presence of phenolics, tannins, flavonoids, saponins, terpenes, steroidal compounds and volatile oil. Steam distillation followed by GC–MS resulted in the identification of 47 volatiles in its aerial parts: 28 were in common including various bioactive volatiles. Major constituents of leaves were phytol (43.63%), palmitic acid (8.55%) and geranyl linalool (6.95%) and stem were palmitic acid (14.95%), eicosane (9.67%) and phytol (9.31%), and hence, a significant difference in the percentage composition of aerial parts was observed. To study seasonal changes, similarity analysis was carried out by calculating correlation coefficient (r) and vector angle cosine (z) that were more than 0.91 for stem-to-stem and leaf-to-leaf batches indicating considerable consistency.

General parametrization of axisymmetric black holes in metric theories of gravity

Following previous work of ours in spherical symmetry, we here propose a new parametric framework to describe the spacetime of axisymmetric black holes in generic metric theories of gravity. In this case, the metric components are functions of both the radial and the polar angular coordinates, forcing a double expansion to obtain a generic axisymmetric metric expression. In particular, we use a continued-fraction expansion in terms of a compactified radial coordinate to express the radial dependence, while we exploit a Taylor expansion in terms of the cosine of the polar angle for the polar dependence. These choices lead to a superior convergence in the radial direction and to an exact limit on the equatorial plane. As a validation of our approach, we build parametrized representations of Kerr, rotating dilaton, and Einstein-dilaton-Gauss-Bonnet black holes. The match is already very good at lowest order in the expansion and improves as new orders are added. We expect a similar behavior for any stationary and axisymmetric black-hole metric.

Evaluating the utility of mid-infrared spectral subspaces for predicting soil properties

We propose four methods for finding local subspaces in large spectral libraries. The proposed four methods include (a) cosine angle spectral matching; (b) hit quality index spectral matching; (c) self-organizing maps and (d) archetypal analysis methods. Then evaluate prediction accuracies for global and subspaces calibration models. These methods were tested on a mid-infrared spectral library containing 1907 soil samples collected from 19 different countries under the Africa Soil Information Service project. Calibration models for pH, Mehlich-3 Ca, Mehlich-3 Al, total carbon and clay soil properties were developed for the whole library and for the subspace. Root mean square error of prediction was used to evaluate predictive performance of subspace and global models. The root mean square error of prediction was computed using a one-third-holdout validation set. Effect of pretreating spectra with different methods was tested for 1st and 2nd derivative Savitzky–Golay algorithm, multiplicative scatter correction, standard normal variate and standard normal variate followed by detrending methods. In summary, the results show that global models outperformed the subspace models. We, therefore, conclude that global models are more accurate than the local models except in few cases. For instance, sand and clay root mean square error values from local models from archetypal analysis method were 50% poorer than the global models except for subspace models obtained using multiplicative scatter corrected spectra with which were 12% better. However, the subspace approach provides novel methods for discovering data pattern that may exist in large spectral libraries.

TOPOGRAPHIC EFFECT ON SPECTRAL VEGETATION INDICES FROM LANDSAT TM DATA: IS TOPOGRAPHIC CORRECTION NECESSARY?

The full potentiality of spectral vegetation indices (VIs) can only be evaluated after removing topographic, atmospheric and soil background effects from radiometric data. Concerning the former effect, the topographic effect was barely investigated in the context of VIs, despite the current availability correction methods and Digital elevation Model (DEM). In this study, we performed topographic correction on Landsat 5 TM spectral bands and evaluated the topographic effect on four VIs: NDVI, RVI, EVI and SAVI. The evaluation was based on analyses of mean and standard deviation of VIs and TM band 4 (near-infrared), and on linear regression analyses between these variables and the cosine of the solar incidence angle on terrain surface (cos i). The results indicated that VIs are less sensitive to topographic effect than the uncorrected spectral band. Among VIs, NDVI and RVI were less sensitive to topographic effect than EVI and SAVI. All VIs showed to be fully independent of topographic effect only after correction. It can be concluded that the topographic correction is required for a consistent reduction of the topographic effect on the VIs from rugged terrain.

Production of heavy Higgs bosons and decay into top quarks at the LHC

We investigate the production of heavy, neutral Higgs boson resonances and their decays to top-quark top-antiquark ($t\bar{t}$) pairs at the Large Hadron Collider (LHC) at next-to-leading order (NLO) in the strong coupling of quantum chromodynamics (QCD). The NLO corrections to heavy Higgs boson production and the Higgs-QCD interference are calculated in the large $m_t$ limit with an effective K-factor rescaling. The nonresonant $t\bar{t}$ background is taken into account at NLO QCD including weak-interaction corrections. In order to consistently determine the total decay widths of the heavy Higgs bosons, we consider for definiteness the type-II two-Higgs-doublet extension of the standard model and choose three parameter scenarios that entail two heavy neutral Higgs bosons with masses above the $t\bar{t}$ threshold and unsuppressed Yukawa couplings to top quarks. For these three scenarios we compute, for the LHC operating at 13 TeV, the $t\bar{t}$ cross section and the distributions of the $t\bar{t}$ invariant mass, of the transverse top-quark momentum and rapidity, and of the cosine of the Collins-Soper angle with and without the two heavy Higgs resonances. For selected $M_{t\bar{t}}$ bins we estimate the significances for detecting a heavy Higgs signal in the $t\bar{t}$ dileptonic and lepton plus jets decay channels.

Technological distance measures: new perspectives on nearby and far away

Understanding the competitive environment of one's company is crucial for every manager. One tool to quantify the technological relationships between companies, evaluate industry landscapes and knowledge transfer potential in collaborations is the technological distance. There are different methods and many different factors that impact the results and thus the conclusions that are drawn from distance calculation. Therefore, the present study derives guidelines for calculating and evaluating technological distances for three common methods, i.e. the Euclidean distance, the cosine angle and the min-complement distance. For this purpose, we identify factors that influence the results of technological distance calculation using simulation. Subsequently, we analyze technological distances of cross-industry collaborations in the field of electric mobility. Our findings show that a high level of detail is necessary to achieve insightful results. If the topic in scope of the analysis does not represent the core business of the companies, we recommend filters to focus on the respective topic. Another key suggestion is to compare the calculated results to a peer group in order to evaluate if a distance can be evaluated as `near' or `far'.

Investigation of terrain illumination effects on vegetation indices and VI-derived phenological metrics in subtropical deciduous forests

We used RapidEye and Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra data to study terrain illumination effects on 3 vegetation indices (VIs) and 11 phenological metrics over seasonal deciduous forests in southern Brazil. We applied TIMESAT for the analysis of the Enhanced Vegetation Index (EVI) and the Normalized Difference Vegetation Index (NDVI) derived from the MOD13Q1 product to calculate phenological metrics. We related the VIs with the cosine of the incidence angle i (Cos i) and inspected percentage changes in VIs before and after topographic C-correction. The results showed that the EVI was more sensitive to seasonal changes in canopy biophysical attributes than the NDVI and Red-Edge NDVI, as indicated by analysis of non-topographically corrected RapidEye images from the summer and winter. On the other hand, the EVI was more sensitive to terrain illumination, presenting higher correlation coefficients with Cos i that decreased with reduction in the canopy background L factor. After C-correction, the RapidEye Red-Edge NDVI, NDVI, and EVI decreased 2%, 1%, and 13% over sunlit surfaces and increased up to 5%, 14%, and 89% over shaded surfaces, respectively. The EVI-related phenological metrics were also much more affected by topographic effects than the NDVI-derived metrics. From the set of 11 metrics, the 2 that described the period of lower photosynthetic activity and seasonal VI amplitude presented the largest correlation coefficients with Cos i. The results showed that terrain illumination is a factor of spectral variability in the seasonal analysis of phenological metrics, especially for VIs that are not spectrally normalized.

Quantitative analysis of highly similar salvianolic acids with 1H qNMR for quality control of traditional Chinese medicinal preparation Salvianolate Lyophilized Injection

Salvianolate Lyophilized Injection (SLI), a traditional Chinese medicinal (TCM) preparation which is used to treat stroke, is composed of multiple salvianolic acids from the aqueous extracts of Salvia miltiorrhiza, and includes mainly protocatechualdehyde, rosmarinic acid, salvianolic acid B, salvianolic acid D, salvianolic acid E, diastereomer of salvianolic acid E, salvianolic acid Y, lithospermic acid and diastereomer of lithospermic acid. It is difficult to quantitatively control the quality of SLI using traditional high performance liquid chromatography due to the highly similar structure of these constituents including three pairs of diastereomers and the lack of commercial resources for most of these constituents as standards. Thus, a highly reproducible, fast, accurate and simple 1H quantitative nuclear magnetic resonance (qNMR) method without the need for calibration curves and complex computation was established by optimizing the solvent system and acquisition parameters to simultaneously determine the nine salvianolic acids and mannitol in SLI. This method was validated and successfully used to determine 10 batches of SLI and the qNMR data were further analyzed with a vector including angle cosine and the partial least squares method for the quality control of SLI. The results indicated that qNMR can be used as a routine method for the quality control of SLI and may have potential in the quantification of diastereomers in other TCM preparations.

Color images steganalysis using rgb channel geometric transformation measures

AbstractIn recent years, information security has received a great deal of attention. To give an example, steganography techniques are used to communicate in a secret and invisible way. Digital color images have become a good medium for digital steganography because of their easy manipulation as carriers via Internet, e‐mails, or used on websites. The main goal of steganalysis is to detect the presence of hidden messages in a digital media. The proposed method is a further extension of the authors' previous work: steganalysis based on color feature correlation and machine learning classification. Fusing features with those obtained from color‐rich models allows increasing the detectability of hidden messages in the color images. Our new proposition uses two types of features, computed between color image channels. The first type of feature reflects local Euclidean transformations, and the second one reflects mirror transformations. These geometric measures are obtained by the sine and cosine of gradient angles between all the color channels. Features are extracted from co‐occurrence correlation matrices of measures. We demonstrate the efficiency of the proposed framework on three steganography algorithms designed to hide messages in images represented in the spatial domain: S‐UNIWARD, WOW, and Synch‐HILL. For each algorithm, we applied a range of different payload sizes. The efficiency of the proposed method is demonstrated by the comparison with the previous authors work and the spatial color‐rich model and color filter array‐aware features for steganalysis. Copyright © 2016 John Wiley & Sons, Ltd.

Flow angle dependent photoacoustic Doppler power spectra under intensity-modulated continuous wave laser excitation

Photoacoustic Doppler (PAD) power spectra showing an evident Doppler shift represent the major characteristics of the continuous wave-excited or burst wave-excited versions of PAD flow measurements. In this paper, the flow angle dependences of the PAD power spectra are investigated using an experiment setup that was established based on intensity-modulated continuous wave laser excitation. The setup has an overall configuration that is similar to a previously reported configuration, but is more sophisticated in that it accurately aligns the laser illumination with the ultrasound detection process, and in that it picks up the correct sample position. In the analysis of the power spectra data, we find that the background power spectra can be extracted by combining the output signals from the two channels of the lock-in amplifier, which is very useful for identification of the PAD power spectra. The power spectra are presented and analyzed in opposite flow directions, at different flow speeds, and at different flow angles. The power spectra at a 90° flow angle show the unique properties of symmetrical shapes due to PAD broadening. For the other flow angles, the smoothed power spectra clearly show a flow angle cosine relationship.

DesR349P Mutation Results in Ultrastructural Disruptions and Compromise of Skeletal Muscle Biomechanics Already at Preclinical Stages in Young Mice before the Onset of Protein Aggregation

Mutations in the extra-sarcomeric protein desmin give rise to protein aggregate myopathies that interfere with cellular function. The R350P mutation is the most common human desminopathy. A murine DesR349P model of human DesR350P desminopathy was previously engineered. In this study, we followed the hypothesis of whether muscle cytoarchitecture and biomechanical properties are already altered in preclinical stages of R349P desminopathy. We applied non-linear second harmonic generation (SHG) and 2-photon fluorescence morphometry analysis to single fibres of fast- and slow-twitch muscle to analyze sarcomeric architecture and nuclear morphology. We found a vast disruption of the lateral sarcomere lattice (indicated by large ‘vernier’ densities) and myofibrillar angular deviations (indicated by cosine angle sums) in single DesR349P fibres. Homozygous fibres were more severely affected than heterozygous. The former showed a marked nuclear pathology with spheric nuclei, increased nuclear number and density. For biomechanics, we assessed active and passive properties. Axial stiffness of DesR349P muscle was much increased, both in fibre and myofibrillar homozygous bundles. Lateral stiffness of the membrane complex in myoblasts was increased in cells heterozygous for DesR349P while homozygous were similar to wt. Caffeine-induced force was compromised in heterozygous fibre bundles but increased in homozygous bundles. This was explained by a marked shift of myofibrillar Ca2+ sensitivity: heterozygous bundles being less sensitive but homozygous more sensitive to Ca2+. This points towards a specific compensatory mechanism in homozygous mutation to cope with otherwise compromised force production already at preclinical stages of desminopathy. Our findings provide a very first complete picture to explain compromised force in heterozygous DesR350P patients: increased muscle stiffness, higher susceptibility towards stretch-induced injury, disrupted myofibrillar alignment and compromised active force production.

Effect of solar zenith angle specification in models on mean shortwave fluxes and stratospheric temperatures

AbstractMany weather and climate models call their radiation schemes only every 3 h, which we show can lead to a stratospheric temperature overestimate of 3–5 K and wavenumber 8 fluctuations in top‐of‐atmosphere (TOA) net shortwave flux around the tropics of amplitude 1.6 W m−2. Solving this problem while retaining a 3h radiation time step requires careful treatment of the cosine of the solar zenith angle, μ0, which appears twice in the calculation of shortwave fluxes, scaling the following: (1) TOA incident flux and (2) the path length of the direct solar beam through the atmosphere. If μ0 is calculated as the average over the radiation time step, rather than at the central time, then the fluctuations are removed, but the stratosphere is still too warm by 2–3 K. It is only if the second μ0 is averaged only over the sunlit part of the radiation time step that the temperature bias is removed.

Evaluation of the accuracy of downward radiative flux observations at the sea surface

Observations of downward radiative flux at the sea surface generally contain uncertainty due to limited numbers of observations and limitations of auxiliary equipment. The lack of shading from direct solar radiation and ventilation systems causes bias or random errors. To evaluate the error of radiation measurements at buoys, downward shortwave and longwave radiative fluxes are compared with International Satellite Cloud Climatology Project (ISCCP), Japanese 55-year Reanalysis (JRA55), and Moderate Resolution Imaging Spectroradiometer (MODIS) retrieved model calculations of 3-h and daytime averages. Cloud masking is evaluated by a combination of MTSAT-1R and in situ observations. Coincident observations from a land-surface station located near the buoy observatories are compared with satellite and reanalysis products. The bias at buoys, compared with retrievals, approximately over- and under-estimate for longwave and shortwave fluxes, respectively. The bias at buoys is larger and smaller than the land by 23–34 W m−2 for longwave and 13–51 W m−2 for shortwave radiation using 3-h averages under clear-sky conditions. The differences in bias decrease when using daytime averages for longwave, but the difference for shortwave increases with daytime averages. To evaluate the effect of environmental factors on buoy observations, we compared rainfall, wind speed, and solar zenith angle with the biases. We found that rainfall and wind speed affect buoy pyrgeometers such that they overestimate the longwave flux. The cosine of solar zenith angle does not cause overestimation for longwave flux, and the effect of dome heating is small. The strong wind causes underestimation of the shortwave radiative flux due to tilting. The effect of wind is reduced when daily averages are used.

Effects of the rotation angle on surface plasmon coupling of nanoprisms.

We studied the effects of relative orientation of bowtie nanostructures on the plasmon resonance both experimentally and theoretically in this work. Specifically, we fabricated gold bowtie nanoantennas with rotated nanoprisms, measured the near-field and the far-field resonance behaviors using Raman spectroscopy and scattering microspectroscopy, and simulated the effects of the rotation angle on the localized surface plasmonic resonance using finite-difference time-domain simulations. In addition to the widely-discussed dipolar resonance in regular bowtie nanostructures, defined as tip-mode resonance in the present study, the excitations of edge-mode resonance were discovered under certain rotation angles of nanoprisms. Because of the resonances of different modes at different wavelengths, two different incident laser sources were used to measure the Raman spectra to provide evidence for the evolution of different resonance modes. Also, both the tip-mode and edge-mode resonances were verified by the simulated charge density distribution and their trends were discussed. Based on the discovered trend, a plasmon protractor was created with a near-exponential decay relationship between the relative resonance wavelength shift and cosine of the rotation angle. A plasmon hybridization model was also proposed for rotated bowties to explain the coupling between nanoprisms during rotation.

Direct and accurate measurement of size dependent wetting behaviors for sessile water droplets.

The size-dependent wettability of sessile water droplets is an important matter in wetting science. Although extensive studies have explored this problem, it has been difficult to obtain empirical data for microscale sessile droplets at a wide range of diameters because of the flaws resulting from evaporation and insufficient imaging resolution. Herein, we present the size-dependent quantitative change of wettability by directly visualizing the three phase interfaces of droplets using a cryogenic-focused ion beam milling and SEM-imaging technique. With the fundamental understanding of the formation pathway, evaporation, freezing, and contact angle hysteresis for sessile droplets, microdroplets with diameters spanning more than three orders of magnitude on various metal substrates were examined. Wetting nature can gradually change from hydrophobic at the hundreds-of-microns scale to super-hydrophobic at the sub-μm scale, and a nonlinear relationship between the cosine of the contact angle and contact line curvature in microscale water droplets was demonstrated. We also showed that the wettability could be further tuned in a size-dependent manner by introducing regular heterogeneities to the substrate.

GPS as a solar observational instrument: Real‐time estimation of EUV photons flux rate during strong, medium, and weak solar flares

AbstractIn this manuscript, the authors show how the Global Navigation Satellite Systems, GNSS (exemplified in the Global Positioning System, GPS), can be efficiently used for a very different purpose from that for which it was designed as an accurate Solar observational tool, already operational from the open global GPS measurements available in real‐time, and with some advantages regarding dedicated instruments onboard spacecraft. The very high correlation of the solar extreme ultraviolet (EUV) photon flux rate in the 26–34 mm spectral band, obtained from the solar EUV monitor instrument onboard the SOHO spacecraft during Solar flares, is shown with the GNSS solar flare activity indicator (GSFLAI). The GSFLAI is defined as the gradient of the ionospheric vertical total electron content rate versus the cosine of the Solar zenith angle in the day hemisphere (which filters out nonsolar over ionization), and it is measured from data collected by a global network of dual frequency GPS receivers (giving in this way continuous coverage). GSFLAI for 60 X class flares, 320 M class flares, and 300 C class flares, occurred since 2001, were directly compared with the EUV solar flux rate data to show existing correlations. It was found that the GSFLAI and EUV flux rate present the same linear relationship for all classes of flares, not only the strong and medium intensity ones, X and M class, as in previous works, but also for the weakest C class solar flares, which is a remarkable result.

Robust object tracking based on local region sparse appearance model

We propose a robust object tracking algorithm based on local region sparse appearance model in this paper. In this algorithm, the object is divided into several sub-regions, and the sparse dictionaries are obtained by clustering in each sub-region. Therefore spatial structure information of the object can be captured well, and the change of object appearance can be also resisted effectively. First, the object is divided into many small patches. Then the object is divided into several sub-regions according to patch distribution again. The establishment of object dictionary base is based on combination of the dictionaries from all the sub-regions, and then space alignment between different parts of the object can be achieved. Meanwhile, noise removal and other operations in the existing sparse reconstruction error maps are performed to retain valuable information. In the updating framework, a novel flexible template set update mechanism is introduced in this paper. In this update mechanism, valuable object samples will be put into the template set. If samples are not valuable, they should not be put into the template set, even when the template set is not full. Then we use patch sparse coefficient histogram of updated templates to extract time domain information of the object in the form of weighted sum. Therefore, it can provide a reliable template basis for obtaining good candidate object. In addition, when tracking result deviates from the actual position of the object, we use a dynamic sub-region resampling method based on cosine angle to correct the position deviation timely. Therefore this method can effectively prevent the object from being completely lost. Both qualitative and quantitative evaluations on challenging video sequences demonstrate that the proposed tracking algorithm performs favorably against several state-of-the-art methods.