 



FRINGE 96An Interferometric QuickLook ProcessorAbstract
IntroductionWith the huge amount of interferometric SAR data that are now and potentially will be available from satellite SAR missions (ERS1, ERS2, RADARSAT, SIRC/XSAR III), the computing time to obtain interferograms and coherence maps in a semioperative scenario should be kept as short as possible. Thus, to select SAR pairs that actually can be exploited for interferometric applications the generation of fast, lowresolution surveys (e.g. fringes and coherence maps) should be carried out routinely on those pairs that show a suitable perpendicular baseline. A ``quicklook`` processor is presented in this paper. The processor computes interferogram and coherence maps in blocks of 33 x 100 km, at a geometric resolution 45 x 45 m on the ground (for flat earth). A good altimetric accuracy has been achieved by averaging the interferograms of 5 azimuth looks. Several adjacent ``blocks`` are then mosaicked to get a continuous azimuth strip. The processor performs near to real time on common lowcost machine: approximately 6 min are requested to compute a single block (32 x 100 km, 5 looks averaged) interferogram & coherence map, by using a PC PENTIUM PRO 200 running under Unix. Real time processing, i.e. 8 x the survey time (including some margin for data ingestion from HDDT), would require a machine 5 times faster. Efficiency has been obtained, at a cost of some quality loss, by designing fast techniques for:
Raw data presummingThe look formation is performed, before range and azimuth focusing, by bandpass filtering and subsampling raw data. This approach gives the following advantages:
The cost of presumming stage is minimized by using polyphase filters and implementing an ``integer'' format that handles one complex (I+Q) sample as a single integer word (16 bits) [2]. The processing time is thus reduced by a factor of 2.2 times due to range presumming (half band is dropped), 2 times due azimuth presumming (3 of the possible 8 looks are discarded) and 1.3 due to approximate focusing. These improvements are achieved at the cost of some information loss: nominally 50 % of the signal (for very low baselines) is lost due to range presumming and almost nothing for azimuth baseline in the tandem case due to the 300 Hz Doppler shift. Coherence estimateA fast algorithm has been designed (named quick & dirty) to compute coherence maps. The estimator achieves a computational time gain greater than 100, with respect of the conventional estimator, at the cost of a reduced statistical confidence (its accuracy being 3 to 5 times worse for coherence values > 0.4). The proposed estimator exploits the images amplitudes, and it is based upon the relation between the absolute value of the complex coherence,
(being the averaged interferometric phase), and the normalized crosscorrelation r of the detected SAR images:
The following expression [3]
can be exploited to compute the complex coherence, given an estimate of the normalized crosscorrelation of the amplitude image, r. The amplitude based coherence estimator is more efficient, since it does not require the knowledge of the interferometric phase (or local frequency). Moreover, it is not affected by possible local frequency estimation errors being insensitive to the interferometric phase. It can be shown that the proposed estimator is strongly biased by amplitude non stationarities, however this bias can be effectively reduced by applying a sort of automatic gain control (AGC) to both detected images before computing the coherence estimate [3]. The proposed estimator can be efficiently implemented to compute a coherence map by exploiting overlapped Bartlett or Boxcar windows at a cost of 50 flops for each image pixel, independently of the window size. Image coregistering & oversamplingThe estimate of coregistering parameters has been improved by means of the ``quick & dirty'' coherence estimator. Large windows can be explored since the estimator is independent of fringes & phases nonstationarities. Image coregistration is performed by two 1D resampling, since, for ERS, image rotation is well approximated by means of two 1D skews. Short space domain kernels are used for resampling & oversampling at the same time. Kernels are optimized for minimal phase noise and tabulated in steps of 1/100 pixel. A 67 samples kernel gives a coherence loss less than 2%. This approach is faster than frequency domain processing. Interferogram generationThe looks image are aligned in phase (to recover the effect of approximate resampling) and then coherently averaged. The multiple look interferogram is performed as usual by conjugate crossmultiplication. Finally the lowresolution interferogram (fringes) is generated by filtering and subsampling. Here again short kernels (6 samples) are used. After mosaicking the blocks to get a large strip, a residual range flattening is given to compensate for flat earth. Notice that these steps are not requested if the generation of coherence maps is the sole requirement. Computational complexityThe computing time achieved by the interferometric quicklook (for a 33 x 100 Km image block) is summarized in table 1.
The most timeconsuming steps have been highlighted in italics. Note the advantage achieved with respect to the standard  public domain interferometric processor, due to a much more efficient way to compute coherence. ConclusionsAn algorithm for generating interferograms & coherence strips has been presented. It is intended for data screening / browsing. All the processing steps requested for generating the interferogram have been revisited  efficiency has been achieved by keeping coherence loss due to processing less then 5%. A ``quickanddirty'', phaseindependent coherence estimator has been introduced. It requires 14'' to compute a 33 x 100 km, 5 looks averaged coherence map. A 5 looks averaged, 33 x 100 km ERS interferogram can be computed in 6' with a 30 Mflops/s Workstation. A prototype processor, installed at ESRIN, has produced several interferograms and coherence stripmaps more than thousand km long. AcknowledgmentsThe authors would like to thank the European Space Agency and Advanced Computer Systems (ACS) for sponsoring the work. Figure 1: Low resolution interferogram from ERS1/ERS2 tandem mission (jan '96), Greenland. Images' size is approximately 300 x 100 km. Left: absolute value; center: coherence map; and right: fringes. Figure 2: Enlargement of the previous image: an area of approximately 7 x 15 km has been shown. Note the fast fringes due to the movement of ice near the coast (iceberg ?). References
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. 