ALBERT

Description: The L-band passive and active microwave geophysical model functions (GMFs) of ocean surface winds from the Aquarius data are derived. The matchups of Aquarius data with the Special Sensor Microwave Imager (SSM/I) and National Centers for Environmental Prediction (NCEP) winds were performed and were binned as a function of wind speed and direction. The radar HH GMF is in good agreement with the PALSAR GMF. For wind speeds above 10 $hbox{m}cdothbox{s}^{-1}$ , the L-band ocean backscatter shows positive upwind–crosswind (UC) asymmetry; however, the UC asymmetry becomes negative between about 3 and 8 $hbox{m}cdothbox{s}^{-1}$ . The negative UC (NUC) asymmetry has not been observed in higher frequency (above C-band) GMFs for ASCAT or QuikSCAT. Unexpectedly, the NUC symmetry also appears in the L-band radiometer data. We find direction dependence in the Aquarius $T_{rm BV}$ , $T_{rm BH}$ , and third Stokes data with peak-to-peak modulations increasing from about a few tenths to 2 K in the range of 10–25- $hbox{m}cdothbox{s}^{-1}$ wind speed. The validity of the GMFs is tested through application to wind and salinity retrieval from Aquarius data using the combined active–passive algorithm. Error assessment using the triple collocation analyses of SSM/I, NCEP, and Aquarius winds indicates that the retrieved Aquarius wind speed accuracy is excellent, with a random error of about 0.75 $hbox{m}cdothbox{s}^{-1}$ . The wind direction retrievals also appear reasonable and accurate above 10 $hbox{m}cdothbo- {s}^{-1}$ . The results of the error analysis indicate that the uncertainty of the GMFs for the wind speed correction of vertically polarized brightness temperatures is about 0.14 K for wind speed up to 10 $hbox{m}cdothbox{s}^{-1}$ .

Description: The local oscillators (LOs) of the Soil Moisture and Ocean Salinity mission payload are used to shift the operating frequency of the 72 receivers to an optimal intermediate frequency needed for the signal processing. The LO temperature variations produce phase errors in the visibility, which result in a blurring of the reconstructed brightness temperature (Tb) image. At the end of the commissioning phase, it was decided to calibrate the LO every 10 min while waiting for a more in-depth analysis. During short periods of time, the LO calibration has been performed every 2 min to assess the impact of a higher calibration rate on the quality of the data. In this paper, by means of a decimation experiment, the relative errors of 6- and 10-min calibration interval data sets are estimated using the 2 min as a reference. A noticeable systematic across- and along-track pattern of amplitude $pm$ 0.3 K is observed for Tb differences between 10 and 2 min, whereas this is reduced between 6 and 2 min. A simulation experiment confirms that the nature of such systematic pattern is due to the visibility phase errors induced by the LO calibration rate. Such pattern is propagated into the sea surface salinity (SSS) retrievals. Overall, the SSS error increase (relative to the 2 min SSS data) is about 0.39 and 0.14 psu for the 10- and 6-min data sets, respectively. This paper shows that a LO calibration rate of at least 6 min would noticeably improve the SSS retrievals.

Description: Radiometric measurements could provide continuous information about atmospheric conditions. In this paper, a sky status indicator (SSI) is proposed as a real-time recognition criterion for the detection, in particular, of the presence of rain events along the propagation path. The computation of the SSI is based on ground-based brightness temperature measurements, at 23.8 and 31.4 GHz, collected in Cabauw, Netherlands, in 2009 by the ESA Atmospheric Propagation and Profiling System (ATPROP) multichannel radiometer. A validation analysis is carried out between simulated data, which are computed by applying the radiative transfer equation to a database of radiosonde profiles collected in De Bilt, Netherlands, by the Royal Netherlands Meteorological Institute, and two data sets of radiometric observations at two elevation angles ( $theta$ equal to 90 $^{circ}$ and $theta$ equal to 69.6 $^{circ}$ ). The analysis based on SSI probability distribution functions has allowed for calculation of the boundary threshold values that are able to discriminate the status of the sky. Furthermore, performances of the SSI were validated against rainfall measurements collected at the ground by a rain gauge located near the ESA ATPROP multichannel radiometer.

Description: A new cloud dynamics and radiation database (CDRD) precipitation retrieval algorithm for satellite passive microwave (PMW) radiometer measurements has been developed. It represents a modification to and an improvement upon the conventional cloud radiation database (CRD) algorithms, which have always been prone to ambiguity. This part 2 paper of a series describes the methodology of the algorithm and the modeling verification analysis involved in creating a synthetic CDRD database for the Europe/Mediterranean basin region. This is followed by a proof-of-concept analysis, which demonstrates that the underlying CDRD theory based on use of meteorological parameters for reducing retrieval ambiguity is valid. This paper uses a regional/mesoscale model, applied in cloud resolving model (CRM) mode, to produce a large set of numerical simulations of precipitating storms and extended precipitating systems. The simulations are used for selection of millions of meteorological/microphysical vertical profiles within which surface rainfall is identified. For each of these profiles, top-of-atmosphere brightness temperature (TB) vectors are calculated (the vector dimension associated with the number of relevant cm–mm wavelengths and polarizations), based on an elaborate radiative-transfer equation (RTE) model system (RMS) coupled to the CRM. This entire body of simulation information is organized into the CDRD database, then used as a priori knowledge to guide a physical Bayesian retrieval algorithm in obtaining rainfall and associated precipitation parameters from the PMW satellite observations. We first prove the physical validity of our CRM-RMS simulations, by showing that the simulated TBs are in close agreement with observations. Agreement is demonstrated using dual-channel-frequency TB manifold sections, which quantify the degree of overlap between the simulated and observed TBs extracted from the full manifolds. Nevertheless, the salient result of this paper is a pro- f that the underlying CDRD theory is valid, found by combining subdivisions of the invoked meteorological parameter ranges of values and showing that such meteorological partitioning associates itself with distinct microphysical profiles. It is then shown that these profiles give rise to similar TB vectors, proving the existence of ambiguity in a CRD-type algorithm. Finally, we show that the CDRD methodology provides significant improvements in reducing retrieval ambiguity and retrieval error, especially for land surface backgrounds where contrasts are typically small between the rainfall TB signatures and surface emission signatures.

Description: The results of remote sensing temperature profiles measurements within a 0–600-m altitude range and total water content measurements during total (Kislovodsk, 2006; Novosibirsk, 2008) and partial (Moscow, 2011) solar eclipses, using microwave radiometers are presented. Initially, continuous data on temperature profiles are obtained at different altitudes before, during, and after total solar eclipses, using two single channel elevation scanning microwave temperature profilers. Terrestrial consequences of solar eclipses (especially total ones) are quite noticeable and important. Solar eclipses support unique, specific conditions, which gives the opportunity for various meteorological research. The most important indicator of thermodynamic processes occurring during solar eclipses is air temperature at different altitudes in the atmospheric boundary layer (ABL). The ABL temperature depends, in general, on the flux of solar radiation and some features of the ground (albedo, absorptivity, and emissivity) and the air (humidity). Temperature profile measurements are accompanied by solar radiation (with net-radiometer) and total water vapor (with microwave radiometers) measurements. The observation results of this paper will contribute detailed model calculations for clarifying meteorological effects of solar eclipses. Observations of the next total solar eclipse over Russia (August 12, 2026) can be used to verify our observational results.

Description: Snow grain size is the snowpack parameter that most affects the microwave snow emission. The specific surface area (SSA) of snow is a metric that allows rapid and reproducible field measurements and that well represents the grain size. However, this metric cannot be used directly in microwave snow emission models (MSEMs). The aim of this paper is to evaluate the suitability and the adaptations required for using the SSA in two MSEMs, i.e., the Dense Media Radiative Theory-Multilayer model (DMRT-ML) and the Helsinki University of Technology model (HUT n-layer), based on in situ radiometric measurements. Measurements of the SSA, using snow reflectance in the short-wave infrared, were taken at 20 snowpits in various environments (e.g., grass, tundra, and dry fen). The results show that both models required a scaling factor for the SSA values to minimize the root-mean-square error between the measured and simulated brightness temperatures. For DMRT-ML, the need for a scaling factor is likely due to the oversimplified representation of snow as spheres of ice with a uniform radius. We hypothesize that the need for a scaling factor is related to the grain size distribution of snow and the stickiness between grains. For HUT n-layer, using the SSA underestimates the attenuation by snow, particularly for snowpacks with a significant amount of depth hoar. This paper provides a reliable description of the grain size for DMRT-ML, which is of particular interest for the assimilation of satellite passive microwave data in snow models.

Description: Ice lens formation, which follows rain on snow events or melt-refreeze cycles in winter and spring, is likely to become more frequent as a result of increasing mean winter temperatures at high latitudes. These ice lenses significantly affect the microwave scattering and emission properties, and hence snow brightness temperatures that are widely used to monitor snow cover properties from space. To understand and interpret the spaceborne microwave signal, the modeling of these phenomena needs improvement. This paper shows the effects and sensitivity of ice lenses on simulated brightness temperatures using the microwave emission model of layered snowpacks coupled to a soil emission model at 19 and 37 GHz in both horizontal and vertical polarizations. Results when considering pure ice lenses show an improvement of 20.5 K of the root mean square error between the simulated and measured brightness temperature (Tb) using several in situ data sets acquired during field campaigns across Canada. The modeled Tbs are found to be highly sensitive to the vertical location of ice lenses within the snowpack.

Description: In recent years, there has been growing interest on the part of the remote sensing community in using the Antarctic area for calibrating and validating data of low-frequency satellite-borne microwave radiometers. In particular, the East Antarctic Plateau appears to be suited for this purpose. The reasons for this interest are the size, structure, spatial homogeneity, and thermal stability of this area. This is particularly interesting for low-frequency microwave radiometers since, due to the low extinction of dry snow, the upper ice-sheet layer is almost transparent and the brightness temperature variability is therefore extremely small. In the context of calibration and validation activities of the European Space Agency's Soil Moisture and Ocean Salinity (SMOS) satellite, an experiment called DOMEX-2, which included radiometric L-band measurements, was carried out at the Italian–French base of Concordia located at Dome C in the East Antarctic Plateau from December 2008 to December 2010. Ground measurements (i.e., snow temperature at different depths, snow structure, meteorological data, etc.) were also collected during the experiment. This paper presents information on the experimental campaign, the characteristics of the radiometric measurements, and the main results. A comparison with SMOS data is also presented.

Description: The 12th Specialist Meeting on Microwave Radiometry and Remote Sensing of the Environment (MicroRad 2012) was held at Villa Mondragone, University of Rome "Tor Vergata," near Frascati, Italy, on March 5-9, 2012. The objective of MicroRad 2012 was to provide an open forum to report and discuss recent advances in the field of microwave radiometry, particularly with application to remote sensing of the environment. The meeting was highly successful, with more than 120 attendees representing 20 countries. There were 76 oral presentations and more than 40 posters. From the papers presented at MicroRad 2012 and others submitted specifically for this special issue, 12 were selected for inclusion in the special issue. The papers were carefully peer reviewed with the usual standards of the IEEE TGRS. As is evident from the table of contents, these papers span a broad range of microwave radiometry and remote sensing applications and reflect the interest in MicroRad and the vitality of research in this area.

Description: Vicarious cold calibration in the frequency range of 85–92 GHz is analyzed. Vicarious cold calibration cannot be applied at these frequencies as easily as at lower frequencies due to greater sensitivity to water vapor and hydrometeor scattering. The effects of that sensitivity are mitigated by selective filtering of the high-frequency brightness temperatures (TBs) to remove those data where large amounts of water vapor and/or hydrometeor scattering are present. Potential filtering algorithms are presented, and the performance of each with respect to vicarious cold calibration TB stability is characterized. A scattering-based precipitation filter that utilizes a combination of both the lower frequencies from 19 to 37 GHz and the frequencies from 85 to 92 GHz is shown to be the most effective and easily implemented filter. For horizontal polarization, the theoretical minimum TB at the higher frequencies occurs at an unphysically high sea surface temperature (SST), which makes the vicarious cold statistic more sensitive to the population of actual SST values as well as the higher amounts of water vapor associated with warm SSTs. The statistic is stabilized in this case by considering the difference between observed and simulated vicarious cold TBs. Intercalibration between two radiometers using the vicarious cold calibration double difference method at high frequencies is shown to be greatly improved when using the precipitation filter.

Description: After 2.5 years of the Soil Moisture and Ocean Salinity (SMOS) mission, the characterization of residual instrumental systematic errors in the measured brightness temperatures $(T_{B})$ is still rather poor. This, in turn, negatively impacts the sea surface salinity retrievals and, as such, notably limits the mission's success. The error mitigation methodology currently used operationally, the so-called Ocean Target Transformation (OTT), mixes both instrumental and model-induced errors. In this paper, it is proposed to distinguish errors by their type of impact on the $T_{B}$ images: mean brightness level, incidence angle dependence, and azimuth angle dependence. A new approach to characterize the azimuth-dependent errors is proposed. First, a careful data selection strategy is applied. Then, an empirically fitted model, which only accounts for the $T_{B}$ incidence angle dependence, is subtracted from the mean $T_{B}$ images of the selected data sets to estimate the systematic antenna-frame errors. The robustness of this methodology is assessed through the estimated anomaly pattern stability when computed for different geophysical conditions, periods of time, and latitudinal bands. The residual variability ranges from 0.03 K to 0.14 K, whereas the OTT variability is about 0.5 K. The new method is forward model independent and generic. It can therefore be applied to estimate the antenna-frame systematic errors over land and ice. Moreover, it proves to be very effective in separating different sources of error and can therefore be used to further characterize other error components and improve the various SMOS forward model terms.

Description: Terrain Observation by Progressive Scans (TOPS) synthetic aperture radar (SAR) and spotlight SAR are advanced SAR imaging modes for wide range swath and high resolution. In order to obtain a wider range coverage, azimuth multichannel is introduced in the literature. Since the azimuth bandwidth of beam steering SAR (BS-SAR; spotlight SAR, sliding spotlight SAR, or TOPS SAR) is much greater than that of a stripmap SAR, a signal reconstruction algorithm used for multichannel stripmap SAR may not be effective for multichannel BS-SAR. In this paper, a multichannel full-aperture azimuth processing algorithm is proposed for a BS-SAR. The key of this algorithm lies in the beam and the azimuth bandwidth compressions of multichannel signals in the Doppler–array and slow time–angle planes, respectively. Through compression processing, the beamwidth and the azimuth bandwidth are smaller than the available angle and equivalent pulse repeating frequency , respectively. Then, an improved post-Doppler STAP method is proposed to recover a 2-D spectrum. With the recovered signal, further processing can be utilized to focus the multichannel signal. Simulation and real data results show the effectiveness of the proposed algorithm.

Description: Remote sensing image fusion can integrate the spatial detail of panchromatic (PAN) image and the spectral information of a low-resolution multispectral (MS) image to produce a fused MS image with high spatial resolution. In this paper, a remote sensing image fusion method is proposed with sparse representations over learned dictionaries. The dictionaries for PAN image and low-resolution MS image are learned from the source images adaptively. Furthermore, a novel strategy is designed to construct the dictionary for unknown high-resolution MS images without training set, which can make our proposed method more practical. The sparse coefficients of the PAN image and low-resolution MS image are sought by the orthogonal matching pursuit algorithm. Then, the fused high-resolution MS image is calculated by combining the obtained sparse coefficients and the dictionary for the high-resolution MS image. By comparing with six well-known methods in terms of several universal quality evaluation indexes with or without references, the simulated and real experimental results on QuickBird and IKONOS images demonstrate the superiority of our method.

Description: The multichromatic analysis (MCA) uses interferometric pairs of SAR images processed at range subbands and explores the phase trend of each pixel as a function of the different central carrier frequencies to infer absolute optical path difference. This approach allows retrieving unambiguous height information on selected pixels, potentially solving the problem of spatial phase unwrapping, which is instead critical in the standard monochromatic processing. The method, based on concepts originally introduced by Madsen and Zebker, has been developed in previous work both theoretically and through simulations. This paper presents the first MCA experimental validation of the procedure, through application to a wideband SAR single-pass interferometric data set acquired by the AES-1 airborne sensor. An evaluation of the impact of the MCA processing parameters on the height estimation performances is obtained through a parametric analysis. The results confirm the indications derived by the theoretical analysis, demonstrating the feasibility of the MCA absolute phase measurement, provided that a sufficient bandwidth is available.

Description: This paper proposes a new semisupervised dimension reduction (DR) algorithm based on a discriminative locally enhanced alignment technique. The proposed DR method has two aims: to maximize the distance between different classes according to the separability of pairwise samples and, at the same time, to preserve the intrinsic geometric structure of the data by the use of both labeled and unlabeled samples. Furthermore, two key problems determining the performance of semisupervised methods are discussed in this paper. The first problem is the proper selection of the unlabeled sample set; the second problem is the accurate measurement of the similarity between samples. In this paper, multilevel segmentation results are employed to solve these problems. Experiments with extensive hyperspectral image data sets showed that the proposed algorithm is notably superior to other state-of-the-art dimensionality reduction methods for hyperspectral image classification.

Description: This paper presents a new framework for the development of generalized composite kernel machines for hyperspectral image classification. We construct a new family of generalized composite kernels which exhibit great flexibility when combining the spectral and the spatial information contained in the hyperspectral data, without any weight parameters. The classifier adopted in this work is the multinomial logistic regression, and the spatial information is modeled from extended multiattribute profiles. In order to illustrate the good performance of the proposed framework, support vector machines are also used for evaluation purposes. Our experimental results with real hyperspectral images collected by the National Aeronautics and Space Administration Jet Propulsion Laboratory's Airborne Visible/Infrared Imaging Spectrometer and the Reflective Optics Spectrographic Imaging System indicate that the proposed framework leads to state-of-the-art classification performance in complex analysis scenarios.

Description: For typical scanning microwave radiometers, a significant source of calibration error arises from thermal gradients on the hot load. Even when direct or reflected solar illumination is blocked, hot load gradients arise from thermal coupling between the target and the surface facing the target which is heated and cooled as the instrument orbits the earth. For the GlobalL Precipitation Measurement (GPM) Microwave Imager (GMI), a rotating metal annular ring called the “hot load tray” serves to guard the hot load against solar intrusion, and is the surface immediately facing the hot load during the majority of the scan. The planned GMI calibration algorithm corrects for the target gradients induced by thermal coupling between the hot load tray and hot load. The correction uses an empirically derived relationship between the target gradient and the temperature differential between the target and the tray. The correction is derived using target-level and GMI system-level calibration testing. The dual calibration of GMI, in connection with thermal vacuum calibration measurements, is a key aid to determining and correcting the hot load gradients.

Description: Satellite images have long been used to study surface manifestations of internal waves (IWs). More recently, marine X-band radar data have been employed to retrieve IW packet parameters. Marine radars have the advantage over satellite systems that their high temporal resolution enables the study of the IW evolution. Until today, no method to automatically detect IW surface signatures in marine radar data has been suggested. In this paper, we present a new fully automated tool to retrieve IW signatures from marine radar image sequences. First, after various preprocessing steps, the IW packet velocity is determined using a combination of localized Radon transform and cross-correlation techniques. Temporal averaging of the marine radar data significantly enhances the IW signatures. The knowledge of the IW packet velocity is used to correct for the IW motion, enabling us to extend the averaging period, which further enhances the IW signal. An IW-motion correction is necessary because, otherwise, the IW signal would become smeared if the averaging period were much longer than the time it takes the IW to propagate between radar resolution cells. The IW-enhanced images are then utilized for the IW signature analysis. Here, we identify local backscatter peaks and exploit the marine radar's high temporal resolution to distinguish signal from noise. The resulting series of IW soliton maps provides information on changes in soliton wavelength, velocity, and backscatter intensity. Our marine radar IW signature analysis tool therefore offers a great opportunity of studying the spatiotemporal evolution of IWs as they grow and decay.

Description: The potential of satellite passive microwave sensors to provide quantitative information about near-source volcanic ash cloud parameters is assessed. To this aim, ground-based microwave weather radar and spaceborne microwave radiometer observations are used together with forward-model simulations. The latter are based on 2-D simulations with the numerical plume model Active Tracer High-Resolution Atmospheric Model (ATHAM), in conjunction with the radiative transfer model Satellite Data Simulator Unit (SDSU) that is based on the deltaEddington approximation and includes Mie scattering. The study area is the Icelandic subglacial volcanic region. The analyzed case study is that of the Grímsvötn eruption in May 2011. ATHAM input parameters are adjusted using available ground data, and sensitivity tests are conducted to investigate the observed brightness temperatures and their variance. The tests are based on the variation of environmental conditions like the terrain emissivity, water vapor, and ice in the volcanic plume. Quantitative correlation analysis between ATHAM/SDSU forward-model columnar content simulations and available microwave radiometric brightness temperature measurements, derived from the Special Sensor Microwave Imager/Sounder (SSMIS), are encouraging in terms of both dynamic range and correlation coefficient. The correlation coefficients are found to vary from $-$ 0.37 to $-$ 0.63 for SSMIS channels from 91 to 183 $pm$ 1 GHz, respectively. The larger sensitivity of the brightness temperature at 183 $pm$ 1 GHz to the columnar content, with respect to other channels, allowed us to consider this channel as the basis for a model-based polynomial relationship of volcanic plume hei- ht as a function of the measured SSMIS brightness temperature.

Description: Bright curvilinear features arising from the geometry of man-made structures are characteristic of synthetic aperture radar (SAR) images of urban areas, particularly due to double-reflection mechanisms. An approach to urban earthquake damage detection using double-reflection line amplitude change in single-look images has been established in previous literature. Based on this method, this paper introduces an automated tool for fast, unsupervised damage detection in urban areas. Ridge-based curvilinear features are extracted from a preevent SAR image, and double-reflection candidates are selected using prior probability distributions derived from a simple geometrical building model. The candidate features are then used with the ratio of a pair of single preevent and postevent SAR single-look amplitude images to estimate damage levels. The algorithm is very efficient, with overall computational complexity of $O(Nlog k)$ for an $N$ -pixel image containing features of mean length $k$ . The technique is demonstrated using COSMO-SkyMed data covering L'Aquila, Italy, and Port-au-Prince, Haiti.

Description: It is necessary to measure the sharpness of distributions in many situations. A class of functions is investigated in this paper. First, the relation between this class and sharpness is clarified, and this justifies this class as sharpness measures. Then, we analyze the performance of different sharpness measures and present a guide to select the sharpness measure. In addition, the relation of this class to the sparsity measure is addressed, which leads to a deeper understanding about sparsity. Finally, we show and discuss the application of this class in inverse synthetic aperture radar imaging.

Description: A detection of radio-frequency interference (RFI) in the space-borne microwave radiometer data is difficult under snow and sea ice-covered conditions. The existing methods such as a spectral difference technique or a principal component analysis (PCA) of RFI indices produce many false RFI signals near the boundary of Greenland and Antarctic ice sheets. In this paper, a double PCA (DPCA) method is developed for RFI detection over Greenland and Antarctic regions. It is shown that the new DPCA method is effective in detecting RFI signals in the C- and X-band radiometer channels of WindSat while removing the false RFI signals over Greenland and Antarctic. It also worked well in other snow-free or snow-rich regions such as winter data over the United States. The proposed DPCA can be applied to satellite radiometer data orbit-by-orbit or granule-by-granule and is thus applicable in an operational environment for fast processing and data dissemination.

Description: Advertisement: Now over three million documents in IEEE Xplore. Thank you for your authorship.

Description: A super-resolution (SR) method based on compressive sensing (CS), structural self-similarity (SSSIM), and dictionary learning is proposed for reconstructing remote sensing images. This method aims to identify a dictionary that represents high resolution (HR) image patches in a sparse manner. Extra information from similar structures which often exist in remote sensing images can be introduced into the dictionary, thereby enabling an HR image to be reconstructed using the dictionary in the CS framework. We use the K-Singular Value Decomposition method to obtain the dictionary and the orthogonal matching pursuit method to derive sparse representation coefficients. To evaluate the effectiveness of the proposed method, we also define a new SSSIM index, which reflects the extent of SSSIM in an image. The most significant difference between the proposed method and traditional sample-based SR methods is that the proposed method uses only a low-resolution image and its own interpolated image instead of other HR images in a database. We simulate the degradation mechanism of a uniform 2 $times$ 2 blur kernel plus a downsampling by a factor of 2 in our experiments. Comparative experimental results with several image-quality-assessment indexes show that the proposed method performs better in terms of the SR effectivity and time efficiency. In addition, the SSSIM index is strongly positively correlated with the SR quality.

Description: This paper presents a method and experimental results for near-surface wind sensing using reflected Global Navigation Satellite Systems (GNSS) signals received on a spacecraft. The estimation method proposed involves four steps. First, the bistatic radar cross section (BRCS) of the received signal is estimated from the measurements. Second, the BRCS measurements are calibrated to agree with existing theoretical and empirical wind–wave models. Next, a geometric optics-based scattering model is used to estimate the sea surface slopes, based on the reflection geometry and the measured BRCS. Finally, the surface winds are estimated using an empirically derived function relating the surface mean square slopes to near-surface wind speed. The accuracy of the proposed inversion technique is then tested using a set of 25 space-based GNSS reflection measurements over a range of wind speeds. These measurements were all taken in the proximity of ocean buoys which provided in situ ocean wind speed information. The wind estimates from the buoys were then compared with the wind retrievals made from the measurements and found to be accurate to a root-mean-square error of 1.84 m/s. Additionally, the potential error sources in the measurements are analyzed, including a simulation of the effects of wind direction on the BRCS measurements. This first demonstration of space-based GNSS scatterometry using a small set of sample measurements will hopefully provide a benchmark and example for future experiments.

Description: There is an error in the above-named article [ibid.,vol. 51, no. 4, pp. 2119??2127, Apr. 2013] regarding the definition and the implementation of equation (3), defining the proposed temporal smoothing index (TSI). The correct formula is provided. These corrections do not change any of the general conclusions of the paper, but some of the comments regarding the interpretation of this table are revised.

Description: We present a novel method for ground moving target imaging using a synthetic aperture radar system transmitting ultranarrowband continuous waveforms (CW). Our method exploits the high Doppler resolution provided by ultranarrowband CW signals to image both the scene reflectivity and to determine the velocity of multiple moving targets. We develop a new forward model based on the temporal Doppler induced by the movement of antennas and moving targets. The forward model relates reflectivity and velocity information at each location to a correlated received signal. We form the reflectivity images of the moving targets and estimate their motion parameters using a filtered-backprojection (FBP) technique combined with the contrast or gradient optimization method. The method results in focused reflectivity images of moving targets and their velocity estimates, regardless of the target location, speed, and velocity direction. We show that the amplitude and visible edges of the targets can be correctly reconstructed when the correct target velocity estimate is used in the FBP imaging. We present the resolution analysis of the reflectivity images. Extensive numerical simulations demonstrate the performance of our method and validate the theoretical results.

Description: The author has written this short message to introduce himself as the new Editor-in-Chief (EiC) of the IEEE JOURNAL OF OCEANIC ENGINEERING. He took on the role last autumn, following the untimely death of his predecessor, Bill Carey. The new EiC gives an indication of his priorities and plans; as they become fleshed out, he will provide more details in subsequent editorial messages.

Description: In recent years, in an attempt to maximize performance, machine learning approaches for event-related potential (ERP) spelling have become more and more complex. In this paper, we have taken a step back as we wanted to improve the performance without building an overly complex model, that cannot be used by the community. Our research resulted in a unified probabilistic model for ERP spelling, which is based on only three assumptions and incorporates language information. On top of that, the probabilistic nature of our classifier yields a natural dynamic stopping strategy. Furthermore, our method uses the same parameters across 25 subjects from three different datasets. We show that our classifier, when enhanced with language models and dynamic stopping, improves the spelling speed and accuracy drastically. Additionally, we would like to point out that as our model is entirely probabilistic, it can easily be used as the foundation for complex systems in future work. All our experiments are executed on publicly available datasets to allow for future comparison with similar techniques.

Description: Compressive sensing has shown significant promise in biomedical fields. It reconstructs a signal from sub-Nyquist random linear measurements. Classical methods only exploit the sparsity in one domain. A lot of biomedical signals have additional structures, such as multi-sparsity in different domains, piecewise smoothness, low rank, etc. We propose a framework to exploit all the available structure information. A new convex programming problem is generated with multiple convex structure-inducing constraints and the linear measurement fitting constraint. With additional a priori information for solving the underdetermined system, the signal recovery performance can be improved. In numerical experiments, we compare the proposed method with classical methods. Both simulated data and real-life biomedical data are used. Results show that the newly proposed method achieves better reconstruction accuracy performance in term of both L1 and L2 errors.

Description: Transcranial magnetic stimulation (TMS) is a tool for noninvasive stimulation of neuronal tissue used for research in cognitive neuroscience and to treat neurological disorders. Many TMS applications call for large electric fields to be sharply focused on regions that often lie deep inside the brain. Unfortunately, the fields generated by present-day TMS coils diffuse and decay rapidly as they penetrate into the head. As a result, they tend to stimulate relatively large regions of tissue near the brain surface. Earlier studies suggested that a focused TMS excitation can be attained using multiple nonuniformly fed coils in a multichannel array. We propose a systematic, genetic algorithm-based technique for synthesizing multichannel arrays that minimize the volume of the excited region required to achieve a prescribed penetration depth and maintain realistic values for the input driving currents. Because multichannel arrays are costly to build, we also propose a method to convert the multichannel arrays into single-channel ones while minimally materially deteriorating performance. Numerical results show that the new multi- and single-channel arrays stimulate tissue 2.4 cm into the head while exciting 3.0 and 2.6 times less volume than conventional Figure-8 coils, respectively.

Description: In difficult epileptic patients, the brain structures are explored by means of depth multicontact electrodes [stereoelectroencephalography (SEEG)]. Recently, a novel diagnostic technique allows an accurate definition of the epileptogenic zone using deep brain stimulation (DBS). The stimulation signal propagates in the brain and thus it appears on most of the other SEEG electrodes, masking the local brain electrophysiological activity. The objective of this paper is the DBS-SEEG signals detrending and denoising in order to recover the masked physiological sources. We review the main filtering methods and put forward an approach based on the combination of filtering with generalized eigenvalue decomposition (GEVD). An experimental study on simulated and real SEEG shows that our approach is able to separate DBS sources from brain activity. The best results are obtained by an original singular spectrum analysis-GEVD approach.

Description: This paper presents a novel palpation probe based on optical fiber technology. It is designed to measure stiffness distribution of a soft tissue while sliding over the tissue surface in a near frictionless manner. A novelty of the probe is its ability to measure indentation depth for nonplanar tissue profiles which are commonly experienced during surgery. Since tumors are often harder than the surrounding tissue, the proposed probe can intraoperatively aid the surgeon to rapidly identify the presence, location, and size of the tumors through the generation of a tissue stiffness map. The probe can concurrently measure tissue reaction force, indentation depth, and the orientation of the probe with respect to the tissue surface. Hence, it can generate an elasticity model of the tissue with minimum measurement inaccuracies caused by surface profile variations. Further, the probe has a tunable force range and the indentation force can be adjusted externally to match tissue limitations. The performance of the probe developed was validated using simulated soft tissues samples. Our tumor identification experiments showed that the probe can accurately identify the location and size of tumors hidden inside nonflat tissue surfaces. Further, the probe has clearly demonstrated its potential to identify tumors with tumor–tissue stiffness ratios as low as 2.1.

Description: We introduce a novel method for the automatic detection of atrial fibrillation (AF) using time-varying coherence functions (TVCF). The TVCF is estimated by the multiplication of two time-varying transfer functions (TVTFs). The two TVTFs are obtained using two adjacent data segments with one data segment as the input signal and the other data segment as the output to produce the first TVTF; the second TVTF is produced by reversing the input and output signals. We found that the resultant TVCF between two adjacent normal sinus rhythm (NSR) segments shows high coherence values (near 1) throughout the entire frequency range. However, if either or both segments partially or fully contain AF, the resultant TVCF is significantly lower than 1. When TVCF was combined with Shannon entropy (SE), we obtained even more accurate AF detection rate of 97.9% for the MIT-BIH atrial fibrillation (AF) database $(n$ = 23) with 128 beat segments. The detection algorithm was tested on four databases using 128 beat segments: the MIT-BIH AF database, the MIT-BIH NSR database ( $n$ = 18), the MIT-BIH Arrhythmia database ( $n$ = 48), and a clinical 24-h Holter AF database ( $n$ = 25). Using the receiver operating characteristic curves from the combination of TVCF and SE, we obtained a sensitivity of 98.2% and specificity of 97.7% for the MIT-BIH AF database. For the MIT-BIH NSR database, we found a specificity of 99.7%. For the MIT-BIH Arrhythmia database, the sensitivity and specificity were 91.1% and 89.7%, respectively. For the clinical database (24-h Holter data), the sensitivity and specificity were 92.3% and 93.6%, respectively. We - lso found that a short segment (12 beats) also provided accurate AF detection for all databases: sensitivity of 94.7% and specificity of 90.4% for the MIT-BIH AF, specificity of 94.4% for the MIT-BIH-NSR, the sensitivity of 92.4% and specificity of 84.1% for the MIT-BIH arrhythmia, and sensitivity of 93.9% and specificity of 84.4% for the clinical database. The advantage of using a short segment is more accurate AF burden calculation as the timing of transitions between NSR and AF are more accurately detected.

Description: In functional MRI, it is often desirable to reduce the readout duration to make the acquired data less prone to $T_2^*$ susceptibility artifacts. In addition, a shorter readout length allows for a shorter minimum TE, which is important for optimizing SNR. This can be achieved by undersampling the $k$ -space. However, the conventional Fourier transform-based reconstruction method suffers from under-sampling artifacts such as high-frequency ringing and loss of resolution. To address this problem, we revisit the constrained-model approach using the generalized-series (GS) which has been proposed to address the undersampling problem for dynamic MRI. We propose a modification to the conventional use of the model in order to reflect small hemodynamic signal changes typical in fMRI. Specifically, while realizing that having high model order is necessary to capture missing information, we found that it is not necessary to span all frequencies of GS basis functions uniformly. Instead, having $k$ -space and GS “sampling” trajectories covering low-frequencies uniformly while spanning high-frequencies sparsely, was observed to be an efficient strategy. The ability of the method over the conventional GS approach in improving resolution of functional images and activation maps while reducing undersampling ringing is demonstrated by simulations and experiments at $3$ T. Reduction in the readout time allowed an increase of statistical signal power as compared to the fully sampled acquisition. Unlike compressed sensing approaches, the proposed method is linear and hence has lower computational complexity. The method could prove useful for other imaging modalities where the signal change is smaller than the ba- eline component.

Description: Interbreath interval (IBI), the time interval between breaths, is an important measure used to analyze irregular breathing patterns in neonates. The discrete bursts of neural activity generate the IBI time series, which exhibits stochastic as well as deterministic dynamics. To quantify the irregularity of breathing, we propose a point process model of IBI using a comprehensive stochastic dynamic modeling framework. The IBIs of immature breathing patterns exhibit a long tail distribution and within a point process model, we have considered the lognormal distribution to represent the stochastic IBI characteristics. An autoregressive (AR) function is embedded within the model to capture the short-term IBI dynamics including abrupt IBI prolongations related to sporadic and periodic apneas that are common in neonates. We tested the utility of our paradigm for depicting the respiratory dynamics in neonatal rats and in preterm infants. Kolmogorov–Smirnov (KS) and independence tests reveal that the model accurately tracks the dynamic characteristics of the signals. In preterm infants, our model-derived indices of IBI instability strongly correlate with clinically derived indices of maturation. Our results validate a new class of algorithms, based on the point process theory, for defining instantaneous measures of breathing irregularity in neonates.

Description: Remote photoplethysmography (rPPG) enables contactless monitoring of the blood volume pulse using a regular camera. Recent research focused on improved motion robustness, but the proposed blind source separation techniques (BSS) in RGB color space show limited success. We present an analysis of the motion problem, from which far superior chrominance-based methods emerge. For a population of 117 stationary subjects, we show our methods to perform in 92% good agreement ( $pm hbox{1.96}sigma$ ) with contact PPG, with RMSE and standard deviation both a factor of 2 better than BSS-based methods. In a fitness setting using a simple spectral peak detector, the obtained pulse-rate for modest motion (bike) improves from 79% to 98% correct, and for vigorous motion (stepping) from less than 11% to more than 48% correct. We expect the greatly improved robustness to considerably widen the application scope of the technology.

Description: A system for capturing habitual, in-home gait measurements using an environmentally mounted depth camera, the Microsoft Kinect, is presented. Previous work evaluating the use of the Kinect sensor for in-home gait measurement in a lab setting has shown the potential of this approach. In this paper, a single Kinect sensor and computer were deployed in the apartments of older adults in an independent living facility for the purpose of continuous, in-home gait measurement. In addition, a monthly fall risk assessment protocol was conducted for each resident by a clinician, which included traditional tools such as the timed up a go and habitual gait speed tests. A probabilistic methodology for generating automated gait estimates over time for the residents of the apartments from the Kinect data is described, along with results from the apartments as compared to two of the traditionally measured fall risk assessment tools. Potential applications and future work are discussed.

Description: Visualization of ex vivo human patellar cartilage matrix through the phase contrast imaging X-ray computed tomography (PCI-CT) has been previously demonstrated. Such studies revealed osteoarthritis-induced changes to chondrocyte organization in the radial zone. This study investigates the application of texture analysis to characterizing such chondrocyte patterns in the presence and absence of osteoarthritic damage. Texture features derived from Minkowski functionals (MF) and gray-level co-occurrence matrices (GLCM) were extracted from 842 regions of interest (ROI) annotated on PCI-CT images of ex vivo human patellar cartilage specimens. These texture features were subsequently used in a machine learning task with support vector regression to classify ROIs as healthy or osteoarthritic; classification performance was evaluated using the area under the receiver operating characteristic curve (AUC). The best classification performance was observed with the MF features perimeter (AUC: $0.94 pm 0.08$ ) and “Euler characteristic” (AUC: $0.94 pm 0.07$ ), and GLCM-derived feature “Correlation” (AUC: $0.93 pm 0.07$ ). These results suggest that such texture features can provide a detailed characterization of the chondrocyte organization in the cartilage matrix, enabling classification of cartilage as healthy or osteoarthritic with high accuracy.

Description: We have developed a pupil-corneal reflection method-based gaze detection system, which allows large head movements and achieves easy gaze calibration. This system contains two optical systems consisting of components such as a camera and a near-infrared light source attached to the camera. The light source has two concentric LED rings with different wavelengths. The inner and outer rings generate bright and dark pupil images, respectively. The pupils are detected from a difference image created by subtracting the bright and dark pupil images. The light source also generates the corneal reflection. The 3-D coordinates of the pupils are determined by the stereo matching method using two optical systems. The vector from the corneal reflection center to the pupil center in the camera image is determined as ${bm r}$ . The angle between the line of sight and the line passing through the pupil center and the camera (light source) is denoted as $theta$ . The relationship $theta = k leftvert {bm r} rightvert$ is assumed, where $k$ is a constant. The theory implies that head movement of the user is allowed and facilitates the gaze calibration procedure. In the automatic calibration method, $k$ is automatically determined while the user looks around on the PC screen without fixating on any specific calibration target. In the one-point calibration method, the user is asked to fixate on one calibration target at the PC screen in order to correct the difference between the optical and visual axes. In the two-point calibration method, in order to correct the nonlinear relationship between $theta$ and $leftvert {bm r} rightvert$ , the user is asked to fixate on two targets. The experimental results show that the three proposed calibration methods improve the precision of gaze detection step by step. In addition, the average gaze error in the visual angle is less than 1 $^circ$ for the seven head positions of the user.

Description: We present an encapsulation scheme that combines atomic layer deposited (ALD) Al $_{2}$ O $_{3}$ and Parylene C for the encapsulation of implantable devices. The encapsulation performances of combining alumina and Parylene C was compared to individual layers of Parylene C or alumina and the bilayer coating had superior encapsulation properties. The alumina–Parylene coated interdigitated electrodes (IDEs) soaked in PBS for up to nine months at temperatures from 37 to 80 °C for accelerated lifetime testing. For 52-nm alumina and 6-μm Parylene C, leakage current was ∼20 pA at 5 VDC, and the impedance was about 3.5 MΩ at 1 kHz with a phase near −87° from electrochemical impedance spectroscopy for samples soaked at 67 °C for equivalent lifetime of 72 months at 37 °C. The change of impedance during the whole soaking period (up to 70 months of equivalent soaking time at 37 °C) over 1 to 10 $^{6}$  Hz was within 5%. The stability of impedance indicated almost no degradation of the encapsulation. Bias voltage effect was studied by continuously applying 5 VDC, and it reduced the lifetime of Parylene coating by ∼75% while it showed no measurable effect on the bilayer coating. Lifetime of encapsulation of IDEs with topography generated by attaching a coil and surface mount device (SMD) capacitor was about half of that of planer IDEs. The stable long-term insulation impedance, low leakage current, and better lifetime under bias voltage and topography made this double-layer encapsulation very promising for chronic implantable devices.

Description: Muscle sympathetic nerve activity is a primary source of cardiovascular control in humans. Traditional analyses smooth away the fine temporal structure of the sympathetic recordings, limiting our understanding of sympathetic activation mechanisms. We use multifiber spike trains extracted from standard microneurography voltage trace to characterize the sympathetic spiking at rest and during sympathoexcitation. Our analysis corroborates known features of sympathetic activity, such as bursting behavior, cardiac rhythmicity, and long conduction delays. It also elucidates new features such as large heartbeat-to-heartbeat variability of firing rates and precise pattern of spiking within cardiac cycles. We find that at low firing rates, spikes occur uniformly throughout the cardiac cycle, but at higher rates, they tend to cluster in bursts around a particular latency. This latency shortens and the clusters tighten as the firing rates grow. Sympathoexcitation increases firing rates and shifts the burst latency later. Negative rate/latency correlation and the sympathoexcitatory shift suggest that spike production of the individual fibers contributes significantly to the control of the sympathetic bursts strength. Access to fine scale temporal information, more physiologically accurate description of nerve activity, and new hypotheses about the nervous outflow control establishes sympathetic spiking as a valuable tool for the cardiovascular research.

Description: In this paper, we report the development of a flexible base array of penetrating electrodes which can be used to interface with the spinal cord. A customizable and feasible fabrication protocol is described. The flexible base arrays were fabricated and implanted into surrogate cords which were elongated by 12%. The resulting strains were optically measured across the cord and compared to those associated with two types of electrodes arrays (one without a base and one with a rigid base connecting the electrodes). The deformation behavior of cords implanted with the flexible base arrays resembled the behavior of cords implanted with individual microwires that were not connected through a base. The results of the strain test were used to validate a 2-D finite element model. The validated model was used to assess the stresses induced by the electrodes of the three types of arrays on the cord, and to examine how various design parameters (thickness, base modulus, etc.,) impact the mechanical behavior of the electrode array. Rigid base arrays induced higher stresses on the cord than the flexible base arrays which in turn imposed higher stresses than the individual microwire implants. The developed flexible base array showed improvement over the rigid base array; however, its stiffness needs to be further reduced to emulate the mechanical behavior of individual microwire arrays without a base.

Description: This paper presents a novel energy management strategy for a hybridized power source small urban electric vehicle (EV). First, an analysis of load requirements for typical urban driving cycles is presented. Thereafter, the energy and power management issues are addressed for a multisource EV to define an improved management architecture. A dynamically restricted search space strategy coupled with a simulated annealing technique is exploited to accomplish the global optimization of the energy management system (EMS). The control of the dc/dc converter operations based on this EMS is also presented. The multiple sources have been simulated using an overall model implemented in Matlab/Simulink. A reduced-scale prototype has been built to analyze the performance of the energy management strategy. The results obtained show that energy management has been enhanced leading to an increase of the vehicle performance with reduced size embarked sources.

Description: Prediction intervals (PIs) are a promising tool for quantification of uncertainties associated with point forecasts of wind power. However, construction of PIs using parametric methods is questionable, as forecast errors do not follow a standard distribution. This paper proposes a nonparametric method for construction of reliable PIs for neural network (NN) forecasts. A lower upper bound estimation (LUBE) method is adapted for construction of PIs for wind power generation. A new framework is proposed for synthesizing PIs generated using an ensemble of NN models in the LUBE method. This is done to guard against NN performance instability in generating reliable and informative PIs. A validation set is applied for short listing NNs based on the quality of PIs. Then, PIs constructed using filtered NNs are aggregated to obtain combined PIs. Performance of the proposed method is examined using data sets taken from two wind farms in Australia. Simulation results indicate that the quality of combined PIs is significantly superior to the quality of PIs constructed using NN models ranked and filtered by the validation set.

Description: This paper presents a supervisory control unit (SCU) combined with short-term ahead wind speed prediction for proper and effective management of the stored energy in a small capacity flywheel energy storage system (FESS) which is used to mitigate the output power fluctuations of an aggregated wind farm. Wind speed prediction is critical for a wind energy conversion system since it may greatly influence the issues related to effective energy management, dynamic control of wind turbine, and improvement of the overall efficiency of the power generation system. In this study, a wind speed prediction model is developed by artificial neural network (ANN) which has advantages over the conventional prediction schemes including data error tolerance and ease in adaptability. The proposed SCU-based control would help to reduce the size of the energy storage system for minimizing wind power fluctuation taking the advantage of prediction scheme. The model for prediction using ANN is developed in MATLAB/Simulink and interfaced with PSCAD/EMTDC. Effectiveness of the proposed control system is illustrated using real wind speed data in various operating conditions.

Description: Due to its reduced communication overhead and robustness to failures, distributed energy management is of paramount importance in smart grids, especially in microgrids, which feature distributed generation (DG) and distributed storage (DS). Distributed economic dispatch for a microgrid with high renewable energy penetration and demand-side management operating in grid-connected mode is considered in this paper. To address the intrinsically stochastic availability of renewable energy sources (RES), a novel power scheduling approach is introduced. The approach involves the actual renewable energy as well as the energy traded with the main grid, so that the supply–demand balance is maintained. The optimal scheduling strategy minimizes the microgrid net cost, which includes DG and DS costs, utility of dispatchable loads, and worst-case transaction cost stemming from the uncertainty in RES. Leveraging the dual decomposition, the optimization problem formulated is solved in a distributed fashion by the local controllers of DG, DS, and dispatchable loads. Numerical results are reported to corroborate the effectiveness of the novel approach.

Description: The present paper investigates the effect of high penetration of photovoltaic (PV) systems on the small signal stability of a large power system. Reduced system inertia and altered power flow patterns as a result of the addition of the utility scale and residential rooftop PVs that replace a portion of conventional generation resources, may lead to decreased damping of the critical modes of the system. To identify the critical modes of the system and the effect of the high PV penetration on these modes, eigenvalue analysis is carried out on the aforementioned system under various PV penetration levels. To substantiate the results observed from the small signal analysis, transient analysis is carried out on the system under various PV penetration levels. The simulation results effectively identify the impact of high PV penetration on small signal stability of the studied system.

Description: A speed exclusion zone algorithm to prevent the rotor rotational frequency from exciting the side-side tower mode in a variable-speed wind turbine was designed. The algorithm was first tested on a simple wind turbine model developed in Simulink and then implemented and tested on an experimental test rig. Different speed exclusion zone widths were tested to observe how the performance of the wind turbine is affected. Results obtained from both simulations and experiments showed the effectiveness of the speed exclusion zone algorithm in preventing excitation of tower resonance. Results from different speed exclusion zone widths showed that with a wider exclusion zone, the tower mode is excited less at the expense of a larger power loss.

Description: A control strategy for the quasi-Z-source inverter (qZSI) with a battery-based photovoltaic (PV) power conversion system is proposed. A battery-assisted qZSI can buck/boost PV panel voltage by introducing shoot-through states, and make full use of PV power by the energy-stored battery paralleled to the quasi-Z-source capacitor. A dynamic small-signal model of the battery-assisted qZSI is established to design a closed-loop controller for regulating shoot-through duty ratio and managing the battery's energy storage. A modified space vector modulation (SVM) technique for the qZSI is applied to achieve low harmonics, high voltage utilization, and high efficiency. A P-Q decoupled grid-tie power injection is fulfilled with the maximum power capture from PV panels and the unity power factor. The validity of the proposed PV system is proved by experimental results, showing an efficient method for the energy-stored PV power generation.

Description: Doubly-fed induction generator (DFIG)-based wind turbines (WTs) are reported to suffer reliability problems due to the presence of slip-rings, brushes, and the gearbox. These disadvantages encouraged several research groups to investigate the viability of employing single or double-frame brushless cascade DFIGs (BCDFIGs) in grid-connected multimegawatt WTs especially offshore and in solar chimney power plants. In this regard, this paper tackles three main issues. Initially, the expected reduction in gear ratio when BCDFIGs are used instead of DFIGs is examined from a steady state perspective. A reduction would lead to less frequent maintenance and an improved return on investment. Next, a detailed comparison between DFIG-WTs and BCDFIG-WTs under unbalanced grid voltage is presented. Finally, the extent to which a multimegawatt BCDFIG-WT is grid code compliant in terms of fault ride-through capability is studied. This is illustrated by comparing and quantifying the response of a DFIG-WT and a BCDFIG-WT to a severe three-phase voltage dip. Simulation and experimental results indicate promising behavior for BCDFIGs during disturbances.

Description: We analyze the potential benefits of colocating wind and concentrating solar power (CSP) plants in the southwestern U.S. Using a location in western Texas as a case study, we demonstrate that such a deployment strategy can improve the capacity factor of the combined plant and the associated transmission investment. This is because of two synergies between wind and CSP: 1) the negative correlation between real-time wind and solar resource availability and 2) the use of low-cost high-efficiency thermal energy storage in CSP. The economic tradeoff between transmission and system performance is highly sensitive to CSP and transmission costs. We demonstrate that a number of deployment configurations, which include up to 67% CSP, yield a positive net return on investment.

Description: Renewable energy resources, especially wind power, are expected to provide a considerable portion of the world energy requirements in the near future. Large-scale wind power penetration impacts the electricity industry in many aspects and raises a number of technical challenges for the electricity network. A day-ahead network-constrained market clearing formulation is proposed which considers demand side resources. The proposed approach can provide flexible load profile and reduce the need for ramp up/down services by the conventional generators. This method can potentially facilitate a large penetration of wind power by shifting the wind power generation from the off-peak periods to the high-peak hours. The validity of the proposed approach has been verified using the IEEE 30 bus and 57 bus test systems.

Description: The small generation units including cogeneration units and renewables have been widely spread during the last decades. These units raise several questions regarding the power system, e.g., their integration into power system control or their role on the supply-side of smart grids. A multiagent model was established in order to facilitate the investigations of the emerged questions. Several viewpoints are taken into account and built into the model considering technical specifications, proper wind-speed and temperature simulations, constraints of heating service, fuel consumption, legal regulation, outages, and services. A state-based method is elaborated in order to generalize different technologies. The different aspects of operation control are represented by agent strategies. Applying these methods resulted in an easy-to-follow and flexible model. This paper is dealing with the developed multiagent-based model applied to study the distributed generation system from many viewpoints like power system balance, virtual power plant, or smart grid scheduling.

Description: The electromagnetic spectrum is a valued shared resource. Its scientific use allows us to learn about our universe, measure and monitor our planet, and communicate scientific data. The use of the spectrum is managed by national, regional, and global regulatory frameworks. There are increasing demands for new or extended allocations because of vast technological advances in the past few years. Understanding spectrum management is important in the successful planning and execution of missions and instruments, as well as in determining the potential source of radio frequency interference in existing data and instruments, and in working to ameliorate its impact. This paper provides a summary of this framework for radio scientists and engineers.

Description: For highly squinted synthetic aperture radar (SAR) imaging, the wavenumber domain SAR processing algorithm is commonly accepted as an ideal solution to SAR focusing in the case of an ideal straight sensor trajectory. However, airborne SAR is very sensitive to atmospheric turbulence that causes serious trajectory deviations. In this paper, we propose a robust autofocusing approach for highly squinted airborne SAR imagery using the extended wavenumber algorithm, being capable of estimating the range-dependent phase errors. To apply the proposed autofocusing scheme, a detailed analysis of the motion error model in the conical reference system is presented, where the formulation of range-dependent phase errors for squinted SAR is given. The proposed autofocusing approach is performed by a three-step process: 1) referring to the inevitable residual phase after deramping for highly squinted SAR, a modified squinted phase gradient autofocusing (SPGA) algorithm is put forward to retrieve the range-independent phase errors; 2) based on the established motion error model, the residual range-dependent phase errors are estimated using a local maximum likelihood-weighted SPGA algorithm; and 3) motion compensation is executed by a two-step approach to reach the range-independent and range-dependent corrections, respectively. Experiments based on measured data have shown that the proposed autofocusing approach performs well for highly squinted SAR imaging.

Description: Hyperspectral radiances from the Infrared Atmospheric Sounding Interferometer (IASI) and Atmospheric Infrared Sounder (AIRS) are used as a reference to improve the calibration accuracy for FengYun-2 (FY-2) infrared (IR) channel radiances. It is shown that the previous FY-2 operational calibration for IR bands produces significant bias in brightness temperatures that can exceed 1.1 K. In particular, the FY-2 IR3 band (6.7 $muhbox{m}$ ) has the largest bias of 2.0 K. The daytime double-difference temperature (DDT) between AIRS and IASI using FY-2 imagers as a transfer medium showed an excellent consistency, is within 0.2 K at 290 K, and is stable over time for FY-2C/2D/2E. This only indicates the robust calibrations applied for both the AIRS and IASI measurements. During the nighttime of the Earth observation, stray light in space affects the long-term stability of the FY-2 DDT, particularly for the Earth scene at 220 K. FY-2E satellite which was launched in 2009 has an instrument design improvement. Intercalibrating FY-2 four times using AIRS and IASI data can reveal the diurnal features of the FY-2 instrument calibration. The temporal DDT appears very large during the spring and autumn eclipse times. Not only can the global-space-based-intercalibration-system intercalibration method provide an excellent operational calibration for the FY-2 imager, but it can also help improve the design of future instruments and onboard blackbody calibration.

Description: Advertisement: Internet television gets a mobile makeover. A mobile version of IEEE.tv is now available for convenient viewing. Plus a new app for IEEE.tv can also be found in your app store.

Description: Compared with traditional remote sensing, multiangular observation provides 3-D structural information of a forest through different directional observations. The MGeoSAIL model, suitable for multiangular observations, was developed based on the single-angle model GeoSAIL. The MGeoSAIL model combines the geometric-optic model with the radiation transfer model and has the advantages of both models. Thus, it is more accurate and feasible. The geometric-optic model calculates the amount of shadowed and illuminated components within a forest scene, while the radiation transfer model [Scattering by Arbitrarily Inclined Leaves (SAIL)] calculates the reflectance and transmittance of tree crowns. The uniform index is introduced to characterize the relationship quantitatively between tree distribution pattern and the bidirectional reflectance distribution function (BRDF). The simulation results show that the MGeoSAIL model could simulate the “hot” spot in red and near-infrared bands, as well as the “bowl” shape in the near-infrared band. The relationship between the uniform index and BRDF is negatively exponential. Finally, the look-up table was calculated using the MGeoSAIL model, and leaf area index (LAI) was inversed from compact high-resolution imaging spectrometry data. The results compared well with the measured LAI in Changbai Mountain area, China.

Description: Provides a listing of current staff, committee members and society officers.

Description: The paper presents a new endoscope motion tracking method that is based on a novel external endoscope tracking device and our modified stochastic optimization method for boosting endoscopy navigation. We designed a novel tracking prototype where a 2-D motion sensor was introduced to directly measure the insertion–retreat linear motion and also the rotation of the endoscope. With our developed stochastic optimization method, which embeds traceable particle swarm optimization in the Condensation algorithm, a full six degrees-of-freedom endoscope pose (position and orientation) can be recovered from 2-D motion sensor measurements. Experiments were performed on a dynamic bronchial phantom with maximal simulated respiratory motion around 24.0 mm. The experimental results demonstrate that our proposed method provides a promising endoscope motion tracking approach with more effective and robust performance than several current available tracking techniques. The average tracking accuracy of the position improved from 6.5 to 3.3 mm, which further approaches the clinical requirement of 2.0 mm in practice.

Description: Radical prostatectomy is performed on approximately 40% of men with organ-confined prostate cancer. Pathologic information obtained from the prostatectomy specimen provides important prognostic information and guides recommendations for adjuvant treatment. The current pathology protocol in most centers involves primarily qualitative assessment. In this paper, we describe and evaluate our system for automatic prostate cancer detection and grading on hematoxylin & eosin-stained tissue images. Our approach is intended to address the dual challenges of large data size and the need for high-level tissue information about the locations and grades of tumors. Our system uses two stages of AdaBoost-based classification. The first provides high-level tissue component labeling of a superpixel image partitioning. The second uses the tissue component labeling to provide a classification of cancer versus noncancer, and low-grade versus high-grade cancer. We evaluated our system using 991 sub-images extracted from digital pathology images of 50 whole-mount tissue sections from 15 prostatectomy patients. We measured accuracies of 90% and 85% for the cancer versus noncancer and high-grade versus low-grade classification tasks, respectively. This system represents a first step toward automated cancer quantification on prostate digital histopathology imaging, which could pave the way for more accurately informed postprostatectomy patient care.

Description: We present DevStaR, an automated computer vision and machine learning system that provides rapid, accurate, and quantitative measurements of C. elegans embryonic viability in high-throughput (HTP) applications. A leading genetic model organism for the study of animal development and behavior, C. elegans is particularly amenable to HTP functional genomic analysis due to its small size and ease of cultivation, but the lack of efficient and quantitative methods to score phenotypes has become a major bottleneck. DevStaR addresses this challenge using a novel hierarchical object recognition machine that rapidly segments, classifies, and counts animals at each developmental stage in images of mixed-stage populations of C. elegans. Here, we describe the algorithmic design of the DevStaR system and demonstrate its performance in scoring image data acquired in HTP screens.

Description: Atrial fibrillation is an increasingly prevalent cardiovascular disease; changes in atrial structure and function induced by atrial fibrillation and its treatments are often spatially heterogeneous. However, spatial heterogeneity of function is difficult to assess with standard imaging techniques. This paper describes a method to assess global and regional mechanical function by combining cardiac magnetic resonance imaging and finite-element surface fitting. We used this fitted surface to derive measures of left atrial volume, regional motion, and spatial heterogeneity of motion in 23 subjects, including healthy volunteers and atrial fibrillation patients. We fit the surfaces using a Newton optimization scheme in under 1 min on a standard laptop, with a root mean square error of $2.3pm 0.5$ mm, less than 9% of the mean fitted radius, and an inter-operator variability of less than 10%. Fitted surfaces showed clear definition of the phases of left atrial motion (filling, passive emptying, active contraction) in both volume-time and regional radius-time curves. Averaged surfaces of healthy volunteers and atrial fibrillation patients provided evidence of substantial regional variation in both amount and timing of regional motion, indicating spatial heterogeneity of function, even in healthy adults.

Description: Magnetic resonance spectroscopic imaging (MRSI) is a powerful tool capable of providing spatially localized maps of metabolite concentrations. Its utility, however, is often depreciated by spectral leakage artifacts resulting from low spatial resolution measurements through an effort to reduce acquisition times. Though model-based techniques can help circumvent these drawbacks, they require strong prior knowledge, and can introduce additional artifacts when the underlying models are inaccurate. We introduce a novel scheme in which a generative model is estimated from the raw MRSI data via a regularized variational framework that minimizes the model approximation error within a measurement-prescribed subspace. As additional a priori information, our approach relies only upon a measured field inhomogeneity map at high spatial resolution. We demonstrate the feasibility of our approach on both synthetic and experimental data.

Description: The absorbed energy distribution (AED) in X-ray imaging detectors is an important factor that affects both energy resolution and image quality through the Swank factor and detective quantum efficiency. In the diagnostic energy range (20–140 keV), escape of characteristic photons following photoelectric absorption and Compton scatter photons are primary sources of absorbed-energy dispersion in X-ray detectors. In this paper, we describe the development of an analytic model of the AED in compound X-ray detector materials, based on the cascaded-systems approach, that includes the effects of escape and reabsorption of characteristic and Compton-scatter photons. We derive analytic expressions for both semi-infinite slab and pixel geometries and validate our approach by Monte Carlo simulations. The analytic model provides the energy-dependent X-ray response function of arbitrary compound materials without time-consuming Monte Carlo simulations. We believe this model will be useful for correcting spectral distortion artifacts commonly observed in photon-counting applications and optimal design and development of novel X-ray detectors.

Description: We investigate an approach to evaluation of emission-tomography (ET) imaging systems used for region-of-interest (ROI) estimation tasks. In the evaluation we employ the concept of “emission counts” (EC), which are the number of events per voxel emitted during a scan. We use the reduction in posterior variance of ROI EC, compared to the prior ROI EC variance, as the metric of primary interest, which we call the “posterior variance reduction index” (PVRI). Systems that achieve a higher PVRI are considered superior to systems with lower PVRI. The approach is independent of the reconstruction method and is applicable to all photon-limited data types including list-mode data. We analyzed this approach using a model of 2-D tomography, and compared our results to the classical theory of tomographic sampling. We found that performance evaluations using the PVRI index were consistent with the classical theory. System evaluation based on EC posterior variance is an intuitively appealing and physically meaningful method that is useful for evaluation of system performance in ROI quantitation tasks.

Description: This paper presents a novel probabilistic voxel selection strategy for medical image registration in time-sensitive contexts, where the goal is aggressive voxel sampling (e.g., using less than 1% of the total number) while maintaining registration accuracy and low failure rate. We develop a Bayesian framework whereby, first, a voxel sampling probability field (VSPF) is built based on the uncertainty on the transformation parameters. We then describe a practical, multi-scale registration algorithm, where, at each optimization iteration, different voxel subsets are sampled based on the VSPF. The approach maximizes accuracy without committing to a particular fixed subset of voxels. The probabilistic sampling scheme developed is shown to manage the tradeoff between the robustness of traditional random voxel selection (by permitting more exploration) and the accuracy of fixed voxel selection (by permitting a greater proportion of informative voxels).

Description: Segmentation of the spinal column from computed tomography (CT) images is a preprocessing step for a range of image-guided interventions. One intervention that would benefit from accurate segmentation is spinal needle injection. Previous spinal segmentation techniques have primarily focused on identification and separate segmentation of each vertebra. Recently, statistical multi-object shape models have been introduced to extract common statistical characteristics between several anatomies. These models can be used for segmentation purposes because they are robust, accurate, and computationally tractable. In this paper, we develop a statistical multi-vertebrae shape+pose model and propose a novel registration-based technique to segment the CT images of spine. The multi-vertebrae statistical model captures the variations in shape and pose simultaneously, which reduces the number of registration parameters. We validate our technique in terms of accuracy and robustness of multi-vertebrae segmentation of CT images acquired from lumbar vertebrae of 32 subjects. The mean error of the proposed technique is below 2 mm, which is sufficient for many spinal needle injection procedures, such as facet joint injections.

Description: In this paper, we report on the first edition of the HEp-2 Cells Classification contest, held at the 2012 edition of the International Conference on Pattern Recognition, and focused on indirect immunofluorescence (IIF) image analysis. The IIF methodology is used to detect autoimmune diseases by searching for antibodies in the patient serum but, unfortunately, it is still a subjective method that depends too heavily on the experience and expertise of the physician. This has been the motivation behind the recent initial developments of computer aided diagnosis systems in this field. The contest aimed to bring together researchers interested in the performance evaluation of algorithms for IIF image analysis: 28 different recognition systems able to automatically recognize the staining pattern of cells within IIF images were tested on the same undisclosed dataset. In particular, the dataset takes into account the six staining patterns that occur most frequently in the daily diagnostic practice: centromere, nucleolar, homogeneous, fine speckled, coarse speckled, and cytoplasmic. In the paper, we briefly describe all the submitted methods, analyze the obtained results, and discuss the design choices conditioning the performance of each method.

Description: Tissue displacements required for mechanical property reconstruction in magnetic resonance elastography (MRE) are acquired in a magnetic resonance imaging (MRI) scanner, therefore, anatomical information is available from other imaging sequences. Despite its availability, few attempts to incorporate prior spatial information in the MRE reconstruction process have been reported. This paper implements and evaluates soft prior regularization (SPR), through which homogeneity in predefined spatial regions is enforced by a penalty term in a nonlinear inversion strategy. Phantom experiments and simulations show that when predefined regions are spatially accurate, recovered property values are stable for SPR weighting factors spanning several orders of magnitude, whereas inaccurate segmentation results in bias in the reconstructed properties that can be mitigated through proper choice of regularization weighting. The method was evaluated in vivo by estimating viscoelastic mechanical properties of frontal lobe gray and white matter for five repeated scans of a healthy volunteer. Segmentations of each tissue type were generated using automated software, and statistically significant differences between frontal lobe gray and white matter were found for both the storage modulus and loss modulus $(p 〈 10^{-6})$ . Provided homogeneous property assumptions are reasonable, SPR produces accurate quantitative property estimates for tissue structures which are finer than the resolution currently achievable with fully distributed MRE.

Description: Recent research has demonstrated that improved image segmentation can be achieved by multiple template fusion utilizing both label and intensity information. However, intensity weighted fusion approaches use local intensity similarity as a surrogate measure of local template quality for predicting target segmentation and do not seek to characterize template performance. This limits both the usefulness and accuracy of these techniques. Our work here was motivated by the observation that the local intensity similarity is a poor surrogate measure for direct comparison of the template image with the true image target segmentation. Although the true image target segmentation is not available, a high quality estimate can be inferred, and this in turn allows a principled estimate to be made of the local quality of each template at contributing to the target segmentation. We developed a fusion algorithm that uses probabilistic segmentations of the target image to simultaneously infer a reference standard segmentation of the target image and the local quality of each probabilistic segmentation. The concept of comparing templates to a hidden reference standard segmentation enables accurate assessments of the contribution of each template to inferring the target image segmentation to be made, and in practice leads to excellent target image segmentation. We have used the new algorithm for the multiple-template-based segmentation and parcellation of magnetic resonance images of the brain. Intensity and label map images of each one of the aligned templates are used to train a local Gaussian mixture model based classifier. Then, each classifier is used to compute the probabilistic segmentations of the target image. Finally, the generated probabilistic segmentations are fused together using the new fusion algorithm to obtain the segmentation of the target image. We evaluated our method in comparison to other state-of-the-art segmentation methods. We demonstrated that our new fusio- algorithm has higher segmentation performance than these methods.

Description: Computing Now House Advertisement

Description: The dc component of biological sources is difficult to measure with body surface electrodes. This is due to the existence of an electrical polarization in the interface between the electrode and the electrolyte. A system to acquire a dc component of the biological origin is proposed by means of differentiating the polarization potential from the biological one. The acquisition system is comprised of an instrumentation amplifier, an integrated semiconductor analog switch, and ordinary silver/silver-chloride body surface electrodes. Biopotentials, electrooculograms from the eyes, were noninvasively recorded by the system with the analog switch which was inserted between the electrodes to make the two inputs electrically open following them being shorted. Asymptotes of a triple exponential function were subsequently extracted by regression model analysis on the output signals that were acquired immediately after the switch opened. One of the distinctive asymptotes of the function showed that the magnitude which corresponded to eye movements stayed stable as long as the eyes were fixed. Analysis of the time constants involved in the acquisition system elucidated that the asymptote was the biological origin.

Description: The trend towards multimegawatt (multi-MW) wind turbines and the increasing interest in direct-drive variable-speed wind energy systems have made multilevel converters a promising candidate for large wind energy conversion systems. This paper presents a new hybrid modular multilevel converter for interfacing a full-scale, permanent magnet synchronous generator (PMSG)-based direct-drive variable-speed wind energy conversion system (WECS). The proposed hybrid converter, which is used on the grid side of the system, consists of a three-level modular multilevel converter (MMC) in series connection with three H-bridge modules. The generator-side converter is based on a conventional three-level neutral-point-clamped converter. The proposed hybrid converter, as opposed to the existing full-scale multilevel converter-based wind energy systems, provides structural modularity and a higher dc-bus voltage utilization. This paper formulates and analyzes the internal dynamics of the proposed hybrid converter including the circulating currents and the capacitor voltage ripples. The ac components of the circulating currents, if not properly reduced, increase the amplitude of the capacitor voltage ripples, rating values of the converter components, and losses. Based on the analysis, closed-loop circulating current and capacitor voltage ripple reduction techniques are developed. The reduction of capacitor voltage ripples help reduce the capacitor value. A mathematical model is also developed for the overall WECS. Performance of the overall WECS, under the proposed multilevel converter-based topology and controls, is evaluated based on time domain simulations in the PSCAD/EMTDC environment.

Description: Electric vehicles (EVs) are likely to have a continued presence in the light-vehicle market in the next few decades. As a result, EV charging will put an extra burden on the distribution grid and adjustments need to be made in some cases. On the other hand, EVs have the potential to support the grid as well. This paper presents a single-phase bidirectional charger topology which pairs up a photovoltaic (PV) source with an EV charger resulting in production cost reduction. The presented topology is then used for vehicle-to-grid (V2G) services. The main focus of this paper is on power quality services which only slightly discharge the battery. Among these services, it studies the possibility of local reactive injection of EVs connected to the grid through a single-phase charger to compensate for voltage drops caused by motor startup or inductive loads. It also studies the possibility of active power injection of EVs for short time periods during PV transients in cloudy weather to keep the system stable. It also studies the potential of EVs to help during low voltage ride-through of the PV sources. The studies are performed using Simulink simulations and a real-time implementation in Real Time Digital Simulator (RTDS). The results demonstrate the effectiveness of power quality V2G services with small wear on the EV battery.

Description: The short-term wind power scenarios have a significant impact on the operation cost and power system reliability due to the stochastic generation scheduling of wind-integrated power systems. In order to obtain the scenarios containing the information of forecast error distribution and fluctuation distribution for short-term wind power, a scenario generation method is proposed. This paper characterizes forecast error via empirical distributions of a set of forecast bins and assumes that wind power fluctuations over unit interval follow $t$ location-scale distribution. An inverse transform sampling from a multivariate normal distribution is adopted to generate a large number of wind power scenarios. The covariance matrix of the multivariate normal distribution is estimated to fit the distribution of historical wind power fluctuations. The proposed scenario generation method is applied to the actual aggregate wind power data in the whole regions of Ireland's Power System. The results indicate that the variability of wind power scenarios can be adjusted by estimating the key range parameter in the exponential covariance structure of a multivariate normal distribution.

Description: Tissue microarrays (TMAs) facilitate the survey of very large numbers of tumors. However, the manual assessment of stained TMA sections constitutes a bottleneck in the pathologist’s work flow. This paper presents a computational pipeline for automatically classifying and scoring breast cancer TMA spots that have been subjected to nuclear immunostaining. Spots are classified based on a bag of visual words approach. Immunohistochemical scoring is performed by computing spot features reflecting the proportion of epithelial nuclei that are stained and the strength of that staining. These are then mapped onto an ordinal scale used by pathologists. Multilayer perceptron classifiers are compared with latent topic models and support vector machines for spot classification, and with Gaussian process ordinal regression and linear models for scoring. Intraobserver variation is also reported. The use of posterior entropy to identify uncertain cases is demonstrated. Evaluation is performed using TMA images stained for progesterone receptor.

Description: The removal of spatially correlated noise is an important step in processing multichannel recordings. Here, a technique termed the adaptive common average reference (ACAR) is presented as an effective and simple method for removing this noise. The ACAR is based on a combination of the well-known common average reference (CAR) and an adaptive noise canceling (ANC) filter. In a convergent process, the CAR provides a reference to an ANC filter, which in turn provides feedback to enhance the CAR. This method was effective on both simulated and real data, outperforming the standard CAR when the amplitude or polarity of the noise changes across channels. In many cases, the ACAR even outperformed independent component analysis. On 16 channels of simulated data, the ACAR was able to attenuate up to approximately 290 dB of noise and could improve signal quality if the original SNR was as high as 5 dB. With an original SNR of 0 dB, the ACAR improved signal quality with only two data channels and performance improved as the number of channels increased. It also performed well under many different conditions for the structure of the noise and signals. Analysis of contaminated electrocorticographic recordings further showed the effectiveness of the ACAR.

Description: This paper proposes a multimodal approach for vessel segmentation of macular optical coherence tomography (OCT) slices along with the fundus image. The method is comprised of two separate stages; the first step is 2-D segmentation of blood vessels in curvelet domain, enhanced by taking advantage of vessel information in crossing OCT slices (named feedback procedure), and improved by suppressing the false positives around the optic nerve head. The proposed method for vessel localization of OCT slices is also enhanced utilizing the fact that retinal nerve fiber layer becomes thicker in the presence of the blood vessels. The second stage of this method is axial localization of the vessels in OCT slices and 3-D reconstruction of the blood vessels. Twenty-four macular spectral 3-D OCT scans of 16 normal subjects were acquired using a Heidelberg HRA OCT scanner. Each dataset consisted of a scanning laser ophthalmoscopy (SLO) image and limited number of OCT scans with size of 496 × 512 (namely, for a data with 19 selected OCT slices, the whole data size was 496 × 512 × 19). The method is developed with least complicated algorithms and the results show considerable improvement in accuracy of vessel segmentation over similar methods to produce a local accuracy of 0.9632 in area of SLO, covered with OCT slices, and the overall accuracy of 0.9467 in the whole SLO image. The results are also demonstrative of a direct relation between the overall accuracy and percentage of SLO coverage by OCT slices.

Description: A brain–computer interface (BCI) based on steady-state visual-evoked potentials (SSVEP) has two difficulties: limitation of the number of commands and uneven probabilities of command execution. To address these problems, the present paper proposes a paradigm of BCI using frequency-modulated visual stimuli. The commands are translated into code words consisting of binary digits, to which visual stimuli with distinct frequencies are assigned. Frequencies of SSVEP are recognized to detect bits, and a command to be executed is determined from the sequence of detected bits. Experimental results show that the proposed paradigm achieves a reliable BCI with higher accuracies and balanced command executing probabilities.

Description: A method for brain monitoring based on measuring the electroencephalogram (EEG) from electrodes placed in-the-ear (ear-EEG) was recently proposed. The objective of this study is to further characterize the ear-EEG and perform a rigorous comparison against conventional on-scalp EEG. This is achieved for both auditory and visual evoked responses, over steady-state and transient paradigms, and across a population of subjects. The respective steady-state responses are evaluated in terms of signal-to-noise ratio and statistical significance, while the qualitative analysis of the transient responses is performed by considering grand averaged event-related potential (ERP) waveforms. The outcomes of this study demonstrate conclusively that the ear-EEG signals, in terms of the signal-to-noise ratio, are on par with conventional EEG recorded from electrodes placed over the temporal region.

Description: Increasing the freedom of communication using conventional row/column (RC) P300 paradigm by naive way (increasing matrix size) may deteriorate inherent distraction effect and interaction speed. In this paper, we propose a two-level predictive (TLP) paradigm by integrating a $3times 3$ two-level matrix paradigm with a statistical language model. The TLP paradigm is evaluated using offline and online data from ten healthy subjects. Significantly larger event-related potentials (ERPs) are evoked by the TLP paradigm compared with the classical $6times 6$ RC. During an online task (correctly spell an English sentence with 57 characters), accuracy and information transfer rate for the TLP are increased by 14.45% and 29.29%, respectively, when compared with the $6times 6$ RC. Time to complete the task is also decreased by 24.61% using TLP. In sharp contrast, an $8times 8$ RC (naive extension of the $6times 6$ RC) consumed $19.18%$ more time than the classical $6times 6$ RC. Furthermore, the statistical language model is also exploited to improve classification accuracy in a Bayesian approach. The proposed Bayesian fusion method is tested offline on data from the online spelling tasks. The results show its potential improvement on single-trial ERP classification.

Description: This paper presents a class of hydration monitor that uses ultrathin, stretchable sheets with arrays of embedded impedance sensors for precise measurement and spatially multiplexed mapping. The devices contain miniaturized capacitive electrodes arranged in a matrix format, capable of integration with skin in a conformal, intimate manner due to the overall skin-like physical properties. These “epidermal” systems noninvasively quantify regional variations in skin hydration, at uniform or variable skin depths. Experimental results demonstrate that the devices possess excellent uniformity, with favorable precision and accuracy. Theoretical models capture the underlying physics of the measurement and enable quantitative interpretation of the experimental results. These devices are appealing for applications ranging from skin care and dermatology, to cosmetology and health/wellness monitoring, with the additional potential for combined use with other classes of sensors for comprehensive, quantitative physiological assessment via the skin.

Description: Microwave radiometers are very sensitive passive sensors that measure the power of the thermal noise within a determined bandwidth. Therefore, any other signal present in the band modifies the value of the measured power, and the corresponding estimated antenna temperature, from which the geophysical parameters are retrieved. Due to the high sensitivity and accuracy required for these instruments, radio frequency interference (RFI) is becoming more and more a serious problem. On one hand, ground-based or global RFI surveys are helping to understand the occurrence and types of RFI sources. If RFI does not necessarily affect the whole bandwidth, or it is not present during the whole integration time, the application of either frequency blanking, time blanking or signal spectrogram techniques can be applied. However, it would be desirable to apply techniques to estimate the RFI signal so that it can be subtracted from the received signal itself so that some useful measurements are still possible. Such a real-time system is currently being developed for RFI detection and mitigation. This work focuses however in the description and performance of a wavelet-based RFI-mitigation technique implemented in a FPGA hardware back-end. The interfering signal is estimated by using the powerful denoising capabilities of the wavelet transform, and it is then subtracted from the total received signal to obtain a RFI-mitigated noise signal.

Description: The narrow-band interference (NBI) is a common jamming signal against synthetic aperture radar (SAR), which can degrade the imaging quality severely. This paper proposes a new method for NBI suppression in the data domain based on the independent component analysis (ICA). In this method, echoes contaminated by the NBI are identified in the frequency domain. Next, time filtering and whitening are performed to the identified echoes. Then, the ICA is carried out to decompose the echoes into a series of basis signals, and the jamming components are selected by thresholding. Finally, the NBI is reconstructed and subtracted from the echoes, and the well-focused SAR imagery is obtained by conventional imaging methods. The proposed method copes well with the time-varying NBI with little signal loss. Results of simulated and measured data have proved the validity of the proposed method.

Description: The Soil Moisture and Ocean Salinity (SMOS) satellite is strongly affected by radio-frequency interference (RFI). A detection algorithm has been developed to accurately obtain the coordinates of the interfering source emitters from the SMOS images. The results obtained from this detection algorithm are regularly used to locate the on-ground sources of interference. This has allowed the identification and termination of over 200 RFI sources observed by SMOS. In the majority of cases, the accuracy of the coordinates provided was better than 4 km, which is a very important achievement considering that SMOS spatial resolution is larger than 35 km and that the contamination of a single RFI can extend to several thousands of kilometers in some cases.

Description: The Soil Moisture and Ocean Salinity (SMOS) radiometer operates within the Earth Exploration Satellite Service passive band at 1400–1427 MHz. Since its launch in November 2009, SMOS images are strongly impacted by radio frequency interference (RFI). So far ${>}{500}$ RFI sources distributed worldwide have been detected. Up to 42% of these RFIs could be suppressed thanks to the co-operation of the National Spectrum Management Authorities. Some of the strongest RFI sources might mask other weaker sources underneath, hence it is expected the total number of RFI detected may increase as strong ones are progressively identified and switched off. Most RFIs are located in Asia and Europe, which together hold ${sim}{73%}$ of the active sources and ${>}{90%}$ of the strongest interference. The areas affected by RFI may experience either an underestimation in the retrieved values of soil moisture and ocean salinity or data loss, with the associated detrimental impact on the scientific return. ESA and the teams participating in SMOS mission have put in place different strategies to alleviate this RFI situation.

Description: Radio-frequency interference (RFI) is probably today's most serious limitation to the accurate retrieval of geophysical parameters from microwave radiometric measurements. Strong RFI inducing a change in the detected power larger than the natural variability is simple to detect. Moderate or weak RFI can be masked by the natural variability of the measurements, passing undetected and corrupting them. A number of techniques have been devised in the past years to detect and, eventually, mitigate RFI present in microwave radiometry measurements: 1) time domain; 2) frequency domain; 3) spectrogram techniques looking for anomalously high power peaks; 4) statistical techniques testing the hypothesis of Gaussianity of the received signal; 5) polarimetric techniques looking for anomalous signatures in the third and fourth Stokes parameters; or 6) wavelet techniques to estimate the RFI signal and cancel it (if any). In this paper, the first four techniques are evaluated with real data gathered with a multifrequency microwave radiometer. It will be shown how spectrogram techniques can detect RFI signals concentrated in narrow frequency bands and/or time intervals that may pass undetected with time-domain and/or frequency-domain techniques alone or with statistical methods. A combined approach is proposed to take advantage of the best performance of each technique. On one side, for strong localized RFI, the approach is spectrogram blanking or, if it is too demanding in terms of computational resources, simple time- and frequency-domain blanking. On the other side, for weak RFI, the approach is the Kurtosis statistical test, which exhibits the best performance among the ten normality tests evaluated, in conjunction with the Anderson–Darling test to detect potential RFI in the blind spots of the Kurtosis test.

Description: The Soil Moisture Active–Passive (SMAP) mission will launch in late 2014 and will carry a combined L-band radiometer/radar instrument for the retrieval of global soil moisture and surface freeze-thaw state. Radio frequency interference (RFI) is a known challenge for Earth remote sensing in the L-band portion of the spectrum. This paper addresses efforts to characterize and mitigate RFI for the SMAP radar. A model for the RFI environment due to surface-based emitters is developed, and is shown to agree well with the observations of currently operating L-band radar systems. An analysis of the environment due to space-based emitters is also presented. Techniques to mitigate RFI in the radar band are described, and are shown to perform sufficiently well to meet the stringent SMAP measurement requirements. A companion paper addresses the different issues encountered with RFI in the radiometer band.

Description: Radio frequency interference (RFI) is a known issue in low-frequency radar remote sensing. In synthetic aperture radar (SAR) image processing, RFI can cause severe degradation of image quality, distortion of polarimetric signatures, and an increase of the SAR phase noise level. To address this issue, a processing system was developed that is capable of reliably detecting, characterizing, and mitigating RFI signatures in SAR observations. In addition to being the basis for image correction, the robust RFI-detection algorithms developed in this paper are used to retrieve a wealth of RFI-related information that allows for mapping, characterizing, and classifying RFI signatures across large spatial scales. The extracted RFI information is expected to be valuable input for SAR-system design, sensor operations, and the development of effective RFI-mitigation strategies. The concepts of RFI detection, analysis, and mapping are outlined. Large-scale RFI mapping results are shown. In case studies, the benefit of detailed RFI information for customized RFI filtering and sensor operations is exemplified.

Description: Radio-frequency interference (RFI) is a persistent threat to Earth-observing microwave radiometers. A number of test statistics are used for radiometric RFI detection. This paper presents a new RFI detection method that uses the information theoretic quantity known as negentropy. In particular, we study six negentropy-based test statistics and compare their performance against kurtosis, Jarque–Bera, Anderson–Darling, and Shapiro–Wilk normality tests for specific RFI signal models. The Neyman–Pearson decision rule is used to develop receiver operating characteristic curves for each test statistic. We show that although negentropy can be used to detect RFI, it does not outperform kurtosis, except for the kurtosis blind-spot case.

Description: Measurements of both the Soil Moisture and Ocean Salinity (SMOS) and Aquarius L-band microwave radiometers show a significant presence of radio frequency interference (RFI), although they operate in a protected frequency band where transmission is prohibited. RFI detection and mitigation remain a challenging problem for both missions, especially for low or moderate (i.e., on the order of 10 K or less) amplitude contributions. An algorithm for low-level source detection and mitigation is already included in Aquarius data sets, and both Aquarius and SMOS have distinct attributes that can potentially enable further improvements in detection and mitigation of these sources to some degree. The combination of SMOS and Aquarius data sets may enable further future improvements as well. Initial efforts toward this goal are reported in this paper. Similarities and differences in RFI effects on SMOS and Aquarius are examined, with a particular focus on instrument properties that cause differences in received RFI power in SMOS and Aquarius observations of a specific source. A study is also performed of SMOS observations for regions reported by Aquarius to contain “low-level” RFI. It is shown that the detection of these sources in the SMOS data set is challenging and that the dependence of the SMOS third and fourth Stokes parameters on incidence angle makes the polarimetric features of SMOS difficult to utilize for low-level source detection. However, an angular fitting procedure suggested previously in the literature can, in some cases, detect such sources in horizontal and vertical polarizations.