For highest possible spatial resolution, multiple orbit passes are combined. There is thus a tradeoff between temporal and spatial resolution. Generally orbit passes of a given point on the Earth occur at nearly the same time of day.
Investigator: | Dr. David G. Long |
Title: | Director,
BYU Center for Remote Sensing Professor, Department of Electrical & Computer Engineering |
Address: | 459 Clyde Building Brigham Young University Provo, UT 84602 |
Email: | long@ee.byu.edu |
Fax: | (801)378-6586 |
...the AMI ground processors ... take the digital radar echo samples and produce a radar image with the required geometrical and radiometric properties. In such imagery, the different gray-tones represent different levels of radar backscattering, expressed by the dimensionless quantity sigma0... In a later processing step, sigma0 values are related to geophysical quantities such as wind speed over the ocean, ice type, soil moisture, agricultural crop type, stage of growth, forest type, etc.
The AMI is a multimode radar operating at a frequency of 5.3 GHz (C-band), using vertically polarized antennas for both transmission and reception...In the wind mode, the AMI is configured as a wind scatterometer and provides three radar images of the ocean surface with a spatial resolution of 50 km and a swath width of 500 km. The three images are acquired by three different antennas: one, the mid-beam, looking to the right side of the satellite, perpendicular to the ERS-1 ground track, one looking forward at 45 degrees azimuth projection angle, with respect to the mid-beam, and one looking backwards at 45 degrees azimuth projection angle with respect to mid-beam.
SIR, average files: 8.9 km pixel gridA field in the header also identifies the resolution.
Gridded files: 44.5 km pixel grid
Antarctic, Arctic files: polar stereographic projectionA field in the header also identifies the projection.
All others: Lambert Equal Area projection
ers-T-reg.rcn
T | image type | x = longitude, y = latitude |
reg | region | Ala = Alaska, Ant = Antarctica, ... |
rcn | reconstruction technique | sir = SIR, grd = gridded |
Two other types of auxiliary files for each region at each spatial resolution contain topography and land mask information. The naming scheme for these files is:
ers-reg.rcn.info
reg | region | Ala = Alaska, Ant = Antarctica, ... |
rcn | reconstruction technique | sir = SIR, grd = gridded |
info | type | topo = topography, lmask = land mask |
ersp-T-regyy-dd1-dd2.rcn[.lmsk]
p | spacecraft |
T | image type
| reg
| geographical | region - see Spatial Coverage above
yy |
two-digit year, always 78
| dd1
| three-digit day of year, start of imaging
| dd2
| three-digit day of year, end of imaging
| rcn |
reconstruction |
technique |
Each file also has header information. The program xv printed the following sample output as it displayed the file ers1-a-Ala92-001-006.sir.lmsk:
SIR file header: 'ers1-a-Ala92-001-006.sir.lmsk' Title: 'SIR image of alaska' Sensor: 'ERS-1/2' Type: 'A image (ers1-a-Ala92-001-006.sir)' Tag: '(c) 2001 BYU MERS Laboratory' Creator: 'BYU MERS:ers_meta_sir v4.0 Ai=-20.00 Bi=-0.130 Bacc= 50.0 RefInc= 40.00 It=27' Created: '04:26:00 05/18/01' Size: 410 x 320 Total:131200 Offset: -33 Scale: 1000 Year: 1992 JD range: 1-6 Region Number: 2 Type: 1 Form: 2 Polarization: 2 Frequency: 5.300000 MHz Datatype: 2 Headers: 1 Ver:31 Nodata: -33.000000 Vmin: -32.000000 Vmax: 0.000000 Lambert form: (local radius) Center point: -155.000000 , 61.500000 Lon, Lat scale: 8.900000 , 8.900000 (km/pix) Lower-Left Corner: -1800.000000 , -1300.000000 Image Min, Max: -32.000000 , 0.000000 Greyscale conversion range: Min: -32.000000, Max:0.000000
The EOSDIS Glossary describes data granularity generally as it applies to the IMS.
A SIR format file consists of one or more 512-byte headers followed by the image data and additional zero padding to insure that the file is a multiple of 512 bytes long. The file header record contains all of the information required to read the remainder of the file and the map projection information required to map pixels to lat/lon on the Earth surface. The image pixel values generally represent floating point values and may be stored in one of three ways. The primary way is as 2 byte integers (with the high order byte first), though the pixels may be stored as single bytes or IEEE floating point values. Scale factors are stored in the header to convert the integer or byte pixel values to native floating point units.
The image is stored in row-scanned (left to right) order from the lower left corner (the origin of the image) up through the upper right corner. By default, the location of a pixel is identified with its lower-left corner. The origin pixel (1,1) is the lower left corner of the image. The array index n of the (i,j)th pixel where i is horizontal and j is vertical is given by
n = (j - 1) × Nx + iwhere Nx is the horizontal dimension of the image. The last pixel stored in the file is at (Nx, Ny).
The sir file header contains various numerical values and strings which describe the image contents. For example, the value for a no-data flag is set in the header as well as a nominal display range and the minimum and maximum representable value. Optional secondary header records (512 bytes) can be used to store additional, non-standard information.
The standard SIR file format supports a variety of image projections including:
Any of the programs described in Software below decodes SIR headers.
In general, sir data files are generated using the scatterometer image reconstruction (SIR) resolution enhancement algorithm or one of its variants for radiometer processing. The multivariate SIR algorithm is a non-linear resolution enhancement algorithm based on modified algebraic reconstruction and maximum entropy techniques [Long, Hardin, and Whiting, 1993]. The singlevariate SIR algorithm was developed originally for radiometers [Long and Daum, 1997] but also used for SeaWinds [Early and Long, 2001]. The SIR w/filtering (SIRF) algorithm has been successfully applied to SASS and NSCAT measurements to study tropical vegetation and glacial ice (e.g. Long and Drinkwater, 1999). Variants of SIR have been successfully applied to the ERS-1/2 scatterometer and various radiometers (SSM/I and SMMR). (SIRF is used for SASS, NSCAT, and SeaWinds slice data processing. SIR is used for ERS-1/2 and SeaWinds egg data. The modified median filter [SIRF] is not used with ERS-1/2 data and SeaWinds egg data.)
For scatterometers, the multivariate form of the SIR algorithm models the dependence of sigma0 on incidence angle as sigma0 (in dB) = A + B * (Inc Ang - 40 deg) over the incidence angle range of 15 to 60 deg. The output of the SIR algorithm is images of the A and B coefficients. See the Data Characteristics section.
A represents the "incidence angle normalized sigma0" (effectively the sigma0 value at 40 deg incidence angle). The units of A are dB. Typically, +2 < A < -45 dB. However, in the SIR images A is typically clipped to a minimum -32 dB with values of A < -32 used to indicate 'no data'.
B describes the incidence angle dependence of sigma0 and has units of dB/deg. At Ku-band the global average of B is approximately -0.13 dB/deg. Typically, -0.2 < B < -0.1. B is clipped to a minimum value of -3 dB/deg. This value is used to denote 'no data' as well.
Single variable SIR or SIRF algorithms are used for radiometers and produce only an A (in this case, the brightness temperature) image. Typically, this can range from 165 to 320. Single variable SIR and SIRF algorithms are used for SeaWinds egg and slice images, respectively. In both cases the A images are at the nominal measurement incidence angle for the sensor and in the sensor measurement units.
Enhanced resolution images made from ERS-1/2 data use the Scatterometer Image Reconstruction (SIR) algorithm. This version of the algorithm does not incoporate a median filter and, for ERS-1/2, uses a 2-d Hamming window as the spatial response function for each beam. In the processing, a linear model relating sigma-0 and incidence angle is assumed, i.e. sigma-0(db) = A + B (theta - 40) where A is the "incidence angle normalized sigma-0" at 40 deg incidence in dB, B is the effective incidence slope of sigma-0 versus incidence angle in dB/deg, and theta is the incidence angle of the observation. This simple linear model is used in place of the gamma=sigma-0/cos(theta) as it more accurately represents the sigma-0 versus backscatter response over a wider range of surface and volume scattering conditions. The SIR algorithm makes images of A and B on an 8.9 km pixel grid. The effective resolution is estimated to be 20-30 km resolution, depending on region and sampling conditions. Raw ERS measurements have a quoted nominal resolution of 50 km on a 25 km sampling grid.
Language | Program Name | Description |
---|---|---|
C | csir_dump.c | dump SIR file to text output |
csir_dump_small.c | ||
csirexample.c | read SIR file, print values of corner pixels | |
sir_ez_example.c | ||
sir2bmp.c | convert SIR file to BMP | |
sir2byte.c | convert SIR file to raw, unsigned byte file | |
sir2gif.c | convert SIR file to GIF | |
sir2gif_ez.c | ||
sirmask.c | mask one SIR file with another | |
Fortran | fsir_dump.f | dump SIR file to text file |
fsir_dump_small.f | ||
fsir_locmap.f | read SIR file, create latitude and longitude maps like the auxiliary files | |
fsir_locmap_EZ.f | ||
fsirexample.f | read SIR file, create an unsigned byte file | |
fsirexample_EZ.f | ||
sir2byte.f | ||
sirmask.f | mask one SIR file with another | |
IDL | xsir_idl.pro | load SIR file, save to file, display image, do forward/inverse transforms |
PV-WAVE | xsir.pro, xsir_pvwave.pro | load SIR file, save to file, display image, do forward/inverse transforms |
MATLAB | loadsir.m, writesir.m, showimage.m, ... | load SIR file, save to file, display image, do forward/inverse transforms |
The IDL and PV-WAVE programs reside in one directory due to the similarity between the languages. xsir_idl.pro, xsir.pro, and xsir_pvwave.pro call the same functions, though the file loadsir.pro must be modified for PV-WAVE.
email: | podaac@podaac.jpl.nasa.gov |
url: |
http://podaac.jpl.nasa.gov/dataset/ERS-1_BYU_L3_OW_SIGMA0_ENHANCED
|
Dr. David Long of BYU is the source of this dataset. Please contact him with more detailed questions. See Investigator for contact information.
This data set is publicized courtesy of the PO.DAAC at JPL.
Early, D.S. and D.G. Long, Feb 2001. "Image Reconstruction and Enhanced Resolution Imaging From Irregular Samples," IEEE Transactions on Geoscience and Remote Sensing, Vol. 39, No.2, pp. 291-302.
Long, D.G. and D. Daum, 1997. "Spatial Resolution Enhancement of SSM/I Data," IEEE Transactions on Geoscience and Remote Sensing, Vol. 36, pp. 407-417.
Long, D.G. and M.R. Drinkwater, 1999. "Cryosphere Applications of NSCAT Data," IEEE Transactions on Geoscience and Remote Sensing, Vol. 37, No. 3, pp. 1671-1684.
Long, D.G., P. Hardin, and P. Whiting, 1993. "Resolution Enhancement of Spaceborne Scatterometer Data," IEEE Transactions on Geoscience and Remote Sensing, Vol. 31, pp. 700-715.