Frequency band selection of radars for buried object detection

Choice of the operational frequency is one of the most responsible parts of any radar design process. Parameters of radars for buried object detection (BOD) are very sensitive to both carrier frequency and ranging signal bandwidth. Such radars have a specific propagation environment with a strong frequency-dependent attenuation and, as a result, short operational range. This fact dictates some features of the radar's parameters: wideband signal-to provide a high range resolution (fractions of a meter) and a low carrier frequency (tens or hundreds megahertz) for deeper penetration. The requirement to have a wideband ranging signal and low carrier frequency are partly in contradiction. As a result, low-frequency (LF) ultrawide-band (UWB) signals are used. The major goal of this paper is to examine the influence of the frequency band choice on the radar performance and develop relevant methodologies for BOD radar design and optimization. In this article, high-efficient continuous wave (CW) signals with most advanced stepped frequency (SF) modulation are considered; however, the main conclusions can be applied to any kind of ranging signals.

Cloud-radar observations of insects in the UK convective boundary layer

Radar has been applied to the study of insect migration for almost 40 years, but most entomological radars operate at X-band (9.4 GHz, 3.2 cm wavelength), and can only detect individuals of relatively large species, such as migratory grasshoppers and noctuid moths, over all of their flight altitudes. Many insects (including economically important species) are much smaller than this, but development of the requisite higher power and/or higher frequency radar systems to detect these species is often prohibitively expensive. In this paper, attention is focussed upon the uses of some recently-deployed meteorological sensing devices to investigate insect migratory flight behaviour, and especially its interactions with boundary layer processes. Records were examined from the vertically-pointing 35 GHz ‘Copernicus’ and 94 GHz ‘Galileo’ cloud radars at Chilbolton (Hampshire, England) for 12 cloudless and convective occasions in summer 2003, and one of these occasions (13 July) is presented in detail. Insects were frequently found at heights above aerosol particles, which represent passive tracers, indicating active insect movement. It was found that insect flight above the convective boundary layer occurs most often during the morning. The maximum radar reflectivity (an indicator of aerial insect biomass) was found to be positively correlated with maximum screen temperature.

Automotive Radars:towards smarter vehicles

Radars are expected to become the main sensors in various civilian applications, especially for autonomous driving. Their success is mainly due to the availability of low cost integrated devices, equipped with compact antenna arrays, and computationally efficient signal processing techniques. This thesis focuses on the study and the development of different deterministic and learning based techniques for colocated multiple-input multiple-output (MIMO) radars. In particular, after providing an overview on the architecture of these devices, the problem of detecting and estimating multiple targets in stepped frequency continuous wave (SFCW) MIMO radar systems is investigated and different deterministic techniques solving it are illustrated. Moreover, novel solutions, based on an approximate maximum likelihood approach, are developed. The accuracy achieved by all the considered algorithms is assessed on the basis of the raw data acquired from low power wideband radar devices. The results demonstrate that the developed algorithms achieve reasonable accuracies, but at the price of different computational efforts. Another important technical problem investigated in this thesis concerns the exploitation of machine learning and deep learning techniques in the field of colocated MIMO radars. In this thesis, after providing a comprehensive overview of the machine learning and deep learning techniques currently being considered for use in MIMO radar systems, their performance in two different applications is assessed on the basis of synthetically generated and experimental datasets acquired through a commercial frequency modulated continuous wave (FMCW) MIMO radar. Finally, the application of colocated MIMO radars to autonomous driving in smart agriculture is illustrated.

Applicabilité de la texture couleur à la différentiation des classes d’occupation du territoire sur des images satellitales multispectrales

La texture est un élément clé pour l’interprétation des images de télédétection à fine résolution spatiale. L’intégration de l’information texturale dans un processus de classification automatisée des images se fait habituellement via des images de texture, souvent créées par le calcul de matrices de co-occurrences (MCO) des niveaux de gris. Une MCO est un histogramme des fréquences d’occurrence des paires de valeurs de pixels présentes dans les fenêtres locales, associées à tous les pixels de l’image utilisée; une paire de pixels étant définie selon un pas et une orientation donnés. Les MCO permettent le calcul de plus d’une dizaine de paramètres décrivant, de diverses manières, la distribution des fréquences, créant ainsi autant d’images texturales distinctes. L’approche de mesure des textures par MCO a été appliquée principalement sur des images de télédétection monochromes (ex. images panchromatiques, images radar monofréquence et monopolarisation). En imagerie multispectrale, une unique bande spectrale, parmi celles disponibles, est habituellement choisie pour générer des images de texture. La question que nous avons posée dans cette recherche concerne justement cette utilisation restreinte de l’information texturale dans le cas des images multispectrales. En fait, l’effet visuel d’une texture est créé, non seulement par l’agencement particulier d’objets/pixels de brillance différente, mais aussi de couleur différente. Plusieurs façons sont proposées dans la littérature pour introduire cette idée de la texture à plusieurs dimensions. Parmi celles-ci, deux en particulier nous ont intéressés dans cette recherche. La première façon fait appel aux MCO calculées bande par bande spectrale et la seconde utilise les MCO généralisées impliquant deux bandes spectrales à la fois. Dans ce dernier cas, le procédé consiste en le calcul des fréquences d’occurrence des paires de valeurs dans deux bandes spectrales différentes. Cela permet, en un seul traitement, la prise en compte dans une large mesure de la « couleur » des éléments de texture. Ces deux approches font partie des techniques dites intégratives. Pour les distinguer, nous les avons appelées dans cet ouvrage respectivement « textures grises » et « textures couleurs ». Notre recherche se présente donc comme une analyse comparative des possibilités offertes par l’application de ces deux types de signatures texturales dans le cas spécifique d’une cartographie automatisée des occupations de sol à partir d’une image multispectrale. Une signature texturale d’un objet ou d’une classe d’objets, par analogie aux signatures spectrales, est constituée d’une série de paramètres de texture mesurés sur une bande spectrale à la fois (textures grises) ou une paire de bandes spectrales à la fois (textures couleurs). Cette recherche visait non seulement à comparer les deux approches intégratives, mais aussi à identifier la composition des signatures texturales des classes d’occupation du sol favorisant leur différentiation : type de paramètres de texture / taille de la fenêtre de calcul / bandes spectrales ou combinaisons de bandes spectrales. Pour ce faire, nous avons choisi un site à l’intérieur du territoire de la Communauté Métropolitaine de Montréal (Longueuil) composé d’une mosaïque d’occupations du sol, caractéristique d’une zone semi urbaine (résidentiel, industriel/commercial, boisés, agriculture, plans d’eau…). Une image du satellite SPOT-5 (4 bandes spectrales) de 10 m de résolution spatiale a été utilisée dans cette recherche. Puisqu’une infinité d’images de texture peuvent être créées en faisant varier les paramètres de calcul des MCO et afin de mieux circonscrire notre problème nous avons décidé, en tenant compte des études publiées dans ce domaine : a) de faire varier la fenêtre de calcul de 3*3 pixels à 21*21 pixels tout en fixant le pas et l’orientation pour former les paires de pixels à (1,1), c'est-à-dire à un pas d’un pixel et une orientation de 135°; b) de limiter les analyses des MCO à huit paramètres de texture (contraste, corrélation, écart-type, énergie, entropie, homogénéité, moyenne, probabilité maximale), qui sont tous calculables par la méthode rapide de Unser, une approximation des matrices de co-occurrences, c) de former les deux signatures texturales par le même nombre d’éléments choisis d’après une analyse de la séparabilité (distance de Bhattacharya) des classes d’occupation du sol; et d) d’analyser les résultats de classification (matrices de confusion, exactitudes, coefficients Kappa) par maximum de vraisemblance pour conclure sur le potentiel des deux approches intégratives; les classes d’occupation du sol à reconnaître étaient : résidentielle basse et haute densité, commerciale/industrielle, agricole, boisés, surfaces gazonnées (incluant les golfs) et plans d’eau. Nos principales conclusions sont les suivantes a) à l’exception de la probabilité maximale, tous les autres paramètres de texture sont utiles dans la formation des signatures texturales; moyenne et écart type sont les plus utiles dans la formation des textures grises tandis que contraste et corrélation, dans le cas des textures couleurs, b) l’exactitude globale de la classification atteint un score acceptable (85%) seulement dans le cas des signatures texturales couleurs; c’est une amélioration importante par rapport aux classifications basées uniquement sur les signatures spectrales des classes d’occupation du sol dont le score est souvent situé aux alentours de 75%; ce score est atteint avec des fenêtres de calcul aux alentours de11*11 à 15*15 pixels; c) Les signatures texturales couleurs offrant des scores supérieurs à ceux obtenus avec les signatures grises de 5% à 10%; et ce avec des petites fenêtres de calcul (5*5, 7*7 et occasionnellement 9*9) d) Pour plusieurs classes d’occupation du sol prises individuellement, l’exactitude dépasse les 90% pour les deux types de signatures texturales; e) une seule classe est mieux séparable du reste par les textures grises, celle de l’agricole; f) les classes créant beaucoup de confusions, ce qui explique en grande partie le score global de la classification de 85%, sont les deux classes du résidentiel (haute et basse densité). En conclusion, nous pouvons dire que l’approche intégrative par textures couleurs d’une image multispectrale de 10 m de résolution spatiale offre un plus grand potentiel pour la cartographie des occupations du sol que l’approche intégrative par textures grises. Pour plusieurs classes d’occupations du sol un gain appréciable en temps de calcul des paramètres de texture peut être obtenu par l’utilisation des petites fenêtres de traitement. Des améliorations importantes sont escomptées pour atteindre des exactitudes de classification de 90% et plus par l’utilisation des fenêtres de calcul de taille variable adaptées à chaque type d’occupation du sol. Une méthode de classification hiérarchique pourrait être alors utilisée afin de séparer les classes recherchées une à la fois par rapport au reste au lieu d’une classification globale où l’intégration des paramètres calculés avec des fenêtres de taille variable conduirait inévitablement à des confusions entre classes.

Atmospheric electricity coupling between earthquake regions and the ionosphere

We propose a mechanism to explain suggested links between seismic activity and ionospheric changes detected overhead. Specifically, we explain changes in the natural extremely low-frequency (ELF) radio noise recently observed in the topside ionosphere aboard the DEMETER satellite at night, before major earthquakes. Our mechanism utilises increased electrical conductivity of surface layer air before a major earthquake, which reduces the surface-ionosphere electrical resistance. This increases the vertical fair weather current, and (to maintain continuity of electron flow) lowers the ionosphere. Magnitudes of crucial parameters are estimated and found to be consistent with observations. Natural variability in ionospheric and atmospheric electrical properties is evaluated, and may be overcome using a hybrid detection approach. Suggested experiments to investigate the mechanism involve measuring the cut-off frequency of ELF “tweeks”, the amplitude and phase of very low frequency radio waves in the Earth–ionosphere waveguide, or medium frequency radar, incoherent scatter or rocket studies of the lower ionospheric electron density.

Retrieval of mesospheric electron densities using an optimal estimation inverse method

We present a new method to determine mesospheric electron densities from partially reflected medium frequency radar pulses. The technique uses an optimal estimation inverse method and retrieves both an electron density profile and a gradient electron density profile. As well as accounting for the absorption of the two magnetoionic modes formed by ionospheric birefringence of each radar pulse, the forward model of the retrieval parameterises possible Fresnel scatter of each mode by fine electronic structure, phase changes of each mode due to Faraday rotation and the dependence of the amplitudes of the backscattered modes upon pulse width. Validation results indicate that known profiles can be retrieved and that χ2 tests upon retrieval parameters satisfy validity criteria. Application to measurements shows that retrieved electron density profiles are consistent with accepted ideas about seasonal variability of electron densities and their dependence upon nitric oxide production and transport.

G band atmospheric radars: new frontiers in cloud physics

Clouds and associated precipitation are the largest source of uncertainty in current weather and future climate simulations. Observations of the microphysical, dynamical and radiative processes that act at cloud scales are needed to improve our understanding of clouds. The rapid expansion of ground-based super-sites and the availability of continuous profiling and scanning multi-frequency radar observations at 35 and 94 GHz have significantly improved our ability to probe the internal structure of clouds in high temporal-spatial resolution, and to retrieve quantitative cloud and precipitation properties. However, there are still gaps in our ability to probe clouds due to large uncertainties in the retrievals. The present work discusses the potential of G band (frequency between 110 and 300 GHz) Doppler radars in combination with lower frequencies to further improve the retrievals of microphysical properties. Our results show that, thanks to a larger dynamic range in dual-wavelength reflectivity, dual-wavelength attenuation and dual-wavelength Doppler velocity (with respect to a Rayleigh reference), the inclusion of frequencies in the G band can significantly improve current profiling capabilities in three key areas: boundary layer clouds, cirrus and mid-level ice clouds, and precipitating snow.

Near real time surface ocean velocity, Tyrrhenian Sea and Ligurian Sea in September 2016

The data set consists of maps of total velocity of the surface current in the North-Western Tyrrhenian Sea and Ligurian Sea averaged over a time interval of 1 hour around the cardinal hour. Surface ocean velocities estimated by HF Radar are representative of the upper 0.3-2.5 meters of the ocean. Total velocities are derived using least square fit that maps radial velocities measured from individual sites onto a cartesian grid. The final product is a map of the horizontal components of the ocean currents on a regular grid in the area of overlap of two or more radar stations.

Near real time surface ocean velocity in the Bay of Biscay - Basque Country from Mai to September 2016

The data set consists of maps of total velocity of the surface current in the Southeastern Bay of Biscay averaged over a time interval of 1 hour around the cardinal hour. Surface ocean velocities estimated by this HF Radar(4.65 MHz) are representative of the upper 2-3 meters of the ocean. The main objective of near real time processing is to produce the best product from available data at the time of processing. Total velocities are derived using least square fit that maps radial velocities measured from individual sites onto a cartesian grid. The final product is a map of the horizontal components of the ocean currents on a regular grid in the area of overlap of two or more radar stations.

Real-time surface current data in the Ibiza Channel from January to October 2016

The data set consists of maps of total velocity of surface currents in the Ibiza Channel, derived from HF radar measurements.

New Possibilities for Geophysical Parameter Retrievals Opened by GCOM-W1 AMSR2

A new approach to retrieve sea surface wind speed (SWS) in tropical cyclones (TCs) from the Advanced Microwave Scanning Radiometer 2 (AMSR2) data is presented. Analysis of all six AMSR2 C- and X-band channel measurements over TCs is shown to efficiently help to separate the rain contribution. Corrected measurements at 6.9 and 10.65 GHz are then used to retrieve the SWS. Spatial and temporal collocation of AMSR2 and tropical rain measurement mission (TRMM) microwave instrument (TMI) data is then further used to empirically relate TMI rain rate (RR) product to RR estimates from AMSR2 in hurricanes. SWS estimates are validated with measurements from the stepped frequency microwave radiometer (SFMR). As further tested, more than 100 North Atlantic and North Pacific TCs are analyzed for the 2012–2014 period. Despite few particular cases, most SWS fields are in a very good agreement with TC center data on maximum wind speeds, radii of storm, and hurricane winds. As also compared, very high consistency between AMSR2 and L-band SMOS wind speed estimates are obtained, especially for the super typhoon Haiyan, to prove the high potential of AMSR2 measurements in TCs.

Global oceanic precipitation: A joint view by TOPEX and the TOPEX microwave radiometer

The TOPEX/POSEIDON mission offers the first opportunity to observe rain cells over the ocean by a dual-frequency radar altimeter (TOPEX) and simultaneously observe their natural radiative properties by a three-frequency radiometer (TOPEX microwave radiometer (TMR)). This work is a feasibility study aimed at understanding the capability and potential of the active/passive TOPEX/TMR system for oceanic rainfall detection. On the basis of past experiences in rain flagging, a joint TOPEX/TMR rain probability index is proposed. This index integrates several advantages of the two sensors and provides a more reliable rain estimate than the radiometer alone. One year's TOPEX/TMR TMR data are used to test the performance of the index. The resulting rain frequency statistics show quantitative agreement with those obtained from the Comprehensive Ocean-Atmosphere Data Set (COADS) in the Intertropical Convergence Zone (ITCZ), while qualitative agreement is found for other regions of the world ocean. A recent finding that the latitudinal frequency of precipitation over the Southern Ocean increases steadily toward the Antarctic continent is confirmed by our result. Annual and seasonal precipitation maps are derived from the index. Notable features revealed include an overall similarity in rainfall pattern from the Pacific, the Atlantic, and the Indian Oceans and a general phase reversal between the two hemispheres, as well as a number of regional anomalies in terms of rain intensity. Comparisons with simultaneous Global Precipitation Climatology Project (GPCP) multisatellite precipitation rate and COADS rain climatology suggest that systematic differences also exist. One example is that the maximum rainfall in the ITCZ of the Indian Ocean appears to be more intensive and concentrated in our result compared to that of the GPCP. Another example is that the annual precipitation produced by TOPEX/TMR is constantly higher than those from GPCP and COADS in the extratropical regions of the northern hemisphere, especially in the northwest Pacific Ocean. Analyses of the seasonal variations of prominent rainy and dry zones in the tropics and subtropics show various behaviors such as systematic migration, expansion and contraction, merging and breakup, and pure intensity variations, The seasonality of regional features is largely influenced by local atmospheric events such as monsoon, storm, or snow activities. The results of this study suggest that TOPEX and its follow-on may serve as a complementary sensor to the special sensor microwave/imager in observing global oceanic precipitation.

Quantifying errors due to frequency changes and target location uncertainty for radar refractivity retrievals

Radar refractivity retrievals can capture near-surface humidity changes, but noisy phase changes of the ground clutter returns limit the accuracy for both klystron- and magnetron-based systems. Observations with a C-band (5.6 cm) magnetron weather radar indicate that the correction for phase changes introduced by local oscillator frequency changes leads to refractivity errors no larger than 0.25 N units: equivalent to a relative humidity change of only 0.25% at 20°C. Requested stable local oscillator (STALO) frequency changes were accurate to 0.002 ppm based on laboratory measurements. More serious are the random phase change errors introduced when targets are not at the range-gate center and there are changes in the transmitter frequency (ΔfTx) or the refractivity (ΔN). Observations at C band with a 2-μs pulse show an additional 66° of phase change noise for a ΔfTx of 190 kHz (34 ppm); this allows the effect due to ΔN to be predicted. Even at S band with klystron transmitters, significant phase change noise should occur when a large ΔN develops relative to the reference period [e.g., ~55° when ΔN = 60 for the Next Generation Weather Radar (NEXRAD) radars]. At shorter wavelengths (e.g., C and X band) and with magnetron transmitters in particular, refractivity retrievals relative to an earlier reference period are even more difficult, and operational retrievals may be restricted to changes over shorter (e.g., hourly) periods of time. Target location errors can be reduced by using a shorter pulse or identified by a new technique making alternate measurements at two closely spaced frequencies, which could even be achieved with a dual–pulse repetition frequency (PRF) operation of a magnetron transmitter.