Model for the interpretation of hyperspectral remote-sensing reflectance

Remote-sensing reflectance is easier to interpret for the open ocean than for coastal regions because the optical signals are highly coupled to the phytoplankton (e.g., chlorophyll) concentrations. For estuarine or coastal waters, variable terrigenous colored dissolved organic matter (CDOM), suspended sediments, and bottom reflectance, all factors that do not covary with the pigment concentration, confound data interpretation. In this research, remote-sensing reflectance models are suggested for coastal waters, to which contributions that are due to bottom reflectance, CDOM fluorescence, and water Raman scattering are included. Through the use of two parameters to model the combination of the backscattering coefficient and the Q factor, excellent agreement was achieved between the measured and modeled remote-sensing reflectance for waters from the West Florida Shelf to the Mississippi River plume. These waters cover a range of chlorophyll of 0.2-40 mg/m(3) and gelbstoff absorption at 440 nm from 0.02-0.4 m(-1). Data with a spectral resolution of 10 nm or better, which is consistent with that provided by the airborne visible and infrared imaging spectrometer (AVIRIS) and spacecraft spectrometers, were used in the model evaluation.     

Effects of a nonuniform vertical profile of chlorophyll concentration on remote-sensing reflectance of the ocean

Numerical simulations of radiative transfer were used to examine the effects of a nonuniform vertical profile of the inherent optical properties of the water column associated with the vertical profile of chlorophyll concentration, Chl(z), on the spectral remote-sensing reflectance, Rrs(gamma), of the ocean. Using the Gaussian function that describes the Chl(z) profile, we simulated a relatively broad range of open-ocean conditions characterized by the presence of a subsurface Chl maximum at depths greater than or equal to 20 m. The simulations for a vertically nonuniform Chl(z) were compared with reference simulations for a homogeneous ocean whose Chl was identical to the surface Chl of inhomogeneous cases. The range of values for the Gaussian parameters that produce significant differences in Rrs(gamma) (> 5%) was determined. For some vertical structures of Chl(z) considered, the magnitude of Rrs(gamma) and the blue-to-green band ratios of Rrs(gamma) differ significantly from the reference values of homogeneous ocean (> 70% in extreme cases of low surface chlorophyll of 0.02 mg m(-3) and shallow pigment maximum at 20 m). The differences are small or negligible when the nonuniform profiles are characterized by a surface Chl greater than 0.4 mg m(-3) or a depth of Chl maximum greater than 45 m (65 m in extremely clear waters with a surface Chl of 0.02 mg m(-3) or less). The comparison of modeling results with the current algorithm for retrieving the global distribution of chlorophyll from satellite imagery of ocean color suggests that strong effects of the subsurface chlorophyll maximum on reflectance at low surface chlorophyll concentrations can lead to a severalfold overestimation in the algorithm-derived surface chlorophyll. Examples of field data from the Sea of Japan and the north polar Atlantic Ocean are used to illustrate various nonuniform pigment profiles and their effect on the blue-to-green ratio of Rrs(gamma).     

Estimation of the remote-sensing reflectance from above-surface measurements

The remote-sensing reflectance R(rs) is not directly measurable, and various methodologies have been employed in its estimation. I review the radiative transfer foundations of several commonly used methods for estimating R(rs), and errors associated with estimating R(rs) by removal of surface-reflected sky radiance are evaluated using the Hydrolight radiative transfer numerical model. The dependence of the sea surface reflectance factor rho, which is not an inherent optical property of the surface, on sky conditions, wind speed, solar zenith angle, and viewing geometry is examined. If rho is not estimated accurately, significant errors can occur in the estimated R(rs) for near-zenith Sun positions and for high wind speeds, both of which can give considerable Sun glitter effects. The numerical simulations suggest that a viewing direction of 40 deg from the nadir and 135 deg from the Sun is a reasonable compromise among conflicting requirements. For this viewing direction, a value of rho approximately 0.028 is acceptable only for wind speeds less than 5 m s(-1). For higher wind speeds, curves are presented for the determination of rho as a function of solar zenith angle and wind speed. If the sky is overcast, a value of rho approximately 0.028 is used at all wind speeds.     

Retrieval of chlorophyll from remote-sensing reflectance in the china seas

The East China Sea is a typical case 2 water environment, where concentrations of phytoplankton pigments, suspended matter, and chromophoric dissolved organic matter (CDOM) are all higher than those in the open oceans, because of the discharge from the Yangtze River and the Yellow River. By using a hyperspectral semianalytical model, we simulated a set of remote-sensing reflectance for a variety of chlorophyll, suspended matter, and CDOM concentrations. From this simulated data set, a new algorithm for the retrieval of chlorophyll concentration from remote-sensing reflectance is proposed. For this method, we took into account the 682-nm spectral channel in addition to the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) channels. When this algorithm was applied to a field data set, the chlorophyll concentrations retrieved through the new algorithm were consistent with field measurements to within a small error of 18%, in contrast with that of 147% between the SeaWiFS ocean chlorophyll 2 algorithm and the in situ observation.     

Model of remote-sensing reflectance including bidirectional effects for case 1 and case 2 waters

A remote-sensing reflectance model based on a lookup table is proposed for use in analyzing satellite ocean color data in both case 1 and case 2 waters. The model coefficients are tabulated for grid values of three angles--solar zenith, sensor zenith, and relative azimuth--to take account of directional variation. This model also requires, as input, a phase function parameter defined by the contribution of suspended particles to the backscattering coefficient. The model is generated from radiative transfer simulations for a wide range of inherent optical properties that cover both case 1 and 2 waters. The model uncertainty that is due to phase function variability is significantly reduced from that in conventional models. Bidirectional variation of reflectance is described and explained for a variety of cases. The effects of wind speed and cloud cover on bidirectional variation are also considered, including those for the fully overcast case in which angular variation can still be considerable (approximately 10%). The implications for seaborne validation of satellite-derived water-leaving reflectance are discussed.     

Diffuse reflectance of oceanic waters. III. Implication of bidirectionality for the remote-sensing problem

The upwelling radiance field beneath the ocean surface and the emerging radiance field are not generally isotropic. Their bidirectional structure depends on the illumination conditions (the Sun's position in particular) and on the optical properties of the water body. In oceanic case 1 waters, these properties can be related, for each wavelength λ, to the chlorophyll (Chl) concentration. We aim to quantify systematically the variations of spectral radiances that emerge from an ocean with varying Chl when we change the geometric conditions, namely, the zenith-Sun angle, the viewing angle, and the azimuth difference between the solar and observational vertical planes. The consequences of these important variations on the interpretation of marine signals, as detected by a satelliteborne ocean color sensor, are analyzed. In particular, the derivation of radiometric quantities, such as R (λ), the spectral reflectance, or [ L(w)(λ)](N), the normalized water-leaving radiance that is free from directional effects, is examined, as well as the retrieval of Chl. We propose a practical method that is based on the use of precomputed lookup tables to provide values of the f/Q ratio in all the necessary conditions[ f relates (R(λ) to the backscattering and absorption coefficients, whereas Q is the ratio of upwelling irradiance to any upwelling radiance]. The f/Q ratio, besides being dependent on the geometric configuration (the three angles mentioned above), also varies with λ and with the bio-optical state, conveniently depicted by Chl. Because Chl is one of the entries for the lookup table, it has to be derived at the beginning of the process, before the radiometric quantities R(λ) or [L(W)(λ)](N) can be produced. The determination of Chl can be made through an iterative process, computationally fast, using the information at two wavelengths. In this attempt to remove the bidirectional effect, the commonly accepted view relative to the data-processing strategy is somewhat modified, i.e., reversed, as the Chl index becomes a prerequisite parameter that must be identified prior to the derivation of the fundamental radiometric quantities at all wavelengths.     

Analytic remote-sensing optical algorithms requiring simple and practical field parameter inputs

Spectral in-water measurements of downward irradiance (E(d)), upward irradiance (E(u)), and nadir radiance (L(u)) are sufficient to calculate the scalar irradiances E(0), E(0d), and E(0u), the average cosines mu, mu(d), and mu(u), the light absorption coefficient a, the backscattering coefficient b(b), and the so-called f factor that relates to R, a, and b(b). The solar elevation of 42 degrees is a special case in which mu(d) is independent of all variables except solar elevation. The algorithms are valid for solar elevations between 12 degrees and 81 degrees for horizontally stratified clear and turbid deep waters.     

Multiparameter retrieval of water optical properties from above-water remote-sensing reflectance using the simulated annealing algorithm

Water-leaving radiance, measured just above the ocean surface, contains important information about near-surface or subsurface processes that occur on or below the deep ocean and coastal water. As such, retrieving seawater inherent optical properties (IOPs) is an important step to determining water type, subsurface light field, turbidity, pigment concentration, and sediment loading. However, the retrieval (or inversion) of seawater IOPs from just above water radiance measurements is a multiparameter nonlinear problem that is difficult to solve by conventional optimization methods. The applicability of the simulated annealing algorithm (SA) is explored as a nonlinear global optimizer to solve this multiparameter retrieval problem. The SA algorithm is combined with widely known semianalytical relations for seawater's IOPs to parameter invert these properties from simulated and measured water-leaving reflectance spectra. Furthermore, given the versatility of the SA algorithm, the scheme is extended to retrieve water depth from input reflectance data. Extensive tests and comparisons with in situ and simulated data sets compiled by the International Ocean-Color Coordinating Group are presented. Field data include reflectance spectra acquired with a handheld GER 1500 spectroradiometer and absorption measurements, performed with the AC-9 instrument on waters around Singapore's nearby islands.     

Effects of electron scatterings on thermal conductivity of thin metal films

In the framework of the statistical conduction models for polycrystalline and monocrystalline metal films, the thermal conductivity due to electron transport is calculated under the further assumption that the electron relaxation time in the bulk material depends on electron energy. It is shown that the same function describes the size effect in thermal and electrical conductivities, and that the Wiedemann-Franz law holds for any electron scattering in thin metal films.     

Directional statistics-based reflectance model for isotropic bidirectional reflectance distribution functions

We introduce a novel parametric bidirectional reflectance distribution function (BRDF) model that can accurately encode a wide variety of real-world isotropic BRDFs with a small number of parameters. The key observation we make is that a BRDF may be viewed as a statistical distribution on a unit hemisphere. We derive a novel directional statistics distribution, which we refer to as the hemispherical exponential power distribution, and model real-world isotropic BRDFs as mixtures of it. We derive a canonical probabilistic method for estimating the parameters, including the number of components, of this novel directional statistics BRDF model. We show that the model captures the full spectrum of real-world isotropic BRDFs with high accuracy, but a small footprint. We also demonstrate the advantages of the novel BRDF model by showing its use for reflection component separation and for exploring the space of isotropic BRDFs.     

A feature parameter of optical diffuse reflectance on normal and malignant human breast tissue

Light that is delivered to the tissue surface undergoes multiple elastic scattering and absorption; part of it returns to the surface as optical diffuse reflectance, carrying quantitative information on the structure and composition of the measured tissue. In this paper, we utilized a well tested and publicly available Monte Carlo simulation software to simulate optical diffuse reflectance on normal and malignant human breast tissue in the visible wavelength region 450–650 nm. Based on the Monte Carlo simulation results, we discovered a feature parameter of optical diffuse reflectance on the simulated tissue surface. Normal and malignant human breast tissue may be discriminated by the value of the feature parameter. The values of the feature parameter are shown for normal and malignant human breast tissue in the visible wavelength region 450–650 nm. Experiments with tissue phantoms showed that the feature parameter is correlated to the component of the phantom. So the feature parameter is useful for the non-invasive optical diagnosis of biological tissue.     

Effects of model scale and particle size on micro-mechanical properties and failure processes of rocks—A particle mechanics approach

A numerical procedure to determine the equivalent micro-mechanical properties of intact rocks is presented using a stochastic representative elementary volume (REV) concept and a particle mechanics approach. More than 200 models were generated in square regions with side lengths varying from 1 to 10 cm, using the Monte Carlo simulation technique. Generated particle models were then used for the calculation of equivalent micro-mechanical properties. Results with a core sample of diorite from Äspö, Sweden, show that the variance of the calculated values of mechanical properties decrease significantly as the side lengths of particle models increase, reaching a REV side length about 5 cm with an acceptable variation of 5%, which is equal to the minimum diameter of rock specimen for uniaxial compressive tests suggested by ISRM. The complete stress–strain curve of the diorite rock sample was predicted under uniaxial compression, as the basis for evaluating the damage and failure processes. The unique contribution of this paper is its study on impacts of sample size and particle size distributions on mechanical behaviour of rocks when particle mechanics approaches are used.     

基于粒子群算法的改进模态参数识别方法

针对一类多模态振动衰减信号的模态参数识别,结合奇异值分解(singular value decomposition,SVD)、解析模态分解(analytical mode decomposition,AMD)、自回归功率谱和粒子群优化(particle swarm optimization,PSO)算法,提出了一种改进...     

Bidirectional reflectance of oceanic waters: accounting for Raman emission and varying particle scattering phase function

The bidirectionality of the upward radiance field in oceanic case 1 waters has been reinvestigated by incorporation of revised parameterizations of inherent optical properties as a function of the chlorophyll concentration (Chl), considering Raman scattering and making the particle phase function shape (beta(rho)) continuously varying along with the Chl. Internal consistency is thus reached, as the decrease in backscattering probability (for increasing Chl) translates into a correlative change in beta(rho). The single particle phase function (previously used) precluded a realistic assessment of bidirectionality for waters with Chl > 1 mg m(-3). This limitation is now removed. For low Chl, Raman emissions significantly affect the radiance field. For moderate Chl (0.1-1 mg m(-3)), new and previous bidirectional parameters remain close. The ocean reflectance anisotropy has implications in ocean color remote-sensing problems, in derivation of coherent water-leaving radiances, in associated calibration-validation activities, and in the merging of data obtained under various geometrical configurations.     

A finite deformation method for discrete modeling: particle rotation and parameter calibration

We present a finite deformation method for 3-D discrete element modeling. In this method particle rotation is explicitly represented using quaternion and a complete set of interactions is permitted between two bonded particles, i.e., normal and tangent forces, rolling and torsional torques. Relative rotation between two particles is decomposed into two sequence-independent rotations, such that an overall torsional and rolling angle can be distinguished and torques caused by relative rotations are uniquely determined. Forces and torques are calculated in a finite deformation fashion, rather than incrementally. Compared with the incremental methods our algorithm is numerically more stable while it is consistent with the non-commutativity of finite rotations. We study the macroscopic elastic properties of a regularly arranged 2-D and 3-D lattice. Using a micro-to-macro approach based on the existence of a homogeneous displacement field, we study the problem of how to choose the particle-scale parameters (normal, tangent, rolling and torsional stiffness) given the macroscopic elastic parameters and geometry of lattice arrangement. The method is validated by reproducing the wing crack propagation and the fracture patterns under uniaxial compression. This study will provide a theoretical basis for the calibration of the DEM parameters required in engineering applications.     

A discrete model for particle deposition

Simulation of deposit growth on a two-dimensional substrate was studied based on a new model that tracks individual cubic particles as they form a deposit structure. The present model is an extension of the classical ballistic deposition model. Effects of three different parameters were studied. These include an attraction parameter that is a measure of the particle to particle attractions, an interaction length within which the particles are assumed to influence and be influenced by surrounding particles, and allowed sticking positions (face-face, edge-edge and corner-corner) that favor particular growth directions. Structures with widely varying properties were obtained using this model. The three parameters were found to have considerable effect on the structure including indications of morphological phase transformations. A new property of the system (saturated roughness/deposit growth rate) was identified that can classify the different types of growth into a single type.     

Random particle model for concrete based on Delaunay triangulation

In this paper, an efficient simulation method for obtaining a random particle model for concrete is outlined. First, the ‘take-and-place method’ and its extension, the ‘directed searching process’, are discussed briefly and the shortcomings are indicated. A new method, the ‘divide-and-fill method’, appears to be more convenient, especially when only a small computer is used. In this simulation method the available space (two-dimensional) is divided in separate areas, using a Delaunay triangulation. These areas are filled with particles taking into account a given grading curve and gravel content. Comparison with physical concrete sections, obtained by means of image analysis, shows that the results of this method closely represent reality. The divide-and-fill method also yields a finite-element mesh in a quasi-automatic way.     

Laser remote-sensing system analysis for search and rescue

We develop a general model of a laser remote-sensing system for search and rescue using targets marked with fluorescent dye. The dye fluoresces at a longer peak wavelength than the incident radiation, enabling a dye-covered target to be distinguished from the unshifted ground echo by the search system. The principal result is a simple expression derived for the average laser power required to search at a particular rate given a required ground energy density. A similar expression is applicable to imaging lidar systems. The example system shown indicates that active probing for lost planes may be practical.