REFLECT : logiciel de restitution des réflectances au sol pour l’amélioration de la qualité de l'information extraite des images satellitales à haute résolution spatiale

RÉSUMÉ - Les images satellitales multispectrales, notamment celles à haute résolution spatiale (plus fine que 30 m au sol), représentent une source d’information inestimable pour la prise de décision dans divers domaines liés à la gestion des ressources naturelles, à la préservation de l’environnement ou à l’aménagement et la gestion des centres urbains. Les échelles d’étude peuvent aller du local (résolutions plus fines que 5 m) à des échelles régionales (résolutions plus grossières que 5 m). Ces images caractérisent la variation de la réflectance des objets dans le spectre qui est l’information clé pour un grand nombre d’applications de ces données. Or, les mesures des capteurs satellitaux sont aussi affectées par des facteurs « parasites » liés aux conditions d’éclairement et d’observation, à l’atmosphère, à la topographie et aux propriétés des capteurs. Deux questions nous ont préoccupé dans cette recherche. Quelle est la meilleure approche pour restituer les réflectances au sol à partir des valeurs numériques enregistrées par les capteurs tenant compte des ces facteurs parasites ? Cette restitution est-elle la condition sine qua non pour extraire une information fiable des images en fonction des problématiques propres aux différents domaines d’application des images (cartographie du territoire, monitoring de l’environnement, suivi des changements du paysage, inventaires des ressources, etc.) ? Les recherches effectuées les 30 dernières années ont abouti à une série de techniques de correction des données des effets des facteurs parasites dont certaines permettent de restituer les réflectances au sol. Plusieurs questions sont cependant encore en suspens et d’autres nécessitent des approfondissements afin, d’une part d’améliorer la précision des résultats et d’autre part, de rendre ces techniques plus versatiles en les adaptant à un plus large éventail de conditions d’acquisition des données. Nous pouvons en mentionner quelques unes : - Comment prendre en compte des caractéristiques atmosphériques (notamment des particules d’aérosol) adaptées à des conditions locales et régionales et ne pas se fier à des modèles par défaut qui indiquent des tendances spatiotemporelles à long terme mais s’ajustent mal à des observations instantanées et restreintes spatialement ? - Comment tenir compte des effets de « contamination » du signal provenant de l’objet visé par le capteur par les signaux provenant des objets environnant (effet d’adjacence) ? ce phénomène devient très important pour des images de résolution plus fine que 5 m; - Quels sont les effets des angles de visée des capteurs hors nadir qui sont de plus en plus présents puisqu’ils offrent une meilleure résolution temporelle et la possibilité d’obtenir des couples d’images stéréoscopiques ? - Comment augmenter l’efficacité des techniques de traitement et d’analyse automatique des images multispectrales à des terrains accidentés et montagneux tenant compte des effets multiples du relief topographique sur le signal capté à distance ? D’autre part, malgré les nombreuses démonstrations par des chercheurs que l’information extraite des images satellitales peut être altérée à cause des tous ces facteurs parasites, force est de constater aujourd’hui que les corrections radiométriques demeurent peu utilisées sur une base routinière tel qu’est le cas pour les corrections géométriques. Pour ces dernières, les logiciels commerciaux de télédétection possèdent des algorithmes versatiles, puissants et à la portée des utilisateurs. Les algorithmes des corrections radiométriques, lorsqu’ils sont proposés, demeurent des boîtes noires peu flexibles nécessitant la plupart de temps des utilisateurs experts en la matière. Les objectifs que nous nous sommes fixés dans cette recherche sont les suivants : 1) Développer un logiciel de restitution des réflectances au sol tenant compte des questions posées ci-haut. Ce logiciel devait être suffisamment modulaire pour pouvoir le bonifier, l’améliorer et l’adapter à diverses problématiques d’application d’images satellitales; et 2) Appliquer ce logiciel dans différents contextes (urbain, agricole, forestier) et analyser les résultats obtenus afin d’évaluer le gain en précision de l’information extraite par des images satellitales transformées en images des réflectances au sol et par conséquent la nécessité d’opérer ainsi peu importe la problématique de l’application. Ainsi, à travers cette recherche, nous avons réalisé un outil de restitution de la réflectance au sol (la nouvelle version du logiciel REFLECT). Ce logiciel est basé sur la formulation (et les routines) du code 6S (Seconde Simulation du Signal Satellitaire dans le Spectre Solaire) et sur la méthode des cibles obscures pour l’estimation de l’épaisseur optique des aérosols (aerosol optical depth, AOD), qui est le facteur le plus difficile à corriger. Des améliorations substantielles ont été apportées aux modèles existants. Ces améliorations concernent essentiellement les propriétés des aérosols (intégration d’un modèle plus récent, amélioration de la recherche des cibles obscures pour l’estimation de l’AOD), la prise en compte de l’effet d’adjacence à l’aide d’un modèle de réflexion spéculaire, la prise en compte de la majorité des capteurs multispectraux à haute résolution (Landsat TM et ETM+, tous les HR de SPOT 1 à 5, EO-1 ALI et ASTER) et à très haute résolution (QuickBird et Ikonos) utilisés actuellement et la correction des effets topographiques l’aide d’un modèle qui sépare les composantes directe et diffuse du rayonnement solaire et qui s’adapte également à la canopée forestière. Les travaux de validation ont montré que la restitution de la réflectance au sol par REFLECT se fait avec une précision de l’ordre de ±0.01 unités de réflectance (pour les bandes spectrales du visible, PIR et MIR), même dans le cas d’une surface à topographie variable. Ce logiciel a permis de montrer, à travers des simulations de réflectances apparentes à quel point les facteurs parasites influant les valeurs numériques des images pouvaient modifier le signal utile qui est la réflectance au sol (erreurs de 10 à plus de 50%). REFLECT a également été utilisé pour voir l’importance de l’utilisation des réflectances au sol plutôt que les valeurs numériques brutes pour diverses applications courantes de la télédétection dans les domaines des classifications, du suivi des changements, de l’agriculture et de la foresterie. Dans la majorité des applications (suivi des changements par images multi-dates, utilisation d’indices de végétation, estimation de paramètres biophysiques, …), la correction des images est une opération cruciale pour obtenir des résultats fiables. D’un point de vue informatique, le logiciel REFLECT se présente comme une série de menus simples d’utilisation correspondant aux différentes étapes de saisie des intrants de la scène, calcul des transmittances gazeuses, estimation de l’AOD par la méthode des cibles obscures et enfin, l’application des corrections radiométriques à l’image, notamment par l’option rapide qui permet de traiter une image de 5000 par 5000 pixels en 15 minutes environ. Cette recherche ouvre une série de pistes pour d’autres améliorations des modèles et méthodes liés au domaine des corrections radiométriques, notamment en ce qui concerne l’intégration de la FDRB (fonction de distribution de la réflectance bidirectionnelle) dans la formulation, la prise en compte des nuages translucides à l’aide de la modélisation de la diffusion non sélective et l’automatisation de la méthode des pentes équivalentes proposée pour les corrections topographiques.

Ajuste radiométrico de imagens aéreas digitais

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Utilizzo delle bande dell’Infrarosso Termico in tecniche di remote sensing per applicazioni in ambito urbano

Il presente studio si concentra sulle diverse applicazioni del telerilevamento termico in ambito urbano. Vengono inizialmente descritti la radiazione infrarossa e le sue interazioni con l’atmosfera terrestre, le leggi principali che regolano lo scambio di calore per irraggiamento, le caratteristiche dei sensori e le diverse applicazioni di termografia. Successivamente sono trattati nel dettaglio gli aspetti caratteristici della termografia da piattaforma satellitare, finalizzata principalmente alla valutazione del fenomeno dell'Urban Heat Island; vengono descritti i sensori disponibili, le metodologie di correzione per gli effetti atmosferici, per la stima dell'emissività delle superfici e per il calcolo della temperatura superficiale dei pixels. Viene quindi illustrata la sperimentazione effettuata sull'area di Bologna mediante immagini multispettrali ASTER: i risultati mostrano come sull'area urbana sia riscontrabile la presenza dell'Isola di Calore Urbano, anche se la sua quantificazione risulta complessa. Si procede quindi alla descrizione di potenzialità e limiti della termografia aerea, dei suoi diversi utilizzi, delle modalità operative di rilievo e degli algoritmi utilizzati per il calcolo della temperatura superficiale delle coperture edilizie. Tramite l’analisi di alcune esperienze precedenti vengono trattati l’influenza dell’atmosfera, la modellazione dei suoi effetti sulla radianza rilevata, i diversi metodi per la stima dell’emissività. Viene quindi introdotto il progetto europeo Energycity, finalizzato alla creazione di un sistema GeoWeb di supporto spaziale alle decisioni per la riduzione di consumi energetici e produzione di gas serra su sette città dell'Europa Centrale. Vengono illustrate le modalità di rilievo e le attività di processing dei datasets digitali per la creazione di mappe di temperatura superficiale da implementare nel sistema SDSS. Viene infine descritta la sperimentazione effettuata sulle immagini termiche acquisite nel febbraio 2010 sulla città di Treviso, trasformate in un mosaico georiferito di temperatura radiometrica tramite correzioni geometriche e radiometriche; a seguito della correzione per l’emissività quest’ultimo verrà trasformato in un mosaico di temperatura superficiale.

Airborne hyperspectral image data of Heron Reef, Australia

This airborne hyperspectral (19 bands) image data of Heron Reef, Great Barrier Reef, Australia is derived from Compact Airborne Spectrographic Imager (CASI) data acquired on 1st and 3rd of July 2002, latitude -23.45, longitude 151.92. Processing and correction to at-surface data was completed by Karen Joyce (Joyce, 2004). Raw imagery consisted several images corresponding to the number of flight paths taken to cover the entire Heron Reef. Spatial resolution is one meter. Radiometric corrections converted the at-sensor digital number values to at surface spectral radiance values using sensor specific calibration coefficients and CSIRO's c-WomBat-c atmospheric correction software. Geometric corrections were done using field collected coordinates of features identified in the image. Projection used was Universal Transverse Mercator Zone 56 South and Datum used was WGS 84. Image data is in TIFF format.

Along-scan radiometric gradient in Landsat images: understanding the relation between floristic and remote sensing information in tropical forests

The along-scan radiometric gradient causes severe interpretation problems in Landsat images of tropical forests. It creates a decreasing trend in pixel values with the column number of the image. In practical applications it has been corrected assuming the trend to be linear within structurally similar forests. This has improved the relation between floristic and remote sensing information, but just in some cases. I use 3 Landsat images and 105 floristic inventories to test the assumption of linearity, and to examine how the gradient and linear corrections affect the relation between floristic and Landsat data. Results suggest the gradient to be linear in infrared bands. Also, the relation between floristic and Landsat data could be conditioned by the distribution of the sampling sites and the direction in which images are mosaicked. Additionally, there seems to be a conjunction between the radiometric gradient and a natural east-west vegetation gradient common in Western Amazonia. This conjunction might have enhanced artificially correlations between field and remotely-sensed information in previous studies. Linear corrections may remove such artificial enhancement, but along with true and relevant spectral information about floristic patterns, because they can´t separate the radiometric gradient from a natural one.

Modelos de regressão lineares para estimativa de produtividade da soja no oeste do Paraná, utilizando dados espectrais

The main objective of this work was to evaluate the linear regression between spectral response and soybean yield in regional scale. In this study were monitored 36 municipalities from the west region of the states of Parana using five images of Landsat 5/TM during 2004/05 season. The spectral response was converted in physical values, apparent and surface reflectances, by radiometric transformation and atmospheric corrections and both used to calculate NDVI and GVI vegetation indices. Those ones were compared by multiple and simple regression with government official yield values (IBGE). Diagnostic processing method to identify influents values or collinearity was applied to the data too. The results showed that the mean surface reflectance value from all images was more correlated with yield than individual dates. Further, the multiple regressions using all dates and both vegetation indices gave better results than simple regression.

Hyperpolarizabilities of the methanol molecule: A CCSD calculation including vibrational corrections

In this work we present the results for hyperpolarizabilities of the methanol molecule including vibrational corrections and electron correlation effects at the CCSD level. Comparisons to random phase approximation results previously reported show that the electron correlation is in general important for both electronic contribution and vibrational corrections. The role played by the anharmonicities on the calculations of the vibrational corrections has also been analyzed and the obtained results indicate that the anharmonic terms are important for the dc-Pockels and dc-Kerr effects. For the other nonlinear optical properties studied the double-harmonic approximation is found to be suitable. Comparison to available experimental result in gas phase for the dc-second harmonic generation second hyperpolarizability shows a very good agreement with the electronic contribution calculated here while our total value is 14% larger than the experimental value.

Finite-size corrections to entanglement in quantum critical systems

We analyze the finite-size corrections to entanglement in quantum critical systems. By using conformal symmetry and density functional theory, we discuss the structure of the finite-size contributions to a general measure of ground state entanglement, which are ruled by the central charge of the underlying conformal field theory. More generally, we show that all conformal towers formed by an infinite number of excited states (as the size of the system L -> infinity) exhibit a unique pattern of entanglement, which differ only at leading order (1/L)(2). In this case, entanglement is also shown to obey a universal structure, given by the anomalous dimensions of the primary operators of the theory. As an illustration, we discuss the behavior of pairwise entanglement for the eigenspectrum of the spin-1/2 XXZ chain with an arbitrary length L for both periodic and twisted boundary conditions.

Design and optimization of na ultra wideband and compact microwave antenna for radiometric monitoring of brain temperature

We present the modeling efforts on antenna design and frequency selection to monitor brain temperature during prolonged surgery using noninvasive microwave radiometry. A tapered log-spiral antenna design is chosen for its wideband characteristics that allow higher power collection from deep brain. Parametric analysis with the software HFSS is used to optimize antenna performance for deep brain temperature sensing. Radiometric antenna efficiency (eta) is evaluated in terms of the ratio of power collected from brain to total power received by the antenna. Anatomical information extracted from several adult computed tomography scans is used to establish design parameters for constructing an accurate layered 3-D tissue phantom. This head phantom includes separate brain and scalp regions, with tissue equivalent liquids circulating at independent temperatures on either side of an intact skull. The optimized frequency band is 1.1-1.6 GHz producing an average antenna efficiency of 50.3% from a two turn log-spiral antenna. The entire sensor package is contained in a lightweight and low-profile 2.8 cm diameter by 1.5 cm high assembly that can be held in place over the skin with an electromagnetic interference shielding adhesive patch. The calculated radiometric equivalent brain temperature tracks within 0.4 degrees C of the measured brain phantom temperature when the brain phantom is lowered 10. C and then returned to the original temperature (37 degrees C) over a 4.6-h experiment. The numerical and experimental results demonstrate that the optimized 2.5-cm log-spiral antenna is well suited for the noninvasive radiometric sensing of deep brain temperature.

QED and relativistic corrections in superheavy elements

In this paper we review the different relativistic and QED contributions to energies, ionic radii, transition probabilities and Landé g-factors in super-heavy elements, with the help of the MultiConfiguration Dirac-Fock method (MCDF). The effects of taking into account the Breit interaction to all orders by including it in the self-consistent field process are demonstrated. State of the art radiative corrections are included in the calculation and discussed. We also study the non-relativistic limit of MCDF calculation and find that the non-relativistic offset can be unexpectedly large.

Radiometric studies on the oxidation of (1-14c) fatty acids by drug-susceptible and drug-resistant mycobacteria

A radiometric assay system has been used to study oxidation patterns of (1-14C) fatty acids by drug-susceptible and drug-resistant organisms of the genus Mycobacterium. Two strains of M. tuberculosis susceptible to all drugs, H37Rv and Erdman, were used. Drug-resistant organisms included in this investigation were M. tuberculosis H37Rv resistant to 5 ug/ml isoniazid, M. bovis, M. avium, M. intracellular, M. kansasii and M. chelonei. The organisms were inoculated in sterile reaction vials containing liquid 7H9 medium, 10% ADC enrichment and 1.0 uCi of one of the (1-14C) fatty acids (butyric, hexánoic, octanoic, decanoic, lauric, myristic, palmitic, stearic, oleic, linoleic, linolenic). Vials were incubated at 37°C and the 14CO2 envolved was measured daily for 3 days with a Bactec R-301 instrument. Although each individual organism displayed a different pattern of fatty oxidation, these patterns were not distinctive enough for identification of the organism. No combination of fatty acids nor preferential oxidation of long chain or of short chain fatty acids were able to separate susceptible from resistant organisms. Further investigation with a larger number of drug susceptible mycobacteria including assimilation studies and oxidation of other substrates may be required to achieve a distinction between drug-susceptible and drug-resistant mycobacteria.

Radiometric studies of mycobacterium tuberculosis

An in vitro assay system that included automated radiometric quantification of 14CO2 released as a result of oxidation of 14C- substrates was applied for studying the metabolic activity of M. tuberculosis under various experimental conditions. These experiments included the study of a) mtabolic pathways, b) detection times for various inoculum sizes, c) effect of filtration on reproducibility of results, d) influence of stress environment e) minimal inhibitory concentrations for isoniazid, streptomycin, ethambutol and rifampin, and f) generation times of M. tuberculosis and M. bovis. These organisms were found to metabolize 14C-for-mate, (U-14C) acetate, (U-14C) glycerol, (1-14C) palmitic acid, 1-14C) lauric acid, (U-14C) L-malic acid, (U-14C) D-glucose, and (U-14C) D-glucose, but not (1-14C) L-glucose, (U-14C) glycine, or (U-14C) pyruvate to 14CO2. By using either 14C-for-mate, (1-14C) palmitic acid, or (1-14C) lauric acid, 10(7) organisms/vial could be detected within 24 48 hours and as few as 10 organisms/vial within 16-20 days. Reproducible results could be obtained without filtering the bacterial suspension, provided that the organisms were grown in liquid 7H9 medium with 0.05% polysorbate 80 and homogenized prior to the study. Drugs that block protein synthesis were found to have lower minimal inhibitory concentrations with the radiometric method when compared to the conventional agar dilution method. The mean generation time obtained for M. bovis and different strains of M. tuberculosis with various substrates was 9 ± 1 hours.

Radiometric detection of metabolic activity of Paracoccidioides brasiliensis and its susceptibility to amphotericin B and diethylstilbestrol

Paracoccidioidomycosis (South American blastomycosis) is a systemic disease, strikingly more frequent in males, caused by the dimorphic fungus Paracoccidioides brasiliensis. A radiometric assay system has been applied to study the metabolic activity and the effect of drugs on this fungus "in vitro". The Y form of the yeast, grown in liquid Sabouraud medium was inoculated into sterile reaction vials containing the 6B aerobic medium along with 2.0 μCi of 14C-substrates. Control vials, prepared in the same way, contained autoclaved fungi. To study the effects of amphotericin B (AB) (0.1 and 10 μg/ml) and diethylstilbestrol (DSB) (1.0, 5.0 and 10 μg/ml) extra controls with live fungi and no drug were used. All vials were incubated at 35°C and metabolism measured daily with a Bactec instrument. 14CO2 production by P. brasiliensis was slow and could be followed for as long as 50 days. AB at 10mg/ml and DSB at 5 μg/ml inhibited the metabolism and had a cidal effect on this fungus. The results with DSB might explain the low incidence of the disease in females. This technique shows promise for studying metabolic pathways, investi gating more convenient 14C-substrates to expedite radiometric detection and for monitoring the effects of other drugs and factors on the metabolism of P. brasiliensis "in vitro".

Radiometric studies on the oxidation of (U-14C) L-amino acids by drug-susceptible and drug-resistant mycobacteria

A radiometric assay system has been used to study oxidation patterns of (U-14C) L-amino acids by drug-susceptible and drug-resistant mycobacteria. Drug-susceptible M. tuberculosis (H37Rv TMC 102 and Erdman) along with the drug-resistant organism M. tuberculosis (H37 Rv TMC 303), M. bovis, M. avium, M. intracellulare, M. kansasii and M. chelonei were used. The organisms were inoculated into a sterile reaction system with liquid 7H9 medium and one of the (U-14C) L-amino acids. Each organism displayed a different pattern of amino acid oxidation, but these patterns were not distinctive enough for identification of the organism. Complex amino acids such as proline, phenylalanine and tyrosine were of no use in identification of mycobacteria, since virtually all organisms failed to oxidize them. There was no combination of substrates able to separate susceptible from resistant organisms.