| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
AATSR Frequently Asked QuestionsPDF VersionDocument Ref: AEP.REP.001 Issue Date: 20 December 2005 Issue: 1
TABLE OF CONTENTS Q1. What does AATSR stand for? Q2. What is AATSR and what does it do? Q4. What orbit does ENVISAT use? Q5. What can AATSR data be used for? Q6. If this FAQ does not answer my question, what should I do? Q7. Where can I find a ‘user guide’ to receiving and utilising AATSR data? Q8. How can I be kept up to date with events that might affect my use of AATSR data? Q9. What are the main differences between ATSR-1, ATSR-2 and AATSR? Q10. How can I order AATSR data? Q13. Where can I find articles/papers on the exploitation of (A)ATSR data? Q14. What tools are available for reading AATSR products? 4.1 Instrument Characteristics Q15. What is the range of nadir and forward view incidence angles applicable to (A)ATSR? Q16. Where can I obtain the full set of AATSR visible channel spectral response functions? Q17. How can I obtain information about any interruptions to AATSR data acquisition? Q18. Which version of the AATSR IPF was used to produce my data? Q19. What differences will there be between data processed in Near-Real-Time and Off-Line? Q22. How do I interpret the information in the filename of my AATSR data? Q23. Why does my AATSR image appear to contain negative temperatures? Q35. Why do the following problems occur with the AATSR LST retrieval? Q36. Where can I find examples of AATSR Level 3 products? APPENDIX A. Interpolation of pixel geolocation in AATSR Full resolution products
AMENDMENT POLICY This document shall be amended by releasing a new edition of the document in its entirety. The Amendment Record Sheet below records the history and issue status of this document. AMENDMENT RECORD SHEET
IntroductionPurpose and ScopeThe purpose of this FAQ is to act as a repository for all of the common questions raised concerning the AATSR project. It includes questions raised with the AATSR Expert Support Laboratory, questions reported to the Quality Working Group and questions of particular interest to new users. The information contained within this document is supplementary to that contained within the AATSR Handbook, available from http://envisat.esa.int/dataproducts/aatsr/ .The questions fall into four basic categories:
Structure of the DocumentAfter this introduction, the document is divided into a number of major sections that are briefly described below: 2 GENERAL QUESTIONS This section addresses general background information relating to the AATSR instrument and the ENVISAT satellite. 3 COMMON QUESTIONS FROM NEW USERS This section addresses issues regularly raised by new users of AATSR data. 4 QUESTIONS ABOUT THE INSTRUMENT This section addresses issues concerning the AATSR instrument and is divided into sub-sections on instrument characteristics and instrument operations. 5 QUESTIONS ON DATA PROCESSING AND PRODUCTS This section addresses issues associated with the processing applied to AATSR data and the format and contents of the resulting products. APPENDIX A Interpolation of pixel geolocation in AATSR Full resolution products This section sets out the method of interpolating the latitude and longitude values from the AATSR full resolution product geolocation Annotation Data Set, to give the co-ordinates of the image pixels in these products. GENERAL QUESTIONSAATSR stands for ‘Advanced Along-Track Scanning Radiometer’. Q What is AATSR and what does it do? AATSR is one of the Announcement of Opportunity (AO) instruments on board the European Space Agency (ESA) satellite ENVISAT. It is the most recent in a series of instruments designed primarily to measure Sea Surface Temperature (SST), following on from ATSR-1 and ATSR-2 on board ERS-1 and ERS-2. AATSR data have a resolution of 1 km at nadir, and are derived from measurements of reflected and emitted radiation taken at the following wavelengths: 0.55 µm, 0.66 µm, 0.87 µm, 1.6 µm, 3.7 µm, 11 µm and 12 µm. Special features of the AATSR instrument include its use of a conical scan to give a dual-view of the Earth’s surface, on-board calibration targets and use of mechanical coolers to maintain the thermal environment necessary for optimal operation of the infrared detectors. The AATSR programme is funded jointly by the UK Department for Environment, Food and Rural Affairs , the Australian Department of Industry, Science and Resources (DISR), and the UK Natural Environment Research Council (NERC).ENVISAT was launched in March 2002 by ESA, and is an advanced polar-orbiting Earth observation satellite providing measurements of the atmosphere, ocean, land, and ice (see http://envisat.esa.int/ ).ENVISAT continues the work of the ERS satellites, and its data supports earth science research and allows monitoring of the evolution of environmental and climatic changes. Q What orbit does ENVISAT use? ENVISAT flies in a sun-synchronous polar orbit of about 800-km altitude. The repeat cycle of the reference orbit is 35 days. Q What can AATSR data be used for? AATSR primarily measures SST to the high levels of accuracy and precision required for monitoring and detecting climate change. Together with its predecessors, ATSR-1 and ATSR-2, AATSR will establish a unique SST data set spanning at least 15 years, supporting oceanographic and climate research. AATSR data can also be used for a number of land surface, cryosphere and atmospheric applications. A link to further information on the applications of (A)ATSR data will be added in the next issue of this FAQ. Q If this FAQ does not answer my question, what should I do? All queries regarding ENVISAT and AATSR should be directed to the ESA Earth Observation Help Desk Team ( EOHelp@esa.int ) in the first instance.COMMON QUESTIONS FROM NEW USERSQ Where can I find a ‘user guide’ to receiving and utilising AATSR data? The AATSR Handbook can be downloaded from http://envisat.esa.int/dataproducts/ . Section 1 of the Handbook is designed as a Product User Guide, intended to help users familiarise themselves with AATSR and "get started" with using AATSR data. This is supported by later sections containing detailed information on the instrument, the products it produces and the algorithms used to generate them.Q How can I be kept up to date with events that might affect my use of AATSR data? AATSR Cyclic Reports are made available by ESRIN to keep the AATSR community informed of any modifications to the processor, updates of auxiliary products, instrument anomalies, the status of data acquisition and processing, and the status of the calibration, validation, and quality control activities. For a full list of the reports, see http://earth.esa.int/pcs/envisat/aatsr/reports/cyclic .The AATSR roject team also maintain an "AATSR Contacts List" and issue Bulletins to users several times per year. These Bulletins are prepared by the VEGA Group PLC, on behalf of ESA and the AATSR Quality Working Group (QWG). The objectives of the bulletins are to keep users informed of instrument operations and data processing issues that may affect the exploitation of the data and which may not be reported in the Cyclic Reports. News and announcements from the AATSR programme are also included. To be added to the contacts list and to receive the AATSR Bulletins, contact EOHelp ( eohelp@esa.int ).ESA also issue ENVISAT news items at http://envisat.esa.int/news/ and http://envisat.esa.int .Q What are the main differences between ATSR-1, ATSR-2 and AATSR? ATSR-1 measured in the infra-red at 1.6 µm, 3.7 µm, 11 µm and 12.0 µm and had no visible channels. ATSR-2 and AATSR include the same four infra-red channels of ATSR-1 and three additional channels at 0.55 µm, 0.67 µm and 0.87 µm. The ATSR-2 instrument is largely the same as ATSR-1, differing only in the inclusion of the extra channels and an on-board visible calibration system. AATSR is functionally largely similar to ATSR-2, but components have been redesigned to match the new platform environment of ENVISAT. The details of the scan, calibration systems, spatial resolution and swath have been kept as close as possible to the earlier instruments to ensure data continuity. The major advantage of AATSR over ATSR-2 is in the downlinking of data. ATSR-2 had restrictions on the amount of data that could be downlinked, due to platform constraints, which particularly affected the visible channels, whereas AATSR can send down data from all the channels continuously at full 12-bit radiometric resolution. This also simplifies the data processing for AATSR, since the processor does not have to cope with the wide range of instrument data formats that could be used for ATSR-2. The specific differences between the ATSR-2 and AATSR source packets are discussed in section 2.5.2 of the AATSR handbook (see Q7). To order AATSR data, please visit the How to apply page. A summary of the relationship between the AATSR and ATSR products is given in Section 2.2.3 of the AATSR Handbook (see Q7). Defra, NERC and ESA are currently supporting the development of an (A)ATSR-series archive in which data from ATSR-1, -2 and AATSR are reprocessed and stored together in an identical format (the ENVISAT product format). This is expected to become operational progressively from 2006 onwards. Further information on access to this archive will be added in future issues of this FAQ. In the meantime, all requests for ATSR-1 and ATXR-2 data should be sent to ESA in the first instance ( EOHelp@esa.int) .Q Where can I find articles/papers on the exploitation of (A)ATSR data? MERIS and (A)ATSR Workshop (26 – 30 September 2005) http://envisat.esa.int/workshops/meris_aatsr2005/ENVISAT/ERS Symposium (6-10 September 2004) http://earth.esa.int/salzburg04/ENVISAT Meris and AATSR Validation Team Workshop - MAVT-2003 (20 - 24 October 2003) http://envisat.esa.int/workshops/mavt_2003/Envisat Validation Workshop (9-13 December 2002) http://envisat.esa.int/pub/ESA_DOC/envisat_val_1202/proceedings/Envisat Calibration Review (9-13 September 2002) http://envisat.esa.int/calval/proceedings/ERS-Envisat Symposium "Looking down to Earth in the New Millennium" (16-20 October 2000) http://earth.esa.int/pub/ESA_DOC/gothenburg/start.pdfATSR Project Web Site – Documentation Section http://www.atsr.rl.ac.uk/documentation/index.shtmlESA Data User Element (DUE) Projects http://dup.esrin.esa.it/index.aspESA EO PI Projects http://eopi.esa.int/esa/esaAATSR PI Team at Leicester University – Publications List http://www.leos.le.ac.uk/aatsr/ENVISAT Library http://envisat.esa.int/services/esa_doc/doc_env.htmlN.B. An AATSR Reference Documents web site is due to be created in 2006. Further details will be provided in subsequent issues of this FAQ.Q What tools are available for reading AATSR products? Most users write their own routines in C, IDL or other languages to read and process AATSR data. The detailed AATSR product specifications (PO-RS-MDA-GS-2009, Volume 7 available at http://envisat.esa.int/services/esa_doc/doc_env.html ) will assist with this task.The EnviView tool provided by ESA (see below) can also be used to convert ENVISAT data into hdf file format. This allows the data to be automatically read by any third-party software supporting this format. The providers of a number of COTS image processing packages such as ERDAS, ENVI and Geomatica also planned to extend their products to directly support the AATSR product format. The individual packages/providers should be consulted for details of each tool’s capabilities. EnviView is a free application that allows Envisat data users to open any Envisat data file and examine its contents. It provides simple visualisation capabilities, and allows data to be exported to HDF for use in other software packages. The Basic ERS & Envisat (A)ATSR and Meris Toolbox (BEAM) is a collection of executable tools and an application programming interface (API) which has been developed to facilitate the utilisation, viewing and processing of ESA MERIS, (A)ATSR and ASAR data. 1 tape_reader.ksh is a Korn shell script which may be used to extract Envisat products supplied on exabyte tapes. All three of the above tools are available from: http://www.envisat.esa.int/services/tools_table.html .There are a number of important concepts associated with the AATSR products that should be understood before trying to read and interpret the data. These are discussed in Section 2.12.1 of the AATSR handbook (see Q7). QUESTIONS ABOUT THE INSTRUMENTQ What is the range of nadir and forward view incidence angles applicable to (A)ATSR? The range of incidence angles between the swath centre and edge is an input to the SST coefficients generation process. The values that are used in the current Radiative Transfer Modeling code, in degrees, are: nadir 0.0 (swath centre) 21.732 (swath edge) forward 55.587 (swath centre) 53.009 (swath edge) Strictly the value depends on which instrument is being considered (because the cone angles differ slightly), but in practice the differences are small. Q Where can I obtain the full set of AATSR visible channel spectral response functions? There is a link to these data from Section 3.2.1.1.2 of the AATSR handbook (see Q7). Q How can I obtain information about any interruptions to AATSR data acquisition? An Instrument Status History log is maintained by the AATSR Flight Operations Support Team and is available from the AATSR Operations web site http://www.aatsrops.rl.ac.uk . The Weekly Reports published on this site also contain a section on Forthcoming Activities listing forthcoming outgassings, Black Body cross-over tests, Envisat orbital maneuvers etc. as and when the FOS team are made aware of them (see Weekly Reports\Latest).Instrument availability is also logged by ESA at http://envisat.esa.int/instruments/availability/ and both instrument and data unavailability is reported in the AATSR Cyclic Reports http://earth.esa.int/pcs/envisat/aatsr/reports/cyclic .Background information on issues relating to AATSR Flight Operations can be found in Sections 1.1.3.2 and 3.2.2.2.1 of the AATSR handbook (see Q7). Questions on data processing and productsData ProcessingQ Which version of the AATSR IPF was used to produce my data? The AATSR IPF version number is a field within the Main Product Header of the data. Volume 5 of the ENVISAT Product Specification (Product Structures) details the contents of the product headers. This document is due to be added to the archive at http://envisat.esa.int/services/esa_doc/doc_env.html shortly. Q What differences will there be between data processed in Near-Real-Time and Off-Line? NRT and off-line products will be identical, but for their level of consolidation and the auxiliary files used in their production (see section 2.10.1 of the AATSR Handbook (see Q7)). Updates to most AATSR auxiliary files will be infrequent, with the exception of the ATS_VC1_AX file containing visible channel calibration coefficients. The version of this file used for NRT processing will not contain the exact coefficients for that orbit. A file from one or two days previous to the acquisition of the orbit in question will be used instead. Under normal circumstances VC1 files used for off-line processing should be no more than 1 day old. For more information on this, see Section 2.11.2 of the AATSR Handbook (see Q7). NRT data are also processed with predicted rather than restituted Orbit State Vectors, but this does not have a significant effect on geolocation/collocation for AATSR.Disclaimers addressing issues affecting AATSR product quality are published at http://envisat.esa.int/dataproducts/availability/ .The auxiliary file of SST coefficients used by the Level 2 processor includes distinct coefficients for each of 38 across-track distances. The process by which these 38 sets of coefficients are generated involves first deriving two sets of coefficients, for the swath centre and swath edge, and then computing intermediate values by linear interpolation with respect to air mass. (Here the air mass refers to the normalised length of the atmospheric path traversed by the line of sight - essentially the secant of the angle of incidence at the pixel.) The air mass at the swath centre is unity (that is the normalisation), but that at the edge depends on factors such as the cone angle of the AATSR scan and the satellite height. A study of the sensitivity of the edge of swath air mass to changes in the cone angle and the satellite height has been conducted. In general, the results are very insensitive to cone angle; satellite height is a slightly more significant parameter, but as it varies around the orbit there is not much that can be done about it other than to use a suitable average value. Further details can be obtained from the AATSR Expert Support Laboratory team at RAL (initial contact to be made through EOHelp@esa.int) . ProductsQ How do I interpret the information in the filename of my AATSR data? This is explained in Section 2.2.1 of the AATSR handbook (see Q7). Volume 5 of the ENVISAT Product Specification (Product Structures) also addresses this and other issues related to ENVISAT product conventions ( http://envisat.esa.int/services/esa_doc/doc_env.html ).Q Why does my AATSR image appear to contain negative temperatures? Small negative values (-1 to -8) are 'exception values' used to flag errors and other special conditions in the data. These are listed in Sections 2.6.2.1.1 of the AATSR Handbook (see Q7) for L1b data. Most topographic correction tie points over sea are set to zero (as the sea is at sea level). However, exception values (notably –999999) can arise in the current algorithm if no pixel at a given across-track position regrids to a particular topographic correction tie point. This occurs because there is no cosmetic fill in the topographic correction and the array is initialised to –999999 before calculations commence. The orbit number in the product file name is that in which the first data contributing to the product falls. This is not a problem for NRT data, as NRT products are not expected to begin exactly at ANX. However, in the case of consolidated data (orbits defined from ANX to ANX), if the first data slightly precedes the ascending node, then the orbit number of the file is one lower than you might expect. Taking the following file as an example: ATS_TOA_1PPLRA20040605_001448_000062202027_00259_11836_1706.N1.gz The start of data (00:14:48 on 2004 June 5) is very close to the ascending node at the start of absolute orbit 11837. Therefore, although the product predominantly represents orbit 11837, the start of data falls (just) in the preceding orbit. There is no unambiguous answer to this. Both the AST and GST products are derived from the same set of L1b data, just using different algorithms, so they are, in effect, both ‘right’. The real issue is the exact definition of the 'number of pixels' field in the AST product.
The field of the AST record designated 'Number of pixels in dual view average' contains the smallest of the numbers of pixels contributing to the averaged brightness temperatures (in the 11 and 12 micron channels and in the nadir and forward views) that were used in the SST retrieval. In an ideal world the same number of pixels would contribute to each of the averaged brightness temperatures in the cell, and this number would appear in the field, but in practice the numbers may differ for different channel/view combinations. In particular, although the numbers of 11 and 12 micron values contributing are likely to be the same in each view, the nadir and forward views may differ because of e.g. missing scans. As the number of pixels is supplied as an indicator of SST quality, the smallest is chosen. Similarly the field 'Number of pixels in nadir-only average' would contain the smaller of the numbers of pixels in the 11 and 12 micron averages where these differ. The 3.7 micron channel numbers are ignored because this channel may or may not contribute to the SST retrieval. The units for the visible channel bandwidth in the IODD, in the ENVISAT Product Specifications and in EnviView are given as microns, but the actual numerical value is in nm. The units for the integrated solar irradiance values are microwatts/cm2. The information required for the interpolation of pixel geolocation can be found in the AATSR Technical note entitled "Interpolation of Pixel Geolocation in AATSR Full Resolution Products" referenced in section 2.12.1.3.1 of the AATSR handbook (see Q7). The contents of this note are also reproduced in Appendix A of this FAQ. The pixel co-ordinate resulting from the interpolation (as described in the above-referenced technical note) corresponds to the lower left corner of the pixel. In the across-track direction the co-ordinate corresponds to the left-hand edge, while in the along track direction, the sampling of the tie points is such that the tie points correspond to the lower edge of the pixel (i.e. the edge corresponding to the lower Y co-ordinate).
The short field names, such as pix_nad, used in EnviView and the Formats sections of the AATSR handbook (see Q7) have not been assigned uniquely; the same names have been used, in different products and MDS, for different things. The data definition document for the AATSR processing software defines a unique parameter ID for each AATSR product field. The unique parameter names of the values referred to are as follows:
The Level 2 processor generates the averaged brightness temperatures by averaging all of the valid image pixels from the Level 1b product that fall in each 10 arc minute cell. Image pixels in the Level 1b product are classifiable into three mutually exclusive types: NATURAL, COSMETIC, and UNFILLED; every pixel is one of these types, and this classification is purely geometrical, and independent of whether or not the pixel contains valid data in any channel. Essentially, NATURAL pixels are those to which a measured instrument pixel has been mapped during regridding. (The majority of image pixels should be of this type.) A pixel to which no instrument pixel maps during regridding is either COSMETIC or UNFILLED. A COSMETIC pixel has been filled by the cosmetic fill algorithm; it will have at least one natural pixel as neighbour. Otherwise the pixel is UNFILLED. In the BT/TOA MDS, the quantity [AST-MDS15-7] (pix_nad) is the total number of filled (i.e. NATURAL or COSMETIC) pixels in the 10 arc minute sub-cell, regardless of whether or not the brightness temperatures are valid. Similarly [AST-MDS15-8] (pix_ss_nad) is the number of filled pixels that are over sea, and [AST-MDS15-9] (clpix_ss_nad) is the fraction of these that have been flagged as cloudy. (Since the cloud clearing algorithms depend upon the pixel data, the last of these quantities has some dependence on whether or not the pixel data is valid.) [AST-MDS15-7] (pix_nad) should equal [AST-MDS15-8] (pix_ss_nad) if the sub-cell is entirely sea, but not if it is intersected by coastline. The quantities [ECM-MDS1-8] and [AST-MDS3-8], both designated pix_nad in EnviView, are identical; that in the Meteo product, [ECM-MDS1-8], is a copy of [AST-MDS3-8]. They represent the number of valid pixels that contribute to the nadir SST in the cell. This number cannot exceed, but may be less than, pix_ss_nad, since some of the filled sea pixels may contain invalid data. Its detailed definition is the smaller of the number of valid pixels that contribute to the averaged brightness temperatures at 11 and 12 microns; this definition reflects the theoretical possibility that a given pixel might have a valid BT in one channel but not in the other, although this is unlikely over clear sea.
To be addressed in the next issue of this FAQ. At night, noise in the visible channels cannot be distinguished from pixel exception values, and can therefore be masked into the confidence flags. This is not strictly an error in the AATSR processing scheme, but rather an unexpected result of adopting the ATSR-2 processing scheme for AATSR (ATSR-2 visible channel data were not available at night). The AATSR pixel exception values are small negative numbers (ranging from -1 to -8). At night, the noise in the visible channels can mimic these exception values, so that the interpretation of negative visible channel reflectances is ambiguous. Also, because the confidence flags are derived from the union of the exception values across all channels, this can cause the confidence flags to be wrongly set. Note that the latter effect affects any user who is using the confidence flags to filter the data, regardless of whether he is interested in the visible channels or not. Users are advised to take note of this effect and to take into account these effects if using the confidence flags or visible channel reflectances in night-time data. This is also due to the presence of noise in the visible channel data at night. NDVI values are calculated whenever visible channel data are present, and as visible channels are present in the AATSR telemetry at all times, the NDVI is calculated for both day and night data. There is no trap for night-time data. However, at night, the visible channels are dominated by noise, making these NDVI values essentially meaningless. Users are therefore advised to disregard the values in the AATSR NDVI field at night. Q Why do the following problems occur with the AATSR LST retrieval? Poor cloud clearing The current AATSR cloud clearing scheme is optimised for use over oceans and its performance cannot be guaranteed over land. Land pixels may be wrongly flagged as cloud and when a pixel is flagged as cloud no LST retrieval is applied. Pixels flagged as cloud will be set to Cloud Top Temperature (currently the 11µm BT) and Cloud Top Height (currently zero). Some AATSR images may contain pixels flagged as cloud, even when ground measurements show that the conditions were clear. It was always appreciated that the imperfections of cloud detection over land would be a limitation of the existing LST implementation, and ultimately a revised cloud clearing scheme optimised for performance over land needs to be defined.
A "blocky" appearance with discontinuities at surface classification cell boundaries. No spatial smoothing of the LST has been included in the current implementation of the LST algorithm. Therefore the outline of the 0.5 degree cells used in the underlying land surface type classification map may be visible in the retrieved LST field. While there may be merit in using an interpolation scheme to smooth such discontinuities at the cell boundaries, it might also distort the land surface temperatures in an unquantifiable way so has not been included in the current LST algorithm. No Lake Surface Temperature retrieval The ATBD for the AATSR LST algorithm indicates that temperature retrieval coefficients are available for use over lakes. However, the geophysical parameter returned over lakes is Sea Surface Temperature. The land surface type file contains a class for lakes and lake surface temperature coefficients are included in the coefficients set. However, most large and medium-sized lakes are flagged as sea in the Level 1b land/sea mask so SSTs are calculated over these regions and the lake surface temperature calculation is never invoked. This may be modified in the future. Regardless of the above problem, retrieving lake surface temperatures over medium-sized lakes via the LST algorithm will always be problematic owing to the fact that the land surface type file is constructed at 0.5 degree resolution. This means that lakes that lie on the boundary of two or more 0.5 degree cells may not be properly identified. The margins of the lakes will also not be correctly identified at this resolution. Consequently it is currently recommended that the lake coefficients are applied manually to the L1b data – Dr Fred Prata can provide further information on how to do this (initial contact through EOHelp@esa.int). Problems with the calculation of LST in coastal regions Clear land pixels over validation sites close to the coast have been seen to contain the 11 micron BT, rather than LST. This occurs where pixels are identified as land in the Level 1b land/sea mask, but are flagged as sea in the Level 2 LST surface type file used by the LST algorithm. An AATSR pixel is identified as land or sea during Level 1b processing with reference to the AATSR land/sea mask. This is an auxiliary file classifying points on the Earth’s surface as land or sea at 1 km resolution. When a Level 1b product is presented for processing to Level 2, the SST retrieval algorithms are invoked for pixels flagged as sea and the LST/NDVI algorithms are invoked for pixels flagged as land. As part of the Level 2 LST algorithm land pixels are then further categorised into 13 biomes with reference to a land surface type file. The Level 2 land surface type file is provided at 0.5 degree resolution. All of the AATSR (1 km) pixels that fall within each 0.5 degree cell are deemed to be either land (including lakes) or sea according to the surface type associated with these cells. However, it is possible for some 0.5 degree cells that cross the boundary between land and sea to be classified as sea when they also contain land. For these cells no surface type or biome is defined for the land pixels that they contain. In this instance the LST algorithm is unable to make a valid LST retrieval; it therefore fills the pixel with the 11 micron BT and sets the pixel status to invalid in the confidence word. [Note that bits 0 and 2 in the confidence word are currently labeled "Nadir-only SST is valid" and "Dual-view SST is valid" respectively, but they are also applicable to LST. They should really be re-named "Nadir field is valid" and "Combined field is valid".] No average LST retrieval in the AST product Although a processing scheme has been defined, the retrieval of average LST has not yet been prototyped and tested, and is therefore not yet ready to become operational. N.B. An initiative is currently underway to address the following issues with the LST retrieval:
These improvements are expected to become operational in Spring 2006. Q Where can I find examples of AATSR Level 3 products? Example L3 products can be found at http://envisat.esa.int/level3/ .This site currently contains global SST Level 3 products from AATSR from September 2002 to January 2005, generated at two spatial resolutions, 10 arcminutes and 30 arcminutes. Further examples may be added in the future.
APPENDIX A: Interpolation of pixel geolocation in AATSR Full resolution products Introduction The AATSR full resolution products are the Gridded Brightness Temperature and Reflectance (GBTR) product ATS_TOA_1P, and the Gridded Surface Temperature (GST) product ATS_NR__2P. In both of these products the information from which the geolocation of the image pixels can be determined is given in an Annotation Data Set (ADS), the Grid Pixel Latitude and Longitude and Topographic Corrections ADS. This note sets out the method of interpolating the latitude and longitude values from the ADS to give the co-ordinates of the image pixels in these products. A.2 Notation The notation used in this document is summarised below.
A.3 X and Y Co-ordinates The Cartesian X and Y co-ordinates used for AATSR are instrument related and are measured in the across-track and along-track directions respectively. The X axis is at right angles to the satellite ground track, and the X co-ordinate increases towards the right (looking forwards from above). The Y axis is parallel to the satellite direction of motion, and the Y co-ordinate increases in the forward direction.
A.4 Product Structure The Measurement Data Sets (MDS) of each full resolution (GBTR or GST) product consist of a series of pixels arranged on a rectangular grid. Each record of the product contains 512 pixels in a line across track, at a sampling interval of 1 km. Let the records in the MDS be numbered by i (i = 0, 1, 2, …) and let the pixels in each record be indexed by j (j = 0, ... 511). The pixels in each record are arranged in order of increasing X co-ordinate, and so the X co-ordinate of the pixel indexed by j is
(As in the case of ATSR-2, a pixel can be thought of as a small quadrilateral, and its nominal co-ordinates are those of its lower left hand corner, and so this co-ordinate refers to the left-hand edge of the pixel.) The along-track co-ordinate Y(i) of the pixels on record i is specified on the record, in the field designated ‘image scan y co-ordinate’, in units of metres. The increment Y(i + 1) - Y(i) is approximately equal to 1 km, but varies around the orbit owing to the equal time interval sampling adopted in AATSR products. The geolocation information can be found in the Grid Pixel Latitude and Longitude ADS, and comprises latitude and longitude at a series of tie points spaced by 25 km across track and by 1 product granule (32 records, or approximately 32 km) along-track. Each record of the ADS contains latitude and longitude values at 23 tie points whose X co-ordinates are X = -275, -250, … +275 km. There is one ADS record for every 32 MDS records, such that the first ADS record corresponds to the first MDS record (that is, they have the same Y co-ordinate). For simplicity we adopt the convention that records in an MDS or ADS are numbered so that the first record in the data set is record zero (although EnviView numbers records from 1). Thus if k (k = 0, 1, 2, …) is the number of the record in the ADS, the ADS record identified by k corresponds to the MDS record numbered
The y co-ordinates of these two records should be identical. A.5 Bilinear Interpolation of Geolocation Information Given the (geodetic) latitudes and longitudes of the tie points from the Grid Pixel Latitude and Longitude ADS, the latitude and longitude of each image pixel may be derived by interpolation between these tie points in two dimensions. In the case of longitude account must be taken of the possibility that the 180 degree meridian intersects the image row. Suppose we wish to find the latitude and longitude of a point whose image co-ordinates are X and Y. Then the index of the tie point to the left of the point (X, Y) is
where int[x] represents the integer part of x. Calculate
The point with along-track co-ordinate Y lies between the tie rows indexed by ig, ig + 1, where
and where Y(ig) is the image scan y co-ordinate of ADS record ig. The interpolation weight for use in the Y direction is then
In the particular case that X and Y are the co-ordinates of the pixel indexed by i, j, then substitution of Equation (1) in Equations (3) and (4) respectively gives
Similarly for the Y co-ordinate, from Equation (2) with k = ig,
The geocentric latitude of the pixel is then calculated by interpolation as follows:
where
and where the tie point latitudes Longitude is treated similarly unless the 180 degree meridian is present:
where
and the tie point longitudes However, if the position to be interpolated is close to longitude 180 degrees, special treatment is required because the longitude values in the ADS are specified in the range -180 to +180 degrees. Thus if longitude 180 degrees falls between the tie points it will give rise to a discontinuity in the longitude. This is the case if
where
In this case 360.0 should be added to each of the tie point
longitudes
[1] At the time of writing the AATSR PI Team report that the BEAM IDL API produces geolocation co-ordinates that are half a pixel out in the along-track direction. It is believed that this has been corrected in BEAM itself, but not yet in the API. Keywords: ESA European Space Agency - Agence spatiale europeenne, observation de la terre, earth observation, satellite remote sensing, teledetection, geophysique, altimetrie, radar, chimique atmospherique, geophysics, altimetry, radar, atmospheric chemistry |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Copyright 2000 - European Space Agency. All rights reserved. | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||